Global evolutionary epidemiology, phylogeography and resistome dynamics of Citrobacter species, Enterobacter hormaechei, Klebsiella variicola, and Proteeae clones: A One Health analyses
=========================================================================================================================================================================================

* John Osei Sekyere
* Melese Abate Reta

## Abstract

**Background.** The global epidemiology and resistomes dynamics of multidrug-resistant *Citrobacter spp., Enterobacter hormaechei, Klebsiella variicola, morganella morganii, Proteus mirabilis* and *Providencia spp*. have not been described, despite their importance as emerging opportunistic clinical pathogens.

**Methods.** The genomes of the above-mentioned organisms were curated from PATRIC and NCBI and used for evolutionary epidemiology, phylogeography and resistome analyses. The phylogeny trees were drawn using RAXmL and edited with Figtree. The resistomes were curated from GenBank and the phylogeography was manually mapped.

**Results and conclusion.** *Mcr-9* and other *mcr* variants were highly prevalent in *E. hormaechei subsp*. and substantial in *C. freundii* whilst KPC, OXA-48, NDM, IMP, VIM, TEM, OXA and SHV were abundant in global *E. hormaechei subsp., Citrobacter freundii, P. mirabilis, P. stuartii* and *P. rettgeri* clones/clades. Species-specific ampCs were highly conserved in respective species whilst fluoroquinolones, aminoglycosides, macrolides, fosfomycin, chloramphenicol, tetracycline, sulphamethoxazole and trimethoprim resistance mechanisms were abundantly enriched in almost all clades of most of the species, making them extensively and pandrug resistant; *K. variicola, C. amalonaticus* and *C, koseri* had relatively few resistance genes. Vertical and horizontal resistome transmissions as well as local and international dissemination of strains evolving from common ancestors were observed, suggesting the anthroponotic, zoonotic, and food-/water-borne infectiousness of these pathogens. There is a global risk of pandrug resistant strains escalating local and international outbreaks of antibiotic-insensitive infections, initiating the dawn of a post-antibiotic era.

Keywords
*   Resistome
*   epidemiology
*   *Citrobacter*
*   *Klebsiella variicola*
*   *Morganella*
*   *Proteus*
*   *Providencia*
*   *Enterobacter hormaechei*

## Introduction

Antibiotic resistance is mainly disseminated via horizontal and vertical transmission through mobile genetic elements such as plasmids and transposons and through clonal and multiclonal expansion of same species 1–6. Conjugative plasmids have been implicated in the transmission of several resistance determinants within and across species, resulting in the presence of same or very similar resistomes in same and different species and clones 1,4,7–9. Thus, the emergence of plasmid-borne resistance genes is always a cause for concern as they help breach the species barrier and shuttle resistance genes (ARGs) from commensals and non-pathogenic bacteria to pathogenic ones or vice versa 3,4,10,11. Such has been the case with the emergence and rapid spread of extended-spectrum β-lactamases (ESBLs) viz., TEM, SHV, OXA and CTX-M, carbapenemases such as NDM, IMP, VIM, KPC and GES, the mobile colistin resistance gene *mcr-1* (to *mcr-10*) and recently, the mobile tigecycline resistance gene, *tet*(*X*) 12–17 Thus, such conjugative plasmids influence the genomic plasticity of several related and unrelated species and genera of bacteria 8,11,18–21.

Coupled with plasmid-borne dissemination of ARGs is the selection and expansion of specific drug-resistant clones 11,21, which quickly spread under antibiotic pressure to overpopulate their environments, facilitating their survival and subsequent spread to other environments 14,22,23. In cases where such clones harbour resistance plasmids, their expansion almost always lead to the concomitant replication and intra-clonal as well as inter-clonal spread of such plasmids 4,5,7,8,24,25 Thus, as such clones are disseminated through contact, food, water, farms, hospitals, and the environment, they carry with them these resistance plasmids to colonize new hosts and environments 6,26–28. It is thus not surprising to have same clones hosting the same plasmids, contain the same resistomes 4–6,21,28. This explains the presence of multi-drug resistance (MDR) in particular international clones such as *Klebsiella pneumoniae* ST258 and *E. coli* ST113 21,25,29,30.

Hence, tracing the phylogeography of clones and their associated resistance genes is highly critical in epidemiology and public health as it provides necessary data to contain the further spread of ARGs 26,28,31. In this work, the global evolutionary epidemiology and resistome dynamics of clinically important but relatively less isolated Enterobacteriaceae pathogens are described 8. It is notable that most of the recently emerged or novel resistance genes in bacteria have occurred in Enterobacteriaceae more than in any other family of bacteria, making Enterobacteriaceae particularly important medically 28,32,33. These include ESBL-, carbapenemase-, *mcr-* and *tet*(*X*)*-producing* producing Enterobacteriaceae, which have been classified by the WHO as high and critical priority pathogens due to their implication in high mortalities and morbidities 32,33. Although *Citrobacter spp., Enterobacter hormaechei, Klebsiella variicola, Morganella morganii, Proteus spp*. and *Providencia spp*. are not mostly reported as *Escherichia coli, Klebsiella pneumoniae*, and *Salmonella enterica*, they have been associated with multiple resistance and clinical fatalities 3,7,14,25,28,34 Due to the transferability of resistance plasmids between members of the Enterobacteriaceae, the global resistome epidemiology of these six genera is important as they could be eventually transferred to commonly isolated Enterobacteriaceae species 11,14,21,25,35,36.

## Results

### Included genomes

A total of 2,377 genomes from *C. freundii* (n=569 genomes), *C. koseri* (n=82 genomes), *C. amalonaticus* (n=35 genomes), *E. steigerwaltii* and *E oharae* (n=121 genomes), *E. xiangfangensis* (n= 90 genomes), *E. hormaechei* (n=563 genomes), *K. variicola* (n=574 genomes), *M. morganii* (n=59 genomes), *P. mirabilis* (n= 156 genomes), and *Providencia spp*. (n=128 genomes) were obtained from PATRIC and NCBI databases as at January 2020 and used for downstream analyses (Tables S1-S3). Included in *C. freundii, C. koseri, E. xiangfangesis, K. variicola* and *Providencia spp*. genomes were genomes of other Gram-negative bacterial species or genera that were initially classified within these respective species but later reclassified by NCBI’s ANI (average nucleotide identity) analysis. The genomes of these reclassified species were however maintained and included in the resistome analyses to serve as controls for comparison (Tables S1-S3).

The genomes were mainly isolated from human specimens, followed by animal (including food animals), plants (including food crops) and environmental specimens. The human and animal specimens used included urine, blood, stool, catheter tip, swabs etc. whilst the environmental specimens used included soils, hospital environments, water, wastewater, sinks etc. (Tables S1-S3). In all, these genomes were obtained from 67 countries globally, with the USA having the most genomes for all species: *C. freundii/spp*. (USA=233, China=34, France=14, Spain=11); *C. amalonaticus* (USA=17); *C. koseri* (USA=55); *E. steigerwaltii/oharae* (USA=40, Japan=15); *E. hormaechei* (USA=142, China=78, Japan=72, France=13, UK=10); *E. xiangfangensis* (USA=50, China=14, India=12); *K. variicola* (USA=289, Germany=50, China=24, Bangladesh=23); *M. morganii* (USA=16); *P. mirabilis* (USA=38, France=27, China=17); *Providencia spp*. (USA=61).

*C. freundii*, the only species among the species included in this analysis to have an MLST scheme, had 84 different clones or sequence types (STs). ST100 (n=51), ST22 (n=51), ST62 (n=18), ST11 (n=14), ST299 (n=11), ST8 (n=10), ST114 (n=8), and ST98 (n=8) were the commonest clones.

### Species epidemiology

#### Citrobacter species

Amongst the *C. freundii* genomes were other *Citrobacter spp*. such as *C. werkmanii, C. youngae, C. brakii, C. portucalensis* etc. *C. werkmanii* were isolated from humans, sprouts, and sinks whilst *C. brakii* were obtained from humans, vegetables, and hospital environments. Some C. *werkmanii* and *C. brakii* strains were found in USA, Germany, and India (Fig. S1A-B). *C. koseri* were mainly found from humans and from a mouse (clade B2) (Fig. S1A-B). A *Citrobacter spp*. (clade A) outbreak was observed in the US (Texas, Boston) and Mexico (Fig. S1A-B). *C. portucalensis* of the same clade (but different clones) were found in Nigeria (vegetables) and Brazil (turtle) as well as from effluents (UK), humans (China), chives, carrots and salad (Germany); these clustered in *C. freundii* clade B1. *C. freundii* strains of different STs and countries clustered together into clades, showing the wide distribution of strains from the same ancestor (or of close evolutionary distance); they were isolated from humans, plants, animals, and the environment (Fig. 1 & S1A-B).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F1)

