ORIGINAL RESEARCH article
Sec. Antimicrobials, Resistance and Chemotherapy
Nanopore MinION Sequencing Reveals Possible Transfer of blaKPC–2 Plasmid Across Bacterial Species in Two Healthcare Facilities
- 1Wadsworth Center, New York State Department of Health, Albany, NY, United States
- 2Maine Center for Disease Control and Prevention, Department of Health and Human Services, Augusta, ME, United States
Carbapenemase-producing Enterobacteriaceae are a major threat to global public health. Klebsiella pneumoniae carbapenemase (KPC) is the most commonly identified carbapenemase in the United States and is frequently found on mobile genetic elements including plasmids, which can be horizontally transmitted between bacteria of the same or different species. Here we describe the results of an epidemiological investigation of KPC-producing bacteria at two healthcare facilities. Using a combination of short-read and long-read whole-genome sequencing, we identified an identical 44 kilobase plasmid carrying the blaKPC–2 gene in four bacterial isolates belonging to three different species (Citrobacter freundii, Klebsiella pneumoniae, and Escherichia coli). The isolates in this investigation were collected from patients who were epidemiologically linked in a region in which KPC was uncommon, suggesting that the antibiotic resistance plasmid was transmitted between these bacterial species. This investigation highlights the importance of long-read sequencing in investigating the relatedness of bacterial plasmids, and in elucidating potential plasmid-mediated outbreaks caused by antibiotic resistant bacteria.
Carbapenem-resistant Enterobacteriaceae (CRE) are an urgent global health threat, and have been categorized by the World Health Organization (World Health Organization [WHO], 2017) and the United States Centers for Disease Control and Prevention (CDC) (Centers for Disease Control and Prevention [CDC], 2019a) as top priorities for research, drug discovery, surveillance, and control. CRE that produce carbapenemases are particularly concerning epidemiologically because carbapenemase genes can be transferred among bacteria via mobile genetic elements, including plasmids (Bonomo et al., 2017). In the United States, the most commonly identified carbapenemase is Klebsiella pneumoniae carbapenemase (KPC), which has become endemic in parts of the country since it was first described in 1996 (Woodworth et al., 2018; Castanheira et al., 2019).
Though first identified in K. pneumoniae, the gene encoding KPC has been detected across multiple genera of gram-negative bacteria (Woodworth et al., 2018; Brandt et al., 2019). Over twenty variants of blaKPC have been described, of which blaKPC–2 and blaKPC–3 are the most frequently detected (Castanheira et al., 2019). blaKPC genes are commonly carried inside transposons, in particular Tn4401, a 10kb self-mobilizing transposon in the Tn3 family (Cuzon et al., 2011; Partridge et al., 2018). Tn4401 and blaKPC are often encoded on plasmids (Brandt et al., 2019), but can also be found integrated into the bacterial chromosome (Mathers et al., 2017).
Though dissemination of KPC-producing bacteria can result from epidemic spread of clonal lineages (Kitchel et al., 2009; Hargreaves et al., 2015), horizontal transfer of KPC-encoding plasmids across unrelated bacteria of the same or different species has also been described (Conlan et al., 2014, 2019; Sheppard et al., 2016; Li et al., 2018; Evans et al., 2020). Detecting such instances of plasmid transfer requires the use of whole-genome sequencing (WGS) to characterize the genomic context of antibiotic resistance genes (Boolchandani et al., 2019). As short-read sequencing has been shown to be limited in its ability to resolve the highly repetitive regions common in plasmids (Arredondo-Alonso et al., 2017), combined short-read and long-read WGS of antibiotic resistant bacteria has become an increasingly common method to better elucidate plasmid structures, and to detect plasmid-mediated outbreaks (Lemon et al., 2017; Conlan et al., 2019; Decano et al., 2019; Van Dorp et al., 2019; Wyres et al., 2019). In this study, we describe a possible case of plasmid transfer detected by hybrid analysis of short read Illumina MiSeq and long-read Oxford Nanopore Technologies (ONT) MinION sequencing, in which an identical plasmid carrying blaKPC–2 was identified across three bacterial species isolated from epidemiologically linked patients in two healthcare facilities.
Materials and Methods
The Maine Center for Disease Control and Prevention utilizes a combination of required and voluntary reporting, along with isolate submission to the state public health laboratory (SPHL), to identify CRE that produce carbapenemases. In 2018, the SPHL identified KPC-producing bacteria in two unique clinical specimens in the same month. In response, and in accordance with guidance from the Centers for Disease Control and Prevention [CDC], 2019b, investigators arranged for colonization screens of epidemiologically linked patients who had been admitted at two associated healthcare facilities by rectal swab. Specimens from patients who had been discharged home were collected by walk-in clinics at one associated facility and by patients’ primary care providers.
The Wadsworth Center, New York State Department of Health’s public health laboratory and the Northeast Regional Laboratory for the CDC funded Antimicrobial Resistance Laboratory Network (ARLN), analyzed all rectal swab specimens from colonization screens with the Cepheid Xpert® Carba-R test. This test is a real-time PCR assay for rapid detection and differentiation of five genes (blaKPC, blaVIM, blaOXA–48, blaIMP–1, and blaNDM) responsible for carbapenem resistance.
This screening identified three additional patients with KPC-producing bacteria representing three species: Citrobacter freundii, Klebsiella pneumoniae, and Escherichia coli. Traditionally, outbreak definitions include a requirement for matching organisms; however, as carbapenemase genes are transferrable among bacterial species and epidemiological links were present, further investigation of possible plasmid transfer was undertaken.
Isolates were recovered from rectal swab specimens that were positive for blaKPC by the Cepheid Xpert® Carba-R test by streaking on MacConkey agar. For clinical and colonization isolates, identification to species was confirmed by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). Antibiotic susceptibility testing was performed by broth microdilution with the Thermo Scientific Sensititre Aris 2X using the GNX2F panel, and by ETEST® (bioMérieux) for ceftazidime-avibactam (all specimens) and ertapenem, meropenem, and imipenem (Isolate 5 only). In addition, disk diffusion was used to test the susceptibility of Isolate 5 to doripenem, ertapenem, imipenem, and meropenem. Antibiotic susceptibilities were interpreted using CLSI M100-ED29 breakpoints for Enterobacteriales (Clinical and Laboratory Standards Institute [CLSI], 2019). The modified carbapenem inactivation method (mCIM) (Clinical and Laboratory Standards Institute [CLSI], 2019) was used to detect carbapenemase production. Molecular characterization of resistance mechanisms was performed using New York State Clinical Laboratory Evaluation Program (CLEP)-approved multiplex real-time PCR assays to detect blaKPC and blaNDM (developed at the Wadsworth Center), as well as blaVIM, blaIMP (all variants), and blaOXA–48–like genes (developed at the CDC).
Illumina WGS and Analysis
Genomic DNA was extracted from isolates using the DNeasy Blood & Tissue Kit on a QIAcube (QIAGEN). DNA was quantified using the Qubit dsDNA BR assay system. Sequence libraries were prepared using the Nextera XT DNA Sample Preparation Kit and sequenced on the Illumina MiSeq system at the Wadsworth Center Applied Genomic Technologies Core.
Raw Illumina reads were processed with Trimmomatic v0.38 (Bolger et al., 2014) and bacterial species identification was confirmed in-silico using Kraken v1.0 (Wood and Salzberg, 2014) with the MiniKraken 8GB database; paired, 250 bp reads were then de novo assembled into contigs with SPAdes v3.12.0 (Bankevich et al., 2012). Assembly quality was assessed using quantitative measurements, including BUSCO v3.1.0 (Simão et al., 2015; Waterhouse et al., 2017), prior to multilocus sequence typing analysis (MLST) with mlst v2.16.21 and AR gene identification with ABRicate v0.8.13.2 Final analysis of the antibiotic resistance genes in the genome assembly compared gene identification between the NCBI Bacterial Antimicrobial Resistance Reference Gene Database (NCBI) (Feldgarden et al., 2019), ResFinder (Zankari et al., 2012), and Comprehensive Antibiotic Resistance Database (CARD) (Jia et al., 2016) databases to determine the best matches.
For mutation event analysis, Mash v1.1 (Ondov et al., 2016) was used to select the best possible reference genome from a number of candidates prior to mapping with BWA/Samtools [v0.7.17 and v1.9, respectively (Li and Durbin, 2009, 2010; Li et al., 2009)]. Mutation events (ME), defined as the number of single nucleotide polymorphisms (SNPs) and insertion/deletion events, were called using FreeBayes v1.0.2 (Garrison and Marth, 2012). A ME matrix was constructed from all isolates within a cluster by pairwise comparison of all reference-aligned sequences to count MEs while ignoring ambiguous or missing bases.
ONT MinION WGS
For genomic DNA extraction, all isolates were sub-cultured from frozen stocks twice on blood agar plates. Colonies were resuspended in 2 ml sterile water to 4 McFarland concentration and harvested as pellets. High molecular weight genomic DNA was extracted from Isolates 1, 3, 4, and 5 using Genomic-tip 20/G (QIAGEN) and Genomic DNA buffers (QIAGEN). Genomic DNA from Isolate 2 was extracted using the Nanobind CBB Big DNA kit (Circulomics). The protocol for genomic DNA isolation for gram negative bacteria was followed for both methods as suggested by the respective manufacturer. The genomic DNA was quantified using a Qubit fluorometer (ThermoFisher Scientific). Quality of the genomic DNA was assessed using the TapeStation (Agilent).
MinION sequencing libraries were prepared from 1.5 μg of input DNA. Genomic DNA was sheared in 50 μl total volume in Covaris G tubes using an Eppendorf 5425 centrifuge at 6000 rpm. The sequencing library was prepared according to manufacturer’s instructions (Oxford Nanopore) and multiplexed using 1D Native barcoding kits (EXP-NBD104, EXP-NBD114) followed by ligation and sequencing kit (SQK-LSK109). The library was loaded on a SpotON flowcell R9.4.1 FLO-MIN106 and sequenced for 72 h on the MinION device. The fast5 data from MinKNOW was converted to fastq format using the Guppy basecaller in fast mode on a MinIT device (Oxford Nanopore, United Kingdom). The fastq reads were demultiplexed using qcat v18.104.22.168
Hybrid Genome Assembly and Annotation
MinION reads were quality filtered using filtlong v0.2.04 with a minimum read length of 1000 and the target number of bases set to 500,000,000 (to provide approximately 100X coverage of the target species’ genomes). Paired-end MiSeq reads were trimmed to remove adapters and low-quality ends (<q10) using trim_galore v0.6.4.5
Genomes were assembled using two methods: Unicycler v0.4.8 (Wick et al., 2017) hybrid assembly using the filtered MinION reads and trimmed Illumina reads, with default settings; and Flye v2.6 (Kolmogorov et al., 2019) using the filtered MinION reads, with plasmids and meta options enabled. Flye assemblies were polished with racon v1.4.7 (Vaser et al., 2017) (-m 8 -x -6 -g -8 -w 500 –no-trimming) using MinION reads mapped to the assembly with minimap2 v2.11 (Li, 2018), followed by medaka v0.8.16 and two rounds of pilon v1.23 (Walker et al., 2014) using Illumina reads mapped to the assembly with bwa-mem (Li, 2013). Discrepancies between the assemblies were assessed by comparing the percentage of Illumina reads that aligned to the assembly using bwa-mem, and the percentage of filtered MinION reads that aligned to the assembly using minimap2. In addition, Assembly Likelihood Evaluation (ALE) scores (Clark et al., 2013) were calculated using both the Illumina and MinION alignment files, and Nanovar v0.1.2 (Tham et al., 2019) was used to identify structural discrepancies between MinION reads and the assemblies. All Illumina and MinION sequences were deposited in the NCBI sequence read archive (SRA), and final assemblies in GenBank (BioProject ID PRJNA636827).
Plasflow v1.1 (Krawczyk et al., 2018) was used to classify contigs as plasmids or chromosomes. Final genomes were annotated with prokka v1.14.0 (Seemann, 2014). Abricate v0.9.8 (see text footnote 2) detected antimicrobial resistance genes using the CARD (Jia et al., 2016), ResFinder (Zankari et al., 2012), and NCBI (Feldgarden et al., 2019) databases, and plasmid replicon genes using the PlasmidFinder (Carattoli et al., 2014) database. Identified plasmids were queried against the PLSDB database v2019_10_07 (Galata et al., 2018) using mash (-S 42 -k 21 -s 1000) (Ondov et al., 2016) to identify related plasmids.
Following the identification of two patients with KPC-producing bacteria associated with Facility A (Patient 1, Isolate 1: C. freundii; Patient 2, Isolate 2: K. pneumoniae), a retrospective investigation including colonization screening of patients who had been roommates with or who had overlapped on the same unit for three or more days with either of the two source patients was performed. Rectal swabs from two additional patients who had been on the same unit at Facility A tested positive for KPC-producing bacteria after these patients had been discharged. Isolate recovery from one rectal swab specimen was unsuccessful (Patient X), but a KPC-producing K. pneumoniae (Isolate 3 from Patient 3) was recovered from the second specimen. Following the initial investigation, routine point prevalence surveys were done at Facility A on the affected unit to identify any additional KPC-producing bacteria. 3 months after the initial investigation, Isolate 4 (C. freundii) was recovered from a rectal swab collected from Patient 4, who had recently been admitted to Facility A (Figure 1).
Figure 1. Timeline of KPC-producing bacterial isolates collected during epidemiological investigations at Facilities A and B. The admission dates of Patients 1-5 and Patient X at Facility A (blue) and Facility B (orange) are shown as arrows. The bacterial species, specimen type, and date of specimen collection are shown for each bacterial isolate.
Prior to the initiation of this investigation, Patients 2, X, and 3 had been transferred from Facility A to Facility B. Patients 2 and X were placed in isolation at the time of admission to Facility B, but Patient 3 was assigned to a double occupancy room with a shared bathroom, sharing the room with Patient 5 for 5 days. After discharge from Facility B, a rectal swab obtained from Patient 5 yielded Isolate 5 (E. coli). All other rectal swabs collected from patients overlapping on the same unit for three or more days at Facility B tested negative for KPC-producing bacteria.
Isolates 1–5 differed in their antibiotic susceptibility profiles (Table 1), but all isolates were resistant to aztreonam and ticarcillin/clavulanic acid; Isolates 3 and 4 were non-susceptible to all beta-lactam antibiotics tested. Isolate 5 was susceptible to all carbapenems by broth microdilution. By disk diffusion, Isolate 5 had intermediate susceptibility to ertapenem (21 mm) and imipenem (22.5 mm), and was susceptible to doripenem (23 mm) and meropenem (23 mm). By gradient diffusion, Isolate 5 was susceptible to ertapenem (MIC 0.25 μg/ml), meropenem (MIC 0.19 μg/ml), and imipenem (MIC 1.0 μg/ml). All isolates tested positive for carbapenemase production by mCIM except for Isolate 4, which was not tested. Real-time PCR detected only blaKPC in all isolates.
By Illumina WGS, the two C. freundii isolates (Isolates 1 and 4) were not able to be assigned to a sequence type, but had distinct alleles at six of seven MLST loci, and therefore were determined to be unrelated. The two K. pneumoniae isolates (Isolates 2 and 3) were both ST37 and were closely related, differing by only seven mutation events. Between five and twenty antibiotic resistance genes were identified from each of the five isolates (Supplementary Table 1). Two of these genes, both encoding beta-lactamases (blaKPC–2 and blaTEM–1) were identified at >90% coverage and identity from all isolates.
Hybrid genome assembly methods combining ONT MinION and Illumina reads were used to determine the location of blaKPC–2 within the genomes of the five isolates. Complete genome assemblies consisting of only circular contigs were generated for all isolates with the exception of Isolate 2, which consisted of four circular contigs and a single small (8.5 kilobase) linear contig (Supplementary Methods). Assemblies for all five isolates consisted of closed chromosomes and between two and six additional small contigs. Details on assembly evaluation and rationales for the final choice of assembly method for each isolate are provided (Supplementary Methods). Each isolate harbored at least one plasmid that carried one or more antibiotic resistance gene (Table 2).
Isolate 1 (C. freundii) carried copies of the blaKPC–2 gene on two unique plasmids: a 43,621 bp plasmid that also carried blaTEM–1 (p1C44), and a 73,366 bp plasmid that did not carry other known resistance genes (p1C73). p1C44 is 100% identical across its entire length to the plasmid pKPC_UVA01 (Genbank accession no. CP017937.1), first described in a K. pneumoniae isolated from an abdominal abscess at the University of Virginia Health System in 2007 (Mathers et al., 2015).
Isolate 2 (K. pneumoniae), Isolate 3 (K. pneumoniae), and Isolate 5 (E. coli) also carried the blaKPC–2 gene on a plasmid identical (p5E44) or nearly identical (p2K44 and p3K44) to pKPC_UVA01 (Figure 2). Compared with pKPC_UVA01, there were two single base indels in p2K44 (a deleted thymine corresponding to position 12,990 of pKPC_UVA01 and an inserted guanine corresponding to position 30,639 of pKPC_UVA01, both located in five-nucleotide homopolymers). Similarly, there was a deleted cytosine in p3K44 in a five nucleotide homopolymer corresponding to position 13,910 of pKPC_UVA01. As the five plasmids were otherwise identical across their entire lengths, and as these indels all occurred in homopolymers, the length of which is known to be commonly mis-identified by MinION sequencers (Wick et al., 2019), they likely represent sequencing, assembly, or polishing errors and not true differences in the plasmid sequences. This is supported by mapping of Illumina reads from Isolates 2 and 3 to pKPC_UVA01, which shows that the Illumina reads do not support the deletions in p2K44 or p3K44, though they do support the insertion in p2K44 (Supplementary Methods).
Figure 2. Identity of plasmids from Isolates 1, 2, 3, and 5 to pKPC_UVA01. Green, brown, purple, and blue colored circles depict the percent identity of the assembled sequence of each plasmid to the sequence of pKPC_UVA01 based on BLASTN alignment (Camacho et al., 2009). Sequences of all four plasmids were 100% identical to pKPC_UVA01 across their entire length, with the exception of p2K44, which had two single nucleotide indels compared to pKPC_UVA01, and p3K44, which had a single nucleotide deletion compared to pKPC_UVA01. The gene coding sequences (CDS) of pKPC_UVA01 are shown in gray, with resistance genes (blaTEM–1 and blaKPC–2) colored red. The location of Tn4401b is shown in yellow. Image created using the Blast Ring Image Generator (BRIG) (Alikhan et al., 2011).
Unlike the other isolates, Isolate 4 (C. freundii) carried blaKPC–2 on a 12,158bp Col440I plasmid (p4C12). Like pKPC_UVA01, p1C73 and p4C12 carried blaKPC–2 within the transposon Tn4401b, but the plasmids did not otherwise share sequence identity (Figure 3). A PLSDB query found that p4C12 was most similar (mash distance 0.006, 781 of 1000 shared hashes) to a 9,803bp blaKPC–3-containing plasmid isolated from K. pneumoniae in Spain (GenBank accession no. NC_019151.1). No plasmids similar to p1C73 were identified (lowest mash distance 0.04, 291 of 1000 shared hashes for GenBank accession no. NZ_CP039300.1).
Figure 3. Annotated gene maps of unique identified KPC-2 plasmids. BLASTN hits from pairwise comparisons of the plasmids are shown in pink; all plasmids contain the 10 kb transposon Tn4401b (>99% sequence identity for both comparisons). Coding sequences (CDS) encoding resistance genes are colored in red, transposases/integrases/recombinases in purple, other identified genes in blue, and hypothetical proteins in gray. Image generated using the genoplotR package (Guy et al., 2010).
In addition to sharing high chromosomal genomic similarity and near-identical blaKPC–2 plasmids, Isolates 2 and 3 also shared similar 157kb IncFIA(HI1) plasmids carrying eleven antibiotic resistance genes (Table 2). Despite carrying identical replicon and resistance genes, these two plasmids were not identical across their entire sequence (mash distance 0.004, 836 of 1000 shared hashes).
Using a combination of short and long-read WGS, we identified identical plasmids carrying blaKPC–2 in four isolates of three bacterial genera recovered from epidemiologically linked patients associated with two healthcare facilities. The combination of these epidemiological and genetic findings strongly supports the transfer of this antibiotic resistance plasmid among different bacterial species in these facilities. Additionally, K. pneumoniae Isolates 2 and 3 shared high chromosomal genetic similarity, differing by only seven mutation events across their genomes. Thus, it is likely that in this case, the entire bacterial organism, including the blaKPC–2-carrying plasmid, was transferred between patients.
Isolate 4 (C. freundii) was identified during a routine point prevalence study at Facility A, 3 months after the outbreak. Sequencing analysis determined that this isolate was genetically unrelated to the other C. freundii isolate from this investigation (Isolate 1), and that the plasmid harboring blaKPC–2 in Isolate 4 was unrelated to the plasmids harboring the gene in Isolates 1, 2, 3, and 5. Further investigation of Isolate 4 identified that the patient was a resident of Massachusetts, a state that had previously identified cases of KPC-producing bacteria (Centers for Disease Control and Prevention [CDC], 2019c), and that this colonization was likely present on admission.
Upon identification of KPC-producing bacteria through colonization screening, Facility A initiated weekly rectal swabs for all patients on the affected unit. No additional KPC-producing bacteria were identified. One year later, there has been no other healthcare-onset KPC-producing bacteria identified at either Facility A or Facility B.
Prevalence of carbapenemase-producing CRE was low in Maine in 2018. At that time, most healthcare facilities in the state did not routinely conduct active surveillance cultures to identify colonized patients upon admission (Centers for Disease Control and Prevention [CDC], 2015). Facility A did note a housekeeping staffing shortage during 2018. Unrecognized colonization and missed opportunities in environmental cleaning may have played a role in the transmission of gastrointestinal flora from the source patient to other patients on the same unit in Facility A. Unrecognized colonization and a shared bathroom may have led to transmission from the source patient to the roommate at Facility B. It has been shown that toilet flushing generates aerosolized bacteria that can land on nearby surfaces or drift in air currents to land on surfaces further away, which can contribute to the direct and/or indirect transmission of gastrointestinal flora (Barker and Jones, 2005; Johnson D. L. et al., 2013; Johnson D. et al., 2013). As environmental screening was not done in the course of this investigation, it is impossible to determine if horizontal transfer of the plasmid harboring blaKPC–2 may have occurred in the environment at Facility A and/or B, or within a patient, or both.
While the epidemiological and genetic data from this investigation support the transmission of the blaKPC–2-carrying plasmid among three species of bacteria, it is also possible that the bacteria independently acquired this resistance plasmid. The identified plasmid is 100% identical to a plasmid first sequenced from K. pneumoniae more than 10 years earlier in Virginia (Mathers et al., 2015). This plasmid (pKPC_UVA01) has also been subsequently identified at 100% identity in other bacterial species (Kluyvera intermedia, C. freundii) from the same hospital in which it was first isolated (Sheppard et al., 2016; Barry et al., 2019), and highly similar plasmids have been identified in Enterobacter species from New York, Michigan, Maryland, Illinois, and Florida (Chavda et al., 2016). This suggests that, despite the fact that pKPC_UVA01 has been shown to have a relatively low conjugation efficiency into K. pneumoniae and E. coli in vitro (Hardiman et al., 2016), it may be widespread in bacterial populations. More long-read sequencing of KPC-producing bacterial populations globally is necessary to gain more insight into the phylogenetic diversity of pKPC_UVA01 and related plasmids.
The identification of three unique plasmids carrying blaKPC–2 across a small sample of KPC-producing bacteria in patients associated with two healthcare facilities highlights the diversity of the plasmid contexts of this resistance gene. The location of blaKPC genes in transposons such as Tn4401b, which are themselves located on plasmids, allows for a high amount of mobility of these genes among bacterial species, patients, and healthcare facilities (Sheppard et al., 2016; Martin et al., 2017; Stoesser et al., 2017; Brandt et al., 2019; Mathers et al., 2019; Hendrickx et al., 2020). The location of blaKPC genes in highly conserved transposons such as Tn4401b also highlights the necessity of using long-read sequencing to differentiate between these plasmids. Genome assembly methods based on short reads alone may not be able to assemble plasmids to the extent necessary to differentiate between plasmids that contain identical transposons on very different backbones (Figure 3). Similarly, the identification of two unique plasmids carrying blaKPC–2 within Tn4401b in Isolate 1 would have been impossible without the use of long-read sequencing.
Long-read and hybrid assembly of bacterial genomes is still an area of active development (De Maio et al., 2019; Wick and Holt, 2019). The discrepancies between assemblies produced by different methods in the current study, including inconsistencies in the genome assembly of Isolate 1 using Flye and in the genome assemblies of Isolates 3 and 5 using Unicycler (Supplementary Methods), support the use of more than one assembly method as well as thorough inspection and evaluation of genome assemblies.
Despite carrying a plasmid encoding blaKPC–2, and testing positive for carbapenemase production by mCIM, Isolate 5 was not carbapenem resistant by broth microdilution (Table 1) or gradient diffusion, though it showed intermediate susceptibility to ertapenem and imipenem by disk diffusion. Carbapenemase producing bacteria with low carbapenem MICs have been previously reported (Fattouh et al., 2016; Tamma et al., 2016), highlighting the importance of utilizing multiple molecular and phenotypic methods to screen clinical isolates for carbapenemase production.
The ability to identify identical plasmids across bacterial species was instrumental in defining this event as an outbreak. Genomic evidence to support the epidemiological suspicion of an outbreak was of great benefit in acceptance by the facility that an outbreak occurred and for the promotion of infection control and prevention activities needed to respond to these novel organisms statewide.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, PRJNA636827.
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
CP wrote the manuscript with contributions from RO, ES, NS, and WH. RO conducted epidemiological investigation and coordinated sample collection. KM, EN, and KM supervised microbiological laboratory work. NS conducted long-read sequencing and base-calling. CP, ES, and WH conducted data analysis with guidance from PL and EL-N. All authors read and approved the final manuscript.
This work was supported by the New York State Department of Health, Cooperative Agreement Number NU50CK000516 funded by the Centers for Disease Control and Prevention (CDC), and Cooperative Agreement U60OE000103 funded by CDC through the Association of Public Health Laboratories. Its contents are solely the responsibility of the authors and do not necessarily represent the official views of the CDC, the Department of Health and Human Services, or the Association of Public Health Laboratories.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors would like to acknowledge Shannon Morris, Kailee Cummings, Janine Bodnar, Christine Wagner, and Dr. Jennifer Liao for their contributions to this investigation.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.02007/full#supplementary-material
- ^ https://github.com/tseemann/mlst
- ^ https://github.com/tseemann/abricate
- ^ https://github.com/nanoporetech/qcat
- ^ https://github.com/rrwick/Filtlong
- ^ https://github.com/FelixKrueger/TrimGalore
- ^ https://github.com/nanoporetech/medaka
Arredondo-Alonso, S., Willems, R. J., van Schaik, W., and Schürch, A. C. (2017). On the (im)possibility of reconstructing plasmids from whole-genome short-read sequencing data. Microb. Genom. 3:e000128. doi: 10.1099/mgen.0.000128
Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021
Barker, J., and Jones, M. V. (2005). The potential spread of infection caused by aerosol contamination of surfaces after flushing a domestic toilet. J. Appl. Microbiol. 99, 339–347. doi: 10.1111/j.1365-2672.2005.02610.x
Barry, K. E., Wailan, A. M., Sheppard, A. E., Crook, D., Vegesana, K., Stoesser, N., et al. (2019). Don’t overlook the little guy: an evaluation of the frequency of small plasmids co-conjugating with larger carbapenemase gene containing plasmids. Plasmid 103, 1–8. doi: 10.1016/j.plasmid.2019.03.005
Bonomo, R. A., Burd, E. M., Conly, J., Limbago, B. M., Poirel, L., Segre, J. A., et al. (2017). Carbapenemase-producing organisms: a global scourge. Clin. Infect. Dis. 66, 1290–1297. doi: 10.1093/cid/cix893
Brandt, C., Viehweger, A., Singh, A., Pletz, M. W., Wibberg, D., Kalinowski, J., et al. (2019). Assessing genetic diversity and similarity of 435 KPC-carrying plasmids. Sci. Rep. 9:11223. doi: 10.1038/s41598-019-47758-47755
Carattoli, A., Zankari, E., García-Fernández, A., Voldby Larsen, M., Lund, O., Villa, L., et al. (2014). In silico detection and typing of plasmids using plasmidfinder and plasmid multilocus sequence typing. Antimicrob. Agents Chemother. 58, 3895–3903. doi: 10.1128/AAC.02412-14
Castanheira, M., Deshpande, L. M., Mendes, R. E., Canton, R., Sader, H. S., and Jones, R. N. (2019). Variations in the occurrence of resistance phenotypes and carbapenemase genes among Enterobacteriaceae isolates in 20 Years of the SENTRY antimicrobial surveillance program. Open Forum Infect. Dis. 6, S23–S33. doi: 10.1093/ofid/ofy347
Centers for Disease Control and Prevention [CDC] (2015). Facility Guidance for Control of Carbapenem-Resistant Enterobacteriaceae (CRE). Washington, DC: Department of Health and Human Services. Available online at: https://www.cdc.gov/hai/pdfs/cre/CRE-guidance-508.pdf (accessed June 1, 2020).
Centers for Disease Control and Prevention [CDC] (2019a). Antibiotic Resistance Threats in the United States, 2019. Washington, DC: Department of Health and Human Services. Available online at: www.cdc.gov/DrugResistance/Biggest-Threats.html (accessed June 1, 2020).
Centers for Disease Control and Prevention [CDC] (2019b). Interim Guidance for a Public Health Response to Contain Novel or Targeted Multidrug-resistant Organisms (MDROs). Washington, DC: Department of Health and Human Services. Available online at: https://www.cdc.gov/hai/containment/guidelines.html (accessed June 1, 2020).
Centers for Disease Control and Prevention [CDC] (2019c). Tracking CRE in United States. Washington, DC: Department of Health and Human Services. Available online at: https://www.cdc.gov/hai/organisms/cre/trackingcre.html (accessed June 1, 2020).
Chavda, K. D., Chen, L., Fouts, D. E., Sutton, G., Brinkac, L., Jenkins, S. G., et al. (2016). Comprehensive genome analysis of carbapenemase-producing Enterobacter spp.: new insights into phylogeny, population structure, and resistance mechanisms. mBio 7, 1–16. doi: 10.1128/mBio.02093-2016
Clark, S. C., Egan, R., Frazier, P. I., and Wang, Z. (2013). ALE: a generic assembly likelihood evaluation framework for assessing the accuracy of genome and metagenome assemblies. Bioinformatics 29, 435–443. doi: 10.1093/bioinformatics/bts723
Conlan, S., Lau, A. F., Deming, C., Spalding, C. D., Lee-Lin, S., Thomas, P. J., et al. (2019). Plasmid Dissemination and selection of a multidrug-resistant Klebsiella pneumoniae strain during transplant-associated antibiotic therapy. mBio 10, e00652-19. doi: 10.1128/mbio.00652-619
Conlan, S., Thomas, P. J., Deming, C., Park, M., Lau, A. F., Dekker, J. P., et al. (2014). Single-molecule sequencing to track plasmid diversity of hospital-associated carbapenemase-producing Enterobacteriaceae. Sci. Transl. Med. 6:254ra126. doi: 10.1126/scitranslmed.3009845
Cuzon, G., Naas, T., and Nordmann, P. (2011). Functional characterization of Tn4401, a Tn3-based transposon involved in blaKPC gene mobilization. Antimicrob. Agents Chemother. 55, 5370–5373. doi: 10.1128/AAC.05202-5211
De Maio, N., Shaw, L. P., Hubbard, A., George, S., Sanderson, N. D., Swann, J., et al. (2019). Comparison of long-read sequencing technologies in the hybrid assembly of complex bacterial genomes. Microb. Genom. 5:e000294. doi: 10.1099/mgen.0.000294
Decano, A. G., Ludden, C., Feltwell, T., Judge, K., Parkhill, J., and Downing, T. (2019). 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, e00130–e00139. doi: 10.1128/mSphere.00130-119
Evans, D. R., Griffith, M. P., Sundermann, A. J., Shutt, K. A., Saul, M. I., Mustapha, M. M., et al. (2020). Systematic detection of horizontal gene transfer across genera among multidrug-resistant bacteria in a single hospital. Elife 9:e53886. doi: 10.7554/eLife.53886
Fattouh, R., Tijet, N., McGeer, A., Poutanen, S. M., Melano, R. G., and Patel, S. N. (2016). What is the appropriate meropenem MIC for screening of carbapenemase-producing Enterobacteriaceae in low-prevalence settings? Antimicrob. Agents Chemother. 60, 1556–1559. doi: 10.1128/AAC.02304-2315
Feldgarden, M., Brover, V., Haft, D. H., Prasad, A. B., Slotta, D. J., Tolstoy, I., et al. (2019). Using the NCBI AMRFinder tool to determine antimicrobial resistance genotype-phenotype correlations within a collection of NARMS isolates. bioRxiv [Preprint] doi: 10.1101/550707
Garrison, E., and Marth, G. (2012). Haplotype-based variant detection from short-read sequencing. arXiv [Preprint] Available online at: https://arxiv.org/abs/1207.3907 (accessed June 1, 2020).
Hardiman, C. A., Weingarten, R. A., Conlan, S., Khil, P., Dekker, J. P., Mathers, A. J., et al. (2016). Horizontal transfer of carbapenemase-encoding plasmids and comparison with hospital epidemiology data. Antimicrob. Agents Chemother. 60, 4910–4919. doi: 10.1128/AAC.00014-16
Hargreaves, M. L., Shaw, K. M., Dobbins, G., Snippes Vagnone, P. M., Harper, J. E., Boxrud, D., et al. (2015). Clonal dissemination of Enterobacter cloacae harboring blaKPC-3 in the upper midwestern United States. Antimicrob. Agents Chemother. 59, 7723–7734. doi: 10.1128/AAC.01291-1215
Hendrickx, A. P. A., Landman, F., de Haan, A., Borst, D., Witteveen, S., van Santen, M., et al. (2020). Plasmid diversity among genetically related Klebsiella pneumoniae bla KPC-2 and blaKPC-3 isolates collected in the dutch national surveillance. bioRxiv [Preprint] doi: 10.1101/2020.01.23.917781
Jia, B., Raphenya, A. R., Alcock, B., Waglechner, N., Guo, P., Tsang, K. K., et al. (2016). CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 45, D566–D573. doi: 10.1093/nar/gkw1004
Johnson, D. L., Mead, K. R., Lynch, R. A., and Hirst, D. V. L. (2013). Lifting the lid on toilet plume aerosol: a literature review with suggestions for future research. Am. J. Infect. Control 41, 254–258. doi: 10.1016/j.ajic.2012.04.330
Kitchel, B., Rasheed, J. K., Patel, J. B., Srinivasan, A., Navon-Venezia, S., Carmeli, Y., et al. (2009). Molecular epidemiology of KPC-producing Klebsiella pneumoniae isolates in the United States: clonal expansion of multilocus sequence type 258. Antimicrob. Agents Chemother. 53, 3365–3370. doi: 10.1128/AAC.00126-129
Lemon, J. K., Khil, P. P., Frank, K. M., and Dekker, J. P. (2017). Rapid nanopore sequencing of plasmids and resistance gene detection in clinical isolates. J. Clin. Microbiol. 55, 3530–3543. doi: 10.1128/JCM.01069-1017
Li, B., Feng, J., Zhan, Z., Yin, Z., Jiang, Q., Wei, P., et al. (2018). Dissemination of KPC-2-encoding IncX6 plasmids among multiple Enterobacteriaceae species in a single chinese hospital. Front. Microbiol. 9:478. doi: 10.3389/fmicb.2018.00478
Li, H. (2013). Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv [Preprint]. Available online at: https://arxiv.org/abs/1303.3997 (accessed June 1, 2020).
Martin, J., Phan, H. T. T., Findlay, J., Stoesser, N., Pankhurst, L., Navickaite, I., et al. (2017). 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. 72, 3025–3034. doi: 10.1093/jac/dkx264
Mathers, A. J., Crook, D., Vaughan, A., Barry, K. E., Vegesana, K., Stoesser, N., et al. (2019). Klebsiella quasipneumoniae provides a window into carbapenemase gene transfer, plasmid rearrangements, and patient interactions with the hospital environment. Antimicrob. Agents Chemother. 63, e02513–e02518. doi: 10.1128/AAC.02513-2518
Mathers, A. J., Stoesser, N., Chai, W., Carroll, J., Barry, K., Cherunvanky, A., et al. (2017). Chromosomal integration of the Klebsiella pneumoniae carbapenemase gene, blaKPC, in Klebsiella. Antimicrob. Agents Chemother. 61, e1823–e1826. doi: 10.1128/AAC.01823-1816
Mathers, A. J., Stoesser, N., Sheppard, A. E., Pankhurst, L., Giess, A., Yeh, A. J., et al. (2015). Klebsiella pneumoniae carbapenemase (KPC)-producing K. pneumoniae at a single institution: insights into endemicity from whole-genome sequencing. Antimicrob. Agents Chemother. 59, 1656–1663. doi: 10.1128/AAC.04292-4214
Ondov, B. D., Treangen, T. J., Melsted, P., Mallonee, A. B., Bergman, N. H., Koren, S., et al. (2016). Mash: fast genome and metagenome distance estimation using MinHash. Genom. Biol. 17:132. doi: 10.1186/s13059-016-0997-x
Sheppard, A. E., Stoesser, N., Wilson, D. J., Sebra, R., Kasarskis, A., Anson, L. W., et al. (2016). Nested russian doll-like genetic mobility drives rapid dissemination of the carbapenem resistance gene blaKPC. Antimicrob. Agents Chemother. 60, 3767–3778. doi: 10.1128/AAC.00464-416
Simão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V., and Zdobnov, E. M. (2015). BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210–3212. doi: 10.1093/bioinformatics/btv351
Stoesser, N., Sheppard, A. E., Peirano, G., Anson, L. W., Pankhurst, L., Sebra, R., et al. (2017). Genomic epidemiology of global Klebsiella pneumoniae carbapenemase (KPC)-producing Escherichia coli. Sci. Rep. 7:5917. doi: 10.1038/s41598-017-06256-6252
Tamma, P. D., Huang, Y., Opene, B. N. A., and Simner, P. J. (2016). Determining the optimal carbapenem mic that distinguishes carbapenemase-producing and non-carbapenemase-producing carbapenem-resistant Enterobacteriaceae. Antimicrob. Agents Chemother. 60, 6425–6429. doi: 10.1128/AAC.00838-816
Tham, C. Y., Tirado-Magallanes, R., Goh, Y., Fullwood, M. J., Koh, B. T. H., Wang, W., et al. (2019). NanoVar: accurate characterization of patients’ genomic structural variants using low-depth nanopore sequencing. Genome Biol. 21:56. doi: 10.1101/662940
Van Dorp, L., Wang, Q., Shaw, L. P., Acman, M., Brynildsrud, O. B., Eldholm, V., et al. (2019). Rapid phenotypic evolution in multidrug-resistant Klebsiella pneumoniae hospital outbreak strains. Microb. Genom. 5:e000263. doi: 10.1099/mgen.0.000263
Walker, B. J., Abeel, T., Shea, T., Priest, M., Abouelliel, A., Sakthikumar, S., et al. (2014). Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One 9:e112963. doi: 10.1371/journal.pone.0112963
Waterhouse, R. M., Seppey, M., Simão, F. A., Manni, M., Ioannidis, P., Klioutchnikov, G., et al. (2017). BUSCO applications from quality assessments to gene prediction and phylogenomics. Mol. Biol. Evol. 35, 543–548. doi: 10.1093/molbev/msx319
Wick, R. R., Judd, L. M., Gorrie, C. L., and Holt, K. E. (2017). Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput. Biol. 13:e1005595. doi: 10.1371/journal.pcbi.1005595
Woodworth, K. R., Walters, M. S., Weiner, L. M., Edwards, J., Brown, A. C., Huang, J. Y., et al. (2018). Vital signs: containment of novel multidrug-resistant organisms and resistance mechanisms - United States, 2006-2017. MMWR. Morb. Mortal. Wkly. Rep. 67, 396–401. doi: 10.15585/mmwr.mm6713e1
World Health Organization [WHO] (2017). Guidelines for the Prevention and Control of Carbapenem-Resistant Enterobacteriaceae, Acinetobacter baumannii and Pseudomonas aeruginosa in Health Care Facilities. Geneva: World Health Organization
Wyres, K. L., Wick, R. R., Judd, L. M., Froumine, R., Tokolyi, A., Gorrie, C. L., et al. (2019). Distinct evolutionary dynamics of horizontal gene transfer in drug resistant and virulent clones of Klebsiella pneumoniae. PLoS Genet. 15:e1008114. doi: 10.1371/journal.pgen.1008114
Keywords: carbapenem-resistant enterobacteriaceae, klebsiella pneumoniae carbapenemase, horizontal gene transfer, plasmids, long-read sequencing, hybrid genome assembly, molecular epidemiology
Citation: Prussing C, Snavely EA, Singh N, Lapierre P, Lasek-Nesselquist E, Mitchell K, Haas W, Owsiak R, Nazarian E and Musser KA (2020) Nanopore MinION Sequencing Reveals Possible Transfer of blaKPC–2 Plasmid Across Bacterial Species in Two Healthcare Facilities. Front. Microbiol. 11:2007. doi: 10.3389/fmicb.2020.02007
Received: 03 June 2020; Accepted: 29 July 2020;
Published: 19 August 2020.
Edited by:John W. A. Rossen, University Medical Center Groningen, Netherlands
Reviewed by:Johann Pitout, University of Calgary, Canada
Roberto Gustavo Melano, Public Health Ontario, Canada
Copyright © 2020 Prussing, Snavely, Singh, Lapierre, Lasek-Nesselquist, Mitchell, Haas, Owsiak, Nazarian and Musser. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Catharine Prussing, email@example.com
†Present address: Emily A. Snavely, ARUP Laboratories, University of Utah Health, Salt Lake City, UT, United States