Recognition of two genetic groups in the Klebsiella oxytoca taxon on the basis of chromosomal beta-lactamase and housekeeping gene sequences as well as ERIC-1 R PCR typing

Description of Klebsiella spallanzanii sp. nov. and of Klebsiella pasteurii sp. nov

Klebsiella oxytoca causes opportunistic human infections and post-antibiotic haemorrhagic diarrhea. This Enterobacteriaceae species is genetically heterogeneous and is currently subdivided into seven phylogroups (Ko1 to Ko4 and Ko6 to Ko8). Here we investigated the taxonomic status of phylogroups Ko3 and Ko4. Genomic sequence-based phylogenetic analyses demonstrate that Ko3 and Ko4 formed well-defined sequence clusters related to, but distinct from, Klebsiella michiganensis (Ko1), K. oxytoca (Ko2), K. huaxiensis (Ko8), and K. grimontii (Ko6). The average nucleotide identity (ANI) of Ko3 and Ko4 were 90.7% with K. huaxiensis and 95.5% with K. grimontii, respectively. In addition, three strains of K. huaxiensis, a species so far described based on a single strain from a urinary tract infection patient in China, were isolated from cattle and human feces. Biochemical and MALDI-ToF mass spectrometry analysis allowed differentiating Ko3, Ko4, and Ko8 from the other K. oxytoca species. Based on these results, we propose the names Klebsiella spallanzanii for the Ko3 phylogroup, with SPARK_775_C1^T (CIP 111695^T and DSM 109531^T) as type strain, and Klebsiella pasteurii for Ko4, with SPARK_836_C1^T (CIP 111696^T and DSM 109530^T) as type strain. Strains of K. spallanzanii were isolated from human urine, cow feces, and farm surfaces, while strains of K. pasteurii were found in fecal carriage from humans, cows, and turtles.

Phylogenetic lineages, clones and β-lactamases in an international collection of Klebsiella oxytoca isolates non-susceptible to expanded-spectrum cephalosporins

OBJECTIVES

The objective of this study was to examine Klebsiella oxytoca clonal and phylogenetic diversity, based on an international collection of carriage isolates non-susceptible to expanded-spectrum cephalosporins (ESCs).

METHODS

The study material comprised 68 rectal carriage K. oxytoca isolates non-susceptible to ESCs recovered in 2008-11 from patients in 14 hospitals across Europe and Israel. ESC resistance was tested phenotypically; genes encoding ESBLs, AmpC cephalosporinases and carbapenemases were amplified and sequenced. The isolates were typed by PFGE and MLST, followed by sequencing of blaOXY genes.

RESULTS

MLST and PFGE distinguished 34 STs and 47 pulsotypes among the isolates, respectively. Six STs were split into several pulsotypes each. Five STs were more prevalent (n = 2-9) and occurred in several countries each, including ST2, ST9 and ST141, which belong to a growing international clonal complex (CC), CC2. Four phylogenetic lineages were distinguished, each with another type of chromosomal OXY-type β-lactamase. Three of these, with OXY-1/-5, OXY-2 types and OXY-4, corresponded to previously described phylogroups KoI, KoII and KoIV, respectively. A single isolate from Israel represented a distinct lineage with a newly defined OXY-7 type. The phylogroups showed interesting differences in mechanisms of ESC resistance; KoI strains rarely overexpressed the OXY enzymes but commonly produced ESBLs, whereas KoII strains often were OXY hyperproducers and carried ESBLs much less frequently. AmpCs (DHA-1) and carbapenemases (VIM-1) occurred sporadically.

CONCLUSIONS

The study confirmed the high genetic diversity of the collection of K. oxytoca ESC-non-susceptible isolates, composed of phylogroups with distinct types of OXY-type β-lactamases, and revealed some STs of broad geographical distribution.

The complete nucleotide sequence of the multi-drug resistance-encoding IncL/M plasmid pACM1