Figure 1. 
Evolutionary epidemiology and resistome of global *Citrobacter freundii* isolates. *C. freundii* clustered into four main clades (A, B1, B2 and B3), highlighted with distinct colours. Clade B3 had the most resistome abundance and diversity. Strains from humans (blue labels), animals (red labels), plants (purple/mauve labels) and the environment (green labels) were found in the same clade/cluster. *Bla*CMY was conserved in these genomes.

*C. werkmannii, C. brakii* and *C. youngae* clustered within *C. freundii* clade A whilst *C. portucalensis* and *C. youngae*, were clustered in clade B; *C. koseri* clustered in clade B3 of *C. freundii* (Fig. 1 & S1A-B). *C. freundii* clades B2 and B3 had a richer resistome than clades A and B1, although blaCMY, *qnr, sul1/2*, aac(6’/3’)-like ARGs were common in all the clades; mph(A/E), *catAB, dfrAB, aadA* and blaTEM were common in clades B2 and B3. Important ARGs such as *mcr*, blaKPc, blaCTX-M, blaNDM, blaVIM, and blaOXA-48 were relatively rare and mainly found in clades B2 and B3 than in A and B1. Comparatively, other Enterobacteriaceae species (*E. coli, K. pneumoniae*, and *S. marcescens*) had strains with richer resistomes diversity than *C. freundii* clades A and B1, including chromosomal mutations and MDR efflux pumps. However, Clades B2 and B3 had comparable resistome diversity and abundance to the above-mentioned Enterobacteriaceae species, in which CMY was well-nigh absent except in *E. coli* (Fig. 1).

A local outbreak of *C. braakii* was observed in the UK (in humans), with closely related strains being isolated from biosolids (Canada), carrots (Germany), chicken (China) and beef (Canada) (Fig. S1A-B). *C. amalonaticus* was also found in humans, animals, plants and the environment (Fig. 2): *bla*CMY was not common in this species, but *bla*SED and *OqxA/B* were almost conserved in clades B1 and C, with clade B2 (all from France) having the richest resistome repertoire that included *bla*TEM, *bla*SHV, *bla*NDM and *mcr*-9.1 (Fig. 2). *C. koseri* had a relatively limited resistome diversity, with *bla*MAL-1/2 (clade B2), *bla*CKO (clade B3) and *fos*A7 (clade B2) being the commonest ARGs. Comparatively, the other Enterobacteriaceae species (e.g. *E. coli, K. pneumoniae, Enterobacter spp., Serratia spp*., and *Providencia spp.)* had richer resistome diversity than *C. koseri* (Fig. 3).

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F2)

Figure 2. 
Evolutionary epidemiology and resistome of global *Citrobacter amalonaticus* isolates. *C. amalonaticus* strains clustered into clades A (red highlight), B1 (green highlight), B2 (blue highlight) and C (mauve highlight); clade B2 had very rich resistome repertoire and were all from France, but the other clades had very few resistance genes. Strains from humans (blue labels), animals (red labels), plants (purple/mauve labels) and the environment (green labels) were found in the same clade/cluster. *Bla*SED and *oqxAB* were almost conserved in these genomes.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F3)

Figure 3. 
Evolutionary epidemiology and resistome of global *Citrobacter koseri* isolates. *C. koseri* strains clustered into clades A (grey highlight), B1 (light blue highlight), B2 (orange highlight) and B3 (mauve highlight). Strains from humans (blue labels) and animals (red labels) were found in the same clade/cluster. *Bla*CKO and *bla*MAL were almost conserved in these genomes.

#### Enterobacter species

The resistomes of *E. steigerwaltii* and *E. oharae* were extraordinarily rich and diverse, although clade A of *E. steigerwaltii* had lesser resistome abundance and diversity compared to clades B and C (Fig. 4). *E. oharae* only clustered in clade B of *E. steigerwaltii* and had most *bla*CTX-M, *bla*SHV, *bla*NDM, *bla*OXA-48 and *bla*VIM ARGs than *E. steigerwaltii* clade C. *bla*ACT was present in almost all *Enterobacter spp*. strains whilst other AmpCs such as *bla*LAP, *bla*SCO-1, *bla*SFO, *bla*TMB, *bla*DHA, and *bla*CARB were virtually absent. *E. cloacae* strains clustered closely with *E. steigerwaltii* clades B and C. overall, *E. steigerwaltii* and *E. oharae* strains were richly endowed with clinically important ARGs including carbapenemases, ESBLs, *aac(6’)-like, aac(3’)-like, aph(3’/3”)-like, aadA, dfrA, catA, fosA, oqxAB, qnrA/B/S/D, sul1/2*, and tet(A/B/D). *Mcr-9* genes were particularly abundant in clade B than in clades A and C of *E. steigerwaltii/cloacae* (Fig. 4).

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F4)

Figure 4. 
Evolutionary epidemiology and resistome of global *Enterobacter steigerwaltii/oharae* isolates. *E. steigerwaltii/oharae* isolates clustered into clades A, B and C, which were all obtained from humans from countries distributed across the globe. Clades B and C had very rich resistome repertoire; *bla*ACT was conserved in these genomes..

*E. cloacae* strains clustered with all *E. xiangfangensis* strains except clade A1, from which it was evolutionarily distant (Fig. 5). ARGs such as *bla* ACT, *catA/B, OqxA/B* and *fosA* were universally present in almost all the *E. xiangfangensis* and *E. cloacae* strains. Furthermore, *E. xiangfangensis* clade A1 had the least ARGs diversity and abundance whilst *E. cloacae* A1 and the remaining *E. xiangfangensis* (clades A2 and B) were richly endowed with multiple ARGs such as *bla*KPC, *mcr-9, bla*CTX-M, *bla*OXA, *bla*NDM, *bla*TEM, and genes mediating resistance to fluoroquinolones (aac(6’/3’)-like, *qnr*A/B/S etc.) and aminoglycosides (*armA, rmtC/G, aadA, aph*(4)-, *aph*(3’/3”)- and *aph*-(6’)-like) as well as *tet* (A/B/D), *dfrA*, and *arr* (particularly in clades A2 and B). Thus, the resistomes diversity and abundance of *E. xiangfangensis* and *E. cloacae* were comparable to that of *E. oharae* and *E. steigerwaltii* (Fig. 4-5).

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F5)

Figure 5. 
Evolutionary epidemiology and resistome of global *Enterobacter xiangfangensis* isolates. The *E. xiangfangensis* contained *E. cloacae* genomes and they clustered into clades A1, A2 and B, with clades A2 and B having rich and diverse resistome repertoire; these clades were distributed globally from humans (blue labels) and animals (red labels). *bla*ACT was conserved in these genomes.

An *E. hormaechei* outbreak was observed in Germany in 2017, evolving with the spread (clade B) (Fig. 6 & S2). *E. steigerwaltii* and *oharae* were mainly isolated from humans whilst *E. xiangfangensis* and *E. hormaechei* were from humans, plants (*E. xiangfangensis* from rice in India), animals and the environment. Strains of closely related (i.e. close evolutionary distance) *E. hormaechei* in clade B were from the US, South Africa, Colombia, China, Germany, Australia and Lebanon from humans, animals, and plants (Fig. 4-6). Other closely related *E. hormaechei* strains from humans, animals, plants, and the environment were found in different countries, showing a gradual evolution of strains emanating from a common ancestor and spreading across countries through different hosts. These observations were made in *E. hormaechei* clades B and C, and represented local and international outbreaks spanning UK, USA, China, Serbia, Japan etc. (Fig. 6 & S2).

![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F6.medium.gif)

[Figure 6.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F6)

Figure 6. 
Evolutionary epidemiology and resistome of global *Enterobacter hormaechei* isolates. The *E. hormaechei* isolates clustered into three main clades A, B and C (with distinct highlights) that contained strains distributed globally from humans (blue labels), and animals (red labels), plants (purple/mauve labels) and the environment (green labels). Clades B and C contained diverse and rich resistome repertoire. *bla*ACT was conserved in these genomes.

Moreover, *bla*ACT, *fosA*, OqxA/B and *cat*A/B were almost conserved in almost all the clades of *E. hormaechei*. As observed with the other *E. hormaechei* subsp. (*oharae, steigerwaltii* and *xiangfangensis*), *E. hormaechei* clades B and C were richly endowed with ARGs than some clade A strains (except clade A in Fig. S2A). Specifically, clades B and C strains were relatively enriched with *bla*TEM-1, *bla*OXA, *bla*CTX-M, *aph(3’/3”)-like, aac(3’)-like, aac(6’)-like, dfrA/B, mcr-9, arr, qnrA/B/S, sul-1/2*, and *tet*(A/B/C/D). Further, clade B was mostly rich in *bla*KPC and *bla*NDM whilst clade C was rich in *bla*NDM and *ble*. A substantial number of clade A strains also had *mcr-9* and *bla*TEM (Fig. 6 & S2).

