Genomic evolution of the globally disseminated multidrug-resistantKlebsiella pneumoniaeclonal group 147

Abstract

Background: The rapid emergence of multidrug-resistant Klebsiella pneumoniae (Kp) is largely driven by the spread of specific clonal groups (CG). Of these, CG147 includes 7-gene MLST sequence types ST147, ST273 and ST392. CG147 has caused nosocomial outbreaks across the world, but its global population dynamics remain unknown. Here, we report a pandrug-resistant ST147 clinical isolate from India (strain DJ) and define the evolution and global emergence of CG147. Methods: Antimicrobial susceptibility testing (EUCAST guidelines) and genome sequencing (Illumina and Oxford Nanopore technologies, Unicycler assembly) were performed on strain DJ. Additionally, we collated 217 publicly available CG147 genomes (NCBI, May 2019). CG147 evolution was inferred within a temporal phylogenetic framework (BEAST) based on a recombination-free sequence alignment (Roary/Gubbins). Comparative genomic analyses focused on resistance and virulence genes and other genetic elements (BIGSdb, Kleborate, PlasmidFinder, PHASTER, ICEFinder and CRISPRCasFinder). Results: Strain DJ had a pandrug resistance phenotype. Its genome comprised 7 plasmids and 1 linear phage-plasmid. Four carbapenemase genes were detected: blaNDM-5 and 2 copies of blaOXA-181 in the chromosome, and a second copy of blaNDM-5 on an 84 kb IncFII plasmid. CG147 genomes carried a mean of 13 acquired resistance genes or mutations; 63% carried a carbapenemase gene and 83% harbored blaCTX-M. All CG147 genomes presented GyrA and ParC mutations and a common subtype IV-E CRISPR-Cas system. ST392 and ST273 emerged in 2005 and 1995, respectively. ST147, the most represented phylogenetic branch, was itself divided into two main clades with distinct capsular loci: KL64 (74%, DJ included, emerged in 1994 and disseminated worldwide, with carbapenemases varying among world regions) and KL10 (20%, 2002, predominantly found in Asian countries, associated with carbapenemases NDM and OXA-48-like). Further, subclades within ST147-KL64 differed in the yersiniabactin locus, OmpK35/K36 mutations, plasmid replicons and prophages. The absence of IncF plasmids in some subclades was associated with a possible activity of a CRISPR-Cas system. Conclusions: K. pneumoniae clonal group CG147 comprises pandrug- or extensively-resistant isolates and carries multiple and diverse resistance genes and mobile genetic elements, including chromosomal blaNDM-5. Its emergence is driven by the spread of several phylogenetic clades marked by their own genomic features and specific temporo-spatial dynamics. These findings highlight the need for precision surveillance strategies to limit the spread of particularly concerning CG147 subsets.

Anatomy of an extensively drug-resistant Klebsiella pneumoniae outbreak in Tuscany, Italy

“…To date, internal nodes of interest on the phylogeny and Bayesian phylogenetic inference were performed using BEAST 2 v2.6.5 ( 46 ) under a coalescent constant population model using tip dates and a mean clock rate constrained to 1.45 × 10 −6 substitutions per site per year, based on a previously reported evolution rate for the ST-147 lineage ( 47 ). Run parameters are detailed in SI Appendix .…”

In-host evolution of classic to convergentKlebsiella pneumoniaesequence type 147 isolates and impact of associated capsular changes on different morphotypes

“…The isolates characterized in this study belong to K. pneumoniae ST147, an international high-risk clonal lineage [47,48] suggest in-host evolution and a "transition" from the classic to the convergent pathotype.…”