Convergence and global molecular epidemiology of Klebsiella pneumoniae plasmids harbouring the iuc3 virulence locus: a population genomic analysis.

Abstract

BACKGROUND: Klebsiella pneumoniae is an important pathogen of humans and animals. In the past five years, increasing reports of convergent strains that carry both virulence factors and antimicrobial resistance genes (ARGs) have raised serious public health concerns. The aim of this study is to describe the global diversity of plasmids carrying iuc3 (a key virulence factor in K pneumoniae associated with pigs and clinical isolates) from diverse settings, and their role in the emergence of convergent strains through hybridisation with plasmids carrying ARGs. METHODS: This population genomic analysis study was designed to describe both the global and local diversity of iuc3-carrying plasmids from diverse sources, and the co-occurrence of iuc3 with ARGs. We used all 4148 Klebsiella spp isolates from two large One-Health studies (SpARK, Italy, and OH-DART, Thailand), including 191 Klebsiella isolates from pigs, 635 from clinical isolates, 1040 from hospital and community carriage, and 2282 from other sources. Short-read sequencing of Klebsiella isolates was performed as part of the SpARK study. We sequenced Klebsiella isolates from the OH-DART (MicrobesNG, Birmingham, UK; HiSeq and NovaSeq, Illumina San Diego, CA, USA; GridION, Oxford Nanopore Technologies, Oxford, UK) and SpARK (MinION or GridION, Oxford Nanopore Technologies, Oxford, UK) studies. We also retrieved plasmid sequences carrying iuc3 from the National Centre for Biotechnology Information (NCBI). To ascertain the degree of diversity, evolutionary dynamics, and structuring across ecological and geographical axes, we detected ARGs and virulence loci, analysed clustering patterns and generated approximate maximum-likelihood phylogenetic trees. FINDINGS: We identified 48 K pneumoniae isolates with iuc3 in the SpARK data and 79 in the OH-DART data. Three (2·4%) of these 127 isolates were from clinical sources, 73 (57·5%) were from pig or pork meat. iuc3 isolates corresponded to multiple (n=47) host sequence types (STs), with ST35, ST45, ST881, ST25, and ST967 harbouring iuc3 in both datasets. We generated hybrid assemblies for 44 (SpARK) and 36 (OH-DART) isolates, plus a single iuc3 isolate from Germany. 53 (65·4%) of these isolates were from pigs, three (3·7%) from clinical sources, and 25 (30·9%) from other sources. There were an additional 48 iuc3 positive isolates from our collections for which only short read data was available. A single iuc3-positive Klebsiella oxytoca isolate from a pig farm was detected in the SpARK data, which was also sequenced. We identified 330 iuc3-positive isolates and 58 iuc3-carrying plasmid assemblies from NCBI, of which 83 (21·4%) were from clinical sources, 120 from pigs (30·9%), and 185 (47·7%) from other sources or of unknown provenance. These isolates were from K pneumoniae except two isolates of Klebsiella quasipneumoniae subsp similipneumoniae and one of Enterobacter hormaechei. The combined dataset of 517 iuc3 plasmids ranged in size from 110 375 bp to 365 580 bp and mostly corresponded to multiple IncFIB(K) and IncFII replicon types. We found seven convergent K pneumoniae plasmids in the Thai data: six from fresh markets and one from a neighbouring hospital. These plasmids emerged through the hybridisation of cocirculating iuc3 plasmids and plasmids encoding extended-spectrum β-lactamases (ESBLs), although none of these seven plasmids carried genes encoding carbapenemases. We also identified putative cocirculating parental plasmids carrying iuc3 and ESBL-encoding genes. Clustering and phylogenetic analysis resolved the iuc3 plasmid sequences into three groups, which were consistent using both complete plasmid sequences (n=139) and short-read data (n=517). In the complete plasmid sequence data, 66 strains contained group 1 plasmids, 38 strains contained group 2 plasmids, and 35 strains contained group 3 plasmids. Group 3 plasmids are mostly carried by isolates circulating in hospitals throughout Asia, with occasional examples in Europe and elsewhere, and carry multiple ARGs and potential virulence factors. By contrast, group 1 plasmids are commonly carried by porcine isolates in Europe, and group 2 are a heterogeneous mixture of geographical and ecological sources. INTERPRETATION: Plasmid hybridisation occurs frequently outside of the health-care environment and can lead to the convergence of resistance and virulence traits. Generating complete plasmid sequences from regional population-scale samples facilitates the identification of convergent plasmids and their putative parental plasmids. Three robust groups of iuc3 plasmids were resolved, which show both epidemiological and geographical differences; one of these groups was associated with clinical isolates in Asia and warrants targeted plasmid surveillance. FUNDING: UKRI, JPIAMR, Evolution Education Trust, and a Schlumberger Foundation Fellowship.