*E. oharae* was evolutionarily closer to *E. xiangfangensis*, clustering on the same branches, with a relatively few *E. oharae* strains clustering with *E. steigerwaltii*. However, *E. steigerwaltii* was mostly distant from *E. oharae* and *E. xiangfangensis*, with a few *E. xiangfangensis* strains (from rice in India) clustering within *E. steigerwaltii* clades A and C (Fig. 2SB). Figure S2B summarises the ARGs in *E. hormaechei* and its subspecies: *bla*ACT, *fosA, OqxAB, qnrA/B/D/S*, and *catA/B* were conserved whilst *bla*OXA, ble, *bla*SHV, *bla*TEM, blaKPC, *bla*NDM, *bla*VIM, *mcr, tet, mph(A), sul-1/2, dfrA*, and fluoroquinolone and aminoglycoside ARGs were richly abundant, particularly in clusters IV, V and VII.

#### K. variicola

*K. variicola* strains were isolated mainly from humans, with some being from animals, plants, and the environment. On the individual branches/clades were closely related strains with very close evolutionary distance but disseminated across countries, suggesting international dissemination of that clade and showing little evolution during the spread from host to host, e.g., clades A3, A4, A5, A6, B1, B2 and C (Fig. 7 & S3). As well, local outbreaks of closely related strains in the USA, Bangladesh, Canada, and Germany were observed. A reassessment of the *K. variicola* tree with genomes of *K. pneumoniae* and *K. quasipneumoniae* largely confirmed the initial clustering of the *K. variicola* strains, with only a few rearrangements of some strains within different clades (Fig. S3A and S3B).

![Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F7.medium.gif)

[Figure 7.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F7)

Figure 7. 
Evolutionary epidemiology and resistome of global *Klebsiella variicola* isolates. The *K. variicola* strains clustered into nine clades viz., A1, A2, A3, A4, A5, A6, B1, B2 and C, which were highlighted with distinct colours and were isolated from countries around the globe. The clades contained strains distributed globally from humans (blue labels), animals (red labels), plants (purple/mauve labels) and the environment (green labels). Besides a few strains in clade B2, the other strains contained very few resistance genes. *bla*LEN was conserved in these genomes.

Notably, *K. variicola* had fewer resistome diversity and abundance than *Citrobacter spp., Enterobacter spp*., and *K. pneumoniae spp*. Conserved within the *K. variicola* genomes were *emrD, fosA, OqxAB* and *bla*LEN-2, whilst the other ARGs were sparse. Whilst *bla*SHV was conserved in *K. pneumoniae*, it was virtually absent in *K. variicola;* further, *bla*LEN was present in the latter but absent in the former. *K. variicola* clade B2 strains (on branch VIII), specifically those from Bangladesh, had richer and more diverse resistomes than the other *K. variicola* clades/clusters. Whilst *mcr* and *bla*NDM were virtually absent in these genomes, *bla*KPC occurred in substantial abundance. These suggest that *K. variicola* is least likely to be MDR and/or harbour ESBLs, carbapenemases, *mcr* and other clinically important ARGs compared to other Enterobacteriaceae species (Fig. 7 & S3).

#### M. morganii

There were three*M. morganii* clusters/clades viz., A, B and C, which were mainly from humans with a few in clades A and C being from animals (Fig. 8). Within each clade are closely related strains with very close evolutionary distance that were distributed across several countries; indeed, strains of the same clone were found in different countries, suggesting international dissemination of the same clones. The branching order of the trees within each clade shows the gradual evolution of the strains as they moved from host to host. *bla*DHA and *catA/B* genes were almost conserved in almost all the *M. morganii* genomes. Other highly abundant ARGs in the *M. morganii* strains were *sul-1/-2, tet*(A/B/D/Y), *bla*OXA-1, *ble, dfrA1, aadA, aac(3)-like* and *aph(3’/3”)-like*. Notably, clade B had more ARG diversity and abundance than clade C and A; clade A had the least diversity and abundance of ARGs (Fig. 8).

![Figure 8.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F8.medium.gif)

[Figure 8.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F8)

Figure 8. 
Evolutionary epidemiology and resistome of global *Morganella morganii* isolates. The *M. morganii* strains clustered into three clades, A (red highlight), B (light blue highlight) and C (yellow/gold highlight), containing isolates obtained globally from humans (blue labels) and animals (red labels).

#### P. mirabilis

*P. mirabilis* clustered into three major clades and included isolates from humans, animals, and the environment. Within the three clades were sub-clades consisting of closely related strains from the same as well as different countries, showing local and international outbreaks involving human and animal hosts, and in some cases environmental mediators (Fig. 9). Specifically, local outbreaks were seen in the USA (clades A3, B2 and B3), France (clade C2) and Japan (C3) whilst international dissemination of clades A2, A3, A4, B, C2 and C3 were observed. Clade C2 had the richest abundance of resistomes, with all the members having a uniform/conserved diversity of the same ARGs, except for *bla*CTX-M, *bla*OXA, *bla*TEM, and *Inu*(F/G). Clade C3 had the 2nd most abundant but more diverse ARGs than C2. Notably, clade B and its subclades had lesser ARGs than clades A and C. *CatA* and *tet*(A/B/D/Y) ARGs were virtually conserved in all the clades whilst *aadA, aac(3’)-like, aph(3’/3”)-like, aph(6’)-like, dfrA7, sat2*, and *sul-1/2/3* were substantially prevalent in all the clades. In particular, blaCARB2 and *florR* were highly conserved in clade C2; *florR* was however less abundant than *bla*CARB. Notably, *mcr* was almost absent except in a few strains in B3. As well, *bla*CTX-M, *bla*CMY, *bla*NDM, *bla*OXA-1, *bla*TEM and *ble-O/Sh* were mainly found in clades A and C3, with traces in B and C1. Thus, ESBLs, carbapenemases and *mcr* genes were relatively less common in *P. mirabilis* strains, compared to other Enterobacteiaceae (Fig. 9).

![Figure 9.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F9.medium.gif)

[Figure 9.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F9)

Figure 9. 
Evolutionary epidemiology and resistome of global *Proteus mirabilis* isolates. The *P. mirabilis* isolates clustered into 10 clades, A-A3, B1-B3, and C1-C3 (shown with different highlights), which contained diverse and abundant resistomes with conserved *catA* and *tet* genes. The clades contained strains distributed globally from humans (blue labels), animals (red labels), plants (purple/mauve labels) and the environment (green labels).

#### Providencia species

There were five *Providencia spp*. viz., *stuartii, rustigianii, alcalifaciens, heimbachae*, and *rettgeri*, with *P. rettgeri* being most isolated and branching into two clades, A and B (Fig. 10). The other species had single clades. *Providencia spp*. were isolated from animals, humans, and the environment. Highly similar strains of the same species were found across countries (*P. stuartii, P. alcalifaciens, P. heimbache* and *P. rettgeri*) and within countries (*P. rettgeri* clade A). The distinction between the various species of Providencia was depicted by the clustering patterns on the tree as strains of the same species clustered together; *Providencia spp*. was closest evolutionarily to *S. marcescens* whilst *E. coli* was closest to *Citrobacter spp*.,and *E. hormaechei*. As well, *K. pneumoniae, K. aerogenes*, and *K. michiganensis* clustered together with relatively short evolutionary distance (Fig. 10).

![Figure 10.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F10.medium.gif)

[Figure 10.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F10)

Figure 10. 
Evolutionary epidemiology and resistome of global *Providencia spp*. isolates. The *Providencia spp*. clustered into 12 branches, with *P. stuartii, P. rustigianii, P. alcalifaciens, P. heimbachei* and *P. rettgeri* clustering into branch VIII to XII respectively. *P. rettgeri* clustered into clades A and B, consisting of globally distributed isolates. *P. stuartii* and *P. rettgeri* contained richer and more abundant resistomes than the other *Providencia* strains and contained strains distributed globally from humans (blue labels), animals (red labels), plants (purple/mauve labels) and the environment (green labels).

