Genomic Diversity within the Enterobacter cloacae Complex

Abstract

Background

Isolates of the Enterobacter cloacae complex have been increasingly isolated as nosocomial pathogens, but phenotypic identification of the E. cloacae complex is unreliable and irreproducible. Identification of species based on currently available genotyping tools is already superior to phenotypic identification, but the taxonomy of isolates belonging to this complex is cumbersome.

Methodology/Principal Findings

This study shows that multilocus sequence analysis and comparative genomic hybridization based on a mixed genome array is a powerful method for studying species assignment within the E. cloacae complex. The E. cloacae complex is shown to be evolutionarily divided into two clades that are genetically distinct from each other. The younger first clade is genetically more homogenous, contains the Enterobacter hormaechei species and is the most frequently cultured Enterobacter species in hospitals. The second and older clade consists of several (sub)species that are genetically more heterogonous. Genetic markers were identified that could discriminate between the two clades and cluster 1.

Conclusions/Significance

Based on genomic differences it is concluded that some previously defined (clonal and heterogenic) (sub)species of the E. cloacae complex have to be redefined because of disagreements with known or proposed nomenclature. However, further improved identification of the redefined species will be possible based on novel markers presented here.

About the Authors:

Armand Paauw

* E-mail: [email protected]

Affiliation: Department of Medical Microbiology, University Medical Centre Utrecht, Utrecht, The Netherlands

Martien P. M. Caspers

Affiliation: TNO Department of Microbiology, Zeist, The Netherlands

Frank H. J. Schuren

Affiliation: TNO Department of Microbiology, Zeist, The Netherlands

Maurine A. Leverstein-van Hall

Affiliation: Department of Medical Microbiology, University Medical Centre Utrecht, Utrecht, The Netherlands

Alexis Delétoile

Affiliation: Unité Biodiversité des Bactéries Pathogènes Emergentes, Institut Pasteur, Paris, France

Roy C. Montijn

Affiliation: TNO Department of Microbiology, Zeist, The Netherlands

Jan Verhoef

Affiliation: Department of Medical Microbiology, University Medical Centre Utrecht, Utrecht, The Netherlands

Ad C. Fluit

Affiliation: Department of Medical Microbiology, University Medical Centre Utrecht, Utrecht, The Netherlands

Introduction

Enterobacter cloacae is a facultative anaerobic Gram-negative bacillus belonging to the family of Enterobacteriaceae. The nomenclature (taxonomy) of the E. cloacae complex is mainly based on whole genome DNA-DNA hybridizations and phenotypic characteristics [1]–[3]. Currently, 6 species have been assigned to the Enterobacter cloacae complex, including Enterobacter asburiae, Enterobacter cloacae, Enterobacter hormaechei, Enterobacter kobei, Enterobacter ludwigii, and Enterobacter nimipressuralis. Only Enterobacter isolates that belong to the E. cloacae complex are of clinical significance and are increasingly isolated as nosocomial pathogens [4], [5]. In surveillance studies, Enterobacter species are often not further classified beyond the genus level probably because identification is difficult. Enterobacter spp. causes 7% of nosocomial infections in intensive care units in the USA [6], [7].

Accurate species identification is desirable for determining whether specific species within the E. cloacae complex are more prone to cause infections. More precise identification of E. cloacae complex isolates may permit differentiation between nosocomial species and commensal or transitional species. However, until now phenotypic identification of species and subspecies within the E. cloacae complex have been largely unreliable and irreproducible [8]. In this study, the discriminatory power of four genetically-based approaches was evaluated. These methods were hsp60- and rpoB-genotyping, multi-locus sequence analysis (MLSA), and comparative genomic hybridization (CGH).

The first method entailed genotypic identification via sequencing of a fragment of the heat shock protein 60 gene (hsp60), and was included because hsp60 genotyping appeared to be a promising novel method. Using this approach, the E. cloacae complex was divided into 12 genetic clusters (I–XII) and an unstable sequence cluster (XIII),(Table 1) [9].

[Figure omitted. See PDF.]

Table 1. The different clusters of the Enterobacter cloacae complex and their nomenclature, according hsp60 genotyping by Hoffmann and Roggenkamp [9].

https://doi.org/10.1371/journal.pone.0003018.t001

The second method consisted of sequencing a fragment of rpoB, and possibly represents a good alternative for strain identification because of the high resolution available to differentiate between closely related species [10].

The third approach employed MLSA, a more sophisticated technique. With this method, several widely separated genomic loci are analyzed, which minimizes the effect of recombination events on clustering analysis.

Finally, microarray-based CGH was also used; this approach is a powerful method for performing genome-wide studies on different bacteria [11]–[19].

The E. cloacae complex could thus be divided into two genetically distinct clades that were not previously recognized. The younger first clade contains the E. hormaechei species. This species is possibly more clinically relevant because it represents the most frequently cultured Enterobacter species in hospitals. The second clade of the E. cloacae complex consists of several clonal and heterogeneous (sub)species, leading us to propose a redefinition of species assignment within this genomically diverse complex.

Results

Identification of isolates based on hsp60 sequences

The population structure of the E. cloacae complex in hospitals was studied using 158 isolates. Of these, 120 (including 27 outbreak I isolates and 13 commensal fecal isolates) were obtained in the University Medical Centre Utrecht, Utrecht, The Netherlands (UMCU); 21 isolates (including 10 outbreak I isolates) from 11 other Dutch hospitals; 12 isolates from 12 different European hospitals; 3 isolates from an industrial site, and 2 isolates were E. cloacae ATCC13047. (See Materials and Methods for more information about the isolates). All isolates were identified as E. cloacae using the Phoenix Automated Microbiology System (Becton Dickinson Biosciences, Sparks, MD, U.S.). These results were compared to identification based on sequencing of a previously described 273 bp fragment of hsp60 that appeared to give reliable identification of E. cloacae complex isolates [9]. Twenty-nine different sequences were obtained and compared with sequences described by Hoffman and Roggenkamp [9]. The isolates were identified as E. cloacae III (n = 26), E. cloacae IV (n = 13), E. cloacae IX (n = 1), E. cloacae cloacae (n = 4), E. cloacae dissolvens (n = 1), E. hormaechei oharae (n = 13), E. hormaechei steigerwaltii (n = 78), Enterobacter kobei (n = 7), E. ludwigii (n = 6), and E. asburiae (n = 9) (Table S1). The hsp60 genotyping data indicate that phenotypic identification using the Phoenix Automated Microbiology System cannot discriminate between different species and subspecies of the E. cloacae complex.

Identification of isolates based on rpoB sequences

Because identification of Enterobacteriaceae by sequencing of a fragment of rpoB has been reported as a robust and promising method for identification [19], isolates were also classified based on the rpoB sequence. Thirty-two different sequences were detected in the set of 158 isolates after analysis of a 501 bp fragment of rpoB. Sequences were compared with a database of rpoB sequences that represents the major phylogenetic clades in the family Enterobacteriaceae [20]; this database is managed by the Pasteur Institute in Paris, France. The isolates were identified as E. hormaechei (n = 117), Enterobacter asburiae (n = 10), E. cloacae (n = 5), Enterobacter cloacae IV (n = 20), and Enterobacter ludwigii (n = 6) (Table S1). As with hsp60 genotyping, the results of rpoB typing showed that phenotypic identification by the Phoenix Automated Microbiology System cannot discriminate between different species and subspecies of the E. cloacae complex. Furthermore, the results also show discrepancies between the two single gene genotyping methods (Table S1). These discrepancies may at least in part be explained by recombination events.

