3. Results
Class 1 integron was detected in 92% of isolates (46/50) using IntI1F and IntI1R primers (
Table 4). Screening of integron class 2 and 3 by detecting the presence of
IntI2, and
IntI3 genes, respectively, showed that these were not present in the collection (Gel electrophoresis
supplementary files). The antimicrobial susceptibility patterns of 50
K. pneumoniae isolates was performed (
Figure 1) . All clinical isolates were resistant to ampicillin (100%), most of the isolates were resistant to cefotaxime (98%), cefepime (92%) and ceftazidime (84%), and almost half of them showed resistance to piperacillin/tazobactam (48%). Moreover, 88% of the isolates showed resistance to ciprofloxacin (88%). Furthermore, these isolates showed intermediate resistance to gentamicin (54%), imipenem (52%), and a lesser extent to amikacin and colistin (30% and 6%), respectively.
According to the results obtained by phoenixBD semi-automated system, 54% of
K. pneumoniae isolates were identified as ESBLs, whereas 42 % (n=21) of isolates were identified as XDR and two isolates (4%) as PDR. All XDRs were susceptible to colistin except 3 isolates (Kp 2, Kp 22, and Kp 50) (
Table 5). In addition, 11.8% and 5.6% of the XDRs were susceptible to amikacin and gentamicin, respectively. Some of the XDR isolates were susceptible to amikacin, gentamicin, imipenem, and meropenem along with colistin (
Table 5).
The conjugation was performed between variable region-positive isolates as donors and
E. coli HB101 strain as a recipient. Seven representative samples were selected to be tested (Kp 5, Kp 21, Kp 27, Kp 37, Kp 42, Kp 49 and, Kp 50). The trans-conjugant PCR products are shown in (
Supplementry file). All samples were positive for
IntI1 gene, which confirmed the successful occurrence of conjugation via means of horizontal transfer. In the tested strains, the constituents of gene cassettes, plasmids, and sequencing types were analyzed by sequencing. WGS data showed that in our tested strains, plasmids were detected in all trans-conjugants which conforms that these conjugative plasmids are highly transmissible. The efficiency of trans-conjugation for each strain was calculated and ranged between 0.07% to 17.3% (
Figure 2.). This exemplifies the role of these plasmids as the main vehicle for the transmission of integrons, and subsequently the transport of various types of gene cassettes from the donor to the recipient strains.
Out of the total 50 isolates of K. pneumoniae, 24 isolates were analysed by WGS, including the 7 isolates which were successful in conjugation. Sixteen of these isolates were positive for the variable regions (Kp 5, 7, 10, 11, 15, 16, 21, 28, 37, 40, 42, 43, 44, 45, 46, 49) and 2 isolates showed bands with unexpected size for intI3 gene (Kp 4 and KP 22). Among the 24 samples, 3 were integron-negative isolates by PCR method (Kp 4, Kp 27, and Kp 30). However, WGS showed that two of these isolates were positive for class 1 integron (Kp 27 and Kp 30), whereas Kp 4 was still negative. Therefore, Kp 4 isolate was excluded from the subsequent molecular analysis.
MLST of all sequenced isolates was determined by the online MLST database from the Center of Genomics Epidemiology (CGE) [
25]. Nine of the isolates belong to sequence type (ST-231) (Kp 5, 6, 7, 10, 11, 15, 28, 30, and 45), six isolates were ST- 395 (Kp16, 22, 41, 43, 44, 50), and 2 isolates were ST-405 (Kp 42 and Kp 46). Only one isolate belongs to each of the following miscellaneous sequence types: ST- 37 (Kp 40), ST-45 (Kp 21), ST-147 (Kp 49), ST-280 (Kp 27), ST-1710 (Kp37) and ST-1741(Kp 25).
WGS data of 23 integron-positive isolates were used to construct a whole genome phylogenetic single nucleotide polymorphism (SNP) (
Figure 3), and the output from the analysis showed three main clusters of strains. The first cluster is composed of 9 isolates belonging to ST-231. The second cluster of strains is in the ST-395 group (n=6). The clinical data of the above mentioned two STs, indicate frequent moving of patients from Intensive Care Unit (ICU) to the male medical ward or vice versa. However, the third cluster consisted of two strains only belonging to ST-405. The remaining strains belong to miscellaneous sequence types (ST-37, ST-147, ST-280, ST-1741, ST-45, ST-13, ST-17, and ST-10).
A close-up analysis of the two predominant STs (ST-231 and ST-395) is done to predict if there any possible recent outbreak through clustering (Figure 4.). SNP tree for ST-231 and ST-395 individually was created and it was observed that the PDR isolates (Kp6 and Kp11) were branching from the same clade suggesting possible vertical transmission among patients. Similarly, among the other ST-231 and ST-395 stains are very colsely related.
Figure 4.
SNPs trees of K. pneumoniae strains (ST-231 and ST-395).
Figure 4.
SNPs trees of K. pneumoniae strains (ST-231 and ST-395).
A phylogeny tree showing relatedness of K. pneumoniae strains that belonged to A) ST-231 (n=9) and, B) ST-395 (n=6) . Isolates are labeled according to their numbers and sequence type (ST). Green: ST-231 and red:ST-395. K. pneumoniae HS11286 was used as a reference isolate (accession number CP003200). The online tool (iTOL) was used to draw and edit the tree.
The demographic data for thirteen samples that belong to ST-231 (n=9) and ST-395 (n=4) which accounts for most samples were list in
Table 6. Ninety-two percent of the patients were males with age more than 50 years (69.2%). The length of stay varied between patients with a range from 1 to 134 days. Most of the patients were admitted to the ICU (61.5%) and the rest of the patients were admitted to intermediate care wards (male medical wards). Highly critical ICU patients usually require broad spectrum antibiotics which exterts tremendous selection pressure, thus driving antimicrobial-resistance bacteria to thrive further in these settings.
At the level of ST, more than two-third (77.8%) of ST-231 and all ST-395 strains were isolated from patients admitted either in ICU or male medical wards or moved between the two places interchangeably. Almost all patients have various risk factors contributing to their likelihood of acquiring highly resistant
K. pneumoniae strains including ICU admission, ventilation, urinary catheters, central venous catheters, and hospitalization for more than 7 days [
28]. Unfortunately, the mortality rate among these patients is very high (84.6%).
All sequenced isolates carried a wide range of acquired antimicrobial resistance-conferring genes. 14 isolates harbored
blaOXA-232 gene, all of which were XDR and PDR isolates except KP 7 which was the only XDR isolate that did not carry
blaOXA-232 gene (
Figure 3). The
bla-NDM5 was detected in only one isolate (Kp 49). Based on the WGS data and Carba-R, 87% of XDR and PDR isolates (20/23) have
blaOXA-48-like, 8.7% (2/23) have
blaNDM, and one isolate (4.3%) has both
blaNDM and
blaOXA48-like genes. There was a high level of agreement between WGS and GeneXpert results for the tested samples (8/9). The
blaOXA-1 allele was detected in eight isolates, most of which belong to ST-395 (n=6). The ESBL gene
blaCTX-M-15 was detected in all isolates, while
blaTEM-1B was seen in 15 isolates. Moreover, 10 isolates were positive for
blaSHV-1, 9 isolates for
blaSHV-11, 2 for
blaSHV-76, and 2 for
blaSHV-27.
A number of various genes conferring resistance against aminoglycosides were detected including aac(6')-Ib-cr, aac(6')-Ib Hangzhou, aadA1,2, aac(3)-lId, aph(3')-Ia, strA, strB, armA and rmtB. Fosfomycin-resistance conferring gene (fosA6) was detected in all isolates except Kp 49 which harbored fosA5. In addition, 6 different genes encoding quinolone resistance were detected in the isolates which are aac(6')-Ib-cr, oqxA, oqxB, qnrS1, B1and B66. Moreover, rifampicin (arr-2) and trimethoprim resistance encoding genes (dfr) were found in 10 and 11 isolates, respectively.
Tetracycline resistance gene, tetA, was detected in four ESBL isolates, whereas tetD in 2 XDR isolates. According to the phoenix automated system, all of these were tigecycline-resistant except one ESBL isolate that carried the tetA gene (Kp 46). The gene of acrAB efflux pumps’ regulator (marA), which contributes to the resistance against tigecycline and other antibiotic classes, was detected in all sequenced isolates. In addition, K. pneumoniae efflux pumps (Kpn E, F, G, H) that confer resistance to different antibiotic classes like macrolide, aminoglycoside, cephalosporin, tetracycline, rifamycin, and colistin, were detected in all isolates.
The four colistin-resistance isolates belonging to ST-231 and ST-395 (Kp 6, Kp 11, Kp 22, and Kp 50) were negative for plasmid-mediated
mcr genes, therefore, further analysis of the SNPs in PhoPQ operons was performed. No chromosomal mutations were detected in the regulatory two-component systems (TCSs) PmrAB and crrAB, and in
mgrB gene, a negative regulator gene of TCSs. However, all isolates harbored the regulatory TCSs (PhoPQ). Our analysis of the PhoPQ operons in the colistin-resistant isolates identified four novel undescribed SNP in the PhoP genes with the amino acid substitutions as follows: (Val129Glu), (Gln147His), (Gln131Glu), and (Pro129Thr) (
Table 7.). The significance of these SNPs is yet to be determined as none of these substitutions have been described.
Most of the tested isolates showed an agreement between their resistance phenotypic patterns to different antibiotic families and the presence of resistance genes as shown in
Table 8. Carbapenems, cephalosporin, and quinolone antibiotics showed a complete level of agreement between phenotype and genotype. However, aminoglycoside antibiotics (AK, CN) showed a low level of agreement with 21 isolates carrying the resistance genes and only 9 and 15 isolates with phenotypic resistance patterns to amikacin and gentamicin, respectively. This observation could be because these resistance-conferring genes might not be expressed. The expression of these genes needs to be further investigated which is not within the scope of this study.
Using WGS data, the variable region of 23 samples including the trans-conjugant strains (Kp 5, Kp 21, Kp 27, Kp 37, Kp 42, Kp 49 and Kp 50) were analyzed to detect gene cassettes (
Table 9 and Figure 3.8). Ten different gene cassettes were identified in 16 strains with the varible region including those encoding resistance to aminoglycoside (
aadA, aacA4, aac(6')-Ib ,
APH(3”)-Ia, ant1) trimethoprim (
dfrA5, dfrA12, dfrA14), rifampin (
arr2 and arr3), chloramphenicol (
catB3, catB8 and cmlA1), macrolide (
ereA2 and erm ) and quaternary ammonium compound
(emrE). The most frequent gene cassette was composed of
dfrA14 gene alone, which was present in 8 isolates (34.8%). However, 3 samples showed empty integrons (In0) without any gene cassette insertions (Kp 16, Kp 22, and Kp 25). One of them had a variable region (Kp 16), while the other two did not (Kp 22 and Kp 25).
Comparative analysis of integrons was performed by locating the gene cassettes within the variable region. Mapping was done based on the genes between the primers of the variable region which varied among strains (
Table 9). It was observed that there was a contig break in the middle of the integron cassettes in some strains and the downstream sequence (3’-end) was found in a separate contig. Therefore, the genetic maps of the integrons were drawn after manual reassembly of the various components (
Figure 3).
Figure 4.
Schematic map of the composition of the integrons in K. pneumoniae strains. Genes and ORFs are denoted by arrowhead indicating the direction of transpcription and colored based on the gene function classifications as shown in the key in the lower part of the figure. Shaded areas denote regions of homology (>95% nucleotide sequence identity). Grouping of strains is based on the sequence similarity as follows: ST-231 strains (Kp 41, Kp 43, Kp 44, and Kp 50), (Kp 27, Kp 37, Kp 42, and Kp 46), (Kp 40, and Kp 21), (Kp 5, and Kp 15) and (Kp 49, Kp 28, and Kp10).
Figure 4.
Schematic map of the composition of the integrons in K. pneumoniae strains. Genes and ORFs are denoted by arrowhead indicating the direction of transpcription and colored based on the gene function classifications as shown in the key in the lower part of the figure. Shaded areas denote regions of homology (>95% nucleotide sequence identity). Grouping of strains is based on the sequence similarity as follows: ST-231 strains (Kp 41, Kp 43, Kp 44, and Kp 50), (Kp 27, Kp 37, Kp 42, and Kp 46), (Kp 40, and Kp 21), (Kp 5, and Kp 15) and (Kp 49, Kp 28, and Kp10).
Among the detected plasmids, eight of them harbored antibiotic resistance genes (Figure 4. And Table 10). The most prevalent plasmid is pKPQIL-IT, which is present in 15 samples that belong to ST-395 and ST-231. It is a 115,300 bp in size IncFIB(QIL) replicon carrying genes associated with resistance to β-lactams (blaTEM-1, and blaKPC-3 ).
Figure 4.
Bar chart representing the distribution of plasmid replicon types identified among K. pneumoniae (n=23). Plasmid types were detected based on WGS data by using ResFinder online tool.
Figure 4.
Bar chart representing the distribution of plasmid replicon types identified among K. pneumoniae (n=23). Plasmid types were detected based on WGS data by using ResFinder online tool.
Plasmid pKP3-A is a 7,605 bp linear ColKP3 replicon, carrying blaOXA-181 gene. ST-231 and ST-405 sharing pKP3-A uniqueness that the integron is interrupted by insertion sequences (Figure 3.8). These isolates exhibited highly similar genotypes. All produced OXA-181, and the majority also have plasmid mediated- blaTEM-1, blaKPC-3 genes. Unlike Xpert Carba-R, the plasmid finder tool detected kpc gene from WGS data in all ST-231 and ST-405 isolates. This could suggest that the copy numbers of the plasmid might be low or the gene might have not been expressed and hence could not be detected by the GeneXpert system.
Isolates Kp 41 and Kp 50 are the only two strains carrying IncFIB (pNDM-Mar) and IncHI1B (pNDM-Mar) plasmids. In this isolate, blaNDM-1 was within pNDM-Mar plasmid. Whereas, Isolate Kp 49 (ST-147) had blaNDM-1 as part of (ble-blaNMD-1) operon where it was flanked by bleomycin resistance (ble) gene and N-(5'-phosphoribosyl) anthranilate isomerase (trpF) gene. Moreover, the plasmid composition of Kp 49 had IncR, IncFII, and pKP3-A. The plasmid pAMA1167-NDM-5 is present only in 7 of the ST-231 isolates. It is a 11,310 bp size IncFII(pAMA1167-NDM-5) replicon with genes encoding resistance against aminoglycosides [aph(3’’)- Ib, aph(6)-Id, aadA2, aadA5, aac(3)-IIa, aac(6’)-Ib-cr5 ], β-lactams (blaNDM-5, blaOXA-1, blaCTX-M-15-1, blaTEM-1), chloramphenicol (cat), sulfonamides (Sul1 and Sul2), trimethoprim (dfrA12 and dfrA17), tetracycline (tet(b) and tet(c) ), and macrolides (emrE and mp(A)).Sequence analysis of pKPN-IT plasmid revealed that it is a 208,191-bp IncFIB(K) replicon carrying (aadA2, cat, Mph(A), Sul1, dfrA12) genes. It is found in KP 10, 25, 27, 37, 42, and 46 isolates. However, plasmid pBK30683 is a 139,941-bp FIA replicon present only in one isolate (KP 37).
Table 11.
plasmids existent in K. pneumoniae isolates and there features.
Table 11.
plasmids existent in K. pneumoniae isolates and there features.
isolate (KP) |
Plasmid |
Size (bp) |
Replicon |
Resistance genes |
5, 6, 7, 10, 11, 15, 16, 22, 28, 30, 41, 43, 44, 45, 50 |
pKPQIL-IT |
115300 |
IncFIB(QIL) |
blaTEM-1, blaKPC-3
|
5, 6, 11, 15, 16, 22, 28, 30, 41, 43, 44, 45, 49, 50 |
pKP3-A |
7605 |
ColKP3 |
BlaOXA-181
|
5, 6, 7, 15, 28, 30, 45 |
pAMA1167-NDM-5 |
11310 |
IncFII(pAMA1167-NDM-5) |
aadA5, aadA2, aac(3)-IId, aph(6)-Id, aph(3’’)-Ib, aac(6’)-Ib-cr5, blaNDM-5,blaOXA-1, blaCTX-M-15-1, blaTEM-1 dfrA17, dfrA12, Mph(A), Sul1, Sul2, emrE, tet(b), tet(C), cat
|
10, 25, 27, 37, 42, 46 |
pKPN-IT |
208191 |
IncFIB(K) |
aadA2, cat, Mph(A), Sul1, dfrA12 |
21, 40 |
pCAV1099-14 |
113992 |
IncFIB(K)(pCAV1099-14) |
dfrA19, APH(3’’)-Ib, APH(3’)-Ia, QnrB52 |
41, 50 |
pNDM-MAR |
267242 |
IncFIB (pNDM-Mar) IncHI1B (pNDM-Mar) |
aac(6’)-Ib, blaOXA-1 , blaNDM-1, cat, QnrB1 |
27, 49 |
pK245 |
98264 |
IncR |
aacC2, strA, strB, dfrA14, catA2, Qnrs, blaSH2A , blaTEM |
49 |
pC15-1a |
92353 |
IncFII |
aac(6’)-Ib, aac(3)-II, blaTEM-1 , blaOXA-1, blaCTX-M-15-1 , tet(A)
|
37 |
pBK30683
|
139941 |
FIA(pBK30683) |
dfrA14, StrA, StrB, blaTEM-1, blaOXA-9, , blaKPC-3 Sul2,ant(3’’)-Ia
|
4. Discussion
The present study focused on molecular characterizations of integrons as a common class of mobile elements, and their significance in the dissemination of multi-drug resistance genes among
K. pneumoniae. Whole genome sequence (WGS) data were analyzed to determine the sequence type as well as charecterisation of integron-carrying plasmids. In addition, WGS was used for studying the correlation between antimicrobial resistance genotype and phenotype, and the level of agreement was analyzed. Moreover, the phylogenetic relatedness of isolates was associated with the patient’s demographic data to explore the possible spread of mobile elements. Most of the studies in the region investigate the epidemiology of AMR with a few studies concerning the prevalence of mobile genetic elements such as genomic islands and their molecular characteristics in Gram negative bacteria (GNB) such as
A. baumannii and
K. pneumoniae [
16,
17,
29,
30,
31,
32,
33]. Moreover, another study was conducted in Gram-negative bacteria isolates in Palestinian hospitals focused on integrons and their role in AMR dissemination [
34]. To our knowledge, this is the first study in Oman that assesses that
transmissibility of plasmids by examining the efficiency of conjugation
in-vitro thus emphasizing the importance of this mechanism in the spontaneous transfer of integron-carried plasmids intra-specially.
In the current study, 96% (48/50) of our
K. pneumoniae isolates carried class 1 integron, whereas none of the isolates contained either class 2 or class 3 integrons. This finding is similar to the results of previous studies that were conducted in Iran [
35,
36]. In addition, a previous study showed that 100% of the MDR isolates were found to be positive for class 1 integron, 36% for class 2, and none for class 3 integron [
10]. However, other studies detected class 2 integron at a low prevalence level of 1.7% and 8.3% [
27,
34].
High prevalence of class 1 integron among our MDR isolates was observed and therefore it is the main focus of this study. The presence of class 1 integrons in most of the isolates suggest that these genetic elements confer a significant advantage to their hosts where exposure to antibiotic overuse creates a selective pressure in hospital environments [
4,
10]. The absence of both class 2 and class 3 integrons in our isolates is expected since these classes are relatively uncommon [
6,
37,
38,
39]. Moreover, this might be due to the fact that our selection is limited to representative MDR isolates which were unusually difficult to treat over a short period of time, which might not reflect the actual prevalence. Therefore, a larger sample size including screening specimens is more likely to reveal the other classes of integrons.
In our study, all isolates were resistant to ampicillin which corresponds to the report of previous studies [
2,
40,
41]. Moreover, our findings showed a higher prevalence of carbapenem resistance (imipenem and meropenem) than other studies [
1,
16]. However, the resistance rate to aminoglycosides in our study was lower than previous reports where amikacin and gentamicin were resistant to 51.6% and 45.2%, respectively (Derakhshan et al., 2013). These results were expected since our isolates were selected as MDRO from the beginning of this study, whereas in the above-mentioned studies, isolates were chosen regardless of their antibiotic resistance profile and were mostly screening samples. This can explain the difference in the antimicrobial susceptibility profile of our isolates against some of the classes of antibiotics.
Our findings showed that 37.5% (18/48) of integron-positive isolates were carrying variable regions. By the analysis of WGS data, a total of 10 different gene cassettes were detected in these isolates. The latter isolates encode for aminoglycoside, trimethoprim, rifampin, chloramphenicol, macrolide, and quaternary ammonium compounds. Trimethoprim resistant gene cassettes (
dfrA14),
which encodes for dihydrofolate reductase enzymes were found to be predominant in our isolates (34.7%). A number of studies showed that
dfr genes were the most frequent gene cassettes present in integrons [
2,
14,
34,
42]. The stability of this gene cassette in class 1 integrons might indicate the overuse of the trimethoprim for a long period of time specially to treat urinary tract infections [
14]. Moreover, despite the initial results of the variable region which was detected by employing conventional PCR using primers that bind in the 5’ and 3’ conserved regions, further WGS analysis showed that 13% of isolates with positive variable regions carried empty gene cassettes. Moreover, previous studies reported the presence of
sul1 at the 3’CS region of the integrons [
43,
44,
45]. However, it was only observed in (3/20) isolates, despite the amplification of 5’CS-3’CS region by PCR. This has been described in
Salmonella enterica [
46,
47]. The non-classical structure of integrons where
intI1 and/or 5’ end of integrons is truncated were also observed in two strains (Kp 16 and Kp 25), which might suggest that amelioration of
IntI1 confers a selective advantage for the host bacterium.
Simple PCR assays tend to overestimate the presence of integron cassettes, and cannot provide further information on the subsequent changes on the integron cassettes such as deletions and insertions. In fact, gene cassettes undergo continuous changes in their compositions which means that these cassettes can integrate antimicrobial resistant genes once conditions are favorable. Truncations and inversions within the integron structures are not uncommon, and could not be further detected by simple PCR. Therefore, WGS data were key in bridging the gap on the various components within these integrons.
The gene cassette arrays (dfrA12-ant1) and (dfra12-emrE-ant1) counted for 8.7%. Kp 49 was the only isolate that harbored a cassette composed of four genes (Arr3- ereA2- aadA- cmlA1). In addition, some of the gene cassettes like aac, cat, emerE, erm, arr, and ant occurred in combination with different cassettes, indicating that new gene cassettes can still be generated. In this study, an attempt was made to classify the integrons by looking for the integron types in the INTEGRALL database, and it was found that some isolates were assigned to the same number of gene cassettes in INTEGRALL despite having different compositions. Upon further manual analysis, it was found that some of the old deposited sequences in the INTEGRALL database have different annotations for some genes. To avoid confusion, INEGRALL classification numbers were omitted in our analysis.
Plasmids are considered as main carriers for antibiotic resistant genes through HGT mechanisms. Spontaneous plasmid conjugation is possible when the plasmid has a compatible mechanism of transfer. A plasmid is classified as conjugative when it helps other conjugative elements like integrons to move between different strains [
48]. In the current study, the trans-conjugant colonies were positive for
intI1 gene in all the representative samples with trans-conjugation efficiency ranged between 0.05% to 17.3%. This confirmed the successful transfer of integrons between cells via HGT. WGS data for the tested isolates supports our hypothesis, as plasmids were identified in all of the tested samples which is similar to the result of a previous study [
14]. Being carried on plasmids, these integrons might have a significant fitness cost since more copies of plasmids can be present compared to chromosomally located integrons. This observation has been reported in previous studies on
K. pneumoniae as well as other Gram negative organisms [
49,
50].
Most of our sequenced isolates showed a hierarchical relationship between different types of MGEs. The detected gene cassettes were embedded inside integrons which in turn inserted inside transposons which are flanked by insertion sequences (IS) that can confer mobility to transposons. Some of these composite transposons were incorporated into conjugative plasmids that offer transfer to other cells [
51]. Most of the isolates in this study (Kp 5, 21, 27, 28, 37, 40, 41, 42, 44, 46, 49 and, Kp 50) have the integrons inserted in transposons that were flanked by IS (Figure 4.). All above-mentioned strains had similar IS which is IS
6 and IS
6100 (Kp 42 and Kp 46). IS
6100 belongs to the IS
6 family of transposable elements forming co-integrates as an endpoint of transposition, which was originally isolated from
Mycobacterium fortuitum [
52]. In general, different previous studies showed the same relationship, were the integrons were located in transposons that imbedded in a conjugative plasmid [
53,
54]. However, MGEs including Tns and ISs had disrupted the integron structure. IS are usually small MGEs which are ubiquitous within most bacterial genomes. These IS act as facilitators for horizontal gene transfer and play a major role in the evolution of microbes thus increasing their diversity and adaptive capacity.
The presence of IS
6 on both sides of the Dcm-methylation operon which is a type II Restriction-Modification (RM) System and next to the integron denotes that RM systems are mobile and located in plasmids. These systems are involved in genome rearrangements and enhance virulence and resistance plasmids dissemination by carriage on other MGEs like in our study on integrons [
55,
56,
57]. The existence of these IS
6 in our isolates entails their importance in disseminating antibiotic resistance genes among various genera and species of bacteria [
58,
59]. In this study
, the specific capture of resistance cassettes by the integrons signifies apparent bias driven by the selective pressure of antibiotic therapy regimes. Other studies have also shown that integron could carry catabolic genes which prove that integrons play a broader role in bacterial evolution [
8,
60].
According to the MLST results in this study, 23 isolates were classified into nine different STs. The most common sequence type was ST-231(39%) followed by ST-395 (26%). Moreover, two isolates belonged to ST-405 was observed in patients from the same ward which is also suggestive of possible spread, while only one isolate in each of the following STs: 37, 45, 147, 280, 1710, and 1741 which were reported in previous studies [
61,
62,
63,
64].
ST-231; the most abundant STs found in this study are widely distributed in South-West Asia with clonal dissemination in Singapore, Brunei, and Darussalam between 2013 and 2015 [
65,
66]. In addition, in India, ST-231 strains were reported as a predominant ST in 22 isolates (45%) [
67]. In Europe, the first occurrence of the MDR
K. pneumoniae ST-231 clone was confirmed in Switzerland [
68]. This dissemination may represent a global public threat toward a new epidemic clone. ST-395 has been reported as the most common MDR
K. pneumoniae clone (69%) in North-Eastern France as well as an outbreak in ICU in Italy [
69,
70].
In a previous study conducted in isolates from Arabian Peninsula, different ST types were detected in MDR
K. pneumoniae in gulf countries including Oman, however, ST-231 was not among these STs [
71]. Similarly, previous studies were conducted in Oman and Saudi Arabia which have concluded that no ST-231 strain was found [
72,
73]. The absence of ST-231 possibly indicates a recent emergence of ST-231 in Oman. However, ST-147 was one of the commonly detected STs from
K. pneumoniae isolates in the two previous studies [
71,
72], whereas in the present study only one strain belongs to ST-147.
Both of these STs, ST-395 carrying strains and 77.8% of ST-231, were clustered in the same wards (mainly medical and ICU wards) with a significant mortality rate (84.6%). This observation highly suggests that there is a horizontal transfer of resistance-conferring genes, which is alarming. The heavy use of broad-spectrum antibiotics for critically-ill patients in high dependency areas and ICU patients selects for highly resistant strains and enhances the spread of resistance determinants. Furthermore, this may indicate inefficient infection control practices that allowed the dissemination of the plasmid-mediated resistance in the hospital. This finding corresponds to a study that was conducted in carbapenem-resistant GNB in SQUH, where 87% of infections were being healthcare-associated with a 62% mortality rate [
17].
Over the last 15 to 20 years, carbapenem-hydrolyzing β-lactamases including OXA-48-like and NDM type carbapenemases disseminated and emerged in Enterobacteriaceae all over the world [
74,
75]. Carbapenemase OXA-48-like differs from the classical OXA-48 by one to five amino acids in which they hydrolyze carbapenems and penicillins but do not affect extended-spectrum cephalosporin [
75].
OXA-232-producing
K. pneumoniae (OXA232Kp) was identified for the first time on a 6.1-kb ColE-type non-conjugative plasmid in France in 2013 from a patient who returned from India and since then it has spread worldwide [
74,
75,
76]. OXA-232 is considered as a point mutation derivative from OXA-181 with one amino acid difference [
77]. In this study, WGS and Xpert Carba-R data detected the presence of the
blaOXA-48-like gene in 87% (20/23) of our XDR and PDR isolates. These genes were identified by WGS as
blaOXA-232 in 93.3% (14/15) of our XDR and PDR sequenced isolates. All of them belonged to ST-231 and ST-395. The level of agreement between WGS and GeneXpert tested samples was highly significant (n=8/9) with only one isolate that was negative for OXA gene at GeneXpert and positive in WGS (Kp 43). Illumina WGS has good coverage for sequencing the whole genome, whereas, GeneXpert is based on real-time PCR, even with 97% specificity might miss some allelic variants. Plasmid analysis of our strains showed that the OXA-181 gene (or OXA-232) was carried on pKP3-A (7605 bp) replicon type ColKP3. It was expressed only in the XDR and PDR OXA-232 positive strains. This is similar to a previous study where 33% and 100% of OXA-232 isolates were carried on ColKP3 plasmid in the USA and China respectively [
77,
78].
On the other hand, NDM was first detected in 2008 from a Swedish patient of Indian origin and then distributed all over the world [
79]. In this study, Xpert Carba-R results showed that only 4% (2/50) of our isolates had the NDM gene. Whereas, only one isolate (4.3%) that belongs to ST-147 harbor both OXA-48-like (
blaOXA-232) and NDM gene (
blaNDM-5) according to WGS data. A study was conducted in UAE to characterize carbapenem-resistant Enterobacteriaceae in the Arabian Peninsula, showed that NDM and OXA-48-like are the most commonly detected genes with a rate of 46.5% and 32.5%, respectively. In the samples collected from Oman, 46% and 44.4% of the isolates were carrying NDM and OXA-48-like (
blaOXA-181), whereas, only 1.6% carry both genes [
71].
Along with the carbapenem genes, ESBL genes including
bla(CTX-M),
bla(TEM), and
bla(SHV) were detected, and have been reported in different studies [
66,
67]. All
bla(CTX-M) genes were identified including
blaCTX-M-15, and it was the dominant ESBL amongst our strains since it was produced by all sequenced isolates. One main factor that contributes to the CTX-M-type ESBLs distribution in
K. pneumoniae is the conjugative plasmids, particularly that belong to IncF with specific insertion sequences [
80]. Moreover, all bla(TEM) genes belonged to blaTEM-1 and detected in 69.6%, whereas, the SHV types were characterized as
blaSHV-1, blaSHV-11, blaSHV-27, and blaSHV-76. All ST-231 strains carrying
blaSHV-1, while ST-395 strains carrying
blaSHV-11. Besides, OXA-1 was positive in all ST-395 strains only with an association with other ESBLs genes
(blaSHV-11, blaCTX-M-15, blaTEM-1). The association between OXA-1 and CTX-M-15 genes will render isolates resistant to β-lactam-β-lactamase inhibitor combinations. Note worthily, a large variety of GNB carried OXA-1 gene in plasmid and integron locations [
81]. Plasmid analysis in this study found that
blaOXA-1 gene was located along with
blaCTX-M-15 in IncF plasmids IncFII and IncFII (pAMA1167-NDM-5).
Resistance to colistin, a drug that is used as the last line drug in the treatment of extensively resistant pathogens, is increasingly reported in Enterobacteriaceae, particularly
K. pneumoniae [
82]. Several mechanisms are involved in colistin resistance, mostly the emergence of mobilized colistin resistance (
mcr) genes via plasmid and the mutations in the chromosomal gene (
mgrB) and operons (PmrAB, and PhoPQ) which are associated with the biosynthesis and modification of lipopolysaccharide (LPS). The PhoPQ regulatory system is activated at a low concentration of Mg2+ or Ca2+ and acidic PH [
83]. The PhoQ (sensor kinase) activates PhoP (regulator protein) by phosphorylation, which in turn activates
pmrFHIJKLM operons. These operons lead to LPS modification by adding 4-amino-4-deoxy-L-arabinose (L-Ara4N) and phosphoethanolamine (PETN) to lipid A. This modification neutralizes the negative charge of LPS leading to the low affinity of LPS to positively charged colistin [
83,
84,
85].
Four sequenced isolates were resistant to colistin via BMD, however, all of these strains were negative for the plasmid-mediated colistin resistance gene
(mcr) and harbored a truncated
mgrB gene. Bioinformatics analysis revealed the presence of four SNPs in the PhoP genes with the following amino acid substitutions: (Val77Glu), (Gln147His), (Gln131Glu), and (Pro129Thr). Further analysis of the promoter region of PhoPQ operon did not reveal any SNPs. A previous study reported different SNPs in PhoP (Val3Phe and Ser86Leu) and PhoQ
(Leu26Pro) that were found to have a role in colistin resistance when tested by real-time PCR [
86]. However, the detected SNPs in this study have not been described before, and therefore, further expression analysis by real-time PCR is needed to show their role a long with the truncated
mgrB gene in colistin resistance. Moreover, colistin resistance might arise
due to other resistance mechanisms such as capsule overproduction that causes a reduction in the interaction between colistin and its target site at
K. pneumoniae (LPS) [
83,
87]. Furthermore, all isolates in this study carried
K. pneumoniae efflux pumps (Kpn E, F, G, H) which might contribute to colistin resistance.
Tigecycline has broad-spectrum activity against Gram-positive and Gram-negative bacteria and has demonstrated efficacy and safety as salvage therapy for MDR/XDR bacteria [
29,
32,
33,
88,
89]. However, the tigecycline resistance is becoming more common in ESBL-producing, MDR, XDR, and carbapenem-resistant isolates [
90,
91]. In this study, the majority of the sequenced isolates (17/19) were tigecycline-resistant, while the remaining four isolates were not tested. The WGS data showed that the efflux pumps’ regulator gene (
marA) was present in all isolates (n=23), including the two sensitive ones (Kp37 and Kp46). In this context, a study reported that the
marA gene, along with other regulatory pathways, could mediate the resistance to tigecycline through upregulating
acrAB efflux pump even in the absence of
ram A [
92]. In another study, the development of tigecycline resistance was attributed to the
tetA gene [
93]. In our study, the
tetA gene was detected in four ESBL isolates, from which only two were resistant
. Thus, it’s more likely that
marA regulator is engaged in the tigecycline resistance.
Different genes that confer resistance against aminoglycosides were detected including
aac(6')-Ib-cr, aac(6')-Ib Hangzhou,
aadA1,2,
aac(3)-lId,
aph(3')-Ia,
strA, strB, armA and rmtB. The latter two genes are plasmid-encoded 16S rRNA methylase genes encoding for enzymes that have the ability to alter the antibiotic target [
94]. Several studies were reported that
armA gene was located in conjugative plasmids of the IncL/M or IncFIA plasmid, while
rmtB in IncFI or IncFIA [
94,
95]. Interestingly, these two genes were only identified in Kp 49, which belongs to ST-147, and it is OXA-232, and NDM-5 producing isolate. In addition, IncFIA(HI1) and IncFII plasmids were uniquely present in this isolate only, and a further search of
rmtB and
armA genes by BLAST using WGS data confirmed the localization of these two genes in the (IncFII) plasmid.
Among carbapenems, cephalosporin, and quinolones, all isolates possessed expressed genes with a complete agreement between AMR genotypic findings and phenotypic expression. These results are close to what has been detected in a previous study [
96]. In contrast, aminoglycosides antibiotics (Amikacin and gentamicin) had high non-expressed genes (57.1% and 28.6%). Since there are other mechanisms behind AMR, the presence or absence of certain genes might not be a sufficient indicator of the isolate’s resistance profile. Furthermore, genotypic analysis has to be correlated with phenotypic findings of various mechanisms involved.
The versatility and ubiquity of integrons in bacterial genomes indicate the key role of these mobile elements in bacterial adaptation. The abundance of repeat sequences of ISs and Tns within the integron structures in our WGS data pertains to the tremendous selective pressure of used antibiotic regimens in the hospital microbial environment. These integrons are very potent capture systems with a limitless capacity to exchange antibiotic gene cassettes as well as other genes that increase their fitness. Integron with combinations of various antibiotic cassettes had been termed as multi-resistance integrons (MRIs) in Gram negative bacteria [
97]
.
The genomic plasticity was apparent in a few cases where there was a deletion of non-essential genes or empty integrons which confirms that these integron systems are performing adaptive rather than housekeeping functions. In support with this hypothesis, two of the isolates we examined had integrons with no antibiotic gene cassettes and the intI1 gene was partially truncated (Kp 16 and Kp25). The different attachment sites and the variety of open reading frames of hypothetical or unknown functions within these integron systems are compelling. Further detailed analysis is required to unravel the importance of integrons among other MGEs towards the adaptation of bacterial evolution.