The other Enterobacteriaceae species had more ARGs than *Providencia spp*. Within *Providencia spp., P. stuartii* and *P. rettgeri* were most endowed with richer and more diverse resistomes whilst *P. alcalifaciens* and *P. rustigianii* were almost bereft of ARGs. Common ARGs within *P. stuartii*, were *aac(2’)-Ia*, aph(3)-like, *bla*OXA, *bla*TEM, *aadA1, dfrA/B, sul1, fosA, ampC, tet* and *catA/B* whilst *P. rettgeri* had *aac(6’)-like, aph(3’/3”)-like, bla*SRT, *bla*OXA, *bla*NDM, *aadA, dfrA/B, OqxAB, sul1, qnrA/B/D/E/S, ampC, catA/B, tet, arr*, mph(A/E), *cmlA*, msr(E), and Inu(F/G). Hence, *P. rettgeri* has the most abundant and diverse ARGs than all other *Providencia spp*. (Fig. 10).

#### Phylogeography

North America (particularly USA) and Europe (particularly Western and Southern Europe) had the highest concentration of the various species, followed by South East Asia, South America (particularly Brazil and Colombia) and South Africa. There were sparse reports on these species from Australasia, the Middle East and Africa (except South Africa) and the Caribbean (Fig. 11).

![Figure 11.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F11.medium.gif)

[Figure 11.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F11)

Figure 11. 
Global phylogeography of *Citrobacter freundii, Citrobacter amalonaticus, Citrobacter koseri, Enterobacter hormaechei subsp. hormaechei, xiangfangensis, steigerwaltii* and *oharae, Klbesiella variicola, Morganella morganii, Proteus mirabilis*, and *Providencia spp*. Most of these genomes were obtained from USA, Europe, South-East Asia and South America in a descending order of frequency. *C. freundii, Enterobacter spp., K. variicola*, and *P. mirabilis* had more diverse distribution across the globe. Each species is designated with a different colour code.

*C. freundii* clade A was distributed mainly in North America, Europe, and South-East Asia whilst clade B was found worldwide on almost all continents. *C. amalonaticus* clade A was found in North America and South Korea; clade B was found in Malawi, USA, Malaysia, and France whilst clade C was only found in Malaysia, USA and Switzerland. *C. koseri* clade A (USA) and clade B (USA, Spain, UK, France, Canada, China and Malaysia) were relatively less reported with clade B being more widely distributed globally (Fig. 11).

*E. oharae* strains and *E. hormaechei* clades A, B and C were globally disseminated, with clade C being most widely distributed, followed by clades B and A. *E. steigerwaltii* clade C was more globally distributed than clades A (USA, Japan and Europe) and B (Argentina, USA, France and Germany). Although *E. xiangfangensis* was globally disseminated, the respective clades A (China, India, USA, S. Africa, and France) and B (France, USA, Egypt, and China) were reported from very few countries (Fig. 11).

*K. variicola* strains, particularly clades B and C, were of wide geographical distribution; clade A was found in relatively fewer countries globally. *M. morganii* strains were found in North America, Europe (including Russia), and South-East Asia, with clades B and C being found in South Africa. *P. mirabilis* clades A, B and C were found globally, with clade C being reported in most countries. *P. rettgeri* strains had the widest global distribution and reports from most countries among the *Providencia spp*., followed by *P. stuartii. P. rettgeri* clade A (China, Brazil, South Africa, Colombia, and USA) was found in fewer countries than clade B. *P. stuartii* strains were of global distribution whilst *P. rustigianii* (UK and USA), *P. alcalifaciens* (USA and India) and *P. heimbachae* (France, China, and Germany) were not (Fig. 11).

### Frequency distribution of ARGs per species

#### MCR ARGs

*Mcr-9.1* ARGs were the commonest *mcr* variants identified, with very few *mcr-1* and *mcr-3* being found in *E. hormaechei* and *C. freundii (mcr-1, -3* and *-10)* and a single *mcr-4.3* gene being found in *P. rettgeri*. Notably, the highest prevalence of *mcr-9* was in *E. hormaechei* (n=67 *mcr-9* genes), *E. steigerwaltii/oharae* (n=32 *mcr-9* genes) and *E. xiangfangensis/cloacae* (n=19 *mcr-9* genes), followed by *C. freundii* (n=19), *C. amalonaticus* (n=5) and other *Citrobacter spp*., some of which had very few or no *mcr* genes (Fig. S4-S9).

#### Carbapenemases

One of the most prevalent carbapenemase among the species was KPC, with KPC-2 (n=123), KPC-3 (n=97), KPC-4 (n=14), and KPC-6 (n=1) being common. KPC-2 was higher in all the species except in *E. xiangfangensis* for which KPC-3 was more abundant (n=35) than KPC-2 (n=9). KPC was most prevalent in *E. hormaechei/oharae/steigerwaltii/xiangfangensis* and *C. freundii* strains whilst *C. amalonaticus, C. koseri*, and *Providencia spp. had no* KPC ARGs. After Ambler class A KPC, Ambler class D OXA-48-like serine carbapenemases (n=65) were also very prominent in all species except *C. amalonaticus, E. xiangfangensis* and *K. variicola*, with other 0XA-48 variants such as OXA-181 (n=1, *M. morganii)* and OXA-396 (n=1, *P. rettgeri)* being relatively scarce. *C. freundii* (n=43), *E. hormaechei* (n=10), *E. oharae/steigerwaltii* (n=9) and *C. koseri* (n=3) had OXA-48 genes whilst OXA-58 (n=2) and OXA-23 (n=23) were only found in *P. mirabilis*. Other class A serine carbapenemases i.e., GES-5 (n=1, *M. morganii)* and IMI-2 (n=1, *E. hormaechei)*, were also rare (Fig. S4-S9).

NDM was the commonest class B carbapenemase (n=159), followed by IMP (n=97) and VIM (n=83). NDM-1 (n=137) was the most prevalent variant and was found in *E. hormaechei* (n=68), *C. freundii* (n=15), *E. steigerwaltii* (n=13), *Providencia spp*. (n=12), *P. mirabilis* (n=10), *E. xiangfangensis* (n=6), *C. amalonaticus* (n=5), *M. morganii* (n=5), and *K. variicola* (n=3), with NDM-5 (n=16; 12 in *E. hormaechei*, 3 in *E. xiangfangensis*, and 1 in *P. mirabilis*), NDM-7 (n=3 in *E. hormaechei*) and NDM-9 (n=3 in *K. variicola*) being less prevalent. IMP-1 (n=71; 67 in *E. hormaechei* and 4 in *C. freundii)*, IMP-8 (n=14; 13 in *C. freundii*, 1 in *E. steigerwaltii)*, IMP-4 (n=9; 5 in *C. freundii*, 4 in *E. hormaechei)*, IMP-27 (n=2 in *P. mirabilis)* and IMP-13 (n=1 in *E. oharae)* were the identified variants. VIM-1 (n=65; 26 in *C. freundii*, 14 in *E. steigerwaltii*, 21 in *E. hormaechei*, 1 in *E. xiangfangensis*, 2 in *Providencia spp*., 1 in *P. mirabilis)*, VIM-4 (n=12; 5 in in *E. steigerwaltii*, 5 in *E. hormaechei*, 2 in *C. freundii)*, VIM-2 (n=2 in *Providencia spp*.), VIM-31 (n=2; *E. steigerwaltii* and *E. hormaechei)*, VIM-5 (n=1 in *E. hormaechei)*, and VIM-67 (n=1 in *E. hormaechei*) were identified in the strains (Fig. S4-S9).

#### ESBLs and ampCs

TEM was the commonest ESBL and TEM-1 was the most common variant to be identified in all species; particularly, TEM-1 was most abundant in *Enterobacter spp*. and *C. freundii*. OXA (particularly OXA-1, -9, and -10), CTX-M (particularly CTX-M-15) and SHV (particularly SHV-12) were also common in almost all species in a descending order of prevalence, but were very abundant in *Enterobacter spp*., specifically *E. hormaechei; C. koseri, C. amalonaticus, P. mirabilis*, and *M. morganii* had relatively low abundance of SHV, OXA and CTX-M genes. Other ESBLs genes such as *bla*SCO, *bla*LAP, *bla*VEB, *bla*TMB, *bla*SFO, *bla*SMB, *bla*CARB, *bla* etc. were also rare in the various species (Fig. S4-S9).

AmpC ARGs were basically strain-specific, with ACT, LEN, CMY, DHA, CKO/MAL and SED being conserved in *Enterobacter spp., K. variicola, C. fruendii, M. morganii, C.koseri* and *C. amalonaticus* (except clade A) respectively; FOX was rare (Fig. S4-S9).

#### Aminoglycoside ARGs