Reliability of DNA fragments for MLSA analysis of the Enterobacter cloacae complex

To minimize the effects of potential recombination events between isolates, MLSA was used. Gene fragments of five additional genes (fusA, gyrB, leuS, pyrG, and rplB) of 50 E. cloacae complex isolates were sequenced. These isolates where selected based on their differences in hsp60 sequence and a preliminary analysis of the CGH data. For unique hsp60 sequence, one isolate was selected. In addition, isolates with similar hsp60 sequences but different preliminary CGH clustering results were included with the limitation that only one isolate was selected from each cluster (data not shown). The accession numbers of the sequenced gene fragments are listed in Table S2. Initial analysis of the sequence data revealed that the E. cloacae complex consists of 2 clades. To determine whether the sequence data was appropriate for phylogenetic studies, we assessed whether positive selection might have affected our data using two population genetic tests: Tajima's D statistic [21] and Fu's Fs statistic [22], [30] (Table 2). Tajima's D statistic only produced a significant (p<0.05) result for the rplB gene from clade 1. This indicates that the rplB gene might be under positive selection. This finding is supported by a significant result in Fu's F test (p<0.01). In addition, the transition/transversion ratio of rplB was estimated by plotting the transversion and transition rates as a function of the genetic distances (data not shown). The data show no saturation for either curve. However, the transition/transversion ratio was estimated to be 1, indicating a non-random distribution of mutation frequency. This supports the notion that the rplB gene is under positive selection. The only other significant Fu's F test was for the rpoB fragment from clade 2, indicating that these isolates are possibly under positive selection or that this population expands. However, this result was not supported by Tajima's D test (Table 2) or the Ka/Ks ratio (data not shown). Therefore, positive selection of rpoB is not likely. These results indicate that all gene fragments except for rplB are suitable for phylogenetic analysis of the E. cloacae complex.

[Figure omitted. See PDF.]

Table 2. Polymorphisms observed in seven protein-encoding genes among 50 E. cloacae complex isolates.

https://doi.org/10.1371/journal.pone.0003018.t002

The quality of the phylogenetic information obtained for the other remaining genes was also analyzed by plotting the transversion and transition rates as a function of the genetic distance (Figure 1). Transversion and transition rates increase concomitantly with genetic distance, indicating that there is little or no saturation. The transition/transversion ratio was estimated to be 2. These data further support that the sequences of the remaining genes can be used for phylogenetic studies.

[Figure omitted. See PDF.]

Figure 1. Transitions/transversion vs genetic distance.

From the six concatenated sequences used for phylogentic studies the transitions (blue crosses) and transversions (green triangles) are plotted against genetic distance according to the K80 model. No saturation is detected and the transitions/transversion ratio is estimated at 2.

https://doi.org/10.1371/journal.pone.0003018.g001

Phylogenetic analyses

The genetic heterogeneity of the E. cloacae complex was analyzed by generating phylogenetic trees based on the MLSA data. Noncongruence between trees based on the sequences of single DNA fragments indicated that recombination events had occurred (data not shown) [23]. Genetic relationships between the isolates could not be determined by E-burst or minimal spanning tree, because the concatenated sequences of the six gene fragments contained more than two mutations per sequence. Therefore, ClonalFrame was used to determine the genetic relationship between the tested isolates. The resulting tree contained the sequences of 51 isolates, including 50 E. cloacae complex isolates and one Escherichia coli K12 that was used as the outgroup (Figure 2A). For reasons mentioned previously, the sequences of the rplB gene were not included in the final analysis. The nodes of several branches were different when using maximum likelihood clustering with concatenated sequences of the six gene fragments (Figure 2B). However, isolates that clustered together in the same branches of the ClonalFrame tree also clustered in the maximum likelihood-tree. This indicates that the genetic relationship between the clusters is perhaps uncertain but that isolates in each cluster are genetically similar to each other because the same isolates cluster together in both trees. The recombination/mutation ratio was 1.04 (95% confidence interval 0.72–1.45), indicating that recombination events are involved in the evolution of the E. cloacae complex. This may explain the differences between the phylogenetic trees. Based on the ClonalFrame tree the isolates, except one were divided into 7 separate clusters (Figure 2A+B).

[Figure omitted. See PDF.]

Figure 2. Phylogenetic trees of E. cloacae complex.

A) Phylogenetic tree based on results of MLSA clustered with Clonalframe. Numbers indicate confidence values of the branches. * Solitaire isolate. K12: E. coli K12 used as the outgroup. B) Phylogenetic tree based on concatenated sequences of six gene fragments clustered with the Maximum Likelihood. Each color depicts a cluster according to the CGH results. * Solitaire isolate. K12: E. coli K12 used as the outgroup.

https://doi.org/10.1371/journal.pone.0003018.g002

Comparative genome hybridization analysis

The six DNA fragments used for MLSA present only a small fraction of the genome. To evaluate whether MLSA results are representative of population structure on the genomic level, comparative genome hybridization (CGH) was performed. CGH was performed with 3,072 DNA fragments from a shotgun library obtained from 8 E. cloacae complex isolates (7 strains); these fragments functioned as probes (mixed genome array or MGA). After analyzing all 158 slides, a total of 2,614 spots (85%) met the quality criteria and were included in the study (Table S3). The average size of the DNA fragments spotted on the slides was 1202 bp. Eighteen randomly selected DNA fragments of the core genome (core genome was defined as genes present in all isolates with a likelihood of >95%) were sequenced (Table S4). The average size of fragments representing the core genome was 1178 bp. No redundancy was detected in these sequenced DNA fragments. The overall redundancy of the MGA was calculated to be less than 5.3%, thus a maximum of 138 of the 2,614 DNA fragments have a sequence overlap with another insert. The total genetic coverage of each of the seven different strains used for construction of the array was calculated [16], [24] (Table 3). The gene coverage varied between 55% and 60% for the 7 isolates and the pan-genomic coverage of the E. cloacae complex was estimated to be 55% to 60%. The calculated minimal nucleotide coverage of the core genome was 58%. Results show that more than half of the genes, but also non-encoding or other coding DNA sequences, are represented in the constructed MGA. For inserts EnCl.014E12 and EnCl.020H03, 98% and 97% of the PCR results correlated with the results of the CGH array, respectively. The high similarity between CGH and PCR results indicate that MGA hybridization results were consistent. All replicates were found in the same cluster as the originals. Clearly, CGH array data are reliable for comparative genomic studies.

[Figure omitted. See PDF.]

Table 3. Characteristics of isolates used for construction of DNA library and the calculated gene coverage of each isolate.

https://doi.org/10.1371/journal.pone.0003018.t003

