Within-patient plasmid dynamics in Klebsiella pneumoniae during an outbreak of a carbapenemase-producing Klebsiella pneumoniae

Abstract

Introduction

Knowledge of within-patient dynamics of resistance plasmids during outbreaks is important for understanding the persistence and transmission of plasmid-mediated antimicrobial resistance. During an outbreak of a Klebsiella pneumoniae carbapenemase-producing (KPC) K. pneumoniae, the plasmid and chromosomal dynamics of K. pneumoniae within-patients were investigated.

Methods

During the outbreak, all K. pneumoniae isolates of colonized or infected patients were collected, regardless of their susceptibility pattern. A selection of isolates was short-read and long-read sequenced. A hybrid assembly of the short-and long-read sequence data was performed. Plasmid contigs were extracted from the hybrid assembly, annotated, and within patient plasmid comparisons were performed.

Results

Fifteen K. pneumoniae isolates of six patients were short-read whole-genome sequenced. Whole-genome multi-locus sequence typing revealed a maximum of 4 allele differences between the sequenced isolates. Within patients 1 and 2 the resistance gene- and plasmid replicon-content did differ between the isolates sequenced. Long-read sequencing and hybrid assembly of 4 isolates revealed loss of the entire KPC-gene containing plasmid in the isolates of patient 2 and a recombination event between the plasmids in the isolates of patient 1. This resulted in two different KPC-gene containing plasmids being simultaneously present during the outbreak.

Conclusion

During a hospital outbreak of a KPC-producing K. pneumoniae isolate, plasmid loss of the KPC-gene carrying plasmid and plasmid recombination was detected within the isolates from two patients. When investigating outbreaks, one should be aware that plasmid transmission can occur and the possibility of within- and between-patient plasmid variation needs to be considered.

About the Authors:

Joep J. J. M. Stohr

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Software, Validation, Visualization, Writing – original draft

* E-mail: [email protected]

Affiliations Laboratory for Medical Microbiology and Immunology, Elisabeth-TweeSteden Hospital, Tilburg, The Netherlands, Department of Infection Control, Amphia Hospital, Breda, The Netherlands

ORCID logo http://orcid.org/0000-0002-4021-9635

Jaco J. Verweij

Roles Conceptualization, Data curation, Methodology, Supervision, Writing – review & editing

Affiliation: Laboratory for Medical Microbiology and Immunology, Elisabeth-TweeSteden Hospital, Tilburg, The Netherlands

Anton G. M. Buiting

Roles Supervision, Writing – review & editing

Affiliation: Laboratory for Medical Microbiology and Immunology, Elisabeth-TweeSteden Hospital, Tilburg, The Netherlands

John W. A. Rossen

Roles Methodology, Supervision, Validation, Writing – review & editing

Affiliations Department of Medical Microbiology and Infection Prevention, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands, Department of Pathology, University of Utah School of Medicine, Salt Lake City, UT, United States of America

Jan A. J. W. Kluytmans

Roles Conceptualization, Formal analysis, Methodology, Supervision, Validation, Writing – review & editing

Affiliations Department of Infection Control, Amphia Hospital, Breda, The Netherlands, Julius Center for Health Sciences and Primary Care, University Medical Center Utrecht, Utrecht University, Utrecht, The Netherlands, Amphia Academy Infectious Disease Foundation, Amphia Hospital, Breda, The Netherlands

Introduction

Recent years have shown a rise in Klebsiella Pneumoniae carbapenemase (KPC)-producing Klebsiella pneumoniae worldwide [1]. Infections with KPC-producing K. pneumoniae are associated with increased mortality and an increased length of hospital stay [2,3]. Moreover, nosocomial infections and colonization with KPC-producing K. pneumoniae are known to be an important source for its transmission within and between health care facilities [4,5]. Prolonged carriage of KPC-producing K. pneumoniae has been described and several risk factors associated with an increased duration of colonization have been identified [6]. The gene encoding the KPC enzyme in K. pneumoniae (blaKPC) is generally located on large conjugative plasmids which can undergo multiple rearrangements during long-term patient colonization [7,8]. Studies investigating the dynamics of blaKPC- plasmids in K.pneumoniae isolates during colonization only include the KPC-producing (or carbapenem-resistant) isolates [8,9]. However, when only blaKPC containing K.pneumoniae isolates (KPC-KP) are included in the analysis, loss of the KPC enzyme encoding plasmid itself cannot be detected. Moreover, studies on blaKPC-plasmid dynamics within-patients during an outbreak remain limited, especially in countries with a low prevalence of blaKPC-plasmid carriage [10].