ARGs mediating resistance to aminoglycosides such as *aadA, aph(2”)-like, aph(3’)-like, aph(4’)-like, aac(6’)-like, aac(3)-like, aph(6)-like*, and *ant(2”)-like* were abundantly prevalent in almost all the clades of *C. freundii, Enterobacter spp*., and *P. mirabilis* and sparsely abundant in the other species; *aadA, aph(4’)-like, aac(6’)-like* and *aac(3’)-like* ARGs were most common. 16S rRNA Methyltransferases such as *rmtB1 (E. steigerwaltii* and *P. mirabilis*,), *rmtC/G (E. hormaechei, C. freundii* and *M. morganii)*, and *armA (C. freundii, E. xiangfangensis, P. mirabilis, P. stuartii, M. morganii* and *P. rettgeri*) were rare (Fig. S4- S9).

#### Fluoroquinolone ARGs

*Aac(6’)-like, OqxAB, QnrA/B/D/S* and *qepA* ARGs were identified, albeit *qepA (M. morganii)* was rare and *Oqx*A was less prevalent than *Oqx*B in all but one species. Notably, *OqxAB* were virtually absent in *C. fruendii* clades. whilst chromosomal mutations in *gyr*AB and *parCE* were only observed in *K. variicola* (Fig. S4-S9).

#### Other ARGs

Chloramphenicol ARGs, *cmlA, catA* and *catB* were found in all the species, albeit *cmlA* was relatively rare in all the species whilst *catA/B* were conserved in *Enterobacter spp*. and *M. morganii; catA* was more prevalent than *catB. catA/B* were also abundant in *Citrobacter spp., P. mirabilis* and *Providencia spp*. (Fig. S4-S9). *Sulphamethoxazole-trimethoprim* ARGs, *sul-1/2/3* and *dfrA*, were enriched in *Providencia spp., P. mirabilis, M. morganii, Enterobacter spp., C. koseri (dfrA was virtually absent), C. amalonaticus* clade B2, and *C. freundii. Sul1* was more prevalent than *Sul2*, with *Sul3* being relatively rare whilst the *dfrA* variants were remarkably diverse. Indeed, both *sul1, sul2* and/or *sul3* as well as several *dfrA* variants were present concurrently in some single strains (Fig. S4-S9; Tables S2).

Several tetracycline ARGs such as *tet*(A), *tet*(B), *tet*(C), *tet*(D), *tet*(G), *tet*(*J*), *tet*(S), *tet*(Y), and *tet*(41), were present in all the species. Notably, *tet*(A), *tet*(B), and *tet*(D), were highly enriched in the various genomes with a descending order of frequency; *tet*(J) was most prevalent in *P. mirabilis* (Fig. S4-S9).*fosA* variants were present in all the species except *M. morganii*, but were most enriched and conserved in *Enterobacter spp. K. variicola, P. stuartii* and *P. rettgeri*. As well, *ere*(A), *emr*(D), *erm*(B), *msr*(E), *mph*(A) and *mph*(E) macrolide ARGs were common in the various species, with *mph*(A) being richly abundant; *mph*(E) and *msr*(E) were enriched in *C. freundii, Providencia spp*. and *P. mirabilis, emr*(D) was abundant in *K. variicola* whilst *ere*(A) was abundantly enriched in *E. hormaechei* and *E. xiangfangensis*. Rifamycin ARG, *arr*, was identified in the various species, represented by *arr-S arr-2* and *arr-1* in all the species (Fig. S4-S9).

## Discussion

Among Gram-negative bacterial species, *Pseudomonas aeruginosa, Acinetobacter baumannii*, and Enterobacteriaceae such as *K. pneumoniae, E. coli, S. enterica, Vibrio cholerae*, and *Shigella spp*. are commonly isolated and implicated in nosocomial (anthroponotic), zoonotic and water/food-borne infections with MDR, extensively and pandrug resistant (XDR and PDR respectively) phenomes 4,26,28,37,38. However, other less isolated Enterobacteriaceae species such as *Citrobacter spp., Enterobacter hormaechei subsp., K. variicola, P. mirabilis, M. morganii and Providencia spp*., are increasingly being implicated in MDR, XDR and PDR infections globally as opportunistic pathogens 7,8,28,34,39–41. We show herein, that *C. freundii, Enterobacter hormaechei subsp. hormaechei, xiangfangensis, oharae* and *steigerwaltii*, and Proteeae strains harbour multiple resistance mechanisms that can make them MDR, XDR and PDR pathogens. More concerning is the global distribution and multiple (human, animal, plants and environmental) specimen sources of these strains, which suggest that they can cause anthroponotic, zoonotic and food- and water-borne infections 42

Hence, these opportunistic pathogens demand more attention than they have been given hitherto as the rich resistome repertoire identified in their genomes makes them reservoirs of ARGs 11,14,21,43. Moreover, being intestinal denizens and commensals, they can easily share these ARGs with facultative and obligate pathogens of humans and animals 22,44–46. Further, their presence on plants and the environment further suggests that they can share their ARGs with food-borne and water-borne pathogens 44–46. Fortunately, the *E. xiangfangensis* strains found in rice from India had very few ARGs, albeit a few had multiple ARGs (Fig. S2B).

Of greater concern is the rich resistome repertoire and abundance of globally distributed *E. hormaechei subsp*. strains. Specifically, *E. hormaechei subsp*. contained clinically important ARGs such as *mcr-9*, carbapenemases, and ESBLs, alongside fluoroquinolones, aminoglycoside, tetracycline, macrolide, fosfomycin, chloramphenicol, rifamycin and sulphomethoxazole-trimethoprim resistance mechanisms. This resistome repertoire was also seen in *C. freundii, P. rettgeri P. mirabilis, P. stuartii*, and *M. morganii*, to a relatively lesser degree in a descending order. This is a worrying observation as colistin, carbapenems and tigecycline are last resort antibiotics used to treat fatal bacterial infections 47,48. In most cases, these antibiotics are used in combination with fosfomycin, fluoroquinolones and aminoglycosides to treat carbapenem-resistant Enterobacteriaceae (CRE) infections 47,48. Evidently, the presence of all these resistance mechanisms to these antibiotics, could make these species PDR and automatically qualify them as critical priority pathogens per the WHO criteria32,33.

As well, members of the tribe Proteeae viz., *M. morganii, P. mirabilis*, and *Providencia spp*., are known to have intrinsic resistance to colistin, tigecycline, aminopenicillins, amikacin, tobramycin, lincosamides, macrolides, fosfomycin and first- and second-generation cephalosporins 1,7,49,50. Thus, the presence of additional resistance determinants in this tribe is especially worrying. Already, there are increasing reports on the isolation of Proteeae species in recurrent urinary tract infections (UTIs) infections, which is facilitated by the increasing use of colistin to treat MDR infections; their broad intrinsic resistance mechanisms enable them to flourish during antibiotic chemotherapy 1,7,8,49,51,52 These observations evince the growing threat of antimicrobial resistance globally and its associated after-effects, supporting the need for efficient antibiotic stewardship to safeguard current antibiotic arsenals as well as discover novel ones 3,14.

As shown in the phylogenomic and phylogeographic analyses, local and international transmission, or outbreaks of strains within the same clone, clade and subclade i.e., of very close evolutionary distance, were observed. Notably, these closely related strains were isolated from humans, animals, plants, and the environment and they harboured important resistance determinants as described above. The phylogenomics showed the gradual evolution of a single strain during dissemination from host to host and depict the fact that antibiotic resistance respects no boundaries. Notably, a large part of *E. hormaechei* clade C consisted of strains from human stools/urine in Japan; these were closely related strains evolving from Spain, Taiwan, China, and Greece (Fig. 6 & S2). A similar observation was made with strains of close evolutionary distance from humans in Nigeria, France, Spain, Portugal, and Lebanon as well as with strains from humans, animals, and the environment in several countries such as France, Germany, USA, Pakistan, Morocco, Lebanon, and Poland in clade C. *E. hormaechei* strains were also found in desert sands in Morocco, showing their broad and diverse niches (Fig. 6 & S2).

Uniform and non-uniform resistome patterns were seen between strains of the same clade/clone in almost all the species. For instance, the same resistome was seen in *C. amalonaticus* clade B2 (Fig. 2), *E. cloacae* clade B (Fig. 5) and *P. mirabilis* clade C2 (Fig. 9) whilst differing resistome patterns were observed in the other species and clades. This observation supports two phenomena: firstly, the clonal expansion of strains harbouring the same resistome repertoire on both chromosomes and plasmids and secondly, the horizontal transmission of genetic elements bearing the same of different ARGs across clones and species. During clonal expansion of strains, there is a concomitant replication of resistance plasmids alongside chromosomal replication, leading to daughter cells with the same resistome diversity 11,12,21,53. As well, horizontal gene transfer of ARGs between bacteria can alter the resistome diversity and composition of daughter clones emanating from the same ancestor 11,12,21,53. In this case, both phenomena are being observed, showing that both vertical and horizontal transmission of ARGs are ongoing during the evolutionary epidemiology of the various clades and species across the globe.

