Original Research ARTICLE
The Genetic Structures of an Extensively Drug Resistant (XDR) Klebsiella pneumoniae and Its Plasmids
- 1State Key Laboratory of Microbial Technology, Microbial Technology Institute, Shandong University, Qingdao, China
- 2Department of Gastroenterology, Qilu Hospital of Shandong University, Jinan, China
- 3Department of Stomatology, Qilu Hospital of Shandong University, Qingdao, China
- 4Department of Medical Biology, Faculty of Health Sciences, University of Tromsø, Tromsø, Norway
Multi-, extensively-, and pan-drug resistant bacteria are a threat to our health today, because their wide resistance spectra make their infections difficult to cure. In this work, we isolated an extensively drug resistant (XDR) Klebsiella pneumoniae 2-1 strain from the stool sample of a patient diagnosed of colorectal cancer. K. pneumoniae 2-1 was found to be resistant to all the antibiotics tested except for cefepime, tigecycline, and ceftazidime-avibactam. By sequencing the complete genome of K. pneumoniae 2-1, we found it contains a chromosome of 5.23 Mb and two circular plasmids with the size of 246 and 90 kb. The larger plasmid, pKP21HI1 was found to be a new conjugation-defective plasmid belonging to incompatibility group HI1B and a new sequence type. Further comparative genomics analysis and antimicrobial resistance gene analysis showed that although a great deal of changes took place on the chromosome of K. pneumoniae 2-1 in comparison with the reference genome, the extensively drug resistance phenotype of K. pneumoniae 2-1 is primarily due to the two multidrug resistant plasmids it contains. This work explains the genetic and mechanistic basis of the extensive drug resistance of K. pneumoniae 2-1, and found that plasmids play key roles in the strong antibiotic resistance of bacteria.
Klebsiella pneumoniae is the most significant clinical species of the Klebsiella genus, and also one of the most frequently observed Gram negative opportunistic pathogens in humans (Podschun and Ullmann, 1998; Navon-Venezia et al., 2017). K. pneumoniae can lead to a variety of diseases including urinary tract infections, pneumonia, bacteremia, and liver abscess (Podschun and Ullmann, 1998; Navon-Venezia et al., 2017).
A large catalog of antibiotics such as β-lactams and aminoglycosides are effective in controlling and curing infections caused by K. pneumoniae (Bush and Jacoby, 2010; Krause et al., 2016). However, antimicrobial resistance (AMR) caused by the overuse and misuse of antibiotics severely reduces the effectiveness of these antibiotics, leading to increasing difficulties in the treatment of K. pneumoniae. For instance, according to the report of a state-wide AMR surveillance program in China, CHINET, the resistance rate of imipenem and meropenem in K. pneumoniae increased significantly (from 3.0 to 10.5% for imipenem and from 2.9 to 13.4% for meropenem) between 2005 and 2014 (Hu et al., 2016). Therefore, antibiotic resistance has become a focus in combating K. pneumoniae infections in the last decade.
To make the AMR scenario worse, in recent years, infections caused by multidrug-resistant (MDR, non-susceptibility to three or more antimicrobial agent categories), extensively drug-resistant (XDR, susceptible to at most two antimicrobial agent categories), and even pandrug-resistant (PDR, non-susceptibility to all drugs) K. pneumoniae have been frequently reported (Falagas et al., 2005, 2008; Leavitt et al., 2010; Karaiskos and Giamarellou, 2014; Huang et al., 2018; Krapp et al., 2018). These K. pneumoniae strains are resistant to most or even all the antibiotics used (Magiorakos et al., 2012). The infection of these “superbugs” often leads to immense difficulties in treatment, and sometimes even the death of infected patients (Giske et al., 2008; Hersh et al., 2012). Oftentimes, these strains emerge in ICU wards or on patients treated with organ transplant procedures due to the long-term and extensive use of various antibiotics (Maseda et al., 2014).
Several mechanisms are exploited by bacteria for resistance against antibiotics: mutations on targets of antibiotics may reduce the effectiveness of antibiotics, the alteration and inactivation of antibiotics may be catalyzed by proteins encoded by chromosomal or plasmid-borne genes, mutated porins and efflux pumps may lead to lower levels of antibiotics in the cytoplasm and subsequently reduces their effective dosages (Zhao et al., 2009; Li et al., 2015; Sharma et al., 2016). Antimicrobial resistance genes (ARGs) that take advantage of these resistance mechanisms may be disseminated between bacteria via horizontal gene transfer mechanisms and spread from environmental bacteria to pathogens (Sharma et al., 2016). Therefore, they pose a greater threat to overall human health. To date, a large consortium of ARGs has been discovered that impacts essentially every class of antibiotics. For instance, 51 ARGs have been registered in the Comprehensive Antibiotic Resistance Database (CARD) for carbapenems, one of the last-line antibiotics against Gram-negative bacteria (Nordmann et al., 2011a,b). This list of ARGs include frequently observed blaOXA−48, blaKPC−2, blaNDM−1, blaIMP, and blaVIM in K. pneumoniae, and is constantly growing (Kliebe et al., 1985; Bradford et al., 2004; Yu et al., 2012).
The multi-, extensively-, and even pan-drug resistance phenotypes of bacteria is usually the result of horizontal gene transfer that gathers ARGs into one cell (Juhas, 2013). Many mobile genetic elements, including plasmids, integrons, and transposons, are capable of expressing genetically linked and co-expressed ARG arrays that, upon dissemination, bring several ARGs all at once, leading to multidrug resistance (Salabi et al., 2013). In recent years, XDR K. pneumoniae turned into an emerging and dangerous pathogen, particularly with the emergence of carbapenem-resistant XDR K. pneumoniae (Santino, 2013; Karaiskos and Giamarellou, 2014; Pontikis et al., 2014; Lim et al., 2015). Several therapeutic options are available for XDR K. pneumoniae. Fosfomycin has been attempted to treat urinary tract infections and gastrointestinal infections caused by XDR and PDR Enterobacteriaceae (Falagas et al., 2010; Leavitt et al., 2010; Braun et al., 2018). Colistin and tigecycline are the last resort treatments for serious carbapenem-resistant K. pneumoniae infections (Olaitan et al., 2014; Piedra-Carrasco et al., 2018). However, AMR against these antibiotics arose, such as fosA3 for fosfomycin, mcr-1 for polymyxin, and tetA mutation for tigecycline (Jiang et al., 2015; Liu et al., 2015; Du et al., 2018). Further successful attempts for the treatment of XDR and PDR K. pneumoniae were made, including antibiotic combination therapies such as using double carbapenems (Piedra-Carrasco et al., 2018), tigecycline plus meropenem (Lim et al., 2015), and β-lactam/β-lactamase inhibitor combos such as ceftazidime plus avibatam (Schimmenti et al., 2018).
Plasmids are the most important carriers for ARGs in MDR K. pneumoniae, and are found in almost all antimicrobial resistant K. pneumoniae isolates (Dolejska et al., 2013; Freire Martín et al., 2014; Pitout et al., 2015; Navon-Venezia et al., 2017). A total of 306 complete K. pneumoniae plasmid sequences are currently available at Genbank (Navon-Venezia et al., 2017). In particular, pKpQIL and pKPN3 that belong to incompatibility group FIIK are commonly found in K. pneumoniae belonging to the epidemic clonal group CG258 (Chen et al., 2014; Navon-Venezia et al., 2017). These two plasmids contain important ARGs such as blaKPC−2 and blaKPC−3, and are frequently reported globally (Leavitt et al., 2010; García-Fernández et al., 2012; Chen et al., 2013, 2014).
In this study, K. pneumoniae 2-1, an XDR K. pneumoniae, was isolated from the stool sample of a patient with colorectal cancer. Two large plasmids were found in this strain. The full genomic sequence of this strain, including both chromosomal DNA sequence and the DNA sequence of two plasmids it harbors, were obtained. Comparative genomic sequence analysis showed one of the plasmids, pKP21HI1, is a new plasmid belonging to IncHI1B and forms a new sequence type. Further analysis of both chromosomal and plasmid-borne ARGs showed the XDR phenotype in this strain is primarily contributed by the ARGs harbored by the two multidrug resistance plasmids. Findings of this work suggest that attention should be focused on multidrug resistance plasmids, as they could determine the XDR phenotype of K. pneumoniae, are highly mobile, and are difficult to contain.
Materials and Methods
K. pneumoniae 2-1 used in this work was isolated from the stool sample of a colorectal cancer patient in Qilu Hospital, Jinan. The patient was treated for lasting abdominal pain and constipation, and was initially prescribed with levofloxacin because he was misdiagnosed of enteritis. Strain identification was performed by analyzing the 16S rDNA sequence of this strain. K. pneumoniae HS11286 strain is a previously identified model clinical strain that was generously gifted from Prof. Hongyu Ou from Shanghai Jiao Tong University (Bi et al., 2015).
Antibiotic Susceptibility Testing
All antibiotic susceptibility assays were performed following CLSI standard (Wayne, 2018b). Disk diffusion assays were performed as previously documented (Vading et al., 2011). Minimum Inhibition Concentrations (MICs) was determined with the agar dilution method as previously reported (Wiegand et al., 2008). Escherichia coli ATCC 25922 was used as a reference strain as recommended by CLSI standard (Wayne, 2018a).
Transfer of Plasmids to E. coli DH5α
To transfer plasmids of K. pneumoniae 2-1 to E. coli DH5α, plasmids were extracted from K. pneumoniae 2-1, and transferred to the recipient E. coli DH5α strain by either chemical transformation or electrotransformation. Successful transfer was determined by trimethoprim resistance (for pKP21HI1) or tetracycline resistance (for pKP21AC2). Further validation was performed by extracting plasmids from plasmid-containing E. coli DH5α and subsequent PCR amplification of aadA2, sul1, dfrA12 for pKP21HI1 or sul2, floR for pKP21AC2.
Conjugations assays were performed in order to test whether plasmids in K. pneumoniae 2-1 can be transmitted via conjugal transfer, following previously published procedures (Borgia et al., 2012). E. coli J53 was used as the recipient strain. Successful conjugal transfer of plasmid was determined by the dual resistance of sodium azide and trimethoprim/streptomycin/tetracycline/chloramphenicol.
Whole Genome Sequencing
The total DNA of K. pneumoniae 2-1 was extracted using the SDS method as previously reported (Natarajan et al., 2016). Total DNA samples were shredded into 10 kb fragments using g-TUBE (Covaris Inc., Woburn, MA, US) to construct a 10 kb SMRT Bell library for PacBio sequencing. Total DNA samples were shredded with ultrasonication to construct a 350 bp library for Illumina sequencing. Libraries were, respectively, sequenced with a PacBio Sequel system (Pacific Biosciences of California, Inc., Menlo Park, CA, US) and an Illumina HiSeq X 10 system (Illumina Inc., San Diego, CA, US) at PE150 mode. DNA sequences were assembled with SMRT Link v5.1.0 software (Ardui et al., 2018). Gaps were manually closed by PCR amplifying gap-containing DNA using primers targeting each end of the gap and subsequent sequencing. Sequence data were deposited in Genbank, with accession numbers of CP031562, CP031563, and CP031564.
Gene model prediction was performed using GeneMarkS Version 4.1.7 (Besemer et al., 2001). The predicted gene models were subsequently annotated with Gene Ontology (GO) (Blake et al., 2015), Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa and Goto, 2000), Cluster of Orthologous Groups of proteins (COG) (Natale et al., 2000), Non-Redundant Protein Database (NR) (Li et al., 2002), Pfam (Punta et al., 2011), and Swiss-Prot (Boeckmann et al., 2003) databases.
For ARG analysis, the Resistance Gene Identifier (RGI) v4.1.0 tool of The CARD was used (Jia et al., 2017).
For classification of K. pneumoniae 2-1, the multilocus sequencing typing method was used as previously reported (Diancourt et al., 2005).
For classification of plasmids, the PlasmidFinder and pMLST webtools were used (Carattoli et al., 2014).
Phylogenetic analysis was performed using the maximum likelihood (ML) method with the MEGA 7.0.21 software (Hall, 2013).
Comparative genomics analysis was performed using the MUMmer v3.2.3 software (Kurtz et al., 2004). InDels were identified using the LASTZ v1.03.54 software.
This study was carried out in accordance with the recommendations of Scientific Ethics Committee of Qilu Hospital of Shandong University. The protocol was approved by the Scientific Ethics Committee of Qilu Hospital of Shandong University. All subjects gave written informed consent in accordance with the Declaration of Helsinki.
Isolation of an XDR K. pneumoniae
A K. pneumoniae 2-1 strain was isolated from the stool sample of a 36-year-old male patient with colorectal cancer. Antimicrobial susceptibility testing of 23 antibiotics was performed using both K-B disk diffusion assay and agar dilution method, showing it is an XDR strain that is resistant to nearly all antibiotics tested (Table 1). A more conservative approach was undertaken to determine the susceptibility of K. pneumoniae 2-1: if the strain is sensitive to an antibiotic using either K-B disk diffusion assay or agar dilution method, the strain is considered sensitive to this antibiotic. The antibiotics that K. pneumoniae 2-1 is non-susceptible to include β-lactams (ampicillin, amoxicillin, ceftazidime, cefotaxime, cefoxitin, cefoperazone), carbapenem (imipenem), aminoglycosides (spectinomycin, kanamycin, streptomycin), quinolones (gatifloxacin, ciprofloxacin, nalidixic acid), diaminopyrimidine (trimethoprim), sulfonamide (sulfisoxazole), rifampicin (rifampicin), macrolide (erythromycin), phenicol (chloramphenicol), tetracycline (tetracycline), and polymyxins (polymyxin B and polymyxin E). K. pneumoniae 2-1 was found to be susceptible to only two antibiotics tested: a fourth generation cephalosporin cefepime and a last-line antibiotic tigecycline. Out of the three β-lactam/β-lactamase inhibitor combos, K. pneumoniae 2-1 was found to be sensitive to only one combo: ceftazidime-avibactam. The antibiotic susceptibility of the reference clinical K. pneumoniae HS11286 strain was also performed as a reference (Table 1).
The Genetic Features of K. pneumoniae 2-1
Total DNA was extracted from K. pneumoniae 2-1 and sequenced on a PacBio system. A circular genomic DNA with the size of 5,230,042 basepairs was assembled. In addition, two circular plasmids, respectively, with the size of 245,534 and 89,508 basepairs were found. The plasmids were subsequently denominated pKP21HI1 and pKP21AC2. The GC content for the genomic DNA, pKP21HI1, and pKP21AC2 are, respectively, 57.57, 51.43 and, 52.92%. The numbers of genes encoded by the genomic DNA, pKP21HI1 and pKP21AC2 are respectively 4972, 268, and 114. Functional annotation was performed with a variety of databases including GO, KEGG, COG, NR, Pfam and Swiss-Prot (Figure 1, Table S1).
Figure 1. Genome map of K. pneumoniae 2-1 chromosome and plasmids. For chromosome, from outside to inside, each ring represents encoding genes, COG, KEGG, GO, ncRNA, GC content, and genomic GC skew value. For plasmids, from outside to inside, each ring represents COG, GC content, and genomic GC skew value. Different color represents different annotation.
Classification of K. pneumoniae 2-1 and its Plasmids
MLST classification of K. pneumoniae 2-1 was performed using the PubMLST database based on genetic variation in seven housekeeping genes (rpoB, gapA, mdh, pgi, phoE, infB, and tonB) (Diancourt et al., 2005), confirming a 2-9-2-1-13-1-4 configuration and supporting the subsequent classification of K. pneumoniae 2-1 as ST726.
The classification of pKP21HI1 and pKP21AC2 was performed based on the typing of its replicon sequences, using the PlasmidFinder database (https://cge.cbs.dtu.dk/services/PlasmidFinder/). These two plasmids were, respectively, classified as IncHI1B and IncA/C2 plasmids, suggested by the presence of repH1B and repA. This classification is further confirmed by phylogenetic analysis of replicon sequences of pKP21HI1 and pKP21AC2, together with replicon sequences of known plasmids of different incompatible groups (Figure 2). This phylogenetic analysis suggested that pKP21HI1 forms a clade with other IncHI1B plasmids, while pKP21AC2 forms a clade with other IncA/C2 plasmids. The sequence types of these two plasmids were further determined with the plasmid multilocus sequence typing (pMLST) analysis with the pMLST webtool (Carattoli et al., 2014). The pKP21AC2 plasmid was confirmed to be a ST-3 plasmid belonging to the IncA/C incompatibility group based on a repA-parA-parB-A053 configuration of 2-2-2-1. However, only one (hcm1043) of the six (hcm1043, hcm1099, hcm1064, hcm1116, hcm1178ac, hcm1259) genes used for sequence typing of IncHI1 plasmids was found in pKP21HI1, suggesting it belongs to a new sequence type of IncHI1 plasmids. Indeed, pMLST analysis was unable to assign pKP21HI1 to any sequence type, and BLAST search of pKP21HI1 in Genbank showed the most similar plasmid (Klebsiella aerogenes strain AR_0161 plasmid unnamed, Genbank accession number CP028952.1) is 206 kbp larger than pKP21HI1.
Figure 2. Phylogenetic analysis of pKP21AC2 and pKP21HI1. Phylogenetic tree of replicon sequence of pKP21HI1, pKP21AC2, and plasmids belonging to different incompatible group was shown. The phylogenetic tree was calculated from 500 individual trees. Bootstrap values were shown on each node. Bar: evolutionary distance of 0.5. Circled are the two plasmids found in K. pneumoniae 2-1 and the most similar plasmid to pKP21HI1.
Comparative Genomic Analysis of K. pneumoniae 2-1 and the Reference K. pneumoniae Strain HS11286
A comparison between K. pneumoniae 2-1 and the reference clinical K. pneumoniae HS11286 chromosomal DNA found significant differences between the two genomes. A total of 22,360 SNPs were found between the two genomes, including 19,659 SNPs that fall into the CDS region, of which 3,789 are non-synonymous. In particular, 25 SNPs lead to premature termination of the gene products (Table S2). A total of 116 InDels were found in K. pneumoniae 2-1 in comparison with K. pneumoniae HS11286, leading to the frameshift of 102 gene products and the premature termination of 1 CDS (Table S2). One hundred and sixteen structural variations were found in K. pneumoniae 2-1, including 34 complex InDels, 34 deletions, 36 insertions, and 1 inversion (Table S2). These changes do not appear to affect the antimicrobial resistance of K. pneumoniae 2-1.
pKP21HI1 as a new Conjugation-Defective Multidrug Resistant Plasmid
The analysis of the presence of ARGs in both pKP21HI1 and its most similar plasmid CP028952.1 showed the presence of, respectively, 18 and 13 ARGs (Table 2), suggesting both plasmids are multidrug resistant plasmids. Comparison between these two plasmids showed significant differences. In comparison with CP028952.1, a large fragment (206 kbp) is missing in pKP21HI1. pKP21HI1 codes for 50 genes that are absent in CP028952.1, while CP028952.1 codes for 278 genes that are absent in pKP21HI1 (Table S3). Of particular interest, the genes absent in pKP21HI1 include 29 transposase/intergrase-coding genes and 13 conjugation/ plasmid transfer protein-coding genes (traN, traX, traG, traF, traI, traD, traE, traB, traV). This is a strong suggestion that CP028952.1 is a more mobile plasmid than pKP21HI1. Indeed, no conjugal protein-coding genes could be found in pKP21HI1, and conjugation assays in attempt to transfer pKP21HI1 to E. coli J53 were unsuccessful. These results suggest that pKP21HI1 is a conjugation-defective multidrug resistant plasmid that limits its mobility between bacteria.
ARGs Harbored by K. pneumoniae 2-1 and Its Plasmids
K. pneumoniae 2-1 and its plasmids are annotated for ARGs with the CARD database (Table 3). A total of 33 ARGs were found, of which 23 are plasmid-borne while the rest are chromosomal. By comparing the ARG catalog of K. pneumoniae 2-1 and K. pneumoniae HS11286, we found little difference except that K. pneumoniae 2-1 encodes two additional quinolone resistance genes oqxA and oqxB. The mutations of several genes that lead to antibiotic resistance were found in the chromosome of K. pneumoniae 2-1, including mutation on elongation factor Tu that leads to resistance to pulvomycin, mutation in uhpT that confers resistance to fosfomycin, parC and gyrA mutation that leads to resistance to quinolones, marR mutation that leads to multidrug resistance, and mutation in ftsI that confers resistance to β-lactams. Therefore, K. pneumoniae 2-1 chromosome encodes genes responsible for quinolone and β-lactam resistance phenotypes found in this work (Table 1). pKP21AC2 hosts 5 ARGs while pKP21HI1 hosts 18 ARGs, including genes for aminoglycoside resistance, quinolone resistance, β-lactam resistance, sulfonamide resistance, chloramphenicol resistance, tetracycline resistance, macrolide resistance, rifamycin resistance, and trimethoprim resistance. The genes encoded by these two plasmids were found responsible for aminoglycoside, sulfonamide, chloramphenicol, tetracycline, macrolide, rifamycin, and trimethoprim resistance that were found in phenotypic analysis but were not encoded by the chromosome. Two class I integrons were found in pKP21HI1, each containing an ARG-rich gene cassette array in the organization of dfrA12-aadA2 and aac-(6')-Ib-cr-blaOXA−1-catB3-arr-3 (Figure S1). These findings suggest that the XDR phenotype is a result of the presence of the two multidrug resistance plasmids harbored by K. pneumoniae 2-1.
To further confirm the roles of the plasmids in antimicrobial resistance, we extracted the plasmids and transferred them to E. coli DH5α. The antibiotic susceptibility tests were performed, showing that E. coli DH5α containing pKP21HI1 is resistant to β-lactams, aminoglycosides, rifamycin, macrolides, quinolones, trimethoprim, and sulfonamides, while E. coli DH5α containing pKP21AC2 is resistant to aminoglycosides, quinolones, rifamycin, macrolides, sulfonamides and tetracycline.
XDR K. pneumoniae is an eminent threat to human health in an era when the discovery of new antibiotics lags behind the emergence and dissemination of antimicrobial resistance. Therefore, understanding the genetic and mechanistic basis for XDR becomes crucial, which may provide a hint on finding solutions to prevent its spread. In this work, with high throughput sequencing, we obtained the complete sequence of an XDR K. pneumoniae 2-1 strain isolated from a patient with colorectal cancer, and identified the genetic basis for its XDR phenotype.
K. pneumoniae 2-1 was found resistant to nearly all antibiotics tested except for cefepime, tigecycline, and ceftazidime-avibactam. A total of 33 ARGs and 6 gene mutations responsible for AMR were found on the chromosome and plasmids of K. pneumoniae 2-1. ARGs for all the resistant antibiotics were found except for carbapenems and polymyxins. The presence of blaOXA−1 would have been responsible for carbapenem resistance, as blaOXA−1-containing strains that are resistant to carbapenems have been previously reported (Sugumar et al., 2014). However, the blaOXA−1-harboring pKP21HI1 didn't increase the resistance of E. coli DH5α to imipenem (Table 4), making this hypothesis unlikely true. Considering neither plasmid was able to increase the resistance of E. coli DH5α (Table 4), it is likely that a yet unknown chromosomal mechanism is responsible for the carbapenem resistance of K. pneumoniae 2-1. K. pneumoniae 2-1 also encodes a series of pmr phosphoethanolamine transferases (basR, arnA, eptA, ugd, pmrF) that involve in polymyxin resistance, when we set the search parameter for CARD database as “Loose”. Therefore, all the resistance phenotypes found in this work are accounted for except for carbapenems, all of which involve one or several ARGs.
Table 4. The antibiotic susceptibility of E. coli DH5α (DH5α), E. coli DH5α containing pKP21HI1 (DH5α-pHI1) and E. coli DH5α containing pKP21AC2 (DH5α-pAC2).
One particularly interesting finding in this work is that the majority of ARGs are encoded by the two plasmids of K. pneumoniae 2-1. The resistance to aminoglycosides, sulfonamides, chloramphenicol, tetracycline, macrolides, rifamycin, and trimethoprim are not encoded by the chromosome, and completely depend on the presence of the two plasmids. Therefore, the two multidrug plasmids play a key role in the XDR phenotype of K. pneumoniae 2-1, which is worrisome as plasmids can easily transfer between bacteria, disseminating ARGs.
The analysis of E. coli DH5α strains, respectively, containing pKP21HI1 and pKP21AC2 further confirmed the role of the two plasmids in antimicrobial resistance. pKP21HI1 but not pKP21AC2 contains ARGs for β-lactams and diaminopyrimidines, while the pKP21HI1-containing E. coli DH5α but not the pKP21AC2-containing E. coli DH5α is resistant to β-lactams and trimethoprim. pKP21AC2 but not pKP21HI1 contains ARGs for tetracycline, while the pKP21AC2-containing E. coli DH5α but not the pKP21HI1-containing E. coli DH5α is resistant to tetracycline. Both plasmids contain ARGs for aminoglycosides and sulfonamides, while both the pKP21HI1-containing E. coli DH5α and the pKP21AC2-containing E. coli DH5α are resistant to these classes of antibiotics. Neither plasmids contain ARGs for carbapenems, while neither plasmid-containing E. coli DH5α strain is resistant to carbapenems, further agreeing with the hypothesis that a chromosomal mechanism could be responsible for carbapenem resistance. Although both plasmids carry ARGs for chloramphenicol, neither plasmid-containing E. coli DH5α strain was shown to be resistant to chloramphenicol. Therefore, the resistance to chloramphenicol requires the presence of both plasmids in the cell. The disagreement between the presence of ARGs in the plasmids and the antimicrobial resistance phenotype is that although pKP21AC2 lacks ARGs for quinolones, macrolides and rifamycin, pKP21AC2-harboring E. coli DH5α is resistant to nalidixic acid, erythromycin, and rifamycin. The resistance of pKP21AC2-harboring E. coli DH5α to nalidixic acid is due to the innate resistance of E. coli DH5α to this antibiotic (Table 4). By further investigations of pKP21AC2 sequences, we found that pKP21AC2 contains a gene that may encode a SoxR mutant conferring resistance to rifamycin (sequence identity 39.71%) and a gene that may encode a mutated repressor (NalD) that could lead to the overexpression of the MexAB-OprM efflux pump for resistance to macrolides (sequence identity 36.17%). Therefore, we hypothesize that these distant homologs to existing antimicrobial resistance determinants may be responsible for rifamycin and macrolide resistance in pKP21AC2. It needs to be addressed that in general, plasmid-harboring E. coli DH5α showed lower levels of antimicrobial resistance in comparison with K. pneumoniae 2-1. This phenomenon suggests that the combination of MDR plasmids and chromosomal mutations can significant increase the overall antibiotic resistance level.
The findings in this work showed the importance of plasmids in XDR. Out of the two multidrug resistance plasmids identified in K. pneumoniae 2-1, pKP21HI1 is a new plasmid that differs significantly to all previously known plasmids. Indeed, the most similar plasmid CP028952.1 from Klebsiella aerogenes is 206 kb larger than pKP21HI1, and analysis of its genetic content revealed that it belongs to a new sequence type of incompatibility group HI1B. This plasmid encodes 18 ARGs that leads to resistance of 9 major classes of antibiotics, as well as two intact Class 1 integrons, suggesting its significant role in AMR and its dissemination. Fortunately, this plasmid appears to be conjugation-defective, therefore limiting its mobility. Nevertheless, it can still disseminate between bacteria via other processes such as natural transformation.
The identification of the genetic and mechanistic basis of XDR in K. pneumoniae 2-1 stresses the role of multidrug resistance plasmids in XDR. Comparison of K. pneumoniae 2-1 and the reference genome of K. pneumoniae HS11286 showed only minor differences in antimicrobial resistance on the chromosomal level. Therefore, although the possibility of chromosomal MDR K. pneumoniae is present and reported (Huang et al., 2017; Mathers et al., 2017), proposal can be made that the external antibiotic pressure favors concentrating ARGs to existing strains in the form of multidrug plasmids, rather than leading to the evolvement of strains with chromosome-encoded XDR. With this proposal, further suggestion can be made that cutting off transmission of plasmids may be a good approach in delaying the appearance of XDR and PDR strains.
In this work, with the analysis of the genetic features of an XDR K. pneumoniae 2-1 strain isolated from the stool sample of a patient diagnosed of colorectal cancer, we found and identified a new conjugation-defective multidrug resistant pKP21HI1 plasmid that belongs to incompatibility group HI1B and a new sequence type. With further analysis of the antimicrobial resistance phenotype of K. pneumoniae 2-1 and the distribution of ARGs on its chromosome and two plasmids, conclusions can be made that the XDR phenotype should be contributed to primarily its plasmids rather than chromosome. Research from this work stressed the importance of multidrug resistant plasmids in leading to extensively- and pan-drug resistance, and lead to the proposal that efforts should be made to control plasmid dissemination in order to delay the appearance of the XDR or PDR pathogens.
LL, YM, and WW performed experiments. TY performed strain isolation. LL, ZY, XS, YS, TG, JK, MW, and HX performed bioinformatical analysis. LL, MW, and HX wrote the manuscript. MW and HX conceived of the study. HX oversaw the project. All authors read and approved the manuscript.
This work was supported by the National Key Research and Development Program of China (2017YFD0400301), the National Natural Science Foundation of China (31770042, 31770043), Shandong Province Key Research and Development Program (2016GSF121040, 2018GSF118008), the Fundamental Research Funds of Shandong University (2017JC028, 2018JC013, 2018JC027), the State Key Laboratory of Microbial Technology Open Project Funds, Shandong University (M2017-07, M2018-07), and Jinan Cultural Industry Development Fund.
Conflict of Interest Statement
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.
We wish to thank Novogene Technology Co. Ltd., Beijing, China, for high throughput sequencing. We would like to thank Prof. Hongyu Ou from Shanghai Jiao Tong University for generously providing K. pneumoniae HS11286 strain.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2018.00446/full#supplementary-material
Ardui, S., Ameur, A., Vermeesch, J. R., and Hestand, M. S. (2018). Single molecule real-time (SMRT) sequencing comes of age: applications and utilities for medical diagnostics. Nucleic Acids Res. 46, 2159–2168. doi: 10.1093/nar/gky066
Besemer, J., Lomsadze, A., and Borodovsky, M. (2001). GeneMarkS: a self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory regions. Nucleic Acids Res. 29, 2607–2618. doi: 10.1093/nar/29.12.2607
Bi, D., Jiang, X., Sheng, Z. K., Ngmenterebo, D., Tai, C., Wang, M., et al. (2015). Mapping the resistance-associated mobilome of a carbapenem-resistant Klebsiella pneumoniae strain reveals insights into factors shaping these regions and facilitates generation of a ‘resistance-disarmed' model organism. J. Antimicrob. Chemother. 70, 2770–2774. doi: 10.1093/jac/dkv204
Boeckmann, B., Bairoch, A., Apweiler, R., Blatter, M. C., Estreicher, A., Gasteiger, E., et al. (2003). The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 31, 365–370. doi: 10.1093/nar/gkg095
Borgia, S., Lastovetska, O., Richardson, D., Eshaghi, A., Xiong, J., Chung, C., et al. (2012). Outbreak of carbapenem-resistant Enterobacteriaceae containing blaNDM−1, Ontario, Canada. Clin. Infect. Dis. 55, e109–e117. doi: 10.1093/cid/cis737
Bradford, P. A., Bratu, S., Urban, C., Visalli, M., Mariano, N., Landman, D., et al. (2004). Emergence of carbapenem-resistant Klebsiella species possessing the class A carbapenem-hydrolyzing KPC-2 and inhibitor-resistant TEM-30 β-lactamases in New York City. Clin. Infect. Dis. 39, 55–60. doi: 10.1086/421495
Braun, G., Cay,ô, R., Matos, A. P., de Mello Fonseca, J., and Gales, A. C. (2018). Temporal evolution of polymyxin B-resistant Klebsiella pneumoniae clones recovered from blood cultures in a teaching hospital during a 7-year period. Int. J. Antimicrob. Agents. 51, 522–527. doi: 10.1016/j.ijantimicag.2017.08.031
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
Chen, L., Chavda, K. D., Melano, R. G., Jacobs, M. R., Koll, B., Hong, T., et al. (2014). Comparative genomic analysis of KPC-encoding pKpQIL-like plasmids and their distribution in New Jersey and New York Hospitals. Antimicrob. Agents Chemother. 58, 2871–2877. doi: 10.1128/AAC.00120-14
Chen, L., Chavda, K. D., Melano, R. G., Jacobs, M. R., Levi, M. H., Bonomo, R. A., et al. (2013). Complete sequence of a blaKPC−2 -harboring IncFIIK1 plasmid from a Klebsiella pneumoniae sequence type 258 Strain. Antimicrob. Agents Chemother. 57, 1542–1545. doi: 10.1128/AAC.02332-12
Diancourt, L., Passet, V., Verhoef, J., Grimont, P. A. D., and Brisse, S. (2005). Multilocus sequence typing of Klebsiella pneumoniae nosocomial isolates. J. Clin. Microbiol. 8, 4178–4182 doi: 10.1128/JCM.43.8.4178-4182.2005
Dolejska, M., Villa, L., Dobiasova, H., Fortini, D., Feudi, C., and Carattoli, A. (2013). Plasmid content of a clinically relevant Klebsiella pneumoniae clone from the Czech Republic producing CTX-M-15 and QnrB1. Antimicrob. Agents Chemother. 57, 1073–1076. doi: 10.1128/AAC.01886-12
Du, X., He, F., Shi, Q., Zhao, F., Xu, J., Fu, Y., et al. (2018). The rapid emergence of tigecycline resistance in blaKPC−2 harboring Klebsiella pneumoniae, as mediated in vivo by mutation in tetA during tigecycline treatment. Front. Microbiol. 9,648. doi: 10.3389/fmicb.2018.00648
Falagas, M. E., Bliziotis, I. A., Kasiakou, S. K., Samonis, G., Athanassopoulou, P., and Michalopoulos, A. (2005). Outcome of infections due to pandrug-resistant (PDR) Gram-negative bacteria. BMC Infect. Dis. 5,24. doi: 10.1186/1471-2334-5-24
Falagas, M. E., Maraki, S., Karageorgopoulos, D. E., Kastoris, A. C., Mavromanolakis, E., and Samonis, G. (2010). Antimicrobial susceptibility of multidrug-resistant (MDR) and extensively drug-resistant (XDR) Enterobacteriaceae isolates to fosfomycin. Int. J. Antimicrob. Agents 35, 240–243. doi: 10.1016/j.ijantimicag.2009.10.019
Falagas, M. E., Rafailidis, P. I., Matthaiou, D. K., Virtzili, S., Nikita, D., and Michalopoulos, A. (2008). Pandrug-resistant Klebsiella pneumoniae, Pseudomonas aeruginosa and Acinetobacter baumannii infections: characteristics and outcome in a series of 28 patients. Int. J. Antimicrob. Agents 32, 450–454. doi: 10.1016/j.ijantimicag.2008.05.016
Freire Martín, I., AbuOun, M., Reichel, R., La Ragione, R. M., and Woodward, M. J. (2014). Sequence analysis of a CTX-M-1 IncI1 plasmid found in Salmonella 4,5,12:i:–, Escherichia coli and Klebsiella pneumoniae on a UK pig farm. J. Antimicrob. Chemother. 69, 2098–2101. doi: 10.1093/jac/dku098
García-Fernández, A., Villa, L., Carta, C., Venditti, C., Giordano, A., Venditti, M., et al. (2012). Klebsiella pneumoniae ST258 producing KPC-3 identified in Italy carries novel plasmids and OmpK36/OmpK35 porin variants. Antimicrob. Agents Chemother. 56, 2143–2145. doi: 10.1128/AAC.05308-11
Giske, C. G., Monnet, D. L., Cars, O., and Carmeli, Y. (2008). Clinical and economic impact of common multidrug-resistant Gram-negative bacilli. Antimicrob. Agents Chemother. 52, 813–821. doi: 10.1128/AAC.01169-07
Hu, F. P., Guo, Y., Zhu, D. M., Wang, F., Jiang, X. F., Xu, Y. C., et al. (2016). Resistance trends among clinical isolates in China reported from CHINET surveillance of bacterial resistance, 2005–2014. Clin. Microbiol. Infect. 22, S9–S14. doi: 10.1016/j.cmi.2016.01.001
Huang, W., Wang, G., Sebra, R., Zhuge, J., Yin, C., Aguero-Rosenfeld, M. E., et al. (2017). Emergence and evolution of multidrug-resistant Klebsiella pneumoniae with both blaKPC and blaCTX−M integrated in the chromosome. Antimicrob. Agents Chemother. 61, e00076–e00017. doi: 10.1128/AAC.00076-17
Huang, Y. H., Chou, S. H., Liang, S. W., Ni, C. E., Lin, Y. T., Huang, Y. W., et al. (2018). Emergence of an XDR and carbapenemase-producing hypervirulent Klebsiella pneumoniae strain in Taiwan. J. Antimicrob. Chemother. 73, 2039–2046. doi: 10.1093/jac/dky164
Jia, B., Raphenya, A. R., Alcock, B., Waglechner, N., Guo, P., Tsang, K. K., et al. (2017). CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 45, D566–D573. doi: 10.1093/nar/gkw1004
Jiang, Y., Shen, P., Wei, Z., Liu, L., He, F., Shi, K., et al. (2015). Dissemination of a clone carrying a fosA3-harbouring plasmid mediates high fosfomycin resistance rate of KPC-producing Klebsiella pneumoniae in China. Int. J. Antimicrob. Agents 45, 66–70. doi: 10.1016/j.ijantimicag.2014.08.010
Karaiskos, I., and Giamarellou, H. (2014). Multidrug-resistant and extensively drug-resistant Gram-negative pathogens: current and emerging therapeutic approaches. Expert Opin. Pharm. 15, 1351–1370. doi: 10.1517/14656566.2014.914172
Kliebe, C., Nies, B. A., Meyer, J. F., Tolxdorff-Neutzling, R. M., and Wiedemann, B. (1985). Evolution of plasmid-coded resistance to broad-spectrum cephalosporins. Antimicrob Agents Chemother. 28, 302–307. doi: 10.1128/AAC.28.2.302
Krapp, F., Ozer, E. A., Qi, C., and Hauser, A. R. (2018). Case report of an extensively drug-resistant Klebsiella pneumoniae infection with genomic characterization of the strain and review of similar cases in the United States. Open Forum Infect. Dis. 5,ofy074. doi: 10.1093/ofid/ofy074
Kurtz, S., Phillippy, A., Delcher, A. L., Smoot, M., Shumway, M., Antonescu, C., et al. (2004). Versatile and open software for comparing large genomes. Genome Biol. 5,R12. doi: 10.1186/gb-2004-5-2-r12
Leavitt, A., Chmelnitsky, I., Ofek, I., Carmeli, Y., and Navon-Venezia, S. (2010). Plasmid pKpQIL encoding KPC-3 and TEM-1 confers carbapenem resistance in an extremely drug-resistant epidemic Klebsiella pneumoniae strain. J. Antimicrob. Chemother. 65, 243–248. doi: 10.1093/jac/dkp417
Lim, T., Cai, Y., Hong, Y., Chun Yong Chan, E., Suranthran, S., Teo, Q., et al. (2015). In vitro pharmacodynamics of various antibiotics in combination against extensively drug-resistant Klebsiella pneumoniae. Antimicrob. Agents Chemother. 59, 2515–2524. doi: 10.1128/AAC.03639-14
Liu, Y., Wang, Y., Walsh, T. R., Yi, L., Zhang, R., Spencer, J., et al. (2015). Emergence of plasmid-mediated colistin resistance mechanism MCR-1 in animals and human beings in China: a microbiological and molecular biological study. Lancet Infect. Dis. 16, 161–168. doi: 10.1016/S1473-3099(15)00424-7
Magiorakos, A. P., Srinivasan, A., Carey, R. B., Carmeli, Y., Falagas, M. E., Giske, C. G., et al. (2012). Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin. Microbiol. Infect. 18, 268–281. doi: 10.1111/j.1469-0691.2011.03570.x
Maseda, E., Mensa, J., Valía, J., Gomez-Herreras, J., Ramasco, F., Samso, E., et al. (2014). Bugs, hosts and ICU environment: countering pan-resistance in nosocomial microbiota and treating bacterial infections in the critical care setting. Rev. Esp. Anestesiol. Reanim. 61, e1–e19. doi: 10.1016/j.redar.2013.11.012
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 species is elusive but not rare. Antimicrob. Agents Chemother 61, e01823–e01816. doi: 10.1128/AAC.01823-16
Natale, D. A., Shankavaram, U. T., Galperin, M. Y., Wolf, Y. I., Aravind, L., and Koonin, E. V. (2000). Towards understanding the first genome sequence of a crenarchaeon by genome annotation using clusters of orthologous groups of proteins (COGs). Genome Biol. 1:RESEARCH0009. doi: 10.1186/gb-2000-1-5-research0009
Natarajan, V. P., Zhang, X., Morono, Y., Inagaki, F., and Wang, F. (2016). A modified SDS-based DNA extraction method for high quality environmental DNA from seafloor environments. Front. Microbiol. 7,986. doi: 10.3389/fmicb.2016.00986
Navon-Venezia, S., Kondratyeva, K., and Carattoli, A. (2017). Klebsiella pneumoniae: a major worldwide source and shuttle for antibiotic resistance. FEMS Microbiol. Rev. 41, 252–275. doi: 10.1093/femsre/fux013
Nordmann, P., Poirel, L., Toleman, M. A., and Walsh, T. R. (2011a). Does broad spectrum beta-lactam resistance due to NDM-1 herald the end of the antibiotic era for treatment of infections caused by Gram-negative bacteria? J. Antimicrob. Chemother. 66, 689–692. doi: 10.1093/jac/dkq520
Olaitan, A. O., Diene, S. M., Kempf, M., Berrazeg, M., Bakour, S., Gupta, S. K., et al. (2014). Worldwide emergence of colistin resistance in Klebsiella pneumoniae from healthy humans and patients in Lao PDR, Thailand, Israel, Nigeria and France owing to inactivation of the PhoP/PhoQ regulator mgrB: an epidemiological and molecular study. Int. J. Antimicrob. Agents 44, 500–507. doi: 10.1016/j.ijantimicag.2014.07.020
Piedra-Carrasco, N., Miguel, J., Fàbrega, A., Viñado, B., Campany, D., Mir, A., et al. (2018). Effectiveness of a double-carbapenem regimen in a KPC-producing Klebsiella pneumoniae infection in an immunocompromised patient. Microb. Drug Resist. 24, 199–202. doi: 10.1089/mdr.2017.0129
Pitout, J. D. D., Nordmann, P., and Poirel, L. (2015). Carbapenemase-producing Klebsiella pneumoniae, a key pathogen set for global nosocomial dominance. Antimicrob. Agents Chemother. 59, 5873–5884. doi: 10.1128/AAC.01019-15
Podschun, R., and Ullmann, U. (1998). Klebsiella spp. as nosocomial pathogens: epidemiology, taxonomy, typing methods, and pathogenicity factors. Clin. Microbiol. Rev. 11, 589–603. doi: 10.1128/CMR.11.4.589
Pontikis, K., Karaiskos, I., Bastani, S., Dimopoulos, G., Kalogirou, M., Katsiari, M., et al. (2014). Outcomes of critically ill intensive care unit patients treated with fosfomycin for infections due to pandrug-resistant and extensively drug-resistant carbapenemase-producing Gram-negative bacteria. Int. J. Antimicrob. Agents 43, 52–59. doi: 10.1016/j.ijantimicag.2013.09.010
Salabi, A. E., Walsh, T. R., and Chouchani, C. (2013). Extended spectrum beta-lactamases, carbapenemases and mobile genetic elements responsible for antibiotics resistance in Gram-negative bacteria. Crit. Rev. Microbiol. 39, 113–122. doi: 10.3109/1040841X.2012.691870
Santino, S. B. A. N. (2013). Microbiological and molecular characterization of extreme drug-resistant carbapenemase-producing Klebsiella pneumoniae isolates. Int. J. Immunopathol. Pharmacol. 26, 785–790. doi: 10.1177/039463201302600325
Schimmenti, A., Brunetti, E., Seminari, E., Mariani, B., Cambieri, P., and Orsolini, P. (2018). Prosthetic joint infection from carbapenemase-resistant Klebsiella pneumoniae successfully treated with ceftazidime-avibactam. Case Rep. Infect. Dis. 2018:1854805. doi: 10.1155/2018/1854805
Sharma, V. K., Johnson, N., Cizmas, L., McDonald, T. J., and Kim, H. (2016). A review of the influence of treatment strategies on antibiotic resistant bacteria and antibiotic resistance genes. Chemosphere 150, 702–714. doi: 10.1016/j.chemosphere.2015.12.084
Sugumar, M., Kumar, K. M., Manoharan, A., Anbarasu, A., and Ramaiah, S. (2014). Detection of OXA-1 beta-lactamase gene of Klebsiella pneumoniae from blood stream infections (BSI) by conventional PCR and In-Silico analysis to understand the mechanism of OXA mediated resistance. PLoS ONE 9,e91800. doi: 10.1371/journal.pone.0091800
Vading, M., Samuelsen, O., Haldorsen, B., Sundsfjord, A. S., and Giske, C. G. (2011). Comparison of disk diffusion, Etest and VITEK2 for detection of carbapenemase-producing Klebsiella pneumoniae with the EUCAST and CLSI breakpoint systems. Clin. Microbiol. Infect. 17, 668–674. doi: 10.1111/j.1469-0691.2010.03299.x
Wayne, P. A. (2018b). Clinical and Laboratory Standards Institute. Methods for Dilution Antimicrobial Susceptibility Tests for Bacteria That Grow Aerobically; Approved Standard-Eleventh Edition. Document M07-A11. CLSI, USA
Wiegand, I., Hilpert, K., and Hancock, R. E. (2008). Agar and broth dilution methods to determine the minimal inhibitory concentration (MIC) of antimicrobial substances. Nat. Protoc. 3, 163–175. doi: 10.1038/nprot.2007.521
Yu, F., Ying, Q., Chen, C., Li, T., Ding, B., Liu, Y., et al. (2012). Outbreak of pulmonary infection caused by Klebsiella pneumoniae isolates harbouring blaIMP−4 and blaDHA−1 in a neonatal intensive care unit in China. J. Med. Microbiol. 61, 984–989. doi: 10.1099/jmm.0.043000-0
Zhao, S., Duncan, M., Tomberg, J., Davies, C., Unemo, M., and Nicholas, R. A. (2009). Genetics of chromosomally mediated intermediate resistance to ceftriaxone and cefixime in Neisseria gonorrhoeae. Antimicrob. Agents Chemother. 53, 3744–3751. doi: 10.1128/AAC.00304-09
Keywords: Klebsiella pneumoniae, extensively drug resistance, antimicrobial resistance, multidrug resistant plasmid, high throughput sequencing, antimicrobial resistance gene
Citation: Li L, Yu T, Ma Y, Yang Z, Wang W, Song X, Shen Y, Guo T, Kong J, Wang M and Xu H (2019) The Genetic Structures of an Extensively Drug Resistant (XDR) Klebsiella pneumoniae and Its Plasmids. Front. Cell. Infect. Microbiol. 8:446. doi: 10.3389/fcimb.2018.00446
Received: 20 August 2018; Accepted: 13 December 2018;
Published: 04 January 2019.
Edited by:Alexandra Aubry, Sorbonne Universités, France
Reviewed by:Daniel Angel Ortiz, Laboratory Corporation of America Holdings (LabCorp), United States
Beiwen Zheng, Zhejiang University, China
Copyright © 2019 Li, Yu, Ma, Yang, Wang, Song, Shen, Guo, Kong, Wang and Xu. 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.