Next Article in Journal
Inanimate Surfaces and Air Contamination with Multidrug Resistant Species of Staphylococcus in the Neonatal Intensive Care Unit Environment
Next Article in Special Issue
Correction: Theilacker et al. PCV13 Vaccination of Adults against Pneumococcal Disease: What We Have Learned from the Community-Acquired Pneumonia Immunization Trial in Adults (CAPiTA). Microorganisms 2022, 10, 127
Previous Article in Journal
Biodegradation of Tetracycline Antibiotics by the Yeast Strain Cutaneotrichosporon dermatis M503
Previous Article in Special Issue
23-Valent Pneumococcal Polysaccharide Vaccination Does Not Prevent Community-Acquired Pneumonia Hospitalizations Due to Vaccine-Type Streptococcus pneumoniae
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

WGS-Based Phenotyping and Molecular Characterization of the Resistome, Virulome and Plasmid Replicons in Klebsiella pneumoniae Isolates from Powdered Milk Produced in Germany

1
Friedrich-Loeffler-Institut, Institute of Bacterial Infections and Zoonoses (IBIZ), Naumburger Str. 96a, 07743 Jena, Germany
2
Institute for Infectious Diseases and Infection Control, Jena University Hospital, Jena University, 07743 Jena, Germany
3
Department of Bacteriology, Immunology and Mycology, Faculty of Veterinary Medicine, Benha University, Toukh 13736, Egypt
4
Department of Microbiology and Biotechnology, Max Rubner-Institut, 24103 Kiel, Germany
5
Research Campus Infectognostics, Philosophenweg 7, 07743 Jena, Germany
*
Author to whom correspondence should be addressed.
Microorganisms 2022, 10(3), 564; https://doi.org/10.3390/microorganisms10030564
Submission received: 17 December 2021 / Revised: 3 March 2022 / Accepted: 3 March 2022 / Published: 5 March 2022

Abstract

:
The emergence of Klebsiella pneumoniae (K. pneumoniae) in German healthcare is worrying. It is not well-investigated in the veterinary world and food chains. In the current study, antibiotic susceptibility profiles of 24 K. pneumoniae strains isolated from powdered milk samples produced in Germany were investigated by a microdilution test. Next-generation sequencing (NGS) was applied to identify genomic determinants for antimicrobial resistance (AMR), virulence-associated genes and plasmids replicons. All isolates were susceptible to the majority (14/18) of tested antibiotics. Resistance to colistin, fosfomycin, chloramphenicol and piperacillin was found. The ambler class A ß-lactamase, blaSHV variants were identified in all isolates, of which blaSHV-187 was most prevalent and found in 50% of isolates. Single-nucleotide-variants of oqxA and oqxB conferring resistance to phenicol/quinolone were found in all isolates, and the oqxB17 was the most prevalent found in 46% of isolates. 67% of isolates harbored fosA genes; however, only one was fosfomycin-resistant. Two isolates harbored genes conferring resistance to colistin, despite being susceptible. The majority of identified virulome genes were iron uptake siderophores. Two enterobactins (entB, fepC), six adherence-related genes belonging to E. coli common pilus (ECP) and one secretion system (ompA gene) were found in all isolates. In contrast, yersiniabactin was found in two isolates. One ST23 strain was susceptible to all tested antibiotics, and harbored determinants discriminatory for hypervirulent strains, e.g., aerobactin, salmochelin, yersiniabactin, enterobactin and regulator of mucoid phenotype A genes that are highly associated with hypervirulent K. pneumoniae. The IncF plasmid family was found in all strains, while almost half of the isolates harbored Col440I-type plasmids and nine isolates harbored various Inc-type plasmids. The presence of K. pneumoniae carrying different resistomes and major virulent specific virulomes in powdered milk samples is alarming. This could threaten public health, particularly of neonates and infants consuming dried milk.

1. Introduction

Klebsiella pneumoniae (K. pneumoniae) is a Gram-negative, non-motile, usually capsule-forming bacterium. The opportunistic pathogen is associated with severe infections and high mortality rates in intensive care unit (ICU) patients [1]. It is one of the ‘ESKAPE’ pathogens, the most common and serious multidrug-resistant (MDR) pathogen in hospitalized patients worldwide [2]. K. pneumoniae has been found to be involved in a wide variety of nosocomial infections, such as sepsis, respiratory and urinary tract infections and wound infections [3]. It has also been isolated from a wide variety of foods, environmental sources such as soil, water and vegetation and dairy farm milk [4]. In recent years, a high prevalence of MDR K. pneumoniae (MDR-KP) and carbapenem-resistant K. pneumoniae (CRKP) has been observed in hospitalized persons worldwide, and has been associated with unprecedented public health problems [5]. K. pneumoniae has been frequently associated with mastitis in bovines [6], leading to high milk production loss, decreased quality and even high mortalities among affected cows [7,8]. The prevalence of this potentially lethal pathogen has been gradually increasing in dairy farms, the environment and the food chain [9].
In Germany, K. pneumoniae isolates developed resistance to third-generation cephalosporins and aminoglycosides in the late 20th century [10]. The high incidence of MDR strains has been confirmed in most German federal states [11,12]. The emergence of extended-spectrum beta-lactamases (ESBL) producing K. pneumoniae strains is on the rise in Germany, where blaCTX-M-15, blaCTX-M-3, blaSHV-11 and blaSHV-5 are the most frequently detected genes [13,14,15,16]. K. pneumoniae has been extensively investigated in humans and hypervirulent (hvKp) and extensively drug-resistant (XDR) strains have been isolated from patients [17,18]. However, research on K. pneumoniae is still an exception in veterinary and environmental health, and the risk of human infection regarding food consumption is only just under investigation [19]. Different virulence and pathogenicity factors contribute to K. pneumoniae pathogenesis and mediate infectivity, e.g., by means of adherence factors, siderophore activity, O-lipopolysaccharide (LPS) and K-capsular antigens [20]. However, many aspects of its pathogenicity are not yet clear.
The current study aimed at the phenotype and genomic characterization of antimicrobial resistance determinants, pathogenicity-associated genes and plasmid replicons in K. pneumoniae strains isolated from milk powder produced in Germany utilizing Next-generation sequence (NGS) technology.

2. Materials and Methods

2.1. Bacterial Isolates and Identification

All Klebsiella spp. isolates used in this study originated from dried milk produced via spray drying (20%) or rotating drum drying (80%) from two companies in Germany. Samples were collected from the end of the production line between 2005 and 2011, as proposed previously [21]. Samples of 10 g were briefly taken by the quality control laboratory of the producers, and tested for the presence of enterobacteria: Bacterial enrichment was carried out in double-buffered peptone water [peptone water 20.0 g L−l (Oxoid, Wesel, Germany), 3.5 g L−l Na2HPO4 (Merck), 1.5 g L−l KH2PO4 (Merck)] for 24 h at 30 °C with subsequent plating onto violet red bile dextrose (VRBD) (Merck KGaA, Darmstadt, Germany) agar and incubation for 24 h at 30 °C. VRBD plates showing bacterial growth were transferred to the Max Rubner-Institute. Colonies with different morphotypes were picked and phenotypically characterized by Gram-staining and testing by API 20E (BioMérieux, Nürtingen, Germany). Twenty-four K. pneumoniae strains isolated from powdered milk were sent to the Institute of Bacterial Infections and Zoonoses (IBIZ, Jena, Germany) for confirmation and typing. Samples metadata are shown in the Supplementary Table S1. All isolates were re-identified using Matrix-Assisted Laser Desorption/Ionization Time-of-Flight Mass Spectrometry (MALDI-TOF MS) at the species level with a log value > 2.300. Protein extraction from pure colonies of each sample was done as previously described [22]. The MALDI-TOF measurements were carried out using a Microflex LT instrument (Bruker Daltonics, Bremen, Germany), following the MALDI Biotyper manufacturer’s recommendation on the log score value of 0–3 for species identification. Score values between 2.300 and 3.000 were considered ‘highly probable for species identification’; values between 2.000 and 2.290 were considered ‘secure genus identification; values between 1.700 and 1.990 were considered ‘probable genus identification’ and values between 0 and 1.690 were not considered for identification. Whole-genome sequencing data were used to confirm the identity of the genera and the species of each strain using Kraken (v2.0.7_beta) [23]. Kraken 2 (v2.0.7_beta) in combination with the database MiniKraken (v2) were used to classify reads and assemblies and to check for possible contamination. At the genus and species levels, the first match (largest percentage proportion) was considered for identification.