Analysing resistance plasmids encoding the KPC enzyme is typically performed using a combination of short- and long-read whole-genome sequencing of an isolate [8,9,11]. Repeat sequences prohibit the complete assembly of the bacterial chromosome and plasmids using short-read sequence data only, resulting in separate contigs of which the origin, plasmid or chromosome, is unknown. [12]. Current automated algorithms aiming to reconstruct plasmids from short-read sequence data are not able to correctly construct large resistance plasmids [12].

In 2017 an outbreak occurred of a KPC-KP in a teaching hospital in Tilburg, the Netherlands. During this period, in all patients colonized or infected with a KPC-producing K. pneumoniae, K. pneumoniae isolates were collected. To investigate the within-patient plasmid and chromosomal dynamics during this outbreak a selection of isolates was sequenced and a plasmid analysis was performed using a hybrid assembly of short- and long-read sequence data.

Method

Klebsiella pneumoniae collection

From 22 October 2017 until 31 December 2017 an outbreak of a KPC-KP occurred in the intensive care unit and surgical ward of a 796-bed teaching hospital in Tilburg, the Netherlands (Fig 1). The outbreak was recognized on the 22nd of October 2017, when a KPC-KP was detected in a urine sample of a patient (patient 2) admitted on the surgical ward. This event followed the earlier repatriation of a patient (the index patient 1) from an Italian hospital on the 11th of September 2017, who was found to be colonized with KPC-KP two days after admission to the intensive care unit (Fig 1). Because carbapenem-resistant Enterobacteriaceae have been practically absent in this hospital so far and because patient 2 was also admitted on the intensive care unit previous to the detection of a KPC-KP in the patient’s urine sample (3–7 October), the finding was considered suspect for nosocomial transmission (Fig 1). An outbreak management team was formed and in both the surgical ward and intensive care unit patient contacts were screened for KPC-KP carriage (Fig 1). During this outbreak, a total of 6 patients were colonized (n = 4) or infected (n = 2) with a KPC-KP. In all patients, for every specimen wherein a K. pneumoniae isolate was obtained, an isolate, regardless of the susceptibility pattern, was collected and stored at -80C° using MicrobankTM. Species identification was performed using the Bruker MALDI BiotyperTM (BD Diagnostics, MD, USA), and antimicrobial susceptibility testing was performed using the PhoenixTM platform (BD Diagnostics, MD, USA) and EUCAST breakpoints v.9.0. [13]. Every specimen, from which a K. pneumoniae isolate was obtained that was measured susceptible to meropenem and/or imipenem, was additionally inoculated on a CHROMagarTM KPC plate (CHROMagar, Paris, France).

[Figure omitted. See PDF.]

Fig 1. Timeline graph with the ward each patient was admitted on previous to and during the outbreak period and the day of the first cultured KPC-KP of each patient detected in the outbreak.

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

Short-read whole-genome sequencing

A selection of K. pneumoniae isolates was sequenced on an Illumina MiSeq using Nextera XT chemistry (Illumina, San Diego, United States) and assembled with SPAdes v. 3.9.1 [14]. The selection was made in a way that: at least each K. pneumoniae isolate with a distinct susceptibility pattern and at least one K. pneumoniae isolate per patient per specimen type was sequenced. A distinct susceptibility pattern was defined as a four-fold difference in minimal inhibitory concentration (MIC) in any of the following antibiotics: amoxicillin-clavulanic acid, ceftriaxone, ceftazidime, meropenem, ciprofloxacin, and gentamicin. Before sequencing isolates were regrown and plated on a CHROMagarTM KPC plate (CHROMagar, Paris, France) when measured resistant to meropenem and on sheep blood agar when measured susceptible to meropenem. Plates were incubated for 18 to 24 hours at 35 to 37C°. The DNA isolation and sequencing protocol are described in the S1 Data. The following quality control criteria for acceptable assemblies were used: coverage: ≥20; number of scaffolds: ≤1000; N50: ≥15.000 bases and maximum scaffold length: ≥50.000 bases.

