The present survey revealed that the picture of the antibiotic resistance, virulence and genetic relationship between clinical strains of K. pneumoniae recovered from the patients admitted to the 12th Bahman Hospital and other laboratories in Shahr-e-Qods. K. pneumoniae causes several types of infections in humans, including respiratory, bloodstream and urinary tract infections (UTIs), which are commonly seen in hospitalized or immunocompromised patients [36,37]. These infections are often treated with beta-lactams and other effective antibiotics against Enterobacteriaceae. Nevertheless, the antibiotic-resistant and highly pathogenic species of K. pneumoniae are rapidly spreading around the world [38]. Bacterial resistance depends vastly on populational and geographical factors. Thus, the body of different people can provide different environments for the bacteria to grow, multiply and be affected by the drugs [39]. Therefore, the study of bacterial resistance in a specific population can provide an appropriate overview of effective drugs for the healthcare staff to provide an effective antibiotic regimen to ensure improved recovery of patients.
Our findings demonstrated that out of 84 isolates identified as K. pneumoniae, more than 50% of the samples had multi-drug resistance, with the highest resistance against ceftriaxone, cefotaxime and Fosfomycin, respectively. Notably, in the tested isolates we observed a high susceptibility to colistin, nitrofurantoin, azithromycin, imipenem, and tetracycline. In line with the present study, Rocha et al., [40] reported that K. pneumoniae ESBL-positive isolates sampled from ICU patients showed the highest resistance to ceftriaxone. Kim Det al., [41] found that the resistance rates of Klebsiella pneumoniae to cefotaxime, cefepime, and carbapenem were 38-41%, 33-41%, and <0.1-2%, respectively, from 2013 to 2015. According to them, all isolates are sensitive to imipenem. The highest susceptibility was reported to gentamicin and piperacillin/tazobactam, respectively. In a different report compared to our results, the highest sensitivity was observed against ceftriaxone, ciprofloxacin and gentamicin, respectively. However, in the present study, the highest resistance was to ceftriaxone [42].
Moreover, our results of the combined disk test indicated that out of 52 resistant samples, 27 (52%) cases were positive for ESBL, almost corroborating evidence with the findings of a study by Rupinder et al., [43] who stated that ESBL production was observed in 48% of E. coli, 44% of K. pneumoniae and 50% of P. aeruginosa isolates in a tertiary hospital in Patiala, Punjab. While according to studies conducted in Iraq, 81.39% of K. pneumoniae isolates were ESBL producers [44]. This discrepancy between the reported frequencies can be due to the geographical diversity, the type of strains and the isolation under study or other different reasons.
The results of the genes involved in ESBL resistance showed that the prevalence of blaCTX-M, blaSHV and blaTEM genes out of 84 samples were 81 (96%), 79 (94%) and 77(91%), respectively. Therefore, the highest frequency was related to blaCTX-M gene. In the study of Pishtiwan et al., [45] the frequency of blaTEM, blaSHV and blaCTX-M genes was 64.7%, 35.2% and 41.1%, respectively. Also, according to Ugbo et al., [46] the prevalence of blaSHV, blaTEM and blaCTX-M was identified to be 55%, 35% and 45% respectively. These values were lower than the percentages of the present study and these results disagreed with the reports of studies by Bajpai et al., [47] in which the gene that predominated was blaTEM (48.7%), followed by blaTEM (7.6%) and blaSHV (5.1%). It is worth mentioning that the genotypic approach can be the method of choice for distinguishing ESBL-producing strains from Enterobacteriaceae because phenotypic tests for ESBL detection only confirm ESBL production.
Analysis of the detection of KPC resistance genes showed that blaKPC (71%) was ampler than blaNDM (60%) without any significant differences in resistant species. Liu et al., [48] stated that among the tested isolates of K. pneumoniae,blaNDM and blaKPC genes had the highest frequency, respectively. This discrepancy may be attributed to genetic variations among strains associated with human populations and antibiotic regimens. In accordance with our data, Xiufeng et al., [49] reported that blaKPC and blaNDM were highly detected in carbapenems-resistant K. pneumoniae samples isolated from a Chinese hospital.
In the case of genes involved in aminoglycoside resistance, the Aac6-Ib and armA genes had a frequency of 90% and 71%, respectively, while among the resistant samples, Aac6-Ib was the most frequently detected gene. Aligned with our results, Cirit et al., [50] also reported Aac6-Ib to have the highest frequency among the genes involved in aminoglycoside resistance in nosocomial K. pneumoniae isolates. A study of highly aminoglycoside-resistant K. pneumoniae and Klebsiella oxytoca from an inpatient in Okinawa, Japan, with no known history of travelling overseas was accomplished. Genome sequencing analysis showed that these isolates harbored armA, which encodes a 16S rRNA methylase, ArmA, that confers pan-aminoglycoside resistance [51]. In a study conducted by Li et al., [52] in China, out of 223 isolates of K. pneumoniae, 13 isolates (5.8%) contained armA and 8 isolates (3.6%) contained rmtB. 110 isolates of K. pneumoniae were phenotypically resistant and after PCR, 11.8% of the isolates contained armA gene and 7.3% of the strains contained rmtB gene. These values are lower than our findings about the armA gene.
According to the results of the modified Hodge test, out of 24 samples resistant to imipenem and meropenem, 15 (62.5%) samples were positive for carbapenems. Examination of virulence genes showed that among 84 samples, mrkD, rmpA, Irp-1 and magA genes showed the frequency of 94%, 83%, 60% and 11%, respectively; among 36 resistant samples, mrkD, rmpA, Irp-1 and magA genes have the highest frequency, respectively, while mrkD was significantly more abundant compared to the other genes. Based on these data, mrkD is evidently highly associated with multidrug resistance. In line with our findings, Liu et al., [53] identified mrkD as the most common virulence gene with the prevalence of 100%. In 2004, an initial study of the magA gene was performed by Fang et al., [54] who identified the gene as a virulence factor in the pathogenesis of K. pneumoniae. In this study, the magA gene was observed in 52 invasive strains (liver abscesses) and 15 non-invasive strains. El Fertas-Aissani et al., [55] obtained the opposite result from the present report. In this regard none of the studied strains carried the magA gene. They also examined the rmpA gene, which contained 3.7% of the 54 strains of K. pneumoniae isolated from different clinical specimens, which is a lower percentage than the present study. In another study, Liu et al., [48] reported that rmpA and magA were the most abundant genes among the 117 isolates of K. pneumoniae, respectively, and Fu et al., [56] observed the same frequency of genes in their study.
Since the presence of some virulence factors can be involved in the pathogenicity of bacteria, knowledge of the existence of these factors and their prevalence can be a good way to identify and treat the studied strain. The phylogenetic tree was drawn based on the UPGMA algorithm and the genetic relationship between the isolates was identified, which the isolates were differentiated into four clusters, G1- G4, with 70% similarity. Moreover, in the strains that were studied, the maximum number, 40, belonged to the G4 cluster, and the minimum, 8, belonged to the G3 cluster.
According to the research of Ferreira et al., [57] the dendrogram obtained from ERIC-PCR results showed a genetic relationship between 25 K. pneumoniae studied. In their study, the clusters were determined using the method (UPGMA) and dice’s similarity coefficient. Based on this, it was found that although the bacteria were isolated from different patients, K. pneumoniae in the bloodstream had a high genetic relationship with each other. In studies by Firmo et al., [34] ERIC-PCR analysis of 16 isolates of K. pneumoniae showed that only three isolates (18.8%) (K4-R3, K6-R3 and K7-R3) present a single band pattern, which is 100% genetically similar, i.e. there was a clonal relationship. The other 13 isolates of K. pneumoniae presented heterogeneous profiles with a maximum of 40% similarity. El-Badawy et al., [58] showed that most isolates have different origins by genotyping K. pneumoniae isolates using ERIC-PCR method. 32 isolates belonged to 18 different single roots, indicating that the prevalence of K. pneumoniae in different parts of the hospital was due to poor infection control.