The presence of multiple ARGs in a single strain might not necessarily mean they are all being expressed in the bacteria’s phenome as antibiotic-susceptible strains have been found to harbour ARGs. For instance, colistin- and fosfomocyin-sensitive Enterobacteriaceae have have been reported in strains harbouring the *mcr-9* and *fosA* genes 1,5,7,8. Nevertheless, the ability of these ARGs to be expressed in the presence of strong promoter or transferred to another host with a stronger promoter for subsequent expression cannot be gainsaid 54 Indeed, antibiotic abuse could serve as an inducer to trigger the transfer and expression of these vast resistomes in the microbial phenomes 49,55, necessitating the importance for judicious antibiotic use.

It is revealing to note that *K. variicola, C. amalonaticus* and *C. koseri* strains had very few ARGs except for *K. variicola* clade B2 and *C. amalonaticus* clade B2, despite the global distribution of *K. variicola* (Fig. 2-3, 7, & S3). Notably, *C. amalonaticus* clade B2, which were all from France, were remarkably enriched with ARGs including *bla*NDM, representing a local outbreak of XDR *C. amalonaticus* strains (Fig. 2). Thus, even in species with fewer resistome diversity and abundance, there are MDR, XDR and PDR strains causing local outbreaks.

Notably, most of the genomes included in this analysis were from the USA, Europe, and South East Asia. This may be due to the fact that these regions have higher prevalence and incidence of infections resulting from these pathogens or that these areas have enough financial and technical means to undertake genomic sequencing of these isolates in periodic surveillance studies. Specifically, genomes of these species were relatively scarce from a large part of Russia, Middle and North-West Asia, Africa, the Caribbean and parts of South America and Canada. Given the alarming resistome diversity and composition realised in this analyses, it is incumbent for all nations to intensify and adopt genome-based epidemiological studies to quickly identify the sources and reservoirs of ARGs to pre-empt outbreaks of MDR, XDR and PDR pathogens.

### Conclusion

In conclusion, less described Enterobacteriaceae species viz., *Enterobacter hormaechei subsp. hormaechei, xiangfangensis, steigerwaltii* and *oharae, C. freundii, M. morganii, P. mirabilis, P. stuartii* and *P. rettgeri*, contain globally distributed MDR, XDR and PDR strains that cause local and international outbreaks, transmitting through humans, animals, food, water and other environmental media or sources. Notably, the resistome repertoire of these relatively rare species were equally or more abundant and diverse as those of commonly isolated species. Hence, intensified efforts should be made to increase education on antibiotic stewardship to safeguard the potency of available antibiotics and reduce the selection and dissemination antibiotic-resistant Enterobacteriaceae. Infection prevention and control as well periodic genomic surveillance of communities, hospitals, farms, water bodies and the general environment (One Health) should be undertaken to pre-empt outbreaks of MDR strains and inform infection control measures.

Notwithstanding the revealing details obtained in this study, strains with clinically important ARGs whose genomes are not deposited in NCBI/PATRIC or whose genomes are not sequenced before January 2020 will be missed; hence, the information contained herein are true up to January 2020. Nevertheless, the global phylogeography and resistome epidemiology of these emerging opportunistic pathogens provide an important picture of the ARGs, sources and transmission patterns of these species.

## Methods

### Included genomes

Genomes of *Citrobacter spp*. (including *amalonaticus, freundii, koseri, werkmanii, brakii, portucalensis* and *youngae), Enterobacter hormaechei subsp. hormaechei, xiangfangensis, steigerwaltii*, and o*harae, Providencia spp*. (including *alcalifaciens, burhodogranariea, heimbachae, rettgeri, rustigianii*, and *stuartii)* and *Proteus mirabilis* deposited at GenBank ([https://www.ncbi.nlm.nih.gov/genbank/](https://www.ncbi.nlm.nih.gov/genbank/)) and PATRIC ([https://www.patricbrc.org/](https://www.patricbrc.org/)) up to January 2020 were pooled and filtered to remove poor genome sequences. These were used for the downstream phylogenetics, phylogeography and resistome analyses.

### Phylogenetics and evolutionary epidemiology analyses

Briefly, plasmid sequences, phages, poor genomes i.e., genomes with sizes below the average genome size of each species viz., 3-4Mb, and genomes of strains that could not share at least 1000 core protein genes with all the included genomes were removed. The remaining genomes were aligned and run through RAXmL in batches of 200 genomes to draw phylogenetic trees using the maximum-likelihood method. A minimum of 1000 genes were used for the alignment and a bootstrap resampling of 1000x was used.The trees were annotated using Figtree to show their sequence type (ST), host (species), country and year of isolation. The various clades and sub-clades within each species or genera were visually identified based on their clustering distance and accordingly labelled.

### Phylogeography

The various clades and subclades per species were manually drawn unto maps to show their phylogeographic distribution using Paint 3D. Different colour codes were used to distinguish between the various species and clades.

### Resistome analyses

The resistomes of the included genomes were individually obtained from the NCBI Pathogen Detection database ([https://www.ncbi.nlm.nih.gov/pathogens/isolates#/search/](https://www.ncbi.nlm.nih.gov/pathogens/isolates#/search/)). The resistomes were aligned per strain and colour-coded per clade or species to show their association per species, clone, or clade. These were then associated with the phylogenomic trees to ascertain the resistome dynamics per clone, clade, species, and geographical location.

## Data Availability

All data are included as supplemental files

## Tweet

Global *E. hormaechei, C. freundii, P. mirabilis, P. stuartii* & *P. rettgeri* strains contain multiple resistance mechanism to important and reserved antibiotics in globally circulating clones. These portend the dawn of pandrug resistance and a return to the pre-antibiotic era.

## Highlights/Importance

*Citrobacter spp., Enterobacter hormaechei subsp., Klebsiella variicola* and *Proteae* tribe members are rarely isolated Enterobacteriaceae increasingly implicated in nosocomial infections. The global phylogenomics, evolution, geographical distribution and resistome repertoire of these species found them to be globally distributed, being isolated from human, animal, plant, and environmental sources. *E. hormaechei subsp., C. freundii, Proteus mirabilis, Providencia stuartii, Providencia rettgeri* and *Morganella morganii* contained multidrug-resistant clades that harboured resistance to clinically important reserved antibiotics, portending the dawn of pandrug resistance and a potential acceleration of the post-antibiotic era.

## Funding

None

## Transparency declaration

None (Authors declare no conflict of interest)

![Supplemental dataset 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F12.medium.gif)

[Supplemental dataset 1.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F12)

Supplemental dataset 1. 
Raw metadata of downloaded genomes from PATRIC containing all the data associated with each genome.

![Supplemental dataset 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F13.medium.gif)

[Supplemental dataset 2.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F13)

Supplemental dataset 2. 
Species by species tabulation and analyses of the resistomes, specimen sources, country of isolation, MLST, Biosample accession number, and strain name of all the genomes according to their order on the phylogeny trees.

![Supplemental dataset 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F14.medium.gif)

[Supplemental dataset 3.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F14)

Supplemental dataset 3. 
Colour-coded species by species tabulation of the resistomes, specimen sources, country of isolation, MLST, Biosample accession number, and strain name of all the genomes according to their order on the phylogeny trees.

![Figure S1 (A and B).](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F15.medium.gif)

[Figure S1 (A and B).](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F15)

Figure S1 (A and B). 
Evolutionary epidemiology and resistome of global *Citrobacter freundii* isolates, A and B. *C. freundii* clustered into four main clades (A, B1, B2 and B3), highlighted with distinct colours. Clade B3 had the most resistome abundance and diversity. Strains from humans (blue labels), animals (red labels), plants (purple/mauve labels) and the environment (green labels) were found in the same clade/cluster. *Bla*CMY was conserved in these genomes.

![Figure S2 (A and B).](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F16.medium.gif)

[Figure S2 (A and B).](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F16)

Figure S2 (A and B). 
Evolutionary epidemiology and resistome of global *Enterobacter hormaechei* isolates, A and B. The *E. hormaechei* isolates clustered into three main clades A, B and C (with distinct highlights) that contained strains distributed globally from humans (blue labels), and animals (red labels), plants (purple/mauve labels) and the environment (green labels). Clades B and C contained diverse and rich resistome repertoire. *bla*ACT was conserved in these genomes.

![Figure S3 (A and B).](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F17.medium.gif)

[Figure S3 (A and B).](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F17)

Figure S3 (A and B). 
Evolutionary epidemiology and resistome of global *Klebsiella variicola* isolates, A and B. The *K. variicola* strains clustered into nine clades viz., A1, A2, A3, A4, A5, A6, B1, B2 and C, which were highlighted with distinct colours and were isolated from countries around the globe. The clades contained strains distributed globally from humans (blue labels), animals (red labels), plants (purple/mauve labels) and the environment (green labels). Besides a few strains in clade B2, the other strains contained very few resistance genes. *bla*LEN was conserved in these genomes.

![Figure S4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F18.medium.gif)

[Figure S4.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F18)

Figure S4. 
Frequency distribution of antibiotic resistance genes in *Citrobacter freundii* (A), and *Citrobacter species* (B and C).

![Figure S5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F19.medium.gif)

[Figure S5.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F19)

Figure S5. 
Frequency distribution of antibiotic resistance genes in *Citrobacter amalonaticus* (A), and *Citrobacter koseri* (B).

![Figure S6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F20.medium.gif)

[Figure S6.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F20)

Figure S6. 
Frequency distribution of antibiotic resistance genes in *Enterobacter steigerwaltii* and *oharae* (A), and *Enterobacter xiangfangensis* (B).

![Figure S7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F21.medium.gif)

[Figure S7.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F21)

Figure S7. 
Frequency distribution of antibiotic resistance genes in *Enterobacter hormaechei* (A, B and C).

![Figure S8.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F22.medium.gif)

[Figure S8.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F22)

Figure S8. 
Frequency distribution of antibiotic resistance genes in *Klebsiella variicola* (A, B and C).

![Figure S9.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.21.20109504/F23.medium.gif)

[Figure S9.](http://medrxiv.org/content/early/2020/05/26/2020.05.21.20109504/F23)

Figure S9. 
Frequency distribution of antibiotic resistance genes in *Morganella morganii* (A), *Proteus mirabilis* (B) and *Providencia species*(C).

## Acknowledgements

None

*   Received May 21, 2020.
*   Revision received May 21, 2020.
*   Accepted May 26, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Mbelle, N. et al. Genomic analysis of two drug-resistant clinical Morganella morganii strains isolated from UTI patients in Pretoria, South Africa. Lett. Appl. Microbiol. 70, 21–28 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/lam.13237&link_type=DOI) 

2.  2.Osei Sekyere, J. Genomic insights into nitrofurantoin resistance mechanisms and epidemiology in clinical Enterobacteriaceae. Futur. Sci. OA 4, FSO293 (2018).
    
    

3.  3.Osei Sekyere, J. & Asante, J. Emerging mechanisms of antimicrobial resistance in bacteria and fungi: advances in the era of genomics. Future Microbiol. 13, 241–262 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2217/fmb-2017-0172&link_type=DOI) 

4.  4.Pedersen, T. et al. Spread of Plasmid-Encoded NDM-1 and GES-5 Carbapenemases among Extensively Drug-Resistant and Pandrug-Resistant Clinical Enterobacteriaceae in Durban, South Africa. Antimicrob. Agents Chemother. 62, e02178–17 (2018).
    
    

5.  5.Mbelle, N. M. et al. Pathogenomics and Evolutionary Epidemiology of Multi-Drug Resistant Clinical Klebsiella pneumoniae Isolated from Pretoria, South Africa. Sci. Rep. 10, 1–17 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-020-59121-0&link_type=DOI) 

6.  6.Osei Sekyere, J. & Mensah, E. Molecular epidemiology and mechanisms of antibiotic resistance in Enterococcus spp., Staphylococcus spp., and Streptococcus spp. in Africa: a systematic review from a One Health perspective. Ann. N. Y. Acad. Sci. 1465, 29–58 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/nyas.14254&link_type=DOI) 

7.  7.Mbelle, N. M. et al. Genomic analysis of a multidrug-resistant clinical Providencia rettgeri (PR002) strain with the novel integron ln1483 and an A/C plasmid replicon. Ann. N. Y. Acad. Sci. 1462, 92–103 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/nyas.14237&link_type=DOI) 