Of the 2,614 spots analyzed, 1,358 were considered to contain core genome DNA fragments because they were present in all isolates tested. These fragments were omitted for further analysis because they do not contribute to the speciation of the E. cloacae complex. Principal component analysis (PCA) showed that the isolates were mainly divided on the basis of the absence or presence of the pQC plasmid (data not shown). Based on the PCA data and because pQC is approximately 300 kb in size, it was expected that pQC-derived DNA fragments in the MGA would influence the assignment of isolates to different species or clades. Statistical and clustering analysis (described in Materials and Methods) showed that 106 DNA fragments were likely derived from the pQC plasmid. These 106 plasmid-derived DNA fragments were removed from further analysis. PCA of the hybridization results of the remaining 1,150 DNA fragments showed that the E. cloacae complex is divided into two genetically distinct clades (Figure 3A). Subsequently, these results where confirmed with a phylogenetic tree created with hierarchical clustering (not shown). Using the Pars program of the PHYLIP software package, a parsimony tree was computed (Figure 3B) that confirmed the genetic distinction between the two clades. Isolates that were not previously typed by MLSA were clustered based on the results of the constructed parsimony tree (Figure 4). Clusters 1–2 represent isolates from the first clade, whereas isolates from the second clade belong to clusters 3–7 and the solitaire isolate. Clusters 1–7 contained 25, 91, 6, 6, 9, 7, and 13 isolates, respectively.

[Figure omitted. See PDF.]

Figure 3. Analyze of the comparative genome hybridization results.

Comparative genome hybridization data clustered with PCA and maximum parsimony. Comparison of CGH data of 180 slides tested for 1150 DNA fragments A) with PCA; B) with the Parsimony method. Both figures show that isolates of the E. cloacae complex are divided in two separate clades. Red: first clade with 116 different isolates; Green: second clade with 42 different isolates.

https://doi.org/10.1371/journal.pone.0003018.g003

[Figure omitted. See PDF.]

Figure 4. Tree constructed on basis of the CGH results of the Parsimony method.

The slanted cladogram depicts genetic relationships between isolates each color represents a different MLSA cluster. *Solitaire isolate in MLSA.

https://doi.org/10.1371/journal.pone.0003018.g004

Comparison of the results of the different typing techniques

Cluster 1 isolates were identified as E. hormaechei using rpoB sequencing, but according to hsp60 genotyping, all isolates were E. cloacae III (Table S1). This cluster was a separate branch in the first clade of both trees on CGH and MLSA results. Cluster 2 isolates were all E. hormaechei isolates according to rpoB, but E. hormaechei steigerwaltii and E. hormaechei oharae according to hsp60 genotyping. These isolates were highly similar according to MLSA, but genetically heterogeneous according to CGH. MLSA-based further subdivision of cluster 2 was not reliable. However, the tree based on CGH results showed that the outbreak I isolates consisted of a separate branch in cluster 2. This clonality of outbreak I isolates is also supported by the similarity of their rpoB and hsp60 sequences. The results from the second clade were less congruent in the CGH tree. Cluster 5 isolates formed a distinct cluster in CGH and MLSA and were all identified as E. asburiae by rpoB and hsp60. Cluster 4 isolates were all genotyped as E. ludwigii and clustered together both in MLSA and CGH. Cluster 6 as well as cluster 7 were separate clusters in MLSA and CGH but all belonged to E. cloacae IV according to rpoB genotyping, whereas according to hsp60 genotyping all cluster 7 isolates were E. cloacae IV and all cluster 6 isolates were E. kobei. Isolates from cluster 3 were mixed on CGH analysis, indicating a heterogeneous group consisting of multiple species. The solitaire isolate belonged to E. asburiae according to rpoB but E. cloacae III according to hsp60. Furthermore, this isolate also clustered separately in the CGH tree and the maximum likelihood tree, indicating that this specific isolate may represent yet another species.

Except for a lack of congruence with cluster 3 the results of CGH and MLSA were identical. This indicates that MLSA results are also highly representative of genomic relationships between E. cloacae complex isolates. However, the findings with CGH or MLSA are inconsistent with those obtained using fragments of rpoB or hsp60 (Table S1).

Evolution of the Enterobacter cloacae complex

The mean pairwise nucleotide diversity (π) of the 50 E. cloacae complex isolates was higher within clade 2 than clade 1 (0.047+/−0.003 vs. 0.014+/−0.002), suggesting that isolates from clade 2 are older than isolates from clade 1. Furthermore, the number of alleles was relatively lower in clade 1 compared to clade 2 (25/30 vs 20/20), which supports this finding (Table 2).

Selective markers

MLSA and CGH are costly and time-consuming techniques, whereas effective patient therapy requires inexpensive and rapid results. Therefore, new approaches that rapidly identify isolates are in demand. Genetic markers specific for each species may facilitate quick and inexpensive tests. Specific markers were identified within the set of 1,150 genes that were variably present among the E. cloacae complex isolates. Thirteen DNA fragments were identified that could potentially serve as selective markers for the first clade, and 2 potential selective markers were detected for the second clade (sensitivity and specificity of 100%) (Table 4). Four DNA fragments were identified for cluster 1 (sensitivity of 100% and specificity of ≥96%) and four specific DNA fragments were identified for outbreak I isolates (sensitivity of ≥95% and specificity of ≥95%) (Table 4).

[Figure omitted. See PDF.]

Table 4. Comparative results from the specific DNA fragments.

https://doi.org/10.1371/journal.pone.0003018.t004

Discussion

The results of this study show that the CGH approach with a random library of 7 different strains (8 isolates) is a powerful method for studying new or poorly characterized species such as those that belong to the E. cloacae complex. Combined with MLSA, this technique can help to accurately assign isolates to species or subspecies. In this study, the CGH array had a high resolution because it was based on the presence or absence of 1,150 different genetic elements derived from different isolates, which were all evaluated on an individual basis. Therefore, CGH was able to determine genetic differences between different isolates and allowed unsupervised comparison of large sets of isolates. However, CGH had some limitations. One potential problem was the large number of mobile elements represented on the MGA. These mobile elements may bias genetic relationships between isolates. If these limitations are taken into account, the CGH approach can be used to determine the genetic relationship between (sub)species and, importantly, to identify specific markers for these (sub)species. Although, the E. cloacae complex is genetically heterogeneous, this study clearly shows that clinical isolates are evolutionary divided into two genetically distinct clades, where clade 1 appears to be younger than clade 2. Clade 1 isolates were mostly obtained from clinical settings with hospitalized patients, which is considered to be an environment with high specific (antibiotic) pressure, whereas the second clade contains a relatively large number of fecal isolates from patients without gastro-intestinal complaints.

Retrospectively, seven of the eight isolates used to construct the MGA belonged to the first clade. Therefore, the CGH results for the second clade should be interpreted with caution. After removal of the pQC sequences, which represented the most important mobile element present in the isolates, a similar clustering of isolates, with the exception of cluster 3, was obtained using either CGH or MLSA (Figure 2 and 4). The diversity of the second clade was larger when using MLSA than, but as mentioned previously, this may have been due to the limited number of isolates of this clade on the MGA. With CGH, specific genetic markers were identified that enabled differentiation between the two clades. In addition, three different genetic markers were identified that enabled specific identification of isolates belonging to cluster 1. Four different genetic markers were identified that allowed discrimination of outbreak I isolates from all tested isolates. MLSA was only performed on one outbreak isolate; therefore it is unknown whether all outbreak isolates cluster together after MLSA analysis. However, it is likely that such clustering occurs because of the clonality of the isolates. No specific genetic markers were detected for clusters of the second clade; this may have been due to the low number of isolates of clade 2 represented on the MGA. In the future, a small array could be constructed that contains specific markers for different (sub)species in order to correctly speciate isolates belonging to the E. cloacae complex. This new method to define and identify new species may be extended to other micro-organisms.