The 89,977 bp nucleotide sequence of pACM1, isolated from a 1993 outbreak strain of cephalosporin-resistant Klebsiella oxytoca, has been completed and assigned GenBank accession number KJ541681. The plasmid has a single 31,842 bp mosaic multi-drug resistance-encoding (MDR) region comprising the mer resistance module of Tn1696, two integrons with a total of seven cassettes, one complete copy each of IS1R and IS26, and the bla(SHV-5)-carrying Tn2003 (with defective IS26 termini), all within a Tn1721-like element inserted into the mucB gene of the IncL/M plasmid backbone. The Tn1721-Tn1696 combination resembles sequence found in the chromosomal MDR islands of some Acinetobacter baumannii isolates. Among the completely sequenced IncL/M resistance plasmids, the Tn1721-based MDR region is unique, but data from older studies suggest that this type of plasmid was widespread in the 1990s. Since resistance gene dosage is affected by plasmid copy number (PCN), we used a relatively simple new "efficiency-corrected" qPCR assay to measure the PCN of pACM1. There are approximately three copies per chromosome in an Escherichia coli DH5α host, and two in the original Klebsiella oxytoca isolate. We could not find similar PCN data for other medically important plasmids for comparison. The study of this plasmid property and its effect on resistance levels should be facilitated in the future by the availability of qPCR instruments and complete genome sequences.

Comparative genomics of an IncA/C multidrug resistance plasmid from Escherichia coli and Klebsiella isolates from intensive care unit patients and the utility of whole-genome sequencing in health care settings

The IncA/C plasmids have been implicated for their role in the dissemination of β-lactamases, including gene variants that confer resistance to expanded-spectrum cephalosporins, which are often the treatment of last resort against multidrug-resistant, hospital-associated pathogens. A bla(FOX-5) gene was detected in 14 Escherichia coli and 16 Klebsiella isolates that were cultured from perianal swabs of patients admitted to an intensive care unit (ICU) of the University of Maryland Medical Center (UMMC) in Baltimore, MD, over a span of 3 years. Four of the FOX-encoding isolates were obtained from subsequent samples of patients that were initially negative for an AmpC β-lactamase upon admission to the ICU, suggesting that the AmpC β-lactamase-encoding plasmid was acquired while the patient was in the ICU. The genomes of five E. coli isolates and six Klebsiella isolates containing bla(FOX-5) were selected for sequencing based on their plasmid profiles. An ∼ 167-kb IncA/C plasmid encoding the FOX-5 β-lactamase, a CARB-2 β-lactamase, additional antimicrobial resistance genes, and heavy metal resistance genes was identified. Another FOX-5-encoding IncA/C plasmid that was nearly identical except for a variable region associated with the resistance genes was also identified. To our knowledge, these plasmids represent the first FOX-5-encoding plasmids sequenced. We used comparative genomics to describe the genetic diversity of a plasmid encoding a FOX-5 β-lactamase relative to the whole-genome diversity of 11 E. coli and Klebsiella isolates that carry this plasmid. Our findings demonstrate the utility of whole-genome sequencing for tracking of plasmid and antibiotic resistance gene distribution in health care settings.

Multidrug and heavy metal-resistant Raoultella planticola isolated from surface water

A surface water isolate of Raoultella sp. having both multidrug- and multimetal-resistant ability was isolated and identified as Raoultella planticola. R. planticola displayed resistance to 15 drugs like ampicillin, amoxicillin/clavulanic acid, aztreonam, erythromycin, imipenem, oxacillin, pefloxacin, penicillin, piperacillin, piperacillin/tazobactam, rifampin, sulbactam/cefoperazone, ticarsillin, ticarsillin/clavulanic acid, vancomycin, and to 11 heavy metals like aluminum, barium, copper, iron, lead, lithium, manganese, nickel, silver, strontium, and tin. The multidrug and multi-metal-resistant R. planticola may remain present in the environment for a long time. Due to a possible health risk of these pathogenic bacteria, a need exists for an accurate assessment of their acquired resistance to multiple drugs and metals.

Evaluation of three phenotypic identification systems for clinical isolates of Raoultella ornithinolytica

