Genomic and functional characterization of carbapenem-resistant Klebsiella pneumoniae from hospital wastewater
BMC Microbiology volume 23, Article number: 115 (2023)
The emergence of carbapenem-resistant Klebsiella pneumoniae (CRKP) attracted extensive attention. Information on CRKP from hospital wastewater (HWW) is limited. The aims of this study were to investigate the genomic characteristics and to evaluate the survivability characteristics of 11 CRKP from HWW in a Chinese teaching hospital in Fujian province.
A total of 11 CRKP from HWW were recovered in this study. All CRKP from HWW were resistant to most antibiotics. Comparative genetic analysis demonstrated that all CRKP isolates were clustered into the three distinct phylogenetic clades and clade 2 and clade 3 were mixtures of samples collected from both HWW and clinical settings. Varieties of resistance genes, virulence genes and plasmid replicon types were detected in CRKP from HWW. In vitro transfer of blaKPC-2 was successful for 3 blaKPC-2-positive CRKP from HWW with high conjugation frequency. Our study demonstrated that the genetic environments of blaKPC−2 shared core structure with ISKpn27-blaKPC−2-ISKpn6. Group analysis showed that CRKP from HWW had a lower survivability in serum compared to clinical CRKP (p < 005); and CRKP from HWW had no significant difference in survivability in HWW compared to clinical CRKP (p > 005).
We analyzed the genomic and survivability characteristics of CRKP from HWW in a Chinese teaching hospital. These genomes represent a significant addition of genomic data from the genus and could serve as a valuable resource for future genomic studies about CRKP from HWW.
The emergence of carbapenem-resistant Klebsiella pneumoniae (CRKP) attracted extensive attention, which was a significant public health challenge worldwide . CRKP can cause many kinds of infections, such as pneumonia, bloodstream infections and urinary tract infections, and responsible for substantial patient morbidity and mortality .
Notably, many CRKP isolates were recovered from the environments, especially the water environment, such as the wastewater, rivers and lakes [3,4,5]. Hospital wastewater (HWW) was recognized as a putative reservoir for CRKP . Current evidence concurs that HWW was an important source for antibiotic resistance in aquatic environments, mainly multidrug-resistant Gram-negative bacteria. Compared to other wastewaters, antibiotic-resistant Gram-negative bacteria were more prevalent in HWW, with higher concentrations of extended spectrum beta-lactamase (ESBLs)-producing pathogens and carbapenemase-producing Enterobacterales in hospital sources .
There are substantial studies regarding the global prevalence and trends of clinical CRKP . Many previous studies of HWW mainly focused on the prevalence of carbapenem-resistant genes, molecular characteristic of bacteriophage in CRKP and CRKP infection or colonization associated with intensive care unit sewage in China [3, 4, 8, 9]. Although there was increased research on CRKP from HWW epidemiology and phylogenomic worldwide , large gaps remain in our understanding regarding the genomic and functional characterization of CRKP from HWW, especially in China. Therefore, the aims of this study were to investigate the genomic characteristics using whole genome sequencing and to evaluate the survivability characteristics of CRKP from HWW, which could help us to obtain a comprehensive understanding of CRKP from HWW.
Bacterial isolates and antimicrobial susceptibility testing
In this study, hospital wastewaters contain CRKP approximately at 1.39 × 104 CFU/ml using serial diluted and a total of 11 CRKP isolates were recovered. All CRKP from HWW were resistant or intermediate to most antibiotics, including ertapenem (100%, MIC values ranged from 8 to 128 μg/mL), imipenem (100%, MIC values ranged from 8 to 64 μg/mL), meropenem (100%, MIC values ranged from 8 to 128 μg/mL) and tigecycline (100%, MIC values ranged from 4 to 16 μg/mL). 91.8% of CRKP isolates (from HWW) were interpreted as intermediate (≤ 2 μg/mL, MIC values ranged from 1 to 64 μg/mL) for colistin (Table 1).
MLST, virulence genes and plasmid replicon types
Figure 1 showed the general genome features of CRKP from HWW and clinical CRKP sequenced genomes in this study. Genetic analysis demonstrated that 11 CRKP from HWW belonged to five sequence types (STs), including ST11, ST15, ST824, ST2433 and ST3184. 11 clinical CRKP belonged to five sequence types, including ST11, ST34, ST37, ST437and ST789.
The distributions of virulence genes (VFs) among all CRKP from HWW were shown in Fig. 1. At least one VF was detected in each isolate. Among them, iutA (11/11, 100%) was the most prevalent VF in this study, following by ccI (4/11, 36.36%), fyuA (4/11, 36.36%), irp2 (4/11, 36.36%), terC (3/11, 27.27%) and traT (3/11, 27.27%). Additionally, iucC was only detected in clinical CRKP isolates.
Overall, we found that all CRKP from HWW contained at least one plasmid replicon, covering seven different known plasmid incompatibility types in this study. F-type plasmid was highly prevalent across CRKP from HWW, including IncFIA, IncFIB and IncFII. In addition, Col440I, IncFII(Yp), IncHI1B(pNDM-MAR), IncX1, IncX3 and repB were only identified in clinical CRKP, and IncHI2A was only detected in CRKP from HWW, respectively.
Determination of resistance genes and porin-associated genes
The distributions of antimicrobial resistance genes (ARGs) differed among CRKP from HWW (Fig. 1). Some ARGs were identified more than 70% among CRKP from HWW isolates, such as fosA (fosfomycin resistance gene, 100%), tetA (tetracycline resistance gene, 100%), oqxA and oqxB (fluoroquinolone resistance genes, 100%).
Only four CRKP from HWW isolates (W25A96, W25A97, W25A158 and W25A58) were detected to carry carbapenems resistance gene, blaKPC-2. The results of carbapenems resistance gene were verified by PCR. Other seven strains carried at least one non-carbapenemase β-lactamase gene targeted for detection in this study. Resistance genes found among the seven isolates without carbapenemase resistance gene were blaDHA, blaTEM, blaSHV, blaCTX-M, blaLAP and blaOXA genes (Table 2). In total, deletion of porin-associated genes was detected in 100% of the CRKP from HWW isolates without carbapenemase resistance gene (Fig. 1 and Tables 2).
Notably, we found two types of carbapenems resistance gene (seven isolates carrying blaKPC-2 and one isolates carrying blaNDM-5) in clinical CRKP isolates (8/11, 72.7%), which was a significant difference between CRKP from HWW and clinical settings (p < 0.05).
A total of 22 recoverable CRKP isolates were whole genome sequenced. The core-genome alignments were applied for phylogenetic tree reconstruction using maximum likelihood estimation (Fig. 1). All the 22 CRKP isolates were clustered into the three distinct phylogenetic clades (clade 1 to 3) in this study. Notably, clade 2 and clade 3 were mixtures of samples collected from both HWW and clinical settings. Among those clades, clade 1 included two isolates clinical CRKP isolates (belonged to ST37 and ST437), clade 2 included 10 CRKP isolates from HWW and three clinical CRKP isolates (belonged to ST15, ST34, ST789, ST824, ST2433 and ST3184), and clade 3 included six clinical CRKP isolates and one CRKP isolate from HWW (all isolates belonged to ST11). Furthermore, same sequence type was clustered in a subbranch relatively in clade 2.
Conjugation experiments and plasmid replicon type analysis
In this study, three of four blaKPC-2-positive isolates (W25A96, W25A97 and W25A158) could transfer blaKPC-2 into EC600 successfully. Among them, the conjugation frequency of W25A96 was 1.58 × 10–4, W25A97 was 1.99 × 10–4, and W25A158 was 4.62 × 10–4. The presence of blaKPC-2 and plasmid replicon types was confirmed by PCR. Only IncFII pHN7A8 was found in those transconjugants. The results of antimicrobial susceptibility testing among recipient cell and transconjugants were shown in Table 1.
Genetic Environment of the bla KPC-2 Gene
Horizontal gene transfer has been a major reason for rapid evolution in the resistance of bacteria, which leads to the capture of resistance genes . The rapid spread of blaKPC-2 between different bacteria is due to the propagation of the blaKPC-2-positive plasmid with transposons and sequence elements (IS).
Here, we further detected the genetic environment of blaKPC-2. Four blaKPC-2-positive CRKP isolates from HWW shared same core structure with ISKpn27-blaKPC−2-ISKpn6 in this study (Fig. 2). A ISKpn27 element was located in upstream of blaKPC-2. Furthermore, a ISKpn6 element was located downstream of blaKPC-2, followed by a hypothetical protein, an antirestriction protein, a hypothetical protein and a replication protein. Furthermore, the genetic environments of blaKPC-2 of CRKP isolates from HWW showed highly similar structure compared to blaKPC-2-positive clinical CRKP strains.
Survivability-associated phenotypic characteristics
Serum resistance assay
The results of serum resistance assay were shown in Fig. 3. All CRKP were tested in 50% human serum to investigate their capacity to resist the serum bactericidal activity. CRKP from HWW and clinical CRKP isolates showed significant growth compared to negative control, E. coli DH5α (p < 0.05). Group analysis revealed that clinical CRKP isolates had a stronger survivability compared to CRKP from HWW isolates (p < 005) in our study.
The results of survival assay were presented in Fig. 4. 72.7% (8/11) of CRKP from HWW and all clinical CRKP survived in HWW for 50 days of monitoring. 81.8% 9/11) of CRKP from HWW and all clinical CRKP had an increase within seven days of incubation. 72.7% (8/11) of CRKP from HWW and 72.7% (8/11) of clinical CRKP showed a slight decrease after 14 days. After 50 days of incubation, 54.5% (6/11) of CRKP from HWW and all of clinical CRKP remained survive more than 50% of their original colony forming units. Group analysis revealed that CRKP from HWW had no significant difference in survival capability in autoclaved HWW compared to clinical CRKP (p > 005) within 21 days in our study. Clinical CRKP showed significant difference in survival capability in 28 days, 35 days, 42 days and 50 days compared to CRKP from HWW (p < 0.05, Fig. 4).
CRKP posed a severe clinical problem given the lack of therapeutic options available . Nowadays, CRKP was detected in many environments, such as clinical settings, hospital wastewater, municipal wastewater and so on . HWW was recognized as a putative reservoir for CRKP due to its selective pressure. Nowadays, there were many reports of CRKP from HWW around the world, such as Germany , United States , China , Romania , South Africa  and so on. Those findings indicated their possible spread in the kinds of water environments and the potential risk of human, livestock and wildlife colonization and/or infection associated with exposure to contaminated water sources. However, there was little research on the genomic diversification and functional characteristics of CRKP from HWW, especially in China . Here, we reported the whole genome sequences and comparative functional analysis of CRKP from HWW in a Chinese teaching hospital. These genomes represent a significant addition of genomic data from the genus and could serve as a valuable resource for future genomic studies about CRKP from HWW.
In the present study, 1.39 × 104 CFU/ml of CRKP were detected from HWW approximately and a total of 11 strains were recovered. The detection rate was higher than that in South Africa and lower than that in urban wastewater in eastern India [6, 15]. Antibiotic selection pressure could affect the fraction of bacterial genera potentially . So, this discrepancy might be related to the different antibiotic selection pressure in different regions, which deserved to investigate further.
Molecular typing is usually used to identify clonal relationships among bacteria and applied to detect the population structure of CRKP strains with respect to their specificity for infecting humans or animals . Previous studies reported that ST11-CRKP was disseminated widely in clinical settings across China [18, 19]. Meanwhile, ST11 was identified as common STs in this hospital (Fig. 1) and Fujian province as described previously [20, 21]. In this study, our results revealed that CRKP from HWW possessed high diversity at the level of sequence types (Fig. 1). As similar to previous studies, there is diversity of ST types of CRKP from HWW and their distribution and dominant groups are different from those in the clinical environment [6, 22]. Interestingly, our findings of comparative genomic analysis showed that clinical CRKP and CRKP from HWW have an intersection in same clade (clade 2 and clade 3) and the blaKPC-2 were all found in those STs (Fig. 1). Those similar circulating STs might suggest the dissemination of this sequence type from hospital settings to aquatic environments and the possibility of spread clonally of blaKPC-2-carrying CRKP strains. Similar to previous report, the ST of non-carbapenemase-producing CRKP (NC-CRKP) also showed diversity, which distinguished from the dominant STs of CRKP carrying carbapenemase genes . In addition, the phylogenetic tree showed that CRKP from HWW possessed high diversity at the level of phylogeny that associated with sequence types (Fig. 1), which was similar to previous study . Combined with the results of profiles of resistance genes, virulence genes and plasmid replicon types, these findings suggested that the most CRKP from HWW and clinical CRKP share variant genetic backgrounds with different STs and might belong to the distance ancestral lineage(s).
Carbapenems resistance was mainly linked to the acquisition of carbapenems resistance genes (such as blaKPC, blaNDM, blaVIM), efflux pumps genes and mutations in outer membrane porins [24, 25]. In this study, CRKP from HWW exhibited a diversity of antibiotics resistance gene profiles, particularly including carbapenems resistance genes in four isolates (Fig. 1). But most of CRKP from HWW (7/11, 63.6%) in current study did not carry any carbapenems resistance genes (Fig. 1). This is an interesting phenomenon. NC-CRKP confers carbapenem resistance association with chromosomal mutations (such as porin gene mutation, and overproduction of efflux pump) and acquired non-carbapenemase resistance mechanisms (acquisition or upregulation of a β-lactamase such as ESBL or AmpC β-lactamase) [26, 27]. Meanwhile, loss or deficiency of outer membrane porins were associated with carbapenem resistance . In our study, some ESBLs genes (such as blaSHV), AmpC β-lactamase genes (such as blaDHA-1) and outer membrane porins genes could be detected in those NC-CRKP strains from HWW (Fig. 1 and Table 2). Therefore, future studies are needed to verify whether it’s a feature of CRKP in HWW by expanding the sample size and warranted to confirm the contribution of these factors to carbapenem resistance for NC-CRKP from HWW. . Additionally, 90.9% CRKP from HWW had resistance to tigecycline, which was not consistent with the previous report . Tigecycline resistance has inevitably emerged over the recent years, mostly identified among extensively drug-resistant and carbapenem-resistant isolates . However, any gene associated with tigecycline resistance has not been found in our study. There may be other undiscovered mechanisms of tigecycline resistance, which are worth further exploration. Interestingly, the transconjugants carried blaKPC-2 were all sensitive to cefepime, which was similar to previous studies [32, 33]. The potential mechanism of cefepime-resistance change through conjugation is worth further study.
The dissemination of carbapenems resistance is associated with horizontal transfer of carbapenems resistance genes . In this study, our study demonstrated that most of blaKPC-2-harboring plasmid in CRKP from HWW (3/4) was transferrable with high conjugation frequency. Given the high transferability of blaKPC-2-harboring-IncFII pHN7A8 plasmid in CRKP from HWW, the spread of blaKPC-2-harboring plasmids into other bacteria had the potential to cause the spread of carbapenems resistance in HWW. Meanwhile, many studies characterized the genetic environments of blaKPC-2 and various mobile genetic elements played a critical role in the rapid spread of blaKPC-2 [35, 36]. In this study, we found that the genetic environments of blaKPC−2 shared the core structure with ISKpn27-blaKPC−2-ISKpn6, which was similar to the previous study . Bacteria could acquire preexisting resistance determinants to promote drug resistance, which was achieved through the concerted activities of mobile genetic elements (such as IS, transposons) and plasmids and integrative conjugative elements . The main genetic structure containing transposons and IS could enhance the spread of the blaKPC-type genes in different plasmid scaffolds . Thus, our results indicated that these mobile elements played a key role in the dissemination of blaKPC−2 in CRKP from HWW. Current study indicated the semblable genetic environments of blaKPC-2 between CRKP from HWW and clinical settings, and our previous study also found similar plasmids carrying blaKPC−2 in this hospital in clinical CRKP , which was also evidence that the possibility of the dissemination of CRKP from hospital settings to aquatic environments.
To gain understanding of CRKP from HWW survivability, the 11 CRKP from HWW along with 11 clinical CRKP were assessed for their survivability to serum and hospital wastewater in vitro. CRKP from HWW had a weaker survivability in serum and similar survivability in HWW compared to clinical CRKP in this study (Fig. 3 and Fig. 4). Serum resistance was one of the major survival mechanisms of bacteria enabling them to survive in the bloodstream of the host . Our study indicated that clinical CRKP might more potentially contribute to the bloodstream infections than CRKP from HWW. Additionally, previous study suggested some abiotic and biotic ecological factors were shown to determine CRKP survival and multiplication when CRKP isolates came in the environment . Our results indicated that CRKP from HWW and clinical CRKP might have the potential to persist in HWW for a long time and CRKP from HWW was a potential risk of carbapenems resistance dissemination. Meanwhile, more research to identify mitigating measures is needed to effectively track the dissemination of these bacteria from different environmental resistomes to humans.
There was a limitation to our study. The CRKP from HWW were isolated from a single hospital with small sample size. Further study needs to increase the sample size, including new CRKP from HWW from other hospital in this area and other provinces in China to make the results of comparative genomics analysis more representatives in the future.
In this study, we analyzed the molecular characteristic of CRKP isolates from HWW in a Chinese teaching hospital using whole genome sequencing. High diversity at the level of sequence types and phylogeny were found. Only one type of carbapenems resistance genes – blaKPC-2 was detected and most of them located on a transferrable plasmid with high conjugation frequency. All of blaKPC−2 in this study shared same genetic environment’s structure. Moreover, in our study, CRKP from HWW showed lower survivability in serum and similar survivability in hospital wastewater compared to clinical CRKP. These findings could add to the data on the genomic and functional characterization of CRKP from HWW in China and provide a good pilotage for further in-depth studies of CRKP’s from HWW properties and careful surveillance of its emergence across different ecologies.
The hospital wastewater was collected from the raw sewage influent of on-site hospital wastewater treatment from Fujian medical university union hospital (Fujian province, China), and the collection points were only at locations specific to this hospital . Briefly, a total of three hospital wastewater samples were taken between November 3 and November 25, each ten days apart. The samples were collected and then transported to laboratory on ice for immediate processing during half-hour. After rinsing bottles twice with target samples, 1 L wastewater samples were filled into sterile plastic bottles. After being thoroughly mixed, wastewater samples were serially diluted. Then 1 mL serial diluted wastewater samples were plated on Luria–Bertani (LB) agar (Sangon Biotech (Shanghai) Co., Ltd, China) (supplemented with Imipenem at 4 µg/mL) to select carbapenem resistance bacteria followed by incubation at 37 ℃ for 24 h. Wastewater samples of each dilution were plated on LB agar in triplicates.
The recovered bacterial species were identified using MALDI-TOF (Bruker Daltonics, Germany). Briefly, every single clone was grown on Columbia blood agar (Autobio, China) and incubated for 18 h at a temperature of 36 °C. Then, fresh bacterial colonies were coated on the MALDI target plates (Bruker Daltonik, Germany) and covered with 1 μL chemical matrix (α-cyano-4-hydroxycinnamic acid). Collection of spectra (2–20 kDa) was performed in linear and positive mode of each point by FlexControl software (Bruker Daltonik). If a colony was unidentified, 16S rRNA sequencing was performed . The primers of 16S rRNA were listed in Table S3 (Supplementary Material). The strains identified to Klebsiella pneumoniae using MALDI-TOF and 16 s rRNA sequencing were collected and tested for carbapenem susceptibility performed by broth dilution method. The isolate, whose MIC value of imipenem was 4 μg/mL or higher, was recorded as a positive confirmation result for CRKP.
For comparative analysis, 11 clinical CRKP isolates were collected from Fujian Medical University Union Hospital between November to December, 2019. The clinical CRKP strains were then randomly assigned in a ratio into matched arm using random number generators (SPSS, version.20) to avoid the selection bias. The microbiological characteristics and antibiotics profiles of 11 clinical CRKP isolates were listed in Table S1 and Table S2 (Supplementary Material).
Antimicrobial susceptibility testing
The antimicrobial susceptibility testing of the CRKP from HWW, recipient and transconjugants was performed by Vitek-2 (Vitek-AST-GN16) (BioMerieux, France) according to the Clinical and Laboratory Standards Institute guidelines (CLSI, 2022). The antibiotics included ampicillin, amoxicillin/clavulanic acid, piperacillin/tazobactam, cefazolin, cefoxitin, ceftriaxone, cefepime, aztreonam, amikacin, gentamicin, tobramycin, ciprofloxacin, levofloxacin, nitrofurantoin, and trimethoprim/sulfamethoxazole. Among them, MICs of ertapenem (Macklin, China), imipenem (Solarbio, China), meropenem (Solarbio, China), tigecycline (Macklin, China) and colistin (Solarbio, China) were determined using broth dilution method. E. coli ATCC 25,922 was used as quality control.
Whole genome sequencing and analysis
The genomic DNA samples were extracted from the CRKP from HWW and clinical CRKP isolates by TIANamp Bacteria DNA Kit (Tiangen, Beijing), sequenced using the Illumina NovaSeq platform in Shanghai Personal Biotechnology Co., Ltd (Shanghai, China), with a paired-end mode of 2 × 150 bp. Raw data was checked for quality with FastQC (version.0.11.7). AdapterRemoval (version.2.2.2) was used for trimming of adapter sequences and low-quality bases and SOAPec (version.2.03) was used for quality correction based on K-mer frequency before data assembly. The filtered raw data were assembled to contigs and scaffolds using SPAdes (version.3.12.0) and A5-MiSeq (version.20160825). The base correction was performed using Pilon (version.1.18). The depth of coverage was ranged from 224-fold for 10,459,852 reads to 292-fold for 10,256,762 reads and the average depth of coverage was 260-fold.
Gene content analysis was performed using tools available through the Center for Genomic Epidemiology (https://cge.cbs.dtu.dk/services/). Specifically, acquired resistance genes and known chromosomal mutations conferring antibiotic resistance were identified using ResFinder . Plasmid types were identified using PlasmidFinder . Virulence genes were identified using VirulenceFinder . The multilocus-sequence typing was identified on MLST (https://cge.food.dtu.dk/services/MLST/).
The genetic environment of blaKPC-2 was predicted in RAST tool (https://rast.nmpdr.org/rast.cgi). Comparative analysis of the genetic environment of blaKPC-2 in CRKP from HWW and clinical CRKP was performed using Easyfig (version.2.1).
Detection of carbapenemase genes
To verify the prediction results of carbapenemase genes using whole genome sequencing, the presence of carbapenemase genes was identified by PCR as described previously , including blaKPC, blaNDM, blaVIM, blaIMP, blaOXA-48 and blaGES. The primers of carbapenemase genes are listed in Table S3 (Supplementary Material).
Core genomes were defined as described previously . Based on the comparative analysis of 11 CRKP from HWW and 11 clinical CRKP isolates, core genomes were selected as single-copy core-genome for sequence alignment using mafft (version.7.429) . The set of single-copy core-genome clusters detected was entered into the pipeline to identify high-quality phylogenetic markers and infer a core genome phylogeny through a maximum-likelihood tree search. Maximum likelihood phylogeny of core of 22 CRKP isolates were inferred from the alignment using FastTree (version.2.1.11) with Whelan-and-Goldman 2001 model and 1000 bootstrap replicates, respectively . The resulting phylogenetic tree was visualized using iTOL.
Conjugation experiments and plasmid replicon types analysis
Conjugation experiments were performed to test the transferability of blaKPC-2-harboring plasmids using filter mating . Rifampin-resistant E. coli EC600 (stored in our laboratory) was used as the recipient. The mating mixture was washed from the filter and spread onto MH agar (Sangon Biotech (Shanghai) Co., Ltd, China) containing sodium rifampicin at 600 μg/mL and imipenem at 1 μg/mL. The transconjugants harbored blaKPC-2 were selected and identified by PCR . blaKPC-2 was amplified using primers listed in Table S3 (Supplementary Material). The conjugation frequency (CF) was calculated as follow :
The transconjugants harbored blaKPC-2 were examined for the presence of three plasmid replicons, combining the results of gene prediction, using the PCR as mentioned above . Three plasmid replicons were amplified using primers listed in Table S3 (Supplementary Material).
Survivability-associated phenotypic characteristics
Serum resistance assay
Serum resistance assay of 11 CRKP from HWW and 11 clinical CRKP was tested using 50% human serum (Precision BioMedicals Co., Ltd, China) for a 3-h incubation as previous studies [41, 51]. Briefly, 5 μL of an overnight bacterial culture was added to 495 μL of LB broth and centrifuged at 9000 rpm for 3 min. Then the pellet obtained was re-suspended in 1 × phosphate-buffered saline. 30 μL of the above bacterial inoculum and 270 μL of the 50% human serum were added and mixed in a 96-well plate for a 3-h incobation. Then the mixture was taken for serial dilutions, plated on LB agar and cultured overnight. The growth in the serum of strain was detected as the total number of CFU/mL recovered per well . The serum resistance assay was performed twice in triplicates. E. coli EC505 (stored in our laboratory) was used as the positive control and E. coli DH5α was used as the negative control [41, 52].
The survival assay of CRKP from HWW and clinical CRKP isolates was monitored for 50 days at room temperature (22 ± 2 ℃). HWW was filtered using membrane filter with 0.25 µm size. HWW was partitioned per 200 mL into bottles and autoclaved (121 °C for 20 min). 1 mL overnight bacterial culture (in LB broth) was suspended in duplicate in test tubes containing 40 mL of the autoclaved HWWs. Tubes were rotated at 180 rpm using a StuartTube Rotator THZ-100 (YiHeng, Shanghai, China). The number of bacteria was determined on LB agar in duplicate, which incubated at 35 ℃ for 24 h at 1-week intervals, and bacterial concentration was expressed as log10 CFU/mL . For group comparison, the number of every strain was counted by growth proportion in 7, 14, 21, 28, 35, 42 and 50 days compared to the original colony forming units (0 day) .
Graphpad prism v.7 was used for statistical analysis. Unpaired two-tailed Student’s t-tests and Wilcoxon sign rank test (unequal variances) were performed to analyze statistical significance for serum resistance assay and survival assay. Only p < 0.05 was considered statistically significant.
Availability of data and materials
The genome sequences of 11 CRKP from HWW and 11 clinical CRKP were deposited in GenBank under BioProject PRJNA854359 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA854359).
Carbapenem-resistant Klebsiella pneumoniae
Extended spectrum beta-lactamase
Multilocus sequence typing
Minimum inhibitory concentration
Yang X, Dong N, Chan EW, Zhang R, Chen S. Carbapenem resistance-encoding and virulence-encoding conjugative plasmids in Klebsiella pneumoniae. Trends Microbiol. 2021;29(1):65–83.
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(1):18.
Li M, Guo M, Chen L, Zhu C, Xiao Y, Li P, Guo H, Chen L, Zhang W, Du H. Isolation and characterization of novel lytic bacteriophages infecting epidemic carbapenem-resistant Klebsiella pneumoniae strains. Front Microbiol. 2020;11:1554.
Feng Y, Wei L, Zhu S, Qiao F, Zhang X, Kang Y, Cai L, Kang M, McNally A, Zong Z. Handwashing sinks as the source of transmission of ST16 carbapenem-resistant Klebsiella pneumoniae, an international high-risk clone, in an intensive care unit. J Hosp Infect. 2020;104(4):492–6.
Runcharoen C, Moradigaravand D, Blane B, Paksanont S, Thammachote J, Anun S, Parkhill J, Chantratita N, Peacock SJ. Whole genome sequencing reveals high-resolution epidemiological links between clinical and environmental Klebsiella pneumoniae. Genome medicine. 2017;9(1):6.
Ekwanzala MD, Dewar JB, Kamika I, Momba MNB. Tracking the environmental dissemination of carbapenem-resistant Klebsiella pneumoniae using whole genome sequencing. Sci Total Environ. 2019;691:80–92.
Hassoun-Kheir N, Stabholz Y, Kreft JU, de la Cruz R, Romalde JL, Nesme J, Sørensen SJ, Smets BF, Graham D, Paul M. Comparison of antibiotic-resistant bacteria and antibiotic resistance genes abundance in hospital and community wastewater: a systematic review. Sci Total Environ. 2020;743: 140804.
Gao M, Wang C, Qiang X, Liu H, Li P, Pei G, Zhang X, Mi Z, Huang Y, Tong Y, et al. Isolation and characterization of a novel bacteriophage infecting carbapenem-Resistant Klebsiella pneumoniae. Curr Microbiol. 2020;77(5):722–9.
Mills MC, Lee J. The threat of carbapenem-resistant bacteria in the environment: evidence of widespread contamination of reservoirs at a global scale. Environ Pollut. 2019;255(Pt 1): 113143.
Lopatkin AJ, Sysoeva TA, You L. Dissecting the effects of antibiotics on horizontal gene transfer: analysis suggests a critical role of selection dynamics. BioEssays. 2016;38(12):1283–92.
Savin M, Bierbaum G, Mutters NT, Schmithausen RM, Kreyenschmidt J, García-Meniño I, Schmoger S, Käsbohrer A, Hammerl JA. Genetic characterization of carbapenem-resistant Klebsiella spp. from municipal and slaughterhouse wastewater. Antibiotics. 2022;11(4):435.
Loudermilk EM, Kotay SM, Barry KE, Parikh HI, Colosi LM, Mathers AJ. Tracking Klebsiella pneumoniae carbapenemase gene as an indicator of antimicrobial resistance dissemination from a hospital to surface water via a municipal wastewater treatment plant. Water Res. 2022;213: 118151.
Zhang L, Ma X, Luo L, Hu N, Duan J, Tang Z, Zhong R, Li Y. The prevalence and characterization of extended-spectrum β-lactamase- and carbapenemase-producing bacteria from hospital sewage, treated effluents and receiving rivers. Int J Environ Res Public Health. 2020;17(4):1183.
Surleac M, Czobor Barbu I, Paraschiv S, Popa LI, Gheorghe I, Marutescu L, Popa M, Sarbu I, Talapan D, Nita M, et al. Whole genome sequencing snapshot of multi-drug resistant Klebsiella pneumoniae strains from hospitals and receiving wastewater treatment plants in Southern Romania. PLoS One. 2020;15(1): e0228079.
Sahoo S, Sahoo RK, Gaur M, Behera DU, Sahu A, Das A, Dey S, Dixit S, Subudhi E. Urban wastewater contributes to the emergence of carbapenem-resistant Klebsiella pneumoniae (CRKP) in an urban receiving river in eastern India. Lett Appl Microbiol. 2023;76(1):ovac005.
Tello A, Austin B, Telfer TC. Selective pressure of antibiotic pollution on bacteria of importance to public health. Environ Health Perspect. 2012;120(8):1100–6.
Qin X, Wu S, Hao M, Zhu J, Ding B, Yang Y, Xu X, Wang M, Yang F, Hu F. The colonization of carbapenem-resistant Klebsiella pneumoniae: Epidemiology, resistance mechanisms, and risk factors in patients admitted to intensive care units in China. J Infect Dis. 2020;221(Suppl 2):S206-s214.
Liao W, Liu Y, Zhang W. Virulence evolution, molecular mechanisms of resistance and prevalence of ST11 carbapenem-resistant Klebsiella pneumoniae in China: A review over the last 10 years. J Global Antimicrob Resist. 2020;23:174–80.
Jin X, Chen Q, Shen F, Jiang Y, Wu X, Hua X, Fu Y, Yu Y. Resistance evolution of hypervirulent carbapenem-resistant Klebsiella pneumoniae ST11 during treatment with tigecycline and polymyxin. Emerg Microbes Infect. 2021;10(1):1129–36.
Li B, Xu XH, Zhao ZC, Wang MH, Cao YP. High prevalence of metallo-β-lactamase among carbapenem-resistant Klebsiella pneumoniae in a teaching hospital in China. Can J Microbiol. 2014;60(10):691–5.
Yu X, Zhang W, Zhao Z, Ye C, Zhou S, Wu S, Han L, Han Z, Ye H. Molecular characterization of carbapenem-resistant Klebsiella pneumoniae isolates with focus on antimicrobial resistance. BMC Genomics. 2019;20(1):822.
Tafoukt R, Touati A, Leangapichart T, Bakour S, Rolain JM. Characterization of OXA-48-like-producing Enterobacteriaceae isolated from river water in Algeria. Water Res. 2017;120:185–9.
Al Fadhli AH, Jamal WY, Rotimi VO: Elucidating the virulence genes harboured by carbapenemase- and non-carbapenemase-producing carbapenem-resistant Klebsiella pneumoniae rectal isolates from patients admitted to intensive care units using whole-genome sequencing in Kuwait. J Med Microbiol 2022, 71(7):001554.
Potter RF, D’Souza AW, Dantas G. The rapid spread of carbapenem-resistant Enterobacteriaceae. Drug Resist Updates. 2016;29:30–46.
Hamidian M, Nigro SJ. Emergence, molecular mechanisms and global spread of carbapenem-resistant Acinetobacter baumannii. Microbial genomics. 2019;5(10): e000306.
Lee NY, Tsai CS, Syue LS, Chen PL, Li CW, Li MC, Ko WC. Treatment outcome of bacteremia due to non-carbapenemase-producing carbapenem-resistant Klebsiella pneumoniae bacteremia: role of carbapenem combination therapy. Clin Ther. 2020;42(3):e33–44.
Lee YQ. Sri La Sri Ponnampalavanar S, Chong CW, Karunakaran R, Vellasamy KM, Abdul Jabar K, Kong ZX, Lau MY, Teh CSJ: Characterisation of non-carbapenemase-producing carbapenem-resistant Klebsiella pneumoniae based on their clinical and molecular profile in Malaysia. Antibiotics. 2022;11(11):1670.
Hamzaoui Z, Ocampo-Sosa A, Fernandez Martinez M, Landolsi S, Ferjani S, Maamar E, Saidani M, Slim A, Martinez-Martinez L, Boutiba-Ben Boubaker I. Role of association of OmpK35 and OmpK36 alteration and bla(ESBL) and/or bla(AmpC) genes in conferring carbapenem resistance among non-carbapenemase-producing Klebsiella pneumoniae. Int J Antimicrob Agents. 2018;52(6):898–905.
Ma J, Song X, Li M, Yu Z, Cheng W, Yu Z, Zhang W, Zhang Y, Shen A, Sun H, et al. Global spread of carbapenem-resistant Enterobacteriaceae: Epidemiological features, resistance mechanisms, detection and therapy. Microbiol Res. 2023;266: 127249.
Lin L, Xiao X, Wang X, Xia M, Liu S. In vitro antimicrobial susceptibility differences between carbapenem-resistant KPC-2-producing and NDM-1-producing Klebsiella pneumoniae in a teaching hospital in northeast China. Microb Drug Resist. 2020;26(2):94–9.
Zhang R, Dong N, Huang Y, Zhou H, Xie M, Chan EW, Hu Y, Cai J, Chen S. Evolution of tigecycline- and colistin-resistant CRKP (carbapenem-resistant Klebsiella pneumoniae) in vivo and its persistence in the GI tract. Emerg Microbes Infect. 2018;7(1):127.
Wu Q, Liu Q, Han L, Sun J, Ni Y. Plasmid-mediated carbapenem-hydrolyzing enzyme KPC-2 and ArmA 16S rRNA methylase conferring high-level aminoglycoside resistance in carbapenem-resistant Enterobacter cloacae in China. Diagn Microbiol Infect Dis. 2010;66(3):326–8.
Su S, Li C, Zhao Y, Yu L, Wang Y, Wang Y, Bao M, Fu Y, Zhang J, Zhang X. Outbreak of KPC-2-producing Klebsiella pneumoniae ST76 isolates in an intensive care unit and neurosurgery unit. Microb Drug Resist. 2020;26(9):1009–18.
Chen L, Mathema B, Chavda KD, DeLeo FR, Bonomo RA, Kreiswirth BN. Carbapenemase-producing Klebsiella pneumoniae: molecular and genetic decoding. Trends Microbiol. 2014;22(12):686–96.
Ageevets V, Sopova J, Lazareva I, Malakhova M, Ilina E, Kostryukova E, Babenko V, Carattoli A, Lobzin Y, Uskov A, et al. Genetic environment of the blaKPC-2 gene in a Klebsiella pneumoniae isolate that may have been imported to russia from southeast Asia. Antimicrob Agents Chemother. 2017;61(2):e01856-e1916.
Wang LH, Wei DD, Wan LG, Yu Y, Deng Q, Liu Y. Diversity of the genetic environment of the blaKPC-2 gene among Klebsiella pneumoniae clinical isolates in a Chinese hospital. Microb Drug Resist. 2016;22(1):15–21.
Chen R, Liu Z, Xu P, Qi X, Qin S, Wang Z, Li R. Deciphering the epidemiological characteristics and molecular features of blaKPC-2- or blaNDM-1-positive Klebsiella pneumoniae isolates in a newly established hospital. Front Microbiol. 2021;12: 741093.
Partridge SR, Kwong SM, Firth N, Jensen SO. Mobile genetic elements associated with antimicrobial resistance. Clin Microbiol Rev. 2018;31(4):e00088-e117.
Zhang X, Li F, Cui S, Mao L, Li X, Awan F, Lv W, Zeng Z. Prevalence and distribution characteristics of blaKPC-2 and blaNDM-1 genes in Klebsiella pneumoniae. Infect Drug Resist. 2020;13:2901–10.
Huang J, Zhang S, Zhao Z, Chen M, Cao Y, Li B. Acquisition of a stable and transferable blaNDM-5-positive plasmid with low fitness cost leading to ceftazidime/avibactam resistance in KPC-2-producing Klebsiella pneumoniae during treatment. Front Cell Infect Microbiol. 2021;11: 658070.
Shaik S, Ranjan A, Tiwari SK, Hussain A, Nandanwar N, Kumar N, Jadhav S, Semmler T, Baddam R, Islam MA, et al. Comparative genomic analysis of globally dominant ST131 clone with other epidemiologically successful extraintestinal pathogenic Escherichia coli (ExPEC) lineages. Bio. 2017;8(5):e01596-17.
Marano RBM, Zolti A, Jurkevitch E, Cytryn E. Antibiotic resistance and class 1 integron gene dynamics along effluent, reclaimed wastewater irrigated soil, crop continua: elucidating potential risks and ecological constraints. Water Res. 2019;164: 114906.
Zhang S, Xu B, Chen M, Zhang Q, Huang J, Cao Y, Li B. Profile and actual transmissibility of carbapenem resistance genes: Intracellular and extracellular DNA in hospital wastewater. J Environ Manage. 2023;329: 117085.
Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, Aarestrup FM, Larsen MV. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother. 2012;67(11):2640–4.
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(7):3895–903.
Joensen KG, Scheutz F, Lund O, Hasman H, Kaas RS, Nielsen EM, Aarestrup FM. Real-time whole-genome sequencing for routine typing, surveillance, and outbreak detection of verotoxigenic Escherichia coli. J Clin Microbiol. 2014;52(5):1501–10.
McNally A, Oren Y, Kelly D, Pascoe B, Dunn S, Sreecharan T, Vehkala M, Välimäki N, Prentice MB, Ashour A, et al. Combined analysis of variation in core, accessory and regulatory genome regions provides a super-resolution view into the evolution of bacterial populations. PLoS Genet. 2016;12(9): e1006280.
Cazares A, Moore MP, Hall JPJ, Wright LL, Grimes M, Emond-Rhéault JG, Pongchaikul P, Santanirand P, Levesque RC, Fothergill JL, et al. A megaplasmid family driving dissemination of multidrug resistance in Pseudomonas. Nat Commun. 2020;11(1):1370.
Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009;26(7):1641–50.
Shintani M, Ohkuma M, Kimbara K: High-resolution comparison of bacterial conjugation frequencies. J Visual Exp 2019(143):e57812.
Nandanwar N, Janssen T, Kuhl M, Ahmed N, Ewers C, Wieler LH. Extraintestinal pathogenic Escherichia coli (ExPEC) of human and avian origin belonging to sequence type complex 95 (STC95) portray indistinguishable virulence features. Int J Med Microbiol. 2014;304(7):835–42.
Huang J, Zhang S, Zhang S, Zhao Z, Cao Y, Chen M, Li B. A Comparative study of fluoroquinolone-resistant Escherichia coli lineages portrays indistinguishable pathogenicity- and survivability-associated phenotypic characteristics between ST1193 and ST131. Infect Drug Resist. 2020;13:4167–75.
Jelić M, Hrenović J, Dekić S, Goić-Barišić I, Tambić Andrašević A. First evidence of KPC-producing ST258 Klebsiella pneumoniae in river water. J Hosp Infect. 2019;103(2):147–50.
We thank Department of Laboratory Medicine, Medical Technology and Engineering College, Fujian Medical University for kindly help.
This study was supported by Fujian Research and Training Grants for Young and Middle-aged Leaders in Healthcare.
Ethics approval and consent to participate
All procedures of this study involving humans (individuals, medical records, human samples, and clinical isolates) were reviewed and approved by the Medical Ethics Committee of Fujian Medical University Union Hospital (No.2020KY088). We confirm that this study was conducted in accordance with the Declaration of Helsinki. Informed consent was obtained from all subjects and legal guardians to participate in this study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1.Main microbiological characteristics of 11 clinical CRKP isolates. Table S2. AntimicrobialSusceptibility Testing of 11 clinical CRKP isolates. Table S3. Nucleotide sequences ofprimers used in this study.
About this article
Cite this article
Xie, Z., Huang, J., zhang, S. et al. Genomic and functional characterization of carbapenem-resistant Klebsiella pneumoniae from hospital wastewater. BMC Microbiol 23, 115 (2023). https://doi.org/10.1186/s12866-023-02862-5