Phylogenomics of Globally Spread Clonal Groups 14 and 15 of Klebsiella pneumoniae

Klebsiella pneumoniae sequence type 14 (ST14) and ST15 caused outbreaks of CTX-M-15 and/or carbapenemase producers worldwide, but their phylogeny and global dynamics remain unclear. We clarified the evolution of K. pneumoniae clonal group 14 (CG14) and CG15 by analyzing the capsular locus (KL), resistome, virulome, and plasmidome of public genomes (n = 481) and de novo sequences (n = 9) representing main sublineages circulating in Portugal. CG14 and CG15 evolved independently within 6 main subclades defined according to the KL and the accessory genome. The CG14 (n = 65) clade was structured in two large monophyletic subclades, CG14-I (KL2, 86%) and CG14-II (KL16, 14%), whose emergences were dated to 1932 and 1911, respectively. Genes encoding extended-spectrum β-lactamase (ESBL), AmpC, and/or carbapenemases were mostly observed in CG14-I (71% versus 22%). CG15 clade (n = 170) was segregated into subclades CG15-IA (KL19/KL106, 9%), CG15-IB (variable KL types, 6%), CG15-IIA (KL24, 43%) and CG15-IIB (KL112, 37%). Most CG15 genomes carried specific GyrA and ParC mutations and emerged from a common ancestor in 1989. CTX-M-15 was especially prevalent in CG15 (68% CG15 versus 38% CG14) and in CG15-IIB (92%). Plasmidome analysis revealed 27 predominant plasmid groups (PG), including particularly pervasive and recombinant F-type (n = 10), Col (n = 10), and new plasmid types. While blaCTX-M-15 was acquired multiple times by a high diversity of F-type mosaic plasmids, other antibiotic resistance genes (ARGs) were dispersed by IncL (blaOXA-48) or IncC (blaCMY/TEM-24) plasmids. We first demonstrate an independent evolutionary trajectory for CG15 and CG14 and how the acquisition of specific KL, quinolone-resistance determining region (QRDR) mutations (CG15), and ARGs in highly recombinant plasmids could have shaped the expansion and diversification of particular subclades (CG14-I and CG15-IIA/IIB). IMPORTANCE Klebsiella pneumoniae represents a major threat in the burden of antibiotic resistance (ABR). Available studies to explain the origin, the diversity, and the evolution of certain ABR K. pneumoniae populations have mainly been focused on a few clonal groups (CGs) using phylogenetic analysis of the core genome, the accessory genome being overlooked. Here, we provide unique insights into the phylogenetic evolution of CG14 and CG15, two poorly characterized CGs which have contributed to the global dissemination of genes responsible for resistance to first-line antibiotics such as β-lactams. Our results point out an independent evolution of these two CGs and highlight the existence of different subclades structured by the capsular type and the accessory genome. Moreover, the contribution of a turbulent flux of plasmids (especially multireplicon F type and Col) and adaptive traits (antibiotic resistance and metal tolerance genes) to the pangenome reflect the exposure and adaptation of K. pneumoniae under different selective pressures.

IS26-mediated plasmid reshuffling results in convergence of toxin-antitoxin systems but loss of resistance genes in XDR Klebsiella pneumoniae from a chronic infection

Carbapenem-resistant Enterobacterales pose an urgent threat to human health worldwide. Klebsiella pneumoniae sequence type (ST) 14, initially identified in the Middle East and South-Asia and co-harbouring the carbapenemase genes bla _OXA-232 and bla _NDM-1, is now emerging globally. One such strain was detected in the USA in 2013 from a patient initially treated in India that also carried armA, a 16S rRNA methyltransferase that confers resistance to all clinically relevant aminoglycosides. Genetic and phenotypic changes were observed in 14 serial isolates collected from this chronically infected patient. The index isolate carried five plasmids, including an IncFIB-IncHI1B (harbouring armA and bla _NDM-1), an IncFIA (bla _CTX-M-15) and a ColE-like (bla _OXA-232), and was extensively resistant to antibiotics. Four years later, a subsequent isolate had accumulated 34 variants, including a loss-of-function mutation in romA, resulting in tigecycline non-susceptibility. Importantly, this isolate now only carried two plasmids, including a large mosaic molecule made of fragments, all harbouring distinct toxin-antitoxin systems, from three of the canonical plasmids. Of the original acquired antibiotic resistance genes, this isolate only retained bla _CTX-M-15, and as a result susceptibility to the carbapenems and amikacin was restored. Long-read sequencing of a subset of five representative isolates, collected between 2013 and 2017, allowed for the elucidation of the complex plasmid patterns and revealed the role of IS26-mediated plasmid reshuffling in the evolution of this clone. Such investigations of the mechanisms underlying plasmid stability, together with global and local surveillance programmes, are key to a better understanding of plasmid host range and dissemination.