Genomic and Geographic Context for the Evolution of High-Risk Carbapenem-Resistant Enterobacter cloacae Complex Clones ST171 and ST78
ABSTRACT
Recent reports have established the escalating threat of carbapenem-resistant Enterobacter cloacae complex (CREC). Here, we demonstrate that CREC has evolved as a highly antibiotic-resistant rather than highly virulent nosocomial pathogen. Applying genomics and Bayesian phylogenetic analyses to a 7-year collection of CREC isolates from a northern Manhattan hospital system and to a large set of publicly available, geographically diverse genomes, we demonstrate clonal spread of a single clone, ST171. We estimate that two major clades of epidemic ST171 diverged prior to 1962, subsequently spreading in parallel from the Northeastern to the Mid-Atlantic and Midwestern United States and demonstrating links to international sites. Acquisition of carbapenem and fluoroquinolone resistance determinants by both clades preceded widespread use of these drugs in the mid-1980s, suggesting that antibiotic pressure contributed substantially to its spread. Despite a unique mobile repertoire, ST171 isolates showed decreased virulence in vitro. While a second clone, ST78, substantially contributed to the emergence of CREC, it encompasses diverse carbapenemase-harboring plasmids, including a potentially hypertransmissible IncN plasmid, also present in other sequence types. Rather than heightened virulence, CREC demonstrates lineage-specific, multifactorial adaptations to nosocomial environments coupled with a unique potential to acquire and disseminate carbapenem resistance genes. These findings indicate a need for robust surveillance efforts that are attentive to the potential for local and international spread of high-risk CREC clones.
IMPORTANCE Carbapenem-resistant Enterobacter cloacae complex (CREC) has emerged as a formidable nosocomial pathogen. While sporadic acquisition of plasmid-encoded carbapenemases has been implicated as a major driver of CREC, ST171 and ST78 clones demonstrate epidemic potential. However, a lack of reliable genomic references and rigorous statistical analyses has left many gaps in knowledge regarding the phylogenetic context and evolutionary pathways of successful CREC. Our reconstruction of recent ST171 and ST78 evolution represents a significant addition to current understanding of CREC and the directionality of its spread from the Eastern United States to the northern Midwestern United States with links to international collections. Our results indicate that the remarkable ability of E. cloacae to acquire and disseminate cross-class antibiotic resistance rather than virulence determinants, coupled with its ability to adapt under conditions of antibiotic pressure, likely led to the wide dissemination of CREC.
INTRODUCTION
Carbapenem-resistant Enterobacteriaceae (CRE) represent a substantial threat to modern health care, challenging our present antibiotic armamentarium and increasing mortality and health care costs, particularly among chronically ill and immunocompromised patients (1). In the United States, Klebsiella pneumoniae bacteria, particularly the widespread health care-associated ST258 clone, have accounted for the majority of CRE infections since their initial detection in the early 2000s (2). However, reports of recent studies from throughout the United States have documented the expanding distribution of carbapenem-resistant Enterobacter cloacae complex (CREC) (3–11), while the incidence of carbapenem-resistant K. pneumoniae resistance has remained stable or has declined in some areas (7, 10, 12). Overall, the genomic background of CREC is characterized by high clonal diversity (7, 8). However, accumulating evidence suggests the widespread distribution and epidemic potential of two high-risk CREC clones, ST171 and ST78 (4–7). In most studies, ST171 and ST78 isolates harbored the plasmid-encoded K. pneumoniae carbapenemase (KPC), while diverse resistance mechanisms were detected in other clonal backgrounds (7, 8). Taken together, these data suggest that the emergence of CREC has a complex history, driven in part by its ability to acquire and maintain blaKPC-harboring plasmids.
Whole-genome sequencing (WGS) has demonstrated the development of CREC clonal outbreaks superimposed on a background of diverse multidrug-resistant (MDR) lineages. A WGS study of MDR E. cloacae in the United Kingdom and Ireland revealed substantial heterogeneity among patterns of clonal spread and emergence of antibiotic resistance but included very few carbapenem-resistant isolates (13). WGS analysis of a localized CREC outbreak in Minnesota and North Dakota was consistent with a clonal outbreak of ST171 associated with an IncFIA plasmid harboring blaKPC-3 (6). More recently, Chavda et al. explored the phylogenetic structure of geographically diverse carbapenem-resistant and carbapenem-susceptible Enterobacter spp., demonstrating diverse mechanisms of resistance driven largely by horizontal transfer of blaKPC-harboring plasmids followed by clonal spread (8). However, to date, little has been known about the timeline of the recent evolution of CREC and its dominant clones in the context of the introduction of carbapenems in the mid-1980s. Additional gaps in knowledge, including the unique genetic features of ST171 and ST78 and their potential impact on the success of these two high-risk clones, could have important clinical and infection control implications. Here we aimed to elucidate the genomic epidemiology of CREC at a New York City hospital and to investigate the evolutionary relationships between ST171 and ST78 genomes in New York City and geographically diverse locations. Using whole-genome phylogenetic and phylogeographic analyses, we assessed the spatiotemporal spread and distribution of mobile genetic elements (MGEs) within ST171 and ST78 clonal sublineages. We also considered the potential contribution of resistance and virulence determinants to the recent emergence of CREC in the United States and worldwide.
RESULTS
Population structure of Enterobacter cloacae complex isolates at a New York City hospital.
To assess the diversity of 125 E. cloacae complex isolates collected at a northern Manhattan hospital system, we constructed a phylogenetic tree based on core genome single nucleotide polymorphisms (SNPs) (Fig. 1). Isolates were distributed across 30 different sequence types (STs), predominately ST171 (n = 60) and ST78 (n = 25). Most isolates belonged to hsp60 clusters previously shown to be pathogenic in humans, including cluster VI (ST171), cluster III (ST78), and cluster VIII (14, 15). The remaining isolates were highly diverse and belonged to six different hsp60 clusters, including several rarely associated with human infection. Within each cluster, E. cloacae complex isolates demonstrated variable phenotypic resistance to different antibiotic classes (Fig. 1), although ST171 and, to a lesser extent, ST78 demonstrated remarkable cross-class antibiotic resistance.
FIG 1
Emergence and regional spread of CREC ST171.
In an analysis of the evolutionary history of the ST171 clone, we mapped 144 ST171 sequences to the E. hormaechei 34978 reference genome (GenBank accession number CP012165) (8) and included 84 publically available genomes to provide temporal and geographic context (see Table S2 in the supplemental material). A total of 1,659 core genome SNPs were identified and used to construct a maximum likelihood phylogenetic tree (Fig. 2A). ST171 core genomes were separated by 0 to 403 SNPs, indicating that many isolates were closely related. Using Bayesian phylogenetic analysis of isolates demonstrating evidence of clocklike evolution (n = 106), we estimated a substitution rate of ~2.7 SNPs per year (95% highest posterior density [HPD], 2.5 to 3.0), corresponding to a rate of 6.0 × 10−7 substitutions per nucleotide site per year. This estimate is comparable to that corresponding to a recently described collection of MDR E. cloacae isolates from the United Kingdom and Ireland, for which an average of 1.5 × 10−7 substitutions per nucleotide site per year across the genome (0.5 to 3.0 SNPs per year) was reported (13). ST171 isolates clustered into two major clades (clades A and B) which we estimate diverged prior to 1962 on the basis of the derived substitution rate (Fig. 2B).
FIG 2
Clade A, comprising approximately 10% of ST171 genomes, consisted of isolates originating from diverse geographic sites, including New York City, Boston, Brazil, and the United Kingdom (Fig. 2A). Isolates belonging to clade A either possessed blaKPC-4 (n = 9) or were blaKPC negative (n = 2). blaKPC-4-harboring isolates formed a distinct branch with a common ancestor dating to 1993 (posterior probability [PP] = 0.999) (Fig. 2B), which predates the first detection of CRE in the United States, in 1996 (16). Clade B consisted of two large clusters of closely related isolates from northern Manhattan and the northern Midwest, including Minnesota and North Dakota, as well as isolates from the United Kingdom, Boston, New York City, New Jersey, Pennsylvania, and Michigan. While most ST171 CREC isolates from clade B harbored blaKPC-3 (n = 124), sporadic isolates possessed blaKPC-2 (n = 2) or were KPC negative (n = 4). According to results of analysis of the ancestral nodes, blaKPC-3-harboring isolates shared one recent common ancestor, which emerged prior to 2002 (PP = 1). Reconstruction of whole-genome phylogeny by the use of Gubbins software revealed that ST171 clade B was additionally defined by several small recombination events (see Fig. S1 in the supplemental material). A 70-kb recombination region encoded a total of 80 proteins, including several proteins involved in DNA replication and 7 different conjugal transfer proteins (Table S3). Within each clade, we also found evidence for the development of fluoroquinolone resistance accompanied by the acquisition and establishment of mutations in the gyrA (g.248C>T, g.260A>C) and parC (g.239G>T) genes (17).
Phylogeographic Bayesian analysis of U.S. isolates with defined collection locations and evidence for clocklike evolution (n = 107) suggests that both clade A and clade B originated in Boston and spread to other regions within the United States (PP = 0.89) (Fig. 2C). From Boston, clade A likely spread to New York City whereas clade B spread to the New York-New Jersey area (PP = 0.95) and subsequently to New York City (PP = 0.88) and the northern Midwest (PP = 0.93) between 2007 and 2008 (Fig. 2B). This was followed by local proliferation within the northern Manhattan and northern Midwestern sites with occasional cross-regional links.
ST78 is defined by distinct CREC and extended-spectrum-beta-lactamase (ESBL)-producing sublineages.
Due to the lack of publicly available ST78 reference genomes, we generated a de novo hybrid assembly for NR0276 using Illumina short reads scaffolded with nanopore reads. The resulting NR0276 genome consisted of 2 chromosomal contigs with an average coverage of 83×, in addition to three complete plasmids assembled into single contigs (Whole Genome Shotgun Project accession number PNXT00000000). We mapped sequencing reads from 25 northern Manhattan isolates and 20 published ST78 genomes from Boston to this new reference sequence in order to reconstruct the evolutionary history of ST78 (Table S2). We identified and used 1,235 SNPs in the core chromosome to build a maximum likelihood phylogenetic tree (Fig. 3A). Pairwise SNP distances between ST78 isolates ranged from 0 to 421. We observed two major clades as confirmed through Bayesian phylogenetic analysis of isolates with clocklike evolution (n = 27) (Fig. 3B). On the basis of an estimated substitution rate of ~4.3 SNPs per year (95% HPD, 3.7 to 5.0), these distinct clades diverged well before 1984 (PP = 1), although a lack of earlier isolates precludes improved resolution of the date of divergence (Fig. 3B). Notably, the branch lengths separating these two clades and the length of time since divergence could not be linked to a major recombination event (Fig. S2).
FIG 3
Clade A, consisting of approximately 10% of the ST78 isolates (n = 5), was composed exclusively of CREC isolates from northern Manhattan harboring blaKPC-3 (Fig. 3A). The second, larger clade, which included isolates from both New York and Boston, consisted of several distinct subclades (subclades B, C, and D, Fig. 3A) defined by beta-lactamase carriage. Subclade B was composed of New York isolates containing blaKPC-2 (n = 4), while subclade C (n = 15) was defined by blaKPC-4 carriage and included both New York and Boston isolates. Subclade D, consisting of carbapenem-susceptible isolates (n = 16) harboring the ESBL blaCTX-M-15 gene, was interspersed with three sporadic carbapenem-resistant isolates, of which two were carbapenemase negative and one demonstrated uptake of blaKPC-3. Unlike the isolates in clade A, isolates in subclades B, C, and D were resistant to fluoroquinolones and were found to have missense mutations in the gyrA (g.247_248TC>AT), parC (g.239G>T), and acrB (g.215G>A) genes (17).
blaKPC-3-harboring ST78 CREC clade A had a most recent common ancestor dated at 1985 as suggested through Bayesian reconstruction (PP = 0.999) (Fig. 3B). Low sample size precluded confidence in the dating of the most recent common ancestor of blaKPC-2-containing subclade B, while subclade C, defined by blaKPC-4 carriage, had an ancestral node dated at 1997 (PP = 1). blaCTX-M-15-containing subclade D had a most recent common ancestor which Bayesian reconstruction suggested dates to 2000 (PP = 1) and likely diverged from subclades B and C in 1984 (PP = 0.999).
Repertoire of mobile genetic elements (MGEs) of ST171 and ST78 genomes.
MGEs unique to ST171 and ST78 genomes encoded a variety of adaptive mechanisms potentially enabling their success within hospital environments. Of eight predicted MGEs in ST171 reference genome 34978, three were distributed throughout both clades (Fig. 2A). A 30-kb genomic island (GI) (GI3; Table S3) unique to ST171 comprised 31 coding sequences (CDS), including a putative type IV toxin-antitoxin (TA) gene pair encoding cytoskeleton-binding toxin (CbtA) and cytoskeleton bundling-enhancing factor A (CbeA), which are thought to participate in cellular stress responses (18). Five additional genomic islands present only in ST171 clade B encoded several putative virulence factors, including arsenic resistance genes from the ars arsenic resistance operon (GI1) (19) and a cluster of genes encoding copper-binding and resistance proteins and belonging to the cop operon (GI2) (20). The ST171 reference genome also harbored six different phage regions, five of which were found throughout both clades (Fig. 2A). The 70-kb HK639 bacteriophage (Phage3; Table S3) putatively encoded an alpha-hemolysin, although this phage was limited to a small subset of ST171 isolates from New York City within clade B and was absent from other CREC clonal backgrounds.
Five of 10 genomic islands identified in the NR0276 ST78 reference genome were present in all ST78 isolates. These genomic islands encoded a variety of potential virulence factors, including a gene cluster encoding S-fimbrial adhesin proteins (GI1) (21) and the flagellin protein FliC and associated flagellar export and regulatory proteins (GI5) (22). However, the other five NR0276 genomic islands were present only in blaKPC-3-harboring clade A. This included GI6, a large genomic island (80 kb) composed of 61 CDS encoding the type II TA system PasT/PasI, which has been shown to enable pathogen persistence in limited-nutrient settings and in the face of oxidative and nitrosative stress in Escherichia coli (23). Genes involved in silver and copper sequestration and resistance, including components of the sil, cus, and cop operons, were also identified in a genomic island (GI9) associated with this clade (20, 24). The ST78 reference genome contained two phage regions, HK639 (Phage1) and the widespread 40-kb mEp390 bacteriophage (Phage2), which was present in most ST78 genomes as well as ST171 genomes.
Plasmid-encoded antimicrobial resistance genes.
The presumed mechanism of carbapenem resistance in almost all CREC isolates regardless of clonal background was carriage of blaKPC-harboring plasmids (87%), although 1 CREC isolate harbored blaNDM-1 and 12 isolates lacked detectable carbapenemase genes. In order to determine the genetic context of blaKPC in ST171 and ST78 isolates, we generated de novo hybrid assemblies of plasmids from representative isolates belonging to the major CREC ST171 and ST78 clades (Table S5). However, in non-ST171/ST78 isolates, particularly those harboring blaKPC-2 (n = 9), blaKPC was associated with variable plasmid profiles and in most cases could not be assigned to plasmid backbones using short-read data.
In ST171 clade A isolates, blaKPC-4 was present on an IncHI2A/IncHI2 plasmid (Fig. 2A). De novo assembly identified a large (315-kb) representative plasmid, pNR3082, which harbored blaKPC-4 within Tn4401b, an additional beta-lactamase gene, and aminoglycoside and sulfonamide resistance genes. The majority of ST171 clade B isolates putatively harbored an IncFIA plasmid harboring blaKPC-3, consistent with clonal spread (Fig. 2A). Representative 141-kb IncFIA plasmid pNR3024 contained the blaKPC-3 gene within Tn4401d as well as two additional beta-lactamase genes and several aminoglycoside and sulfonamide resistance genes (Fig. 4A). The backbone structures of pNR3024 were highly similar to those of p34798, an IncFIA plasmid derived from the 34978 reference genome, and pBK30683 in K. pneumoniae, which was previously shown to be widely distributed in hospitals in New York and New Jersey (Fig. 4; see also Table S5) (25). pNR3024 differed from these plasmids by 80 SNPs distributed in methyltransferase genes and genes associated with hypothetical proteins, a 98-bp insertion in the DNA polymerase III theta subunit of pNR3024, and an inverted 7.5-kb region distal to blaTEM-1 harboring an aminoglycoside resistance gene (Fig. 4). IncFIA plasmid pMNCRE44, which was identified among ST171 isolates in the northern Midwest (6), also shared a similar plasmid backbone consisting of three collinear blocks with close identity that resulted in approximately 60% homology with pNR3024 (Fig. 4), suggesting recent diversification.
FIG 4
Although ST78 included fewer CREC isolates, we found evidence for multiple distinct blaKPC-encoding plasmids in this clonal lineage (Fig. 3A). Most clade A CREC isolates harbored an IncN plasmid represented by pNR0276, a 69-kb plasmid containing blaKPC-3 within a Tn4401b transposon (Fig. 4B; see also Table S5). pNR0276 was also detected in three isolates belonging to unique STs (ST93, ST133, and ST269) and hsp60 clusters, indicating a broad host range. pNR0276 closely resembled the IncN pKm38_N plasmid from a Klebsiella michiganensis strain harboring blaKPC-2 collected in 1997 (Fig. 4B) (26), differing by only 22 SNPs and a 210-bp transposase deletion. In subclades B and C, respectively, blaKPC-2 and blaKPC-4 were putatively located on distinct and yet related plasmids which shared an IncN replicon type and contained conserved backbone structures (Fig. 3A) (5, 8). The evidence for conserved elements in IncN plasmid harboring different blaKPC subtypes was further demonstrated by alignment of pNR0276 to pBK31551 and blaKPC-4-harboring IncN plasmid derived from a K. pneumoniae isolate from a New Jersey hospital (27), which revealed the presence of three large collinear blocks with a total length of approximately 50 kb demonstrating >95% identity (Fig. 4B). We also found evidence in our collection for horizontal transfer of a pNR3024 plasmid from ST171 to ST78, as two closely related ST78 isolates obtained from the same patient 1 month apart (patient A, Table S1) differed only by the presence of this IncFIA plasmid.
To determine the in vitro transmissibility of plasmids from distinct CREC backgrounds, we determined transformation efficiencies for IncFIA and IncN plasmids. IncFIA plasmids were derived from a representative isolate within the major northern Manhattan cluster of ST171 clade A (pNR0011) as well as from outlying branches (pNR3040 and pNR3041). We also selected IncN plasmids from ST78 clade A (pNR3055) and from different STs (pNR1247 and ST454). We transformed these plasmids into E. coli DH10B cells and observed a significant difference in mean transformation efficiencies between the IncN pNR1247 plasmid and the IncFIA pNR0011_1 plasmid (7.2 × 104 CFU/µl [standard deviation {SD}, 6.5 × 103] versus 1.2 × 103 [SD, 1.4 × 103]; P = 0.01) (Fig. S3). Intermediate values were observed for IncFIA plasmids pNR3040 and pNR3041 and IncN plasmid pNR3055.
CREC is a low-virulence pathogen.
To evaluate differences in virulence between CREC clones and between CREC and ESBL isolates, we employed a cytotoxicity assay using immortalized, bone marrow-derived murine macrophages. Macrophages exposed to cell-free bacterial supernatant derived from CREC isolates compared to ESBL isolates demonstrated a significant difference in the levels of cell killing (median versus maximum lactate dehydrogenase [LDH] release, 7.6% versus 21.0%, P = 0.04) (Fig. 5A). Compared to ST78 isolates, ST171 isolates showed reduced cytotoxicity, although the difference was not significant (6.5% versus 21.0%, P =0.052) (Fig. 5B). However, the subset of CREC ST171 isolates showed significantly lower toxicity than the ESBL ST78 isolates (6.0% versus 21.5%, P = 0.007). There were no significant differences in the levels of cell killing in comparisons of isolates from different ST171 and ST78 clades (overall P = 0.2). We also noted no significant differences in the levels of cytotoxicity between ST171 and non-ST171/ST78 isolates (6.5% versus 10.6%, P = 0.8), although the latter produced variable levels of cell killing (Fig. 5B). There were no significant differences in cytotoxicity among the isolates based on culture source or blaKPC subtype (Fig. 5C and D).
FIG 5
DISCUSSION
In this comprehensive collection of newly sequenced CREC isolates and published genomes, we found evidence for five putative mechanisms of CREC evolution and dissemination (summarized in Fig. 6). These included (i) independent acquisition of blaKPC in ST171 by at least two distinct clonal sublineages with subsequent clonal spread and regional outbreaks; (ii) multiple instances of blaKPC uptake into a dominant ESBL clone, ST78; (iii) distribution of an IncN plasmid among several different STs, indicating a possible hypertransmissible plasmid; (iv) uptake of a variety of carbapenemase gene-harboring plasmids into diverse STs; and (v) sporadic emergence of non-carbapenemase-producing isolates. Overall, the introduction of blaKPC into multiple E. cloacae complex clones, including several distinct ST171 and ST78 sublineages, suggests that a combination of species-specific and external factors, such as antibiotic pressure, has contributed to the widespread emergence of CREC.
FIG 6
Our study provided robust evidence for the spread of two distinct ST171 clades across the United States, including a successful blaKPC-3-harboring sublineage demonstrating multiple foci of rapid, regionalized clonal proliferation in the Northeastern and northern Midwestern United States. Within the two clades, we also identified links between isolates from different regional and international locations, highlighting the potential for transmission of high-risk clones to contribute to the global spread of CRE. Although results of in vitro experiments suggested that ST171 is a low-virulence organism relative to other E. cloacae clones, ST171 isolates demonstrated remarkable cross-class resistance to antibiotics used commonly to treat Gram-negative infections, with the exception of amikacin. Notably, the acquisition of carbapenem and fluoroquinolone resistance determinants prior to widespread use of carbapenems and fluoroquinolones suggests that ST171 was already present in hospital settings prior to the mid-1980s and may have proliferated in the setting of rising antibiotic pressure.
Conversely, the emergence of CREC ST78 was characterized by multiple instances of sporadic uptake of blaKPC-harboring plasmids without evidence of extensive clonal spread. We identified three different CREC subclades defined by different blaKPC subtypes and a large sublineage of isolates from northern Manhattan demonstrating update of the plasmid-mediated blaCTX-M-15 ESBL gene. Moreover, ST78 isolates harbored unique genetic factors that may make this clone particularly adept at thriving in nosocomial environments, such as metal-binding and resistance proteins and genes putatively enabling niche-specific colonization. Given that previous population analyses of MDR E. cloacae demonstrated that ST78 is a widespread, globally dominant ESBL clone (28, 29), this suggests that ST78 is a highly successful hospital-associated clone with a unique ability to accept plasmids harboring resistance genes.
In addition to high-risk CREC clones, our findings suggest that increased attention to plasmid-mediated transmission of blaKPC is also warranted. We found evidence for recent horizontal transfer of a hypertransmissible IncN plasmid harboring blaKPC, which was identified in multiple CREC clones and demonstrated relatively high in vitro transformation efficiency. Among previous studies, IncN family plasmids were frequently implicated in dissemination of blaKPC among multiple bacterial strains (30, 31). Interestingly, pNR0276 was found to be closely related to blaKPC-harboring IncN plasmids collected in the 1990s (26), further suggesting that it is a well-established plasmid backbone with the potential to enable multispecies spread of blaKPC. However, we also detected a wide diversity of CREC isolates harboring heterogeneous plasmids, particularly among blaKPC-harboring isolates. This may be driven by a growing hospital reservoir of blaKPC-harboring plasmids as well as intrinsic factors that enable rapid uptake and acquisition of novel plasmids in E. cloacae complex.
Our study had several notable limitations. Our sample size was limited to available clinical CREC isolates at our hospital system and to published CREC genomes and may have reflected oversampling of clonal outbreaks. In order to perform Bayesian phylogenetic and phylogeographic analyses, several assumptions were required, and inherent limitations of these techniques may have introduced uncertainty into our findings. The accuracy of the evolutionary timeline for both ST171 and ST78 would be improved by the availability of more-recent clinical isolates or of additional earlier isolates collected prior to the early 2000s. Although we demonstrated high node support at all primary branch points in our analyses, the availability of additional isolates may lead to improved resolution of both recent branches with low support and long branches which could not be resolved due to a lack of earlier isolates. Moreover, our phylogeographic analysis was limited by the availability of isolates from only a few health care facilities and, more broadly, by the inclusion of isolates from only the Northeastern and northern Midwestern United States. Ultimately, phylogeographic analysis will produce a tree rooted in one of the sampled locations. Therefore, despite our high statistical support for the idea of an origin of ST171 in Boston, the future inclusion of isolates from additional locations may either support or refute this finding. Lastly, our cytotoxicity model was limited to demonstrating cytolytic activity of bacterial supernatant and may not fully reflect the in vivo virulence profile of CREC clones. In vivo models and functional studies are needed to further explore the role of specific genomic determinants in shaping the distribution and success of ST171 and ST78.
In summary, here we elucidated diverse evolutionary pathways of high-risk and sporadic E. cloacae clones in a large data set of genomes from across the United States and international sites. The inclusion of long-read sequencing and de novo hybrid assemblies, including new genomic and plasmid reference sequences, minimized reference bias and enabled detailed analysis of lineage-specific MGE and resistance gene repertoires. Our results indicate that, rather than heightened virulence, it was likely the proclivity of E. cloacae to acquire and disseminate multidrug resistance determinants, coupled with its ability to adapt to nosocomial environments, that led to the wide dissemination of CREC. The clonal expansion of ST171 across the United States and its subsequent local proliferation point to the need to pay increased attention to this high-risk clone. However, ongoing uptake of blaKPC-harboring plasmids has the potential to generate new CREC sublineages in ST78 and other E. cloacae complex subtypes. This indicates a need for heightened surveillance efforts attentive to the potential for both local and international spread of high-risk clones and the emergence of new CREC sublineages.
MATERIALS AND METHODS
Selection of bacterial isolates.
We sequenced 125 MDR E. cloacae complex isolates collected between November 2009 and July 2016 at a tertiary care hospital in northern Manhattan (see Table S1 in the supplemental material). This collection encompassed all available CREC isolates, including 46 isolates for which clinical and genotyping data were described previously (7), and a subset of ESBL isolates belonging to ST78 or ST171 for comparison. We obtained clinical and microbiological data from a review of electronic health records and ascertained colonization status versus infection status using CDC National Healthcare Safety Network surveillance definitions (32). Antimicrobial susceptibilities were determined using Vitek2 or Etest (BioMérieux) as part of routine clinical care and interpreted using Clinical and Laboratory Standards Institute breakpoints (33). We identified an additional 104 publicly available CREC genomes for inclusion in the analysis (Table S2).
Whole-genome sequencing and hybrid assembly.
Genomic DNA was extracted using a DNeasy UltraClean microbial DNA isolation kit (Qiagen). Index-tagged whole-genome libraries were sequenced using a HiSeq 2500 system (Illumina) and a MiSeq system (Illumina). We also conducted long-read sequencing for representative ST171 and ST78 isolates. DNA library preparation was performed using a Rapid Barcoding sequencing kit (Oxford Nanopore Technologies). Sequencing was then performed using a MinION DNA sequencer (Oxford Nanopore Technologies), and preprocessing of long reads with Poretools was performed as previously described (34, 35). Hybrid de novo assembly performed with SPAdes v3.10.1 (36) was followed by contig reordering using progressiveMauve in Geneious v10.0.5 (37) and annotation using Prokka v1.12 (38) and BLAST as previously described (35). Mobile genetic elements (MGEs) and prophage regions were identified by IslandViewer 3 (39) and PHAST (40), respectively. Gubbins v2.3.1 was used to identify possible recombination sites (41).
Isolates were classified by cluster membership (in clusters I to XI) based on hsp60 housekeeping gene typing (14) as well as on multilocus sequence typing (MLST) (42), antimicrobial resistance gene profiling (43), and plasmid replicon typing (44) using SRST2 (45). To identify sublineage-specific MGEs, we used bowtie2 v2.3.1 (46) to align short reads against MGEs present in the ST171 and ST78 reference genomes, with no ambiguous matches permitted. Alignment lengths were tabulated using SAMtools v1.5 (47), and MGE presence was defined as alignment of ≥95% of MGE lengths. Putative plasmid presence was determined based on identification of replicon type and blaKPC subtype using BLAST, as well as mapping of Illumina reads to reference plasmid sequences in Geneious v8.1.7 (48). Plasmids demonstrating at least 99% identity over 90% of the sequence length, including the replicon regions and the entire Tn4401 transposon (49), were classified as closely related.
Phylogenetic analysis.
We constructed phylogenetic trees by mapping individual genomes to ST171 and ST78 reference genomes to identify core genome single nucleotide polymorphisms (SNPs) using Snippy v3 (https://github.com/tseemann/snippy). Maximum likelihood phylogenetic analysis based on core genome concatenated SNPs was performed using RAxML v8.0.0 (50). Support for nodes was assessed using 100 rapid bootstrap inferences, and a final tree was selected through a maximum likelihood search under Gamma (50). Phylogenetic trees were visualized using iTOL v3 (51).
We performed Bayesian analysis to estimate the time to the most recent ancestor for ST78 and ST171. After all outliers were removed, a subset of ST171 and ST78 isolates displaying evidence of clock-like evolution as determined through TempEst (52) were used to estimate the substitution rate, time of emergence, and divergence dates within each clonal lineage using BEAST2 v2.4.7 (53–55). Using the HKY substitution Gamma site model with a relaxed lognormal molecular clock and coalescent Bayesian skyline prior, BEAST2 was run for 100 million generations, with sampling every 1,000 states and with isolate collection dates defined as available. Tracer v1.6 was used to inspect chain convergence and ensure an effective sample size (ESS) of ≥100 for all key estimated parameters (56). BEAST2’s TreeAnnotator v2.4.7 was used to identify the maximum clade credibility tree after 10% burn-in.
We also used Bayesian analysis to estimate the origins and phylogeographic spread of ST171 within the United States, with tip dates defined as described above and location assigned to NYC, New Jersey, Boston, or Minnesota/North Dakota. The BEAUti MultiTypeTree template was used to prepare the BEAST2 input using "location" as a discrete trait (57). The HKY substitution model and relaxed lognormal molecular clock were used, with the clock rate set to the ST171 BEAST analysis estimate (2e−3). The MultiTypeTree structured coalescent prior was used with lognormal distributions of kappa, population size, and rate matrix priors. The phylogeographic BEAST analysis was run for 400 million generations, with sampling every 10,000 states, and was repeated three times to ensure convergence of chains and sufficient ESS as calculated. LogCombiner v2.4.7 and Tree Annotator v2.4.7 were used to obtain a clade tree with maximum credibility after combining trees from all three runs with 10% burn-in. FigTree v1.4.3 was used for visualization of all BEAST trees.
Microbiological characterization.
The cytotoxicity of cell-free culture supernatant obtained from a randomly selected subset of previously sequenced E. cloacae isolates (n = 90) was assessed in vitro using an immortalized murine bone marrow-derived macrophage (iBMDM) cell line. Once iBMDM cells were grown to confluence in a 100-mm-diameter petri dish, they were subjected to trypsinization and seeded the day before intoxication at 50,000 cells/well. Cultures grown overnight in Trypticase soy broth at 37°C with shaking at 2,469 × g were diluted to 1:100 and incubated for an additional 8 h. Bacterial supernatant was then collected by centrifugation and incubated with murine macrophages for 2 h at 37°C in 96-well flat-bottom plates. Cell viability was assessed by measuring release of lactate dehydrogenase (LDH), corrected for maximum LDH release by addition of Triton X to a control well. Cell viability assays were performed on 10 different colonies on 4 different days and pooled for analysis. Median differences in percentages of the maximum LDH release by Triton X were compared between isolate groups using the Mann-Whitney test or the Kruskal-Wallis nonparametric test with Dunn’s correction for multiple comparisons as appropriate, and data were visualized using GraphPad Prism v7.02 (GraphPad Software, Inc.).
IncFIA and IncN plasmids, including plasmids derived from major CREC ST171 and ST78 sublineages, were isolated and purified using a plasmid minikit (Qiagen) and were electroporated into competent E. coli DH10B cells (Invitrogen) using a Gene Pulser Xcell system (Bio-Rad Laboratories). Transformants were selected using LB agar plates containing 100 µg/ml carbenicillin. Gel electrophoresis was used to confirm plasmid sizes, and the presence of blaKPC was confirmed by PCR using established primers (58). Transformation experiments were performed in triplicate, and efficiencies were calculated as the number of transformants per microgram of competent cells.
Accession number(s).
Sequence data were deposited in the NCBI sequence read archive (SRA) (https://www.ncbi.nlm.nih.gov/sra) under accession numbers SRP099597 and SRP126514 (BioProjects PRJNA374677 and PRJNA421733). Individual BioSample numbers for each isolate are provided in Table S1. The NR0276, NR3024, and NR3082 hybrid assembly whole-genome shotgun projects have been deposited at DDBJ/ENA/GenBank under accession numbers PNXT00000000, PNXS00000000, and PNXR00000000, respectively. The versions described in this paper are PNXT01000000, PNXS01000000, and PNXR01000000. Novel housekeeping gene allelic variants were reported to the E. cloacae MLST database curator (http://pubmlst.org/ecloacae/).
ACKNOWLEDGMENTS
We thank the Columbia University Medical Center Clinical Microbiology Laboratory staff for their assistance in collecting the MDR isolates included in this study. We also thank Alexis Sommerfield, Rita Chan, and William Sause for assistance with streaking and growing the strains used for the cytotoxicity assay and Jon Kagan for providing the iBMDM cell line.
This work was supported by grants from the National Institutes of Health National Institute of Allergy and Infectious Diseases (R01 AI116939 [A.-C.U.]) and the Irving Institute of Columbia University Medical Center (A.G.-S., A.-C.U.). The laboratory of V.J.T. was supported by funds from NYU School of Medicine’s Office of Therapeutic Alliance (V.J.T.). The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Supplemental Material
File (mbo003183911sf1.pdf)
- Download
- 14.90 KB
File (mbo003183911sf2.pdf)
- Download
- 5.71 KB
File (mbo003183911sf3.pdf)
- Download
- 35.30 KB
File (mbo003183911st1.xlsx)
- Download
- 28.95 KB
File (mbo003183911st2.xlsx)
- Download
- 13.17 KB
File (mbo003183911st3.xlsx)
- Download
- 11.47 KB
File (mbo003183911st4.xlsx)
- Download
- 10.07 KB
File (mbo003183911st5.xlsx)
- Download
- 11.71 KB
REFERENCES
1.
CDC. 2013. Antibiotic resistance threats in the United States. http://www.cdc.gov/drugresistance/threat-report-2013/. Accessed 1 September 2016.
2.
CDC. 2013. Vital signs: carbapenem-resistant Enterobacteriaceae. MMWR Morb Mortal Wkly Rep 62:165–170.
3.
Kiedrowski LM, Guerrero DM, Perez F, Viau RA, Rojas LJ, Mojica MF, Rudin SD, Hujer AM, Marshall SH, Bonomo RA. 2014. Carbapenem-resistant Enterobacter cloacae isolates producing KPC-3, North Dakota, USA. Emerg Infect Dis 20:1583–1585.
4.
Ahn C, Syed A, Hu F, O’Hara JA, Rivera JI, Doi Y. 2014. Microbiological features of KPC-producing Enterobacter isolates identified in a U.S. hospital system. Diagn Microbiol Infect Dis 80:154–158.
5.
Pecora ND, Li N, Allard M, Li C, Albano E, Delaney M, Dubois A, Onderdonk AB, Bry L. 2015. Genomically informed surveillance for carbapenem-resistant enterobacteriaceae in a health care system. MBio 6:e01030-15.
6.
Hargreaves ML, Shaw KM, Dobbins G, Snippes Vagnone PM, Harper JE, Boxrud D, Lynfield R, Aziz M, Price LB, Silverstein KA, Danzeisen JL, Youmans B, Case K, Sreevatsan S, Johnson TJ. 2015. Clonal dissemination of Enterobacter cloacae harboring blaKPC-3 in the upper midwestern United States. Antimicrob Agents Chemother 59:7723–7734.
7.
Gomez-Simmonds A, Hu Y, Sullivan SB, Wang Z, Whittier S, Uhlemann AC. 2016. Evidence from a New York City hospital of rising incidence of genetically diverse carbapenem-resistant Enterobacter cloacae and dominance of ST171, 2007–14. J Antimicrob Chemother 71:2351–2353.
8.
Chavda KD, Chen L, Fouts DE, Sutton G, Brinkac L, Jenkins SG, Bonomo RA, Adams MD, Kreiswirth BN. 2016. Comprehensive genome analysis of carbapenemase-producing Enterobacter spp.: new insights into phylogeny, population structure, and resistance mechanisms. MBio 7:e02093-16.
9.
Sheppard AE, Stoesser N, Wilson DJ, Sebra R, Kasarskis A, Anson LW, Giess A, Pankhurst LJ, Vaughan A, Grim CJ, Cox HL, Yeh AJ, Modernising Medical Microbiology (MMM) Informatics Group, Sifri CD, Walker AS, Peto TE, Crook DW, Mathers AJ. 2016. Nested Russian doll-like genetic mobility drives rapid dissemination of the carbapenem resistance gene blaKPC. Antimicrob Agents Chemother 60:3767–3778.
10.
Wilson BM, El Chakhtoura NG, Patel S, Saade E, Donskey CJ, Bonomo RA, Perez F. 2017. Carbapenem-resistant Enterobacter cloacae in patients from the US Veterans Health Administration, 2006–2015. Emerg Infect Dis 23:878–880.
11.
Kanamori H, Parobek CM, Juliano JJ, van Duin D, Cairns BA, Weber DJ, Rutala WA. 2017. A prolonged outbreak of KPC-3-producing Enterobacter cloacae and Klebsiella pneumoniae driven by multiple mechanisms of resistance transmission at a large academic burn center. Antimicrob Agents Chemother 61:e00912-17.
12.
Park SO, Liu J, Furuya EY, Larson EL. 2016. Carbapenem-resistant Klebsiella pneumoniae infection in three New York City hospitals trended downwards from 2006 to 2014. Open Forum Infect Dis 3:ofw222.
13.
Moradigaravand D, Reuter S, Martin V, Peacock SJ, Parkhill J. 2016. The dissemination of multidrug-resistant Enterobacter cloacae throughout the UK and Ireland. Nat Microbiol 1:16173.
14.
Hoffmann H, Roggenkamp A. 2003. Population genetics of the nomenspecies Enterobacter cloacae. Appl Environ Microbiol 69:5306–5318.
15.
Kremer A, Hoffmann H. 2012. Prevalences of the Enterobacter cloacae complex and its phylogenetic derivatives in the nosocomial environment. Eur J Clin Microbiol Infect Dis 31:2951–2955.
16.
Yigit H, Queenan AM, Anderson GJ, Domenech-Sanchez A, Biddle JW, Steward CD, Alberti S, Bush K, Tenover FC. 2001. Novel carbapenem-hydrolyzing beta-lactamase, KPC-1, from a carbapenem-resistant strain of Klebsiella pneumoniae. Antimicrob Agents Chemother 45:1151–1161.
17.
Ruiz J. 2003. Mechanisms of resistance to quinolones: target alterations, decreased accumulation and DNA gyrase protection. J Antimicrob Chemother 51:1109–1117.
18.
Masuda H, Tan Q, Awano N, Wu KP, Inouye M. 2012. YeeU enhances the bundling of cytoskeletal polymers of MreB and FtsZ, antagonizing the CbtA (YeeV) toxicity in Escherichia coli. Mol Microbiol 84:979–989.
19.
Andres J, Bertin PN. 2016. The microbial genomics of arsenic. FEMS Microbiol Rev 40:299–322.
20.
Bondarczuk K, Piotrowska-Seget Z. 2013. Molecular basis of active copper resistance mechanisms in Gram-negative bacteria. Cell Biol Toxicol 29:397–405.
21.
Girardeau JP, Bertin Y, Callebaut I. 2000. Conserved structural features in class I major fimbrial subunits (Pilin) in gram-negative bacteria. Molecular basis of classification in seven subfamilies and identification of intrasubfamily sequence signature motifs which might Be implicated in quaternary structure. J Mol Evol 50:424–442.
22.
Soutourina OA, Bertin PN. 2003. Regulation cascade of flagellar expression in Gram-negative bacteria. FEMS Microbiol Rev 27:505–523.
23.
Norton JP, Mulvey MA. 2012. Toxin-antitoxin systems are important for niche-specific colonization and stress resistance of uropathogenic Escherichia coli. PLoS Pathog 8:e1002954.
24.
Randall CP, Gupta A, Jackson N, Busse D, O’Neill AJ. 2015. Silver resistance in Gram-negative bacteria: a dissection of endogenous and exogenous mechanisms. J Antimicrob Chemother 70:1037–1046.
25.
Chen L, Chavda KD, Melano RG, Hong T, Rojtman AD, Jacobs MR, Bonomo RA, Kreiswirth BN. 2014. Molecular survey of the dissemination of two blaKPC-harboring IncFIA plasmids in New Jersey and New York Hospitals. Antimicrob Agents Chemother 58:2289–2294.
26.
Eilertson B, Chen L, Chavda KD, Kreiswirth BN. 2017. Genomic characterization of two KPC-producing Klebsiella isolates collected in 1997 in New York City. Antimicrob Agents Chemother 61:e02458-16.
27.
Chen L, Chavda KD, Fraimow HS, Mediavilla JR, Melano RG, Jacobs MR, Bonomo RA, Kreiswirth BN. 2013. Complete nucleotide sequences of blaKPC-4- and blaKPC-5-harboring IncN and IncX plasmids from Klebsiella pneumoniae strains isolated in New Jersey. Antimicrob Agents Chemother 57:269–276.
28.
Girlich D, Poirel L, Nordmann P. 2015. Clonal distribution of multidrug-resistant Enterobacter cloacae. Diagn Microbiol Infect Dis 81:264–268.
29.
Izdebski R, Baraniak A, Herda M, Fiett J, Bonten MJ, Carmeli Y, Goossens H, Hryniewicz W, Brun-Buisson C, Gniadkowski M; MOSAR WP2, WP3 and WP5 Study Groups. 2015. MLST reveals potentially high-risk international clones of Enterobacter cloacae. J Antimicrob Chemother 70:48–56.
30.
Conlan S, Thomas PJ, Deming C, Park M, Lau AF, Dekker JP, Snitkin ES, Clark TA, Luong K, Song Y, Tsai YC, Boitano M, Dayal J, Brooks SY, Schmidt B, Young AC, Thomas JW, Bouffard GG, Blakesley RW, NISC Comparative Sequencing Program, Mullikin JC, Korlach J, Henderson DK, Frank KM, Palmore TN, Segre JA. 2014. Single-molecule sequencing to track plasmid diversity of hospital-associated carbapenemase-producing Enterobacteriaceae. Sci Transl Med 6:254ra126.
31.
Adler A, Khabra E, Paikin S, Carmeli Y. 2016. Dissemination of the blaKPC gene by clonal spread and horizontal gene transfer: comparative study of incidence and molecular mechanisms. J Antimicrob Chemother 71:2143–2146.
32.
CDC/NHSN. 2016. CDC/NHSN surveillance definitions for specific type of infections. In National Healthcare Safety Network (NHSN) patient safety component manual. http://www.cdc.gov/nhsn/pdfs/pscmanual/pcsmanual_current.pdf. Accessed 15 September 2016.
33.
CLSI. 2015. Performance standards for antimicrobial susceptibility testing. Twenty-fifth informational supplement, M100-S25 edition. Clinical and Laboratory Standards Institute, Wayne, PA.
34.
Loman NJ, Quinlan AR. 2014. Poretools: a toolkit for analyzing nanopore sequence data. Bioinformatics 30:3399–3401.
35.
Giddins MJ, Macesic N, Annavajhala MK, Stump S, Khan S, McConville TH, Mehta M, Gomez-Simmonds A, Uhlemann AC. 2018. Successive emergence of ceftazidime-avibactam resistance through distinct genomic adaptations in blaKPC-2-harboring Klebsiella pneumoniae sequence type 307 isolates. Antimicrob Agents Chemother 62.
36.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, Pyshkin AV, Sirotkin AV, Vyahhi N, Tesler G, Alekseyev MA, Pevzner PA. 2012. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol 19:455–477.
37.
Darling AE, Mau B, Perna NT. 2010. ProgressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One 5:e11147.
38.
Seemann T. 2014. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30:2068–2069.
39.
Dhillon BK, Chiu TA, Laird MR, Langille MG, Brinkman FS. 2013. IslandViewer update: improved genomic island discovery and visualization. Nucleic Acids Res 41:W129–W132.
40.
Zhou Y, Liang Y, Lynch KH, Dennis JJ, Wishart DS. 2011. PHAST: a fast phage search tool. Nucleic Acids Res 39:W347–W352.
41.
Croucher NJ, Page AJ, Connor TR, Delaney AJ, Keane JA, Bentley SD, Parkhill J, Harris SR. 2015. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res 43:e15.
42.
Miyoshi-Akiyama T, Hayakawa K, Ohmagari N, Shimojima M, Kirikae T. 2013. Multilocus sequence typing (MLST) for characterization of Enterobacter cloacae. PLoS One 8:e66358.
43.
Gupta SK, Padmanabhan BR, Diene SM, Lopez-Rojas R, Kempf M, Landraud L, Rolain JM. 2014. ARG-ANNOT, a new bioinformatic tool to discover antibiotic resistance genes in bacterial genomes. Antimicrob Agents Chemother 58:212–220.
44.
Carattoli A, Zankari E, García-Fernández A, Voldby Larsen M, Lund O, Villa L, Møller Aarestrup F, Hasman H. 2014. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother 58:3895–3903.
45.
Inouye M, Dashnow H, Raven LA, Schultz MB, Pope BJ, Tomita T, Zobel J, Holt KE. 2014. SRST2: rapid genomic surveillance for public health and hospital microbiology labs. Genome Med 6:90.
46.
Langmead B, Salzberg SL. 2012. Fast gapped-read alignment with Bowtie 2. Nat Methods 9:357–359.
47.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R; 1000 Genome Project Data Processing Subgroup. 2009. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25:2078–2079.
48.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, Thierer T, Ashton B, Meintjes P, Drummond A. 2012. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28:1647–1649.
49.
Naas T, Cuzon G, Truong HV, Nordmann P. 2012. Role of ISKpn7 and deletions in blaKPC gene expression. Antimicrob Agents Chemother 56:4753–4759.
50.
Stamatakis A, Hoover P, Rougemont J. 2008. A rapid bootstrap algorithm for the RAxML Web servers. Syst Biol 57:758–771.
51.
Letunic I, Bork P. 2016. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res 44:W242–W245.
52.
Rambaut A, Lam TT, Max Carvalho L, Pybus OG. 2016. Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen). Virus Evol 2:vew007.
53.
Bouckaert R, Heled J, Kühnert D, Vaughan T, Wu CH, Xie D, Suchard MA, Rambaut A, Drummond AJ. 2014. BEAST 2: a software platform for Bayesian evolutionary analysis. PLoS Comput Biol 10:e1003537.
54.
Hasegawa M, Kishino H, Yano T. 1985. Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol 22:160–174.
55.
Drummond AJ, Ho SY, Phillips MJ, Rambaut A. 2006. Relaxed phylogenetics and dating with confidence. PLoS Biol 4:e88.
56.
Rambaut A, Suchard MA, Xie D, Drummond AJ. 2014. Tracer v1.6. http://tree.bio.ed.ac.uk/software/tracer/. Accessed 15 August 2017.
57.
Vaughan TG, Kühnert D, Popinga A, Welch D, Drummond AJ. 2014. Efficient Bayesian inference under the structured coalescent. Bioinformatics 30:2272–2279.
58.
Poirel L, Walsh TR, Cuvillier V, Nordmann P. 2011. Multiplex PCR for detection of acquired carbapenemase genes. Diagn Microbiol Infect Dis 70:119–123.
Information & Contributors
Information
Published In
mBio
Volume 9 • Number 3 • 5 July 2018
eLocator: 10.1128/mbio.00542-18
Editor: Robert A. Bonomo, Louis Stokes Veterans Affairs Medical Center
PubMed: 29844109
Copyright
© 2018 Gomez-Simmonds et al. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.
History
Received: 7 March 2018
Accepted: 7 May 2018
Published online: 29 May 2018
Keywords
Contributors
Editor
Robert A. Bonomo
Editor
Louis Stokes Veterans Affairs Medical Center
Metrics & Citations
Metrics
Note:
- For recently published articles, the TOTAL download count will appear as zero until a new month starts.
- There is a 3- to 4-day delay in article usage, so article usage will not appear immediately after publication.
- Citation counts come from the Crossref Cited by service.
Citations
Citation
2018.
Genomic and Geographic Context for the Evolution of High-Risk Carbapenem-Resistant Enterobacter cloacae Complex Clones ST171 and ST78. mBio
9:10.1128/mbio.00542-18.
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. For an editable text file, please select Medlars format which will download as a .txt file. Simply select your manager software from the list below and click Download.