2.2. Antibiotics Susceptibility Testing (AST)

The minimum inhibitory concentration (MIC) was determined by the broth microdilution method using an automated MICRONAUT-S system (Micronaut, MERLIN Diagnostics GmbH, Bornheim-Hersel, Germany) and MICRONAUT-S MDR MRGN-Screening MIC plates (Catalog Nr. E1-114-040) as previously described [24]. Strains were automatically classified as susceptible, intermediate and resistant by the MICRONAUTS software, which applies the MIC values according to the Clinical and Laboratory Standards Institute (CLSI, 2020) breakpoint guidelines available for K. pneumoniae [25]. The sensitivity of isolates was determined from a panel of 18 antibiotics at different concertations on the plates: ciprofloxacin (CIP), levofloxacin (LEV), amikacin (AMK), chloramphenicol (CMP), fosfomycin (FOS), tigecycline (TGC), colistin (COL), trimethoprim/sulfamethoxazole (T/S), piperacillin (PIP), piperacillin/tazobactam (PIT), cefotaxime (CTX), ceftazidime (CAZ), ceftazidime/avibactam (CAA), ceftolozane/tazobactam (CTA), cefepime (CEP), imipenem (IMP), meropenem (MER) and ertapenem (ERT).

2.3. WGS and In-Silico Detection of AMR Determinants, Virulome Genes and Plasmidome

