Emergence of a Multidrug-Resistant Hypervirulent Klebsiella pneumoniae Sequence Type 23 Strain with a Rare blaCTX-M-24-Harboring Virulence Plasmid

Here, we report a multidrug-resistant hypervirulent Klebsiella pneumoniae (MDR-HvKP) strain of sequence type 23 (ST23) with a rare hybrid plasmid harboring virulence genes and blaCTX-M-24, and we analyze the genetic basis for relationship between genotypes and MDR-hypervirulence phenotypes. Further analysis indicates that the hybrid plasmid is formed by IS903D-mediated intermolecular transposition of the blaCTX-M-24 gene into the virulence plasmid.

H ypervirulent Klebsiella pneumoniae (HvKP) is a serious threat to public health, as it could cause severe infection with high mortality and morbidity in young/healthy individuals (1). Previous studies showed that most HvKP strains are sensitive to common antibiotics, whereas multidrug-resistant (MDR) strains possess lower virulence (2). In addition to the differentiated phenotypes, HvKP and MDR K. pneumoniae show different genotypes with different clonal groups (3). However, recent studies found that several HvKP and MDR strains could evolve into MDR-HvKP strains through acquiring multidrug-resistant or hypervirulent plasmids (4)(5)(6)(7). The emergence and spread of MDR-HvKP is undoubtedly one of the most severe threat/challenges to global public health. Understanding of the genetic basis for MDR-HvKP strains is essential to control this deadly infection. In this study, we reported an MDR-HvKP strain of sequence type 23 (ST23) with a rare plasmid harboring virulence and bla CTX-M-24 genes, and we analyzed the genetic basis for relationship between genotypes and MDR-hypervirulence phenotypes.
MDR-HvKP strain 11492 was isolated from the blood of a male community infection patient (in his 30s) with an abscess of the kidney at the Chinese PLA General Hospital in November 2014. The patient was diagnosed with pancreatitis, pancreatic abscess, high fever (ϳ40.8°C), and shock (for detailed clinical information, see Table S1). The patient was subjected to antibiotic treatment (including imipenem-cilastatin sodium, teicoplanin, linezolid, meropenem, and biapenem), and operated on to remove the abscess.
We first evaluated the virulence phenotype of the strains using a string test. Strains 11492 and NTUH-K2044 showed a hypermucoviscous phenotype (viscous strings of Ͼ5 mm). We further analyzed the hypervirulent phenotype by wax moth larvae (Galleria mellonella; Huiyude Biotech Company, Tianjin) infection testing. G. mellonella larvae were maintained on wood chips in the dark at 15°C until use. We then infected the G. mellonella larvae with 10 4 , 10 5 , 10 6 , and 10 7 CFU/ml of K. pneumoniae strains. Experiments were performed in triplicate. The virulence of K. pneumoniae was determined by the survival rate of the G. mellonella larvae, using GraphPad Prism software 5.01 (10,11). The result indicated a significant decrease in survival rates when infected with 11492 and NTUH-K2044 HvKP strains in relative to those when infected with classical K. pneumoniae (cKP) strain 13190 under various infection concentrations at multiple time points (Fig. 1, Fig. S1 and S2). Thus, using HvKP strain NTUH-K2044 as the hypervirulence-positive control and cKP strain 13190 as the negative control, strain 11492 was demonstrated to be hypervirulent.
Drug susceptibility testing was performed by broth microdilution using 23 antibiotics, and the drug resistance phenotypes were determined according to CLSI standards (12 To reveal the genetic basis of the MDR-hypervirulence phenotype, we obtained the complete genome of MDR-HvKP 11492 by single-molecule real-time (SMRT) sequencing (Table S2 and  assembly was performed using the Hierarchical Genome Assembly Process 3 (HGAP3) within SMRT Analysis v2.3.0. Gap closing was completed by PBJelly (13), and circularization was achieved by manual comparison. To correct the polymer errors, we resequenced these isolates using Illumina sequencing (Table S3). Raw reads were trimmed, filtered, and mapped onto assembled genome sequences using BWA v0.5.9. Pilon v1.13 was subsequently employed to polish genome sequences using the obtained alignments (14). Genome sequences were annotated with the Rapid Annotations using Subsystem Technology (RAST) pipeline (15). Bioinformatic analysis provided the general chromosome genome information (Table S4), including GC% content (57.45%), genome size (5.25 Mb), predicted proteincoding genes (5,017), gene length (922 bp), and coding region (88.1%). Strain 11492 contained a 193-kb plasmid, p11492-vir-CTXM, with a lower GC content (50.37%), lower ratio of coding regions (75.71%), and shorter average gene length (724 bp).
Importantly, although most drug resistance and virulence genes are generally located in different plasmids (6), our analysis revealed a rare hybrid drug-resistant virulent plasmid, p11492-vir-CTXM (Table S5). It contains some known virulent genes (such as rmpA and rmpA2) and an extended-spectrum ␤-lactamase (ESBL) gene (bla CTX-M-24 ), contributing a lot to the MDR-hypervirulence phenotype. We then annotated p11492-vir-CTXM in detail. It was classified as a multireplicon IncHI1B/IncFIB plasmid using the PlasmidFinder database. We further performed a full-plasmid BLAST comparative analysis. The result showed that p11492-vir-CTXM exhibited 99% identity with the three IncHI1B/IncFIB virulent plasmids (pRJF999, GenBank accession no. CP014011; pK2044, GenBank accession no. AP006726; and pLVPK, GenBank accession no. AY378100), with high coverages (82%, 83%, and 76%, respectively) (Fig. 3A). Notably, p11492-vir-CTXM carried an ϳ7-kb region harboring a bla CTX-M-24 gene that was unique compared with the other three virulent plasmids. If we excluded the unique region, the backbone region of p11492-vir-CTXM would show high similarity to those of pRJF999 and pK2044, with 99% identity and 99% coverage, We further explored the genetic basis for the integration of unique region into the virulence plasmid. This region is inserted into a conserved backbone gene encoding the permease of the drug/metabolite transporter (DMT), with two IS903D elements at terminal regions (in opposite directions). Each IS903D is in turn linked to an external 8-bp sequence (GCACAGAGA), possibly a product of target site duplications, indicating the insertion event of the bla CTX-M-24 -embedded region. Additionally, this region is located in a novel transposon element, termed Tn6501, with the structure IS903D-bla CTX-M-24 -ISEcp1-ΔTn1721-ΔISKpn19-IS903D. Similar transposon genetic structure was also detected in the 68-kb plasmid pKP69 (GenBank accession no. EU195449.1) from Klebsiella pneumoniae, with 100% nucleotide identity and 84% coverage (Fig. 3B). These provided evidence of insertion of this region into the backbone of the virulence plasmid, based on IS903D-mediated intermolecular replicative transposition.
Although HvKP and MDR strains usually show nonoverlapping genotypes and phenotypes, our research demonstrated that a relatively infrequent MDR-HvKP strain evolved from HvKP through obtaining drug resistance genes. Furthermore, we discovered a rare hybrid plasmid, p11492-vir-CTXM, with both known virulent and bla CTX-M-24 genes , although drug resistance and hypervirulence genes are unlikely within the same plasmid. The emergence of MDR-HvKP strains, especially that of those carrying drug-resistant virulent plasmids, poses unprecedented threats/ challenges to public health. Mobile genetic elements may accelerate the formation of MDR-HvKP through horizontal gene transfer events. This is a dangerous tendency and should be closely monitored.
Accession number(s). Complete sequences of the chromosome of strain 11492 and of plasmid p11492-vir-CTXM have been deposited in the GenBank database under accession no. CP026021 and CP026022.