- Open Access
Molecular epidemiology and antimicrobial resistance features of Acinetobacter baumannii clinical isolates from Pakistan
Annals of Clinical Microbiology and Antimicrobials volume 19, Article number: 2 (2020)
Acinetobacter baumannii is a Gram-negative opportunistic pathogen with a notorious reputation of being resistant to antimicrobial agents. The capability of A. baumannii to persist and disseminate between healthcare settings has raised a major concern worldwide.
Our study investigated the antibiotic resistance features and molecular epidemiology of 52 clinical isolates of A. baumannii collected in Pakistan between 2013 and 2015. Antimicrobial susceptibility patterns were determined by the agar disc diffusion method. Comparative sequence analyses of the ampC and blaOXA-51-like alleles were used to assign the isolates into clusters. The whole genomes of 25 representative isolates were sequenced using the MiSeq Desktop Sequencer. Free online applications were used to determine the phylogeny of genomic sequences, retrieve the multilocus sequence types (ST), and detect acquired antimicrobial resistance genes.
Overall, the isolates were grouped into 7 clusters and 3 sporadic isolates. The largest cluster, Ab-Pak-cluster-1 (blaOXA-66 and ISAba1-ampC-19) included 24 isolates, belonged to ST2 and International clone (IC) II, and was distributed between two geographical far-off cities, Lahore and Peshawar. Ab-Pak-clusters-2 (blaOXA-66, ISAba1-ampC-2), and -3 (blaOXA-66, ISAba1-ampC-20) and the individual isolate Ab-Pak-Lah-01 (ISAba1-blaOXA-66, ISAba1-ampC-2) were also assigned to ST2 and IC II. On the other hand, Ab-Pak-clusters-4 (blaOXA-69, ampC-1), -5 (blaOXA-69, ISAba1-ampC-78), and -6A (blaOXA-371, ISAba1-ampC-3) belonged to ST1, while Ab-Pak-cluster-6B (blaOXA-371, ISAba1-ampC-8) belonged to ST1106, with both ST1 and ST1106 being members of IC I. Five isolates belonged to Ab-Pak-cluster-7 (blaOXA-65, ampC-43). This cluster corresponded to ST158, showed a well-delineated position on the genomic phylogenetic tree, and was equipped with several antimicrobial resistance genes including blaOXA-23 and blaGES-11.
Our study detected the occurrence of 7 clusters of A. baumannii in Pakistan. Altogether, 6/7 of the clusters and 45/52 (86.5%) of the isolates belonged to IC I (n = 9) or II (n = 36), making Pakistan no exception to the global domination of these two clones. The onset of ST158 in Pakistan marked a geographical dispersal of this clone beyond the Middle East and brought up the need for a detailed characterization.
Acinetobacter baumannii is an important opportunistic pathogen that has increasingly been reported worldwide. The World Health Organization (WHO) has recently classified carbapenem resistant A. baumannii as a top priority organism for research and development of new antibiotics against antibiotic-resistant bacteria . The occurrence of A. baumannii in Pakistan was first reported in 2004, where it counted for 4.4% of the Gram-negative bacilli in a collection of 812 isolates obtained from a variety of clinical samples in Rawalpindi in 2002 . Importantly, 72% of the reported A. baumannii isolates were found to be extended-spectrum beta-lactamase (ESBL) producers. Later, eight clusters of carbapenem-resistant isolates of A. baumannii were detected in two intensive care units in Karachi, Pakistan, between November 2007 and August 2008 . Concurrently, A. baumannii was the most frequent bacterial pathogen among a collection of 50 carbapenem-resistant clinical isolates obtained in Rawalpindi in 2009 . According to Kaleem et al., twenty-seven (84%) of the carbapenem-resistant A. baumannii isolates were metallo-beta-lactamase positive by the E-test strip method.
Recently, A. baumannii was isolated in 18% of patients with ventilator-associated pneumonia, making it the 4th most frequent pathogen, at the Khyber Teaching Hospital in Peshawar in 2013 . Notably, A. baumannii was declared as the most resistant pathogen identified in this study. A. baumannii also ranked the second most commonly detected pathogen (24.7%) among a collection of carbapenemase producers obtained from one hospital in Lahore between 2015 and 2016 . Furthermore, A. baumannii was the main isolated bacterium (39.8%) among a total of 113 isolates recovered from 80 ventilator-supported patients at the Foundation Hospital, Rawalpindi, in 2016 .
Nevertheless, only few of these studies have plotted the epidemiology of A. baumannii in Pakistan onto a global map. Most of the 2007–2008 hospital-acquired carbapenem-resistant A. baumannii belonged to either International clone (IC) I or II . Similarly, 1/17 and 7/17 clinical isolates of A. baumannii, obtained in 2008, were linked to ICs I and II, respectively . Two isolates collected in Norway in 2009 with a history of import from Pakistan were assigned to sequence types 2 (ST2) and ST15, using the multilocus sequence typing (MLST) scheme of Pasteur Institute [9, 10]. ST2 is a key member of IC II while ST15 has often been reported in South America [11, 12]. Our search in the A. baumannii MLST databases yielded records of 4 isolates from Pakistan (last accessed in 11 September 2019). The 4 isolates were entered in 2019, which might explain why we could not find published literature on them. Two of these isolates belonged to ST1310 (8, 1, 5, 1, 6, 2, 3) or ST1327 (13, 1, 5, 26, 7, 1, 29), for which there was no linkage to known ICs. The other two isolates belonged to ST1106 (2, 1, 2, 1, 5, 1, 1), a member of clonal complex CC1 corresponding to IC I .
The aim of this study was to investigate the molecular epidemiology and antibiotic resistance features of 52 A. baumannii clinical isolates collected in Pakistan between 2013 and 2015. The clonality among our isolates was crossed over the global population of A. baumannii.
Materials and methods
Two collections of A. baumannii clinical isolates were included in this study. The isolates, each recovered from one patient, were collected on a consecutive basis regardless of their antimicrobial susceptibility features (Table 1). The first collection (n = 16) was obtained between February and September 2013 at the Combined Military Hospital in Lahore (CMH Lahore), while the second one (n = 36) was obtained between February and November 2015 at the Combined Military Hospital in Peshawar (CMH Peshawar). The geographic distance between the two hospitals is about 520 kilometres (Additional file 1: Figure S1). The isolates were cultured from pus (n = 17), respiratory secretions (n = 27), intravenous catheter tip (n = 5), Foley catheter tip (n = 2), and bile (n = 1). The clinical samples were processed according to standard in-house culture protocols. When needed, isolates were revived from frozen glycerol stocks onto Brain Heart Infusion agar plates or slants (Oxoid, Basingstoke, United Kingdom). Species identification was first determined using API 20 NE (BioMerieux, France) and later confirmed by detecting occurrence of the intrinsic blaOXA-51-like gene in all the isolates . The identification was further verified by partial rpoB (zone 1, 352 bp) gene sequence analysis in 25 isolates selected for whole-genome sequencing .
Antimicrobial susceptibility testing
Susceptibility of the isolates to ampicillin/sulbactam, piperacillin/tazobactam, cefepime, cefotaxime, ceftazidime, ceftriaxone, meropenem, imipenem, gentamicin, amikacin, tobramycin, ciprofloxacin, levofloxacin, trimethoprim/sulfamethoxazole, and minocycline was investigated by the agar disc diffusion method, using discs from Oxoid (Basingstoke, United Kingdom). The broth microdilution method was used for doxycycline (Pfizer Global pharmaceuticals), polymyxin B (Glaxosmith Kline pharmaceuticals), and colistin (Forest pharmaceuticals). The tests were performed and susceptibility patterns were interpreted following the guidelines of the Clinical and Laboratory Standards Institute (CLSI) .
The isolates were typed using two single-locus molecular schemes based on the allelic identity of the A. baumannii-intrinsic ampC and blaOXA-51-like genes [16, 17]. PCR amplifications of ampC and blaOXA-51-like were performed using in-house designed primers (Additional file 2: Table S1) and followed by Sanger sequencing of the amplicons, as previously described [16, 17]. This approach was able to detect the occurrence of insertion sequence (IS) elements, such as ISAba1, in the bordering regions of ampC and/or blaOXA-51-like.
Whole-genome sequence analyses
Twenty-five isolates were selected for whole-genome sequencing based on their antimicrobial resistance patterns and strain typing results. Sequencing was performed using the MiSeq Desktop Sequencer and MiSeq Reagent Kit v3 (Illumina, San Diego, CA, USA). DNA preparation, library construction, and genome sequencing were done according to the manufacturers’ instructions. Sequence data were assembled and analysed using the CLC genomics workbench (v7.0.4; CLC bio, Aarhus, Denmark). The web server Reference sequence Alignment-based Phylogeny builder (REALPHY) was used to construct a phylogenetic tree based on multiple alignments of the genomic sequence data . Sequence reads were mapped to the genomes of three well-known reference strains: AYE (IC I, GenBank: NC_010410.1), ACICU (IC II, GenBank: NC_010611.1), and ATCC 17978 (sporadic ST437, GenBank: NZ_CP018664.1).
An In-silico search was performed to detect the existence of 24 A. baumannii plasmid-borne replicase genes, according to the A. baumannii PCR-based replicon typing (AB-PBRT) scheme and a number of other studies on A. baumannii plasmids [19,20,21,22]. The MLST web-based search engine, hosted by the Center for Genomic Epidemiology in Denmark (http://www.genomicepidemiology.org/), was used to assign the isolates into STs according to the Institute Pasteur’s MLST scheme (http://www.pasteur.fr/mlst) . The occurrence of acquired antimicrobial resistance genes was detected using the ResFinder service, also hosted by the Center for Genomic Epidemiology in Denmark . The occurrence of resistance genes was verified, and genetic surroundings were annotated based on the yields of nucleotide similarities obtained using the Basic Local Alignment Search Tool (http://blast.ncbi.nlm.nih.gov/Blast.cgi) against the “Nucleotide collection (nr/nt)” and/or “Whole-genome shotgun contigs (wgs)” databases . The presence of neighbouring IS elements was detected using the ISfinder online application .
Nucleotide sequence accession numbers
Draft genome sequences of all the isolates were deposited in the DDBJ/EMBL/GenBank database under the BioProject accession number: PRJNA482499. The versions described in this paper are QQPR01000000, SMUA01000000, QQPS01000000, QQPT01000000, SMUB01000000, SMUC01000000, SMUD01000000, QQPU01000000, QQPV01000000, QQPW01000000, SMUE01000000, SMUF01000000, SMUG01000000, SMUQ01000000, QQPX01000000, SMUR01000000, QQPY01000000, SMUH01000000, SMUI01000000, QQPZ01000000, SMUJ01000000, QQQA01000000, QQQB01000000, SMUK01000000, and QQQC01000000. The new alleles ampC-78 and ampC-79 were deposited in the ampC database hosted at the A. baumannii MLST website (http://pubmlst.org/abaumannii/) and in the DDBJ/EMBL/GenBank database under the accession numbers MK634301.1 (blaADC-191) and MK690541.1 (blaADC-91), respectively.
Results and discussion
Antimicrobial resistance patterns
The isolates showed resistance rates of 100% (52/52) to ampicillin/sulbactam, piperacillin/tazobactam, cefepime, cefotaxime, ceftazidime, ceftriaxone, meropenem, gentamicin, tobramycin, ciprofloxacin, levofloxacin, and trimethoprim/sulfamethoxazole (Table 1 and Additional file 3: Table S2). In addition, 98.1% (51/52) of the isolates were resistant to amikacin and imipenem, 32.7% (17/52) to minocycline, 30.8% (16/52) to doxycycline, and 0% (0/52) to colistin and polymyxin B (Table 1 and Additional file 3: Table S2). All the isolates were resistant to ≥ 3 classes of antibiotics and were accordingly defined as multidrug-resistant . The results were in line with previous studies describing extensive occurrence of multidrug-resistant strains of A. baumannii in Pakistan [28,29,30]. However, our results could be partially biased by the predominance of few genetically related clusters of isolates, as described in the next paragraph.
Comparative sequence analysis of the blaOXA-51-like and ampC loci enabled us to group 49/52 of the isolates into seven clusters, designated Ab-Pak-cluster-1 to -7 (Table 2). The remaining 3 isolates (Ab-Pak-Lah-01, Ab-Pak-Pesh-04 and Ab-Pak-Pesh-22) were considered as individual strains. Ab-Pak-cluster-1 (blaOXA-66, ISAba1-ampC-19) was the largest cluster with 24 isolates distributing between Lahore, 2013 (n = 14) and Peshawar, 2015 (n = 10). The second largest cluster, Ab-Pak-cluster-2 (blaOXA-66, ISAba1-ampC-2), included 7 isolates, which were all detected in Peshawar, 2015. Ab-Pak-cluster-3 (blaOXA-66, ISAba1-ampC-20) included 4 isolates and was also exclusive to Peshawar, 2015. Individual strain Ab-Pak-Lah-01 (ISAba1-blaOXA-66, ISAba1-ampC-2) was the only isolate in our collection carrying insertion sequence ISAba1 upstream its intrinsic blaOXA-51-like gene (Table 2). Based on their blaOXA-51-like and ampC alleles, Ab-Pak-clusters-1, -2 and -3, and individual strain Ab-Pak-Lah-01 belonged to A. baumannii IC II [11, 16, 17].
Ab-Pak-cluster-4 (blaOXA-69, ampC-1) and Ab-Pak-cluster-5 (blaOXA-69, ISAba1-ampC-78) included 3 and 2 isolates, respectively, and were both collected in Peshawar, 2015. The ampC-78 allele had a novel sequence according to the nucleotide database of ampC alleles in A. baumannii (https://pubmlst.org/abaumannii/). Ab-Pak-cluster-6 (blaOXA-371) included 4 isolates scattering between Lahore, 2013 (1 isolate) and Peshawar, 2015 (3 isolates). The 4 isolates were equally divided into two sub-clusters, namely Ab-Pak-sub-cluster-6A (blaOXA-371, ISAba1-ampC-3) and -6B (blaOXA-371, ISAba1-ampC-8). The nucleotide identity of ampC-8 (1167 bp) is 100% identical to ampC-3 (1152 bp) apart from having a duplicated segment of 15 nucleotides (Additional file 4: File S1). The blaOXA-371 allele was only one nucleotide different in comparison to blaOXA-69 (GenBank accession numbers NG_049662.1 and NG_049809.1). Ab-Pak-clusters-4, -5, -6A, and -6B could conceivably be classified under the umbrella of IC I [16, 17, 31].
Ab-Pak-cluster-7 (blaOXA-65, ampC-43) included 5 isolates that were collected in Peshawar, 2015. The blaOXA-65 allele in Ab-Pak-cluster-7 encoded for ß-lactamase OXA-65 (GenBank: WP_001021782.1). However, this allele had 3 synonymous nucleotide substitutions in comparison to the first GenBank-deposited allele for OXA-65 (GenBank: NG_049805.1). This deviation revealed an inherent shortage in the current numbering system of the blaOXA genes [32, 33]. In this regard, we would restate our view that alleles of well-defined bacterial genes should be numbered based on their nucleotide identifies rather than their amino acid sequences . Lastly, the individual strains Ab-Pak-Pesh-04 (blaOXA-64, ampC-25) and Ab-Pak-Pesh-22 (blaOXA-68, ISAba1-ampC-79) were assigned into CC25 and CC23, respectively [34, 35].
Whole-genome phylogenetic tree
Among the 25/52 isolates selected for whole genome sequence analysis, Ab-Pak-clusters-1, -2, -3, -4, -5, -6A, -6B, -7 were represented by 6 (showing 3 different antimicrobial resistance patterns)/24, 3/7, 2/4, 2/3, 2/2, 2/2, 2/2, and 3/5 isolates, respectively. The three individual isolates were also included. The genome assembly features were presented in Additional file 5: Table S3. The whole-genome phylogenetic tree showed a branching pattern indorsing the assembly of Ab-Pak-cludters-2, -3, -4, -5 and -7 (Fig. 1). The distinction between Ab-Pak-cludters-6A and -6B was also confirmed. The six isolates of Ab-Pak-cluster-1 were distributed on two detached branches and a few sub-branches. Such a high discriminatory power of whole-genome phylogenetic analyses, compared to loci-focused strain typing approaches, has been noted in previous studies [36, 37]. Nevertheless, the internal splits in Ab-Pak-cluster-1 were most likely related to a long-standing presence of this cluster. In addition, the tree showed a noticeably well-demarcated positioning of Ab-Pak-cluster-7 far from ICs I and II. Sequencing the whole genome of more isolates, especially from Ab-Pak-cluster-1, would have increased our capacity to confirm or exclude transmission events between patients . However, such a limitation is generally tolerable when resources are limited.
Plasmid replicase genes
The isolates carried between 1 and 4 plasmid replicase genes (Additional file 6: Table S4). The most commonly detected genes were repAci6 (16/25) and repAci1/repAci2 (15/25), which was in line with earlier studies elsewhere . Interestingly, repAci9 was detected in 12/25 of the isolates. The repAci9-positive plasmids ranged in size between 10,000 and 16,000 base pairs (bp), similarly to what has commonly been described (GenBank accession numbers: LT594096.1, CP024577.1, CP040088.1, CP035673.1, and AY541809.1). However, the individual isolate Ab-Pak-Pesh-22 carried a relatively large repAci9-positive plasmid, with the size of 31,442 bp (GenBank: SMUR01000022.1), on which the tet(39) tetracycline resistance operon was detected.
Other replicase genes, such as repAci4 and the replicase genes of p2ABSDF and p4AYE, were present in 2/25, 5/25, and 4/25 of the isolates, respectively. Two isolates, Ab-Pak-Pesh-20 and Ab-Pak-Pesh-27, from Ab-Pak-cluster-4 carried a replicase gene with novel nucleotide identity, designated repAci24, showing only 58% nucleotide identity to the replicase gene of p3ABSDF (GenBank: CU468233.1). The three Ab-Pak-cluster-7 isolates carried another novel replicase gene, designated repAci25, which will be further discussed below. The occurrence of partial sequences of some replicase genes was noted in the genomes of few isolates (Additional file 6: Table S4). Our in silico analysis was not able to confirm or exclude if these incomplete sequences represented or belonged to actual plasmids.
Multilocus sequence types
The selected representatives from Ab-Pak-clusters-1, -2, and -3 and individual strain Ab-Pak-Lah-01 were all assigned to ST2 (2, 2, 2, 2, 2, 2, 2), as shown in Table 2. Ab-Pak-clusters-4, -5, and -6A belonged to ST1 (1, 1, 1, 1, 5, 1, 1), whereas Ab-Pak-cluster-6B belonged to ST1106 (2, 1, 2, 1, 5, 1, 1). ST1 and ST1106 shared identical alleles in 5/7 of the loci. Ab-Pak-cluster-7 was assigned to ST158 (41, 42, 13, 1, 5, 4, 14). The individual isolates Ab-Pak-Pesh-04 and Ab-Pak-Pesh-22 belonged to ST25 (3, 3, 2, 4, 7, 2, 4) and ST23 (1, 3, 10, 1, 4, 4, 4), respectively. Overall, the MLST results were consistent with the two single-loci, blaOXA-51-like and ampC, typing results and with the whole-genome-based phylogenetic analysis.
Ab-Pak-cluster-7, a new member of clone ST158
Searching for homologous sequences in the GenBank databases enabled us to detect two A. baumannii strains, K50 (GenBank: OHJL00000000.1) and AA-014 (GenBank: AMGA00000000.1), carrying identical ampC and blaOXA-51-like alleles as of Ab-Pak-cluster-7. Both K50 and AA-014 belonged to ST158 [40, 41]. K50 was recovered from a hospitalized patient in Kuwait in 2008, while AA-014 was collected from a wound specimen in Iraq in 2008. Additional two isolates, 2226C and ABC002, were assigned into ST158 according to the MLST database (https://pubmlst.org/abaumannii/). 2226C was isolated in Turkey in 2009 while ABC002 was collected in Egypt in 2011 (https://pubmlst.org/). ST158 isolates were also reported in Lebanon between 2011 and 2013 . The recognition of Ab-Pak-cluster-7 in Pakistan has then expanded the geographic distribution of ST158 beyond the Middle East [43,44,45,46].
Antimicrobial resistance genes
A variety of antimicrobial resistance genes were detected in the genomes of the 25 isolates selected for whole genome sequencing (Additional file 7: Table S5). Class D ß-lactamase gene blaOXA-23 was present in 24/25 isolates. blaOXA-23 was located either in Tn2006 (12 isolates) or Tn2008 (12 isolates). Some isolates, such as Ab-Pak-Lah-04 and Ab-Pak-Lah-08 from Ab-Pak-cluster-1, carried three copies of blaOXA-23. Of note, only one copy was intact while the other two were truncated at their 5′ extremities (data not published). A similar scenario was detected in strain AbPK1 (GenBank: CP024576). Interestingly, AbPK1 was collected in Pakistan in 2012. It also belonged to ST2 and had the ampC-19 and blaOXA-66 alleles . Although it was isolated from infected animals, AbPK1 had a strong genotypic linkage to Ab-Pak-cluster-1.
The individual strain Ab-Pak-Pesh-22 was resistant to meropenem but susceptible to imipenem. Interestingly, Ab-Pak-Pesh-22 did not carry blaOXA-23 or any other commonly known mechanism for carbapenem resistance. On the other hand, Ab-Pak-Lah-01 was equipped with two mechanisms conferring carbapenem resistance, an acquired blaOXA-23 gene and an intrinsic blaOXA-66 gene proceeded by ISAba1 . None of our isolates carried the class D ß-lactamase blaOXA-24-like or blaOXA-58-like genes.
Overall, our results were in line with previous studies nominating OXA-23 as the most commonly detected group of acquired Class D β-lactamases in A. baumannii both in the region and throughout the world [49,50,51]. The occurrence of ISAba1 upstream of blaOXA-23, blaOXA-66, ampC-2, ampC-3, ampC-8, ampC-19, ampC-20, ampC-78, or ampC-79 has the potential to overexpress these genes, which might subsequently confer resistance to carbapenems (blaOXA) or ceftazidime (ampC) as previously proposed [48, 52, 53].
Class A extended-spectrum ß-lactamase genes were detected in 14/25 of the isolates, among which blaPER-1, blaTEM-1D, blaGES-11 and blaPER-7 were carried by 5, 5, 3, and 1 isolate, respectively. There was no co-existence of more than one class A ß-lactamase gene per isolate. All the isolates (25/25) carried aminoglycoside resistance genes. Each isolate carried between 2 to 7 aminoglycoside resistance genes. The occurrence rates ranged from 20/25 for aphA6a and strB, 16/25 for strA, 13/25 for aacC1, 8/25 for aphA1b, 7/25 for aadA1, 4/25 for aphA6b, aadB, and aacA4, and 1/25 for aacC2. The 16S rRNA methylase armA gene, commonly conferring high levels of resistance to aminoglycosides , was detected in 3 isolates, namely Ab-Pak-Pesh-04, Ab-Pak-Pesh-31, and Ab-Pak-Pesh-36. The inhibition zone diameters of gentamicin 10 μg, amikacin 30 μg, and tobramycin 10 μg for these armA-positive isolates did not show a unique resistance pattern compared to the armA-negative isolates (Additional file 3: Table S2).
The macrolide resistance 2′-phophotransferase msr(E)-mph(E) operon was present in 15/25 isolates. The sulphonamide resistance genes sul2 and sul1 were detected in 14/25 and 11/25 isolates, respectively, with 3 isolates having both sul2 and sul1. The tetracycline resistance genes tet(B), tet(A), and tet(39) were detected in 9/25, 2/25, and 1/25 isolates, respectively, with no cases of overlapping. Other resistance genes, such as floR, catA1, cmlA5, and catB8 (chloramphenicol resistance), dfrA7 and dfrA1 (trimethoprim resistance), and arr-2 (rifampicin resistance) had low rates of occurrence ranging between 1/24 and 4/24 (Additional file 7: Table S5).
Antimicrobial resistance features of Ab-Pak-cluster-7
Ab-Pak-cluster-7 was equipped with 8 acquired antimicrobial resistance genes, namely: blaOXA-23, blaGES-11, aphA6a, aacA4, sul1, drfA7, msr(E) and mph(E). The blaOXA-23 gene was located on transposon Tn2008 with genetic surroundings identical to plasmid pK50a carried by isolate K50 . Four genes, aacA4, drfA7, blaGES-11 and sul1, were part of a truncated class one integron (Fig. 2). The integron was bounded by two miniature inverted-repeat transposable elements (MITE), forming a mobile element of 7486 bp. The insertion of this mobile element created a characteristic target site duplication of 5-bp . The same element was previously found on plasmids pK50a , p1AB5075 (CP008707.1), and pAb8098 (KY022424.1). It was also detected in the whole genome sequence of isolates 428 (MDTS01000040.1) and IS-251 (AMEJ01000018.1). The aphA6a gene in Ab-Pak-cluster-7 was located on transposon TnaphA6, inserted at the same location as previously described for plasmids pAb-G7-2 (KF669606.1) and pAbPK1b (CP024578.1). Of note, TnaphA6 was missing in isolate K50.
The msr(E)-mph(E) macrolide resistance operon in Ab-Pak-cluster-7 was surrounded by two pdif sites creating a module with a probable movability mediated by the XerC–XerD system . The msr(E)-mph(E) operon was co-located on a contig carrying the repAci25 gene, encoding for novel plasmid replication initiator protein from the Rep_3 protein family (pfam01051) and superfamily (cl19398). repAci25 had a nucleotide identity of 91% to repAci4 (GenBank: GU978998.1). A close allele to repAci25, with 99.6% nucleotide identity, was detected on plasmids p597A-14.8 (GenBank: CP033871.1), pAb825_36 (GenBank: MG100202.1), pAb244_7 (GenBank: MG520098.1), and pAb242_9 (GenBank: KY984045.1). Among other Acinetobacter species, repAci25 was detected in the whole genome sequence of Acinetobacter bereziniae (100% nucleotide identity, GenBank: CDEL01000266.1) and Acinetobacter radioresistens (99.7% nucleotide identity, GenBank: PXJE01000090.1). Interestingly, plasmid pALWVS1.1 in Acinetobacter lwoffii strain VS15 had a gene showing 96.7% nucleotide identity to repAci25 (GenBank: KX426232.1). The age of the permafrost sediment, from which VS15 was obtained, was estimated to be around 2 to 3 million years .
Our study demonstrated that uncomplicated sequence typing schemes, based on comparative analysis of one or two intrinsic loci, such as ampC and blaOXA-51-like, could be a practical approach for rapid grouping of bacterial clinical isolates, such as A. baumannii. However, further validation of the results is commonly needed. The occurrence of a 100% rate of multidrug-resistant strains is alarming and worth a rapid action plan, including regular follow ups and urgent management procedures. The study enabled us to detect seven clusters of A. baumannii prevailing in two clinical settings in Pakistan. Two of these clusters, Ab-Pak-cluster-1 and -6, lasted for more than 2 years and were able to spread between Lahore and Peshawar. The predominance of IC II in Pakistan was in line with the intensive circulation of this clone worldwide (11). The frequent occurrence of isolates belonging to IC I underlined that this clone is still a key trouble-maker in several parts of the world, including Pakistan. Significantly, Pakistan and the Middle East could be a reservoir for under-detected clones of carbapenem-resistant A. baumannii, including the blaOXA-23- and blaGES-11-positive ST158. The clinical significance and virulence features of ST158, represented by Ab-Pak-cluster-7 in Pakistan, is worth further investigations.
Availability of data and materials
World Health Organisation. Global priority list of antibiotic-resistant bacteria to guide research, discovery, and development of new antibiotics. Geneva: World Health Organisation; 2017.
Ali AM, Rafi S, Qureshi AH. Frequency of extended spectrum beta lactamase producing gram negative bacilli among clinical isolates at clinical laboratories of Army Medical College, Rawalpindi. J Ayub Med Coll Abbottabad. 2004;16:35–7.
Irfan S, Turton JF, Mehraj J, Siddiqui SZ, Haider S, Zafar A, Memon B, Afzal O, Hasan R. Molecular and epidemiological characterisation of clinical isolates of carbapenem-resistant Acinetobacter baumannii from public and private sector intensive care units in Karachi, Pakistan. J Hosp Infect. 2011;78:143–8.
Kaleem F, Usman J, Hassan A, Khan A. Frequency and susceptibility pattern of metallo-beta-lactamase producers in a hospital in Pakistan. J Infect Dev Ctries. 2010;4:810–3.
Ahmad S, Bacha N, Bakht J, Ahmed J. Characterization of pathogens involved in ventilator associated pneumonia in surgical and medical intensive care units—a single center experience. Pak J Pharm Sci. 2017;30(6):2091–9.
Qamar MU, Walsh TR, Toleman MA, Tyrrell JM, Saleem S, Aboklaish A, Jahan S. Dissemination of genetically diverse NDM-1, -5, -7 producing-Gram-negative pathogens isolated from pediatric patients in Pakistan. Future Microbiol. 2019;14:691–704.
Shah AA, Jamil B, Naseem S, Khan AW, Ali Y, Hussain K, Abbasi SA. Susceptibility pattern of tracheal tube isolates from intensive care unit of Fauji Foundation Hospital Rawalpindi. J Pak Med Assoc. 2019;69(7):981–4.
Evans BA, Hamouda A, Abbasi SA, Khan FA, Amyes SG. High prevalence of unrelated multidrug-resistant Acinetobacter baumannii isolates in Pakistani military hospitals. Int J Antimicrob Agents. 2011;37:580–1.
Diancourt L, Passet V, Nemec A, Dijkshoorn L, Brisse S. The population structure of Acinetobacter baumannii: expanding multiresistant clones from an ancestral susceptible genetic pool. PLoS ONE. 2010;5:e10034.
Karah N, Haldorsen B, Hermansen NO, Tveten Y, Ragnhildstveit E, Skutlaberg DH, Tofteland S, Sundsfjord A, Samuelsen O. Emergence of OXA carbapenemase- and 16S rRNA methylase-producing international clones of Acinetobacter baumannii in Norway. J Med Microbiol. 2011;60:515–21.
Karah N, Sundsfjord A, Towner K, Samuelsen O. Insights into the global molecular epidemiology of carbapenem non-susceptible clones of Acinetobacter baumannii. Drug Resist Updat. 2012;15:237–47.
Rodríguez C, Balderrama Yarhui N, Nastro M, Nunez Quezada T, Castro Canarte G, Magne Ventura R, Ugarte Cuba T, Valenzuela N, Roach F, Mota M, Burger N, Velázquez Aguayo G, Ortellado-Canese J, Bruni G, Pandolfo C, Bastyas N, Famiglietti A. Molecular epidemiology of carbapenem-resistant Acinetobacter baumannii in South America. J Med Microbiol. 2016;65(10):1088–91.
Turton JF, Woodford N, Glover J, Yarde S, Kaufmann ME, Pitt TL. Identification of Acinetobacter baumannii by detection of the blaOXA-51-like carbapenemase gene intrinsic to this species. J Clin Microbiol. 2006;44:2974–6.
La Scola B, Gundi VA, Khamis A, Raoult D. Sequencing of the rpoB gene and flanking spacers for molecular identification of Acinetobacter species. J Clin Microbiol. 2006;44:827–32.
Clinical and Laboratory Standards Institute (CLSI). Performance standards for antimicrobial susceptibility testing. 29th ed. Wayne: Clinical and Laboratory Standards Institute (CLSI); 2019.
Karah N, Jolley KA, Hall RM, Uhlin BE. Database for the ampC alleles in Acinetobacter baumannii. PLoS ONE. 2017;12(5):e0176695.
Pournaras S, Gogou V, Giannouli M, Dimitroulia E, Dafopoulou K, Tsakris A, Zarrilli R. Single locus sequence-based typing of blaOXA-51-like gene for rapid classification of Acinetobacter baumannii clinical isolates to international clones. J Clin Microbiol. 2014;52:1653–7.
Bertels F, Silander OK, Pachkov M, Rainey PB, van Nimwegen E. Automated reconstruction of whole-genome phylogenies from short-sequence reads. Mol Biol Evol. 2014;31(5):1077–88.
Bertini A, Poirel L, Mugnier PD, Villa L, Nordmann P, Carattoli A. Characterization and PCR-based replicon typing of resistance plasmids in Acinetobacter baumannii. Antimicrob Agents Chemother. 2010;54:4168–77.
Cameranesi MM, Limansky AS, Morán-Barrio J, Repizo GD, Viale AM. Three novel Acinetobacter baumannii plasmid replicase-homology groups inferred from the analysis of a multidrug-resistant clinical strain isolated in Argentina. J Infect Dis Epidemiol. 2017;3:46.
Hamidian M, Nigro SJ, Hall RM. Variants of the gentamicin and tobramycin resistance plasmid pRAY are widely distributed in Acinetobacter. J Antimicrob Chemother. 2012;67:2833–6.
Lean SS, Yeo CC. Small, enigmatic plasmids of the nosocomial pathogen, Acinetobacter baumannii: good, bad, who Knows? Front Microbiol. 2017;8:1547.
Larsen MV, Cosentino S, Rasmussen S, Friis C, Hasman H, Marvig RL, Jelsbak L, Sicheritz-Pontén T, Ussery DW, Aarestrup FM, Lund O. Multilocus sequence typing of total-genome-sequenced bacteria. J Clin Microbiol. 2012;50(4):1355–61.
Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, Aarestrup FM, Larsen MV. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother. 2012;67:2640–4.
Zhang Z, Schwartz S, Wagner L, Miller W. A greedy algorithm for aligning DNA sequences. J Comput Biol. 2000;7:203–14.
Siguier P, Perochon J, Lestrade L, Mahillon J, Chandler M. ISfinder: the reference centre for bacterial insertion sequences. Nucleic Acids Res. 2006;34:D32–6.
Hujer KM, Hujer AM, Hulten EA, Bajaksouzian S, Adams JM, Donskey CJ, Ecker DJ, Massire C, Eshoo MW, Sampath R, Thomson JM, Rather PN, Craft DW, Fishbain JT, Ewell AJ, Jacobs MR, Paterson DL, Bonomo RA. Analysis of antibiotic resistance genes in multidrug-resistant Acinetobacter sp. isolates from military and civilian patients treated at the Walter Reed Army Medical Center. Antimicrob Agents Chemother. 2006;50:4114–23.
Anwar M, Ejaz H, Zafar A, Hamid H. Phenotypic detection of metallo-β-lactamases in carbapenem resistant Acinetobacter baumannii isolated from pediatric patients in Pakistan. J Pathog. 2016;2016:1–6.
Begum S, Hasan F, Hussain S, Shah AA. Prevalence of multi drug resistant Acinetobacter baumannii in the clinical samples from Tertiary Care Hospital in Islamabad, Pakistan. Pak J Med Sci. 2013;29:1253–8.
Hasan B, Perveen K, Olsen B, Zahra R. Emergence of carbapenem resistant Acinetobacter baumannii in hospitals in Pakistan. J Med Microbiol. 2014;63:50–5.
Holt K, Kenyon JJ, Hamidian M, Schultz MB, Pickard DJ, Dougan G, Hall R. Five decades of genome evolution in the globally distributed, extensively antibiotic-resistant Acinetobacter baumannii global clone 1. Microb Genom. 2016;2(2):e000052.
Hall MR, Schwarz S. Resistance gene naming and numbering: is it a new gene or not? J Antimicrob Chemother. 2016;71(3):569–71.
Jacoby GA, Bonomo RA, Bradford PA, Bush K, Doi Y, Feldgarden M, Haft D, Klimke W, Nordmann P, Palzkill T, Poirel L, Prasad A, Rossolini GM, Walsh T. Comment on: resistance gene naming and numbering: is it a new gene or not? J Antimicrob Chemother. 2016;71(9):2677–8.
Sahl JW, Del Franco M, Pournaras S, Colman RE, Karah N, Dijkshoorn L, Zarrilli R. Phylogenetic and genomic diversity in isolates from the globally distributed Acinetobacter baumannii ST25 lineage. Sci Rep. 2015;5:15188.
Zhu L, Yan Z, Zhang Z, Zhou Q, Zhou J, Wakeland EK, Fang X, Xuan Z, Shen D, Li QZ. Complete genome analysis of three Acinetobacter baumannii clinical isolates in China for insight into the diversification of drug resistance elements. PLoS ONE. 2013;8(6):e66584.
Jia H, Chen Y, Wang J, Xie X, Ruan Z. Emerging challenges of whole-genome-sequencing-powered epidemiological surveillance of globally distributed clonal groups of bacterial infections, giving Acinetobacter baumannii ST195 as an example. Int J Med Microbiol. 2019;19:151339 (Epub ahead of print).
Li H, Liu F, Zhang Y, Wang X, Zhao C, Chen H, Zhang F, Zhu B, Hu Y, Wang H. Evolution of carbapenem-resistant Acinetobacter baumannii revealed through whole-genome sequencing and comparative genomic analysis. Antimicrob Agents Chemother. 2015;59(2):1168–76.
Eigenbrod T, Reuter S, Gross A, Kocer K, Günther F, Zimmermann S, Heeg K, Mutters NT, Nurjadi D. Molecular characterization of carbapenem-resistant Acinetobacter baumannii using WGS revealed missed transmission events in Germany from 2012–15. J Antimicrob Chemother. 2019;74(12):3473–80.
Towner KJ, Evans B, Villa L, Levi K, Hamouda A, Amyes SG, Carattoli A. Distribution of intrinsic plasmid replicase genes and their association with carbapenem-hydrolyzing class D β-lactamase genes in European clinical isolates of Acinetobacter baumannii. Antimicrob Agents Chemother. 2011;55:2154–9.
Chan AP, Sutton G, DePew J, Krishnakumar R, Choi Y, Huang XZ, Beck E, Harkins DM, Kim M, Lesho EP, Nikolich MP, Fouts DE. A novel method of consensus pan-chromosome assembly and large-scale comparative analysis reveal the highly flexible pan-genome of Acinetobacter baumannii. Genome Biol. 2015;16(1):143.
Wibberg D, Salto IP, Eikmeyer FG, Maus I, Winkler A, Nordmann P, Pühler A, Poirel L, Schlüter A. Complete genome sequencing of Acinetobacter baumannii strain K50 discloses the large conjugative plasmid pK50a encoding carbapenemase OXA-23 and extended-spectrum β-lactamase GES-11. Antimicrob Agents Chemother. 2018;62(5):e00212–e00218-18.
Rafei R, Pailhoriès H, Hamze M, Eveillard M, Mallat H, Dabboussi F, Joly-Guillou ML, Kempf M. Molecular epidemiology of Acinetobacter baumannii in different hospitals in Tripoli, Lebanon using blaOXA-51-like sequence based typing. BMC Microbiol. 2015;15:103.
Bonnin RA, Rotimi VO, Al Hubail M, Gasiorowski E, Al Sweih N, Nordmann P, Poirel L. Wide dissemination of GES-type carbapenemases in Acinetobacter baumannii isolates in Kuwait. Antimicrob Agents Chemother. 2013;57(1):183–8.
Castanheira M, Costello SE, Woosley LN, Deshpande LM, Davies TA, Jones RN. Evaluation of clonality and carbapenem resistance mechanisms among Acinetobacter baumannii-Acinetobacter calcoaceticus complex and Enterobacteriaceae isolates collected in European and Mediterranean countries and detection of two novel β-Lactamases, GES-22 and VIM-35. Antimicrob Agents Chemother. 2014;58(12):7358–66.
Hammerum AM, Hansen F, Skov MN, Stegger M, Andersen PS, Holm A, Jakobsen L, Justesen US. Investigation of a possible outbreak of carbapenem-resistant Acinetobacter baumannii in Odense, Denmark using PFGE, MLST and whole-genome-based SNPs. J Antimicrob Chemother. 2015;70(7):1965–8.
Vali L, Dashti K, Opazo-Capurro AF, Dashti AA, Al Obaid K, Evans BA. Diversity of multi-drug resistant Acinetobacter baumannii population in a major hospital in Kuwait. Front Microbiol. 2015;6:743.
Linz B, Mukhtar N, Shabbir MZ, Rivera I, Ivanov YV, Tahir Z, Yaqub T, Harvill ET. Virulent epidemic pneumonia in sheep caused by the human pathogen Acinetobacter baumannii. Front Microbiol. 2018;9:2616.
Turton FJ, Ward ME, Woodford N, Kaufmann ME, Pike R, Livermore DM, Pitt TL. The role of ISAba1 in expression of OXA carbapenemase genes in Acinetobacter baumannii. FEMS Microbiol Lett. 2006;258(1):72–7.
Evans BA, Amyes SG. OXA β-lactamases. Clin Microbiol Rev. 2014;27(2):241–63.
Li P, Niu W, Li H, Lei H, Liu W, Zhao X, Guo L, Zou D, Yuan X, Liu H, Yuan J, Bai C. Rapid detection of Acinetobacter baumannii and molecular epidemiology of carbapenem-resistant A. baumannii in two comprehensive hospitals of Beijing. China. Front Microbiol. 2015;6:997.
Vijayakumar S, Mathur P, Kapil A, Das BK, Ray P, Gautam V, Sistla S, Parija SC, Walia K, Ohri VC, Anandan S, Subramani K, Ramya I, Veeraraghavan B. Molecular characterization & epidemiology of carbapenem-resistant Acinetobacter baumannii collected across India. Indian J Med Res. 2019;149(2):240–6.
Hamidian M, Hall RM. ISAba1 targets a specific position upstream of the intrinsic ampC gene of Acinetobacter baumannii leading to cephalosporin resistance. J Antimicrob Chemother. 2013;68:2682–3.
Nigro SJ, Hall RM. Does the intrinsic oxaAb (blaOXA-51-like) gene of Acinetobacter baumannii confer resistance to carbapenems when activated by ISAba1? J Antimicrob Chemother. 2018;73(12):3518–20.
Cho YJ, Moon DC, Jin JS, Choi CH, Lee YC, Lee JC. Genetic basis of resistance to aminoglycosides in Acinetobacter spp. and spread of armA in Acinetobacter baumannii sequence group 1 in Korean hospitals. Diagn Microbiol Infect Dis. 2009;64:185–90.
Domingues S, Nielsen KM, da Silva GJ. The blaIMP-5-carrying integron in a clinical Acinetobacter baumannii strain is flanked by miniature inverted-repeat transposable elements (MITEs). J Antimicrob Chemother. 2011;66(11):2667–8.
Blackwell GA, Hall RM. The tet39 determinant and the msrE-mphE genes in Acinetobacter plasmids are each part of discrete modules flanked by inversely oriented pdif (XerC-XerD) sites. Antimicrob Agents Chemother. 2017;61(8):e00780-17.
Mindlin S, Petrenko A, Kurakov A, Beletsky A, Mardanov A, Petrova M. Resistance of permafrost and modern Acinetobacter lwoffii strains to heavy metals and arsenic revealed by genome analysis. Biomed Res Int. 2016;2016:3970831.
We thank all colleagues who generously provided isolates for this study. We acknowledge Dr. Raja Kamran Afzal (combined military hospital) for his assistance in collecting clinical isolates. We thank the team of curators of the Institute Pasteur Acinetobacter MLST system for curating the data and making them publicly available at http://pubmlst.org/abaumannii/.
Open access funding provided by Umea University. This work was supported by project grants from the Swedish Research Council (grants 2015-03007, 2015-06824, 2018-02914), Kempe Foundations (grants JCK-1527, JCK-1724), and HEC Pakistan (Grant 8666/Punjab/NRPU/2017). The work was performed as part of the Umeå Centre for Microbial Research (UCMR) Linnaeus Program at The Laboratory for Molecular Infection Medicine Sweden (MIMS) with support from Umeå University and the Swedish Research Council (grants 2007-8673 and 2016-06598). The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Ethics and approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1.
The geographic distance between the Combined Military Hospital in Lahore (CMH Lahore) and the Combined Military Hospital in Peshawar (CMH Peshawar) as drawn by Google Maps https://www.google.com/maps/dir/Combined+Military+Hospital,+Mall+Rd,+Peshawar+Cantonment,+Peshawar,+Khyber+Pakhtunkhwa,+Pakistan/CMH+Lahore,+Abdul+Rehman+Rd,+Saddar+Town,+Lahore,+Punjab,+Pakistanemail@example.com,72.6522947,8z/data=!4m14!4m13!1m5!1m1!1s0x38d917bbe7a07855:0xf8ccf377d1a61673!2m2!1d71.5423502!2d34.0033769!1m5!1m1!1s0x39190517756f887b:0xc49fc1d23f55d0b0!2m2!1d74.3726762!2d31.5408931!3e0.
Additional file 2: Table S1.
Primers used for amplification and sequencing of the blaOXA-51-like and ampC loci in 52 clinical isolates of Acinetobacter baumannii collected in Pakistan between 2013 and 2015.
Additional file 3: Table S2.
Antimicrobial susceptibility results for 52 clinical isolates of Acinetobacter baumannii collected in Pakistan between 2013 and 2015.
Additional file 4: File S1.
CLUSTAL O(1.2.4) multiple sequence alignment of ampC-3 and ampC-8.
Additional file 5: Table S3.
Genome assembly features of 25 clinical isolates of Acinetobacter baumannii collected in Pakistan between 2013 and 2015.
Additional file 6: Table S4.
Plasmid replicon genes in the whole genome sequences of 25 clinical isolates of Acinetobacter baumannii collected in Pakistan between 2013 and 2015.
Additional file 7: Table S5.
Acquired antimicrobial resistance genes retrived from the whole genome sequences of 25 clinical isolates of Acinetobacter baumannii collected in Pakistan between 2013 and 2015.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Karah, N., Khalid, F., Wai, S.N. et al. Molecular epidemiology and antimicrobial resistance features of Acinetobacter baumannii clinical isolates from Pakistan. Ann Clin Microbiol Antimicrob 19, 2 (2020). https://doi.org/10.1186/s12941-019-0344-7
- Strain typing
- International clone