Improved molecular characterization of the Klebsiella oxytoca complex reveals the prevalence of the kleboxymycin biosynthetic gene cluster

Abstract

As part of ongoing studies with clinically relevant Klebsiella spp., we characterized the genomes of three clinical GES-5-positive ST138 strains originally identified as Klebsiella oxytoca. blaOXY gene, average nucleotide identity and phylogenetic analyses showed the strains to be Klebsiella michiganensis. Affiliation of the strains to ST138 led us to demonstrate that the current multi-locus sequence typing scheme for K. oxytoca can be used to speciate members of this genetically diverse complex of bacteria. The strains encoded the kleboxymycin biosynthetic gene cluster (BGC), previously only found in K. oxytoca strains and one strain of Klebsiella grimontii. The finding of this BGC, associated with antibiotic-associated haemorrhagic colitis, in K. michiganensis led us to carry out a wide-ranging study to determine the prevalence of this BGC in Klebsiella spp. Of 7,170 publicly available Klebsiella genome sequences screened, 88 encoded the kleboxymycin BGC. All BGC-positive strains belonged to the K. oxytoca complex, with strains of four (K. oxytoca, K. pasteurii, K. grimontii, K. michiganensis) of the six species of the complex found to encode the complete BGC. In addition to being found in K. grimontii strains isolated from preterm infants, the BGC was found in K. oxytoca and K. michiganensis metagenome-assembled genomes recovered from neonates. Detection of the kleboxymycin BGC across the K. oxytoca complex may be of clinical relevance and this cluster should be included in databases characterizing virulence factors, in addition to those characterizing BGCs. IMPACT STATEMENTMembers of the Klebsiella oxytoca complex are difficult to speciate using phenotypic and chemotaxonomic methods. Consequently, many genomes deposited in public databases are misclassified as K. oxytoca. Here we demonstrate that the current multi-locus sequence typing (MLST) system for the complex can be used to accurately speciate many strains, which will be of use to clinical laboratories in resource-limited settings which rely on the MLST scheme for typing and epidemiological tracking of isolates. In addition, extended analyses of the genomes of Klebsiella spp. have revealed the kleboxymycin biosynthetic gene cluster (BGC) is restricted to species of the Klebsiella oxytoca complex (K. oxytoca, K. michiganensis, K. pasteurii and K. grimontii). Species- and/or gene-specific differences in the clusters sequences may be relevant to virulence of K. oxytoca and related species. The finding of the kleboxymycin BGC in the preterm infant gut microbiota may have implications for disease presentation in a subset of neonates.