Raoultella spp. have recently been separated from the genus Klebsiella based on their molecular characteristics. It was discovered that Raoultella ornithinolytica can be misidentified as Klebsiella oxytoca by commonly used phenotypic identification systems. Therefore, this study evaluated the ability of three phenotypic systems to identify R. ornithinolytica compared with the genotypic methods sequence-specific primer PCR (SSP-PCR), 16S rRNA gene sequence analysis using the MicroSeq 500 system16S rDNA bacterial identification system or comparison with GenBank sequences using blast. The phenotypic systems examined in this study were the VITEK 2 GN ID card, the MicroScan Neg Combo 32 panel and API 20E. The SSP-PCR panel was able to distinguish the R. ornithinolytica reference strain from other Raoultella spp. and K. oxytoca. Of the 27 isolates identified as R. ornithinolytica by SSP-PCR, VITEK 2 identified all of them as R. ornithinolytica. MicroScan and API identified 25 isolates (92.6%) and 24 isolates (88.9%) as K. oxytoca, respectively. These isolates were ornithine decarboxylase (ODC) negative in all three phenotypic systems. MicroSeq 500 identified 24 isolates (88.9%) as R. ornithinolytica, whereas GenBank identification was heterogeneous. Of the 68 isolates identified as K. oxytoca by SSP-PCR, 66 isolates (97.1%) were identified as K. oxytoca by VITEK 2, MicroScan and API. MicroScan and API require additional biochemical tests to differentiate between ODC-negative R. ornithinolytica and K. oxytoca.

Persistence and spread of qnr, extended-spectrum beta-lactamase, and ampC resistance genes in the digestive tract of chickens

The aim of this assay was to develop an experimental model of digestive colonization of chickens with bacteria harboring qnr, extended-spectrum beta-lactamase, or ampC genes. Specific pathogen-free chickens were orally inoculated with two Escherichia coli strains containing either the plasmid pMG252 bearing bla(FOX) and qnrA genes, or pMG298 bearing bla(CTX-M) and qnrB genes. Analysis of strains isolated from fecal samples showed that the two strains were able to persist for several weeks in the digestive flora of inoculated birds and could rapidly spread to noninoculated ones. However, the multi-resistant isolates were maintained as a small proportion of the overall enterobacterial population. The qnr, extended-spectrum beta-lactamase, and ampC resistance genes could be transferred, in vivo, in the absence of selective pressure, to other chicken E. coli or Klebsiella pneumoniae isolates.

Klebsiella oxytoca: opportunistic infections in laboratory rodents

Opportunistic pathogens have become increasingly relevant as the causative agents of clinical disease and pathological lesions in laboratory animals. This study was conducted to evaluate the role of Klebsiella oxytoca as an opportunistic pathogen in laboratory rodents. Therefore, K. oxytoca-induced lesions were studied from 2004 to early 2006 in naturally infected rodent colonies maintained at The Jackson Laboratory (TJL), Bar Harbor, USA, the Animal Research Centre (Tierforschungszentrum, TFZ) of the University of Ulm, Germany and the Central Animal Facility (ZTM) of the Hannover Medical School, Germany. K. oxytoca infections were observed in substrains of C3H/HeJ mice, which carry the Tlr4(Lps-d) allele; in LEW.1AR1-iddm rats, the latter being prone to diabetes mellitus; in immunodeficient NMRI-Foxn1(nu) mice; and in mole voles, Ellobius lutescens. The main lesions observed were severe suppurative otitis media, urogenital tract infections and pneumonia. Bacteriological examination revealed K. oxytoca as monocultures in all cases. Clonality analysis performed on strains isolated at the ZTM and TFZ (serotyping, pulse field gel electrophoresis [PFGE], enterobacterial repetitive intergenic consensus (ERIC) polymerase chain reaction, sequencing of 16S rRNA and rpoB genes) revealed that the majority of bacteria belonged to two clones, one in each facility, expressing the capsule type K55 (ZTM) or K72 (TFZ). Two strains, one isolated at the ZTM and one at the TFZ, showed different PFGE and ERIC pattern than all other isolates and both expressed capsule type K35. In conclusion, K. oxytoca is an opportunistic pathogen capable of inducing pathological lesions in different rodent species.