Short-read whole genome analysis

Whole-genome MLST (wgMLST) (core and accessory genome) was performed for all sequenced isolates using Ridom SeqSphere+, version 4.1.9 (Ridom, Münster, Germany). Species-specific typing schemes were used as described by Kluytmans-van den Bergh et al. [15]. The all-to-all pairwise genetic difference was calculated between the isolates by counting the total number of allele differences in the wgMLST typing scheme and by dividing the total number of allele differences in the wgMLST typing scheme by the total number of shared alleles in the wgMLST typing scheme, ignoring pairwise missing values. The phylogenetic tree was visualized using iTOL v5.5.1 [16]. The genomes of the sequenced isolates were uploaded to the online bioinformatic tools ResFinder v.3.1 and PlasmidFinder v.2.0 (Center for Genomic Epidemiology, Technical University of Denmark, Lingby, Denmark) [17,18]. Acquired resistance genes were called when at least 60% of the length of the best matching gene in the ResFinder database was covered with a sequence identity of at least 90%. Plasmid replicon genes were called when at least 60% of the sequence length of the replicon gene in the PlasmidFinder database was covered with a sequence identity of at least 80%.

Long-read whole genome sequencing

A selection of the isolates was long-read sequenced on a MinION sequencer using the FLO-MIN106D flow cell and the Rapid Barcoding Sequencing Kit SQK RBK004 according to the standard protocol provided by the manufacturer (Oxford Nanopore Technologies, Oxford, United Kingdom). The selection was made in a way that in each patient in which more than one isolate was short-read sequenced all isolates with a unique plasmid replicon content were long-read sequenced. Short-and long-read sequencing was performed from extracted DNA of the same regrown culture.

Hybrid short- and long-read plasmid analysis