When combined, the MLSA and CGH methods have a synergetic effect. Where the CGH-array offers the advantage of being able to detect specific discriminating markers, MLSA is a standardized method whose results can be used for additional comparative studies. MLSA analysis divided the isolates into seven clusters. Hsp60 and rpoB genotyping both did not sufficiently discriminate between the subspecies of E. hormaechei. Because CGH and MLSA also did also not discriminate between E. hormaechei subspecies, it is concluded that the current nomenclature for subspecies of E. hormaechei is inadequate. Furthermore, identification based on sequencing of hsp60 and rpoB fragments did not always give identical outcomes. However, this was mostly based on differences in nomenclature between genotyping methods (e.g. E. cloacae IV / E. kobei and E. hormaechei / E. cloacae III). If differences in nomenclature are discarded, only 1.9% (3/158) of the isolates were misidentified based on sequencing of hsp60 or rpoB fragments. Discrepancies arose for E. asburiae, E. cloacae, and an E. hormaechei according to rpoB genotyping, as well as E. cloacae III, E. cloacae dissolvens, and an E. cloacae IX according to hsp60 genotyping (Table S1). For comparison, the Phoenix Automated Microbiology System (Becton Dickinson Biosciences) identified all isolates as E. cloacae. The present study demonstrates the superiority of genotypic identification over typical phenotypic identification for members of the E. cloacae complex. CGH and MLSA are clearly better methods for identifying isolates of the E. cloacae complex, but these methods are more expensive and labor-intensive. For diagnostic purposes, sequencing the 501 bp fragment of rpoB remains an acceptable alternative.

Correct speciation of isolates is clinically relevant. The current study showed that even when outbreak isolates were not included, isolates of the first clade were isolated nearly two times more often than isolates of the second clade. Furthermore, it is notable that more than half of the 13 fecal E. cloacae isolates from patients without gastro-intestinal complaints were found in the smaller second clade. Two of the DNA fragments specific for the first clade contained a partial operon for pyrroloquinoline quinone (pqq) biosythesis. Pqq is a co-factor of several dehydrogenases and transfers redox equivalents to the respiratory chain [25]. It has also been postulated that pqq reduces oxidative stress [26]. In addition, pqq stimulates bacterial growth [27]. A tempting but speculative hypothesis is that the isolates in the first clade are more pathogenic and isolates in the second clade are more commensal.

In conclusion, MLSA combined with CGH is a powerful method for studying species assignment within the E. cloacae complex. The complex is divided into two genetically distinct clades and consists of several more related and heterogenic (sub)species that must be redefined because of disagreements with known or proposed nomenclature. The CGH approach is a novel tool that could potentially be used to rapidly identify Enterobacteriaceae or even detect specific strains that threaten hospital populations.

Materials and Methods

Bacterial strains

One hundred twenty of the 158 isolates included in this study originated from the UMCU including 13 fecal isolates that were obtained from cultures taken from patients without gastro-intestinal complaints during their admission (Table S1). These 13 fecal isolates were considered commensal or transitional flora. These isolates were a subset of a large collection of isolates typed with Pulsed Field Gel Electrophoresis (PFGE) that represented the E. cloacae complex population in our hospital [28]. Twenty seven of the isolates belonged to PFGE genotype 1 (outbreak I), which is of great interest to us because this strain caused a nationwide outbreak and is now endemic in the UMCU [28]–[30]. Most isolates of PFGE genotype 1 contain a conjugative plasmid, pQC that carries the qnrA1, blaCTX-M-9, and aadB genes, which encode resistance to quinolones, extended-spectrum β-lactamases, and aminoglycosides, respectively, and thus provide this outbreak strain with a multidrug-resistant phenotype [30]. In addition, 21 isolates (10 outbreak I isolates) from 11 other Dutch hospitals were included. Also included were 12 isolates from 12 different European hospitals and 3 isolates from an industrial site. Finally, E. cloacae ATCC13047 and an ATCC 13047 conjugant with the pQC plasmid were included.

All isolates were previously identified as E. cloacae using an automated system and software (Phoenix Automated Microbiology System, Becton Dickinson Biosciences, Sparks, MD, U.S.). The Enterobacter species form a heterogenic group and its nomenclature is still undergoing alterations. Isolates identified with the Phoenix apparatus as E. cloacae were used as a starting point. When hsp60 sequences indicated that an isolate did not belong to the genus Enterobacter it was excluded from the study.

In addition, E. cloacae ATCC 1307 with and without pQC was used for determining the contribution of DNA fragments on the microarray that were derived from this plasmid. Conjugation was performed as previously described [30]. All isolates were cultured overnight at 37°C on tryptic soy agar plates with 5% sheep blood. DNA was extracted with a NucleoSpin Tissue kit (Macherey-Nagel Gmbh & Co. KG, Düren, Germany) according to the manufacturer's instructions. DNA quality and quantity were measured with an ND-1000 Spectrophotometer (Wilmington, DE, U.S.) according to the manufacturer's instructions. A shotgun library was created from eight isolates (seven different strains). These eight isolates contained five isolates from the UMCU, two of which represented the outbreak strain (PFGE genotype 1) harboring the R-plasmid pQC (+/−300 kb), three isolates with different genetic backgrounds based on PFGE (including one with pQC), and three isolates that originated from different European countries (Table 3) [28], [30]. No high copy plasmid DNA was present in the DNA of the eight isolates used for the shotgun library.

Detection of specific genes and genotyping

Target DNA for PCR was extracted by heating bacterial suspensions for 10 min at 95°C. Amplification of DNA fragments from blaCTX-M-9, qnrA1, irp1, hsp60, fusA, gyrB, leuS, pyrG, rplB and rpoB genes was performed with a Hotstart Taq MasterMix kit (Qiagen, Westburg b.v., Leusden, The Netherlands). Primers used, the size of the obtained products, relative primer positions, and annealing temperatures are described in Table S5. PCR products were purified with a Qiaquick PCR purification kit (Qiagen). Purified DNA was bidirectionally sequenced by BaseClear B.V. (Leiden, The Netherlands) according to their quick shot protocol (www.baseclear.com).

Analysis of sequencing and genotyping results