Phylogeny and species identification of the family Enterobacteriaceae based on dnaJ sequences

Phylogenetic relations within the family Enterobacteriaceae were analyzed using partial dnaJ sequences of 165 strains belonging to 93 species from 27 enterobacterial genera. The dnaJ phylogeny was in relative agreement with that constructed by 16S rDNA sequences, but more monophyletic groups were obtained from the dnaJ tree than from the 16S rDNA tree. The degree of divergence of the dnaJ gene was approximately 6 times greater than that of 16S rDNA. Also, the dnaJ gene showed the most discriminatory power in comparison with tuf and atpD genes, facilitating clear differentiation of any 2 enterobacterial species by dnaJ sequence analysis. The application of dnaJ sequences to the identification was confirmed by assigning 72 clinical isolates to the correct enterobacterial species. Our data indicate that analysis of the dnaJ gene sequences can be used as a powerful marker for phylogenetic study and identification at the species level of the family Enterobacteriaceae.

Molecular characterization of beta-lactamase genes blaA and blaB of Yersinia enterocolitica biovar 1A

The beta-lactamase genes blaA and blaB were detected by PCR amplification in strains of Yersinia enterocolitica biovar 1A isolated from India, Germany, France and the USA. Both genes were detected in all strains. Polymerase chain reaction-restriction fragment length polymorphism revealed genetic heterogeneity in blaA but not in blaB. Cluster analysis of blaA restriction profiles grouped the strains into three groups. The blaA gene of Y. enterocolitica biovar 1A showed a high degree of sequence homology to that of Y. enterocolitica 8081 (biovar 1B) and Y. enterocolitica Y-56 (biovar 4), whereas homology was low with class A beta-lactamase genes of other members of the family Enterobacteriaceae. The pI 8.7 of enzyme Bla-A of Y. enterocolitica biovar 1A was similar to that of biovars 2, 3 and 4. The enzyme Bla-B focused at 6.8 and 7.1, indicating that biovar 1A strains produced a 'B-like' enzyme. This is the first study to have investigated the genetic heterogeneity of the beta-lactamase genes of Y. enterocolitica.

Six groups of the OXY beta-Lactamase evolved over millions of years in Klebsiella oxytoca

The diversity and evolution of the class A OXY beta-lactamase from Klebsiella oxytoca were investigated and compared to housekeeping gene diversity. The entire bla(OXY) coding region was sequenced in 18 clinical isolates representative of the four K. oxytoca beta-lactamase gene groups bla(OXY-1) to bla(OXY-4) and of two new groups identified here, bla(OXY-5) (with four isolates with pI 7.2 and one with pI 7.7) and bla(OXY-6) (with four isolates with pI 7.75 and three with pI 8.1). Genes bla(OXY-5) and bla(OXY-6) showed 99.8% within-group nucleotide similarity but differed from each other by 4.2% and from bla(OXY-1), their closest relative, by 2.5% and 2.9%, respectively. Antimicrobial susceptibility to beta-lactams was similar among OXY groups. Nucleotide sequence diversity of the 16S rRNA (1,454 bp), rpoB (940 bp), gyrA (383 bp), and gapDH (573 bp) genes was in agreement with the beta-lactamase gene phylogeny. Strains with bla(OXY-1), bla(OXY-2), bla(OXY-3), bla(OXY-4), and bla(OXY-6) genes formed five phylogenetic groups, named KoI, KoII, KoIII, KoIV, and KoVI, respectively. Isolates harboring bla(OXY-5) appeared to represent an emerging lineage within KoI. We estimated that the bla(OXY) gene has been evolving within K. oxytoca for approximately 100 million years, using as calibration the 140-million-year estimation of the Escherichia coli-Salmonella enterica split. These results show that the bla(OXY) gene has diversified along K. oxytoca phylogenetic lines over long periods of time without concomitant evolution of the antimicrobial resistance phenotype.