Genomic DNA was extracted using the High Pure PCR Template Preparation Kit (Roche Diagnostics GmbH, Mannheim, Germany) according to the manufacturer’s instructions. Sequencing libraries were prepared with the Nextera XT DNA Library Prep Kit (Illumina, Inc., San Diego, CA, USA), followed by paired-end sequencing on an Illumina MiSeq sequencer (Illumina). To analyze raw sequencing data, the Linux-based pipeline WGSBAC (v2.0, https://gitlab.com/FLI_Bioinfo/WGSBAC, accessed on 28 February 2022) was used as previously described [24,26]. The pipeline performs a quality control of raw sequencing data utilizing FastQC (v0.11.7) files [27], and calculates raw sequencing coverage. For assembly, the pipeline uses Shovill (v1.0.4, https://github.com/tseemann/shovill, accessed on 28 February 2022), which is based on the SPAdes assembler [28]. The quality of the assembled genomes is checked with QUAST (v. 45.0.2) [29]. To identify contaminations, WGSBAC uses the sequence classifier Kraken 2 v1.1 [23] and the database Kraken2DB. For the prediction of virulence factors, WGSBAC utilizes ABRIcate (v0.8.10, https://github.com/tseemann/abricate, accessed on 28 February 2022) together with the Virulence Factor Database [30]. For the determination of genetic features leading to AMR, ABRIcate, together with the databases, Comprehensive Antibiotic Resistance Database (CARD) [31] and ResFinder [32], were applied. In addition, NCBI’s AMRFinderPlus tool [33] was used, together with the Klebsiella-specific parameters, to identify point mutations leading to AMR. Abricate against the PlasmidFinder database [34] was used to identify the presence of plasmid replicons. In-silico determination of classical multilocus sequence typing (MLST) was performed by WGSBAC based on the assembled genomes using the software mlst (v2.16.1, https://github.com/tseemann/mlst, accessed on 28 February 2022), which incorporates the species-specific scheme (K. pneumoniae). To construct a phylogenetic tree, Snippy (v4.3.6 https://github.com/tseemann/snippy, accessed on 28 February 2022) was used to identify core-genome single nucleotide polymorphisms (cgSNPs), together with the strain HS11286 (GenBank accession no. ASM24018v2) as a reference genome. WGSBAC finally used RAxML (Randomized Accelerated Maximum Likelihood) v8 [35] to reconstruct a phylogenetic tree based on the cgSNP alignment and the interactive Tree of Life (iTOL) v4 web tool (https://itol.embl.de/login.cgi, accessed on 28 February 2022) for visualization [36].

2.4. Data Availability

All study data are included in the article and supporting information. The data have also been submitted to the European Nucleotide Archive (ENA). The project accession number is PRJEB45776.

3. Results

3.1. Whole Genome Sequencing Data and MLST Analysis

Genome sequencing of twenty-four K. pneumoniae isolates yielded an average total number of 1,855,429 reads per isolate (range 1,597,112–2,155,888 Table S1). The isolates’ mean coverage was 86-fold (range from 73-fold to 95-fold). At the genus level, the first match (largest percentage proportion) for all isolates was always “Klebsiella” which was the classification for 91.52% of the reads on average (max 94.58%, min 86.57%). At the species level, the first match for all 24 isolates was always “K. pneumoniae.” Genome assembly yielded a genome size with a minimum of 5.305.557 bp and a maximum of 5.770.949 bp. The GC content was on average 57.08%. The mean N50 of the 24 assembled genomes was 27,275 bp (range 154,040–399,650 bp) (Table S1).
MLST analyzes of the 24 K. pneumoniae isolates allocated ten strains into a distinct sequence type (ST) and 14 strains failed to be assigned due to new alleles. Three strains were assigned to ST/1322, three were assigned to ST/567 and one isolate each to ST/23, ST/220, ST/502 and ST/1083. Of the 14 strains (58%) which initially could not be assigned to a distinct sequence type, five strains isolated in 2007, 2008, 2010 and 2011 were allocated to novel STs (New 1), and later on assigned to ST/5625. One isolate from 2010 was another novel ST (New 2), which is assigned to ST/6016. The remaining eight isolates, which were isolated in 2005, 2006 and 2009, showed the occurrence of new alleles, and were near to ST/299 with the presence of new alleles with less than 100% identity emerged in the locus pgi (New 3), and later on were assigned to ST6014 (Table 1). The phylogenetic analysis using core-genome SNP calling included additional eight K. pneumoniae isolates from our collection (DE-FLI1 to DE-FLI8). The resulting phylogenetic tree implies that the six strains of unknown ST are closely related to ST/1083 (DE-MRI4), while the eight strains with untyped ST are similar to ST 6 (DE-FLI8), as shown in Figure 1.

3.2. Phenotyping and AMR Determinants in K. pneumoniae Isolates

Antibiotic susceptibility testing of the 24 K. pneumoniae strains revealed that all strains were susceptible to 14 of 18 tested antibiotics, i.e., ciprofloxacin, levofloxacin, amikacin, the tetracycline derivative tigecycline, trimethoprim/sulfamethoxazole, cefepime, cefotaxime, ceftazidime, ceftazidime/avibactam, ceftolozane/tazobactam, piperacillin/tazobactam, imipenem, meropenem and ertapenem in-vitro. Resistance to colistin and fosfomycin was seen for isolates DE-MRI3 and DE-MRI17 obtained in 2005 and 2007, respectively. One-third (8 out of 24) of the isolates displayed resistance to piperacillin, and four isolates showed non-susceptibility (one resistant and three intermediate) to chloramphenicol.
No genes mediating resistance to aminoglycosides were found. All strains harbored at least one variant of the sulfhydryl variable (blaSHV) ambler class A β-lactamase. Seven blaSHV variants were identified, and blaSHV-187 was the most prevalent and found in 12 (50%) isolates, followed by blaSHV-11, blaSHV-1 and blaSHV-145, which were found in ten, nine and eight isolates, respectively. The variants of oqxA and oqxB, conferring resistance to phenicol/quinolone, were found in all isolates, despite only four isolates showing non-susceptibility to chloramphenicol. The oqxB17 variant was the most prevalent found in elven isolates, followed by oqxA7 and oqxA10, found in six isolates, while oqxB25 and oqxB26 were found in the same five isolates. The two genes pmrB_R256G and mcr-9 conferring resistance to colistin were found in two different isolates obtained from milk powder taken in 2005 and 2009, respectively. However, both isolates were susceptible to colistin. In contrast, the colistin-resistant strain did not carry any genes from the mcr family. More than half of the isolates, 67% (n = 16), harbored the fosA gene, which confers reduced susceptibility to fosfomycin; however, only one isolate was resistant to fosfomycin. Moreover, all isolates harbored genes for the antibiotic efflux pump, i.e., ramA the positive regulator of AcrAB-TolC, the emrR and emrD genes, Klebsiella_pneumoniae_KpnE, KpnF, KpnG, KpnH; Klebsiella_pneumoniae_acrA, acrB, acrD and Klebsiella_pneumoniae_OmpK37 (Table S1).

3.3. Characterization of Plasmid Replicons and Virulence-Associated Genes in K. pneumoniae

The plasmid replicons in the K. pneumoniae isolates were analyzed using Plasmid Finder. Ten different types of plasmid replicons were identified in the 24 K. pneumoniae isolates. Four replicons belonged to the IncF plasmid family, three belonged to Col-type plasmids, two belonged to the Inc-type plasmid family and one belonged to the repA. All isolates were found to be carrying plasmid replicons; however, the number of replicons was variable in each strain, and ranged from one to six. The strain ST502 (DE-MRI9) carried the highest number of plasmid replicons (n = 6) among the strains and contained the second-highest number of virulence-associated genes (n = 21) after the ST23 strain, which carried 29 virulence factors. The IncFIB(K)-1-Kpn3 plasmid replicon was predominant in 16 (67%) isolates, followed by IncFII_1_pKP91 and Col440I_1 that were found in 13 (54%) and 11 (46%) isolates, respectively. The IncFIB.Mar.-1-pNDM.Mar and IncHI1B-1-pNDM.MAR replicons were found in seven isolates, while the IncFIB.pKPHS1.-1-pKPHS1 replicon was found in four isolates and IncR-1 was found in three isolates. The ColRNAI-1 was found in two isolates, while Col440II_1 and RepA_1_pKPC.CAV1321 were found once (Table S1).
In the current study, 29 virulence-associated genes were identified in the 24 strains. The majority of identified virulence-associated genes belonged to the Iron uptake siderophores. The prototypic catecholate siderophore enterobactins (entB, fepC) were found in all isolates (100%), while entA was found in 21 (87.5%) isolates. Eleven phenolate siderophore yersiniabactin-related genes (ybtA, ybtE, ybtP, ybtQ, ybtS, ybtT, ybtU, ybtX, irp1, irp2, fyuA) were detected and all were found in the same two (8%) isolates. The citrate-hydroxamate siderophore aerobactins (iucA, iucB, iucC, iutA), the glycosylated salmochelins (iroB, iroC, iroD, iroN) and the regulator of mucoid phenotype A (rmpA2, rmpA, rmpC and rmpD) were found only in an ST23 strain (4%). The outer membrane protein A (ompA), which is related to the T6SS-II secretion system of Klebsiella and further six adherence virulence genes related to ECP from E. coli (ykgK.ecpR, yagV.ecpE, yagW.ecpD, yagX.ecpC, yagY.ecpB, yagZ.ecpA) were found in all isolates (100%) (Table S1). One ST23 strain (DE-MRI13) isolated in 2007 was susceptible to all tested antibiotics and harbored determinants discriminatory for hypervirulent K. pneumoniae (hvKp), e.g., aerobactin, salmochelin, yersiniabactin, enterobactin genes and regulator of mucoid phenotype A, while it was hypermucoviscosity negative and devoid of magA and wcaG genes (Table 2). Plasmids of the ST/23 strain (DE-MRI13) were compared with the pK2044 plasmid of an ST/23 isolate NC_006625 available at https://www.ncbi.nlm.nih.gov/nuccore/NC_006625.1 (accessed on 28 February 2022) regarding the content of rmpA/A2. Our ST/23 isolate (DE-MRI13) contains two contigs harboring INK factors. Contig 24 is 63.276bp long and harbours IncFIB(K)_1_Kpn3 and contig 20 is 84.546 bp long and contains IncHI1B_1_pNDM-MA. The rmpA2 was found only on the contig 20 of DE-MRI13 strain, while rmpA was found on the plasmid pK2044 of strain NC_006625 from NCBI.
The plasmid carrying the mcr-9 gene in the present study was compared with the plasmids pA2483 and pA2504 harboring mcr-9 found in carbapenem-resistant and colistin non-resistant Enterobacter cloacae complex isolates recently published by Kananizadeh et al. in 2020 [37]. The plasmids pA2483 and pA2504 are 288.696 bp and 276.927 bp in size, respectively. Both contain three INK factors: IncHI2A, IncHI2 and RepA. In the present study, the strain (DE-MRI21) harbours mcr-9 on contig 37. This contig exclusively contains the INK factor RepA and is 31.103 bp in size. Regarding genetic factors for AMR, the plasmid carrying the contig 37 contains exclusively mcr-9, while the other plasmids contain additional AMR genes conferring resistance to heavy metals.

4. Discussion

Klebsiella pneumoniae is a notorious pathogen throughout the health care system worldwide. It causes a wide range of severe infections that are often difficult to treat [38]. In animals, it has been implicated in bovine mastitis [6], as well as in mortalities among affected cows [7,8], resulting in substantial economic losses in the dairy industry [39]. K. pneumoniae was among the most frequently isolated bacterial species in milk substitution formulas for infants collected from 35 countries [40]. A high incidence of MDR strains has been reported in pasteurized milk and whole milk powder samples collected from retail shops in Mexico, representing a public health hazard [41]. ESBL-producing K. pneumoniae was isolated from food handlers who consumed unpasteurized milk and raw meat [42]. To assess whether the powdered milk might be an essential source for disseminating K. pneumoniae or resistance genes in Germany, we investigated 24 isolates obtained from the milk powder of powdered milk producers by WGS. Eight of the strains displayed resistance to piperacillin, one to colistin and one to fosfomycin. All isolates were susceptible to ciprofloxacin and levofloxacin (fluoroquinolones), amikacin (aminoglycoside), the tetracycline derivative tigecycline, trimethoprim/sulfamethoxazole, 3rd generation cephalosporins cefotaxime, ceftazidime, ceftazidime/avibactam and ceftolozane/tazobactam, cefepime (4th generation cephalosporin), piperacillin/tazobactam and to the carbapenems (imipenem, meropenem and ertapenem). Colistin is considered a last-resort antibiotic and is used to treat infections caused by MDR Gram-negative bacteria [43]. Therefore, the isolation of a strain resistant to colistin from milk powder is alarming and carries a potential risk for humans, particularly for neonates and infants consuming dried milk. Although none of the isolates recovered from milk powder showed an MDR pattern, several AMR genes mediating resistance to ß-lactams, efflux pumps, phenicol/quinolone, colistin and fosfomycin were identified.
In the current study, all isolates were historical isolates recovered from milk powder samples between 2005–2011 and had not been investigated regarding resistance profile, MLST, AMR and virulence-associated genes at the time. The investigation of historical isolates by WGS provides a large amount of data, which can help to elucidate the development of AMR in Germany and previously unknown resistance mechanisms. All isolates were obtained from the end product of the manufacturing process. However, it was not possible to determine the source of contamination, i.e., human, animal or environmental. As the milk is pasteurized and heat treatment is applied during concentration and drying, the survival of K. pneumoniae and the introduction of strains from animal origin, i.e., cow, is highly unlikely. However, recontamination in the production plant via personnel, circulating air and packing material may cause contamination in the end product. Only a few articles exploring ST in human isolates have been published to date, and none on animal isolates. In the current study, an ST/23 strain, which is associated with a hypervirulent K. pneumoniae type (hvKp) was isolated in 2007. The isolate harbored determinants discriminatory for hypervirulent strains, e.g., aerobactin, salmochelin, yersiniabactin, enterobactin and regulator of mucoid phenotype A genes; however, it was hypermucoviscosity negative and devoid of magA and wcaG genes. An ST/23 strain producing OXA-48 was isolated previously from tracheal secretions of hospitalized patients [44]. A hvKp ST268 strain was also isolated from water in northern Germany [45] and an hvKp ST2398 strain was isolated from a patient with liver abscess and endophthalmitis in 2016 [46]. An MDR ST/11 strain was reported in European mouflons (ovis orientalis musimon) in 2016 [47].
Carbapenems and colistin are considered the drugs of choice for treating MDR bacterial infections. Infections caused by carbapenemase-producing K. pneumoniae are increasing in Germany [48,49], and colistin-resistant strains carrying the mcr-1 gene were reported in clinical samples in 2017 [50], and from wastewater treatment plants [51]. In the current study, blaSHV variants ß-lactamase were found in all strains. Since the late 20th century, blaSHV-11 and blaSHV-5 are among the most frequently detected genes in K. pneumoniae of human origin in Germany [52]. The blaSHV variants were also identified in K. pneumoniae strains isolated from European mouflons [47], pet animals [53] and environmental samples collected from pig farms [54]. Although blaOXA-48 and blaCTX-M-15 are the most frequently identified resistance genes in K. pneumoniae strains isolated from patients [15,18,55] and pet animals [56,57] in the last decade in Germany, none were found in the current study. The presence of pmrB_R256G and mcr-9 genes did not confer colistin resistance in our samples. The mcr-9 is a mobilized and plasmid-mediated colistin resistance gene identified in Salmonella enterica [58]. It is the most widely disseminated gene of the mcr-family after mcr-1, identified in isolates from 40 countries across six continents [59]. Its ability to cause colistin resistance has been investigated in several studies, and the occurrence of two different frameshift mutations in mcr-9 probably leads to non-functional mcr-9 proteins [58]. However, the occurrence of frameshift mutation was not investigated in the mcr-9 positive strain. To the best of our knowledge, the mcr-9 was not reported in K. pneumoniae strains isolated from humans, veterinaries and the environment in Germany. However, colistin-resistant strains carrying the mcr-1 gene have been reported in clinical samples from leukemia patients [50] and wastewater treatment plants [51]. Our results agree with Kananizadeh and coworkers, who found colistin-susceptible Enterobacter cloacae harboring mcr-9 [37]. Comparison of the plasmids carrying mcr-9 in our study and the plasmids from Kananizadeh et al. [37] revealed that our strain harbored exclusively mcr-9 on contig 37. This contig exclusively contains INK factor RepA, while the plasmids described by Kananizadeh contain three INK factors, IncHI2A, IncHI2 and RepA, and contain additional AMR genes such as pcoS conferring resistance to Arsenic and terD, terZ and terW conferring resistance to Tellurium.
Of the 24 isolates, only one was resistant to fosfomycin, despite the presence of fosA in 16 isolates. Chromosomal fosA genes confer a high level of resistance to fosfomycin in Gram-negative bacteria [60]. However, it failed to mediate resistance in the current K. pneumoniae isolates. The fosA gene is widely distributed on the chromosomes of many Gram-negative bacteria. Fosfomycin has become a valid choice against MDR, extended-spectrum ß-lactamases (ESBL) and carbapenem-resistant Enterobacteriaceae, and it is recommended as an emerging treatment for infection caused by MDR bacterial pathogens [61,62]. Recently, it was reported that FosA inhibitors such as sodium phosphonoformate (PPF) could restore the activity of fosfomycin in chromosomally encoded FosA Gram-negative bacterial strains, including K. pneumoniae [63]. Moreover, a small-molecule inhibitor of FosA significantly potentiated fosfomycin activity against Gram-negative pathogens harboring the fosA gene, including K. pneumoniae [64]. Some of these FosA inhibitors are clinically approved antiviral agents. Exposure of K. pneumoniae strains to such inhibitors could explain the susceptibility of the strains to fosfomycin despite harboring the fosA gene.
Four main components are known to mediate infectivity and pathogenesis in K. pneumoniae, including adherence factors, siderophore activity, O-lipopolysaccharide (LPS) and K-capsular antigens [20]. Six adherence-related genes belonging to ECP (E. coli common pilus) were found in all strains; however, no Type I fimbriae was found. ECP is a pilus of EHEC O157:H7, necessary for the initial attachment to host epithelial cells and subsequent colonization [65]. Iron acquisition systems or siderophores are needed in K. pneumoniae to overcome host defense and for deep penetration of the tissues [66]. The majority of identified virulence genes belonged to iron uptake siderophores. Enterobactins were found in almost all isolates, while salmochelin, aerobactin and the regulator of mucoid phenotype A were found in one isolate and yersiniabactin in two. The siderophore aerobactin is a critical and dominant virulence factor produced by hypervirulent K. pneumoniae [67]. hvKp infections are currently increasing, and are considered a global threat [68,69]. Several virulence genes were found to contribute to the hypervirulent phenotype, including mucoviscosity-associated gene A (magA) and (wcaG) and regulator of mucoid phenotype A (rmpA and rmpA2) [70,71]. The rmpA2/A are highly associated with hvKp as they positively regulate the cps locus during the capsular polysaccharide synthesis, resulting in the hypercapsule phenotype [72], and were identified in the ST/23 strain in the current study. However, recent studies showed that aerobactin is the major hvKp-specific virulence factor [73,74] and is used to differentiate hvkp either alone or combined with salmochelin and rmpA2 [75], or with hypermucoviscosity [76]. The role of yersiniabactin in pathogenesis is still unknown [77]. However, it was found in 17.7% of K. pneumoniae strains isolated from blood cultures and urine from hospitalized patients in Munich [78].
Plasmids are the major factors included in the dissemination of pathogenicity and resistance-associated genes in bacteria. Their duplication and survival capacity depends on their replication determinants (replicon) [79]. The plasmid and replicons content in K. pneumoniae isolates were investigated. Plasmids were identified in almost all K. pneumoniae strains with a diverse replicon content. The IncF plasmid family members are the most frequent plasmid types in Enterobacteriaceae and were found in all strains. More than half of the isolates contained IncFIB and IncFII plasmids. IncF plasmids play a significant role in disseminating AMR, and are often associated with virulence genes in Enterobacteriaceae [80]. Although the IncF plasmids have been associated with the emergence of blaCTX-M-15 globally [81,82], none of the isolates in the current study harbored blaCTX variants. Two Inc-type plasmids were found, IncHI1B and IncR. The IncR plasmid is often responsible for the horizontal transmission of blaKPC to other Enterobacteriaceae by conjugation [1]. Half of the isolates contained the Col440I-type plasmid replicon, and two isolates contained ColRNAI. Resistance to carbapenems in XDR K. pneumoniae due to the blaNDM-1 is mediated by the acquisition of multi-replicons containing ColRNAI, IncFIB Col440I, IncFII and IncFIB plasmid [83], and all have been identified in the current study.

5. Conclusions

K. pneumoniae isolates recovered from the end of the production process of powdered milk samples in Germany revealed a wide variety of AMR genes, mostly mediating resistance to ß-lactams and virulence-associated genes, which were mostly siderophores. The isolates showed susceptibility to the majority (14 of 18) of tested antibiotics. The appearance of isolates resistant to colistin and fosfomycin, chloramphenicol and piperacillin is alarming. Plasmids are known to host virulence and antibiotic resistance coding genes. The presence of a wide variety of such plasmids in dried milk samples highlights the threat that can occur from consuming contaminated milk powder, and is considered a threat to public health. To the best of our knowledge, mcr-9 has not been isolated from K. pneumoniae in Germany, and further investigations of the occurrence of frameshift mutation are required in the mcr-9 positive strain. Additional studies on a large number of recently isolated strains applying conjugation or transformation assays are required to evaluate the AMR gene transferability.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/microorganisms10030564/s1, Table S1: The average read length and the coverage, data of genome assembly, the results of antibiotic susceptibility testing, metadata, MLST, AMR and virulence-associated genes and plasmid replicons of 24 K. pneumoniae isolates from milk powder samples produced in Germany.

Author Contributions

Conceptualization, research design and data analysis, G.W., H.N., L.D.S. and M.W.P.; drafting and writing of the manuscript, G.W. and L.D.S.; bioinformatic data analysis, J.L.; and sample collection and strain identification P.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by internal funding of the Friedrich-Loeffler-Institut (FLI) and no external fund has been used.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All study data are included in the article and supporting information. The data have also been submitted to the European Nucleotide Archive (ENA). The project accession number is PRJEB45776.

Acknowledgments

We thank Gernot Schmoock, Johannes Solle, Claudia Grosser and Birgit Schikowski for excellent technical assistance. We thank also the Institut Pasteur teams for the curation and maintenance of BIGSdb-Pasteur databases at http://bigsdb.pasteur.fr/, accessed on 28 February 2022.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Navon-Venezia, S.; Kondratyeva, K.; Carattoli, A. Klebsiella pneumoniae: A major worldwide source and shuttle for antibiotic resistance. FEMS Microbiol. Rev. 2017, 41, 252–275. [Google Scholar] [CrossRef] [PubMed]
  2. Santajit, S.; Indrawattana, N. Mechanisms of antimicrobial resistance in ESKAPE pathogens. Biomed. Res. Int. 2016, 2016, 2475067. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Akinyemi, K.O.; Abegunrin, R.O.; Iwalokun, B.A.; Fakorede, C.O.; Makarewicz, O.; Neubauer, H.; Pletz, M.W.; Wareth, G. The emergence of Klebsiella pneumoniae with reduced susceptibility against third-generation cephalosporins and carbapenems in Lagos hospitals, Nigeria. Antibiotics 2021, 10, 142. [Google Scholar] [CrossRef] [PubMed]
  4. Jindal, P.; Bedi, J.; Singh, R.; Aulakh, R.; Gill, J. Phenotypic and genotypic antimicrobial resistance patterns of Escherichia coli and Klebsiella isolated from dairy farm milk, farm slurry, and water in Punjab, India. Environ. Sci. Pollut. Res. Int. 2021, 28, 28556–28570. [Google Scholar] [CrossRef] [PubMed]
  5. Xu, L.; Sun, X.; Ma, X. Systematic review and meta-analysis of mortality of patients infected with carbapenem-resistant Klebsiella pneumoniae. Ann. Clin. Microbiol. Antimicrob. 2017, 16, 18. [Google Scholar] [CrossRef] [Green Version]
  6. Saishu, N.; Ozaki, H.; Murase, T. CTX-M-type extended-spectrum beta-lactamase-producing Klebsiella pneumoniae isolated from cases of bovine mastitis in Japan. J. Vet. Med. Sci. 2014, 76, 1153–1156. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Grohn, Y.T.; Wilson, D.J.; Gonzalez, R.N.; Hertl, J.A.; Schulte, H.; Bennett, G.; Schukken, Y.H. Effect of pathogen-specific clinical mastitis on milk yield in dairy cows. J. Dairy Sci. 2004, 87, 3358–3374. [Google Scholar] [CrossRef]
  8. Hertl, J.A.; Schukken, Y.H.; Welcome, F.L.; Tauer, L.W.; Grohn, Y.T. Pathogen-specific effects on milk yield in repeated clinical mastitis episodes in Holstein dairy cows. J. Dairy Sci. 2014, 97, 1465–1480. [Google Scholar] [CrossRef] [Green Version]
  9. Hu, Y.; Anes, J.; Devineau, S.; Fanning, S. Klebsiella pneumoniae: Prevalence, reservoirs, antimicrobial resistance, pathogenicity, and infection: A hitherto unrecognized zoonotic bacterium. Foodborne Pathog. Dis. 2021, 18, 63–84. [Google Scholar] [CrossRef]
  10. Ritter, E.; Bauernfeind, A.; Becker-Boost, E.; Fiehn, A.; Stocker, H.; Wirsing von Konig, C.H.; Finger, H. Outbreak of a nosocomial infection of SHV2-beta-lactamase-containing Klebsiella pneumonia strains in an operative intensive care unit. Immun. Infekt. 1992, 20, 3–6. [Google Scholar]
  11. Rodloff, A.C.; Dowzicky, M.J. Antimicrobial susceptibility among European Gram-negative and gram-positive isolates collected as part of the tigecycline evaluation and surveillance trial (2004–2014). Chemotherapy 2017, 62, 1–11. [Google Scholar] [CrossRef] [PubMed]
  12. Wareth, G.; Sprague, L.; Neubauer, H.; Pletz, M. Klebsiella pneumoniae in Germany: An overview on spatiotemporal distribution and resistance development in humans. Ger. J. Microbiol. 2021, 1, 16–25. [Google Scholar] [CrossRef]
  13. Zhou, X.; Garcia-Cobos, S.; Ruijs, G.; Kampinga, G.A.; Arends, J.P.; Borst, D.M.; Moller, L.V.; Holman, N.D.; Schuurs, T.A.; Bruijnesteijn van Coppenraet, L.E.; et al. Epidemiology of extended-spectrum beta-lactamase-producing E. coli and vancomycin-resistant Enterococci in the Northern Dutch-German cross-border region. Front. Microbiol. 2017, 8, 1914. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Both, A.; Büttner, H.; Huang, J.; Perbandt, M.; Belmar Campos, C.; Christner, M.; Maurer, F.P.; Kluge, S.; König, C.; Aepfelbacher, M.; et al. Emergence of ceftazidime/avibactam non-susceptibility in an MDR Klebsiella pneumoniae isolate. J. Antimicrob. Chemother. 2017, 72, 2483–2488. [Google Scholar] [CrossRef] [Green Version]
  15. Kresken, M.; Körber-Irrgang, B.; Pfeifer, Y.; Werner, G. Activity of temocillin against CTX-M-producing Escherichia coli and Klebsiella pneumoniae from Germany. Int. J. Antimicrob. Agents 2018, 51, 159–160. [Google Scholar] [CrossRef]
  16. Gürntke, S.; Kohler, C.; Steinmetz, I.; Pfeifer, Y.; Eller, C.; Gastmeier, P.; Schwab, F.; Leistner, R. Molecular epidemiology of extended-spectrum beta-lactamase (ESBL)-positive Klebsiella pneumoniae from bloodstream infections and risk factors for mortality. J. Infect. Chemother. 2014, 20, 817–819. [Google Scholar] [CrossRef]
  17. Klaper, K.; Wendt, S.; Lübbert, C.; Lippmann, N.; Pfeifer, Y.; Werner, G. Hypervirulent Klebsiella pneumoniae of lineage ST66-K2 caused tonsillopharyngitis in a German patient. Microorganisms 2021, 9, 133. [Google Scholar] [CrossRef]
  18. Heiden, S.E.; Hübner, N.O.; Bohnert, J.A.; Heidecke, C.D.; Kramer, A.; Balau, V.; Gierer, W.; Schaefer, S.; Eckmanns, T.; Gatermann, S.; et al. A Klebsiella pneumoniae ST307 outbreak clone from Germany demonstrates features of extensive drug resistance, hypermucoviscosity, and enhanced iron acquisition. Genome Med. 2020, 12, 113. [Google Scholar] [CrossRef]
  19. Wareth, G.; Neubauer, H. The animal-foods-environment interface of Klebsiella pneumoniae in Germany: An observational study on pathogenicity, resistance development, and the current situation. Vet. Res. 2021, 52, 16. [Google Scholar] [CrossRef]
  20. Brisse, S.; Grimont, F.; Grimont, P.A.D. The genus Klebsiella. In The Prokaryotes: Volume 6: Proteobacteria: Gamma Subclass; Dworkin, M., Falkow, S., Rosenberg, E., Schleifer, K.-H., Stackebrandt, E., Eds.; Springer: New York, NY, USA, 2006; pp. 159–196. [Google Scholar] [CrossRef]
  21. Jacobs, C.; Braun, P.; Hammer, P. Reservoir and routes of transmission of Enterobacter sakazakii (Cronobacter spp.) in a milk powder-producing plant. J. Dairy Sci. 2011, 94, 3801–3810. [Google Scholar] [CrossRef]
  22. Khater, D.F.; Lela, R.A.; El-Diasty, M.; Moustafa, S.A.; Wareth, G. Detection of harmful foodborne pathogens in food samples at the points of sale by MALDT-TOF MS in Egypt. BMC Res. Notes 2021, 14, 112. [Google Scholar] [CrossRef] [PubMed]
  23. Wood, D.E.; Lu, J.; Langmead, B. Improved metagenomic analysis with Kraken 2. Genome Biol. 2019, 20, 257. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Wareth, G.; Linde, J.; Hammer, P.; Nguyen, N.H.; Nguyen, T.N.M.; Splettstoesser, W.D.; Makarewicz, O.; Neubauer, H.; Sprague, L.D.; Pletz, M.W. Phenotypic and WGS-derived antimicrobial resistance profiles of clinical and non-clinical Acinetobacter baumannii isolates from Germany and Vietnam. Int. J. Antimicrob. Agents 2020, 56, 106127. [Google Scholar] [CrossRef] [PubMed]
  25. CLSI. Performance Standards for Antimicrobial Susceptibility Testing, 30th ed.; CLSI supplement M100; Clinical and Laboratory Standards Institute: Malvern, PA, USA, 2020. [Google Scholar]
  26. Linde, J.; Homeier-Bachmann, T.; Dangel, A.; Riehm, J.M.; Sundell, D.; Öhrman, C.; Forsman, M.; Tomaso, H. Genotyping of Francisella tularensis subsp. holarctica from Hares in Germany. Microorganisms 2020, 8, 1932. [Google Scholar] [CrossRef]
  27. Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data. 2010. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc (accessed on 28 February 2022).
  28. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [Green Version]
  29. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality assessment tool for genome assemblies. Bioinformatics 2013, 29, 1072–1075. [Google Scholar] [CrossRef]
  30. Chen, L.; Yang, J.; Yu, J.; Yao, Z.; Sun, L.; Shen, Y.; Jin, Q. VFDB: A reference database for bacterial virulence factors. Nucleic Acids Res. 2005, 33, D325–D328. [Google Scholar] [CrossRef] [Green Version]
  31. Jia, B.; Raphenya, A.R.; Alcock, B.; Waglechner, N.; Guo, P.; Tsang, K.K.; Lago, B.A.; Dave, B.M.; Pereira, S.; Sharma, A.N.; et al. CARD 2017: Expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 2017, 45, D566–D573. [Google Scholar] [CrossRef]
  32. Zankari, E.; Hasman, H.; Cosentino, S.; Vestergaard, M.; Rasmussen, S.; Lund, O.; Aarestrup, F.M.; Larsen, M.V. Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother. 2012, 67, 2640–2644. [Google Scholar] [CrossRef]
  33. Feldgarden, M.; Brover, V.; Haft, D.H.; Prasad, A.B.; Slotta, D.J.; Tolstoy, I.; Tyson, G.H.; Zhao, S.; Hsu, C.H.; McDermott, P.F.; et al. Validating the AMRFinder tool and resistance gene database by using antimicrobial resistance genotype-phenotype correlations in a collection of isolates. Antimicrob. Agents Chemother. 2019, 63, e00483-19. [Google Scholar] [CrossRef] [Green Version]
  34. Carattoli, A.; Zankari, E.; García-Fernández, A.; Voldby Larsen, M.; Lund, O.; Villa, L.; Møller Aarestrup, F.; Hasman, H. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob. Agents Chemother. 2014, 58, 3895–3903. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef]
  36. Letunic, I.; Bork, P. Interactive tree of life (iTOL) v4: Recent updates and new developments. Nucleic Acids Res. 2019, 47, W256–W259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Kananizadeh, P.; Oshiro, S.; Watanabe, S.; Iwata, S.; Kuwahara-Arai, K.; Shimojima, M.; Ogawa, M.; Tada, T.; Kirikae, T. Emergence of carbapenem-resistant and colistin-susceptible Enterobacter cloacae complex co-harboring bla(IMP-1) and mcr-9 in Japan. BMC Infect. Dis. 2020, 20, 282. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Gonzalez-Ferrer, S.; Peñaloza, H.F.; Budnick, J.A.; Bain, W.G.; Nordstrom, H.R.; Lee, J.S.; Van Tyne, D. Finding order in the chaos: Outstanding questions in Klebsiella pneumoniae pathogenesis. Infect. Immun. 2021, 89, e00693-20. [Google Scholar] [CrossRef] [PubMed]
  39. Paulin-Curlee, G.G.; Sreevatsan, S.; Singer, R.S.; Isaacson, R.; Reneau, J.; Bey, R.; Foster, D. Molecular subtyping of mastitis-associated Klebsiella pneumoniae isolates shows high levels of diversity within and between dairy herds. J. Dairy Sci. 2008, 91, 554–563. [Google Scholar] [CrossRef] [PubMed]
  40. Muytjens, H.L.; Roelofs-Willemse, H.; Jaspar, G.H. Quality of powdered substitutes for breast milk with regard to members of the family Enterobacteriaceae. J. Clin. Microbiol. 1988, 26, 743–746. [Google Scholar] [CrossRef] [Green Version]
  41. Diaz, D.E.A.M.E.; Leon Duarte, A.B.; Montes De Oca Canastillo, F. Incidence of multiple antibiotic resistant organisms isolated from retail milk products in Hermosillo, Mexico. J. Food Prot. 1992, 55, 370–373. [Google Scholar] [CrossRef]
  42. Diriba, K.; Awulachew, E., 2nd; Tekele, L.; Ashuro, Z. Fecal carriage rate of extended-spectrum beta-lactamase-producing Escherichia coli and Klebsiella pneumoniae among apparently health food handlers in Dilla University student cafeteria. Infect. Drug Resist. 2020, 13, 3791–3800. [Google Scholar] [CrossRef]
  43. Wareth, G.; Brandt, C.; Sprague, L.D.; Neubauer, H.; Pletz, M.W. Spatio-temporal distribution of Acinetobacter baumannii in Germany—A comprehensive systematic review of studies on resistance development in humans (2000–2018). Microorganisms 2020, 8, 375. [Google Scholar] [CrossRef] [Green Version]
  44. Becker, L.; Kaase, M.; Pfeifer, Y.; Fuchs, S.; Reuss, A.; von Laer, A.; Sin, M.A.; Korte-Berwanger, M.; Gatermann, S.; Werner, G. Genome-based analysis of Carbapenemase-producing Klebsiella pneumoniae isolates from German hospital patients, 2008–2014. Antimicrob. Resist. Infect. Control 2018, 7, 62. [Google Scholar] [CrossRef] [PubMed]
  45. Falgenhauer, L.; Schwengers, O.; Schmiedel, J.; Baars, C.; Lambrecht, O.; Heß, S.; Berendonk, T.U.; Falgenhauer, J.; Chakraborty, T.; Imirzalioglu, C. Multidrug-resistant and clinically relevant Gram-negative bacteria are present in German surface waters. Front. Microbiol. 2019, 10, 2779. [Google Scholar] [CrossRef] [PubMed]
  46. Pichler, C.; Büchsel, M.; Rossen, J.W.; Vavra, M.; Reuter, S.; Kern, W.V.; Thimme, R.; Mischnik, A. First report of invasive liver abscess syndrome with endophthalmitis caused by a K2 serotype ST2398 hypervirulent Klebsiella pneumoniae in Germany, 2016. New Microbes New Infect. 2017, 17, 77–80. [Google Scholar] [CrossRef]
  47. Loncaric, I.; Beiglbock, C.; Fessler, A.T.; Posautz, A.; Rosengarten, R.; Walzer, C.; Ehricht, R.; Monecke, S.; Schwarz, S.; Spergser, J.; et al. Characterization of ESBL- and AmpC-producing and fluoroquinolone-resistant Enterobacteriaceae isolated from Mouflons (Ovis orientalis musimon) in Austria and Germany. PLoS ONE 2016, 11, e0155786. [Google Scholar] [CrossRef] [PubMed]
  48. Kaiser, T.; Finstermeier, K.; Hantzsch, M.; Faucheux, S.; Kaase, M.; Eckmanns, T.; Bercker, S.; Kaisers, U.X.; Lippmann, N.; Rodloff, A.C.; et al. Stalking a lethal superbug by whole-genome sequencing and phylogenetics: Influence on unraveling a major hospital outbreak of carbapenem-resistant Klebsiella pneumoniae. Am. J. Infect. Control 2017, 46, 54–59. [Google Scholar] [CrossRef] [PubMed]
  49. Ducomble, T.; Faucheux, S.; Helbig, U.; Kaisers, U.X.; Konig, B.; Knaust, A.; Lubbert, C.; Moller, I.; Rodloff, A.C.; Schweickert, B.; et al. Large hospital outbreak of KPC-2-producing Klebsiella pneumoniae: Investigating mortality and the impact of screening for KPC-2 with polymerase chain reaction. J. Hosp. Infect. 2015, 89, 179–185. [Google Scholar] [CrossRef] [PubMed]
  50. Tietgen, M.; Semmler, T.; Riedel-Christ, S.; Kempf, V.A.J.; Molinaro, A.; Ewers, C.; Gottig, S. Impact of the colistin resistance gene mcr-1 on bacterial fitness. Int. J. Antimicrob. Agents 2017, 51, 554–561. [Google Scholar] [CrossRef]
  51. Hembach, N.; Schmid, F.; Alexander, J.; Hiller, C.; Rogall, E.T.; Schwartz, T. Occurrence of the mcr-1 colistin resistance gene and other clinically relevant antibiotic resistance genes in microbial populations at different municipal wastewater treatment plants in Germany. Front. Microbiol. 2017, 8, 1282. [Google Scholar] [CrossRef]
  52. Bauernfeind, A.; Rosenthal, E.; Eberlein, E.; Holley, M.; Schweighart, S. Spread of Klebsiella pneumoniae producing SHV-5 beta-lactamase among hospitalized patients. Infection 1993, 21, 18–22. [Google Scholar] [CrossRef]
  53. Ewers, C.; Stamm, I.; Pfeifer, Y.; Wieler, L.H.; Kopp, P.A.; Schonning, K.; Prenger-Berninghoff, E.; Scheufen, S.; Stolle, I.; Gunther, S.; et al. Clonal spread of highly successful ST15-CTX-M-15 Klebsiella pneumoniae in companion animals and horses. J. Antimicrob. Chemother. 2014, 69, 2676–2680. [Google Scholar] [CrossRef] [Green Version]
  54. Garcia-Cobos, S.; Kock, R.; Mellmann, A.; Frenzel, J.; Friedrich, A.W.; Rossen, J.W. Molecular typing of Enterobacteriaceae from pig holdings in North-Western Germany reveals extended-spectrum and AmpC beta-lactamases producing but no carbapenem-resistant ones. PLoS ONE 2015, 10, e0134533. [Google Scholar] [CrossRef] [PubMed]
  55. Borgmann, S.; Pfeifer, Y.; Becker, L.; Rieß, B.; Siegmund, R.; Sagel, U. Findings from an outbreak of carbapenem-resistant Klebsiella pneumoniae emphasize the role of antibiotic treatment for cross transmission. Infection 2018, 46, 103–112. [Google Scholar] [CrossRef] [PubMed]
  56. Schmiedel, J.; Falgenhauer, L.; Domann, E.; Bauerfeind, R.; Prenger-Berninghoff, E.; Imirzalioglu, C.; Chakraborty, T. Multiresistant extended-spectrum beta-lactamase-producing Enterobacteriaceae from humans, companion animals and horses in central Hesse, Germany. BMC Microbiol. 2014, 14, 187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Stolle, I.; Prenger-Berninghoff, E.; Stamm, I.; Scheufen, S.; Hassdenteufel, E.; Guenther, S.; Bethe, A.; Pfeifer, Y.; Ewers, C. Emergence of OXA-48 carbapenemase-producing Escherichia coli and Klebsiella pneumoniae in dogs. J. Antimicrob. Chemother. 2013, 68, 2802–2808. [Google Scholar] [CrossRef] [Green Version]
  58. Carroll, L.M.; Gaballa, A.; Guldimann, C.; Sullivan, G.; Henderson, L.O.; Wiedmann, M. Identification of novel mobilized colistin resistance gene mcr-9 in a multidrug-resistant, colistin-susceptible Salmonella enterica serotype typhimurium isolate. mBio 2019, 10, e00853-19. [Google Scholar] [CrossRef] [Green Version]
  59. Ling, Z.; Yin, W.; Shen, Z.; Wang, Y.; Shen, J.; Walsh, T.R. Epidemiology of mobile colistin resistance genes mcr-1 to mcr-9. J. Antimicrob. Chemother. 2020, 75, 3087–3095. [Google Scholar] [CrossRef]
  60. Ito, R.; Mustapha, M.M.; Tomich, A.D.; Callaghan, J.D.; McElheny, C.L.; Mettus, R.T.; Shanks, R.M.Q.; Sluis-Cremer, N.; Doi, Y. Widespread fosfomycin resistance in Gram-negative bacteria attributable to the chromosomal fosA Gene. mBio 2017, 8, e00749-17. [Google Scholar] [CrossRef] [Green Version]
  61. Wong, M.; Wong, D.; Malhotra, S. Intravenous fosfomycin as salvage therapy for osteomyelitis caused by multidrug-resistant Pseudomonas aeruginosa. Am. J. Health Syst. Pharm. 2021, 78, 2209–2215. [Google Scholar] [CrossRef]
  62. Giurazza, R.; Mazza, M.C.; Andini, R.; Sansone, P.; Pace, M.C.; Durante-Mangoni, E. Emerging treatment options for multi-drug-resistant bacterial infections. Life 2021, 11, 519. [Google Scholar] [CrossRef]
  63. Ito, R.; Tomich, A.D.; McElheny, C.L.; Mettus, R.T.; Sluis-Cremer, N.; Doi, Y. Inhibition of fosfomycin resistance protein FosA by phosphonoformate (foscarnet) in multidrug-resistant Gram-negative pathogens. Antimicrob. Agents Chemother. 2017, 61, e01424-17. [Google Scholar] [CrossRef] [Green Version]
  64. Tomich, A.D.; Klontz, E.H.; Deredge, D.; Barnard, J.P.; McElheny, C.L.; Eshbach, M.L.; Weisz, O.A.; Wintrode, P.; Doi, Y.; Sundberg, E.J.; et al. Small-molecule inhibitor of fosA expands fosfomycin activity to multidrug-resistant Gram-negative pathogens. Antimicrob. Agents Chemother. 2019, 63, e01524-18. [Google Scholar] [CrossRef] [Green Version]
  65. Rendón, M.A.; Saldaña, Z.; Erdem, A.L.; Monteiro-Neto, V.; Vázquez, A.; Kaper, J.B.; Puente, J.L.; Girón, J.A. Commensal and pathogenic Escherichia coli use a common pilus adherence factor for epithelial cell colonization. Proc. Natl. Acad. Sci. USA 2007, 104, 10637–10642. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Bachman, M.A.; Lenio, S.; Schmidt, L.; Oyler, J.E.; Weiser, J.N. Interaction of lipocalin 2, transferrin, and siderophores determines the replicative niche of Klebsiella pneumoniae during pneumonia. mBio 2012, 3, e00224-11. [Google Scholar] [CrossRef] [Green Version]
  67. Russo, T.A.; Olson, R.; MacDonald, U.; Beanan, J.; Davidson, B.A. Aerobactin, but not yersiniabactin, salmochelin, or enterobactin, enables the growth/survival of hypervirulent (hypermucoviscous) Klebsiella pneumoniae ex vivo and in vivo. Infect. Immun. 2015, 83, 3325–3333. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Kochan, T.J.; Ozer, E.A.; Pincus, N.B.; Fitzpatrick, M.A.; Hauser, A.R. Complete genome sequence of Klebsiella pneumoniae strain TK421, a conjugative hypervirulent isolate. Microbiol. Resour. Announc. 2020, 9, e01408-19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Pajand, O.; Darabi, N.; Arab, M.; Ghorbani, R.; Bameri, Z.; Ebrahimi, A.; Hojabri, Z. The emergence of the hypervirulent Klebsiella pneumoniae (hvKp) strains among circulating clonal complex 147 (CC147) harbouring bla(NDM/OXA-48) carbapenemases in a tertiary care center of Iran. Ann. Clin. Microbiol. Antimicrob. 2020, 19, 12. [Google Scholar] [CrossRef] [PubMed]
  70. Zhang, R.; Lin, D.; Chan, E.W.; Gu, D.; Chen, G.X.; Chen, S. Emergence of carbapenem-resistant serotype K1 hypervirulent Klebsiella pneumoniae strains in China. Antimicrob. Agents Chemother. 2016, 60, 709–711. [Google Scholar] [CrossRef] [Green Version]
  71. Catalán-Nájera, J.C.; Garza-Ramos, U.; Barrios-Camacho, H. Hypervirulence and hypermucoviscosity: Two different but complementary Klebsiella spp. phenotypes? Virulence 2017, 8, 1111–1123. [Google Scholar] [CrossRef] [Green Version]
  72. Hsu, C.R.; Lin, T.L.; Chen, Y.C.; Chou, H.C.; Wang, J.T. The role of Klebsiella pneumoniae rmpA in capsular polysaccharide synthesis and virulence revisited. Microbiology 2011, 157, 3446–3457. [Google Scholar] [CrossRef]
  73. Russo, T.A.; Gulick, A.M. Aerobactin synthesis proteins as antivirulence targets in hypervirulent Klebsiella pneumoniae. ACS Infect. Dis. 2019, 5, 1052–1054. [Google Scholar] [CrossRef]
  74. Russo, T.A.; Olson, R.; Macdonald, U.; Metzger, D.; Maltese, L.M.; Drake, E.J.; Gulick, A.M. Aerobactin mediates virulence and accounts for increased siderophore production under iron-limiting conditions by hypervirulent (hypermucoviscous) Klebsiella pneumoniae. Infect. Immun. 2014, 82, 2356–2367. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Shankar, C.; Jacob, J.J.; Vasudevan, K.; Biswas, R.; Manesh, A.; Sethuvel, D.P.M.; Varughese, S.; Biswas, I.; Veeraraghavan, B. Emergence of multidrug-resistant hypervirulent ST23 Klebsiella pneumoniae: Multidrug-resistant plasmid acquisition drives evolution. Front. Cell Infect. Microbiol. 2020, 10, 575289. [Google Scholar] [CrossRef] [PubMed]
  76. Liu, C.; Guo, J. Hypervirulent Klebsiella pneumoniae (hypermucoviscous and aerobactin positive) infection over 6 years in the elderly in China: Antimicrobial resistance patterns, molecular epidemiology and risk factor. Ann. Clin. Microbiol. Antimicrob. 2019, 18, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Koczura, R.; Kaznowski, A. Occurrence of the Yersinia high-pathogenicity island and iron uptake systems in clinical isolates of Klebsiella pneumoniae. Microb. Pathog. 2003, 35, 197–202. [Google Scholar] [CrossRef]
  78. Schubert, S.; Cuenca, S.; Fischer, D.; Heesemann, J. High-pathogenicity island of Yersinia pestis in enterobacteriaceae isolated from blood cultures and urine samples: Prevalence and functional expression. J. Infect. Dis. 2000, 182, 1268–1271. [Google Scholar] [CrossRef] [Green Version]
  79. Bardaji, L.; Añorga, M.; Ruiz-Masó, J.A.; Del Solar, G.; Murillo, J. Plasmid replicons from Pseudomonas are natural chimeras of functional, exchangeable modules. Front. Microbiol. 2017, 8, 190. [Google Scholar] [CrossRef]
  80. Dolejska, M.; Brhelova, E.; Dobiasova, H.; Krivdova, J.; Jurankova, J.; Sevcikova, A.; Dubska, L.; Literak, I.; Cizek, A.; Vavrina, M.; et al. Dissemination of IncFII(K)-type plasmids in multiresistant CTX-M-15-producing Enterobacteriaceae isolates from children in hospital paediatric oncology wards. Int. J. Antimicrob. Agents 2012, 40, 510–515. [Google Scholar] [CrossRef]
  81. Ho, W.S.; Yap, K.P.; Yeo, C.C.; Rajasekaram, G.; Thong, K.L. The complete sequence and comparative analysis of a multidrug-resistance and virulence multireplicon IncFII Plasmid pEC302/04 from an extraintestinal pathogenic Escherichia coli EC302/04 indicate extensive diversity of IncFII plasmids. Front. Microbiol. 2015, 6, 1547. [Google Scholar] [CrossRef] [Green Version]
  82. Villa, L.; García-Fernández, A.; Fortini, D.; Carattoli, A. Replicon sequence typing of IncF plasmids carrying virulence and resistance determinants. J. Antimicrob. Chemother. 2010, 65, 2518–2529. [Google Scholar] [CrossRef] [Green Version]
  83. Ramsamy, Y.; Mlisana, K.P.; Allam, M.; Amoako, D.G.; Abia, A.L.K.; Ismail, A.; Singh, R.; Kisten, T.; Han, K.S.; Muckart, D.J.J.; et al. Genomic analysis of carbapenemase-producing extensively drug-resistant Klebsiella pneumoniae isolates reveals the horizontal spread of p18-43_01 plasmid encoding bla(NDM-1) in South Africa. Microorganisms 2020, 8, 137. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Phylogenetic tree using core-genome SNP of K. pneumoniae strains recovered from milk powder samples produced in Germany between 2005–2011.
