Previous Article | Next Article ![]()
Applied and Environmental Microbiology, September 2003, p. 5306-5318, Vol. 69, No. 9
0099-2240/03/$08.00+0 DOI: 10.1128/AEM.69.9.5306-5318.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Klinikum Grosshadern, Max von Pettenkofer Institute for Hygiene and Medical Microbiology, Ludwig Maximilian University Munich, 81377 Munich, Germany
Received 7 February 2003/ Accepted 6 June 2003
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Over the last few years, five new species (Enterobacter asburiae, Enterobacter cancerogenus, Enterobacter dissolvens, Enterobacter kobei, and Enterobacter nimipressuralis) have been reassigned to or split out of the nomenspecies E. cloacae (4, 31). E. dissolvens is the closest relative to E. cloacae, since the DNA relatedness of the two species was up to 82% in DNA-DNA hybridization experiments and they were phenotypically not differentiable (4, 25). E. nimipressuralis had a DNA relatedness of up to 67% to E. cloacae and differed phenotypically from it by being sucrose negative. E. asburiae had a DNA relatedness of 63% to E. cloacae and differed from it by being Voges-Proskauer, malonate, melibiose, and L-rhamnose negative (4). E. hormaechei also had a DNA relatedness of 63% to E. cloacae and differed from it by being D-sorbitol, melibiose, and esculin negative (31). E. cancerogenus is considered purely phytopathogenic and was transferred from the genus Erwinia to the genus Enterobacter as a senior synonym of Enterobacter taylorae (13). It had a DNA relatedness of 61% to E. cloacae and differed from it by being, inter alia, ornithine decarboxylase negative and D-arabinose positive. E. kobei had a DNA relatedness of 28% to E. cloacae and of 72% to CDC enteric Group 69. It differed from E. cloacae only by being Voges-Proskauer negative (24). Most of these species have a DNA relatedness of over 60% to E. cloacae and form the so-called E. cloacae complex.
Two DNA-DNA hybridization studies have so far partly analyzed the genetic clustering of the E. cloacae complex by examining local isolates with E. cloacae phenotypes (21, 25), delineating the following genogroups. One group contained both the E. cloacae and E. dissolvens type strains. A second formed around reference strain CDC 1347 41R and was not ascribed to a species (42). A third formed around the E. hormaechei type strain and contained seven biogroups. A fourth was found around the E. asburiae type strain and consisted of three biogroups. Finally, one group in the study of Grimont and Grimont (21) and three in that of Lindh and Ursing (25) as well as up to 11 ungroupable strains were not assigned to type or reference strains. However, the hybridization studies of Grimont and Grimont (21) and Lindh and Ursing (25) are difficult to compare. Both authors examined collections of strains recovered at local hospitals in Paris and Copenhagen, respectively. Beside others, this might be the reason for the partly different results of the authors. For example, Grimont and Grimont (21) reported 30% of study strains clustering with the E. asburiae type strain, while none of those of Lindh and Ursing (25) clustered with it. Grimont and Grimont (21) found a genetic group of nitrogen-fixing strains (group 5), whereas Lindh and Ursing (25) did not find such a group but reported almost 10% ungroupable strains. An extensive, population genetic-based classification of E. cloacae remains elusive.
Population genetic studies have proven to be extremely useful for getting insights into the epidemiology and evolution of the pathogenicity of bacterial species, as exemplified in previous studies on the species (3) and on the genus level (14). The objectives of this study were to describe the genetic structure of the nomenspecies E. cloacae and to explore the phylogenetic relationships between its genetic groups and related species. For this purpose, we analyzed 225 strains with E. cloacae phenotypes from 22 European regions and all type strains of the genus Enterobacter with a combination of partial sequencing of three housekeeping genes and PCR-restriction fragment length polymorphism (PCR-RFLP) of a further protein-coding gene.
| MATERIALS AND METHODS |
|---|
|
|
|---|
All 13 type strains of the genus Enterobacter were used as points of reference. They are listed in Table 1. Reference strain CDC 1347-71R (42) was included as a reference for DNA relatedness groups 2 of Lindh and Ursing (25) and Grimont and Grimont (21). Type and reference strains were obtained from the American Type Culture Collection and the Collection de l'Institut Pasteur (CIP). In addition, the published hsp60 sequence of Enterobacter aerogenes (accession number AB008141) and the hemB sequence of E. coli (accession number L44595) were downloaded from GenBank for direct comparison.
|
Biochemical and susceptibility testing.
API20E (BioMérieux, Marcy l'Etoile, France) testing was performed following the manufacturer's instructions. The results were interpreted with the Analytical Profile Index (API) database of the ApiLab Plus software (version 3.3.3; BioMérieux, Marcy l'Etoile, France). Antibiotic susceptibilities to ampicillin, ampicillin plus clavulanic acid, cefoxitin, ceftazidime, and cefotaxime were tested. The expression of the chromosomal AmpC cephalosporinase was interpreted as follows: no expression if the strain was susceptible to all antibiotics tested; basal expression if the strain was resistant to ampicillin, ampicillin-clavulanic acid, and/or cefoxitin, hyperexpression if the strain was resistant to ampicillin, ampicillin-clavulanic acid, cefoxitin, ceftazidime, and cefotaxime.
DNA preparation.
Bacterial DNA was prepared for PCR by quick-heat lyses. Bacteria were grown in LB to an optical density at 600 nm of 1.0. Then 100 µl of the culture was centrifuged (2 min at 10,000 rpm). The pellet was diluted in 1 ml of distilled water and boiled for 10 min.
ERIC-PCR analysis.
ERIC primers annealing to enterobacterial repetitive intergenic consensus regions were used as previously described (18). Amplification reactions were carried out in an Applied Biosystems GeneAmp 2700 PCR thermocycler. Thin-walled PCR tubes were purchased from Applied Biosystems (Branchburg, N.J.). Primers ERIC 1and ERIC 2 (Table 2) were used and synthesized by Metabion, Munich, Germany. Reactions were performed in a total volume of 30 µl under standard conditions (500 pmol of each primer, 200 µM nucleotide mix, 10 mM Tris-Cl [pH 8.3], 5 mM KCl, 1.5 mM MgCl2, and 2.5 U of standard Taq polymerase). After a 10-min heating period at 95°C, 35 amplification cycles were run with denaturation at 95°C for 1 min, annealing at 45°C for 1 min, and extension at 72° for 1 min followed by a final extension time of 15 min at 72°C. ERIC-PCR products were visualized after electrophoresis in 3% agarose gels in 1x TAE buffer, staining with ethidium bromide, and exposure to UV light. In each run, ERIC-PCR was repeated for strain EN114 in order to control the standardization of the conditions. The ERIC patterns of EN114 were identical for all three PCR runs performed. Patterns of different strains were compared by visual inspection. The patterns were interpreted as identical if an identical number of bands of identical sizes was found.
|
Single sets of primers could be used for Hsp60 and RpoB PCRs. Different sets of primers were needed for HemB PCR (sets HemB, HemBu, and HemB_all) and AmpC PCR (sets AmpC and AmpCu) because conserved regions were not identical in all sequences and for some clusters PCR products could not be achieved with one of the primer sets alone. PCR products of HemB, HemBu, AmpC, and AmpCu were sequenced unidirectionally with the forward primers for representatives of each strain cluster in order to verify that the correct target gene had been amplified. The hemB sequences obtained in this first step were aligned. Universal HemB_all primers were designed corresponding to their most homologous regions. PCR products were purified with the Nucleospin PCR cleaning kit (Macherey and Nagel, Dueren, Germany) before further processing.
PCR assays were performed with a Perkin-Elmer Gene Amp PCR System 2400 or an Applied Biosystems GeneAmp PCR System 2700 and Taq Gold DNA polymerase (Perkin-Elmer, Branchburg, N.J.) under standard conditions (see ERIC-PCR). Reactions were carried out in a total volume of 50 µl and started after a preheating time of 10 min at 95°C. Samples were subjected to 30 cycles of 30 s at 95°C, 30 s at annealing temperature, and 30 min of extension at 72°C. Cycles were followed by a final elongation step of 15 min at 72°C. Annealing temperatures were optimized by increasing the annealing temperature stepwise by 2°C, starting at 6°C below the lowest Tm of the primer sets used, until specific, single-banded PCR products were obtained. Due to the relatively high variability of the genes, PCRs of the ampC and hemB genes were only applicable at annealing temperatures of 50°C. Otherwise, the following annealing temperatures were chosen: 57.5°C for hsp60 and 59°C for rpoB. PCR products were visualized in 2.0% agarose gels after electrophoresis and ethidium bromide staining.
DNA sequencing.
PCR-amplified DNA was sequenced with the dye terminator method in both directions. To prepare cycle sequencing reactions, 0.2 µg of purified DNA was added to a 20-µl reaction solution containing 4 µl of BigDye dye terminator sequencing mix and 0.7 µl of a 5-pmol primer solution. Cycle sequencing was performed in the Applied Biosystems GeneAmp PCR System 2700. Cycle parameters were 25 cycles with an initial 96°C denaturation step of 10 s, followed by an annealing step at 57°C of 10 s and an extension step of 4 min at 60°C. Sequences were determined by electrophoresis with the ABI Prism 377 DNA sequencer. About 1% of base pairs sequenced resulted in a nontypeable character. They have been corrected, based on the reverse sequences, with the SeqMan (DNAStar) program.
PCR-RFLP.
In order to determine group-specific restriction patterns of ampC, PCR products were subjected to analysis of the PCR-RFLP patterns after restriction with HaeIII (Invitrogen, Groningen, The Netherlands). This restriction enzyme was chosen after computer analysis of the different ampC sequences of representative strains (EN30, EN117, EN124, E. asburiae type strain, and E. dissolvens type strain) with MapDraw (DNAStar Inc., 1999). Restriction was performed with 1 µg of DNA in a total volume of 20 µl over 4 h at 37°C. DNA was electrophoresed in 2% or 3% agarose gels, stained with ethidium bromide, and visualized with UV light.
Data analysis.
Sequence data were analyzed with MegAlign (DNAStar Inc., 1999) and PAUP version 4.0b10 (43) on a Macintosh G4 personal computer. Multiple alignment was performed with Clustal V, included in the MegAlign program. With MegAlign, pairwise percent sequence divergence was calculated as 100 times the sum of the residue distances from one strain to the other along the tree divided by the sum of all branch lengths. With Microsoft Excel, mean percent sequence divergences between clusters of sequences and the distinctness parameters k were calculated as previously described (32).
With PAUP 4.0b10, neighbor-joining trees (38) were estimated based on pairwise genetic distances on the basis of all substitutions with the Jukes-Cantor distance parameter. A group of strains was called a genetic cluster if the mean distinctness parameter k was above 2.0; otherwise it was called a sequence crowd. The significance of branchings was evaluated by bootstrap analysis of 100 replicates. Only groups with frequencies above 50% were kept. Bootstrap values above 80% were considered significant, and the cluster was designated as being strongly supported.
Nucleotide sequence accession numbers.
The GenBank accession numbers for the sequences reported in this paper are: for hsp60, AJ417108 to AJ417143, AJ543761 to AJ543910, and AJ567846 to AJ567902; for rpoB, AJ543656 to AJ543728 and AJ566948 to AJ566941; for hemB, AJ426490 to AJ426522 and AJ567697 to AJ567724.
| RESULTS |
|---|
|
|
|---|
To eliminate this risk, ERIC-PCRs were performed for all isolates. The good specificity of ERIC-PCR for the detection of identical isolates has recently been demonstrated (18). Comparing the ERIC patterns obtained by visual inspection, 35 isolates showed patterns which were identical to those of other strains (data not shown). They were excluded. In addition to the remaining 206 study strains, all 13 type strains of the genus Enterobacter and reference strain CDC 134771R, which represents genetic group 2 of Grimont and Grimont (21), were included. All 220 strains and their characteristics are listed in Table 1. They were isolated from blood cultures (26 strains), the respiratory tract (59 strains), urogenital tract (38 strains), abdomen (9 strains), superficial swabs (44 strains), artificial devices (7 strains), stool (7 strains), plants (7 strains), and unknown origins (23 strains).
Molecular systematics analysis.
Despite their public health significance, only very few sequence data for housekeeping genes were available for Enterobacter species at public databases. Hence, we tried oligonucleotides for the amplification of the gyrase A gene (gyrA), gyrase B gene (gyrB), or subunit C of topoisomerase IV (parC) published for other species (5, 10), but PCR products were not obtained with any of them for all of the study strains. However, we obtained PCR products for all strains by amplifying hsp60, rpoB, and hemB.
Molecular phylogeny based on hsp60 sequence analysis.
The sequences of a portion of 273 nucleotides of the hsp60 gene were determined for the study set of 220 strains by bidirectional sequencing. All sequences were pairwise aligned; 91 characters (33%) were variable (defined as different in at least one strain). A neighbor-joining tree was constructed with all 220 study strains (Fig. 1). Since E. aerogenes is a homotypic synonym to Klebsiella mobilis and belongs to another genus, its hsp60 sequence was used to root the tree; 13 groups of strains were found.
|
|
In order to get a better understanding of the distances of the clusters, the mean numbers of nucleotide differences along the tree were calculated as parameters of distance (Table 3). Mean percent sequence divergences within the clusters were equal to or below 1.0 (0 to 1.0 ± 0.6), with the exception of cluster I (3.8 ± 2.3), which was due to a relatively high divergence of the E. asburiae type strain to the other strains of the cluster. Sequence crowd xiii had a within group divergence of 6.7 ± 1.3. Divergences between the clusters were above 2.5 (2.7 ± 0.5 to 9.1 ± 0.1) with the exceptions of cluster pairs VI-VII, VI-VIII, and VII-VIII.
Some clusters showed specific amino acid alterations in the 60-kDa heat shock protein coded by the hsp60 gene (Table 1). For example, cluster V had a specific amino acid alteration at position 465, clusters IV and VI at position 468, cluster XI, including the E. cloacae type strain, at positions 430 and 466, and cluster XII, including the E. dissolvens type strain, at position 430 [E. coli amino acid numbering system (M11294)]. Interestingly, the E. asburiae type strain differed from the other strains of cluster I by an alteration at position 466.
Summarizing, hsp60 analysis suggested the existence of at least 12 genetic clusters and one heterogeneous sequence crowd within the E. cloacae complex. Eleven clusters (II to XII) were supported by bootstraps, and eight of them had strong bootstrap support. Clusters VII and VIII were well supported at a higher node.
Molecular phylogeny based on rpoB sequence analysis.
The sequences of a second housekeeping gene were investigated in order to reconfirm the clustering observed after hsp60 analysis. For this purpose, a representative subset of 78 strains was chosen, including all type and reference strains of the E. cloacae complex and the type strain of E. gergoviae. Sequences of a portion of 1,008 nucleotides of the rpoB gene were determined for these strains by bidirectional sequencing. About 700 nucleotides overlapped from both directions, leaving about 150 bp at each end sequenced unidirectionally. Multiple alignment by Clustal V revealed a total of 218 (22%) variable characters. The neighbor-joining relationships of the sequences are presented in Fig. 2 (left side). With the only exception being sequence crowd xiii, the rpoB tree was in strong agreement with the hsp60 tree. All strains were assigned to the same clusters. Moreover, much sharper delineations were observed, and bootstrap values were higher for clusters I, VIII, and IX. All k parameters but those of sequence crowd xiii were well above 2.0 (2.8 to 99), indicating that the rpoB sequence clusters could be considered distinct genetic groups (Table 2). However, in contrast to the hsp60 results, sequence crowd xiii was split into two groups (xiii.a and xiii.b). Correspondingly, its mean k parameter was well below 2.0 (1.6 ± 0.4). Recapitulating, the same genetic clusters observed for hsp60 appeared for rpoB. Sequence crowd xiii was unstable and fell apart.
|
|
The ampC gene was amplified with one of the two PCR assays presented and restricted with HaeIII for all strains included in the rpoB analysis (Fig. 2). By visual inspection of the resulting PCR-RFLP patterns, 10 different clusters were delineable, designated a to k. Three pairs of sequence clusters showed similar restriction patterns (I and IV, III and VII, and VI and VIII). Four strains of cluster VIII showed patterns aberrant from the others, demonstrating a little within-group variability. The similarity of the restriction patterns of clusters III and VII was striking: whereas the distance of cluster VII to cluster III was much higher than that to clusters VI and VIII in the hsp60 and hemB trees, the distances were just the other way around in the ampC PCR-RFLP and rpoB tree. The E. asburiae type strain differed from the other strains of cluster I by an additional restriction site in the ampC pattern, underlining the relatively high within-group divergences of cluster I due to the distance of the E. asburiae type strain from the other members of the group (Fig. 1, Fig. 2, Fig. 3, and Table 3).
Whereas the PCR-RFLP patterns were similar for all strains of cluster XII, no PCR products were obtained for the strains of sequence cluster XI, including the E. cloacae type strain, with both PCR assays used. This supported the clustering of the hsp60 and rpoB trees, in which E. cloacae and E. dissolvens represented two distinct phylogenetic lineages, whereas they were mixed in the hemB tree.
Do different genetic clusters show specific distributions to clinical materials?
To some extent, genetic clusters were not randomly distributed to the different groups of clinical materials from which the strains were recovered. Some features should be outlined. First, E. nimipressuralis (cluster X) was not isolated from a human specimen. Second, 15% of strains of cluster VIII (59 strains) were recovered from blood cultures representing more than a third of all blood culture isolates. In contrast, only 3% of strains of the closely related cluster VI (28 strains) were recovered from blood cultures. Third, whereas none of the 58 strains of cluster III were recovered from abdominal specimens, 40% of studied strains of cluster II were recovered from abdominal specimens, representing almost half of all abdominal isolates. Fourth, whereas about 20% of strains of all clusters derived from superficial skin or wound swabs, none of the strains of clusters I and II did so. These data suggest that the different clusters could have different virulence-associated properties. However, since only strains with E. cloacae phenotypes were considered in this study, the other species might have been underrepresented. This might have distorted the picture. More appropriately designed studies are needed in order to evaluate the different clinical relevances of the clusters.
API20E analysis.
Twenty-one different API20E codes were generated by the 206 study strains included. The predominating codes were 3.305.573 (58%) and 3.305.773 (19%). Only a few group-specific characteristics could be drawn out of the biochemical information from the API20E analysis. First, cluster IX was 100% myoinositol positive, whereas clusters VI and VII were 100% negative. Second, cluster VII was 100% sorbitol negative, but sorbitol-negative strains were also found in clusters I, III, VI, and VIII. Third, cluster IV was 100% rhamnose negative, which was otherwise found only in the E. asburiae type strain. This short analysis suggested that at least some of the clusters could have specific biochemical properties.
| DISCUSSION |
|---|
|
|
|---|
Six species (E. asburiae, E. cancerogenus, E. dissolvens, E. hormaechei, E. kobei, and E. nimipressuralis) which are closely related to E. cloacae are subsumed in the so-called E. cloacae complex; 13% of our study strains clustered around their type strains, whereas 84% formed genetic clusters which were not associated with any type or reference strain. Only 3% of our study strains clustered with the type strain of E. cloacae. In DNA-DNA hybridization studies, Grimont and Grimont (21), who examined 49 clinical isolates, and Lindh and Ursing (25), who examined 123 clinical isolates, reported similarly low frequencies of true E. cloacae strains. It has been known for years that the choice of the type strain of E. cloacae was unfortunate, since its genotype does not represent the majority of what are generally considered E. cloacae strains (21).
For this reason, many authors preferred to use reference strain CDC 1347-71R as representative of E. cloacae for their DNA-DNA hybridization experiments (4, 31, 37); 5% of our study strains grouped around CDC 1347-71R which complied with the rate reported by Lindh and Ursing (25). In our study, the type strain of E. kobei fell into the same cluster. E. kobei was described based on a small group of clinical isolates selected as being Voges-Proskauer negative (24). The species showed close relatedness to CDC enteric group 69 (16). Data on the relatedness of strain CDC 1347-71R and CDC enteric group 69 are not available in the literature. CDC 1347-71R and our study strains clustering with the E. kobei type strain were all Voges-Proskauer positive. The following two problems are raised by these findings and need to be addressed in further DNA-DNA hybridization studies. First, if reference strain CDC 1347-71R in fact represents E. kobei and not E. cloacae, the relatedness of the species E. hormaechei, E. asburiae, E. dissolvens, and E. nimipressuralis to E. cloacae needs to be reevaluated. Second, if CDC1347-71R, CDC enteric group 69, and E. kobei represent the same relatedness group, the species description of E. kobei needs to be emended. Forty percent of our strains of the E. kobei cluster were isolated from intra-abdominal sites, and 45% of all intra-abdominal isolates belonged to the E. kobei cluster. Whether this is just an odd coincidence needs to be investigated in an appropriately designed epidemiological study including more clinical strains.
We found that 2% of our study strains clustered with the E. dissolvens type strain. The species E. dissolvens was reassigned from the genus Erwinia to the genus Enterobacter with the reservation that it might in fact be identical to E. cloacae (4). Strains of E. dissolvens and E. cloacae were 60 to 82% related and had
Tm values far below 5°C in reciprocal DNA. hybridization studies (21, 25). In our study, E. cloacae and E. dissolvens strains were closely related to each other but still represented two distinct phylogenetic lineages in three of the four genes analyzed. Extended multilocus sequencing and DNA-DNA hybridization studies with larger sets of strains are needed in order to finally decide whether both species denominations should be maintained or whether one of them should be abandoned.
We found that 3% of our study strains clustered around the E. hormaechei type strain. The description of the species E. hormaechei was based on 23 strains which were selected on the criterions that they were similar to E. cloacae but did not produce acid from D-sorbitol and melibiose. They formed a distinct DNA-DNA hybridization group and had 63% DNA relatedness to reference strain CDC 1347-71R (31). In the hybridization study of Gimont and Grimont (21), the largest cluster formed around the E. hormaechei type strain and consisted of seven biotypes with only one (3a, one strain) matching the original description of the species. Later, Davin-Regli et al. (11) reported outbreaks with E. hormaechei strains that were D-sorbitol and melibiose positive, suggesting that additional biotypes might exist.
In our study, two study strains together with the E. hormaechei type strain formed the small genetic cluster VII. As far as possible with the information from the API20E system, the biochemical properties of this cluster largely met the species definition of E. hormaechei (31). These data suggest that the species E. hormaechei is actually a valid description (31) and represents only a small phylogenetic lineage within a larger genogroup (21). Our clusters VI and VIII were closely related to E. hormaechei cluster VII. DNA-DNA hybridization studies are needed to verify whether these clusters form a common DNA relatedness group allowing emending and broadening of the species description of E. hormaechei.
We found that 4% of our strains clustered with the type strain of E. asburiae. The species E. asburiae was deduced from the former enteric group 17, which formed a single DNA relatedness group and consisted of Enterobacter strains which were Voges-Proskauer, melibiose, and L-rhamnose negative (4). In previous studies, this biotype represented only one of at least three different biotypes found in the DNA-DNA hybridization group around the E. asburiae type strain (21). In our hsp60 analysis, cluster I was only poorly supported by bootstraps due to the E. asburiae type strain, which grouped slightly apart the cluster. In the ampC PCR-RFLP, the E. asburiae type strain had an additional restriction site. None of the strains of cluster I except the E. asburiae type strain was Voges-Proskauer, melibiose, or L-rhamnose negative. These data reflect the geno- and phenotypic heterogeneity of E. asburiae. In former studies, it played only a minor role as a clinical pathogen (9, 12, 40). In our study, one strain of the E. asburiae cluster (EN373) was isolated from the blood culture of a neonate (18). Many more strains need to be analyzed in order to get a better understanding of the species E. asburiae.
E. nimipressuralis was isolated from elm trees with a disease called wetwood. The species was transferred from the genus Erwinia to the genus Enterobacter and was 52 to 67% related to E. cloacae in reciprocal DNA hybridization reactions (4). Until now, it has never been described as a human pathogen. In our study, only one isolate from a dog's footpad clustered with the E. nimipressuralis type strain, reconfirming the suggested low clinical relevance of this species. Little is known about the species E. cancerogenus. First ascribed to the genus Erwinia, it has recently been transferred to the genus Enterobacter (13) as a senior synonym of Enterobacter taylorae, which was 61% related to E. cloacae reference strain CDC 1371-71R in reciprocal DNA hybridization reactions (17). Human infections with E. cancerogenus are incidental (1), which was also reflected by the fact that none of our strains clustered with the E. cancerogenus type strain.
All our study strains were identified as E. cloacae by their phenotypic characteristics. Nonetheless, 13% of them belonged to other species of the genus. The species of the E. cloacae complex were without exception described based on phenotypically selected strain collections, of which the genetic relatedness was proven in a second step. From our findings, it follows that the species descriptions of presumably all species of the E. cloacae complex need to be revised and potentially emended.
The value of sequence data in population genetic studies has been more and more acknowledged (27). Since its advent, the outstanding role of the analysis of 16S rDNA sequences for the determination of phylogenetic relationships between prokaryotes is indisputable. However, for closely related species such as Enterobacteriaceae, the role of the 16S rDNA sequence analysis in population genetic studies has been questioned (29, 32). For the E. cloacae complex in particular, its limitations have been demonstrated (44). We analyzed the genetic relationships of 16S rDNA sequences of E. cloacae strains in a limited prestudy (data not shown). Only vague genetic clusters were generated. In view of this restricted discriminatory power of the 16S rDNA, we decided to focus on analysis of protein coding genes (45).
The analyses of four protein coding genes, hsp60, rpoB, hemB, and ampC, were included in the study. Sequencing of partial hsp60, also known as groEL or cpn60, has been successfully applied for the classification of many bacteria and turned out to be useful for the phylogenetic analysis of Enterobacter in the present study (19, 20, 28, 33, 45). Mollet et al. (29) demonstrated that phylogenetic trees generated on the bases of rpoB sequences were well compatible with the currently accepted classification of Enterobacteriaceae. In this aspect, it was far superior to 16S rDNA, a fact that we could reproduce in our study. Anyway, housekeeping genesas most of the genes typically sequenced by bacterial systematistsare not likely to be involved in the ecological delimitation of bacterial populations (32). ß-Lactamases evolved as defense mechanisms against competitors. They directly affect the competitive superiority of the bacterial strains in their ecological niches (8). ampC, coding for an Ambler class C cephalosporinase, is coded on the chromosome and is universally present in the genus Enterobacter. The clonality of ampC genes on the species level (including E. asburiae and E. hormaechei) has recently been shown (36). Following our PCR-RFLP results, ampC genes seemed to be clonal even on the level of genetic groups.
The population genetics presented provide a framework of the genetic structure of the E. cloacae complex. This may be of value for further studies examining the epidemiological, phenotypic, and virulence-associated properties of clones of the E. cloacae complex. In Escherichia coli, it has been demonstrated that factors associated with virulence are not randomly distributed among genetic groups (3, 7). Some virulence-associated properties have been described for strains of the E. cloacae complex (23). Clonality of virulence factors might cause the unequal distribution of the genetic clusters to the clinical materials which we noticed in our study. DNA-DNA hybridization experiments and biotyping of our strain collection will help determine to what extent our sequence clusters correspond to distinct species. Additionally, this will allow us to correlate sequence data for protein coding genes with DNA reassociation data, as exemplified in previous studies for 16S rDNA (41).
| ACKNOWLEDGMENTS |
|---|
We thank Kleinhuber and ESCMID for the publication of our call for strains. Special thanks to all colleagues who collected and sent strains for this studyS. Alberti, Palma de Mallorca; F. Allerberger, Innsbruck; J. Bille, Lausanne; D. Bitter-Suermann, Hannover; K. Bottolfsen, Bergen; M. Breuer-Wera, Aachen; A. Dierkes-Kersting and P. Breuer, Gelsenkirchen; H. Erichsen, Kiel; H. K. Geiss, Heidelberg; P. B. Heczko, Cracow; A. Leanord, Lanarkshire, United Kingdom; S. Lukas, Regensburg; C. E. Nord, Stockholm; W. Pfister, Jena; D. Pierard, Brussels; K. Poschinger, Munich; D. P. Roberts, New Brunswick; V. Schäfer, Frankfurt; R. Smyth, Växjö; M. Stark and I. Authenrieth, Tübingen; J. Wagner, Berlin; A. Wenger, Lausanne; Zwilling; and M. Kist, Freiburgand all other colleagues who have not specified their names but provided us with strains. We are grateful to Gudrun Maindok for excellent technical assistance in laboratory work.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||