1. Introduction
The rapid emergence and spread of Antibiotic-Resistant Bacteria (ARB) is a global health concern. Particularly, Carbapenem-Resistant Bacteria (CRB), according to the recent World Health Organization (WHO) report, are an alarming threat to public health worldwide [1]. Resistance against carbapenem, the last-resort antibiotic, in the Enterobacteriaceae and Acinetobacter baumannii complex, has progressively worsened and led to the development of Multi-Drug Resistance (MDR) phenotypes. This is a major challenge for health institutions around the global due to limited antibiotic choices and high mortality rates [2]. In Pakistan, the over-the-counter availability of antibiotics promotes misuse among humans and animals. The use of antibiotic in food-producing animals as prophylactic and growth-promoting agents is poorly controlled [3]. Pharmaceutical manufacturing and agricultural runoffs also contribute to environmental contamination and the persistence of antibiotic compounds in the environment [4]. Inadequate treatment of antibiotic-laden wastewater facilitates their dispersion into aquatic environments such as streams, canals, and rivers, exerting continuous selection pressure on both pathogenic and non-pathogenic bacterial populations. This pressure promotes the transmission of antibiotic resistance genes via plasmids or other mobile genetic elements (MGEs) among different species and even genera. Consequently, this significantly contributes to the emergence and dissemination of Antibiotic-Resistant Genes (ARGs) in diverse environmental settings [5,6]. ARB can potentially infiltrate into drinking water supplies, especially in urban areas of low-resource countries like Pakistan, as reported in a similar developing country, Bangladesh, where the researcher isolated ARB from drinking water facilities [7]. In such cases, wastewater acts as a crucial vector for the transmission of bacterial pathogens to humans in regions with inefficient water quality management systems. Factors such as larger populations, rapid industrialization, and exposure to climatic changes, as seen in Pakistan, pose serious risks to public health and individuals interacting with these environments [8,9].
Integrons are genetic elements identified to facilitate the horizontal gene transfer between bacteria by incorporating exogenous genes [10]. The genetic structure of an integron is composed of the integron integrase (IntI), which is a site-specific recombinase, a recombination site (attI), and a promoter (Pc) that directs the transcription of foreign genes. The IntI mediates the integration that captures, removes, and shuffles circularized genetic components called gene cassettes. These gene cassettes, consisting of a cassette recombination site (attC) and an open reading frame (ORF), can lead to diverse adaptive functions [11]. Integrating new gene cassettes, which are diverse and abundant in the environment, can lead to rapid adaptation and help the bacteria evolve in specific traits. The integrons were first identified due to their significant contribution in the spread of antibiotic resistance genes. Recently, it has been revealed that they can encode gene cassettes, which are of an unknown function but often encompass proteins involved in resistance to disinfectants, heavy metals, cell interactions, and niche adaptation [12]. There are three integron classes based on the amino acid sequences of IntI integrases as (IntI1, IntI2, and IntI3) that are known to be responsible for resistance to antibiotics. The class 1 integrons are the most predominant among bacteria and are mostly found in the Enterobacteriaceae family as well as other clinically significant Gram-negative bacteria [13].
In this study, we investigated thirty-two bacterial isolates from wastewater sources in Pakistan and determined the phenotypic antibiotic resistance profiles. Whole-genome sequencing was performed to identify the ARGs, point mutations responsible for resistance, MLST, plasmid replicons, and integron gene cassettes. Additionally, Oxford nanopore technology was used to sequence one representative isolate to detect the plasmid and the genetic environment of ARGs using the combined Illumina and nanopore reads for that isolate. Although extensive studies have reported antibiotic resistance in bacterial isolates from environmental water sources globally [14,15,16], in Pakistan the majority of the studies focused on the detection of ARGs using the PCR-based amplification approaches from wastewater bacterial isolates [17,18,19,20]. To the best of our knowledge, there are limited studies published regarding the whole-genome sequencing analysis of wastewater samples [21,22]. There are reports regarding the whole-genome sequencing of CRB, but these are from humans, hospital settings (including surfaces, intensive care units, and sinks), animal sources, and food products—not wastewater sources [23,24,25,26,27,28,29]. However, the characterization of ARB and particularly the CRB isolates from wastewater using whole-genome sequencing analysis in Pakistan needs attention because there is little evidence concerning the diversity, spread, and resistome content of the wastewater-borne bacterial population. This investigation aimed to evaluate the prevalence of ARB and CRB in wastewater samples from Pakistan, to characterize the bacterial isolates and their ARGs profiles using phenotypic assay and whole-genome sequencing. Additionally, the plasmid replicons and integron–intergase gene cassettes were investigated in these sequenced isolates. We explored previously overlooked environmental sources and detected the clinically relevant ARGs in these environmental niches, which pose concerns for public and environmental health.
2. Results
2.1. Bacterial Isolation and Species Identification
The thirty-two bacterial isolates were preliminary identified on selective CHROMagar™ mSuperCARBA culture media according to the provided instructions. The reddish pink and uidA positive isolates were identified as E. coli and differentiated from the rest of the 14 isolates. All the thirty-two isolates were whole-genome sequenced on an Illumina Miseq platform with an average coverage of ~70×. The genus and species identities were further confirmed using GAMBIT, fastANI, and BIGSdb identification. The isolates were grouped into 7 different genera and 12 different species. Among them, 18 were E. coli, 4 Klebsiella pneumoniae, and 1 each for Klebsiella quasipneumoniae, Citrobacter telavivensis, Citrobacter portucalensis, Enterobacter cloacae, Raoultella ornithinolytica, Raoultella planticola, Acinetobacter baumannii, Morganella morganii, Serratia ureilytica, and Serratia nevei. The assemblies’ statistics, including specific genome sizes, contig numbers, G + C content, N50 values, ANI (%), and the nearest BIGSdb species identification, are shown in Supplementary File Table S1.
2.2. Phylogroups, MLST, Serotypes, and Capsule Types
Classification of the phylogroups of the 18 sequenced E. coli isolates was performed and they were grouped into 5 phylogroups (Figure 1A). The most prevalent phylogroup was phylogroup group A (50%; n = 9), followed by phylogroup C (28%; n = 5), and phylogroup D (11%; n = 2). One isolate each for phylogroup B1 and F was detected (Figure 1A).
The MLST analysis revealed diverse sequence types for the studied isolates with some belonging to the pathogenic clonal lineages. The 18 E. coli isolates can be grouped into 11 different STs. ST410 was observed as a predominant ST in 5 out of 18 (28%) with all having phylogroup C. ST46 and ST167, both in phylogroup A, were found in 3 and 2 out of 18 isolates (17% and 11%), respectively. The other detected STs belonging to phylogroup A were ST1284, ST361, ST10886, and ST10. The two STs, ST38 and ST405, were the only phylogroup D isolates. The one isolate each for phylogroup B1 and F were ST2351 and ST648, respectively (Figure 1A). The STs reported here, including ST10, ST410, ST167, and ST46, are pathogenic clonal lineages reported globally. ST10, ST167, and ST1284 are classified into Clonal Complex 10 (CC10) according to the EnteroBase, which are the pathogenic lineages [30].
The serotypes were also identified for the E. coli isolates. The most frequent ST, ST410, was found to have different serotypes, including H9 (n = 2), H21 (n = 1), H30 (n = 1) and O1novel:H9 (n = 1). The STs, including ST648 and ST46, have the same serotype O8:H4 (n = 4), while the two ST167 have different serotypes: Onovel31:H5 (n = 1) and Onovel32:H9 (n = 1). The other detected serotypes with STs included ST38 of serotype O153:H30, ST1284 of serotype Onovel32/O9, ST361 of serotype O9:H30, ST405 of serotype O102:H6, ST10886 of serotype H4, ST2351 of serotype O13:H28, and ST10 of serotype O28novel:H21 (Figure 1A).
For Klebsiella, the in silico MLST analysis showed that all four isolates were of different sequence types including ST3830-ILV, ST245-2LV, ST152, ST4291, and ST2534 (Figure 1B). Sequence variations within the capsular polysaccharides, K antigens, and lipopolysaccharides such as O antigens, are another typing method for K. pneumoniae. In the studied isolates, three distinct O capsule types (O1/O2v1, O1/O2v2, and O3b) and four distinct K capsule types (KL62, KL102, KL104 and KL105) were detected (Supplementary File Table S3).
2.3. Antimicrobial Sensitivity Assay (ASA)
The ASA analysis against ten antibacterial agents of thirty-two isolates is presented in Table 1. Overall, the isolates displayed resistance against the cephalosporin, i.e., cefotaxime in 31 out of 32 (97%) isolates, followed by ciprofloxacin in 29 out of 32 (91%) isolates. Regarding carbapenem, 25 out of 32 (78%) isolates were resistant to imipenem. The sensitive drugs identified were gentamicin followed by trimethoprim-sulfamethoxazole. The resistance against these drugs were found in 18 (56%) and 20 (63%) of the 32 isolates as shown in Table 1.
In the E. coli isolates, 11 out of 18 were resistant to all antibiotics tested. Additionally, two isolates were resistant to nine antibiotics, while resistance against eight, seven, and six antibiotics was found in one isolate, three isolates, and one isolate, respectively (Table 1). Resistance to carbapenems (imipenem, meropenem, doripenem, and ertapenem) was noticed in 16 out of 18 (89%) E. coli isolates. Regarding cephalosporins, all isolates (100%) were resistant to cefotaxime and ceftazidime, except for 42AC1, which was sensitive to the latter. Similarly, 42AC1 was the only isolate sensitive to ciprofloxacin and nalidixic acid, while all other E. coli isolates were resistant. Among the 18 E. coli isolates, 13 (72.2%) were resistant to trimethoprim-sulfamethoxazole and gentamicin. The isolate 42AC1 was the least resistant isolate, followed by 76EC1, being sensitive to eight and four antibiotics, respectively.
In the case of Klebsiella, the Klebsiella quasipneumoniae EB7-1 displayed non-susceptibility to all antibiotics, followed by Klebsiella pneumoniae 22EB1, which was only sensitive to trimethoprim-sulfamethoxazole. The isolate EC118 was resistant to meropenem, cefotaxime, and ciprofloxacin. The Klebsiella pneumoniae EC117 was sensitive to ceftazidime and gentamicin, while EC139 was sensitive to imipenem, nalidixic acid, and trimethoprim-sulfamethoxazole in addition to gentamicin (Table 1).
Furthermore, among the nine out of thirty-two isolates with single species, the Raoultella planticola 68EB1 isolate was resistant to all antibiotics tested. Of these nine isolates, seven were resistant to carbapenem, especially imipenem and doripenem. All nine isolates were resistant to cefotaxime, except Raoultella ornithinolytica 17EB1, which was sensitive to cefotaxime but resistant to ciprofloxacin and trimethoprim-sulfamethoxazole. The Citrobacter telavivensis isolate EB1-3 was resistant to meropenem, doripenem, cefotaxime, and ciprofloxacin but sensitive to all other antibiotics. The Citrobacter portucalensis isolate 66EB1 was sensitive to meropenem, ciprofloxacin, and gentamicin while the Enterobacter cloacae EC145 isolate was sensitive to meropenem, nalidixic acid, and trimethoprim-sulfamethoxazole. (Figure 1). Also, the Acinetobacter baumannii 51AC1 isolate was sensitive to nalidixic acid and gentamicin, whereas the Morganella morganii 77AC2 was only sensitive to ertapenem. Both the Serratia isolates, including Serratia ureilytica 78EB1 and Serratia nevei 78EB2, were resistant to imipenem, cefotaxime, and ceftazidime; the former was additionally resistant to imipenem and doripenem (Table 1).
2.4. Detection of ARGs
The whole-genome sequencing of thirty-two ARB recovered from wastewater displayed a diverse resistome, including the two clinically relevant carbapenamases, NDM and OXA. Other important β-lactamases, including CTX-M-15, TEM, SHV, and CMY genes, were also noticed in some of the isolates (Figure 1A).
In E. coli isolates, 15 resistance genes were found in the isolate 23AC, while on average, 8 ARGs were detected in all other isolates. Notably, the critical carbapenamase genes, bla NDM-5 and bla OXA-181, were each present in 6 out of 18 (33%) isolates. The six isolates bearing bla NDM-5 belonged to different STs including ST648, ST1284, ST361, ST405, ST167, and ST10886. On the other hand, bla OXA-181 was found in two STs, including ST46 (n = 2) and ST410 (n = 4). Additional β-lactamases were found, including bla CTX-M-15 in 13 out of 18 (72%) isolates, and bla OXA-1 in 7 out of 18 (38%) isolates. Variants of the bla TEM and bla CMY genes were present in 10 out of 18 (56%) and 3 out of 18 (17%) isolates, respectively, while bla SHV-12 was identified in only one isolate, 68EC2. Regarding quinolone resistance, genes such as qnrS1 and aac(6′)-Ib-cr (which confer resistance to both quinolones and aminoglycosides) were found in 14 out of 18 (77.7%) isolates, and isolate 23AC1 co-harboring both the genes. Sulfonamide and trimethoprim resistance genes (sul and dfrA) were detected in 13 and 11 out of 32 (72% and 61%) isolates, respectively. Isolates 23AC1 and 37EC2 co-harbored two genes each for sulfonamide (sul1 and sul2) and trimethoprim (dfrA12 and dfrA17) on different contigs in their genomes. On the other hand, the isolate 45EC1 contained two dfrA genes (dfrA1 and dfrA12) on separate contigs. Aminoglycosides-resistant genes were detected in 13 out of 18 (72%) isolates, with 10 of these harboring the variants of aadA including the aadA5 and aadA2 genes. Interestingly, the isolate 37EC2 encoded both the aadA5 and aadA2 genes. Other important ARGs, for which phenotypic susceptibility was not tested, included the macrolide resistance gene (mph(A)) in 15 out of 18 (83%) isolates, the chloramphenicol resistance gene (catA1) in 2 out of 18 (11%) isolates, and the tetracycline resistance genes (tet(A) and tet(B)) in 9 out of 18 (50%) isolates. The glycopeptide-resistant gene, ble, was found in all six isolates harboring the bla NDM-5 gene. Additionally, the qacEdelta1 gene associated with quaternary ammonium compound efflux pump transporter, was present in 10 out of 18 (55%) isolates.
Point mutations conferring antibiotic resistance were also detected in the sequenced E. coli isolates. Notably, the pmrB (Y358N) mutation was detected in 6 out of 18 (33%) isolates, known for heteroresistance to colistin. Additionally, 2 chromosomal point mutations conferring resistance to carbapenems were identified in 12 out of 18 (67%) E. coli isolates, including five isolates with the ftsI_N337NYRIN point mutation and seven isolates with the ftsI_I336IKYRI. The glpT_E448K point mutation for fosfomycin resistance was noticed in 15 out of 32 (83%) E. coli isolates. Furthermore, the mutations, blaTEMp_C32T and cyanA_S352T (fosmidomycin resistance) were detected in 3 and 2 out of 18 (17% and 11%) isolates. Mutations in the gyrA, parC, and parE genes responsible for quinolone resistance were also detected in the E. coli isolates. More specifically, the amino acid substitutions, S83L and D837N in gyrA and S80I in the parC gene, were detected in 13 out of 18 (72%) isolates. The S458A mutation in the parE gene was present in 11 of 18 (61%) isolates. Additionally, the S57T and E84G in the parC genes were only detected in one isolate each, as shown in Figure 1A.
Regarding Klebsiella, the Klebsiella quasipneumoniae, EB7-1 isolate harbored 14 ARGs, including the β-lactamases bla OKP-B-15, bla OXA-1, bla CTX-M-15, and bla TEM-1B. Other detected ARGs included the qnrB1, aac(6′)-Ib-cr5, sul2, dfrA14, aph(3″)-Ib, aph(6)-Id, fosA_6, and tet(A). Additionally, isolates were found to harbor genes encoding multi-drug efflux RND (Resistance-Nodulation-Division) transporter periplasmic adaptor subunit (oqxA) and multi-drug efflux RND transporter permease subunit (oqxB). The isolate EC117 included the bla SHV-187, qnrS1, oqxA/B, dfrA14, and fosA_3. The Klebsiella pneumoniae strains EC118 and EC139 harbored the same ARGs, including bla SHV-187, oqxA/B, and fosA_3. Among the Klebsiella isolates, the Klebsiella pneumoniae 22EB1 isolate carried 16 ARGs including β-lactamases such as bla NDM-1, bla OXA-9, bla CTX-M-15, bla DHA-1, bla SHV-145, and bla TEM-1A. This isolate also harbored additional ARGs such as qnrS1, oqxA/B, aadA1, aac(6′)-Ib, aph(3″)-Ib, aph(3′)-VI, aph(6)-Id, ble, and fosA_3 (Figure 1B).
Additionally, amino acid substitutions for tigecycline resistance were detected in the carbapenem-resistant Klebsiella pneumoniae 22EB1 isolate. The mutation in the ramR gene (A19V), which is part of the AcrAB efflux pump, was found in the isolate. However, this has previously been identified to have no significant role in the tigecycline resistance [31].
Moreover, in nine out of the thirty-two isolates from single species, the Citrobacter telavivensis strains EB1-3 harbored bla SED1 and oqxA/B, while Citrobacter portucalensis 66EB1 was found to have bla CMY-119, bla DHA-1, qnrB43, qnrB4, sul1, dfrA7, mph(A), and tet(A). The Raoultella ornithinolytica strains 17EB1 harbored bla LAP-2, bla PLA-1a, qnrS1, sul1, dfrA1, fosA_2, and tet(B), while Raoultella planticola 68EB1 was found to have bla CTX-M-15, bla PLA-5A, bla TEM-1A, qnrB1, oqxB sul2, dfrA14, aac(3)-IId, aph(3″)-Ib, aph(6)-Id, and fosA_2. The Enterobacter cloacae EC145 harbored one gene of bla CMH_3 while the Acinetobacter baumannii 51AC1 harbored bla OXA-70, bla OXA-72, bla ADC_25, and ant(3”)_IIa. In Morganella morganii 77AC2, bla DHA-13, sul1, dfrA10, ant(2”)_Ia, aph(3′)-Ia, catA1, tet(B), and qacEdelta1 were detected. The two Serratia harbored aac(6′)-Ic, and the bla SST-1 and bla SRT-1 in Serratia ureilytica and Serratia nevei, respectively, as shown in Figure 1B. Most of the isolates confirmed the genotypes for a specific phenotype, and a correlation of the genetic resistance mechanism responsible for a phenotype is provided in Table 2.
2.5. Detection of Integrons–Integrases
The isolates from this study were investigated for the presence of integrons, the genetic elements identified for the spread of gene cassettes conferring adaptive functions. Table 3 summarizes the prevalence of integrons and their gene cassettes carrying different genes in the integron-positive isolates. Of the 32 isolates sequenced from this study, 11 harbored a class 1 integron integrase gene. Among these 11, 6 were E. coli, accounting for 55% of the total, while others included Klebsiella quasipneumoniae, Klebsiella pneumoniae, Citrobacter portucalensis, Raoultella planticola, and Morganella morganii. The integron was present on the plasmid-borne contig of all the isolates except for Citrobacter portucalensis, where the integron was found on a chromosomal-borne contig. All the E. coli isolates carried the same array of genes in the gene cassettes: genes encoding resistance to trimethoprim (dfrA), aminoglycosides (aadA), and a semi-functional derivative of a small multi-drug resistance (SMR) protein with a role in biocide resistance (qacEΔ1). The one exception was for 76EC1, which only contained the dfrA gene and a partial IS6 family transposon. The E. coli isolates were found to have three and two variants of the dfrA and the aadA gene, respectively, along with the common qacEΔ1 gene. Among the six E. coli integron-positive isolates, five co-harbor bla NDM-1 while 76EC1 did not (Figure 1A). The integron-positive isolates from the Klebsiella genus, including Klebsiella quasipneumoniae and Klebsiella pneumoniae, were found with the same array of genes: the trimethoprim resistance gene drfA and the bacterial mobilization protein, mobC. The integron with the gene cassette array containing dfrA7 and qacEΔ1was found on the chromosomal contig of a Citrobacter portucalensis. Interestingly, the Raoultella planticola was the only isolate with only the drfA14 gene in an integron cassette on a plasmid. The Morganella morganii isolate 77AC2 was found to carry an aminoglycoside gene, ant(2′)-ia (“2”-aminoglycoside nucleotidyltransferase), along with the qacEΔ1 gene (Table 3).
2.6. Detection of Plasmid Replicons
The PlasmidFinder was used to predict the plasmid replicons in the sequenced isolates. The replicons were grouped into ten types: Col, IncFIA, IncFIB, IncFII, IncI, IncR, IncX, IncY, repB, and p0111. IncFIB was the major replicon detected, found in 17 out of 32 (53%) isolates, followed by Col, IncFIA, and IncFII in 15, 13, and 8 (46%, 40%, and 25%) isolates, respectively.
The IncFIB replicon was also detected in most of the E. coli isolates (12 out of 18 (67%)), followed by IncFIA in 10 out of 18 (56%) isolates. Six replicons were detected in both 23AC1 and 20EC1; in contrast, only one replicon was detected in 23EC1 and 41EC2 (Figure 1A).
In Klebsiella, IncFIB was detected in all isolates except EC139 (no replicon was detected). Klebsiella quasipneumoniae EB7-1 was found to have three replicons (Col, IncFIB, and IncFIA). There were six replicons in 22EB1 and EC117, while the isolate EC118 had only one replicon (Figure 1B).
Moreover, among the 9 out of 32 isolates with single species, there were 3 replicons in Raoultella ornithinolytica and Raoultella planticola isolates. The Citrobacter portucalensis 66EB1 isolate was found to have two Col-type plasmid replicons. No plasmid replicons were predicted by the PlasmidFinder tool in Enterobacter cloacae, Citrobacter telavivensis, Morganella, and the two Serratia isolates (Figure 1B).
2.7. Plasmid Sequence Analysis of E. coli 41EC2 and Localization of Its ARGs
To further investigate the plasmid content, genetic environment, and localization of the bla NDM-5 gene in the sequenced E. coli isolates, a representative 41EC2 isolate was further selected. The combined long and short read sequencing approach was utilized to obtain a full-length complete plasmid sequence from the isolate 41EC2. The combined Illumina and nanopore reads for 41EC2 resulted in a single plasmid, identified as IncFIA/IncFIC/rep_cluster_2244 (p41EC2; GenBank accession: CP155728) with a size of 129,877 bp. Plasmid multilocus sequence typing (pMLST) using the pMLST tool (
3. Discussion
The circulation of Gram-negative ARB in the environment in Pakistan remains poorly understood due to limited data availability. Among ARB, CRB are of concern, having been designated as critical-priority pathogens by the WHO [32]. While most reports of CRB in this region are from clinical settings, their presence in environmental matrices is a growing global concern, with evidence of their spread in community settings, sewage environments, and wastewater sources [33]. The dissemination of Gram-negative ARB isolates into environmental sources poses a serious public health threat. This study documents the presence of clinically relevant ARGs in 32 sequenced bacterial isolates from various sites in Pakistan, including household water, sewage, and canals. These findings suggest the community-level carriage, disposal, and spread of ARB into the environment. The bacterial species were identified using the k-mer-based identification tool, followed by ANI analysis and BIGSdb [34]. The isolates comprised 18 Escherichia coli, 4 Klebsiella pneumoniae and 10 other bacterial species. The phenotypic assay revealed a high resistance rate to cephalosporins, ciprofloxacin, and imipenem with most E. coli isolates showing resistance to carbapenem. These findings are consistent with previous studies on wastewater isolates from Pakistan, which reported high levels of resistance to β-lactams and fluoroquinolones, including cephalosporin, carbapenem, and ciprofloxacin [18,22]. Our results align with those earlier studies, providing further evidence of the escalating burden of antimicrobial resistance in the South Asian region[35,36].
This project also highlights the presence of diverse sequence types (ST) among the bacterial isolates carrying important ARGs, particularly the carbapenemases. Notably, bla OXA-181 and bla NDM-5 were prevalent among the E. coli isolates in this study. The bla OXA-181 gene was found in two STs, the ST361 and ST46, which to our knowledge, are reported for the first time in wastewater isolates from Pakistan. Previous studies have reported the bla NDM-5 gene in E. coli from human and black kites (Milvus migrans) in the same region, with a wide range of STs, including ST167, ST361, ST648, and ST1284 [37,38,39]. Our study extended these findings by identifying these high-risk, human-associated E. coli clones in wastewater isolates. Importantly, this study reports the detection of a rare E. coli ST10886 isolate, 45EC1, belonging to phylogroup A with serotype H9. This isolate harbor the bla NDM-5 and bla CMY-16 genes, along with a pmrB (Y358N) mutation conferring resistance to colistin (Figure 1). According to the Escherichia/Shigella Enterobase database (
Additionally, Extended Spectrum β-lactamases (ESBLs), including the Cefotaximase-Munich (CTX-M), Sulfhydryl Variable (SHV), and Temoneria (TEM) families, were identified in the studied isolates. Our results showed that CTX-M-15 was the most prevalent ESBL, present in 72% (13 out of 18) of the E. coli isolates. Notably, 10 of these 13 isolates co-harbored carbapenemases, specifically the bla NDM-5 and bla OXA-181. The high prevalence of ESBLs, particularly CTM-15, has been consistently reported in wastewater samples from this region [17].
Additionally, resistance to β-lactam was also observed due to the presence of AmpC types β-lactamases. AmpC β-lactamases are important cephalosporinases encoded either chromosomally (intrinsic) or by plasmids (acquired) [45]. In this study, the plasmid mediated AmpC β-lactamase bla CMY, which was previously reported as the first plasmid mediated AmpC β-lactamase [46], was only detected in three E. coli isolates and one Citrobacter portucalensis isolate (66EB1). Other types of AmpC β-lactamases identified in the sequenced isolates included the bla OKP-B-15 gene, a chromosomal β-lactamase gene found in Klebsiella quasipneumoniae, the bla CMH-3 gene in Enterobactor cloaceae, ADC (Acinetobacter-derived cephalosporinases) in Acinetobacter baumannii, and bla DHA-1 in Klebsiella quasipneumoniae, Citrobacter portucalensis 66EB1, and Morganella morganii. Furthermore, the bla PLA-1a and bla PLA-5a genes, identified as class A β-lactamase, were present in Raoultella species. These genes have previously been reported as part of the Raoultella chromosome, providing the bacteria with intrinsic resistance to β-lactamases [47]. Another class A β-lactamase, bla LAP-2, known for its narrow resistance spectrum against β-lactam antibiotics, was identified alongside the qnrS1 gene in Raoultella ornithinolytica 17EB1 on a single contig. Previously, a similar gene arrangement was found in Col- and IncFIB-type plasmids in Klebsiella pneumoniae [48]. The presence of Col- and IncFIB-type replicons in these isolates suggests a similar backbone structure to that found in Klebsiella pneumonia, now reported in Raoultella ornithinolytica from Pakistan.
The phenotypic resistance to quinolones was observed in most isolates, with the exception of 42AC1 and a few others. This resistance was primarily due to acquired genes such as qnrS1, aac(6)-Ib-cr, qnrB1, qnrB43 and qnrB4, along with multi-drug efflux proteins oqxA and oqxB, particularly in Klebsiella and one Raoultella species. Additionally, mutations in quinolone resistance-determining regions (QRDRs) were detected in over 72% of E. coli isolates, specifically the amino acid substitutions D87N, S83L, and S458A in DNA gyrase (gyrA) and topoisomerase IV (parE). These substitutions have been strongly associated with high levels of quinolone resistance in E. coli, consistent with previous reports [49]. Interestingly, although the E. coli isolate 41AC1 carried genetic markers for quinolone resistance, it may be non-susceptible to other quinolones, such as levofloxacin, norfloxacin, and ofloxacin, which were not included in the phenotypic sensitivity tests conducted in this study.
Aminoglycoside resistance genes were in 72% of the E. coli isolates, with the aadA5 gene being most prevalent. Notably, some E. coli isolates (42AC1, 19EC2, and 76EC1) possessed one or more aminoglycoside-modifying enzymes (AMEs) that did not correspond with the phenotypic resistance patterns observed (Table 2). This discrepancy suggests potential resistance to additional aminoglycosides, such as amikacin, streptomycin, kanamycin, and tobramycin, which were not tested in the phenotypic assays conducted in this study.
Additionally, genes conferring resistance to macrolide, fosfomycin and tetracyclines were also identified in the studied isolates. These findings are consistent with previous reports from the same region, which detected similar resistance genes in different bacterial isolates [50]. Interestingly, in our study, the fosA gene was exclusively found in the Klebsiella and Raoultella species. It may be attributed to the presence of the fosA gene in the chromosome, which confers intrinsic resistance to Fosfomycin and occurrence that is rarely identified in E. coli [51].
Previously, colistin resistance due to the mcr-1 gene has been reported in Pakistan in E. coli isolates originating from humans, wild migratory birds, healthy broiler chickens [52], and sewage samples [20]. The later study highlighted a 35% resistance rate to colistin in E. coli, underscoring the significant burden of colistin resistance in environmental settings compared to clinical contexts, which is a serious concern. In our study, we observed a mutation in the pmrB gene (Y358N) in 33% of E. coli isolates. This nucleotide substitution is noteworthy as it has been previously associated with colistin non-susceptibility, exhibiting a high MIC > 32 μg/mL even in the absence of the mcr gene [53]. This suggests that the single amino acid change could contribute to colistin resistance. Given these isolates were derived from wastewater, the presence of colistin in such environments may have driven the emergence of resistance. The potential for these resistant isolates to spread to other sources, such as community water supplies, poses a significant public health risk. Therefore, the detection of this mutation in wastewater-borne E. coli is concerning and underscores the need for ongoing surveillance and proactive measures to prevent the dissemination of resistant strains.
The identified plasmid replicons indicate that most of the E. coli and Klebsiella isolates were prone to plasmid uptake, which may facilitate horizontal gene transfer and contribute to the dissemination of critical ARGs among pathogenic bacteria. The bla NDM-5 isolates exhibited a combination of diverse plasmid replicons, including Col, IncFIA, IncFIB and IncFII types. These replicon types have previously been associated with bla NDM-5 genes from clinical samples [37]. Additionally, our study also identified IncX3 replicons in two isolates carrying the bla OXA-181 gene. These findings are consistent with the previous observation linking InX3 replicons to the bla OXA-181 gene in Enterobacterales [54].
The bla NDM-1 gene in Klebsiella pneumoniae 22EB1 ST2534 is likely carried on an IncFIB(pQil) replicon-type plasmid. Similarly, the ESBL gene in this same isolate is probably located on a repB(R1701) replicon-type plasmid. These replicon-type plasmids have been recently identified in clinical isolates associated with the ESBL antibiotic group in Pakistan [43]. In contrast, there were no plasmid replicons detected in the Acinetobacter baumannii 51AC1 isolate, which harbors the carbapenamase bla OXA-72. This could be because the novel replicon type is not present in the database, was missed due to short read sequencing, and/or was present on the chromosome, all of which merit further investigation.
The representative bla NDM-5-positive E. coli isolate 41EC2 was additionally sequenced by using both Illumina and Oxford nanopore platforms to assemble the full-length plasmid sequence. Combined Illumina and nanopore sequencing reads revealed a single conjugative plasmid in the E. coli isolate 41EC2. BLASTn analysis identified this plasmid as being closely related to plasmids in the NCBI database, with 100% coverage and nucleotide identity to other ST167 E. coli isolates from USA, including the pMB9245_1 plasmid, GenBank accession No. CP103530.1, with a size of 133,836 bp. Other closely related plasmids included pMyNCGM167_1 (GenBank accession No. LC744455.1, 132,545 bp from Myanmar), pM309-NDM5 (GenBank accession No. AP018833.1, 136,947 bp from Myanmar), and a plasmid from the E. coli strain NDM_11.16372 (GenBank accession No. CP082130.1, 136,961 bp from Tanzania), with 100% coverage and 99.99% nucleotide identities, respectively. Concentric ring-based circular comparisons between p41EC2 and similar plasmids (Figure 2) demonstrated that all these plasmids share the IncFIA/IncFIC/InFII-type replicon, consistent with the IncFIA/IncFIC/rep_cluster_2244 type of p41EC2. Previous reports have described similar IncF-type plasmids where the bla NDM-5 is flanked by two copies of IS26, and with ISCR1 located in between [37,55]. These results highlight the critical roles of IS26 and ISCR1 in the transmission of bla NDM-5 and underscore their significance in the dissemination of antibiotic-resistant elements. The presence of bla NDM-5 in this conjugative plasmid likely facilities its spread. The identification of similar plasmids in isolates from various countries suggests a global distribution of such plasmid types.
Furthermore, the intergon–integrase elements were identified in the sequenced isolates. The class 1 integron was found in 55% of the isolates from this study, which is reported for the first time from the wastewater-borne sequenced isolates from Pakistan. Previous studies reported similar percentages of class 1 integrons in clinical isolates [56]. The integrons found in the E. coli isolates shared the same genetic array as those present in a clinical class 1 integron [12]. The presence of qacEΔ1 at the 3′ conserved region suggests it is a clinical class 1 integron that has made its way into the wastewater. Interestingly, this study identified the presence of the mobC gene in the integron cassette. The mobC is a plasmid-encoded protein making up the relaxosome and assists in the conjugative transfer of plasmids [57]. The integrons present in plasmid-borne contig indicate the possibility of disseminating these genes cassettes into different pathogenic and non-pathogenic bacterial isolates in the wastewater system. The co-occurrence of integrons and the bla NDM-5 gene in E. coli isolates underscores the ability of the genomes to adapt and evolve in response to the environmental changes and conditions. The integrons present in these isolates offer a probable route of transmitting these gene cassettes or transmitting them into a new array with other important resistance genes. This may pose a serious threat to drinking water facilities, such as those in Pakistan with limited wastewater treatment options.
This study highlights the detection of clinically relevant ARGs in bacterial isolates from wastewater, emphasizing the need for a One Health-oriented molecular surveillance approach to combat the growing threats of antibiotic resistance. The numerous reservoirs and dissemination routes of AGRs suggest the importance of integrated efforts to address this issue.
4. Materials and Methods
4.1. Sample Collection, Microbiological Analysis, and Antimicrobial Sensitivity Assay (ASA)
Wastewater samples from different sites of Pakistan, including sewage, household wastewater, and canals of different sites, were collected over a two-year period from August 2020 to July 2022. Fifty ml of each water sample was collected in sterile falcon tubes, transported, and processed immediately at the Bacterial Pathogen Laboratory, Health Biotechnology Division, National Institute for Biotechnology and Genetic Engineering College, Pakistan Institute of Engineering and Applied Sciences (NIBGE-C, PIEAS) Faisalabad, Pakistan. Otherwise they were kept at 4 °C to be processed on the next day for isolation. The collected samples were enriched in buffered peptone water (BPW; Oxoid, UK) in a ratio of 1:9 (sample:BPW) and incubated at 37 °C for 18–24 h. One loop-full of the enriched culture was streaked onto CHROMagar™ mSuperCARBA™ plates (CHROMagar, Saint-Denis, France), which is a selective and differential chromogenic culture medium for the isolation and detection of CRB, and then incubated at 37 °C for 18–24 h. Based on the species-enrichment guidelines provided by CHROMagar™ mSuperCARBA™ media, the isolates were grouped into three categories; pink for suspected E. coli, metallic blue colonies for presumptive Klebsiella, Enterobacter, or Citrobacter, and cream, opaque colonies for suspected Acinetobacter. These colored colonies were further streaked on MacConkey agar (Oxoid: CM0007) and CHROMagar™ mSuperCARBA™ plates to purify the individual colonies. The pure presumptive E. coli isolates were further identified for the presence of the uidA gene that encodes a β-D-glucuronidase enzyme used for the identification of E. coli. PCR profile and conditions were followed as published previously [58].
Antibiograms of each of the isolates were detected using standard Kirby–Bauer disc diffusion protocol [59] and the results were interpreted using Clinical and Laboratory Standards Institute (CLSI) [60] guidelines and breakpoints. Briefly, the bacterial cultures were grown overnight in Tryptone Soya Broth (TSB; Oxoid: CM0129, UK), diluted to the equal turbidity of freshly prepared 0.5 MacFarland standard and spread as a lawn culture on Mueller–Hinton agar (MHA; Oxoid: CM0337). The E. coli ATCC 25922 strain was used as quality control for this assay. The antibiotic discs were dispensed on MHA using a 6-cartridge antimicrobial susceptibility disc dispenser (Oxoid, UK). The isolates were tested against a panel of 10 antimicrobial discs (Oxoid, UK), commonly used in human and veterinary medicine including imipenem (10 μg), meropenem (10 μg), ertapenem (10 μg), doripenem (10 μg), ceftazidime (30 μg), cefotaxime (5 μg), ciprofloxacin (5 μg), gentamicin (10 μg), trimethoprim/sulfamethoxazole (25 μg), and nalidixic acid (30 µg). The diameters of inhibition zones were determined, and the isolates were categorized as resistant, intermediate resistant, or susceptible using the CLSI guidelines. The isolates with intermediate resistant were considered as “resistant″ in the Section 2.
4.2. Whole-Genome Sequencing and Bioinformatics Analysis
Thirty-two isolates were selected for whole-genome sequencing based on the following criteria: (a) phenotypic resistance to carbapenem and cephalosporins; (b) samples collected from different provinces/locations of Pakistan; (c) susceptible isolates for the comparative study; and (d) based on the ARGs’ specific PCR results (unpublished data). A single colony of selected isolates was suspended in 1 mL of TSB at 37 °C overnight, and the genomic DNA was extracted using the GeneJET DNA purification kit (Thermo Fisher Scientific, Waltham, Massachusetts, USA). The DNA samples were visualized and confirmed on ethidium bromide-stained 2% agarose gel electrophoresis system. The genomic libraries were prepared and sequenced using a previously published protocol [61]. Briefly, the Nextera XT DNA library preparation kit (Illumina, San Diego, CA, USA) was used to generate paired-end libraries, followed by sequencing on an Illumina MiSeq (Illumina, San Diego, CA, USA) using the MiSeq Reagent V2 500 cycle kit (Illumina, USA). The representative E. coli isolate harboring the bla NDM-5 gene, 41EC2, was additionally sequenced using long-read sequencing approach on a Mk1C device (Oxford Nanopore Technologies, Oxford, UK) using Native Barcoding Kit 24 V14 (SQK-NBD114.24) with R10.4.1 flow cell. The bases were named using Guppy v6.4.6, and the assembly of Illumina and nanopore sequencing reads were performed using Unicycler v0.4.9b to ascertain the whole-genome and plasmid sequences for this isolate [62].
The raw reads were processed using the same protocol as published previously [34]. Briefly, the raw reads obtained were checked for quality using FastQC v0.12.1 and multiQC v1.19. The adapters and low-quality reads were removed to achieve a Phred score of 33 with Trimmomatic v0.39. The quality of trimmed reads was checked thorough FastQc and multiQC. SPAdes v3.15.5 was performed on the trimmed reads for de novo assembly of the genomes. The contigs with length and coverage less than 500 bp and 2X, respectively, were removed from the assembled genomes [63]. The quality of the assemblies was checked with QUAST v 3.9, while the contamination and completeness of the assembled genomes were evaluated with Checkm2 [64,65]. The bacterial isolates were identified using the following approaches: (a) the k-mer approach, using the Genomic Approximation Method for Bacterial Identification and Tracking tool (GAMBIT) [66]; (b), the average nucleotide identity analysis using FastANI v1.33, including all the species reference genome sequences available at NCBI RefSeq database for a genus (list of genomes with accessions available in Supplementary File Table S2) [67]; and (c) cross-checking of species confirmations with the Public databases for molecular typing and microbial genome diversity (BIGSdb) (
The assemblies were searched for integrons using IntegronFinder (v2.0) with options [--local-max --func-annot –gbk] [72,73]. The gene cassettes from IntegronFinder’s GenBank-formatted output file were further validated and identified with Prokka, Resfinder and BLASTp. Additionally, sequence searches were performed using DIAMOND in eggNOG-mapper (v 2.1.12) to functionally annotate the gene cassettes based on the eggNOG orthology data [74]. Additionally, Plasmid Hunter was used to screen for integron-positive contigs and to classify them as plasmid- or chromosomal-borne contigs [75]. The plasmid replicons were predicted using PlasmidFinder in ABRicate (
5. Limitations of the Study
This study has some limitations. The bacterial isolates analyzed may not fully represent the overall resistance in wastewater, as only thirty-two isolates were selected from different locations. This limited sampling may be an underestimation or overestimation of CRB presence. Additionally, the large amount of debris and fungal growth also presents a challenge in obtaining single bacterial colonies, ultimately leading to the underestimation of ARB in the samples. Future studies are needed to assess the resistome content of bacterial isolates nationwide and by incorporating more samples collected at different time points of the year to provide a more comprehensive understanding of the resistome in wastewater.
6. Conclusions
For the first time, the data presented here highlight the role of environment samples as hotspots for WHO-priority bacterial pathogens and ARGs in Pakistan. The detection of clinically relevant ARGs, particularly the carbapenemases variants, bla NDM and bla OXA, detected in this study, underscores the role of wastewater as a significant reservoir of carbapenmases. The presence of carbapenem- and colistin-resistant isolates, the last-resort antibiotics, suggests a large arsenal of mechanisms in these isolates under antimicrobial pressure and is a serious threat to public health and environmental integrity. Therefore, it is crucial to focus on community wastewater to prevent the dissemination of these key ARGs by novel genetic arrangements via mobile genetic elements, especially the integrons, or by horizontal gene transfer of bla NDM-5 on conjugative plasmids. Implementing effective preventive and control measures is essential to reducing the spread of these emerging threats. This research contributes valuable insights into the dynamics of antimicrobial resistance from a One Health perspective.
Conceptualization, M.A.G. and A.A.; data curation, F.S. and G.L.; formal analysis, F.S., X.H. and A.S.; investigation, F.S., K.M. and H.S.; methodology, H.A., M.A.G., Y.S., A.A. and G.L.; project administration, A.A.; resources, H.A., Y.S., A.A. and G.L.; software, X.H. and A.S.; supervision, A.A. and G.L.; validation, F.S., X.H. and G.L.; writing—original draft, F.S.; writing—review and editing, X.H., K.M., H.S., A.A. and G.L. All authors have read and agreed to the published version of the manuscript.
Not applicable.
Not applicable.
The raw reads for all the studied isolates are available under Bio Project PRJNA1106557 and the list of NCBI accessions for these isolates are included in
We thank Iowa State University Research IT for use of computational resources to execute the bioinformatics analysis of this project.
The authors declare no conflicts of interest.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Enlarge this image.Figure 1. Summary of ARGs detected in the sequenced isolates, (A) E. coli and (B) rest of the isolates. The species column indicates the predicted species for each isolate. The phylotypes and serotypes for E. coli were determined along with the MLST for each isolate if the scheme was available. * Pasteur MLST scheme for Acinetobacter baumannii 51AC1. “-” in the MLST indicates that the MLST scheme is not available or an unknown MLST. The total ARGs are summarized by bacterial isolate identified in silico, and conferring resistance as summarized by antibiotic category. (A) Summary of ARGs, point mutations and plasmid replicons in E. coli isolates from this study. (B) The ARGs in the sequenced isolates are those belonging to Klebsiella and other individual species sequenced in this study. BTL, β-lactams; QNL, quinolone; SUL, sulfonamide; TMP, trimethoprim; AMG, aminoglycoside; C, chloramphenicol; FOS, Fosfomycin; M, macrolide; T, tetracycline; G, Glycopeptide; EFP; multi-drug efflux RND transporters; Q; quaternary ammonium compound; T; point mutations for tigecycline. FD, fosmidomycin; Plasmids; the plasmids replicons hits for each isolate.
Enlarge this image.Figure 1. Summary of ARGs detected in the sequenced isolates, (A) E. coli and (B) rest of the isolates. The species column indicates the predicted species for each isolate. The phylotypes and serotypes for E. coli were determined along with the MLST for each isolate if the scheme was available. * Pasteur MLST scheme for Acinetobacter baumannii 51AC1. “-” in the MLST indicates that the MLST scheme is not available or an unknown MLST. The total ARGs are summarized by bacterial isolate identified in silico, and conferring resistance as summarized by antibiotic category. (A) Summary of ARGs, point mutations and plasmid replicons in E. coli isolates from this study. (B) The ARGs in the sequenced isolates are those belonging to Klebsiella and other individual species sequenced in this study. BTL, β-lactams; QNL, quinolone; SUL, sulfonamide; TMP, trimethoprim; AMG, aminoglycoside; C, chloramphenicol; FOS, Fosfomycin; M, macrolide; T, tetracycline; G, Glycopeptide; EFP; multi-drug efflux RND transporters; Q; quaternary ammonium compound; T; point mutations for tigecycline. FD, fosmidomycin; Plasmids; the plasmids replicons hits for each isolate.
Phenotypic antimicrobial susceptibility profile of the sequenced isolates from this study.
| Sample ID | Species | Antimicrobial Susceptibility Profile | |||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| β-lactams | Quinolones | Combination | Aminoglycosides | ||||||||
| Carbapenem | Cephalosporins | ||||||||||
| IPM | MEM | DOR | ETP | CTX | CAZ | CIP | NAL | SXT | GEN | ||
| 19EC2 | Escherichia coli | I | S | S | I | R | R | R | R | R | S |
| 20EC1 | Escherichia coli | R | R | R | R | R | R | R | R | R | R |
| 21EC2 | Escherichia coli | R | R | R | R | R | R | R | R | R | R |
| 22EC1 | Escherichia coli | R | R | R | R | R | R | R | R | I | I |
| 22EC2 | Escherichia coli | R | R | R | R | R | R | R | R | S | R |
| 23EC1 | Escherichia coli | R | R | R | R | R | R | R | R | S | I |
| 23AC1 | Escherichia coli | R | R | R | R | R | R | R | R | R | R |
| 26AC1 | Escherichia coli | R | R | R | R | R | R | R | R | R | I |
| 30EC1 | Escherichia coli | R | R | R | R | R | R | R | R | R | R |
| 36AC2 | Escherichia coli | R | R | R | R | R | R | R | R | R | I |
| 37EC1 | Escherichia coli | S | R | I | R | R | R | R | I | S | S |
| 37EC2 | Escherichia coli | R | R | R | R | R | R | R | R | R | R |
| 41EC2 | Escherichia coli | R | R | R | R | R | R | R | R | R | I |
| 42AC1 | Escherichia coli | S | S | S | R | R | S | S | S | S | S |
| 43EC1 | Escherichia coli | R | I | R | R | R | R | R | R | R | R |
| 45EC1 | Escherichia coli | R | R | R | R | R | R | R | R | R | I |
| 68EC2 | Escherichia coli | I | I | R | S | R | R | R | R | S | S |
| 76EC1 | Escherichia coli | I | S | S | S | R | I | I | I | R | S |
| EB7-1 | Klebsiella quasipneumoniae | R | R | R | R | R | R | R | R | R | I |
| EC117 | Klebsiella pneumoniae | I | I | I | I | R | S | R | I | I | S |
| EC118 | Klebsiella pneumoniae | S | I | S | S | I | S | I | S | S | S |
| EC139 | Klebsiella pneumoniae | S | I | I | I | R | I | R | S | S | S |
| 22EB1 | Klebsiella pneumoniae | R | R | R | R | R | R | R | R | S | I |
| EB1-3 | Citobacter telavivensis | S | R | I | S | R | S | I | S | S | S |
| 66EB1 | Citobacter portucalensis | R | S | I | I | I | R | S | R | R | S |
| 17EB1 | Raoultella ornithinolytica | S | S | S | S | S | S | R | S | R | S |
| 68EB1 | Raoultella planticola | I | I | I | R | R | R | I | I | R | R |
| EC145 | Enterobacter cloacae | R | S | I | I | R | I | R | S | S | R |
| 51AC1 | Acinetobacter baumannii | R | R | R | R | R | I | R | S | I | S |
| 77AC2 | Morganella morganii | R | I | R | S | R | R | R | R | R | R |
| 78EB1 | Serratia ureilytica | R | I | R | S | R | I | S | S | S | S |
| 78EB2 | Serratia nevei | I | S | S | S | R | I | I | S | S | S |
Definitions: IPM: imipenem; MEM: meropenem DOR: doripenem: ETP: ertapenem; CTX: cefotaxime; CAZ: ceftazidime; CIP: ciprofloxacin; NAL: nalidixic Acid; SXT: trimethoprim-sulfamethoxazole; GEN: gentamycin; abbreviation of standard antibiotic codes from (
Correlation between phenotypic resistance profiles and the genetic mechanisms detected in the isolates sequenced in this study.
| Sample ID | Species | βlactam | Quinolones | Combination | Aminoglycosides | No Phenotypes Tested | ||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Carbapenem | Cephalosporins | |||||||||||
| Phenotype | Genetic Mechanism | Phenotype | Genetic Mechanism | Phenotype | Genetic Mechanism | Phenotype | Genetic Mechanism | Phenotype | Genetic Mechanism | Other Genetic Mechanism | ||
| 19EC2 | Escherichia coli | I (IPM, ETP), S (MEM, DOR) | fstI I336IKYR1 | R | OXA-1, CTX-M-15, CMY-16 | R | aac(6′)-Ib-cr_1, gyrA (S83L and D837N), parC (S80I and S57T), parE (S458A) | R | sul2, dfrA1 | S | aph(3″)-Ib_2, aph(6)-Id_1 | tet(A), |
| 20EC1 | Escherichia coli | R | NDM-5, fstI N337NYRIN | R | OXA-1, CTX-M-15, TEM-1B | R | aac(6′)-Ib-cr_1, gyrA (S83L and D837N), parC (S80I), parE (S458A) | R | sul1, dfrA17 | R | aadA5, aac(3)-IIa | catA1, mph(A), tet(B), ble, qacEdelta1, glpT_E448K, cyanA_S352T |
| 21EC2 | Escherichia coli | R | NDM-5, fstI I336IKYR1 | R | OXA-1, CTX-M-15 | R | aac(6′)-Ib-cr_1, gyrA (S83L and D837N), parC (S80I), parE (S458A) | R | sul1, dfrA17 | R | aac(6′)-Ib-cr_1, aadA5, aac(3)-IIa | mph(A), ble, qacEdelta1, glpT_E448K |
| 22EC1 | Escherichia coli | R | OXA-181 | R | CTX-M-15, TEM-35, blaTEMp_C32T | R | qnrS1 | I | I | glpT_E448K | ||
| 22EC2 | Escherichia coli | R | Nil | R | CTX-M-15, TEM-35, blaTEMp_C32T | R | qnrS1 | S | R | glpT_E448K | ||
| 23EC1 | Escherichia coli | R | Nil | R | TEM-1B | R | gyrA (S83L and D837N), parC (S80I), parE (S458A) | S | Nil | I | Nil | mph(A), pmrB_Y358N, glpT_E448K |
| 23AC1 | Escherichia coli | R | OXA-181, fstI N337NYRIN | R | OXA-1, CTX-M-15, CMY-2 | R | qnrS1, aac(6′)-Ib-cr_1, gyrA (S83L and D837N), parC (S80I), parE (S458A), parE (S458A) | R | sul1, sul2 | R | aac(6′)-Ib-cr_1, aadA5, aac(3)-IId, | mph(A), tet(B), qacEdelta1, pmrB_Y358N, glpT_E448K |
| 26AC1 | Escherichia coli | R | OXA-181, fstI I336IKYR1 | R | CTX-M-15 | R | qnrS1, gyrA (S83L and D837N), parC (S80I), parE (S458A) | R | sul1, dfrA17 | I | aadA5, | mph(A), qacEdelta1, pmrB_Y358N, glpT_E448K |
| 30EC1 | Escherichia coli | R | NDM-5, fstI N337NYRIN | R | OXA-1 | R | gyrA (S83L and D837N), parC (S80I), parE (S458A) | R | sul1, dfrA12 | R | aadA2, | glpT_E448K |
| 36AC2 | Escherichia coli | R | OXA-181, fstI I336IKYR1 | R | CTX-M-15, TEM-169 | R | qnrS1, gyrA (S83L and D837N), parC (S80I), parE (S458A) | R | sul1, dfrA17 | I | aadA5, | mph(A), qacEdelta1, pmrB_Y358N, glpT_E448K |
| 37EC1 | Escherichia coli | R (MEM, ETP), I (DOR), S (IMP) | OXA-181 | R | CTX-M-15, TEM-35, blaTEMp_C32T | R (CIP), I (NAL) | qnrS1 | S | Nil | S | Nil | mph(A), glpT_E448K |
| 37EC2 | Escherichia coli | R | NDM-5, fstI I336IKYR1 | R | OXA-1, CTX-M-15 | R | aac(6′)-Ib-cr_1, gyrA (S83L and D837N), parC (S80I), parE (S458A) | R | sul1, dfrA12, dfrA17 | R | aadA2, aadA5, aac(3)-IIa | mph(A), tet(B, ble, qacEdelta1, glpT_E448K, cyanA_S352T |
| 41EC2 | Escherichia coli | R | NDM-5, fstI N337NYRIN | R | OXA-1, CTX-M-15, TEM-1B | R | aac(6′)-Ib-cr_1, gyrA (S83L and D837N), parC (S80I), parE (S458A) | R | sul1, dfrA12 | I | aadA2, aac(6′)-Ib-cr_1 | mph(A), tet(A), ble, qacEdelta1, |
| 42AC1 | Escherichia coli | R (ETP), S (IPM, MEM, DOR) | OXA-181, fstI I336IKYR1 | R (CTX), S (CAZ) | CTX-M-15, TEM-169 | S | qnrS1, gyrA (S83L and D837N), parC (S80I), parE (S458A) | S | sul1 | S | aadA5 | mph(A), qacEdelta1, pmrB_Y358N, glpT_E448K |
| 43EC1 | Escherichia coli | R (IPM, DOR, ETP), I (MEM) | fstI N337NYRIN | R | TEM-35 | R | gyrA (S83L and D837N), parC (S80I and E84G) | R | sul2, dfrA1 | R | aph(3″)-Ib_5 | mph(A), tet(B), |
| 45EC1 | Escherichia coli | R | NDM-5, fstI I336IKYR1 | R | CMY-16 | R | gyrA (S83L and D837N), parC (S80I) | R | sul1, dfrA1, dfrA12 | I | aadA2, | mph(A), tet(A), ble, qacEdelta1, pmrB_Y358N, glpT_E448K |
| 68EC2 | Escherichia coli | R (DOR), I (IPM, MEM), S (ETP) | Nil | R | SHV-12 | R | qnrS1 | S | Nil | S | Nil | glpT_E448K |
| 76EC1 | Escherichia coli | I (IPM), S (ETP, MEM, DOR) | Nil | R (CTX), I (CAZ) | CTX-M-15, TEM-1B | I | qnrS1 | R | sul2, dfrA14 | S | aph(3″)-Ib_5, aph(6)-Id | tet(A) |
| EB7-1 | Klebsiella quasipneumoniae | R | Nil | R | OKP-B-15, OXA-1, CTX-M-15, TEM-1B | R | qnrB1, aac(6′)-Ib-cr5, oqxA/B | R | sul2, dfrA14 | I | aph(3″)-Ib_5, aph(6)-Id | fosA_6, tet(A) |
| EC117 | Klebsiella pneumoniae | I | Nil | R (CTX), | SHV-187 | R (CIP), I (NAL) | qnrS1, oqxA/B | I | dfrA14 | S | Nil | fosA_3 |
| EC118 | Klebsiella pneumoniae | I (MEM), S (IPM, ETP, DOR) | Nil | I (CTX), S (CAZ) | SHV-187 | I (CIP), S (NAL) | oqxA/B | S | Nil | S | Nil | fosA_3 |
| EC139 | Klebsiella pneumoniae | I (MEM, ETP, DOR), S (IPM) | Nil | R (CTX), I (CAZ) | SHV-187 | R (CIP), S (NAL) | oqxA/B | S | Nil | S | Nil | fosA_3 |
| 22EB1 | Klebsiella pneumoniae | R | NDM-1 | I | OXA-9, CTX-M-15, DHA-1, SHV-145, TEM-1A | R | qnrS1, oqxA/B | S | Nil | I | aadA1, aac(6′)-Ib, aph(3″)-Ib_2, aph(3′)-VI, aph(6)-Id | ble, fosA6 |
| EB1-3 | Citobacter telavivensis | R(MEM), I(DOR), S (IPM, ETP) | Nil | R (CTX), | SED1 | I (CIP), S (NAL) | oqxA/B | S | Nil | S | Nil | Nil |
| 66EB1 | Citobacter portucalensis | R (IPM), I (ETP, DOR), S (MEM) | Nil | R (CAZ), | CMY-119, DHA-1 | S (CIP), R (NAL) | qnrB43, qnrB4 | R | Sul1, dfrA7 | S | Nil | mph(A), tet(A) |
| 17EB1 | Raoultella ornithinolytica | S | Nil | S | LAP-2, PLA-1a | R (CIP), S (NAL) | qnrS1 | R | Sul1, dfrA1 | S | Nil | fosA_2, tet(B), qacEdelta1 |
| 68EB1 | Raoultella planticola | R (ETP), I (IMP, MEM, DOR) | Nil | R | CTX-M-15, PLA-5A, TEM-1A | I | qnrB1, oqxB | R | Sul2, dfrA14 | R | aac(3)-IId, aph(3″)-Ib_5, aph(6)-Id | fosA_2 |
| EC145 | Enterobacter cloacae | R (IPM), I (ETP, DOR), S (MEM) | Nil | R (CTX), | CMH_3 | R (CIP), S (NAL) | Nil | S | Nil | R | Nil | Nil |
| 51AC1 | Acinetobacter baumannii | R | OXA-72 | R (CTX), | OXA-70, ADC_25 | R (CIP), S (NAL) | Nil | I | Nil | S | ant(3″)_IIa | Nil |
| 77AC2 | Morganella morganii | R (IMP, DOR), I (MEM), S (ETP) | Nil | R | DHA-13 | R | Nil | R | sul1, dfrA10 | R | ant(2″)_Ia, aph(3′)-Ia | catA1, tet(B), qacEdelta1 |
| 78EB1 | Serratia ureilytica | R (IMP, DOR), I (MEM), S (ETP) | Nil | R (CTX), | SST-1 | S | Nil | S | Nil | S | aac(6′)-Ic | |
| 78EB2 | Serratia nevei | I (IMP), S (MEM, ETP, DOR) | Nil | R (CTX), | SRT-1 | I (CIP), S (NAL) | Nil | S | Nil | S | aac(6′)-Ic | Nil |
Detection of integron elements and their associated gene cassettes in the wastewater isolates using IntegronFinder 2.0.
| Sample ID | Species | Contig Origin | Size of Integron (bp) | Gene Cassette(s) and Order |
|---|---|---|---|---|
| 20EC1 | Escherichia coli | P | 2802 | dfrA17-aadA5-qacEΔ1 |
| 21EC2 | Escherichia coli | P | 2802 | dfrA17- aadA5-qacEΔ1 |
| 30EC1 | Escherichia coli | P | 3366 | dfrA12-aadA2-qacEΔ1 |
| 41EC2 | Escherichia coli | P | 3240 | dfrA12-aadA2-qacEΔ1 |
| 45EC1 | Escherichia coli | P | 3240 | dfrA12-aadA2-qacEΔ1 |
| 76EC1 | Escherichia coli | P | 1952 | dfrA14-Tnp(partial) |
| EB7-1 | Klebsiella quasipneumoniae | P | 2130 | dfrA14-mobC |
| EC117 | Klebsiella pneumoniae | P | 2130 | dfrA14-mobC |
| 66EB1 | Citobacter portucalensis | C | 1691 | dfrA7-qacEΔ1 |
| 68EB1 | Raoultella planticola | P | 1723 | dfrA14 |
| 77AC2 | Morganella morganii | P | 2200 | ant(2″)-Ia- qacEΔ1 |
P = Plasmid, C = Chromosome.
Supplementary Materials
The following supporting information can be downloaded at:
References
1. WHO. List of Drug-Resistant Bacteria Most Threatening to Human Health Available; WHO: Geneva, Switzerland, 2024.
2. Jean, S.S.; Harnod, D.; Hsueh, P.R. Global Threat of Carbapenem-Resistant Gram-Negative Bacteria. Front. Cell. Infect. Microbiol.; 2022; 12, 823684. [DOI: https://dx.doi.org/10.3389/fcimb.2022.823684] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35372099]
3. Mohsin, M.; Azam, M.; ur Rahman, S.; Esposito, F.; Sellera, F.P.; Monte, D.F.; Cerdeira, L.; Lincopan, N. Genomic background of a colistin-resistant and highly virulent MCR-1-positive Escherichia coli ST6395 from a broiler chicken in Pakistan. Pathog. Dis.; 2019; 77, ftz064. [DOI: https://dx.doi.org/10.1093/femspd/ftz064] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31782775]
4. Zafar, R.; Bashir, S.; Nabi, D.; Arshad, M. Occurrence and quantification of prevalent antibiotics in wastewater samples from Rawalpindi and Islamabad, Pakistan. Sci. Total Environ.; 2021; 764, 142596. [DOI: https://dx.doi.org/10.1016/j.scitotenv.2020.142596]
5. Stange, C.; Sidhu, J.P.S.; Tiehm, A.; Toze, S. Antibiotic resistance and virulence genes in coliform water isolates. Int. J. Hyg. Environ. Health; 2016; 219, pp. 823-831. [DOI: https://dx.doi.org/10.1016/j.ijheh.2016.07.015]
6. Hamelin, K.; Bruant, G.; El-Shaarawi, A.; Hill, S.; Edge, T.A.; Fairbrother, J.; Harel, J.; Maynard, C.; Masson, L.; Brousseau, R. Occurrence of virulence and antimicrobial resistance genes in Escherichia coli isolates from different aquatic ecosystems within the St. Clair River and Detroit River areas. Appl. Environ. Microbiol.; 2007; 73, pp. 477-484. [DOI: https://dx.doi.org/10.1128/AEM.01445-06]
7. Talukdar, P.K.; Rahman, M.; Rahman, M.; Nabi, A.; Islam, Z.; Hoque, M.M.; Endtz, H.P.; Islam, M.A. Antimicrobial Resistance, Virulence Factors and Genetic Diversity of Escherichia coli Isolates from Household Water Supply in Dhaka, Bangladesh. PLoS ONE; 2013; 8, e61090. [DOI: https://dx.doi.org/10.1371/journal.pone.0061090]
8. Ahmed, T.; Zounemat-Kermani, M.; Scholz, M. Climate change, water quality and water-related challenges: A review with focus on Pakistan. Int. J. Environ. Res. Public Health; 2020; 17, 8518. [DOI: https://dx.doi.org/10.3390/ijerph17228518] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33212957]
9. Fida, M.; Li, P.; Wang, Y.; Alam, S.M.K.; Nsabimana, A. Water contamination and human health risks in Pakistan: A review. Expo. Health; 2023; 15, pp. 619-639. [DOI: https://dx.doi.org/10.1007/s12403-022-00512-1]
10. Hall, R.M. Integrons and gene cassettes: Hotspots of diversity in bacterial genomes. Ann. New York Acad. Sci.; 2012; 1267, pp. 71-78. [DOI: https://dx.doi.org/10.1111/j.1749-6632.2012.06588.x]
11. Mazel, D. Integrons: Agents of bacterial evolution. Nat. Rev. Microbiol.; 2006; 4, pp. 608-620. [DOI: https://dx.doi.org/10.1038/nrmicro1462]
12. Gillings, M.R.; Gaze, W.H.; Pruden, A.; Smalla, K.; Tiedje, J.M.; Zhu, Y.G. Using the class 1 integron-integrase gene as a proxy for anthropogenic pollution. ISME J.; 2015; 9, pp. 1269-1279. [DOI: https://dx.doi.org/10.1038/ismej.2014.226] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25500508]
13. Ghaly, T.M.; Chow, L.; Asher, A.J.; Waldron, L.S.; Gillings, M.R. Evolution of class 1 integrons: Mobilization and dispersal via food-borne bacteria. PLoS ONE; 2017; 12, e0179169. [DOI: https://dx.doi.org/10.1371/journal.pone.0179169] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28586403]
14. Bleichenbacher, S.; Stevens, M.J.A.; Zurfluh, K.; Perreten, V.; Endimiani, A.; Stephan, R.; Nüesch-Inderbinen, M. Environmental dissemination of carbapenemase-producing Enterobacteriaceae in rivers in Switzerland. Environ. Pollut.; 2020; 265, 115081. [DOI: https://dx.doi.org/10.1016/j.envpol.2020.115081]
15. Gomi, R.; Matsumura, Y.; Yamamoto, M.; Tanaka, M.; Komakech, A.J.; Matsuda, T.; Harada, H. Genomic surveillance of antimicrobial-resistant Escherichia coli in fecal sludge and sewage in Uganda. Water Res.; 2024; 248, 120830. [DOI: https://dx.doi.org/10.1016/j.watres.2023.120830]
16. Falgenhauer, L.; Schwengers, O.; Schmiedel, J.; Baars, C.; Lambrecht, O.; Heß, S.; Berendonk, T.U.; Falgenhauer, J.; Chakraborty, T.; Imirzalioglu, C. Multidrug-Resistant and Clinically Relevant Gram-Negative Bacteria Are Present in German Surface Waters. Front. Microbiol.; 2019; 10, 2779. [DOI: https://dx.doi.org/10.3389/fmicb.2019.02779]
17. Ahsan, A.; Rehman, T.A.U.; Irshad, H.; Shahzad, M.A.; Siddique, A.; Jamil, A.; Ali, A. Antibiotic resistance pattern and molecular detection of ESBL-associated genes in E. coli from surface and wastewater of Islamabad capital territory, Pakistan. J. Water Health; 2022; 20, pp. 601-609. [DOI: https://dx.doi.org/10.2166/wh.2022.216]
18. Saima, S.; Fiaz, M.; Manzoor, M.; Zafar, R.; Ahmed, I.; Nawaz, U.; Arshad, M. Molecular investigation of antibiotic resistant bacterial strains isolated from wastewater streams in Pakistan. 3 Biotech; 2020; 10, 378. [DOI: https://dx.doi.org/10.1007/s13205-020-02366-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32802720]
19. Khan, G.A.; Berglund, B.; Khan, K.M.; Lindgren, P.E.; Fick, J. Occurrence and Abundance of Antibiotics and Resistance Genes in Rivers, Canal and near Drug Formulation Facilities—A Study in Pakistan. PLoS ONE; 2013; 8, e62712. [DOI: https://dx.doi.org/10.1371/journal.pone.0062712]
20. Yasmin, S.; Karim, A.M.; Lee, S.H.; Zahra, R. Temporal Variation of Meropenem Resistance in E. coli Isolated from Sewage Water in Islamabad, Pakistan. Antibiotics; 2022; 11, 635. [DOI: https://dx.doi.org/10.3390/antibiotics11050635]
21. Zahra, R.; Javeed, S.; Malala, B.; Babenko, D.; Toleman, M.A. Analysis of Escherichia coli STs and resistance mechanisms in sewage from Islamabad, Pakistan indicates a difference in E. coli carriage types between South Asia and Europe. J. Antimicrob. Chemother.; 2018; 73, pp. 1781-1785. [DOI: https://dx.doi.org/10.1093/jac/dky109]
22. Mustafa, S.S.; Batool, R.; Kamran, M.; Javed, H.; Jamil, N. Evaluating the Role of Wastewaters as Reservoirs of Antibiotic-Resistant ESKAPEE Bacteria Using Phenotypic and Molecular Methods. Infect. Drug Resist.; 2022; 15, pp. 5715-5728. [DOI: https://dx.doi.org/10.2147/IDR.S368886]
23. Ahmad, K.; Khattak, F.; Ali, A.; Rahat, S.; Noor, S.; Mahsood, N.; Somayya, R. Carbapenemases and extended-spectrum β-lactamase–producing multidrug-resistant Escherichia coli isolated from retail chicken in peshawar: First report from Pakistan. J. Food Prot.; 2018; 81, pp. 1339-1345. [DOI: https://dx.doi.org/10.4315/0362-028X.JFP-18-045] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30019956]
24. Apanga, P.A.; Ahmed, J.; Tanner, W.; Starcevich, K.; VanDerslice, J.A.; Rehman, U.; Channa, N.; Benson, S.; Garn, J.V. Carbapenem-resistant Enterobacteriaceae in sink drains of 40 healthcare facilities in Sindh, Pakistan: A cross-sectional study. PLoS ONE; 2022; 17, e0263297. [DOI: https://dx.doi.org/10.1371/journal.pone.0263297] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35113948]
25. Nieto-Rosado, M.; Sands, K.; Portal, E.A.R.; Thomson, K.M.; Carvalho, M.J.; Mathias, J.; Milton, R.; Dyer, C.; Akpulu, C.; Boostrom, I. et al. Colonisation of hospital surfaces from low- and middle-income countries by extended spectrum β-lactamase- and carbapenemase-producing bacteria. Nat. Commun.; 2024; 15, 2758. [DOI: https://dx.doi.org/10.1038/s41467-024-46684-z]
26. Khawaja, T.; Mäklin, T.; Kallonen, T.; Gladstone, R.A.; Pöntinen, A.K.; Mero, S.; Thorpe, H.A.; Samuelsen, Ø.; Parkhill, J.; Izhar, M. et al. Deep sequencing of Escherichia coli exposes colonisation diversity and impact of antibiotics in Punjab, Pakistan. Nat. Commun.; 2024; 15, 5196. [DOI: https://dx.doi.org/10.1038/s41467-024-49591-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/38890378]
27. Ilyas, S.; Rasool, M.H.; Arshed, M.J.; Qamar, M.U.; Aslam, B.; Almatroudi, A.; Khurshid, M. The Escherichia coli sequence type 131 harboring extended-spectrum beta-lactamases and carbapenemases genes from poultry birds. Infect. Drug Resist.; 2021; 14, pp. 805-813. [DOI: https://dx.doi.org/10.2147/IDR.S296219]
28. Ilyas, S.; Qamar, M.U.; Rasool, M.H.; Abdulhaq, N.; Nawaz, Z. Multidrug-resistant pathogens isolated from ready-to-eat salads available at a local market in Pakistan. Br. Food J.; 2016; 118, pp. 2068-2075. [DOI: https://dx.doi.org/10.1108/BFJ-02-2016-0058]
29. D’Souza, A.W.; Potter, R.F.; Wallace, M.; Shupe, A.; Patel, S.; Sun, X.; Gul, D.; Kwon, J.H.; Andleeb, S.; Burnham, C.A.D. et al. Spatiotemporal dynamics of multidrug resistant bacteria on intensive care unit surfaces. Nat. Commun.; 2019; 10, 4569. [DOI: https://dx.doi.org/10.1038/s41467-019-12563-1]
30. Zhou, Z.; Alikhan, N.-F.; Mohamed, K.; Fan, Y.; Achtman, M.; Brown, D.; Chattaway, M.; Dallman, T.; Delahay, R.; Kornschober, C. The EnteroBase user’s guide, with case studies on Salmonella transmissions, Yersinia pestis phylogeny, and Escherichia core genomic diversity. Genome Res.; 2020; 30, pp. 138-152. [DOI: https://dx.doi.org/10.1101/gr.251678.119]
31. Chiu, S.K.; Huang, L.Y.; Chen, H.; Tsai, Y.K.; Liou, C.H.; Lin, J.C.; Siu, L.K.; Chang, F.Y.; Yeh, K.M. Roles of ramR and tet(A) mutations in conferring tigecycline resistance in carbapenem-resistant Klebsiella pneumoniae clinical isolates. Antimicrob. Agents Chemother.; 2017; 61, e00391-17. [DOI: https://dx.doi.org/10.1128/AAC.00391-17]
32. Tacconelli, E.; Carrara, E.; Savoldi, A.; Harbarth, S.; Mendelson, M.; Monnet, D.L.; Pulcini, C.; Kahlmeter, G.; Kluytmans, J.; Carmeli, Y. et al. Discovery, research, and development of new antibiotics: The WHO priority list of antibiotic-resistant bacteria and tuberculosis. Lancet Infect. Dis.; 2018; 18, pp. 318-327. [DOI: https://dx.doi.org/10.1016/S1473-3099(17)30753-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29276051]
33. Mills, M.C.; 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, 113143. [DOI: https://dx.doi.org/10.1016/j.envpol.2019.113143] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31541827]
34. Puljko, A.; Barišić, I.; Dekić Rozman, S.; Križanović, S.; Babić, I.; Jelić, M.; Maravić, A.; Udiković-Kolić, N. Molecular epidemiology and mechanisms of carbapenem and colistin resistance in Klebsiella and other Enterobacterales from treated wastewater in Croatia. Environ. Int.; 2024; 185, 108554. [DOI: https://dx.doi.org/10.1016/j.envint.2024.108554]
35. Chaurasia, S.; Sivanandan, S.; Agarwal, R.; Ellis, S.; Sharland, M.; Sankar, M.J. Neonatal sepsis in South Asia: Huge burden and spiralling antimicrobial resistance. Br. Med. J.; 2019; 364, k5314. [DOI: https://dx.doi.org/10.1136/bmj.k5314] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30670451]
36. Malchione, M.D.; Torres, L.M.; Hartley, D.M.; Koch, M.; Goodman, J.L. Carbapenem and colistin resistance in Enterobacteriaceae in Southeast Asia: Review and mapping of emerging and overlapping challenges. Int. J. Antimicrob. Agents; 2019; 54, pp. 381-399. [DOI: https://dx.doi.org/10.1016/j.ijantimicag.2019.07.019]
37. Abdullah, S.; Almusallam, A.; Li, M.; Mahmood, M.S.; Mushtaq, M.A.; Eltai, N.O.; Toleman, M.A.; Mohsin, M. Whole genome-based genetic insights of bla NDM producing clinical E. coli isolates in hospital settings of Pakistan. Microbiol. Spectr.; 2023; 11, e00584-23. [DOI: https://dx.doi.org/10.1128/spectrum.00584-23]
38. Ahlstrom, C.A.; Woksepp, H.; Sandegren, L.; Mohsin, M.; Hasan, B.; Muzyka, D.; Hernandez, J.; Aguirre, F.; Tok, A.; Söderman, J. et al. Genomically diverse carbapenem resistant Enterobacteriaceae from wild birds provide insight into global patterns of spatiotemporal dissemination. Sci. Total Environ.; 2022; 824, 153632. [DOI: https://dx.doi.org/10.1016/j.scitotenv.2022.153632]
39. Habib, A.; Lo, S.; Villageois-Tran, K.; Petitjean, M.; Malik, S.A.; Armand-Lefèvre, L.; Ruppé, E.; Zahra, R. Dissemination of carbapenemase-producing Enterobacterales in the community of Rawalpindi, Pakistan. PLoS ONE; 2022; 17, e0270707. [DOI: https://dx.doi.org/10.1371/journal.pone.0270707] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35802735]
40. Hernández-García, M.; Pérez-Viso, B.; León-Sampedro, R.; Navarro-San Francisco, C.; López-Fresneña, N.; Díaz-Agero, C.; Morosini, M.I.; Ruiz-Garbajosa, P.; Cantón, R. Outbreak of NDM-1+CTX-M-15+DHA-1-producing Klebsiella pneumoniae high-risk clone in Spain owing to an undetectable colonised patient from Pakistan. Int. J. Antimicrob. Agents; 2019; 54, pp. 233-239. [DOI: https://dx.doi.org/10.1016/j.ijantimicag.2019.05.021]
41. Bilal, H.; Zhang, G.; Rehman, T.; Han, J.; Khan, S.; Shafiq, M.; Yang, X.; Yan, Z.; Yang, X. First report of blandm-1 bearing incx3 plasmid in clinically isolated st11 klebsiella pneumoniae from Pakistan. Microorganisms; 2021; 9, 951. [DOI: https://dx.doi.org/10.3390/microorganisms9050951]
42. Chaudhry, T.H.; Aslam, B.; Arshad, M.I.; Alvi, R.F.; Muzammil, S.; Yasmeen, N.; Aslam, M.A.; Khurshid, M.; Rasool, M.H.; Baloch, Z. Emergence of blandm-1 harboring klebsiella pneumoniae st29 and st11 in veterinary settings and waste of Pakistan. Infect. Drug Resist.; 2020; 13, pp. 3033-3043. [DOI: https://dx.doi.org/10.2147/IDR.S248091] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32904734]
43. Lascols, C.; Cherney, B.; Conley, A.B.; Rishishwar, L.; Crawford, M.A.; Morse, S.A.; Fisher, D.J.; Anderson, K.; Hodge, D.R.; Pillai, S.P. et al. Investigation of multidrug-resistant plasmids from carbapenemase-producing Klebsiella pneumoniae clinical isolates from Pakistan. Front. Microbiol.; 2023; 14, 1192097. [DOI: https://dx.doi.org/10.3389/fmicb.2023.1192097]
44. Ejaz, H.; Ahmad, M.; Younas, S.; Junaid, K.; Abosalif, K.O.A.; Abdalla, A.E.; Alameen, A.A.M.; Elamir, M.Y.M.; Bukhari, S.N.A.; Ahmad, N. et al. Molecular epidemiology of extensively-drug resistant acinetobacter baumannii sequence type 2 co-harboring blandm and blaoxa from clinical origin. Infect. Drug Resist.; 2021; 14, pp. 1931-1939. [DOI: https://dx.doi.org/10.2147/IDR.S310478]
45. Saikia, S.; Gogoi, I.; Oloo, A.; Sharma, M.; Puzari, M.; Chetia, P. Co-production of metallo-β-lactamase and OXA-type β-lactamases in carbapenem-resistant Acinetobacter baumannii clinical isolates in North East India. World J. Microbiol. Biotechnol.; 2024; 40, 167. [DOI: https://dx.doi.org/10.1007/s11274-024-03977-1]
46. Jacoby, G.A. AmpC β-lactamases. Clin. Microbiol. Rev.; 2009; 22, pp. 161-182. [DOI: https://dx.doi.org/10.1128/CMR.00036-08] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19136439]
47. Walckenaer, E.; Poirel, L.; Leflon-Guibout, V.; Nordmann, P.; Nicolas-Chanoine, M.H. Genetic and Biochemical Characterization of the Chromosomal Class A β-Lactamases of Raoultella (formerly Klebsiella) planticola and Raoultella ornithinolytica. Antimicrob. Agents Chemother.; 2004; 48, pp. 305-312. [DOI: https://dx.doi.org/10.1128/AAC.48.1.305-312.2004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/14693555]
48. Xu, J.; Zhu, Z.; Chen, Y.; Wang, W.; He, F. The Plasmid-Borne tet(A) Gene Is an Important Factor Causing Tigecycline Resistance in ST11 Carbapenem-Resistant Klebsiella pneumoniae Under Selective Pressure. Front. Microbiol.; 2021; 12, 644949. [DOI: https://dx.doi.org/10.3389/fmicb.2021.644949]
49. Moghnia, O.H.; Al-Sweih, N.A. Whole Genome Sequence Analysis of Multidrug Resistant Escherichia coli and Klebsiella pneumoniae Strains in Kuwait. Microorganisms; 2022; 10, 507. [DOI: https://dx.doi.org/10.3390/microorganisms10030507]
50. Bilal, H.; Khan, M.N.; Rehman, T.; Hameed, M.F.; Yang, X. Antibiotic resistance in Pakistan: A systematic review of past decade. BMC Infect. Dis.; 2021; 21, 244. [DOI: https://dx.doi.org/10.1186/s12879-021-05906-1]
51. Ito, R.; Mustapha, M.M.; Tomich, A.D.; Callaghan, J.D.; McElheny, C.L.; Mettus, R.T.; Shanks, R.M.Q.; Sluis-Cremer, N.; Doi, Y. Widespread fosfomycin resistance in gram-negative bacteria attributable to the chromosomal fosA gene. mBio; 2017; 8, e00749-17. [DOI: https://dx.doi.org/10.1128/mBio.00749-17]
52. Lv, J.; Mohsin, M.; Lei, S.; Srinivas, S.; Wiqar, R.T.; Lin, J.; Feng, Y. Discovery of a mcr-1-bearing plasmid in commensal colistin-resistant Escherichia coli from healthy broilers in Faisalabad, Pakistan. Virulence; 2018; 9, pp. 994-999. [DOI: https://dx.doi.org/10.1080/21505594.2018.1462060]
53. Luo, Q.; Yu, W.; Zhou, K.; Guo, L.; Shen, P.; Lu, H.; Huang, C.; Xu, H.; Xu, S.; Xiao, Y. et al. Molecular epidemiology and colistin resistant mechanism of mcr-positive and mcr-negative clinical isolated Escherichia coli. Front. Microbiol.; 2017; 8, 2262. [DOI: https://dx.doi.org/10.3389/fmicb.2017.02262] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29250039]
54. Yu, Z.; Zhang, Z.; Shi, L.; Hua, S.; Luan, T.; Lin, Q.; Zheng, Z.; Feng, X.; Liu, M.; Li, X. In silico characterization of IncX3 plasmids carrying bla OXA-181 in Enterobacterales. Front. Cell. Infect. Microbiol.; 2022; 12, 988236. [DOI: https://dx.doi.org/10.3389/fcimb.2022.988236]
55. Shropshire, W.C.; Konovalova, A.; McDaneld, P.; Gohel, M.; Strope, B.; Sahasrabhojane, P.; Tran, C.N.; Greenberg, D.; Kim, J.; Zhan, X. et al. Systematic Analysis of Mobile Genetic Elements Mediating β-Lactamase Gene Amplification in Noncarbapenemase-Producing Carbapenem-Resistant Enterobacterales Bloodstream Infections. mSystems; 2022; 7, e0047622. [DOI: https://dx.doi.org/10.1128/msystems.00476-22]
56. Wang, L.; Zhu, M.; Yan, C.; Zhang, Y.; He, X.; Wu, L.; Xu, J.; Lu, J.; Bao, Q.; Hu, Y. et al. Class 1 integrons and multiple mobile genetic elements in clinical isolates of the Klebsiella pneumoniae complex from a tertiary hospital in eastern China. Front. Microbiol.; 2023; 14, 985102. [DOI: https://dx.doi.org/10.3389/fmicb.2023.985102]
57. Godziszewska, J.; Kuliñska, A.; Jagura-Burdzy, G. MobC of conjugative RA3 plasmid from IncU group autoregulates the expression of bicistronic mobC-nic operon and stimulates conjugative transfer. BMC Microbiol.; 2014; 14, 235. [DOI: https://dx.doi.org/10.1186/s12866-014-0235-1]
58. Bashir, S.; Haque, A.; Sarwar, Y.; Ali, A.; Anwar, M.I. Virulence profile of different phylogenetic groups of locally isolated community acquired uropathogenic E. coli from Faisalabad region of Pakistan. Ann. Clin. Microbiol. Antimicrob.; 2012; 11, 23. [DOI: https://dx.doi.org/10.1186/1476-0711-11-23]
59. Sambrook, J.; Fritsch, E.F.; Maniatis, T. Molecular Cloning: A Laboratory Manual; 2nd ed. Cold Spring Harbor Laboratory Press: Cold Spring Harbor, NY, USA, 1989.
60. Clinical and Laboratory Standards Institute (CLSI). Performance Standards for Antimicrobial Susceptibility Testing; CLSI Supplement M100-S30 CLSI: Wayne, PA, USA, 2020.
61. Zheng, Z.; Gorden, P.J.; Xia, X.; Zheng, Y.; Li, G. Whole-genome analysis of Klebsiella pneumoniae from bovine mastitis milk in the U.S. Environ. Microbiol.; 2022; 24, pp. 1183-1199. [DOI: https://dx.doi.org/10.1111/1462-2920.15721]
62. Wick, R.R.; Judd, L.M.; Gorrie, C.L.; Holt, K.E. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput. Biol.; 2017; 13, e1005595. [DOI: https://dx.doi.org/10.1371/journal.pcbi.1005595]
63. Chen, X.; Liu, W.; Li, H.; Yan, S.; Jiang, F.; Cai, W.; Li, G. Whole genome sequencing analysis of avian pathogenic Escherichia coli from China. Vet. Microbiol.; 2021; 259, 109158. [DOI: https://dx.doi.org/10.1016/j.vetmic.2021.109158]
64. Chklovski, A.; Parks, D.H.; Woodcroft, B.J.; Tyson, G.W. CheckM2: A rapid, scalable and accurate tool for assessing microbial genome quality using machine learning. Nat. Methods; 2023; 20, pp. 1203-1212. [DOI: https://dx.doi.org/10.1038/s41592-023-01940-w] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37500759]
65. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality assessment tool for genome assemblies. Bioinformatics; 2013; 29, pp. 1072-1075. [DOI: https://dx.doi.org/10.1093/bioinformatics/btt086] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23422339]
66. Lumpe, J.; Gumbleton, L.; Gorzalski, A.; Libuit, K.; Varghese, V.; Lloyd, T.; Tadros, F.; Arsimendi, T.; Wagner, E.; Stephens, C. et al. GAMBIT (Genomic Approximation Method for Bacterial Identification and Tracking): A methodology to rapidly leverage whole genome sequencing of bacterial isolates for clinical identification. PLoS ONE; 2023; 18, e0277575. [DOI: https://dx.doi.org/10.1371/journal.pone.0277575] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36795668]
67. Jain, C.; Rodriguez-R, L.M.; Phillippy, A.M.; Konstantinidis, K.T.; Aluru, S. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun.; 2018; 9, 5114. [DOI: https://dx.doi.org/10.1038/s41467-018-07641-9] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30504855]
68. Seemann, T. Prokka: Rapid prokaryotic genome annotation. Bioinformatics; 2014; 30, pp. 2068-2069. [DOI: https://dx.doi.org/10.1093/bioinformatics/btu153]
69. Waters, N.R.; Abram, F.; Brennan, F.; Holmes, A.; Pritchard, L. Easy phylotyping of Escherichia coli via the EzClermont web app and command-line tool. Access Microbiol.; 2020; 2, e000143. [DOI: https://dx.doi.org/10.1099/acmi.0.000143]
70. Jolley, K.A.; Maiden, M.C.J. BIGSdb: Scalable analysis of bacterial genome variation at the population level. BMC Bioinform.; 2010; 11, 595. [DOI: https://dx.doi.org/10.1186/1471-2105-11-595]
71. Feldgarden, M.; Brover, V.; Gonzalez-Escalona, N.; Frye, J.G.; Haendiges, J.; Haft, D.H.; Hoffmann, M.; Pettengill, J.B.; Prasad, A.B.; Tillman, G.E. et al. AMRFinderPlus and the Reference Gene Catalog facilitate examination of the genomic links among antimicrobial resistance, stress response, and virulence. Sci. Rep.; 2021; 11, 12728. [DOI: https://dx.doi.org/10.1038/s41598-021-91456-0]
72. Haft, D.H.; DiCuccio, M.; Badretdin, A.; Brover, V.; Chetvernin, V.; O’Neill, K.; Li, W.; Chitsaz, F.; Derbyshire, M.K.; Gonzales, N.R. et al. RefSeq: An update on prokaryotic genome annotation and curation. Nucleic Acids Res.; 2018; 46, pp. D851-D860. [DOI: https://dx.doi.org/10.1093/nar/gkx1068]
73. Néron, B.; Littner, E.; Haudiquet, M.; Perrin, A.; Cury, J.; Rocha, E.P.C. IntegronFinder 2.0: Identification and Analysis of Integrons across Bacteria, with a Focus on Antibiotic Resistance in Klebsiella. Microorganisms; 2022; 10, 700. [DOI: https://dx.doi.org/10.3390/microorganisms10040700]
74. Ghaly, T.M.; Gillings, M.R.; Rajabal, V.; Paulsen, I.T.; Tetu, S.G. Horizontal gene transfer in plant microbiomes: Integrons as hotspots for cross-species gene exchange. Front. Microbiol.; 2024; 15, 1338026. [DOI: https://dx.doi.org/10.3389/fmicb.2024.1338026] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/38741746]
75. Tian, R.; Zhou, J.; Imanian, B. PlasmidHunter: Accurate and fast prediction of plasmid sequences using gene content profile and machine learning. Brief. Bioinform.; 2024; 25, bbae322. [DOI: https://dx.doi.org/10.1093/bib/bbae322] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/38960405]
76. Alikhan, N.F.; Petty, N.K.; Ben Zakour, N.L.; Beatson, S.A. BLAST Ring Image Generator (BRIG): Simple prokaryote genome comparisons. BMC Genom.; 2011; 12, 402. [DOI: https://dx.doi.org/10.1186/1471-2164-12-402] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21824423]
77. Robertson, J.; Nash, J.H.E. MOB-suite: Software tools for clustering, reconstruction and typing of plasmids from draft assemblies. Microb. Genom.; 2018; 4, e000206. [DOI: https://dx.doi.org/10.1099/mgen.0.000206]
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Hu, Xiao 2
; Saxena, Anugrah 2
; Mou, Kathy 2 ; Shen, Huigang 2
; Hazrat, Ali 3 ; Muhammad Afzal Ghauri 3 ; Sarwar, Yasra 3 ; Ali, Aamir 3
; Li, Ganwu 2 1 National Institute for Biotechnology and Genetic Engineering College, Pakistan Institute of Engineering and Applied Sciences (NIBGE-C, PIEAS), Faisalabad 38000, Punjab, Pakistan;
2 Department of Veterinary Diagnostic and Production Animal Medicine, College of Veterinary Medicine, Iowa State University, Ames, IA 50011, USA;
3 National Institute for Biotechnology and Genetic Engineering College, Pakistan Institute of Engineering and Applied Sciences (NIBGE-C, PIEAS), Faisalabad 38000, Punjab, Pakistan;