Figure 1. Phylogenetic tree using core-genome SNP of K. pneumoniae strains recovered from milk powder samples produced in Germany between 2005–2011.
Microorganisms 10 00564 g001
Table 1. Description of the typed and novel STs and allelic profiles of the untyped K. pneumoniae strains recovered from milk powder samples produced in Germany between 2005–2011.
Table 1. Description of the typed and novel STs and allelic profiles of the untyped K. pneumoniae strains recovered from milk powder samples produced in Germany between 2005–2011.
STsIDOriginYeargapAinfBmdhpgiphorpoBtonB
1322DE-MRI1Milk powder200525121818016
DE-MRI12Milk powder200725121818016
DE-MRI14Milk powder200725121818016
567DE-MRI10Milk powder20062177117442
DE-MRI17Milk powder20072177117442
DE-MRI18Milk powder20082177117442
23DE-MRI13Milk powder200721119412
220DE-MRI24Milk powder201121214549
502DE-MRI9Milk powder20032533110418
1083DE-MRI4Milk powder2005212113123
New 1
5625
DE-MRI15Milk powder2007211340122
DE-MRI16Milk powder2007211340122
DE-MRI19Milk powder2008211340122
DE-MRI23Milk powder2010211340122
DE-MRI25Milk powder2011211340122
New2
6016
DE-MRI22Milk powder20102631110225
New3
6014
DE- MRI 2Milk powder20052101393562431
DE- MRI 3Milk powder20052101393562431
DE-MRI 5Milk powder20052101393562431
DE-MRI 6Milk powder20062101393562431
DE-MRI 7Milk powder20062101393562431
DE-MRI 8Milk powder20062101393562431
DE-MRI20Milk powder20092101393562431
DE-MRI21Milk powder20092101393562431
ST: Sequence type.
Table 2. Genetic features of 24 K. pneumoniae isolates originating from milk powder samples.
Table 2. Genetic features of 24 K. pneumoniae isolates originating from milk powder samples.
Parameters of Twenty-Four K. pneumoniae Isolates
Average total number of reads18,554,295 reads per isolate
Isolates mean coverage8558-fold
Genome
Minimum genome size (bp)5.305.557 bp
Maximum genome size (bp)5.770.949 bp
The average of the GC content (%)57.08%
The mean N5027,275 bp
Database accession no.
European Nucleotide Archive (ENA)Project accession number: PRJEB45776
Virulence associated genes (% coverage/No. of isolates)
Enterobactin (Iron uptake)entA (99.20/21); entB (99.18/24); fepC (94.61/24)
Salmochelin (Iron uptake)IroB (97.08/1); iroC (99.28/1); iroD (100/1); iroN (98.81/1). All were found in the same isolate.
Aerobactin (Iron uptake)iucA (99.72/1); iucB (100/1); iucC (99.89/1); iutA (99:96/1). All were found in the same isolate.
Yersiniabactin (Iron uptake)ybtA, ybtE, ybtP, ybtQ, ybtS, ybtT, ybtU, ybtX, irp1, irp2, fyuA. All were found with 100% coverage in the same two isolates.
T6SS-IIompA (100/24).
Adherence (Gene related to ECP from E. coli)ykgK.ecpR (98.31/24); yagV.ecpE (99.74/24); yagW.ecpD (100/24); yagX.ecpC (99.96/24); yagY.ecpB (100/24); yagZ.ecpA (99.32/24).
Regulator of mucoid phenotype ArmpA2 (96.67/1); rmpA (98.57/1); rmpC (90.23/1); and rmpD (100/1).
Antibiotic resistance determinants
Resistance genes (%) blaSHV (100); oqxA (100); oqxB (100); KpnE (100); KpnF (100); KpnG (100); KpnH (100); Ompk37 (100); acrA (100); acrB (100); acrD (100); emrR (100); emrD (100); ramA (100); fosA (71); pmrB_R256G (4); mcr-9 (4).
Plasmid repliconsCol family (Col440II; Col440I; ColRNAI); IncF family (IncFIB.K.; IncFIB.Mar.; IncFIB.pKPHS1.; IncFII_1_pKP91); Inc family (IncHI1B and IncR); Rep A family (RepA_1).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wareth, G.; Linde, J.; Hammer, P.; Pletz, M.W.; Neubauer, H.; Sprague, L.D. WGS-Based Phenotyping and Molecular Characterization of the Resistome, Virulome and Plasmid Replicons in Klebsiella pneumoniae Isolates from Powdered Milk Produced in Germany. Microorganisms 2022, 10, 564. https://doi.org/10.3390/microorganisms10030564