Diversity and evolution of the class A chromosomal beta-lactamase gene in Klebsiella pneumoniae

We investigated the diversity of the chromosomal class A beta-lactamase gene in Klebsiella pneumoniae in order to study the evolution of the gene. A 789-bp portion was sequenced in a panel of 28 strains, representative of three phylogenetic groups, KpI, KpII, and KpIII, recently identified in K. pneumoniae and of different chromosomal beta-lactamase variants previously identified. Three groups of sequences were found, two of them corresponding to the families SHV (pI 7.6) and LEN (pI 7.1), respectively, and one, more heterogeneous, corresponding to a new family that we named OKP (for other K. pneumoniae beta-lactamase). Levels of susceptibility to ampicillin, cefuroxime, cefotaxime, ceftazidime, and aztreonam and inhibition by clavulanic acid were similar in the three groups. One new SHV variant, seven new LEN variants, and four OKP variants were identified. The OKP variants formed two subgroups based on nucleotide sequences, one with pIs of 7.8 and 8.1 and the other with pIs of 6.5 and 7.0. The nucleotide sequences of the housekeeping genes gyrA, coding for subunit A of gyrase, and mdh, coding for malate dehydrogenase, were also determined. Phylogenetic analysis of the three genes studied revealed parallel evolution, with the SHV, OKP, and LEN beta-lactamase families corresponding to the phylogenetic groups KpI, KpII, and KpIII, respectively. This correspondence was fully confirmed for 34 additional strains in PCR assays specific for the three beta-lactamase families. We estimated the time since divergence of the phylogenetic groups KpI and KpIII at between 6 and 28 million years, confirming the ancient presence of the beta-lactamase gene in the genome of K. pneumoniae.

Diversity and evolution of the class A chromosomal beta-lactamase gene in Klebsiella pneumoniae

We investigated the diversity of the chromosomal class A beta-lactamase gene in Klebsiella pneumoniae in order to study the evolution of the gene. A 789-bp portion was sequenced in a panel of 28 strains, representative of three phylogenetic groups, KpI, KpII, and KpIII, recently identified in K. pneumoniae and of different chromosomal beta-lactamase variants previously identified. Three groups of sequences were found, two of them corresponding to the families SHV (pI 7.6) and LEN (pI 7.1), respectively, and one, more heterogeneous, corresponding to a new family that we named OKP (for other K. pneumoniae beta-lactamase). Levels of susceptibility to ampicillin, cefuroxime, cefotaxime, ceftazidime, and aztreonam and inhibition by clavulanic acid were similar in the three groups. One new SHV variant, seven new LEN variants, and four OKP variants were identified. The OKP variants formed two subgroups based on nucleotide sequences, one with pIs of 7.8 and 8.1 and the other with pIs of 6.5 and 7.0. The nucleotide sequences of the housekeeping genes gyrA, coding for subunit A of gyrase, and mdh, coding for malate dehydrogenase, were also determined. Phylogenetic analysis of the three genes studied revealed parallel evolution, with the SHV, OKP, and LEN beta-lactamase families corresponding to the phylogenetic groups KpI, KpII, and KpIII, respectively. This correspondence was fully confirmed for 34 additional strains in PCR assays specific for the three beta-lactamase families. We estimated the time since divergence of the phylogenetic groups KpI and KpIII at between 6 and 28 million years, confirming the ancient presence of the beta-lactamase gene in the genome of K. pneumoniae.

Genetic and biochemical characterization of the chromosomal class A beta-lactamases of Raoultella (formerly Klebsiella) planticola and Raoultella ornithinolytica