8.  8.Osei Sekyere, J., Maningi, N. E., Modipane, L. & Mbelle, N. M. Emergence of mcr-9.1 in ESBL-producing Clinical Enterobacteriaceae in Pretoria, South Africa: Global Evolutionary Phylogenomics, Resistome and Mobilome. John. mSystems 5, e00148–20 (2020).
    
    

9.  9.Mbelle, N. M. et al. The Resistome, Mobilome, Virulome and Phylogenomics of Multidrug-Resistant Escherichia coli Clinical Isolates from Pretoria, South Africa. Sci. Rep. 9, 1–43 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-019-40055-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30862893&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F26%2F2020.05.21.20109504.atom) 

10. 10.Chirindze, L. M. et al. Faecal colonization of E. coli and Klebsiella spp. producing extended-spectrum beta-lactamases and plasmid-mediated AmpC in Mozambican university students. BMC Infect. Dis. 18, 244 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12879-018-3154-1&link_type=DOI) 

11. 11.Decano, A. G. et al. Plasmids shape the diverse accessory resistomes of Escherichia coli ST131. *bioRxiv* 2020.05.07.081380 (2020). doi:10.1101/2020.05.07.081380
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4wNS4wNy4wODEzODB2NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA1LzI2LzIwMjAuMDUuMjEuMjAxMDk1MDQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

12. 12.Partridge, S. R., Kwong, S. M., Firth, N. & Jensen, S. O. Mobile Genetic Elements Associated with Antimicrobial Resistance. Clin. Microbiol. Rev. 31, 1–61 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1128/CMR.00020-18&link_type=DOI) 

13. 13.Böhm, M.-E., Razavi, M., Marathe, N. P., Flach, C.-F. & Larsson, D. G. J. Discovery of a novel integron-borne aminoglycoside resistance gene present in clinical pathogens by screening environmental bacterial communities. Microbiome 8, 41 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s40168-020-00814-z&link_type=DOI) 

14. 14.Asante, J. & Osei Sekyere, J. Understanding antimicrobial discovery and resistance from a metagenomic and metatranscriptomic perspective: Advances and applications. Environ. Microbiol. Rep. 11, 62–86 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/1758-2229.12735&link_type=DOI) 

15. 15.Van Boeckel, T. P. et al. Global trends in antimicrobial resistance in animals in low- and middle-income countries. Science (80-.). 365, eaaw1944 (2019).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE3OiIzNjUvNjQ1OS9lYWF3MTk0NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA1LzI2LzIwMjAuMDUuMjEuMjAxMDk1MDQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

16. 16.He, T. et al. Emergence of plasmid-mediated high-level tigecycline resistance genes in animals and humans. Nat. Microbiol. 4, 1450–1456 (2019).
    
    

17. 17.Li, J. et al. Inter-host Transmission of Carbapenemase-Producing *Escherichia coli* among Humans and Backyard Animals. Environ. Health Perspect. 127, 107009 (2019).
    
    

18. 18.Jaidane, N. et al. Genomic analysis of in vivo acquired resistance to colistin and rifampicin in Acinetobacter baumannii. Int. J. Antimicrob. Agents 51, 266–269 (2018).
    
    

19. 19.Olaitan, A. O., Diene, S. M., Assous, M. V. & Rolain, J.-M. Genomic Plasticity of Multidrug-Resistant NDM-1 Positive Clinical Isolate of Providencia rettgeri. Genome Biol. Evol. 8, 723–8 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/gbe/evv195&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26454012&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F26%2F2020.05.21.20109504.atom) 

20. 20.Conlan, S. et al. Plasmid Dynamics in KPC-Positive Klebsiella pneumoniae during Long-Term Patient Colonization. MBio 7, (2016).
    
    

21. 21.Decano, A. G. et al. Complete Assembly of Escherichia coli Sequence Type 131 Genomes Using Long Reads Demonstrates Antibiotic Resistance Gene Variation within Diverse Plasmid and Chromosomal Contexts. mSphere 4, 130–19 (2019).
    
    

22. 22.Osei Sekyere, J., Maningi, N. E. & Fourie, P. B. Mycobacterium tuberculosis, antimicrobials, immunity, and lung–gut microbiota crosstalk: current updates and emerging advances. Ann. N. Y. Acad. Sci. 1467, 21–47 (2020).
    
    

23. 23.Liu, B.-T. et al. Characteristics of Carbapenem-Resistant Enterobacteriaceae in Ready-to-Eat Vegetables in China. Front. Microbiol. 9, 1147 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fmicb.2018.01147&link_type=DOI) 

24. 24.Tshisevhe, V. S., Lekalakala, M. R., Tshuma, N., Janse van Rensburg, S. & Mbelle, N. Outbreak of carbapenem-resistant Providencia rettgeri in a tertiary hospital. S. Afr. Med. J. 107, 31–33 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7196/SAMJ.2017.v107i1.12002&link_type=DOI) 