AMA Style

Wareth G, Linde J, Hammer P, Pletz MW, Neubauer H, Sprague LD. WGS-Based Phenotyping and Molecular Characterization of the Resistome, Virulome and Plasmid Replicons in Klebsiella pneumoniae Isolates from Powdered Milk Produced in Germany. Microorganisms. 2022; 10(3):564. https://doi.org/10.3390/microorganisms10030564

Chicago/Turabian Style

Wareth, Gamal, Jörg Linde, Philipp Hammer, Mathias W. Pletz, Heinrich Neubauer, and Lisa D. Sprague. 2022. "WGS-Based Phenotyping and Molecular Characterization of the Resistome, Virulome and Plasmid Replicons in Klebsiella pneumoniae Isolates from Powdered Milk Produced in Germany" Microorganisms 10, no. 3: 564. https://doi.org/10.3390/microorganisms10030564

APA Style

Wareth, G., Linde, J., Hammer, P., Pletz, M. W., Neubauer, H., & Sprague, L. D. (2022). WGS-Based Phenotyping and Molecular Characterization of the Resistome, Virulome and Plasmid Replicons in Klebsiella pneumoniae Isolates from Powdered Milk Produced in Germany. Microorganisms, 10(3), 564. https://doi.org/10.3390/microorganisms10030564

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop