Next Article in Journal
16S-23S rRNA Internal Transcribed Spacer Region (ITS) Sequencing: A Potential Molecular Diagnostic Tool for Differentiating Lactococcus garvieae and Lactococcus petauri
Next Article in Special Issue
A Review of Carbapenem Resistance in Enterobacterales and Its Detection Techniques
Previous Article in Journal
Integrated Microbiota and Metabolome Analysis to Assess the Effects of the Solid-State Fermentation of Corn–Soybean Meal Feed Using Compound Strains
Previous Article in Special Issue
High-Molecular-Weight Plasmids Carrying Carbapenemase Genes blaNDM-1, blaKPC-2, and blaOXA-48 Coexisting in Clinical Klebsiella pneumoniae Strains of ST39
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 11
Issue 5
10.3390/microorganisms11051314
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Whole-Genome Sequencing Revealed the Fusion Plasmids Capable of Transmission and Acquisition of Both Antimicrobial Resistance and Hypervirulence Determinants in Multidrug-Resistant Klebsiella pneumoniae Isolates

by
Andrey Shelenkov
*,†,
Yulia Mikhaylova
†,
Shushanik Voskanyan
,
Anna Egorova
and
Vasiliy Akimkin
Central Research Institute of Epidemiology, Novogireevskaya Str., 3a, 111123 Moscow, Russia
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Microorganisms 2023, 11(5), 1314; https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms11051314
Submission received: 23 March 2023 / Revised: 11 May 2023 / Accepted: 15 May 2023 / Published: 17 May 2023
(This article belongs to the Special Issue Antimicrobial Resistance Mechanisms in Bacteria)
Browse Figures
Review Reports Versions Notes

Abstract

:
Klebsiella pneumoniae, a member of the Enterobacteriaceae family, has become a dangerous pathogen accountable for a large fraction of the various infectious diseases in both clinical and community settings. In general, the K. pneumoniae population has been divided into the so-called classical (cKp) and hypervirulent (hvKp) lineages. The former, usually developing in hospitals, can rapidly acquire resistance to a wide spectrum of antimicrobial drugs, while the latter is associated with more aggressive but less resistant infections, mostly in healthy humans. However, a growing number of reports in the last decade have confirmed the convergence of these two distinct lineages into superpathogen clones possessing the properties of both, and thus imposing a significant threat to public health worldwide. This process is associated with horizontal gene transfer, in which plasmid conjugation plays a very important role. Therefore, the investigation of plasmid structures and the ways plasmids spread within and between bacterial species will provide benefits in developing prevention measures against these powerful pathogens. In this work, we investigated clinical multidrug-resistant K. pneumoniae isolates using long- and short-read whole-genome sequencing, which allowed us to reveal fusion IncHI1B/IncFIB plasmids in ST512 isolates capable of simultaneously carrying hypervirulence (iucABCD, iutA, prmpA, peg-344) and resistance determinants (armA, blaNDM-1 and others), and to obtain insights into their formation and transmission mechanisms. Comprehensive phenotypic, genotypic and phylogenetic analysis of the isolates, as well as of their plasmid repertoire, was performed. The data obtained will facilitate epidemiological surveillance of high-risk K. pneumoniae clones and the development of prevention strategies against them.

1. Introduction