25. 25.Kopotsa, K., Osei Sekyere, J. & Mbelle, N. M. Plasmid evolution in carbapenemase-producing Enterobacteriaceae: a review. Ann. N. Y. Acad. Sci. 1457, 61–91 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/nyas.14223&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31469443&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F26%2F2020.05.21.20109504.atom) 

26. 26.Weill, F. X. et al. Global phylogeography and evolutionary history of Shigella dysenteriae type 1. Nat. Microbiol. 1, (2016).
    
    

27. 27.Hargreaves, M. L. et al. Clonal Dissemination of Enterobacter cloacae Harboring bla KPC-3 in the Upper Midwestern United States. Antimicrob. Agents Chemother. 59, 7723–7734 (2015).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjEwOiI1OS8xMi83NzIzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMjYvMjAyMC4wNS4yMS4yMDEwOTUwNC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

28. 28.Osei Sekyere, J. & Reta, M. A. Phylogeography and Resistome Epidemiology of Gram-Negative Bacteria in Africa: A Systematic Review and Genomic Meta-Analysis from a One-Health Perspective. *medRxiv* 2020.04.09.20059766 (2020). doi:10.1101/2020.04.09.20059766
    
    

29. 29.Cerqueira, G. C. et al. Multi-institute analysis of carbapenem resistance reveals remarkable diversity, unexplained mechanisms, and limited clonal outbreaks. Proc. Natl. Acad. Sci. U. S. A. 114, 1135–1140 (2017).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiMTE0LzUvMTEzNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA1LzI2LzIwMjAuMDUuMjEuMjAxMDk1MDQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

30. 30.Tijet, N. et al. Lateral dissemination and inter-patient transmission of blaKPC-3: role of a conjugative plasmid in spreading carbapenem resistance. J. Antimicrob. Chemother. 71, 344–347 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jac/dkv356&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26518052&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F26%2F2020.05.21.20109504.atom) 

31. 31.Petty, N. K. et al. Global dissemination of a multidrug resistant Escherichia coli clone. Proc. Natl. Acad. Sci. U. S. A. 111, 5694–5699 (2014).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTExLzE1LzU2OTQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNS8yNi8yMDIwLjA1LjIxLjIwMTA5NTA0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

32. 32.Rello, J. et al. A global priority list of the TOp TEn resistant Microorganisms (TOTEM) study at intensive care: a prioritization exercise based on multi-criteria decision analysis. Eur. J. Clin. Microbiol. Infect. Dis. 38, 319–323 (2019).
    
    

33. 33.Asokan, G. V, Ramadhan, T., Ahmed, E. & Sanad, H. WHO Global Priority Pathogens List: A Bibliometric Analysis of Medline-PubMed for Knowledge Mobilization to Infection Prevention and Control Practices in Bahrain. Oman Med. J. 34, 184–193 (2019).
    
    

34. 34.Mbelle, N. et al. Genomic analysis of two drug-resistant clinical Morganella morganii strains isolated from UTI patients in Pretoria, South Africa. Lett. Appl. Microbiol. 70, 21–28 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/lam.13237&link_type=DOI) 

35. 35.Lv, L. et al. Emergence of a Plasmid-Encoded Resistance-Nodulation-Division Efflux Pump Conferring Resistance to Multiple Drugs, Including Tigecycline, in Klebsiella pneumoniae. MBio 11, (2020).
    
    

36. 36.Ellem, J. A. et al. Locally Acquired mcr-1 in Escherichia coli, Australia, 2011 and 2013. Emerg. Infect. Dis. 23, 1160–1163 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2307.161638&link_type=DOI) 

37. 37.Perez, A. et al. High incidence of MDR and XDR Pseudomonas aeruginosa isolates obtained from patients with ventilator-associated pneumonia in Greece, Italy and Spain as part of the MagicBullet clinical trial. J. Antimicrob. Chemother. 74, 1244–1252 (2019).
    
    

38. 38.D’Onofrio, V. et al. Epidemiology of colistin-resistant, carbapenemase-producing Enterobacteriaceae and Acinetobacter baumannii in Croatia. Infect. Genet. Evol. 81, 104263 (2020).
    
    

39. 39.Lob, S. H. et al. In vitro activity of imipenem-relebactam against resistant phenotypes of Enterobacteriaceae and Pseudomonas aeruginosa isolated from intraabdominal and urinary tract infection samples - SMART Surveillance Europe 2015–2017. J. Med. Microbiol. 69, 207–217 (2020).
    
    

40. 40.Lu, Y., Feng, Y., McNally, A. & Zong, Z. Occurrence of colistin-resistant hypervirulent Klebsiella variicola. J. Antimicrob. Chemother. (2018). doi:10.1093/jac/dky301
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jac/dky301&link_type=DOI) 

41. 41.Rodriguez-Medina, N., Barrios-Camacho, H., Duran-Bedolla, J. & Garza-Ramos, U. Klebsiella variicola: an emerging pathogen in humans. Emerg. Microbes Infect. 8, 973–988 (2019).
    
    

42. 42.Exner, M. et al. Antibiotic resistance: What is so special about multidrug-resistant Gram-negative bacteria? GMS Hyg. Infect. Control 12, Doc05 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3205/dgkh000290&link_type=DOI) 

43. 43.Yang, D. et al. The Occurrence of the Colistin Resistance Gene mcr-1 in the Haihe River (China). Int. J. Environ. Res. Public Health 14, (2017).
    
    

44. 44.Lerner, A., Matthias, T. & Aminov, R. Potential effects of horizontal gene exchange in the human gut. Front. Immunol. 8, (2017).
    
    

45. 45.Boto, L., Pineda, M. & Pineda, R. Potential impacts of horizontal gene transfer on human health and physiology and how anthropogenic activity can affect it. FEBS J. 286, 3959–3967 (2019).
    
    

46. 46.Anderson, M. T. & Seifert, H. S. Opportunity and means: Horizontal gene transfer from the human host to a bacterial pathogen. MBio 2, (2011).
    
    

47. 47.Osei Sekyere, J., Govinden, U., Bester, L. A. & Essack, S. Y. Colistin and tigecycline resistance in carbapenemase-producing Gram-negative bacteria: emerging resistance mechanisms and detection methods. J. Appl. Microbiol. 121, 601–617 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/jam.13169&link_type=DOI) 

48. 48.Osei Sekyere, J., Sephofane, A. K. & Mbelle, N. M. Comparative Evaluation of CHROMagar COL-APSE, MicroScan Walkaway, ComASP Colistin, and Colistin MAC Test in Detecting Colistin-resistant Gram-Negative Bacteria. Sci. Rep. 10, 1–13 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-020-59121-0&link_type=DOI) 

49. 49.Liu, H., Zhu, J., Hu, Q. & Rao, X. Morganella morganii, a non-negligent opportunistic pathogen. International Journal of Infectious Diseases 50, 10–17 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2016.07.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F26%2F2020.05.21.20109504.atom) 

50. 50.Hope, R. et al. Tigecycline activity: low resistance rates but problematic disc breakpoints revealed by a multicentre sentinel survey in the UK. J Antimicrob Chemother 65, 2602–2609 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jac/dkq370&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20930043&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F26%2F2020.05.21.20109504.atom) 

51. 51.Hayakawa, K. et al. Growing prevalence of Providencia stuartii associated with the increased usage of colistin at a tertiary health care center. Int. J. Infect. Dis. 16, e646–8 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2012.05.1029&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22818111&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F26%2F2020.05.21.20109504.atom) 

52. 52.Gajdács, M. & Urbán, E. Comparative epidemiology and resistance trends of Proteae in urinary tract infections of inpatients and outpatients: A 10-year retrospective study. Antibiotics 8, (2019).
    
    

53. 53.Carattoli, A. Plasmids in Gram negatives: Molecular typing of resistance plasmids. Int. J. Med. Microbiol. 301, 654–658 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijmm.2011.09.003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21992746&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F26%2F2020.05.21.20109504.atom) 

54. 54.Carroll, L. M. et al. Identification of novel mobilized colistin resistance gene mcr-9 in a multidrug-resistant, colistin-susceptible salmonella enterica serotype typhimurium isolate. MBio 10, e00853–19 (2019).
    
    

55. 55.Liu, G., Bogaj, K., Bortolaia, V., Olsen, J. E. & Thomsen, L. E. Antibiotic-Induced, Increased Conjugative Transfer Is Common to Diverse Naturally Occurring ESBL Plasmids in Escherichia coli. Front. Microbiol. 10, 2119 (2019).