Editing and analysis of chromatogram traces was performed using BioEdit [31]. Sequences of seven gene fragments (hsp60, fusA, gyrB, leuS, pyrG, rplB and rpoB) were confirmed by at least two chromatograms (forward and reverse). The quality of the phylogenetic information was determined by plotting transition and transversion rates as a function of the genetic distance calculated with a K80 model using DAMBE software [32], [33]. Pairwise nucleotide diversity (π), the number of segregating sides (θ), and tests for selection (Tajima's D, and Fu's F test) were calculated using DNAsp, version 4.20.2 [21], [22], [34]. The results showed (see results) that rplB is not suitable for phylogenetic analysis; thus rplB results were discarded for MLSA. ClonalFrame was used to assess the clonal relationship between E. cloacae complex isolates because we noticed high sequence diversity and putative recombination events in the six remaining gene fragments of the tested population. (http://www2.warwick.ac.uk/fac/sci/statistics/staff/research/didelot/clonalframe/). ClonalFrame determines the genetic relationships of bacteria based on point mutations and homologous recombination [35]. The burn-in length was 50,000 iterations and Markov chain Monte Carlo iterations were set at 50,000. After every 100 iterations, a posterior sample was recorded. Finally, a majority-ruled consensus tree was generated from all posterior samples. Congruence of the obtained tree was tested by phylogenic analyses of the concatenated sequences of the six gene fragments. The tree was constructed using maximum likelihood statistics with the “dnaml” program in PHYLIP. Parameter settings were: transition/transversion ratio = 2 based on data presented in Figure 1 and the transversion and transition rates as a function of the genetic distances, empirical base frequency, one category of rates, constant rate variation, and unweighted sites because sequence fragments are not expected to have different mutation frequencies and no selective pressure is expected (see results section). Prior to tree construction, 50.000 bootstraps of the data were performed using seqboot from the PHYLIP software package. Subsequently, the output was used to construct trees with SplitsTree4 (http://www.splitstree.org) [36].

Microarray development

Equal amounts (10 µg) of chromosomal DNA from eight E. cloacae complex isolates (7 strains) were mixed to create a shotgun library as described by Borucki et al. (Table 3) [16], [37], [38]. For each isolate, equal amounts of genomic DNA were mixed, 10 µg of the pooled DNA were sonicated (Branson 250/450 Sonifier, 6 mm microtip, output intensity 1, Geneva, Switzerland) and fragments of approximately 1.2 kb were extracted from agarose gels (Qiaquick columns, Qiagen) and end-repaired (DNA Terminator End Repair Kit, Lucigen Corporation, Middleton, WI, U.S.). End-repaired fragments were ligated into the pSMART-HC-Kan vector (Clone-SMART, Lucigen). Subsequently, the recombinant plasmids were transformed into E. coli (ElectroMAX DH10B Cells, Invitrogen, Breda, The Netherlands) and plated on tryptone yeast plates with 30 µg/ml kanamycin. A total of 3,072 recombinant clones were arrayed into 96-well plates. Cloned DNA fragments were amplified by PCR using SMART primers (Lucigen) with 5′-C6 aminolinkers to facilitate cross-linking to the aldehyde-coated glass slides. PCR products were ethanol-purified and resuspended in 3×SSC (1×SSC: 0.15 M NaCl, 0.015 M sodium citrate, pH 7). All PCR products were printed on CSS silylated slides (European Biotech, Network, Dolembreux, Belgium) by an ESI three-axis DB-3 robot (ChipWriter Pro, Biorad, Hercules, CA, U.S.) at a controlled humidity of 55%. Slides were printed in four batches, after which they were blocked and denatured (15 sec dip in boiling water) following the manufacturer's instructions.

Labeling, hybridization and data acquisition

Labeling, hybridization and data acquisition were performed as previously described by Leavis et al. [16]. A Bioprime system (Invitrogen) was used to label 0.5 µg DNA. For normalization, equal amounts of the 8 library isolates (7 strains) were mixed as the reference pool and labeled with Cy3 dUTP. Tester isolates were labeled with Cy5 dUTP. To test reproducibility, fifteen isolates were tested in duplicate, two in triplicate and one quadruplicate. For all hybridizations, Cy5 and Cy3 probes were combined with 100 µg yeast tRNA, dried with a Speedvac (30 min, at high temperature), resuspended in 40 µl Easyhyb buffer (Roche, Almere, The Netherlands), and denatured for 2 min at 100°C. Printed slides were pre-hybridized in 0.45 µm-filtered pre-hybridization buffer [1% BSA, 5×SSC, and 0.1% sodium dodecyl sulfate (SDS)] at 42°C for 45 min with rotation, then washed twice with purified water (MilliQ, Millipore, Billerica, MA, U.S.), dried with N2, and pre-warmed at 42°C. The hybridization mix was then pipetted on the printed side of the slide, covered with a hybrislib, placed in hybridization chambers (Corning Life Sciences B.V. Schiphol-Rijk, The Netherlands), and incubated overnight at 42°C in a water bath. Slides were then thoroughly washed sequentially in (a) 1×SSC, 0.2%SDS for 10 sec at 37°C, (b) 0.5×SSC for 10 sec at 37°C, and (c) twice in 0.2×SSC for 10 min at room temperature. Slides were dried with N2 and scanned using a Scanarray Express 680013 Microarray Analysis System (Perkin Elmers Life Analytical Sciences Inc., Wellesley, MA, U.S.). Images were obtained and quantified with ImaGene 4.2 software (Biodiscovery, Marina del Rey, CA, U.S.).

Processing and analysis of MGA data

Inferior spots (empty spots or those when the Cy3 signal was less than 2 times the background), were excluded from normalization and data analysis [16]. Each slide was independently normalized to correct for individual differences. Cy5 minus background/Cy3 minus background ratios were calculated. Ratios were normalized by correcting for the overall signal intensities in both channels. Spots were selected for further analyses when 95% of the slides had a hybridization result for that spot. Next, the data were log2 transformed. The estimated probability of presence (EPP) of each insert was determined using a GACK-transformation [39]. (http://falkow.stanford.edu/whatwedo/software/software.html). With GACK-transformation, it is possible to dynamically choose cut-offs for grouping into present or divergent/absence of DNA-fragments based on the shape of the distribution. For hierarchical clustering and principal component analysis a graded output was selected. In our case the graded output generated a range of values from −0.5 to 0.5 in increments of 0.05, with −0.5 corresponding to 5% or less EPP and 0.5 corresponding to 95% or more EPP. The E. cloacae complex core genome was defined as genomic DNA fragments present in all isolates (GACK data>−0.5)[16]. For statistical analysis, the data were transformed to binary output, where 0 corresponded to 5% or less EPP, 1 to 95% or more EPP, and values in the 5%–95% interval were discarded. For parsimony analyze, the data were transformed into binary output, where 0 corresponded to <50% EPP, 1 to >50% EPP, and blanks were represented with a question mark. http://evolution.genetics.washington.edu/phylip.html [16], [40].

Complete linkage hierarchical clustering with Pearson correlation and principal component analysis (PCA) were performed and visualized with TIGR MeV version 3.1 software (http://www.tm4.org/mev.html). Data were also analyzed using Pars from the data package PHYLIP [40]. (http://evolution.genetics.washington.edu/phylip.html) Pars is a general parsimony program that performs the Wagner parsimony method with multiple states. Isolate 76 was chosen as the outgroup because this isolate appeared to be the most diverse after hierarchical clustering. A consensus tree was constructed from 50,000 trees. All analyses used all DNA fragments minus core and minus pQC-related DNA fragments. The output from the Pars analysis was used to construct trees using SplitsTree4 (http://www.splitstree.org) [36].

Validation of the array

For each strain included in the shotgun library, the number of positive hybridizations (EPP >95%) was determined and gene coverage of the CGH array was calculated using a formula previously described by Akopyants et al. [24]: CG = 1−(1−(T+I−2(RO/G)))N

where CG = gene coverage, T = gene size, I = insert size, RO = required overlap and G = genome size. The average gene size was estimated to be 922 bp, which is the average gene length in Enterobacter sp. 638. The average length of DNA fragments spotted on the array was 1202 bp, based on sequences of 425 DNA fragments (data not shown). In contrast to oligonucleotide-based hybridizations, only overlaps between spotted DNA and test DNA are necessary for a random shotgun library based on DNA-DNA hybridization. The assumption was made that 100 nucleotides would be necessary for a positive hybridization result [24]. Finally, the size of the genomes of the Enterobacter isolates used to construct the CGH-array were estimated to be 4.5 million bp based on the whole genome sequences of Enterobacter sp. 638 and Enterobacter sakazakii ATCC BAA-894 (4,518,712 and 4,368,373 bp, respectively). Next, the minimal percentage of coverage of the genome was defined as:where CM represents the minimal coverage of the core genome, N the number of DNA fragments on the MGA that passed quality control, R the redundancy (1/number of unique DNA fragments sequenced-1), I the DNA fragment size, and G the genome size [16]. To evaluate the robustness of hybridization, the presence or absence in the tested isolates of two DNA fragments located on two different bidirectionally sequenced DNA fragments was confirmed by PCR. The DNA fragments chosen were part of irp1 encoding HMWP1 and part of orf513 and blaCTX-M-9.

Determination and exclusion of R-plasmid pQC-related DNA fragments

Since this study focused on studying the genetic relationships between species of the E. cloacae complex, interference by plasmid DNA should be excluded. The large size of pQC (+/−300 kb) prohibits its separation from chromosomal DNA. However, DNA fragments derived from pQC influence clustering analyses. Therefore, DNA fragments from pQC were removed for analyses. The paired student-T test (p<0.01) was used for three pairs of isolates to determine which DNA fragments were derived from pQC. Pairs of isolates were defined as isolates obtained from one patient with highly similar PFGE patterns (in our case all outbreak strains) where one isolate harbored pQC but the other did not. To verify pQC specificity, all isolates were hierarchically clustered using the 118 DNA fragments potentially derived from pQC. Twelve DNA fragments that were considered not to be specific for pQC were discarded. Isolates that tested positive on qnrA1-specific PCR were considered to carry pQC, as this gene is located on the plasmid [28], [30].

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Isolates used in this study

https://doi.org/10.1371/journal.pone.0003018.s001

(0.05 MB XLS)

Table S2.

Accession numbers of sequenced gene fragments

https://doi.org/10.1371/journal.pone.0003018.s002

(0.04 MB XLS)

Table S3.

Flag filtered, normalized 2log transferred data; * replicate slide

https://doi.org/10.1371/journal.pone.0003018.s003

(6.85 MB XLS)

Table S4.

Sequences of the sequenced DNA fragments used in the study

https://doi.org/10.1371/journal.pone.0003018.s004

(0.06 MB XLS)

Table S5.

Primers used for amplification and sequencing.

https://doi.org/10.1371/journal.pone.0003018.s005

(0.02 MB XLS)

Acknowledgments

We thank Tjomme van der Bruggen for his comments and reading of the manuscript.

Author Contributions

Conceived and designed the experiments: AP MPMC FHJS RCM JV ACF. Performed the experiments: AP MPMC. Analyzed the data: AP MPMC ACF. Contributed reagents/materials/analysis tools: AP MPMC AD ACF. Wrote the paper: AP MPMC FHJS MALvH AD ACF. Contributed to hypothesis generation and overall study design: AP MPMC FHJS MALvH AD RCM JV ACF.

Citation: Paauw A, Caspers MPM, Schuren FHJ, Leverstein-van Hall MA, Delétoile A, Montijn RC, et al. (2008) Genomic Diversity within the Enterobacter cloacae Complex. PLoS ONE3(8): e3018. https://doi.org/10.1371/journal.pone.0003018

References

1. Brenner DJ, McWhorter AC, Knutson JK, Steigerwalt AG (1982) Escherichia vulneris: a new species of Enterobacteriaceae associated with human wounds. J Clin Microbiol 15: 1133–1140.DJ BrennerAC McWhorterJK KnutsonAG Steigerwalt1982Escherichia vulneris: a new species of Enterobacteriaceae associated with human wounds.J Clin Microbiol1511331140

2. Brenner DJ, McWhorter AC, Kai A, Steigerwalt AG, Farmer JJ III (1986) Enterobacter asburiae sp. nov., a new species found in clinical specimens, and reassignment of Erwinia dissolvens and Erwinia nimipressuralis to the genus Enterobacter as Enterobacter dissolvens comb. nov. and Enterobacter nimipressuralis comb. nov. J Clin Microbiol 23: 1114–1120.DJ BrennerAC McWhorterA. KaiAG SteigerwaltJJ Farmer III1986Enterobacter asburiae sp. nov., a new species found in clinical specimens, and reassignment of Erwinia dissolvens and Erwinia nimipressuralis to the genus Enterobacter as Enterobacter dissolvens comb. nov. and Enterobacter nimipressuralis comb. nov.J Clin Microbiol2311141120

3. Mehlen A, Goeldner M, Ried S, Stindl S, Ludwig W, et al. (2004) Development of a fast DNA-DNA hybridization method based on melting profiles in microplates. Syst Appl Microbiol 27: 689–695.A. MehlenM. GoeldnerS. RiedS. StindlW. Ludwig2004Development of a fast DNA-DNA hybridization method based on melting profiles in microplates.Syst Appl Microbiol27689695

4. Hoffmann H, Stindl S, Ludwig W, Stumpf A, Mehlen A, et al. (2005) Reassignment of Enterobacter dissolvens to Enterobacter cloacae as E. cloacae subspecies dissolvens comb. nov and emended description of Enterobacter asburiae and Enterobacter kobei. Systematic and Applied Microbiology 28: 196–205.H. HoffmannS. StindlW. LudwigA. StumpfA. Mehlen2005Reassignment of Enterobacter dissolvens to Enterobacter cloacae as E. cloacae subspecies dissolvens comb. nov and emended description of Enterobacter asburiae and Enterobacter kobei.Systematic and Applied Microbiology28196205

5. Sanders WE Jr, Sanders CC (1997) Enterobacter spp.: pathogens poised to flourish at the turn of the century. Clin Microbiol Rev 10: 220–241.WE Sanders JrCC Sanders1997Enterobacter spp.: pathogens poised to flourish at the turn of the century.Clin Microbiol Rev10220241

6. Jones RN (2003) Global epidemiology of antimicrobial resistance among community-acquired and nosocomial pathogens: a five-year summary from the SENTRY Antimicrobial Surveillance Program (1997–2001). Semin Respir Crit Care Med 24: 121–134.RN Jones2003Global epidemiology of antimicrobial resistance among community-acquired and nosocomial pathogens: a five-year summary from the SENTRY Antimicrobial Surveillance Program (1997–2001).Semin Respir Crit Care Med24121134

7. Streit JM, Jones RN, Sader HS, Fritsche TR (2004) Assessment of pathogen occurrences and resistance profiles among infected patients in the intensive care unit: report from the SENTRY Antimicrobial Surveillance Program (North America, 2001). Int J Antimicrob Agents 24: 111–118.JM StreitRN JonesHS SaderTR Fritsche2004Assessment of pathogen occurrences and resistance profiles among infected patients in the intensive care unit: report from the SENTRY Antimicrobial Surveillance Program (North America, 2001).Int J Antimicrob Agents24111118

8. Kämpfer P, Nienhüser A, Packroff G, Wernicke F, Mehling A, et al. (2007) Molecular identification of coliform bacteria isolated from drinking water reservoirs with traditional methods and the Colilert-18 system. Int J Hyg Environ Health 211: 374–384.P. KämpferA. NienhüserG. PackroffF. WernickeA. Mehling2007Molecular identification of coliform bacteria isolated from drinking water reservoirs with traditional methods and the Colilert-18 system.Int J Hyg Environ Health211374384

9. Hoffmann H, Roggenkamp A (2003) Population genetics of the nomenspecies Enterobacter cloacae. Appl Environ Microbiol 69: 5306–5318.H. HoffmannA. Roggenkamp2003Population genetics of the nomenspecies Enterobacter cloacae.Appl Environ Microbiol6953065318

10. Case RJ, Boucher Y, Dahllöf I, Holmström C, Doolittle WF, et al. (2007) Use of 16S rRNA and rpoB genes as molecular markers for microbial ecology studies. Appl Environ Microbiol 73: 278–288.RJ CaseY. BoucherI. DahllöfC. HolmströmWF Doolittle2007Use of 16S rRNA and rpoB genes as molecular markers for microbial ecology studies.Appl Environ Microbiol73278288

11. Dobrindt U, Agerer F, Michaelis K, Janka A, Buchrieser C, et al. (2003) Analysis of genome plasticity in pathogenic and commensal Escherichia coli isolates by use of DNA arrays. J Bacteriol 185: 1831–1840.U. DobrindtF. AgererK. MichaelisA. JankaC. Buchrieser2003Analysis of genome plasticity in pathogenic and commensal Escherichia coli isolates by use of DNA arrays.J Bacteriol18518311840

12. Fitzgerald JR, Sturdevant DE, Mackie SM, Gill SR, Musser JM (2001) Evolutionary genomics of Staphylococcus aureus: insights into the origin of methicillin-resistant strains and the toxic shock syndrome epidemic. Proc Natl Acad Sci U S A 98: 8821–8826.JR FitzgeraldDE SturdevantSM MackieSR GillJM Musser2001Evolutionary genomics of Staphylococcus aureus: insights into the origin of methicillin-resistant strains and the toxic shock syndrome epidemic.Proc Natl Acad Sci U S A9888218826

13. Fitzgerald JR, Reid SD, Ruotsalainen E, Tripp TJ, Liu M, et al. (2003) Genome diversification in Staphylococcus aureus: Molecular evolution of a highly variable chromosomal region encoding the Staphylococcal exotoxin-like family of proteins. Infect Immun 71: 2827–2838.JR FitzgeraldSD ReidE. RuotsalainenTJ TrippM. Liu2003Genome diversification in Staphylococcus aureus: Molecular evolution of a highly variable chromosomal region encoding the Staphylococcal exotoxin-like family of proteins.Infect Immun7128272838

14. Fukiya S, Mizoguchi H, Tobe T, Mori H (2004) Extensive genomic diversity in pathogenic Escherichia coli and Shigella strains revealed by comparative genomic hybridization microarray. J Bacteriol 186: 3911–3921.S. FukiyaH. MizoguchiT. TobeH. Mori2004Extensive genomic diversity in pathogenic Escherichia coli and Shigella strains revealed by comparative genomic hybridization microarray.J Bacteriol18639113921

15. Israel DA, Salama N, Krishna U, Rieger UM, Atherton JC, et al. (2001) Helicobacter pylori genetic diversity within the gastric niche of a single human host. Proc Natl Acad Sci U S A 98: 14625–14630.DA IsraelN. SalamaU. KrishnaUM RiegerJC Atherton2001Helicobacter pylori genetic diversity within the gastric niche of a single human host.Proc Natl Acad Sci U S A981462514630

16. Leavis HL, Willems RJL, van Wamel WJB, Schuren FHJ, et al. (2007) Insertion sequence-driven diversification creates a globally dispersed emerging multiresistant subspecies of E. faecium. PLoS Pathog 3: e7. doi:10.1371/journal.ppat.0030007.HL LeavisRJL WillemsWJB van WamelFHJ Schuren2007Insertion sequence-driven diversification creates a globally dispersed emerging multiresistant subspecies of E. faecium.PLoS Pathog3e7. doi:10.1371/journal.ppat.0030007

17. Lindsay JA, Moore CE, Day NP, Peacock SJ, Witney AA, et al. (2006) Microarrays reveal that each of the ten dominant lineages of Staphylococcus aureus has a unique combination of surface-associated and regulatory genes. J Bacteriol 188: 669–676.JA LindsayCE MooreNP DaySJ PeacockAA Witney2006Microarrays reveal that each of the ten dominant lineages of Staphylococcus aureus has a unique combination of surface-associated and regulatory genes.J Bacteriol188669676

18. Porwollik S, Wong RMY, McClelland M (2002) Evolutionary genomics of Salmonella: gene acquisitions revealed by microarray analysis. Proc Natl Acad Sci U S A 99: 8956–8961.S. PorwollikRMY WongM. McClelland2002Evolutionary genomics of Salmonella: gene acquisitions revealed by microarray analysis.Proc Natl Acad Sci U S A9989568961

19. Salama N, Guillemin K, McDaniel TK, Sherlock G, Tompkins L, et al. (2000) A whole-genome microarray reveals genetic diversity among Helicobacter pylori strains. Proc Natl Acad Sci U S A 97: 14668–14673.N. SalamaK. GuilleminTK McDanielG. SherlockL. Tompkins2000A whole-genome microarray reveals genetic diversity among Helicobacter pylori strains.Proc Natl Acad Sci U S A971466814673

20. Salerno A, Delétoile A, Lefevre M, Ciznar I, Krovacek K, et al. (2007) Recombining population structure of Plesiomonas shigelloides (Enterobacteriaceae) revealed by multilocus sequence typing. J Bacteriol 189: 7808–7818.A. SalernoA. DelétoileM. LefevreI. CiznarK. Krovacek2007Recombining population structure of Plesiomonas shigelloides (Enterobacteriaceae) revealed by multilocus sequence typing.J Bacteriol18978087818

21. Tajima F (1989) Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123: 585–595.F. Tajima1989Statistical method for testing the neutral mutation hypothesis by DNA polymorphism.Genetics123585595

22. Fu YX (1997) Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147: 915–925.YX Fu1997Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection.Genetics147915925

23. Feil EJ, Spratt BG (2001) Recombination and the population structures of bacterial pathogens. Annu Rev Microbiol 55: 561–590.EJ FeilBG Spratt2001Recombination and the population structures of bacterial pathogens.Annu Rev Microbiol55561590

24. Akopyants NS, Clifton SW, Martin J, Pape D, Wylie T, et al. (2001) A survey of the Leishmania major Friedlin strain V1 genome by shotgun sequencing: a resource for DNA microarrays and expression profiling. Mol Biochem Parasitol 113: 337–340.NS AkopyantsSW CliftonJ. MartinD. PapeT. Wylie2001A survey of the Leishmania major Friedlin strain V1 genome by shotgun sequencing: a resource for DNA microarrays and expression profiling.Mol Biochem Parasitol113337340

25. Velterop JS, Sellink E, Meulenberg JJM, David S, Bulder I, et al. (1995) Synthesis of pyrroloquinoline quinone in vivo and in vitro and detection of an intermediate in the biosynthetic pathway. J Bacteriol 177: 5088–5098.JS VelteropE. SellinkJJM MeulenbergS. DavidI. Bulder1995Synthesis of pyrroloquinoline quinone in vivo and in vitro and detection of an intermediate in the biosynthetic pathway.J Bacteriol17750885098

26. Misra HS, Khairnar NP, Barik A, Indira PK, Mohan H, et al. (2004) Pyrroloquinoline-quinone: a reactive oxygen species scavenger in bacteria. FEBS Lett 578: 26–30.HS MisraNP KhairnarA. BarikPK IndiraH. Mohan2004Pyrroloquinoline-quinone: a reactive oxygen species scavenger in bacteria.FEBS Lett5782630

27. Adamowicz M, Conway T, Nickerson KW (1991) Nutritional complementation of oxidative glucose metabolism in Escherichia coli via pyrroloquinoline quinone-dependent glucose dehydrogenase and the Entner-Doudoroff pathway. Appl Environ Microbiol 57: 2012–2015.M. AdamowiczT. ConwayKW Nickerson1991Nutritional complementation of oxidative glucose metabolism in Escherichia coli via pyrroloquinoline quinone-dependent glucose dehydrogenase and the Entner-Doudoroff pathway.Appl Environ Microbiol5720122015

28. Paauw A, Verhoef J, Fluit AC, Blok HEM, Hopmans TEM, et al. (2007) Failure to control an outbreak of qnrA1-positive multidrug-resistant Enterobacter cloacae infection despite adequate implementation of recommended infection control measures. J Clin Microbiol 45: 1420–1425.A. PaauwJ. VerhoefAC FluitHEM BlokTEM Hopmans2007Failure to control an outbreak of qnrA1-positive multidrug-resistant Enterobacter cloacae infection despite adequate implementation of recommended infection control measures.J Clin Microbiol4514201425

29. Leverstein-van Hall MA, Blok HEM, Paauw A, Fluit AC, Troelstra A, et al. (2006) Extensive hospital-wide spread of a multidrug-resistant Enterobacter cloacae clone, with late detection due to a variable antibiogram and frequent patient transfer. J Clin Microbiol 44: 518–524.MA Leverstein-van HallHEM BlokA. PaauwAC FluitA. Troelstra2006Extensive hospital-wide spread of a multidrug-resistant Enterobacter cloacae clone, with late detection due to a variable antibiogram and frequent patient transfer.J Clin Microbiol44518524

30. Paauw A, Fluit AC, Verhoef J, Leverstein-van Hall MA (2006) Enterobacter cloacae outbreak and emergence of quinolone resistance gene in Dutch hospital. Emerg Infect Dis 12: 807–812.A. PaauwAC FluitJ. VerhoefMA Leverstein-van Hall2006Enterobacter cloacae outbreak and emergence of quinolone resistance gene in Dutch hospital.Emerg Infect Dis12807812

31. Hall TA (1999) BioEdit:a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser 41: 95–98.TA Hall1999BioEdit:a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT.Nucl Acids Symp Ser419598

32. Wirth T, Morelli G, Kusecek B, van Belkum A, van der Schee C, et al. (2007) The rise and spread of a new pathogen: seroresistant Moraxella catarrhalis. Genome Res 17: 1647–1656.T. WirthG. MorelliB. KusecekA. van BelkumC. van der Schee2007The rise and spread of a new pathogen: seroresistant Moraxella catarrhalis.Genome Res1716471656

33. Xia X, Xie Z (2001) DAMBE: software package for data analysis in molecular biology and evolution. J Hered 92: 371–373.X. XiaZ. Xie2001DAMBE: software package for data analysis in molecular biology and evolution.J Hered92371373

34. Rozas J, Sánchez-DelBarrio JC, Messeguer X, Rozas R (2003) DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 19: 2496–2497.J. RozasJC Sánchez-DelBarrioX. MesseguerR. Rozas2003DnaSP, DNA polymorphism analyses by the coalescent and other methods.Bioinformatics1924962497

35. Didelot X, Falush D (2007) Inference of bacterial microevolution using multilocus sequence data. Genetics 175: 1251–1266.X. DidelotD. Falush2007Inference of bacterial microevolution using multilocus sequence data.Genetics17512511266

36. Huson DH, Bryant D (2006) Application of phylogenetic networks in evolutionary studies. Mol Biol Evol 23: 254–267.DH HusonD. Bryant2006Application of phylogenetic networks in evolutionary studies.Mol Biol Evol23254267

37. Borucki MK, Kim SH, Call DR, Smole SC, Pagotto F (2004) Selective discrimination of Listeria monocytogenes epidemic strains by a mixed-genome DNA microarray compared to discrimination by pulsed-field gel electrophoresis, ribotyping, and multilocus sequence typing. J Clin Microbiol 42: 5270–5276.MK BoruckiSH KimDR CallSC SmoleF. Pagotto2004Selective discrimination of Listeria monocytogenes epidemic strains by a mixed-genome DNA microarray compared to discrimination by pulsed-field gel electrophoresis, ribotyping, and multilocus sequence typing.J Clin Microbiol4252705276

38. Vlaminckx BJM, Schuren FHJ, Montijn RC, Caspers MPM, Fluit AC, et al. (2007) Determination of the relationship between group A streptococcal genome content, M type, and toxic shock syndrome by a mixed genome microarray. Infect Immun 75: 2603–2611.BJM VlaminckxFHJ SchurenRC MontijnMPM CaspersAC Fluit2007Determination of the relationship between group A streptococcal genome content, M type, and toxic shock syndrome by a mixed genome microarray.Infect Immun7526032611

39. Kim CC, Joyce EA, Chan K, Falkow S (2002) Improved analytical methods for microarray-based genome-composition analysis. Genome Biol 3: RESEARCH0065.CC KimEA JoyceK. ChanS. Falkow2002Improved analytical methods for microarray-based genome-composition analysis.Genome Biol3RESEARCH0065

40. Felsenstein J (2005) 2005. PHYLIP (Phylogeny Inference Package) version 3.6. J. Felsenstein20052005. PHYLIP (Phylogeny Inference Package) version 3.6.Distributed by the author. Department of Genome Sciences, University of Washington, Seattle. Cladistics. Distributed by the author. Department of Genome Sciences, University of Washington, Seattle. Cladistics.

41. Kosako Y, Tamura K, Sakazaki R, Miki K (1996) Enterobacter kobei sp nov, a new species of the family Enterobacteriaceae resembling Enterobacter cloacae. Current Microbiology 33: 261–265.Y. KosakoK. TamuraR. SakazakiK. Miki1996Enterobacter kobei sp nov, a new species of the family Enterobacteriaceae resembling Enterobacter cloacae.Current Microbiology33261265

42. Hoffmann H, Stindl S, Stumpf A, Mehlen A, Monget D, et al. (2005) Description of Enterobacter ludwigii sp. nov., a novel Enterobacter species of clinical relevance. Syst Appl Microbiol 28: 206–212.H. HoffmannS. StindlA. StumpfA. MehlenD. Monget2005Description of Enterobacter ludwigii sp. nov., a novel Enterobacter species of clinical relevance.Syst Appl Microbiol28206212

43. Hoffmann H, Stindl S, Ludwig W, Stumpf A, Mehlen A, et al. (2005) Enterobacter hormaechei subsp. oharae subsp. nov., E. hormaechei subsp. hormaechei comb. nov., and E. hormaechei subsp. steigerwaltii subsp. nov., three new subspecies of clinical importance. J Clin Microbiol 43: 3297–3303.H. HoffmannS. StindlW. LudwigA. StumpfA. Mehlen2005Enterobacter hormaechei subsp. oharae subsp. nov., E. hormaechei subsp. hormaechei comb. nov., and E. hormaechei subsp. steigerwaltii subsp. nov., three new subspecies of clinical importance.J Clin Microbiol4332973303