Klebsiella pneumoniae is a significant hospital pathogen that can cause infections of varying severity in the bloodstream, respiratory and nervous systems, and urogenital tract of a human or animal body, often leading to serious adverse events [1]. This Gram-negative bacterium is one of the most widespread and dangerous members of the Enterobacteriaceae family, and it has become one of the major causes of healthcare-associated infections worldwide [2].
Two independent evolutionary branches (pathotypes) have been considered to represent the K. pneumoniae population: classical (cKp) and hypervirulent Klebsiella (hvKp) [3,4]. CKp is distributed globally and constitutes a fraction of the opportunistic flora of healthy people. Such strains can rapidly acquire multiple resistance to antibiotics in clinical settings, which makes them a leading causative agent of nosocomial infections. A subgroup of cKp is designated as “high-risk clones” representing distinct genetic lineages that usually are multidrug-resistant (MDR) or extensively drug-resistant (XDR). Such isolates have spread worldwide and have caused large numbers of infections, mostly in clinical settings [5,6]. The best-known and most-characterized K. pneumoniae high-risk clones belong to multilocus sequence typing (MLST)-based types ST11, ST14, ST15, ST17, ST45, ST147, ST258, ST307 and ST512 [3].
HvKp is associated with more severe and aggressive infections, usually in healthy and community-residing individuals [7,8]. HvKp, more often than cKp strains, carries genetic factors associated with a hypermucoid phenotype and with the synthesis of lipo- and polysaccharide capsules, as well as genes encoding the elements of an iron acquisition system, allantoin utilization system, and type I fimbriae [8,9,10]. Such strains are predominantly associated with clones ST23, ST65 and ST86 [11].
CKp and hvKp were previously considered as non-overlapping evolutionary branches due to their significant differences [12]. However, during the last fifteen years, several research studies have demonstrated the acquisition of multidrug resistance by hypervirulent strains of Klebsiella in South Korea [13], Argentina [14], France [15], China [16] and the Russian Federation [17,18]. Moreover, a case of the acquisition of a virulence plasmid by a classical extremely resistant hospital strain of K. pneumoniae in China has been described [19]. This confirmed the concerns of some researchers regarding the possibility of a new “superpathogen” formation combining the pathogenetic properties of the two K. pneumoniae evolutionary branches [20].
Hypervirulence, as well as multidrug resistance, is usually associated with the acquisition of additional genetic material and the formation of genetic lines that effectively support these acquired determinants [21]. The ecology of K. pneumoniae and its genome plasticity designated this species as a key reservoir and distributor of antibiotic resistance genes. This pathogen is more competent to acquire genetic material of plasmid origin from different donors of different ecological niches [22]. This bacterium is also capable of keeping such plasmids long enough to pass them on to new recipients living in humans and animals [5]. K. pneumoniae often contains more than one plasmid in the genome, including both low-molecular-weight, high-copy-number, and high-molecular-weight, low-copy-number, plasmids. High-molecular-weight plasmids belong to different incompatibility groups and usually carry genetic determinants of virulence and antibiotic resistance [23]. In recent years, whole-genome sequencing (WGS), especially long-read sequencing, has become a powerful tool for the determination of plasmid structures, as well as antimicrobial resistance (AMR) and virulence gene locations [17,24,25].
In this study, we described the detailed phenotypic and genomic characteristics of four MDR K. pneumoniae isolates based on antimicrobial susceptibility testing, hybrid whole-genome assembly and phylogenetic comparison. The analysis performed allowed us to assess thoroughly the plasmid composition of the samples studied and to determine the localization of AMR and virulence determinants. We also succeeded in revealing the convergence of hvKP and MDR-cKp traits in one of the isolates, which resulted in the emergence of a strain that could pose a serious threat to the healthcare system. The data obtained will be useful for epidemiological surveillance of this important pathogen, especially in the Russian Federation, where it has become the major source of nosocomial infections in recent years (AMRmap database [26], https://amrmap.ru/, accessed on 1 March 2023).

2. Materials and Methods

2.1. Sample Collection, Susceptibility Testing and DNA Isolation

Three K. pneumoniae isolates were collected from different patients of the same medical department (ICU) during the short intervals (3–4 days) of their hospital stays. One K. pneumoniae sample was isolated from a patient in the pulmonology department during the same period. This patient was discharged from the hospital and the other three patients died. More detailed metadata are presented in Table 1.
The total number of bacterial isolates analyzed in the tertiary care hospital per year in the pilot study was approximately 850, of which 160 were sequenced and about 80% represented K. pneumoniae. The reasons for choosing these particular isolates for investigation were the severity of infections they caused, their resistance to several classes of antibiotics including carbapenems and their possible carriage of hybrid plasmid replicons revealed by PCR.
Species identification was performed using time-of-flight mass spectrometry (MALDI-TOF MS) with the VITEK MS (bioMerieux, Marcy-l’Étoile, France). Antimicrobial susceptibility/resistance of the isolates was determined by the disc diffusion method using the Mueller–Hinton medium (bioMerieux, Marcy-l’Étoile, France) and disks with antibiotics (BioRad, Marnes-la-Coquette, France), and by the Minimum Inhibitory Concentration (MIC) method using a VITEK 2 Compact 30 analyzer (bioMerieux, Marcy-l’Étoile, France).
The antibiotics tested included amikacin, amoxicillin, ampicillin, cefepime, cefotaxime, ceftazidime, ceftriaxone, ceftazidime/avibactam, cefuroxime, ciprofloxacin, ertapenem, gentamicin, imipenem, levofloxacin, meropenem, piperacillin/tazobactam, tobramycin and trimethoprim/sulfamethoxazole. EUCAST clinical breakpoints, version 11.0 (https://www.eucast.org/clinical_breakpoints/, accessed on 20 December 2021), were applied to interpret the susceptibility/resistance results obtained.

2.2. Whole-Genome Sequencing

Genomic DNA was isolated with the DNeasy Blood and Tissue kit (Qiagen, Hilden, Germany); then, the Nextera™ DNA Sample Prep Kit (Illumina®, San Diego, CA, USA) was applied to prepare paired-end libraries for WGS of the isolates on an Illumina® NextSeq 2000 platform (Illumina®, San Diego, CA, USA).
The same extracted genomic DNA was used for library preparation for the Oxford Nanopore MinION sequencing system (Oxford Nanopore Technologies, Oxford, UK) with the Rapid Barcoding Sequencing kit SQK-RBK004 (Oxford Nanopore Technologies, Oxford, UK). The amount of initial DNA was 400 ng for each sample. The libraries were prepared according to the manufacturer protocols and were sequenced on FLO-MIN106 R9.4 flow cell with a standard 24 h sequencing protocol using the MinKNOW software version 22.03 (Oxford Nanopore Technologies, Oxford, UK).

2.3. Genome Assembly, Data Processing and Annotation

Basecalling of the raw MinION data was made using Guppy basecalling software version 6.1.7 (Oxford Nanopore Technologies, Oxford, UK), and demultiplexing of the samples was performed using Guppy barcoding software version 6.1.7 (Oxford Nanopore Technologies, Oxford, UK).
Hybrid short- and long-read assemblies were obtained using Unicycler version 0.5.0 (normal and bold mode) [27]. Additional plasmid assembly was performed using Flye software, version 2.9.1 [28], using default parameters, and assembly polishing and correction was made with Polypolish version 0.5.0 [29] and Medaka version 1.7.3 (https://github.com/nanoporetech/medaka, accessed on 1 March 2023) using the ‘r941_min_fast_g507’ model.
Genome assemblies were submitted to NCBI Genbank under the project PRJNA942929.
The pipeline described earlier [30,31] was used for the raw data filtration and assembled genome processing and annotation. Resfinder 4.3.0 software was used for antimicrobial gene detection (https://cge.cbs.dtu.dk/services/ResFinder/, accessed on 20 February 2023, using default parameters) and point mutation analysis. Virulence factors were searched in VFDB (http://www.mgc.ac.cn/VFs/main.htm, accessed on 20 February 2023, using default parameters). Additional virulence genes (peg-344, etc.) were searched using BLASTn using reference sequences from Genbank (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/genbank/, accessed on 12 March 2023).
Plasmid replicon types, the predicted mobility and other annotation features were detected using MOB-suite [32] using default parameters (mob_typer version 3.1.2). Plasmid visualization was made using BRIG (version 0.95, https://github.com/happykhan/BRIG, accessed on 12 March 2023). CRISPRCasFinder [33] was used to identify the presence of CRISPR/Cas systems in the genomes studied. Spacers were analyzed using the CRISPRTarget tool [34].
Isolate typing was performed by MLST using BIGSdb (https://bigsdb.pasteur.fr/klebsiella/, accessed on 20 February 2023). In addition, the types based on capsule synthesis loci (K-loci) and lipooligosaccharide outer core loci (OCL) were detected using Kaptive software version 2.0.3 with default parameters [35]. Virulence scores for the isolates were assessed with Kleborate [36] version 2.3.2.
Detection of cgMLST profiles was performed using MentaList (https://github.com/WGS-TB/MentaLiST, version 0.2.4, default parameters) using the scheme obtained from cgmlst.org (https://www.cgmlst.org/ncs/schema/schema/2187931/, contained 2358 loci, last update 20 February 2023) [37]. The minimum spanning tree was built using PHYLOViz online (http://online.phyloviz.net, accessed on 20 February 2023).

3. Results

3.1. Isolate Typing and Resistance Profiles

In silico MLST revealed that three isolates (CriePt492, 494 and 495) belonged to ST512 and they were characterized as possessing KL107 capsular type and O2v2 oligosaccharide type. The typing results showed another profile for the isolate CriePt491—ST147/KL20/O2v1. It is worth noting that the determined genetic lines belong to the “high-risk clones” of K. pneumoniae mentioned above.
According to phenotypic analysis, all isolates were susceptible to ceftazidime/avibactam and resistant to the other antimicrobial compounds tested. The isolates CriePt492 and CriePt495 were additionally susceptible to gentamycin. Thus, all the isolates under investigation exhibited multidrug resistance (Figure 1).
CgMLST analysis revealed that the genomic sequences of ST512 isolates were very similar (3 allele differences between CriePt492 and 495, and 29 different alleles between either of these isolates and CriePt494). Thus, according to the criterion described previously (≤18 cgMLST allele differences between the isolates to be considered as belonging to a single strain or clone for K. pneumoniae [38]), the isolates CriePt492 and 495 were highly likely representing a single strain; CriePt494 demonstrated the beginning steps of divergence from its relative samples.
In silico searching for antimicrobial resistance determinants showed that chromosomal genes (oqxAB and fosA) were represented in all the isolates under investigation. In addition, fosA and blaSHV-182 were found on the chromosomes of CriePt492 and 495, and only fosA (without blaSHV, which was found on a plasmid) on the chromosome of CriePt494. In contrast, CriePt491 had the major part of its AMR genes on its chromosome.
Eleven antibiotic-resistance genes (ARGs) were found in both ST512 isolates, CriePt492 and 495, and their resistance gene profiles were completely the same. These isolates carried two aminoglycoside resistance genes, two beta-lactamase genes (blaKPC-3 and blaSHV-182) and other genes associated with resistance to chloramphenicol (cat), macrolides (mph(A)), fosfomycin (fos), fluoroquinolones (oqx), trimethoprim (dfr) and sulfonamides (sul). Another ST512 isolate, CriePt494, differed from its relatives by having additional determinants, namely, two additional aminoglycoside resistance genes, two beta-lactamase genes (blaCTX-M15 and blaNDM-1), as well as one more each of macrolide-, fluoroquinolone- and trimethoprim-resistance genes. The composition of the antimicrobial-resistance genes of the ST147 isolate (CriePt491) was more similar to the isolate CriePt494, and it carried one additional beta-lactamase gene, blaTEM-1A, thus possessing four different beta-lactam resistance genes in the genome. It also had one more gene associated with resistance to macrolides—msr (E).
The differences also included the armA gene, which was revealed only in CriePt491 (chromosome) and CriePt494 (plasmid), and which corresponded to their phenotypic resistance to gentamicin.
Point mutation analysis revealed porin mutations associated with carbapenem resistance in all isolates. OmpK36 mutation A217S and ompK37 mutations I70M and I128M were revealed in all isolates, and ompK37 N230G was additionally found in ST512 isolates. Several mutations conferring resistance to fluoroquinolones were also revealed in the acrR gene for all isolates studied. A summary of point mutations revealed is provided in Table S1.
Annotation results of the resistance determinants for the isolates studied were consistent with the resistance phenotypes for the antimicrobial drugs included in the panel.

3.2. Isolate Typing and Resistance Profiles

The sets of virulence factors in K. pneumoniae isolates under investigation were very similar and included up to 68 virulence genes. The most important of them, which were revealed in all four isolates, comprised acridine efflux pump genes (acr), fimbria involved in the processes of adhesion, formation of biofilms, other regulators of the initial stages of infection (fimABCDEFGHIK, mrkABCDFHIJ, ecpRABCDE, RscAB), genes encoding outer membrane, ATP-binding proteins and regions associated with type IV secretion system and enzymes of catalytic activity (dotU, ompA, clpV, tssF/G, ugd, gnd, etc.), as well as several siderophores (entABCEF, fepABCDG) and genes of iron uptake (salmochelin and aerobactin) clusters (iutA and iroE, respectively).
Two isolates (CriePt492 and 495) belonging to the single strain had exactly the same set of virulence factors. The differences were observed in the genome of the isolate CriePt494. It possessed the aerobactin cluster (iucABCD) and its receptor (iutA), and the capsule upregulation gene prmpA2 (where subscript p indicates its plasmid origin), a regulator of a mucoid phenotype. The plasmid CriePt494_p2 also contained peg-344 metabolite transporter, which is considered one of the markers of hypervirulent phenotype [9]. The isolate CriePt491 differed by the presence of the rfbK1 gene encoding O9 family phosphomannomutase involved in the synthesis of capsular polysaccharide.
The presence/absence for a selected set of the most important virulence genes for the isolates studied is shown in Figure 2. A complete list of virulence factors is given in Table S2.

3.3. Plasmid Typing, Classification, Annotation and Comparison

The isolates studied differed from each other by plasmid number and composition. A brief description of replicon typing is given in Table 2, while the extended information for each plasmid assembled is shown in Table 3. Hybrid assembly allowed us to determine the refined plasmid structure and obtain the data regarding AMR and virulence gene location, which represents an important step in elucidating the possible mechanisms of their transfer.
CriePt492 and CriePt495 both included six plasmids with the same repertoire of replicons, while CriePt494 possessed an additional IncFIA plasmid. CriePt491 had the lowest number of plasmids (five) and did not possess the ColRNAI replicon. However, we revealed more differences in plasmid structure with a more detailed analysis. CriePr492, 494 and 495 shared the four smallest plasmids with lengths ranging from 2963 to 53,292, and their corresponding plasmids had the same or very close sequences. At the same time, the plasmids of CriePt491 were completely different despite having similar replicon types. Notably, the former three isolates had a small non-typeable plasmid with a length of 2963, while CriePt491 included a different, yet also non-typeable, plasmid with a length of 53,009.
Next, the plasmids of CriePt492, 494 and 495 will be discussed together since they shared many features, and the plasmids of CriePt491 will then be discussed separately.

3.3.1. Plasmids of CriePt492, 494 and 495

The smallest non-typeable plasmid (length = 2963) had a 99% similarity with the pR17.4849_3.0k plasmid from Salmonella enterica obtained in Taiwan (Genbank acc. CP100750), which also did not have a commonly known replicon type and was attributed to the AG463 cluster by the MOB-typer tool. This plasmid was also sequenced on MinION, and sequence mismatches probably reflected the differences in the assembly algorithms used. The relaxases of MOBV type were detected in these plasmids, and they were mobilizable, which means that they rely on conjugative plasmids to provide the mating pair formation components, and MPFT plasmids can serve as helpers for their conjugation [39]. No AMR or virulence genes were found in these plasmids.
The IncX3 plasmids contained the gene of class A beta-lactamase blaKPC-3 known to provide resistance to various antibiotics including third-generation cephalosporins and carbapenems for K. pneumoniae, as well as the blaSHV-182 beta-lactamase gene, which was also revealed in the chromosomes. The plasmids were found to be conjugative and included MOBP relaxase. In addition, they were attributed to MPFT type, which is based on the plasmids’ T4SS system, involved in mating pair formation (MPF) during conjugation [40]. The closest Genbank match for this plasmid was the Escherichia coli strain SUISSEKPC3NDM5 plasmid p1606b (Genbank acc. CP083703.1, Switzerland), which shared higher than 99% identity with our plasmids, contained the same AMR genes and had the same total length.
The plasmids with the Col(pHAD28) replicon carried the single AMR gene aac(6′)-Ib, providing resistance to aminoglycosides. Col(pHAD28) plasmids had almost identical sequences in ST512 isolates and carried the aac(6′)-Ib gene providing aminoglycoside resistance. No relaxases were detected, and the plasmids were considered as non-mobilizable. The closest match from Genbank for these plasmids was pCR14_5 from K. pneumoniae (Genbank acc. CP015397.1, USA), which had a length of 9456 and included the same AMR gene.
The ColRNAI plasmids also had nearly identical sequences for the isolates studied. AMR genes and virulence factors were not revealed in them. The plasmids included MOBC relaxase and were considered mobilizable. The K. pneumoniae strain BA4656 plasmid pColRNAI (Genbank acc. CP035911.1, length = 9729, India) had a higher than 99% similarity with this set of plasmids. The reference plasmid also did not include any virulence or resistance determinants.
CriePt494 possessed an additional conjugative IncFIA plasmid with a relaxase of MOBF type. It carried important AMR determinants, such as the gene encoding class A extended-spectrum beta-lactamase (ESBL), blaCTX-M-15, which can provide resistance to various beta-lactam antibiotics including fourth-generation cephalosporins. This gene was not found in other ST512 isolates, and was revealed in CriePt491, but in a plasmid of different type. Additional AMR genes located on this IncFIA plasmid included erm(B), aac(3)-IIa and qnrS1, providing resistance to lincosamides/macrolides, aminoglycosides and fluoroquinolones, respectively. A similarity among reference Genbank plasmids was found in the Shigella sonnei plasmid p6904-27 (CP045525.2, Switzerland, length = 83,273) and Escherichia coli strain SUISSEKPC3NDM5 plasmid p1606d (CP083705.1, Switzerland, length = 146,558), which also included the blaCTX-M-15 gene.
The IncFIB plasmids of CriePt 492 and 495 were very similar and included MOBF relaxase. However, they differed by one AMR gene—aac(6′)-Ib—which was revealed in a plasmid from CriePt495, but not in CriePt492. Other AMR genes were the same for both plasmids and included aadA2, aph(3′)-Ia, catA1, dfrA12, mph(A) and sul1, thus empowering these isolates with a broad range of weapons against various classes of antimicrobial drugs. In contrast, the IncFIB plasmid of CriePt494 was almost twice as short and did not include any AMR genes, while most of the genes given above were revealed in the IncHI1B plasmid of this isolate. In addition, the CriePt494 plasmid did not possess any relaxases. CriePt494_p3 was similar to the pIT-394-FIB K. pneumoniae plasmid (CP110954.1, Italy), which had almost the same length and did not include any AMR genes either. CriePt492_p3 was similar to the pKPN-a68 plasmid (CP009777.1, USA, length = 212,192), which included the same AMR genes.
Meanwhile, the main differences inside the group of ST512 isolates (CriePt 492, 494 and 495) were revealed in the large fusion conjugative plasmids IncHI1B/IncFIB (p2). All three plasmids were attributed to MOBH mobility type and to MPFF type. However, while the plasmids from CriePt492 and 495 were almost completely similar and did not carry any AMR genes or virulence factors, the differences between CriePt494_p2 and its counterparts were dramatic. This plasmid included both AMR genes, including the New Delhi metallo-beta-lactamase-1 gene and other genes providing resistance to several classes of antibiotics, and virulence factors iucABCD and rmpA2. As was mentioned above, CriePt494 was the only isolate carrying these hypervirulence factors. In addition, it also included peg-344. This plasmid was similar to a fusion plasmid, phvKpST147_NDM-1_2566 (MW911671.1), revealed recently in hypervirulent K. pneumoniae from Russia, but it had a slightly different repertoire of virulence/resistance genes. All other virulence genes revealed in the isolates studied were found in the chromosomes.
A circular diagram representing the IncHI1B/IncFIB plasmids from CriePt isolates and the phvKpST147_NDM-1_2566 reference plasmid is shown in Figure 3. It is easily seen that the plasmids of CriePt 491, 492 and 495 lack the parts containing resistance and virulence genes, and phvKpST147_NDM-1_2566 lacks the regions containing aph(3′)-Ia, dfrA12 and aadA2 genes from CriePt494.

3.3.2. Plasmids of CriePt491

Although CriePt491 was also MDR, as were its counterparts, a large fraction of the AMR genes in this isolate were located in its chromosome, while additional resistance determinants were revealed in the IncR plasmid not possessed by ST512 isolates. The plasmids of CriePt491 were mostly non-mobilizable and did not include any relaxases.
The large IncHI1B/IncFIB plasmid did not include any AMR genes or virulence factors. It had a region of more than 100,000 in length almost identical to the phvKpST147_NDM-1_2566 plasmid of K. pneumoniae (MW911671.1, Saint Petersburg, Russia), but multiple AMR genes and virulence factors in this 350k reference plasmid were located outside of this region.
The 100k IncFIB plasmid did not include any pathogenicity determinants either, but it possessed the gene encoding the RepB replication initiation protein. The closest match for this plasmid in Genbank was a phvKpST147_3 plasmid (CP066857.1, Saint Petersburg, Russia) with 99.9% identity, which also did not possess any AMR or virulence determinants.
The 53k p4 was similar to a plasmid from a Norwegian K. pneumoniae isolate of approximately the same length (CP034049.1), for which no known replicons were revealed either.
The small 5k Col(pHAD28) plasmid was similar to a P7 plasmid from a Spanish K. pneumoniae isolate (OW970547.1).
Finally, the plasmid of particular interest was a 47k IncR plasmid, which carried multiple AMR genes. The plasmid included three beta-lactamases, namely, class A blaCTX-M-15 and ESBL blaTEM-1B, as well as class D blaOXA-9, which provided resistance to multiple beta-lactam antibiotics, including fourth-generation cephalosporins. In addition, aminoglycoside- and ciprofloxacin-resistance genes were detected. Together with the chromosome-encoded carbapenemase blaNDM-1 and other resistance genes, these plasmids provided CriePt491 with a multifunctional antibiotic defense system. However, the IncR plasmid did not include any relaxases or mobilization proteins, and thus it requires additional helper plasmids to be transferred via a conjugation process. This plasmid was similar to the EFN 299 plasmid p4 (CP092593.1, Ghana) and phvKpST147_4 (CP066858.1, Saint Petersburg, Russia), which also carried multiple AMR genes including beta-lactamase genes.

3.4. CRISPR-Cas Systems and Anti-CRISPR Genes

A type I-E CRISPR-Cas system was identified in the CriePt491 chromosome (start = 896,591; end = 905,034) that included intact and apparently functional cas2, cas1, cas6, cas5, cas7, cse2, cse1 and cas3 genes. No reliable CRISPR sequences were identified in other isolates. The scheme of the system is shown in Figure 4.
Spacer analysis for the CriePt491 CRISPR-Cas system revealed a large fraction (about 70%) of anti-plasmid sequences, in particular, to IncFII and IncX3 plasmids. However, detailed analysis of these sequences lies beyond the scope of the manuscript.
Remarkably, the genes encoding the AcrIIA7 Cas9 inhibitor protein were found with >90% amino acid sequence identity in the chromosomes of all isolates except CriePt491.

3.5. Phylogenetic Analysis of the Isolates

We compared the ST512 isolates and ST147 CriePt491 separately with the corresponding reference isolates from Genbank (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/genbank/, accessed on 10 February 2023) based on cgMLST analysis. The minimum spanning tree for the ST512 isolates and the best matches from Genbank are shown in Figure 5.
Pairwise comparison of cgMLST profiles according to the threshold set previously (≤18, [38]) revealed two isolates that were within the single strain range with CriePt492 and CriePt495—strain 52/1 (GCA_022555555 from Samara, Russia, 17 allele differences), and 2021CK-01402 (GCA_021904695.1, USA, 11 allele differences), while the closest neighbor of CriePt494 (again, GCA_021904695.1) had 26 allele differences, thus exceeding the threshold value. Both reference isolates contained various AMR genes, including blaKPC-3, but did not possess the hypervirulence determinants iucABCD and rmpA2. No plasmid structures were determined for these reference isolates by their uploaders, which made direct comparison impossible.
The closest neighbor of CriePt491 was the isolate 2566 (GCA_011044915, Saint Petersburg, Russia, 20 allele mismatches), which possessed the fusion plasmids mentioned above. This isolate had blaCTX-M-15 and blaNDM-1 beta-lactamase genes, which were also possessed by CriePt491. Interestingly, 2566 also possessed an apparently functional type I-E CRISPR-Cas system.
CgMLST profiles for the isolates studied and the corresponding reference isolates are presented in Table S3.

4. Discussion

K. pneumoniae has already become a widespread healthcare-associated pathogen constituting a significant threat to public health due to its ability to rapidly acquire and spread AMR determinants and to develop resistance even to novel antibiotics [22,41]. Several lineages of “classical”, often MDR or XDR, K. pneumoniae strains, which have usually infected chronically ill patients residing in hospitals, were designated as high-risk clones due to their increased morbidity and mortality [3]. The best-characterized of these K. pneumoniae clones include, among others, ST147 and ST512 [3,42], to which the isolates described in our study belong. Although previously high-risk clones and hypervirulent clones of K. pneumoniae were considered to be separate evolutional groups, an increasing number of reports worldwide, especially in the last 5 years, have described the convergence of hypervirulence and MDR properties in single isolates of this species [17,19,43,44,45].
In this study, we analyzed three MDR ST512 isolates and one MDR ST147 isolate, which caused severe adverse conditions in patients admitted to a tertiary care hospital in which three out of the four patients died. Two ST512 isolates, CriePt492 and CriePt495, were almost identical both in terms of their genomic structure and in terms of their phenotypic characteristics, and the last ST512 isolate, CriePt494, was close to them in cgMLST profile but had a different plasmid repertoire. CriePt494 was the only isolate containing peg-344, iucA and the plasmid-borne rmpA2 gene, which are considered reliable diagnostic markers for identifying strains in the hvKp-rich cohort [9,46], and it was assigned a virulence score of three on a five-point scale developed earlier [36]. All these genes, together with other virulence factors from the iuc cluster and several AMR genes, including blaNDM-1, were contained in a large 300k fusion conjugative plasmid, IncHI1B/IncFIB. A recent comprehensive in silico analysis revealed that about 70% of the known virulence plasmids from Genbank carried hybrid replicons of this type, and it specifically reported the genes iucA and iutA, which were found in the CriePt494_p2 plasmid, as being associated with high virulence [47]. The structure of this plasmid was similar to phvKpST147_NDM-1_2566 previously found in Russia [43] (the isolates with similar plasmids in that study belonged to ST15, ST147, ST395 and ST874), but their AMR gene content was different. Both CriePt492 and CriePt495 also had a conjugative plasmid IncHI1B/IncFIB each, but in this case, there were no AMR or virulence genes in these plasmids. The same was true for the IncHI1B/IncFIB plasmid of ST147 CriePt491, but it also lacked any relaxases and was predicted to be non-mobilizable. Thus, we can conclude that the possession of a fusion plasmid was not a property of the isolates having particular ST, and we reported its carriage in ST512, which had not been revealed previously by other researchers [43].
The formation of K. pneumoniae genetic lineages exhibiting both high virulence and carbapenem resistance can occur by several pathways associated with horizontal gene transfer [43,48], and the conjugation was considered its most significant mechanism [49]. Hypervirulent isolates can acquire MDR plasmids, and MDR isolates can also acquire plasmids carrying virulence factors. Some reports have supposed that the latter is more likely [21], but others have confirmed that the former occurs more often, or that both are possible [10,50], and the involvement of IncHI1B/FIB was revealed in both situations [51]. In our case, we can suppose the acquisition of additional virulence genes by the MDR isolate CriePt494 to form the fusion virulence/resistance-conferring plasmid since we revealed CriePt492 and 495, which (i) also contained IncHI1B/FIB, but without AMR or virulence factors, and (ii) included another IncFIB containing AMR, but not virulence, genes similar to the ones from IncHI1B/FIB of CriePt495, while the corresponding IncFIB from the latter did not contain any pathogenic factors. However, additional intermediate plasmids should be revealed to confirm the possible evolutionary pathway. Meanwhile, the formation of such mosaic plasmids by possible stepwise acquisition of transposon elements has been hypothesized in previous reports [43,52,53].
In general, the isolates studied showed great correspondence of their phenotypic and genomic resistance profiles. The resistance to gentamicin exhibited by CriePt491 and CriePt494 was highly likely conferred by the presence of the armA gene on the chromosome and the IncHI1B/FIB plasmid, respectively. Carbapenem resistance of the isolates could be attributed to NDM-1 and KPC-3 carbapenemases, as well as to mutations revealed in ompK36 and ompK37 porins, which were reported previously for carbapenem-resistant K. pneumoniae isolates from Turkey [54] also producing KPC-3, and from Nigeria [55]. Notably, CriePt491 had its blaNDM-1 gene in a chromosome, which is quite unusual for K. pneumoniae and was previously associated with the integration from an IncHI1B-like plasmid [56]. All isolates exhibited a high level of AMR, and CriePt491 and CriePt494 were susceptible only to a ceftazidime/avibactam combination from the panel of the antibiotics tested.
CriePt491 was the only isolate not having conjugative plasmids and carrying a large fraction of its AMR genes on the chromosome, and these facts could be explained by its possession of a complete CRISPR-Cas type I-E system. The situation of chromosomal AMR gene location is not typical for K. pneumoniae since this species usually carries most of its acquired resistance genes on the plasmids [2,25,57], which was observed in other isolates from the current study. At the same time, all other isolates included the AcrIIA7 Cas9 inhibitor protein gene, which could be the reason why they did not have CRISPR-Cas, although AcrIIA7 has been supposed to target II-A systems [58]. The type I-E CRISPR system was previously suggested to efficiently limit the acquisition of antibiotic resistance genes and plasmids in K. pneumoniae strains [59,60,61]. Nevertheless, other recent findings have shown that the presence of a putative CRISPR-Cas system in a bacterial genome is not a limitation for the acquisition of common MDR plasmids by K. pneumoniae isolates [17,62], and some plasmids carrying resistance genes can even simultaneously contain a CRISPR-Cas system against other types of plasmids [53,63].
The targeting of CRISPR-Cas systems is achieved by incorporating the DNA fragments (protospacers) from infecting bacteriophages or plasmids to the CRISPR array, where they form a sequence-specific memory against subsequent invasions by the same bacteriophage or plasmid [63]. For example, I-E CRISPR-Cas-matched protospacers have commonly been found in conjugative IncFII-MDR plasmids [60], which play an important role in AMR dissemination within K. pneumoniae populations [64]. InxX3 plasmids also have a great conjugation ability and are more commonly associated with carbapenemase production than, e.g., IncR or IncH plasmids [65]. In our study, it was in fact the IncX3 plasmid that carried the blaKPC-3 gene in ST512 isolates. The analysis of spacer sequences in the CriePt491 CRISPR-Cas system revealed a large number of spacers targeting IncFII and IncX3 plasmids, which this isolate did not possess. However, this purely correlative observation needs to be verified by additional experiments. Interestingly, a recent study has demonstrated the possibility of using the re-engineered conjugative endogenous CRISPR-Cas3 system to cure the AMR gene carrying K. pneumoniae IncFII plasmids in a Galleria mellonella infection model [64]. Thus, future investigations in this field seem to be very promising for tackling AMR in K. pneumoniae and other bacterial species.
The limitation of our study lies in the small number of the isolates investigated. However, the aim of this report was not to give a snapshot of K. pneumoniae spreading in a particular hospital or region, but rather to draw attention to the emergence of fusion plasmids carrying both virulence and resistance genes and to look for clues regarding the possible mechanisms of their formation. Prompt accumulation of such information is very important, and previous reports have revealed that timely and appropriate analysis of bacterial sequencing data could be a bottleneck in epidemiological surveillance [66,67]. We plan to continue the surveillance of MDR K. pneumoniae isolates and their plasmids in the future when more samples become available.
We believe that the plasmid data provided in this manuscript will facilitate the epidemiological surveillance of K. pneumoniae in Russia and other countries and that it presents useful data for developing prevention measures against this prominent pathogen. Such measures could possibly include the development of novel antibiotics [68] or other antimicrobial compounds like antimicrobial peptides [69,70], implementation of bacteriophages [71], and advancing antibiotic stewardship [72].

5. Conclusions

We have performed long- and short-read WGS and conducted comprehensive plasmid analysis of four clinical MDR K. pneumoniae isolates and we revealed a hybrid IncHI1B/FIB plasmid carrying both AMR and virulence genes in one of the isolates, CriePt494, belonging to a high-risk clone, ST512. This plasmid carried peg-344, iucABCD and prmpA2 high-virulence determinants, as well as multiple resistance genes, including blaNDM-1 beta-lactamase. Other isolates also possessed plasmids of this type, but they did not carry any pathogenicity determinants, although their structure allows assuming that they are ‘ready’ to obtain AMR or virulence genes in the future. The fact that such plasmids were revealed in a rather small collection of clinical isolates suggests the possibility of their spread rate increasing and draws attention to the threat posed by the increasing risk of fusion AMR/virulence plasmid formation in clinical settings.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/microorganisms11051314/s1, Table S1: The summary of point mutations revealed in clinical K. pneumoniae isolated under study; Table S2: The list of virulence factors revealed in clinical K. pneumoniae isolates studied; Table S3: cgMLST profiles for clinical K. pneumoniae isolates studied and reference isolates from Genbank.

Author Contributions

Conceptualization, Y.M.; data curation, A.S.; formal analysis, A.S. and Y.M.; funding acquisition, V.A.; investigation, A.S., Y.M., S.V. and A.E.; methodology, A.S. and Y.M.; project administration, V.A.; resources, Y.M., S.V. and A.E.; software, A.S.; supervision, V.A.; validation, A.S. and Y.M.; visualization, A.S., Y.M. and A.E.; writing—original draft, A.S. and Y.M.; writing—review and editing, A.S. and Y.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by a grant from the Ministry of Science and Higher Education of the Russian Federation (agreement No. 075-15-2019-1666).

Institutional Review Board Statement

Ethical review and approval were waived for this study since the human samples were routinely collected and patients’ data remained anonymous.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The assembled genome sequences for all isolates were uploaded to the NCBI Genbank under the project number PRJNA942929.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Magill, S.S.; Edwards, J.R.; Bamberg, W.; Beldavs, Z.G.; Dumyati, G.; Kainer, M.A.; Lynfield, R.; Maloney, M.; McAllister-Hollod, L.; Nadle, J.; et al. Multistate point-prevalence survey of health care-associated infections. N. Engl. J. Med. 2014, 370, 1198–1208. [Google Scholar] [CrossRef] [PubMed]
  2. Gorrie, C.L.; Mirceta, M.; Wick, R.R.; Judd, L.M.; Lam, M.M.C.; Gomi, R.; Abbott, I.J.; Thomson, N.R.; Strugnell, R.A.; Pratt, N.F.; et al. Genomic dissection of Klebsiella pneumoniae infections in hospital patients reveals insights into an opportunistic pathogen. Nat. Commun. 2022, 13, 3017. [Google Scholar] [CrossRef]
  3. Wyres, K.L.; Lam, M.M.C.; Holt, K.E. Population genomics of Klebsiella pneumoniae. Nat. Rev. Microbiol. 2020, 18, 344–359. [Google Scholar] [CrossRef] [PubMed]
  4. Chen, D.; Zhang, Y.; Wu, J.; Li, J.; Chen, H.; Zhang, X.; Hu, X.; Chen, F.; Yu, R. Analysis of hypervirulent Klebsiella pneumoniae and classic Klebsiella pneumoniae infections in a Chinese hospital. J. Appl. Microbiol. 2022, 132, 3883–3890. [Google Scholar] [CrossRef]
  5. Navon-Venezia, S.; Kondratyeva, K.; Carattoli, A. Klebsiella pneumoniae: A major worldwide source and shuttle for antibiotic resistance. FEMS. Microbiol. Rev. 2017, 41, 252–275. [Google Scholar] [CrossRef] [PubMed]
  6. Arcari, G.; Carattoli, A. Global spread and evolutionary convergence of multidrug-resistant and hypervirulent Klebsiella pneumoniae high-risk clones. Pathog. Glob. Health 2022, 1–14. [Google Scholar] [CrossRef]
  7. Lee, C.R.; Lee, J.H.; Park, K.S.; Jeon, J.H.; Kim, Y.B.; Cha, C.J.; Jeong, B.C.; Lee, S.H. Antimicrobial Resistance of Hypervirulent Klebsiella pneumoniae: Epidemiology, Hypervirulence-Associated Determinants, and Resistance Mechanisms. Front. Cell. Infect. Microbiol. 2017, 7, 483. [Google Scholar] [CrossRef]
  8. Choby, J.E.; Howard-Anderson, J.; Weiss, D.S. Hypervirulent Klebsiella pneumoniae—Clinical and molecular perspectives. J. Intern. Med. 2020, 287, 283–300. [Google Scholar] [CrossRef]
  9. Russo, T.A.; Olson, R.; Fang, C.T.; Stoesser, N.; Miller, M.; MacDonald, U.; Hutson, A.; Barker, J.H.; La Hoz, R.M.; Johnson, J.R. Identification of Biomarkers for Differentiation of Hypervirulent Klebsiella pneumoniae from Classical K. pneumoniae. J. Clin. Microbiol. 2018, 56, e00776-18. [Google Scholar] [CrossRef]
  10. Tang, M.; Kong, X.; Hao, J.; Liu, J. Epidemiological Characteristics and Formation Mechanisms of Multidrug-Resistant Hypervirulent Klebsiella pneumoniae. Front. Microbiol. 2020, 11, 581543. [Google Scholar] [CrossRef]
  11. Wyres, K.L.; Nguyen, T.N.T.; Lam, M.M.C.; Judd, L.M.; van Vinh Chau, N.; Dance, D.A.B.; Ip, M.; Karkey, A.; Ling, C.L.; Miliya, T.; et al. Genomic surveillance for hypervirulence and multi-drug resistance in invasive Klebsiella pneumoniae from South and Southeast Asia. Genome. Med. 2020, 12, 11. [Google Scholar] [CrossRef]
  12. Bialek-Davenet, S.; Criscuolo, A.; Ailloud, F.; Passet, V.; Jones, L.; Delannoy-Vieillard, A.S.; Garin, B.; Le Hello, S.; Arlet, G.; Nicolas-Chanoine, M.H.; et al. Genomic definition of hypervirulent and multidrug-resistant Klebsiella pneumoniae clonal groups. Emerg. Infect. Dis. 2014, 20, 1812–1820. [Google Scholar] [CrossRef]
  13. Su, S.C.; Siu, L.K.; Ma, L.; Yeh, K.M.; Fung, C.P.; Lin, J.C.; Chang, F.Y. Community-acquired liver abscess caused by serotype K1 Klebsiella pneumoniae with CTX-M-15-type extended-spectrum beta-lactamase. Antimicrob. Agents Chemother. 2008, 52, 804–805. [Google Scholar] [CrossRef]
  14. Cejas, D.; Fernandez Canigia, L.; Rincon Cruz, G.; Elena, A.X.; Maldonado, I.; Gutkind, G.O.; Radice, M.A. First isolate of KPC-2-producing Klebsiella pneumonaie sequence type 23 from the Americas. J. Clin. Microbiol. 2014, 52, 3483–3485. [Google Scholar] [CrossRef] [PubMed]
  15. Surgers, L.; Boyd, A.; Girard, P.M.; Arlet, G.; Decre, D. ESBL-Producing Strain of Hypervirulent Klebsiella pneumoniae K2, France. Emerg. Infect. Dis. 2016, 22, 1687–1688. [Google Scholar] [CrossRef] [PubMed]
  16. Zhang, R.; Lin, D.; Chan, E.W.; Gu, D.; Chen, G.X.; Chen, S. Emergence of Carbapenem-Resistant Serotype K1 Hypervirulent Klebsiella pneumoniae Strains in China. Antimicrob. Agents Chemother. 2016, 60, 709–711. [Google Scholar] [CrossRef] [PubMed]
  17. Shaidullina, E.; Shelenkov, A.; Yanushevich, Y.; Mikhaylova, Y.; Shagin, D.; Alexandrova, I.; Ershova, O.; Akimkin, V.; Kozlov, R.; Edelstein, M. Antimicrobial Resistance and Genomic Characterization of OXA-48- and CTX-M-15-Co-Producing Hypervirulent Klebsiella pneumoniae ST23 Recovered from Nosocomial Outbreak. Antibiotics 2020, 9, 862. [Google Scholar] [CrossRef] [PubMed]
  18. Khaertynov, K.S.; Anokhin, V.A.; Davidyuk, Y.N.; Nicolaeva, I.V.; Khalioullina, S.V.; Semyenova, D.R.; Alatyrev, E.Y.; Skvortsova, N.N.; Abrahamyan, L.G. Case of Meningitis in a Neonate Caused by an Extended-Spectrum-Beta-Lactamase-Producing Strain of Hypervirulent Klebsiella pneumoniae. Front. Microbiol. 2017, 8, 1576. [Google Scholar] [CrossRef]
  19. Gu, D.; Dong, N.; Zheng, Z.; Lin, D.; Huang, M.; Wang, L.; Chan, E.W.; Shu, L.; Yu, J.; Zhang, R.; et al. A fatal outbreak of ST11 carbapenem-resistant hypervirulent Klebsiella pneumoniae in a Chinese hospital: A molecular epidemiological study. Lancet Infect. Dis. 2018, 18, 37–46. [Google Scholar] [CrossRef]
  20. Hennequin, C.; Robin, F. Correlation between antimicrobial resistance and virulence in Klebsiella pneumoniae. Eur. J. Clin. Microbiol. Infect. Dis. 2016, 35, 333–341. [Google Scholar] [CrossRef]
  21. Wyres, K.L.; Wick, R.R.; Judd, L.M.; Froumine, R.; Tokolyi, A.; Gorrie, C.L.; Lam, M.M.C.; Duchene, S.; Jenney, A.; Holt, K.E. Distinct evolutionary dynamics of horizontal gene transfer in drug resistant and virulent clones of Klebsiella pneumoniae. PLoS Genet. 2019, 15, e1008114. [Google Scholar] [CrossRef] [PubMed]
  22. Wyres, K.L.; Holt, K.E. Klebsiella pneumoniae as a key trafficker of drug resistance genes from environmental to clinically important bacteria. Curr. Opin. Microbiol. 2018, 45, 131–139. [Google Scholar] [CrossRef]
  23. Partridge, S.R.; Kwong, S.M.; Firth, N.; Jensen, S.O. Mobile Genetic Elements Associated with Antimicrobial Resistance. Clin. Microbiol. Rev. 2018, 31, e00088-17. [Google Scholar] [CrossRef] [PubMed]
  24. Peter, S.; Bosio, M.; Gross, C.; Bezdan, D.; Gutierrez, J.; Oberhettinger, P.; Liese, J.; Vogel, W.; Dorfel, D.; Berger, L.; et al. Tracking of Antibiotic Resistance Transfer and Rapid Plasmid Evolution in a Hospital Setting by Nanopore Sequencing. mSphere 2020, 5, e00525-20. [Google Scholar] [CrossRef] [PubMed]
  25. Shelenkov, A.; Petrova, L.; Mironova, A.; Zamyatin, M.; Akimkin, V.; Mikhaylova, Y. Long-Read Whole Genome Sequencing Elucidates the Mechanisms of Amikacin Resistance in Multidrug-Resistant Klebsiella pneumoniae Isolates Obtained from COVID-19 Patients. Antibiotics 2022, 11, 1364. [Google Scholar] [CrossRef] [PubMed]
  26. Kuzmenkov, A.Y.; Trushin, I.V.; Vinogradova, A.G.; Avramenko, A.A.; Sukhorukova, M.V.; Malhotra-Kumar, S.; Dekhnich, A.V.; Edelstein, M.V.; Kozlov, R.S. AMRmap: An Interactive Web Platform for Analysis of Antimicrobial Resistance Surveillance Data in Russia. Front. Microbiol. 2021, 12, 620002. [Google Scholar] [CrossRef]
  27. Phillippy, A.M.; Wick, R.R.; Judd, L.M.; Gorrie, C.L.; Holt, K.E. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput. Biol. 2017, 13, e1005697. [Google Scholar] [CrossRef]
  28. Kolmogorov, M.; Yuan, J.; Lin, Y.; Pevzner, P.A. Assembly of long, error-prone reads using repeat graphs. Nat. Biotechnol. 2019, 37, 540–546. [Google Scholar] [CrossRef]
  29. Wick, R.R.; Holt, K.E. Polypolish: Short-read polishing of long-read bacterial genome assemblies. PLoS Comput. Biol. 2022, 18, e1009802. [Google Scholar] [CrossRef] [PubMed]
  30. Shelenkov, A.; Mikhaylova, Y.; Yanushevich, Y.; Samoilov, A.; Petrova, L.; Fomina, V.; Gusarov, V.; Zamyatin, M.; Shagin, D.; Akimkin, V. Molecular Typing, Characterization of Antimicrobial Resistance, Virulence Profiling and Analysis of Whole-Genome Sequence of Clinical Klebsiella pneumoniae Isolates. Antibiotics 2020, 9, 261. [Google Scholar] [CrossRef] [PubMed]
  31. Shelenkov, A.; Petrova, L.; Zamyatin, M.; Mikhaylova, Y.; Akimkin, V. Diversity of International High-Risk Clones of Acinetobacter baumannii Revealed in a Russian Multidisciplinary Medical Center during 2017–2019. Antibiotics 2021, 10, 1009. [Google Scholar] [CrossRef] [PubMed]
  32. Robertson, J.; Bessonov, K.; Schonfeld, J.; Nash, J.H.E. Universal whole-sequence-based plasmid typing and its utility to prediction of host range and epidemiological surveillance. Microb. Genom. 2020, 6. [Google Scholar] [CrossRef] [PubMed]
  33. Couvin, D.; Bernheim, A.; Toffano-Nioche, C.; Touchon, M.; Michalik, J.; Neron, B.; Rocha, E.P.C.; Vergnaud, G.; Gautheret, D.; Pourcel, C. CRISPRCasFinder, an update of CRISRFinder, includes a portable version, enhanced performance and integrates search for Cas proteins. Nucleic Acids. Res. 2018, 46, W246–W251. [Google Scholar] [CrossRef]
  34. Biswas, A.; Gagnon, J.N.; Brouns, S.J.; Fineran, P.C.; Brown, C.M. CRISPRTarget: Bioinformatic prediction and analysis of crRNA targets. RNA Biol. 2013, 10, 817–827. [Google Scholar] [CrossRef] [PubMed]
  35. Wyres, K.L.; Wick, R.R.; Gorrie, C.; Jenney, A.; Follador, R.; Thomson, N.R.; Holt, K.E. Identification of Klebsiella capsule synthesis loci from whole genome data. Microb. Genom. 2016, 2, e000102. [Google Scholar] [CrossRef]
  36. Lam, M.M.C.; Wick, R.R.; Watts, S.C.; Cerdeira, L.T.; Wyres, K.L.; Holt, K.E. A genomic surveillance framework and genotyping tool for Klebsiella pneumoniae and its related species complex. Nat. Commun. 2021, 12, 4188. [Google Scholar] [CrossRef]
  37. Feijao, P.; Yao, H.T.; Fornika, D.; Gardy, J.; Hsiao, W.; Chauve, C.; Chindelevitch, L. MentaLiST—A fast MLST caller for large MLST schemes. Microb. Genom. 2018, 4. [Google Scholar] [CrossRef]
  38. Schurch, A.C.; Arredondo-Alonso, S.; Willems, R.J.L.; Goering, R.V. Whole genome sequencing options for bacterial strain typing and epidemiologic analysis based on single nucleotide polymorphism versus gene-by-gene-based approaches. Clin. Microbiol. Infect. 2018, 24, 350–354. [Google Scholar] [CrossRef]
  39. Loftie-Eaton, W.; Rawlings, D.E. Diversity, biology and evolution of IncQ-family plasmids. Plasmid 2012, 67, 15–34. [Google Scholar] [CrossRef]
  40. Smillie, C.; Garcillan-Barcia, M.P.; Francia, M.V.; Rocha, E.P.; de la Cruz, F. Mobility of plasmids. Microbiol. Mol. Biol. Rev. 2010, 74, 434–452. [Google Scholar] [CrossRef]
  41. Eger, E.; Schwabe, M.; Schulig, L.; Hubner, N.O.; Bohnert, J.A.; Bornscheuer, U.T.; Heiden, S.E.; Muller, J.U.; Adnan, F.; Becker, K.; et al. Extensively Drug-Resistant Klebsiella pneumoniae Counteracts Fitness and Virulence Costs That Accompanied Ceftazidime-Avibactam Resistance Acquisition. Microbiol. Spectr. 2022, 10, e0014822. [Google Scholar] [CrossRef] [PubMed]
  42. Kochan, T.J.; Nozick, S.H.; Medernach, R.L.; Cheung, B.H.; Gatesy, S.W.M.; Lebrun-Corbin, M.; Mitra, S.D.; Khalatyan, N.; Krapp, F.; Qi, C.; et al. Genomic surveillance for multidrug-resistant or hypervirulent Klebsiella pneumoniae among United States bloodstream isolates. BMC. Infect. Dis. 2022, 22, 603. [Google Scholar] [CrossRef] [PubMed]
  43. Starkova, P.; Lazareva, I.; Avdeeva, A.; Sulian, O.; Likholetova, D.; Ageevets, V.; Lebedeva, M.; Gostev, V.; Sopova, J.; Sidorenko, S. Emergence of Hybrid Resistance and Virulence Plasmids Harboring New Delhi Metallo-beta-Lactamase in Klebsiella pneumoniae in Russia. Antibiotics 2021, 10, 691. [Google Scholar] [CrossRef] [PubMed]
  44. Liu, S.; Ding, Y.; Xu, Y.; Li, Z.; Zeng, Z.; Liu, J. An outbreak of extensively drug-resistant and hypervirulent Klebsiella pneumoniae in an intensive care unit of a teaching hospital in Southwest China. Front. Cell. Infect. Microbiol. 2022, 12, 979219. [Google Scholar] [CrossRef]
  45. Davoudabadi, S.; Goudarzi, H.; Goudarzi, M.; Ardebili, A.; Faghihloo, E.; Sharahi, J.Y.; Hashemi, A. Detection of extensively drug-resistant and hypervirulent Klebsiella pneumoniae ST15, ST147, ST377 and ST442 in Iran. Acta Microbiol. Immunol. Hung. 2021, 69, 77–86. [Google Scholar] [CrossRef]
  46. Liu, C.; Du, P.; Xiao, N.; Ji, F.; Russo, T.A.; Guo, J. Hypervirulent Klebsiella pneumoniae is emerging as an increasingly prevalent K. pneumoniae pathotype responsible for nosocomial and healthcare-associated infections in Beijing, China. Virulence 2020, 11, 1215–1224. [Google Scholar] [CrossRef]
  47. Bolourchi, N.; Naz, A.; Sohrabi, M.; Badmasti, F. Comparative in silico characterization of Klebsiella pneumoniae hypervirulent plasmids and their antimicrobial resistance genes. Ann. Clin. Microbiol. Antimicrob. 2022, 21, 23. [Google Scholar] [CrossRef]
  48. Li, P.; Luo, W.; Xiang, T.X.; Jiang, Y.; Liu, P.; Wei, D.D.; Fan, L.; Huang, S.; Liao, W.; Liu, Y.; et al. Horizontal gene transfer via OMVs co-carrying virulence and antimicrobial-resistant genes is a novel way for the dissemination of carbapenem-resistant hypervirulent Klebsiella pneumoniae. Front. Microbiol. 2022, 13, 945972. [Google Scholar] [CrossRef]
  49. Holmes, A.H.; Moore, L.S.; Sundsfjord, A.; Steinbakk, M.; Regmi, S.; Karkey, A.; Guerin, P.J.; Piddock, L.J. Understanding the mechanisms and drivers of antimicrobial resistance. Lancet 2016, 387, 176–187. [Google Scholar] [CrossRef]
  50. Lan, P.; Jiang, Y.; Zhou, J.; Yu, Y. A global perspective on the convergence of hypervirulence and carbapenem resistance in Klebsiella pneumoniae. J. Glob. Antimicrob. Resist. 2021, 25, 26–34. [Google Scholar] [CrossRef]
  51. Dong, N.; Yang, X.; Chan, E.W.; Zhang, R.; Chen, S. Klebsiella species: Taxonomy, hypervirulence and multidrug resistance. eBioMedicine 2022, 79, 103998. [Google Scholar] [CrossRef] [PubMed]
  52. Liu, C.; Du, P.; Yang, P.; Lu, M.; Shen, N. Fusion plasmid enhanced the endemic extensively drug resistant Klebsiella pneumoniae clone ST147 harbored bla(OXA-48) to acquire the hypervirulence and cause fatal infection. Ann. Clin. Microbiol. Antimicrob. 2023, 22, 11. [Google Scholar] [CrossRef] [PubMed]
  53. Shankar, C.; Vasudevan, K.; Jacob, J.J.; Baker, S.; Isaac, B.J.; Neeravi, A.R.; Sethuvel, D.P.M.; George, B.; Veeraraghavan, B. Hybrid Plasmids Encoding Antimicrobial Resistance and Virulence Traits Among Hypervirulent Klebsiella pneumoniae ST2096 in India. Front. Cell. Infect. Microbiol. 2022, 12, 875116. [Google Scholar] [CrossRef] [PubMed]
  54. Ozad Duzgun, A. From Turkey: First Report of KPC-3- and CTX-M-27-Producing Multidrug-Resistant Klebsiella pneumoniae ST147 Clone Carrying OmpK36 and Ompk37 Porin Mutations. Microb. Drug. Resist. 2021, 27, 1265–1270. [Google Scholar] [CrossRef]
  55. Ngbede, E.O.; Adekanmbi, F.; Poudel, A.; Kalalah, A.; Kelly, P.; Yang, Y.; Adamu, A.M.; Daniel, S.T.; Adikwu, A.A.; Akwuobu, C.A.; et al. Concurrent Resistance to Carbapenem and Colistin Among Enterobacteriaceae Recovered from Human and Animal Sources in Nigeria is Associated with Multiple Genetic Mechanisms. Front. Microbiol. 2021, 12, 740348. [Google Scholar] [CrossRef]
  56. Sakamoto, N.; Akeda, Y.; Sugawara, Y.; Takeuchi, D.; Motooka, D.; Yamamoto, N.; Laolerd, W.; Santanirand, P.; Hamada, S. Genomic Characterization of Carbapenemase-Producing Klebsiella pneumoniae with Chromosomally Carried bla(NDM-1). Antimicrob. Agents Chemother. 2018, 62, e01520-18. [Google Scholar] [CrossRef]
  57. Wang, X.; Zhao, J.; Ji, F.; Chang, H.; Qin, J.; Zhang, C.; Hu, G.; Zhu, J.; Yang, J.; Jia, Z.; et al. Multiple-Replicon Resistance Plasmids of Klebsiella Mediate Extensive Dissemination of Antimicrobial Genes. Front. Microbiol. 2021, 12, 754931. [Google Scholar] [CrossRef]
  58. Trasanidou, D.; Geros, A.S.; Mohanraju, P.; Nieuwenweg, A.C.; Nobrega, F.L.; Staals, R.H.J. Keeping crispr in check: Diverse mechanisms of phage-encoded anti-crisprs. FEMS Microbiol. Lett. 2019, 366, fnz098. [Google Scholar] [CrossRef]
  59. Shen, J.; Lv, L.; Wang, X.; Xiu, Z.; Chen, G. Comparative analysis of CRISPR-Cas systems in Klebsiella genomes. J. Basic Microbiol. 2017, 57, 325–336. [Google Scholar] [CrossRef]
  60. Zhou, Y.; Tang, Y.; Fu, P.; Tian, D.; Yu, L.; Huang, Y.; Li, G.; Li, M.; Wang, Y.; Yang, Z.; et al. The type I-E CRISPR-Cas system influences the acquisition of bla(KPC)-IncF plasmid in Klebsiella pneumonia. Emerg. Microbes Infect. 2020, 9, 1011–1022. [Google Scholar] [CrossRef]
  61. Mackow, N.A.; Shen, J.; Adnan, M.; Khan, A.S.; Fries, B.C.; Diago-Navarro, E. CRISPR-Cas influences the acquisition of antibiotic resistance in Klebsiella pneumoniae. PLoS ONE 2019, 14, e0225131. [Google Scholar] [CrossRef] [PubMed]
  62. Hu, Y.; Jiang, J.; Wang, D.; Guo, Q.; Wang, M. Coexistence of blaKPC-IncFII plasmids and type I-E* CRISPR-Cas systems in ST15 Klebsiella pneumoniae. Front. Microbiol. 2023, 14, 1125531. [Google Scholar] [CrossRef] [PubMed]
  63. Kamruzzaman, M.; Iredell, J.R. CRISPR-Cas System in Antibiotic Resistance Plasmids in Klebsiella pneumoniae. Front. Microbiol. 2019, 10, 2934. [Google Scholar] [CrossRef] [PubMed]
  64. Zhou, Y.; Yang, Y.; Li, X.; Tian, D.; Ai, W.; Wang, W.; Wang, B.; Kreiswirth, B.N.; Yu, F.; Chen, L.; et al. Exploiting a conjugative endogenous CRISPR-Cas3 system to tackle multidrug-resistant Klebsiella pneumoniae. eBioMedicine 2023, 88, 104445. [Google Scholar] [CrossRef]
  65. Guo, X.; Chen, R.; Wang, Q.; Li, C.; Ge, H.; Qiao, J.; Li, Y. Global prevalence, characteristics, and future prospects of IncX3 plasmids: A review. Front. Microbiol. 2022, 13, 979558. [Google Scholar] [CrossRef]
  66. Viehweger, A.; Blumenscheit, C.; Lippmann, N.; Wyres, K.L.; Brandt, C.; Hans, J.B.; Holzer, M.; Irber, L.; Gatermann, S.; Lubbert, C.; et al. Context-aware genomic surveillance reveals hidden transmission of a carbapenemase-producing Klebsiella pneumoniae. Microb. Genom. 2021, 7. [Google Scholar] [CrossRef]
  67. Hodcroft, E.B.; De Maio, N.; Lanfear, R.; MacCannell, D.R.; Minh, B.Q.; Schmidt, H.A.; Stamatakis, A.; Goldman, N.; Dessimoz, C. Want to track pandemic variants faster? Fix the bioinformatics bottleneck. Nature 2021, 591, 30–33. [Google Scholar] [CrossRef]
  68. Karampatakis, T.; Tsergouli, K.; Behzadi, P. Carbapenem-Resistant Klebsiella pneumoniae: Virulence Factors, Molecular Epidemiology and Latest Updates in Treatment Options. Antibiotics 2023, 12, 234. [Google Scholar] [CrossRef]
  69. Shelenkov, A.A.; Slavokhotova, A.A.; Odintsova, T.I. Cysmotif Searcher Pipeline for Antimicrobial Peptide Identification in Plant Transcriptomes. Biochemistry 2018, 83, 1424–1432. [Google Scholar] [CrossRef]
  70. de Souza, C.M.; da Silva, A.P.; Junior, N.G.O.; Martinez, O.F.; Franco, O.L. Peptides as a therapeutic strategy against Klebsiella pneumoniae. Trends. Pharm. Sci. 2022, 43, 335–348. [Google Scholar] [CrossRef]
  71. Gorodnichev, R.B.; Volozhantsev, N.V.; Krasilnikova, V.M.; Bodoev, I.N.; Kornienko, M.A.; Kuptsov, N.S.; Popova, A.V.; Makarenko, G.I.; Manolov, A.I.; Slukin, P.V.; et al. Novel Klebsiella pneumoniae K23-Specific Bacteriophages from Different Families: Similarity of Depolymerases and Their Therapeutic Potential. Front. Microbiol. 2021, 12, 669618. [Google Scholar] [CrossRef] [PubMed]
  72. De Carvalho, F.R.T.; Telles, J.P.; Tuon, F.F.B.; Rabello Filho, R.; Caruso, P.; Correa, T.D. Antimicrobial Stewardship Programs: A Review of Strategies to Avoid Polymyxins and Carbapenems Misuse in Low Middle-Income Countries. Antibiotics 2022, 11, 378. [Google Scholar] [CrossRef] [PubMed]
Microorganisms 11 01314 g001 550
Figure 1. Phenotypic and genomic antibiotic resistance profiles for clinical K. pneumoniae isolates studied. Cft/avb—ceftazidime/avibactam; pip/taz—piperacillin/tazobactam; tmp/sul—trimethoprim/sulfamethoxazole.
Figure 1. Phenotypic and genomic antibiotic resistance profiles for clinical K. pneumoniae isolates studied. Cft/avb—ceftazidime/avibactam; pip/taz—piperacillin/tazobactam; tmp/sul—trimethoprim/sulfamethoxazole.
Microorganisms 11 01314 g001
Microorganisms 11 01314 g002 550
Figure 2. Important virulence factors revealed in the K. pneumoniae isolates studied.
Figure 2. Important virulence factors revealed in the K. pneumoniae isolates studied.
Microorganisms 11 01314 g002
Microorganisms 11 01314 g003 550
Figure 3. Circular diagram of IncHI1B/IncFIB plasmids from CriePt isolates and reference plasmid phvKpST147_NDM-1_2566. Antimicrobial resistance genes are indicated in red and virulence factors in gray.
Figure 3. Circular diagram of IncHI1B/IncFIB plasmids from CriePt isolates and reference plasmid phvKpST147_NDM-1_2566. Antimicrobial resistance genes are indicated in red and virulence factors in gray.
Microorganisms 11 01314 g003
Microorganisms 11 01314 g004 550
Figure 4. The scheme of the CRISPR-Cas system for K. pneumoniae isolate CriePt491.
Figure 4. The scheme of the CRISPR-Cas system for K. pneumoniae isolate CriePt491.
Microorganisms 11 01314 g004
Microorganisms 11 01314 g005 550
Figure 5. Minimum spanning cgMLST-based tree for ST512 K. pneumoniae isolates revealed by us and reference isolates from Genbank. Values of the pairwise allele differences between the isolates are given in red.
Figure 5. Minimum spanning cgMLST-based tree for ST512 K. pneumoniae isolates revealed by us and reference isolates from Genbank. Values of the pairwise allele differences between the isolates are given in red.
Microorganisms 11 01314 g005
Table
Table 1. Metadata for clinical K. pneumoniae isolates obtained from patients of ICU and pulmonology departments of a single tertiary care hospital.
Table 1. Metadata for clinical K. pneumoniae isolates obtained from patients of ICU and pulmonology departments of a single tertiary care hospital.
Isolate IdPatient Age (Years)Patient GenderIsolation DateDepartmentLocusDiagnosis (Main)
CriePt49175female12 December 2021ICU 1BAL 2Thrombosis of the superior mesenteric artery with necrosis of the colon
CriePt49282female12 December 2021ICUBALAcute thrombosis of the femoral artery
CriePt49464male8 December 2021pulmonologysputumInterstitial pneumonitis
CriePt49558male9 December 2021ICUBALIntracerebral parenchymal hemorrhage in the suprasellar region of the brain
1 ICU—intensive care unit; 2 BAL—bronchoalveolar lavage.
Table
Table 2. Plasmid replicon composition for the clinical K. pneumoniae isolates studied.
Table 2. Plasmid replicon composition for the clinical K. pneumoniae isolates studied.
Sample Id/PlasmidsnumCol(pHAD28)ColRNAIIncFIAIncFIBIncHI1BIncRIncX3Unkn1 1Unkn2
CriePt4915+--+++--+
CriePt4926++-++-++-
CriePt4947+++++-++-
CriePt4956++-++-++-
1 Unkn1 and Unkn2 stand for two different types of non-typeable plasmids.
Table
Table 3. Extended description of plasmids for the clinical K. pneumoniae isolates studied.
Table 3. Extended description of plasmids for the clinical K. pneumoniae isolates studied.
IDSize, bpGC ContentReplicon Type (s)Relaxase Typempf_TypePredicted MobilityAMR/Virulence (VIR) Genes
CriePt491
p2179,7070.448IncHI1B/IncFIB-MPF_Fnon-mobilizable-
p3109,6490.493IncFIB--non-mobilizable-
p453,0090.534---non-mobilizable-
p546,9840.534IncR--mobilizableAMR
p649150.431Col(pHAD28)--non-mobilizable-
CriePt492
p2250,6300.461IncHI1B/IncFIBMOBHMPF_Fconjugative-
p3184,7250.533IncFIBMOBF-mobilizableAMR
p453,2920.496IncX3MOBPMPF_TconjugativeAMR
p510,6890.475Col(pHAD28)--non-mobilizableAMR
p697300.532ColRNAIMOBC-mobilizable-
p729630.650-MOBV-mobilizable-
CriePt494
p2308,2080.472IncHI1B/IncFIBMOBHMPF_FconjugativeAMR, VIR
p393,8080.521IncFIB--non-mobilizable-
p489,7530.510IncFIAMOBFMPF_FconjugativeAMR
p553,2920.496IncX3MOBPMPF_TconjugativeAMR
p610,6890.475Col(pHAD28)--non-mobilizableAMR
p797370.532ColRNAIMOBC-mobilizable-
p829630.650-MOBV-mobilizable-
CriePt495
p2251,8300.461IncHI1B/IncFIBMOBHMPF_Fconjugative-
p3190,4760.533IncFIBMOBF-mobilizableAMR
p453,2920.496IncX3MOBPMPF_TconjugativeAMR
p510,6890.475Col(pHAD28)--non-mobilizableAMR
p697300.532ColRNAIMOBC-mobilizable-
p729630.650-MOBV-mobilizable-
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shelenkov, A.; Mikhaylova, Y.; Voskanyan, S.; Egorova, A.; Akimkin, V. Whole-Genome Sequencing Revealed the Fusion Plasmids Capable of Transmission and Acquisition of Both Antimicrobial Resistance and Hypervirulence Determinants in Multidrug-Resistant Klebsiella pneumoniae Isolates. Microorganisms 2023, 11, 1314. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms11051314

AMA Style

Shelenkov A, Mikhaylova Y, Voskanyan S, Egorova A, Akimkin V. Whole-Genome Sequencing Revealed the Fusion Plasmids Capable of Transmission and Acquisition of Both Antimicrobial Resistance and Hypervirulence Determinants in Multidrug-Resistant Klebsiella pneumoniae Isolates. Microorganisms. 2023; 11(5):1314. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms11051314

Chicago/Turabian Style

Shelenkov, Andrey, Yulia Mikhaylova, Shushanik Voskanyan, Anna Egorova, and Vasiliy Akimkin. 2023. "Whole-Genome Sequencing Revealed the Fusion Plasmids Capable of Transmission and Acquisition of Both Antimicrobial Resistance and Hypervirulence Determinants in Multidrug-Resistant Klebsiella pneumoniae Isolates" Microorganisms 11, no. 5: 1314. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms11051314

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop