Clonal spread of carbapenem-resistant Klebsiella pneumoniae among patients at admission and discharge at a Vietnamese neonatal intensive care unit
Antimicrobial Resistance & Infection Control volume 10, Article number: 162 (2021)
The increasing prevalence of carbapenem-resistant Enterobacteriaceae (CRE) is a growing problem globally, particularly in low- to middle-income countries (LMICs). Previous studies have shown high rates of CRE colonisation among patients at hospitals in LMICs, with increased risk of hospital-acquired infections.
We isolated carbapenem-resistant Klebsiella pneumoniae (CRKP) from faecal samples collected in 2017 from patients at admission and discharge at a Vietnamese neonatal intensive care unit (NICU). 126 CRKP were whole-genome sequenced. The phylogenetic relationship between the isolates and between clinical CRKP isolates collected in 2012–2018 at the same hospital were investigated.
NDM-type carbapenemase-(61%) and KPC-2-encoding genes (41%) were the most common carbapenem resistance genes observed among the admission and discharge isolates. Most isolates (56%) belonged to three distinct clonal clusters of ST15, carrying blaKPC-2, blaNDM-1 and blaNDM-4, respectively. Each cluster also comprised clinical isolates from blood collected at the study hospital. The most dominant ST15 clone was shown to be related to isolates collected from the same hospital as far back as in 2012.
Highly resistant CRKP were found colonising admission and discharge patients at a Vietnamese NICU, emphasising the importance of continued monitoring. Whole-genome sequencing revealed a population of CRKP consisting mostly of ST15 isolates in three clonally related clusters, each related to blood isolates collected from the same hospital. Furthermore, clinical isolates collected from previous years (dating back to 2012) were shown to likely be clonally descended from ST15 isolates in the largest cluster, suggesting a successful hospital strain which can colonise inpatients.
The prevalence of antibiotic-resistant bacteria is increasing globally, leading to elevated healthcare costs and patient mortality. The impact of antibiotic resistance is particularly heavy in low- and middle-income countries (LMICs) [1, 2]. Lower health indicators, higher rates of infection and inability to afford adequate treatment contribute to an elevated burden of infections caused by antibiotic-resistant bacteria in these settings. Furthermore, emergence and dissemination of antibiotic resistance is facilitated by unregulated access to antibiotics and a lack of resources to undertake extensive infection control measures . Data on the burden of carbapenem-resistant Enterobacteriaceae (CRE) is currently scarce in LMICs, however, Stewardson et al.  showed in a recent study, with data from ten LMICs, that bloodstream infections caused by CRE were associated with increased length of hospital stays and mortality. In Vietnam, prevalence of hospital-acquired infections (HAIs) at intensive care units (ICUs) has previously been reported to be high. A point prevalence study conducted in 2012 and 2013 found a total prevalence of 30%, with Klebsiella pneumoniae being the most common Enterobacteriaceae to cause HAIs, 15% of which were resistant to carbapenems . Carbapenems were the third most prescribed antibiotics, constituting 14% of total antibiotic prescriptions. At the ICUs of three Vietnamese paediatric referral hospitals during the same period, the rate of HAI was reported to be 33% with K. pneumoniae being the most common causative agent . 55% of K. pneumoniae isolates were reported to be resistant to carbapenems.
Colonisation of the intestinal tract with CRE has been shown to increase the risk of subsequent infection with CRE and to be associated with a higher mortality rate [7, 8]. In a recent study, we reported that 52% of patients admitted to 12 study hospitals in Vietnam were colonised with CRE, most commonly of K. pneumoniae (69%) . Risk factors for CRE colonisation included diagnosis with HAI and treatment with a carbapenem. The length of hospital stay was also found to be an independent risk factor. The mean colonisation rate increased from 13% on admission to 89% at day 15. A sub-study at one of the included neonatal ICUs (NICUs) showed ratios of colonisation of 32% at admission and 87% at discharge; a significantly higher mortality rate was observed among patients colonised with CRE. A study among neonatal children with HAIs caused by multidrug-resistant bacteria showed a significant correlation between antibiotic resistance and mortality . These findings prompted further investigation into carriage of carbapenem-resistant K. pneumoniae (CRKP) among the patients. In the current study, whole-genome sequencing (WGS) was used to characterise the CRKP colonising the intestinal tract of patients at the abovementioned NICU and to infer the phylogenetic relationships between strains colonising the patients at admission contra discharge. Isolates of CRKP from clinical samples were also characterised to compare the colonising bacteria to strains causing infections at the hospital.
Materials and methods
Collection of faecal samples, clinical isolates and patient baseline data
The samples for this study were prospectively collected from patients at a NICU at the largest paediatric hospital in northern Vietnam. The sampling procedure and collection of baseline patient data has been described elsewhere . In brief, between March and June 2017, faecal samples were collected from 326 neonates at admission to the NICU. Because the study hospital does not have a childbirth unit, a few of the admissions are from other wards with the majority of admissions being from other hospitals. Faecal swabs were collected for 248 of these patients at discharge from the NICU. The samples were collected with rectal swabs by using ESwabs (Copan, Brescia, Italy). All patients admitted to the NICU during the study period were included except for patients with anal atresia or other conditions affecting the feasibility of collecting the faecal sample or for whom a caretaker, legal guardian or relative chose to opt out of the study. Patient data collected included sex, date of birth, source of admission, reason for admission, diagnosis, and patient outcome. To compare the findings of CRKP in the faecal samples with the epidemiological situation at the hospital, clinical isolates of CRKP identified in the clinical routine work were retrospectively included in the study. These were constituted of 28 isolates collected from blood from 28 patients in 2017 and 2018, and 13 isolates collected from tracheal aspirates from 12 patients in 2017.
Screening for CRKP and antimicrobial susceptibility testing
The ESwabs were plated on CarbaID agar (bioMérieux, Marcy l’Étoile, France) and incubated at 35 ˚C for 16 to 24 h. Based on growth on the selective medium and morphology, suspected CRKP were species-determined by using a VITEK 2 (bioMérieux). Isolates determined as CRKP were recultivated on UTI Brilliance agar (Oxoid, Basingstoke, UK). Isolates were also cultured on COL-APSE agar (CHROMagar, Paris, France) to screen for colistin resistance. CRKP originating from the admission (n = 34) and discharge (n = 92) faecal samples, in addition to clinical isolates, were subjected to antimicrobial susceptibility testing to a panel of antibiotics consisting of meropenem, amikacin, gentamicin, tobramycin, cefotaxime, ciprofloxacin, fosfomycin, tigecycline, trimethoprim-sulfamethoxazole, chloramphenicol, tetracycline, ceftazidime-avibactam and colistin. Disc diffusion (Oxoid) on Mueller–Hinton agar plates (Oxoid) was used to test susceptibility to all antibiotics except for colistin, for which broth microdilution with a cation-adjusted Mueller–Hinton broth (BD, Franklin Lakes, NJ, USA) in custom made plates from Sensititre (TrekDS/Thermo Fisher Scientific, Oakwood Village, OH) was used. All zone diameters and MICs were interpreted according to the clinical breakpoints suggested by EUCAST  except for tetracycline, for which CLSI guidelines were used  and tigecycline. A wild type distribution of tigecycline zone diameters was defined from data from clinical isolates of extended-spectrum β-lactamase-(ESBL)-producing K. pneumoniae from Sweden, and isolates were considered resistant if zone diameters were outside of this distribution (i.e., < 15 mm). Multidrug resistance (MDR) was defined as non-susceptibility to antibiotics of three or more antimicrobial categories .
Whole-genome sequencing of CRKP, genome assembly and bioinformatic analysis
Total DNA from 149 CRKP isolates from admission and discharge faecal samples, blood and tracheal aspirates was extracted by using the EZ1 DNA Tissue Kit and the EZ1 Advanced XL instrument (Qiagen, Hilden, Germany) and subsequently stored in – 80 °C. The DNA concentration was determined by using a Qubit 2.0 Fluorometer (Thermo Fisher Scientific, Waltham, MA), and 20 ng of DNA was used as input to construct a sequencing library for whole-genome sequencing. The library preparation was performed by using a QIAseq FX DNA Library Kit (Qiagen) following the manufacturer’s instructions. DNA quantity and quality of the sequencing library was assessed by using a Qubit 2.0 Fluorometer and a QIAxcel instrument (Qiagen), followed by paired-end sequencing performed on a MiSeq instrument (Illumina, San Diego, CA). The raw reads generated were uploaded to the Sequence Read Archive at NCBI (accession numbers: SRR13256598-SRR13246708). An additional 18 isolates collected from admission and discharge faecal samples previously sequenced by using the same methods (accession numbers: SRR12149859-SRR12149874) were also analysed. Genome assembly of the whole-genome sequenced isolates was performed by using CLC Genomics Workbench v.9.5.3 (Qiagen). Multilocus sequence typing (MLST) and querying of antibiotic resistance genes were performed via the databases at the Center for Genomic Epidemiology (http://www.genomicepidemiology.org). Single nucleotide polymorphisms (SNPs) were called against a reference genome (NC_022566) if fulfilling the criteria of having a sequencing depth ≥ 10x, a frequency of ≥ 90% and a Phred score ≥ 20. The resulting variants were used to construct a phylogeny of the sequenced isolates by using CLC Genomics Workbench V.9.5.3. To compare the genetic relatedness of the study isolates to CRKP isolated prior to 2017 at the study hospital from clinical specimen, 5 isolates from 2015 previously reported on elsewhere , 4 isolates from 2014, 9 isolates from 2013 and 1 isolate from 2012, all belonging to ST15, were included in the analysis. Isolates differing by ≤ 18 SNPs were considered to constitute a single clone . All raw reads generated and presented in this study are available at the Sequence Read Archive at NCBI (accession numbers: SRR13256598-SRR13246708).
The associations between patient baseline variables sex, treatment outcome and diagnoses with carriage of CRKP were tested with Fisher’s exact tests and χ2-tests. Associations between age at admission and length of stay at the NICU with carriage of CRKP were tested with Mann–Whitney U and Kruskal–Wallis one-way analysis of variance tests. Differences in antibiotic susceptibility were assessed with Fisher’s exact tests. To compensate for multiple comparison errors, P ≤ 0.01 was considered significant. Due to the low number of patients which were positive at admission but negative at discharge (n = 5), this cohort was excluded from the statistical analysis. All statistical tests were performed in Prism for Windows 8.4.3 and SPSS Statistics version 26.
Prevalence of CRKP at admission and discharge and baseline patient statistics
Faecal samples were collected at admission and discharge from a total of 326 patients at the study NICU (Table 1). A total of 69 patients (21%) had faecal samples positive for CRKP at admission and 184 had faecal samples positive for CRKP at discharge (74%). Only one patient which had a faecal sample positive for CRKP at admission had a negative faecal sample at discharge (2.0%). 137 patients (55%) were negative for CRKP at admission but were positive at discharge (i.e., were colonized during their stay at the NICU). A total of 59 patients (24%) had faecal samples which neither tested positive for CRKP at admission nor at discharge.
Overall, 37% of the patients were female (n = 118). The overall median age of patients in the study was 4 days. Patients CRKP-positive at admission were significantly older compared to patients CRKP-negative at admission (P < 0.001) and patients both CRKP-positive at admission and discharge were significantly older (P < 0.01) than those CRKP-negative at both admission and discharge and those positive only at discharge. The median length of stay at the NICU overall was 5 days. Patients which were CRKP-positive at discharge had significantly longer NICU stays compared to patients which were CRKP-negative at discharge (P < 0.001). Furthermore, a significant difference was observed when comparing both admission and discharge status (P < 0.01); patients CRKP-negative at both admission and discharge had shorter NICU stays. Patients transferred from the NICU were found to be significantly likelier to be admission-negative (P < 0.001) and discharge-negative (P < 0.01). Regarding diagnoses, patients with CRKP-positive faecal samples at either admission or discharge were more likely to be afflicted with respiratory distress syndrome (P < 0.01), and patients with congenital heart disease (CHD) were found to be less likely to be admission-positive for CRKP (P < 0.001).
Antimicrobial susceptibility of CRKP isolates from admission and discharge samples
Data on antimicrobial susceptibility is presented in Table 2. Isolates from patients at admission did not significantly differ in antibiotic susceptibility compared to isolates from patients at discharge who were negative at admission (data no shown). None of the admission- or discharge-isolates were susceptible to meropenem (which was an inclusion criterion). Cefotaxime susceptibility rates were also accordingly at 0% for each included isolate. All admission- and discharge-isolates (except one) were MDR. The susceptibility rates for aminoglycosides were low; for gentamicin, amikacin and tobramycin, the susceptibility rates were 41%, 8.8% and 0% for admission-isolates and 36%, 4.3% and 1.1% for discharge-isolates, respectively. Susceptibility ratios were also low for ciprofloxacin and fosfomycin; for admission-isolates, susceptibility ratios were 15% and 5.9% respectively, whereas for discharge-isolates, the corresponding ratios were 6.5% for each antibiotic. Ceftazidime/avibactam susceptibility ratios were also fairly low, 35% of the admission-isolates and 39% of the discharge-isolates were susceptible. The highest susceptibility ratios observed were for tigecycline and colistin. 100% and 96% of the admission- and discharge-isolates were susceptible to tigecycline whereas 82% and 89% were susceptible to colistin, respectively. Notably, for 4 colistin-susceptible admission isolates (12%), heteroresistance was observed (i.e., the main population was susceptible), however, a resistant subpopulation could be isolated from the COL-APSE agar. The colistin-heteroresistant isolates were denoted as resistant in the summary of the antimicrobial susceptibility data.
Antibiotic resistance genes among CRKP from admission and discharge samples
Although blaKPC-2, which was carried by 35% and 44% of the isolates collected at admission and discharge, respectively, was the single most common carbapenemase gene observed among the faecal isolates, genes encoding NDM-type carbapenemases were more common taken together (Fig. 1). For admission-isolates, blaNDM-1 (32%), blaNDM-4 (24%) and blaNDM-5 (8.8%) together constituted 65% if the carbapenemase genes observed whereas for discharge-isolates, the abovementioned NDM-type genes constituted 60% of all carbapenemase genes, with blaNDM-4 being the most common (40%), followed by blaNDM-1 (15.2%) and blaNDM-5 (4.3%). The carbapenemase gene blaOXA-181 was also observed among a few of the isolates (2.9% of the admission-isolates and 3.3% of the discharge-isolates, respectively). Plasmid-mediated ampC-genes were also detected. Among admission isolates, 4 carried blaDHA-1 (12%) and 2 carried blaCMY-6 (5.9%). Among discharge isolates, 1 carried blaDHA-1 (1.1%) and 2 carried blaCMY-6 (2.2%). blaCTX-M-14 and blaCTX-M-15 were the most common CTX-M-genes detected; their prevalences were, respectively, 21% and 29% among admission isolates, and 33% and 29% for discharge isolates.
A wide range of genes conferring resistance to antibiotics other than β-lactams were observed (Fig. 1). Notably, genes encoding 16S rRNA methylases, conferring high grade aminoglycoside resistance, were carried by 38% of admission-isolates and 20% of discharge-isolates. The most common were rmtB, followed by rmtC and armA, which were carried by 27%, 8.8% and 2.9% of the admission-isolates and 15%, 3.3% and 1.1% of the discharge-isolates, respectively. The colistin resistance gene mcr-8 was detected in one admission-isolate; the isolate was resistant to colistin and the MIC for the isolate was 16 mg/L.
Multilocus sequence typing of CRKP isolates from admission and discharge samples
Typing data based on MLST is presented in Table 3. Among both admission- and discharge-isolates, ST15 was found to be the predominant ST. Among the admission-isolates, 14 isolates (41%) were found to belong to ST15. This was followed by ST11, to which 4 isolates (12%) belonged. From the discharge samples, 57 isolates (62%) belonged to ST15. The second most common ST among these isolates was ST14, to which 10 isolates (11%) belonged. This ST was followed in prevalence by ST11 and ST16, to which 7 (7.6%) and 4 isolates (4.3%) belonged, respectively. When the discharge-isolates were stratified on isolates originating from patient with negative admission samples (n = 59), ST15 was found at a significantly higher ratio of 75% (n = 44) (P = 0.0015). The second most common ST among these isolates, ST16 was found only among 5 isolates (8.5%). By comparing the STs of isolates from admission samples to isolates from discharge samples originating from the same patient, 19 patients (59%) had faecal isolates belonging to the same ST.
Characteristics of CRKP isolated from blood
The 28 bloodstream isolates, collected in 2017 and 2018, showed a predominance of NDM-type carbapenemase genes (Additional file 1: Fig. S1), 26 isolates (93%) carried a blaNDM-gene among which blaNDM-1 (43%) and blaNDM-4 (43%) were the most and equally common, followed by blaNDM-5 (7.1%). Only two isolates carried other carbapenemase genes, one carried blaKPC-2 (3.6%) and the other carried blaOXA-181 (3.6%). The plasmid-mediated ampC-gene blaDHA-1 was detected among 21.4% of the blood-isolates. blaCTX-M-15 (50%) and blaCTX-M-14 (39%) were the most common CTX-M-genes although blaCTX-M-9 (3.6%) was also detected in 1 isolate. The high-grade aminoglycoside resistance gene rmtB was carried by 8 isolates (29%). The susceptibility testing showed that all blood isolates were MDR (Table 2). The aminoglycoside susceptibility among the blood isolates was very low. The susceptibility ratio to gentamicin was 18%, which was the highest for any aminoglycoside. The susceptibility ratio for ceftazidime/avibactam was only 7.1%. Ciprofloxacin and fosfomycin susceptibility ratios were also low, only 3 isolates (11%) were susceptible to each antibiotic. The susceptibility ratio was highest to tigecycline (100%) and colistin (93%), both which were the only antibiotics for which susceptibility ratios were over 90%. Notably, two isolates displayed heteroresistance to colistin as discovered by plating the isolates on COL-APSE agar plates. These isolates were denoted as colistin-resistant in the summary of the data. All blood isolates were susceptible to tigecycline. The MLST analysis of the blood isolates showed that ST15 was the most common ST, to which 5 isolates (18%). This ST was followed in prevalence by ST14 and ST397 (14%, respectively), and ST17 (11%).
Characteristics of CRKP isolated from tracheal aspirates
Among the isolates from tracheal aspirates collected in 2017 and 2018 from the study hospital, 7 carried blaKPC-2 (54%) and 6 carried NDM-type carbapenemase-encoding genes (Additional file 1: Fig. S1); 3 carried blaNDM-1 (23%) and 3 carried blaNDM-4 (23%). blaCTX-M-15 (46%) and blaCTX-M-14 (23%) were the only CTX-M-genes found among these isolates. The susceptibility testing showed that all tracheal aspirate isolates were MDR (Table 2). Although none of the tracheal aspirate isolates were susceptibility to amikacin and tobramycin, 54% were susceptible to gentamicin. None of the isolates were susceptible to ciprofloxacin and fosfomycin. The susceptibility ratio was highest for tigecycline (100%) and colistin (85%). The tracheal aspirate isolates were found to belong to either of two STs, 10 isolates (77%) belonged to ST15 and 3 isolates (23%) belonged to ST14.
Genetic relatedness of CRKP isolates
The phylogenetic tree based on the SNP analysis of the isolates belonging to ST15 showed that most isolates clustered in three different clusters (Fig. 2), separated by a minimum distance of > 200 SNPs. Isolates within each cluster differed by ≤ 18 SNPs, indicating that the three clusters comprised three different clones. The largest cluster, consisting of 48 isolates in total, was comprised of 8 admission isolates, 33 discharge isolates, 1 blood isolate and 6 tracheal aspirate isolates. All isolates in this cluster carried blaKPC-2. Seven of the 8 patients from which the admission isolates originated were also positive for CRKP at discharge. 4 of the patients had discharge isolates which clustered together with the admission isolates. The discharge isolate from 1 patient was also ST15 but differed > 1,800 SNPs from the admission isolate whereas 2 patients carried isolates at discharge which had different STs compared to the admission isolate (ST11 and a novel ST, respectively). The clinical ST15 isolates collected 2012, 2013, 2014 and 2015 formed clusters close to the largest ST15 cluster. The differences between the most closely related ST15 isolates from 2017 and 2015 were 27 SNPs, 2015 and 2014 were 2 SNPs, 2014 and 2013 were 47 SNPs, and 2013 and 2012 were 16 SNPs. The second largest cluster was comprised of 21 isolates, all of which carried blaNDM-4. 2 isolates in this cluster were recovered from admission samples, 17 from discharge samples and 2 from blood. The third largest cluster consisting of 8 isolates in total, was comprised of 3 discharge isolates, 2 blood isolates and 3 tracheal aspirate isolates.
SNP-analyses were also carried out on isolates belonging to numerous ST11, ST14, ST16 and ST334 isolates (Additional file 2–5: Fig. S2–S5). The ST11-isolates comprised two clusters (maximum distance within clusters were 1 and 2 SNPs, respectively) separated by > 4,000 SNPs. The isolates comprising the larger cluster (n = 9) all carried blaNDM-4, whereas the 2 isolates comprising the smaller cluster all carried blaNDM-1. The ST14-isolates (n = 18) formed a single cluster with a ≤ 15 SNPs difference, all of which carried blaNDM-4. The ST16-isolates (n = 5) clustered together with an outlier isolate differing 45 SNPs from its closest relative. Three isolates carried both blaKPC-2 and blaNDM-4 whereas the outlier carried blaOXA-181. Interestingly, 1 isolate, differing only 2 SNPs from its closest relative, carried blaNDM-1 and neither blaKPC-2 nor blaNDM-4. The ST334-isolates (n = 8) formed a closely related cluster with isolates carrying blaNDM-4, with an outlier isolate carrying blaNDM-1 which differed 170 SNPs from its closest relative.
In this study, faecal samples collected at admission and discharge from patients at a NICU in a large Vietnamese hospital were screened for CRKP and the isolates were further characterised with WGS. The results showed a 21% prevalence of CRKP at admission and an increase to 74% at discharge. Statistical analyses revealed several significant associations. Admission-positive patients were significantly older compared with admission-negative patients (P < 0.001). Due to the relatively low age of the included patients (median age of all patients was 4 days) this is likely reflecting the longer exposure time of older neonates to the environment. Similarly, significantly longer stays at the NICU (P < 0.001) were observed among patients discharge-positive for CRKP, likely reflecting an increased risk by longer exposure time in the NICU environment. Patients transferred from the NICU to other wards of the hospital were found to be significantly likelier to be admission-negative (P < 0.001) and discharge-negative (P < 0.01) which may reflect more adverse conditions for patients colonized with CRKP, requiring them to stay at the NICU. Respiratory distress syndrome was found to more likely afflict patients who were either positive for CRKP at admission or discharge or both (P < 0.01), indicating a possible connection between colonisation and respiratory infection. A significant association was also found that patients with CHD were less likely to be negative for CRKP at admission (P < 0.001). Although not significant (P = 0.089), patients with CHD tended to be older at admission (median age = 12 days) compared with children without CHD (median age = 4 days), which could explain this association.
A total of 34 isolates from faecal samples collected at admission and 92 collected at discharge were analysed with whole-genome sequencing. MLST and phylogenetic analyses based on SNPs were used to evaluate the phylogenetic relationships between the isolates. ST15 was the most common ST among both admission and discharge isolates. 75% of the discharge isolates from patients which were negative for CRKP at admission belonged to ST15, indicating that the patients colonised during their hospital stay were predominantly colonised with ST15 strains. ST15 is an emerging, global high-risk ST and strains of this ST has been reported on in hospital outbreaks worldwide . Clinical ST15-isolates are often ESBL-producers but have been observed as carbapenemase-producers in several outbreaks. Phylogenetic analysis of the ST15 isolated in the current study showed that most isolates clustered in three different clusters comprised of clones carrying blaKPC-2, blaNDM-1 and blaNDM-4, respectively. To investigate whether the admission- and discharge-isolates were genetically related to clinical isolates, we sequenced and compared blood and tracheal aspirates collected at the study hospital in 2017 and 2018. Each of the clusters consisting of ST15 admission- and discharge-isolates additionally had clinical isolates which clustered together, indicating the potential of the colonising isolates to cause infections. A previous study at the study hospital carried out in 2015 detailed the characteristics of ST15 CRKP which constituted the most common ST among CRKP in clinical specimen . The ST15 from 2015 were found to be comprised of two clones. Representative isolates of these two clones were in the current study shown to be closely related to the ST15 CRKP isolates clustered in the largest cluster among the ST15 isolated in 2017 (all of which carried blaKPC-2). The closest relatives were separated by only 27 SNPs, indicating that these isolates belong to the same clone, although separated by differences accumulated over two years. Phylogenetic analysis of sporadic isolates of ST15 CRKP carrying blaKPC-2 recovered in 2012, 2013 and 2014 from clinical specimen at the same hospital further showed a general clustering of the isolates on sampling year, but the close distance between the clusters and the largest cluster of ST15 CRKP from 2017 could indicate clonal descent from the isolate collected in 2012. The two other major clusters of ST15 CRKP (carrying blaNDM-1 and blaNDM-4, respectively) from the current study were clearly separated both from each other and the major cluster comprised of isolates carrying blaKPC-2, indicating the emergence of these two new clones at the hospital. Notably, the phylogenetic analyses also revealed the occurrence of mostly genetically homogenous isolates of other prevalent STs, including ST11, ST14, ST16 and ST334, carrying various carbapenemase resistance genes. These data show that the CRKP at the study hospital are comprised of a diverse set of clones. The data suggests that the ST15 clone has been established and well-adapted to the hospital environment since at least as early as 2015, which may explain its continued predominance in 2017. However, the influx of new clones, including genetically divergent strains of ST15, indicate the need for continued monitoring of potential changes in the epidemiological situation.
All CRKP included in the study were MDR. Resistance ratios to prevalently used antibiotics including aminoglycosides, ciprofloxacin and fosfomycin were high, indicating the difficulty of empirical and culture-based treatment. Notably, a considerable number of isolates carried 16S rRNA methylase genes conferring high grade resistance to aminoglycosides. In our earlier study on isolates collected at the hospital in 2015 , we observed a predominance of blaKPC-2 with a corresponding high ratio of susceptibility to ceftazidime/avibactam. With the emergence of clones carrying blaNDM-genes, outnumbering blaKPC-2-carrying isolates in the current data, the susceptibility ratio to ceftazidime/avibactam was considerably lower, reaching as low as 7.1% among blood isolates (Table 2). Colistin is another drug regarded as a last-resort antibiotic for MDR CRKP. The relatively high susceptibility rate observed among isolates in this study (> 80%) suggests colistin’s applicability as such in the current context. Regardless, colistin should be used carefully to avoid selection of resistant clones. Aside from chromosomally mediated colistin resistance (which was the most common mechanism observed), one admission isolate carried the plasmid-mediated colistin resistance gene mcr-8 and isolates of CRKP with mcr-1 have previously been reported among clinical isolates at the study hospital . We have previously reported on a high colistin resistance ratio (42%) among the major ST15 clone at the study hospital , mainly due to chromosomal alterations, further showing that use of colistin may result in selection of resistant clones and subclones. Tigecycline was the antibiotic to which the highest susceptibility ratio was observed. However, since no clinical breakpoints are available for tigecycline, susceptibility was defined based on an ecological cut-off determined from a wild type population . By using appropriately determined clinical breakpoints, it is possible that the susceptibility would be considered substantially lower. Furthermore, although tigecycline has been used as salvage therapy in critically ill children infected by MDR bacteria, the safety of the antibiotic has not been clearly established in the patient group . Further studies assessing the efficacy and safety of paediatric use of tigecycline are urgently needed.
Notably, no significant difference in antibiotic susceptibility was observed among isolates from patients at admission compared to isolates from patients at discharge which were negative at admission. This likely reflects that CRKP from the study hospital and many CRKP from other hospitals (from which most patients were admitted) share a high degree of genetic homogeneity. In a recent study conducted in four hospitals in Hanoi, CRKP of the same lineages were isolated from three hospitals with patient transfer . However, in the fourth hospital, which did not transfer patients with other hospitals, CRKP were of different lineages, supporting the hypothesis that patient transfer facilitates the dissemination of CRKP strains between hospitals.
In this study, CRKP from faecal samples of patients at admission and discharge at a large Vietnamese NICU were analysed with whole-genome sequencing and most isolates were shown to belong to ST15. Most of these isolates clustered in three different clonal clusters, each also consisting of blood isolates collected at the same hospital. Furthermore, clinical isolates collected from previous years (dating back to 2012) were shown to likely be clonally descended from ST15 isolates in the largest cluster, suggesting a successful hospital strain which can colonise inpatients. The high levels of resistance among the CRKP isolates emphasize the importance of antimicrobial susceptibility testing of clinical isolates of K. pneumoniae to ascertain an efficient treatment regimen, as the data does not indicate an unequivocal empirical regimen. Aside from prudent usage of antibiotics, infection control measures are also urgently needed to limit the number of HAIs overall. To ensure clinical readiness, continued monitoring of the epidemiological situation is required.
Availability of data and materials
All raw reads generated and presented in this study are available at the Sequence Read Archive at NCBI (accession numbers: SRR13256598-SRR13246708).
Clinical & Laboratory Standards Institute
Carbapenem-resistant Klebsiella pneumoniae
European Committee on Antimicrobial Susceptibility Testing
Intensive care unit
Low- and middle-income countries
Multilocus sequence typing
Neonatal intensive care unit
Limmathurotsakul D, Dunachie S, Fukuda K, Feasey NA, Okeke IN, Holmes AH, et al. Surveillance and epidemiology of drug resistant infections consortium (SEDRIC). Improving the estimation of the global burden of antimicrobial resistant infections. Lancet Infect Dis. 2019;19:e392–8.
Pokharel S, Raut S, Adhikari B. Tackling antimicrobial resistance in low-income and middle-income countries. BMJ Glob Health. 2019;4:e002104.
Li G, Bielicki JA, Ahmed ASMNU, Islam MS, Berezin EN, Gallacci CB, et al. Towards understanding global patterns of antimicrobial use and resistance in neonatal sepsis: insights from the NeoAMR network. Arch Dis Child. 2020;105:26–31.
Stewardson AJ, Marimuthu K, Sengupta S, Allignol A, El-Bouseary M, Carvalho MJ, et al. Effect of carbapenem resistance on outcomes of bloodstream infection caused by Enterobacteriaceae in low-income and middle-income countries (PANORAMA): a multinational prospective cohort study. Lancet Infect Dis. 2019;19:601–10.
Phu VD, Wertheim HF, Larsson M, Nadjm B, Dinh QD, Nilsson LE, et al. Burden of hospital acquired infections and antimicrobial use in Vietnamese adult intensive care units. PLoS ONE. 2016;11:e0147544.
Le NK, Hf W, Vu PD, Khu DT, Le HT, Hoang BT, et al. High prevalence of hospital-acquired infections caused by gram-negative carbapenem resistant strains in Vietnamese pediatric ICUs: a multi-centre point prevalence survey. Medicine (Baltimore). 2016;95:e4099.
Tischendorf J, de Avila RA, Safdar N. Risk of infection following colonization with carbapenem-resistant Enterobacteriaceae: a systematic review. Am J Infect Control. 2016;44:539–43.
McConville TH, Sullivan SB, Gomez-Simmonds A, Whittier S, Uhlemann AC. Carbapenem-resistant Enterobacteriaceae colonization (CRE) and subsequent risk of infection and 90-day mortality in critically ill patients, an observational study. PLoS ONE. 2017;12:e0186195.
Tran DM, Larsson M, Olson L, Hoang NTB, Le NK, Khu DTK, et al. High prevalence of colonization with carbapenem-resistant Enterobacteriaceae among patients admitted to Vietnamese hospitals: risk factors and burden of disease. J Infect. 2019;79:115–22.
Peters L, Olson L, Khu DTK, Linnros S, Le NK, Hanberger H, et al. Multiple antibiotic resistance as a risk factor for mortality and prolonged hospital stay: a cohort study among neonatal intenstive care patients with hospital-acquired infections caused by gram-negative bacteria in Vietnam. PLoS ONE. 2019;14:e0215666.
The European Committee on Antimicrobial Susceptibility Testing. Breakpoint tables for interpretation of MICs and zone diameters. Version 9.0, 2019. http://www.eucast.org.
CLSI. Performance Standards for Antimicrobial Susceptibility Testing. 28th ed. CLSI supplement M100. Wayne, PA: Clinical and Laboratory Standards Institute; 2018.
Magiorakos AP, Srinivasan A, Carey RB, Carmeli Y, Falagas ME, Giske CG, et al. Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin Microbiol Infect. 2012;18:268–81.
Berglund B, Hoang NTB, Tärnberg M, Le NK, Nilsson M, Khu DTK, et al. Molecular and phenotypic characterization of clinical isolates belonging to a KPC-2-producing strain of ST15 Klebsiella pneumoniae from a Vietnamese pediatric hospital. Antimicrob Resist Infect Control. 2019;8:156.
Schürch AC, Arredondo-Alonso S, Willems RJL, Goering RV. Whole genome sequencing options for bacterial strain typing and epidemiologic analysis based on single nucleotide polymorphism versus gene-by-gene-based approaches. Clin Microbiol Infect. 2018;24:350–4.
Navon-Venezia S, Kondratyeva K, Carattoli A. Klebsiella pneumoniae: a major worldwide source and shuttle for antibiotic resistance. FEMS Microbiol Rev. 2017;41:252–75.
Berglund B, Hoang NTB, Tärnberg M, Le NK, Welander J, Nilsson M, et al. Colistin- and carbapenem-resistant Klebsiella pneumoniae carrying mcr-1 and blaOXA-48 isolated at a paediatric hospital in Vietnam. J Antimicrob Chemother. 2018;73:1100–2.
Berglund B, Hoang NTB, Tärnberg M, Le NK, Svartström O, Khu DTK, et al. Insertion sequence transpositions and point mutations in mgrB causing colistin resistance in a clinical strain of carbapenem-resistant Klebsiella pneumoniae from Vietnam. Int J Antimicrob Agents. 2018;51:789–93.
The European Committee on Antimicrobial Susceptibility Testing. MIC and zone diameter distributions and ECOFFs. http://www.eucast.org. Accessed 23 Apr 2020.
Iosifidis E, Violaki A, Michalopoulou E, Volakli E, Diamanti E, Koliouskas D, et al. Use of tigecycline in pediatric patients with infections predominantly due to extensively drug-resistant gram-negative bacteria. J Pediatric Infect Dis Soc. 2017;6:123–8.
Linh TD, Thu NH, Shibayama K, Suzuki M, Yoshida L, Thai PD, et al. Expansion of KPC-producing Enterobacterales in four large hospitals in Hanoi. Vietnam J Glob Antimicrob Resist. 2021;S2213–7165(21):00213–7.
We would like to thank the staff at neonatal and microbiology departments of Vietnam National Children’s Hospital for their hard work, Rogier van Doorn for his suggestions on how to conduct the study, and Newborn Vietnam for their support, teaching and valuable comments regarding the screening procedure during the study.
Open access funding provided by Linköping University. This study was supported financially by Training and Research Academic Collaboration (TRAC)—Sweden—Vietnam, the Swedish Foundation for International Cooperation in Research and Higher Education (STINT) (SG 2015-5972), and ALF grants, Region Östergötland (LIO-793191 and LIO-899871).
Ethics approval and consent to participate
Ethics approval for this study was granted by the Ethics Committee of the Vietnam National Children’s Hospital (reference number: VNCH-RICH-16-014 with approved extension/appendix). Caretakers, legal guardians, and relatives of the patients included in this study were given the opportunity to opt out of the study at any time.
Consent for publication
In this study, no individually identifiable data are included, all data are presented as pooled data. As such, consent for publication was not required. Nonetheless, caretakers, legal guardians, and relatives of all included patients were fully informed of the conditions of the study and were given the opportunity to opt out of the study at any time.
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This paper is dedicated to the memory of the late Hai Thanh Le
: Fig. S1. Antibiotic resistance genes detected in isolates of carbapenem-resistant Klebsiella pneumoniae originating from clinical isolates collected from blood and tracheal fluid samples. For clarity, certain gene variants (specified below) have been aggregated into gene categories. These are, and consist of: blaSHV (non-ESBL): blaSHV-28, blaSHV-33, blaSHV-56, blaSHV-65, blaSHV-67, blaSHV-73, blaSHV-79, blaSHV-96, blaSHV-110, blaSHV-111, blaSHV-133, blaSHV-142, blaSHV-165, blaSHV-172, blaSHV-179, blaSHV-182, blaSHV-194, blaSHV-196; blaSHV (ESBL): blaSHV-12, blaSHV-13, blaSHV-30, blaSHV-38, blaSHV-42, blaSHV-99, blaSHV-106; blaTEM: blaTEM-1B, blaTEM-216; blaOXA: blaOXA-1, blaOXA-9, blaOXA-10; aac: aac(3)-IIa, aac(3)-IId, aac(6’)-Ib, aac(6’)-Ib3, aac(6’)-Ib-cr; aad: aadA1, aadA2b, aadA3, aadA16, aadA5; ant: ant(2″)-Ia, ant(3’)-VI; aph: aph(3″)-Ib, aph(6)-Id; blaOKP: blaOKP-B-2, blaOKP-B-3, blaOKP-B-4, blaOKP-B-8, blaOKP-B-10, blaOKP-B-14; fosA: fosA, fosA3, fosA5; oqxAB: oqxA, oqxB; qnrB: qnrB1, qnrB4, qnrB6, qnrB9; catB: catB3, catB8; arr: arr-2, arr-6; drfA: dfrA1, dfrA7, dfrA12, dfrA14, dfrA23, dfrA27.
: Fig. S2. The genetic relatedness based on single-nucleotide polymorphisms (SNPs) of 11 carbapenem-resistant Klebsiella pneumoniae (CRKP) belonging to ST11 collected from faecal samples from patients at the study hospital was illustrated in a phylogenetic tree. The closest distances between these clusters are denoted by the closest related inter-cluster relatives as indicated along the horizontal line. The maximum number of SNPs along the clusters indicate the largest distance between two isolates in the indicated cluster.
: Fig. S3. The genetic relatedness based on single-nucleotide polymorphisms (SNPs) of 18 carbapenem-resistant Klebsiella pneumoniae (CRKP) belonging to ST14 collected from faecal samples from patients at the study hospital was illustrated in a phylogenetic tree. The closest distances between these clusters are denoted by the closest related inter-cluster relatives as indicated in parenthesis. The maximum number of SNPs along the clusters indicate the largest distance between two isolates in the indicated cluster.
: Fig. S4. The genetic relatedness based on single-nucleotide polymorphisms (SNPs) of 5 carbapenem-resistant Klebsiella pneumoniae (CRKP) belonging to ST16 collected from faecal samples from patients at the study hospital was illustrated in a phylogenetic tree. The closest distances between these clusters are denoted by the closest related inter-cluster relatives as indicated along the vertical line. The maximum number of SNPs along the clusters indicate the largest distance between two isolates in the indicated cluster.
: Fig. S5. The genetic relatedness based on single-nucleotide polymorphisms (SNPs) of 8 carbapenem-resistant Klebsiella pneumoniae (CRKP) belonging to ST334 collected from faecal samples from patients at the study hospital was illustrated in a phylogenetic tree. The closest distances between these clusters are denoted by the closest related inter-cluster relatives as indicated along the vertical line. The maximum number of SNPs along the clusters indicate the largest distance between two isolates in the indicated cluster.
About this article
Cite this article
Berglund, B., Hoang, N.T.B., Lundberg, L. et al. Clonal spread of carbapenem-resistant Klebsiella pneumoniae among patients at admission and discharge at a Vietnamese neonatal intensive care unit. Antimicrob Resist Infect Control 10, 162 (2021). https://doi.org/10.1186/s13756-021-01033-3