Enterobacterial strains of Raoultella spp. display a penicillinase-related beta-lactam resistance pattern suggesting the presence of a chromosomal bla gene. From whole-cell DNA of Raoultella planticola strain ATCC 33531(T) and Raoultella ornithinolytica strain ATCC 31898(T), bla genes were cloned and expressed into Escherichia coli. Each gene encoded an Ambler class A beta-lactamase, named PLA-1 and ORN-1 for R. planticola and R. ornithinolytica, respectively. These beta-lactamases (291 amino acids), with the same pI value of 7.8, had a shared amino acid identity of 94%, 37 to 47% identity with the majority of the chromosome-encoded class A beta-lactamases previously described for Enterobacteriaceae, and 66 to 69% identity with the two beta-lactamases LEN-1 and SHV-1 from Klebsiella pneumoniae. However, the highest identity percentage (69 to 71%) was found with the plasmid-mediated beta-lactamase TEM-1. PLA-1, which displayed very strong hydrolytic activity against penicillins, also displayed significant hydrolytic activity against cefepime and, to a lesser extent, against cefotaxime and aztreonam, but there was no hydrolytic activity against ceftazidime. Such a substrate profile suggests that the Raoultella beta-lactamases PLA-1 and ORN-1 should be classified into the group 2be of the beta-lactamase classification of K. Bush, G. A. Jacoby, and A. A. Medeiros (Antimicrob. Agents Chemother. 39:1211-1233, 1995). The highly homologous regions upstream of the bla(PLA-1A) and bla(ORN-1A) genes comprised a nucleotide sequence identical to the -35 region and another one very close to the -10 region of the bla(LEN-1) gene. From now on, as the bla gene sequences of the most frequent Raoultella and Klebsiella species are available, the bla gene amplification method can be used to differentiate these species from each other, which the biochemical tests currently carried out in the clinical laboratory are unable to do.

New Klebsiella oxytoca beta-lactamase genes bla(OXY-3) and bla(OXY-4) and a third genetic group of K oxytoca based on bla(OXY-3)

The two genetic groups (oxy-1 and oxy-2) previously identified in the Klebsiella oxytoca taxon are recognizable by four independent molecular markers: (i). ERIC-1R profiles, (ii). 16S ribosomal DNA (rDNA) signature sequences, (iii). singular nucleotides in a defined fragment of the rpoB gene, and (iv) the type of the strain's bla(OXY) gene (i.e., bla(OXY-1) or bla(OXY-2)). K. oxytoca strains SG266 and SG271 could not be classified into these genetic groups based on their ERIC-1R profile and bla(OXY) gene sequence. With regard to the gene identity percentages between the bla(OXY-1) and bla(OXY-2) gene groups (86.8% +/- 0.4%) and within a bla(OXY) gene group (>99%), it was concluded that the bla(OXY) gene of strain SG271 was representative of a new bla(OXY) gene group (bla(OXY-3)), since the mean identity percentages between it and the two bla(OXY) gene groups were 85.5% +/- 0.2% and 84.4% +/- 0.4%, respectively. Since the corresponding percentages were 95.0% +/- 0.4% and 86.2% +/- 0.3% for strain SG266, it was impossible to classify its bla(OXY) gene, which was therefore named bla(OXY-4). The 16S rDNA signature sequences of the two strains could be determined only after cloning experiments. The SG266 clones displayed the same signature sequence as that of the genetic group oxy-1, whereas the SG271 clones displayed three different 16S rDNA signature sequences that also differed from those of the two genetic groups. Singular nucleotides were found within the rpoB sequence of the two strains, allowing for their distinction from the two genetic groups. All of these results, combined with those previously obtained by the ERIC-1R PCR method, indicate that strain SG271 is representative of a new K. oxytoca genetic group (oxy-3), whereas strain SG266 could not be classified.

Enterobacterial repetitive intergenic consensus 1R PCR assay for detection of Raoultella sp. isolates among strains identified as Klebsiella oxytoca in the clinical laboratory

The enterobacterial repetitive intergenic consensus 1R PCR method, which provided recognizable profiles for reference strains of the three species of Raoultella and the two genetic groups of Klebsiella oxytoca, was applied to 19 clinical isolates identified as K. oxytoca. By this method, as confirmed by species-specific gene sequencing, two Raoultella ornithinolytica and two unclassifiable K. oxytoca isolates were identified.