A hybrid assembly of long-read and short-read sequence data was performed using Unicycler v.0.8.4 [19]. The genomes created using the hybrid assembly were uploaded to the online bioinformatic tools ResFinder v.3.1 and PlasmidFinder v.2.0 (Center for Genomic Epidemiology, Technical University of Denmark, Lingby, Denmark) [17,18]. Contigs created by the hybrid assembly that were smaller than 1000kb and that contained plasmid replicons were extracted from the assembly graph using BANDAGE v.0.8.1 [20]. All extracted plasmid contigs were annotated using Prokka v. 1.13.3 [21]. A pan-genome was constructed and pairwise comparisons were performed of plasmids between isolates of the same patient using BLAST+ v.2.6.0. (identity cut-off 95%) and Gview v.1.7. via the Gview webserver (https://server.gview.ca/) [22,23].

Accession numbers

Generated raw reads were submitted to the European Nucleotide Archive (ENA) of the European Bioinformatics Institute (EBI) under the study accession number: PRJEB35018 (link to data: https://www.ebi.ac.uk/ena/data/view/PRJEB35018).

Results

Klebsiella pneumoniae collection

During the outbreak period, a total of 35 K. pneumoniae isolates (patient 1: n = 15, patient 2: n = 13, patient 3: n = 4 and patient 4–6: n = 1) with two distinct susceptibility patterns were collected(Fig 2; S1 Table). The two distinct susceptibility patterns were detected in isolates cultured from patient 2 only and were based on differences in MIC for amoxicillin-clavulanic acid, ceftriaxone, ceftazidime, and meropenem (Table 1). In the specimens containing a K. pneumoniae isolate susceptible to meropenem, no growth was detected on the CHROMagarTM KPC plate. Fifteen isolates were selected to be sequenced: patient 1: n = 4, patient 2: n = 7, patient 3–6: n = 1 (Fig 2; Table 1).

[Figure omitted. See PDF.]

Fig 2. Timeline graphic containing the collected isolates per patient during the study period.

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

[Figure omitted. See PDF.]

Table 1. K. pneumoniae isolates used for short-read whole-genome sequencing.

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

Short-read whole genome analysis

Short-read WGS was performed on 15 K. pneumoniae isolates of 6 patients (Fig 2). Despite MIC testing revealing two distinct susceptibility patterns, using wgMLST the maximum number of allele differences detected between the various isolates was 4 (0.09%)(S1 Fig; S2 Table). Moreover, in none of the pairwise comparisons of the sequenced isolates did the number of allele differences exceed the limit of clonal relatedness (smaller or equal to 0.45%) as defined by Kluytmans-van den Bergh et al. (S1 Fig) [15]. The acquired resistance gene content did differ most notably with four isolates not containing a blaKPC gene in the Whole Genome Assembly (WGA) (Table 2), explaining the difference in antimicrobial susceptibility profile seen between the isolates. Moreover, 2 isolates contained a tet(A) gene not detected in any of the other genomes (Table 2). Plasmid replicon content also differed between the isolates: three isolates contained one plasmid replicon gene, three isolates contained two plasmid replicon genes and nine isolates contained three plasmid replicon genes. The difference in plasmid replicon- and acquired resistance gene-content was shown between isolates collected from the same patient (both in isolates from patient 1 and patient 2) (Tables 1 and 2).

[Figure omitted. See PDF.]

Table 2. Acquired resistance gene- and plasmid replicon-content of whole-genome assembly of the sequenced K. pneumoniae isolates.

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

Hybrid short- and long-read plasmid analysis

Four isolates were selected, based on within-patient plasmid replicon content differences, to be long-read sequenced: KP1 and KP3 of patient 1 and KP5 and KP9 of patient 2 (Fig 2). Hybrid assembly of short- and long-read sequences revealed 10 contigs of which, based on size, 6 were assumed to be of plasmid origin (Table 3). Despite the fact that only one IncFII(k) replicon was detected in the short-read whole-genome assembly of isolate KP1, the hybrid assembly revealed that the IncFII(K) plasmid replicon was actually present in two separate plasmids in the KP1 isolate: an IncFII(K) and IncFIB(K) replicon plasmid and an IncFII(K)and IncFIB(pQil) replicon plasmid. Within patient 1, GView BLAST analysis revealed the 3.1 plasmid to be a recombinant plasmid resulting from a recombination event between 18 CDS of plasmid 1.1 (containing the blaKPC gene), 50 CDS present in plasmid 1.1 and plasmid 1.2, 56 CDS of plasmid 1.2 only combined with an introduction of a tet(A) gene containing transposon (Fig 3A). Moreover, it revealed a loss of major parts of the plasmid content between isolates KP1 and KP3 without affecting the isolates susceptibility pattern. In patient 2, the loss of the entire 5.1 blaKPC gene-containing plasmid was observed between isolates KP5 and KP9 (Table 3). Additionally, a 26.786 bp deletion occurred in plasmid 5.2 when compared to plasmid 9.1 resulting in loss of 34 coding sequences among which were the antibiotic resistance genes aac(6')Ib-cr, blaOXA-1, and a catB3-like gene (Fig 3B; Table 3).

[Figure omitted. See PDF.]

Fig 3.

a) Gview BLAST plasmid comparison of all annotated plasmids in patient 1. b) Plasmid comparison of annotated plasmids 5.2 and 9.1 of patient 2. Each arrow represents a coding sequence and not necessarily transcriptional direction; Gene names are depicted as generated by prokka.

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

[Figure omitted. See PDF.]

Table 3. Size, gene-, acquired antimicrobial resistance gene-and plasmid replicon-content of the plasmid contigs created with the hybrid assemblies.

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

Discussion

The plasmid replicon content of the first sequenced blaKPC containing isolate of each patient during the hospital outbreak was similar between the different patients. However, during the outbreak within both patient 1 and patient 2 the plasmid replicon contents highly varied. This variation in plasmid replicon content was partially the result of plasmid loss observed in both patients, leading to a distinct susceptibility pattern in the isolates of one of these patients. A previous study also described plasmid loss during long time colonization in K. pneumoniae [9]. However, the present study also includes isolates with all distinct resistance patterns revealing the loss of a blaKPC containing plasmid within an outbreak setting. The similar replicon content of the plasmids 1.1, 1.2 and 5.1 would suggest that incompatibility between these plasmids [24,25] was the cause of the plasmid loss observed in the isolates of patient 1 and 2. Besides plasmid loss, within-patient 1 acquisition of a tet(A) containing transposon was detected in the plasmid content of the KP3 isolate when compared to the plasmid content of isolate KP1. Thus not only gene loss was revealed but also the acquisition of genetic elements in the plasmid content of the isolates collected during this outbreak. The plasticity of the plasmid content in bacterial isolates observed in this study has been described before both in vitro as in vivo [7,9,26]. However, recent reports also describe plasmids which remain highly stable [27,28]. This suggests that in vivo plasmid stability is likely the result of an interplay between host factors, plasmid content and the different plasmids composing the plasmid content of a bacterial isolate, possibly resulting in either a highly stable or unstable plasmid content.

The hybrid assembly revealed that the 3.1 plasmid was the result of a recombination event between the 1.1 and the 1.2 plasmid occurring in the KP3 isolate only (and possibly the KP4 isolate) and not in the KP5 isolate. These recombination events between different plasmids in the same isolates have also been described in other studies [9,26]. However, this is to the best of our knowledge the first study to describe within patient blaKPC gene-containing plasmid recombination and loss during an hospital outbreak. This recombination event led to two different blaKPC plasmids occurring within one patient and plasmid transmission of these two different plasmids could have occurred during this outbreak. Several studies have already reported blaKPC-plasmid transmission between different isolates during outbreaks [11,29]. Distinguishing outbreak related from non-outbreak related plasmids based on sequence data is essential for using molecular data to confirm blaKPC-containing plasmid transmission in outbreaks. Our findings suggest that when investigating plasmid transmission during outbreaks, the possibility of within-patient plasmid variation needs to be considered. Therefore, it could well be that transmission of blaKPC-containing IncF plasmids within hospital outbreaks cannot be dismissed based on sequence dissimilarity between the different plasmids investigated. Complicating the investigation of plasmid transmission during hospital outbreaks even when the antibiotic susceptibility pattern is not altered.

The present study has some limitations. The first was that in a culture of a specific specimen, colonies of the same morphology were not routinely isolated and stored. Therefore, possible subpopulations were not detected. Despite this, in the cultures in which a K. pneumoniae was isolated that was measured susceptible to carbapenem, no growth was observed on the CHROMagar KPCTM suggesting that no resistant subpopulations were present in these specimens. Based on susceptibility pattern differences, specimen type of isolation, and plasmid replicon content a selection of the isolates were sequenced therefore plasmid variations that did not influence the susceptibility pattern, specimen type of isolation, and plasmid replicon content might go undetected. Moreover, since only 1 isolate was collected in patients 3–6 no conclusions can be drawn regarding longitudinal plasmid variation in these patients.

Concluding, during a hospital outbreak of a blaKPC producing K. pneumoniae isolate plasmid loss of the blaKPC carrying plasmid and plasmid recombination was detected in two patients. When investigating outbreaks wherein plasmid transmission can occur, the possibility of within- and between-patient plasmid variation needs to be considered.

Supporting information

[Figure omitted. See PDF.]

S1 Data. Supplementary materials and methods.

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

(DOCX)

S1 Table. Susceptibility testing results of the collected K. pneumoniae isolates.

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

(XLSX)

S2 Table. Distance matrix of allele differences (% of allele differences) between the isolates sequenced in this study as determined using wgMLST.

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

(XLSX)

S3 Table. Table A response to reviewers.

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

(XLSX)

S1 Fig. Neighbour joining tree based on the wgMLST analysis of the different isolates.

A publicly available genome of a K. pneumoniae sequence type (ST) 256 strain (ATCC® BAA-1705TM) was included in the tree as reference for genetic distance.

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

(EPS)

Citation: Stohr JJJM, Verweij JJ, Buiting AGM, Rossen JWA, Kluytmans JAJW (2020) Within-patient plasmid dynamics in Klebsiella pneumoniae during an outbreak of a carbapenemase-producing Klebsiella pneumoniae. PLoS ONE 15(5): e0233313. https://doi.org/10.1371/journal.pone.0233313

References

1. Munoz-Price LS, Poirel L, Bonomo RA, Schwaber MJ, Daikos GL, Cormican M, et al. Clinical epidemiology of the global expansion of Klebsiella pneumoniae carbapenemases. Lancet Infect Dis [Internet]. 2013 Sep;13(9):785–96. Available from: https://linkinghub.elsevier.com/retrieve/pii/S1473309913701907 pmid:23969216

2. Cassini A, Högberg LD, Plachouras D, Quattrocchi A, Hoxha A, Simonsen GS, et al. Attributable deaths and disability-adjusted life-years caused by infections with antibiotic-resistant bacteria in the EU and the European Economic Area in 2015: a population-level modelling analysis. Lancet Infect Dis. 2019;19(1):56–66. pmid:30409683

3. Schwaber MJ, Klarfeld-Lidji S, Navon-Venezia S, Schwartz D, Leavitt A, Carmeli Y. Predictors of carbapenem-resistant Klebsiella pneumoniae acquisition anions hospitalized adults and effect of acquisition on mortality. Antimicrob Agents Chemother. 2008;52(3):1028–33. pmid:18086836

4. Weterings V, Zhou K, Rossen JW, van Stenis D, Thewessen E, Kluytmans J, et al. An outbreak of colistin-resistant Klebsiella pneumoniae carbapenemase-producing Klebsiella pneumoniae in the Netherlands (July to December 2013), with inter-institutional spread. Eur J Clin Microbiol Infect Dis. 2015;34(8):1647–55. pmid:26067658

5. David S, Reuter S, Harris SR, Glasner C, Feltwell T, Argimon S, et al. Epidemic of carbapenem-resistant Klebsiella pneumoniae in Europe is driven by nosocomial spread. Nat Microbiol [Internet]. 2019 Jul 29; Available from: http://www.nature.com/articles/s41564-019-0492-8

6. Feldman N, Adler A, Molshatzki N, Navon-Venezia S, Khabra E, Cohen D, et al. Gastrointestinal colonization by KPC-producing Klebsiella pneumoniae following hospital discharge: Duration of carriage and risk factors for persistent carriage. Clin Microbiol Infect [Internet]. 2013;19(4):E190–6. Available from: http://dx.doi.org/10.1111/1469-0691.12099 pmid:23331385

7. Hardiman CA, Weingarten RA, Conlan S, Khil P, Dekker J, Mathers A, et al. Horizontal Transfer of Carbapenemase-Encoding Plasmids and. Antimicrob Agents Chemother. 2016;60(8):4910–9. pmid:27270289

8. Ruppé E, Olearo F, Pires D, Baud D, Renzi G, Cherkaoui A, et al. Clonal or not clonal? Investigating hospital outbreaks of KPC-producing Klebsiella pneumoniae with whole-genome sequencing. Clin Microbiol Infect. 2017;23(7):470–5. pmid:28143787

9. Conlan S, Park M, Deming C, Thomas PJ, Young AC, Coleman H, et al. Plasmid Dynamics in KPC-Positive Klebsiella pneumoniae during Long-Term Patient Colonization. MBio [Internet]. 2016 Jul 6;7(3):1–9. Available from: http://mbio.asm.org/lookup/doi/10.1128/mBio.00742-16

10. van den Bunt G, van Pelt W, Hidalgo L, Scharringa J, de Greeff SC, Schürch AC, et al. Prevalence, risk factors and genetic characterisation of extended-spectrum beta-lactamase and carbapenemase-producing Enterobacteriaceae (ESBL-E and CPE): a community-based cross-sectional study, the Netherlands, 2014 to 2016. Euro Surveill. 2019;24(41).

11. Martin J, Phan HTT, Findlay J, Stoesser N, Pankhurst L, Navickaite I, et al. Covert dissemination of carbapenemase-producing Klebsiella pneumoniae (KPC) in a successfully controlled outbreak: Long- and short-read whole-genome sequencing demonstrate multiple genetic modes of transmission. J Antimicrob Chemother. 2017;72(11):3025–34. pmid:28961793

12. Arredondo-Alonso S, Willems RJ, van Schaik W, Schürch AC. On the (im)possibility of reconstructing plasmids from whole-genome short-read sequencing data. Microb Genomics [Internet]. 2017 Oct 1;3(10). Available from: http://www.microbiologyresearch.org/content/journal/mgen/10.1099/mgen.0.000128.v1

13. The European Committee on Antimicrobial Susceptibility Testing. Breakpoint tables for interpretation of MICs and zone diameters. Version 8.1, 2018. http://www.eucast.org.

14. Nurk S, Bankevich A, Antipov D, Gurevich A, Korobeynikov A, Lapidus A, et al. Assembling Genomes and Mini-metagenomes from Highly Chimeric Reads. In 2013. p. 158–70. Available from: http://link.springer.com/10.1007/978-3-642-37195-0_13

15. Kluytmans-van den Bergh MFQ, Rossen JWA, Bruijning-Verhagen PCJ, Bonten MJM, Friedrich AW, Vandenbroucke-Grauls CMJE, et al. Whole-Genome Multilocus Sequence Typing of Extended-Spectrum-Beta-Lactamase-Producing Enterobacteriaceae. J Clin Microbiol [Internet]. 2016;54(12):2919–27. Available from: http://www.ncbi.nlm.nih.gov/pubmed/27629900%0Ahttp://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=PMC5121380 pmid:27629900

16. Letunic I, Bork P. Interactive Tree Of Life (iTOL) v4: recent updates and new developments. Nucleic Acids Res [Internet]. 2019 Jul 2;47(W1):W256–9. Available from: https://academic.oup.com/nar/article/47/W1/W256/5424068 pmid:30931475

17. Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, et al. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother. 2012;67(11):2640–4. pmid:22782487

18. Carattoli A, Zankari E, Garciá-Fernández A, Larsen MV, Lund O, Villa L, et al. In Silico detection and typing of plasmids using plasmidfinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother. 2014;58(7):3895–903. pmid:24777092

19. Wick RR, Judd LM, Gorrie CL, Holt KE. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol. 2017;13(6):1–22.

20. Wick RR, Schultz MB, Zobel J, Holt KE. Bandage: Interactive visualization of de novo genome assemblies. Bioinformatics. 2015;31(20):3350–2. pmid:26099265

21. Seemann T. Prokka: Rapid prokaryotic genome annotation. Bioinformatics. 2014;30(14):2068–9. pmid:24642063

22. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol [Internet]. 1990 Oct;215(3):403–10. Available from: https://linkinghub.elsevier.com/retrieve/pii/S0022283605803602 pmid:2231712

23. Petkau A, Stuart-Edwards M, Stothard P, van Domselaar G. Interactive microbial genome visualization with GView. Bioinformatics. 2010;26(24):3125–6. pmid:20956244

24. Villa L, García-Fernández A, Fortini D, Carattoli A. Replicon sequence typing of IncF plasmids carrying virulence and resistance determinants. J Antimicrob Chemother. 2010;65(12):2518–29. pmid:20935300

25. Datta N, Hedges RW. Compatibility Groups among fi− R Factors. Nature [Internet]. 1971 Nov;234(5326):222–3. Available from: http://www.nature.com/articles/234222a0 pmid:5002028

26. Xie M, Li R, Liu Z, Chan EWC, Chen S. Recombination of plasmids in a carbapenem-resistant NDM-5-producing clinical Escherichia coli isolate. J Antimicrob Chemother. 2018;73(5):1230–4. pmid:29373691

27. Löhr IH, Hülter N, Bernhoff E, Johnsen PJ, Sundsfjord A, Naseer U. Persistence of a pKPN3-like CTX-M-15-encoding IncFIIK plasmid in a Klebsiella pneumonia ST17 host during two years of intestinal colonization. PLoS One. 2015;10(3):1–16.

28. Roer L, Overballe-Petersen S, Hansen F, Johannesen TB, Stegger M, Bortolaia V, et al. ST131 fimH22 Escherichia coli isolate with a blaCMY-2/IncI1/ST12 plasmid obtained from a patient with bloodstream infection: Highly similar to E. coli isolates of broiler origin. J Antimicrob Chemother. 2019;74(3):557–60. pmid:30496481

29. Conlan S, Thomas PJ, Deming C, Park M, Lau AF, Dekker JP, et al. Single-molecule sequencing to track plasmid diversity of hospital-associated carbapenemase-producing Enterobacteriaceae. Sci Transl Med [Internet]. 2014 Sep 17;6(254):254ra126–254ra126. Available from: http://stm.sciencemag.org/cgi/doi/10.1126/scitranslmed.3009845 pmid:25232178