An outbreak of multiply antibiotic-resistant ST49:ST128:KL11:OCL8 Acinetobacter baumannii isolates at a Sydney hospital

Genomic analysis of early ST32 Acinetobacter baumannii strains recovered in US military treatment facilities reveals distinct lineages and links to the origins of the Tn6168 ampC transposon

OBJECTIVES

To study the population structure and genomic characteristics, including antimicrobial resistance genes, plasmid types and surface polysaccharide type, of the globally distributed Acinetobacter baumannii belonging to ST32 (Institut Pasteur scheme).

METHODS

Antibiotic resistance phenotype for 19 antibiotics was determined using Vitek 2. Whole-genome sequencing was performed using the Illumina MiSeq platform. Genomes were assembled using Newbler. Phylogenetic analysis was done by determining the core-genome alignments using Panaroo v1.3, analysed in IQ-Tree2 v2.2.0.3 to construct Maximum Likelihood trees using the RaxML software. Resistance genes and IS were identified using the Abricate programme, and ISFinder databases.

RESULTS

One hundred and thirty-three (n = 133) ST32 A. baumannii isolates were analysed in this study. These genomes originated mainly from US military treatment facilities (n = 113), but also included additional publicly available genomes in GenBank (n = 20) recovered from a broad geographic distribution extending to Asia and South America. Phylogenetic analysis of all 133 genomes revealed at least four clades, with over 80 genomes forming a tightly clustered branch, suggesting they are likely to represent outbreak strains. Analysis of the ampC region showed that ST32 strains played a significant role in the formation of the widely distributed ampC transposon, Tn6168, and supplying DNA segments containing an ISAba1-ampC from ST32s via homologous recombination.

CONCLUSIONS

ST32 strains played a significant role in the evolution of antibiotic resistance in several widely distributed sequence types including ST1 (global clone 1) and ST3.

Rep_3 plasmids found in Acinetobacter species that include orfX also harbour dif modules

Plasmids found in Acinetobacter species are not found in other Gram-negative species. There are many distinct plasmid types and the majority encode a Rep_3 family RepA replication initiation protein. Among these, a number were known to carry dif modules. Here, the representative plasmid for each of the 78 reported R3 types was examined to identity features of the plasmid backbone that are associated with the carriage of dif modules. A conserved open reading frame designated orfX (IPR047783) was found downstream of repA in 35 of them, and the backbones of those 35 plasmids were bounded by recombination sites recognised by XerC and XerD, known as pdif sites. These plasmid backbones are all equivalent to a C-type dif module as the pdif sites are in the orientation D/C at one end and C/D at the other end, i.e. the XerC binding sites are internal. Hence, to provide the XerD binding sites and generate a complete plasmid at least one D-type dif module is needed. Phylogenies of the RepA and OrfX proteins revealed that plasmids with closely-related RepA proteins are not always associated with closely-related OrfX proteins and vice-versa indicating extensive backbone recombination. Folded structures of diverse OrfX proteins predicted using AlphaFold 3 revealed an N-terminal HTH domain followed by a long α-helix that is predicted to promote dimerization and a disordered C-terminus. Given the correlation between the presence of orfX and one or more dif modules, the possibility that OrfX is involved in dif module movement deserves to be investigated.

IS26 and the IS26 family: versatile resistance gene movers and genome reorganizers

SUMMARYIn Gram-negative bacteria, the insertion sequence IS26 is highly active in disseminating antibiotic resistance genes. IS26 can recruit a gene or group of genes into the mobile gene pool and support their continued dissemination to new locations by creating pseudo-compound transposons (PCTs) that can be further mobilized by the insertion sequence (IS). IS26 can also enhance expression of adjacent potential resistance genes. IS26 encodes a DDE transposase but has unique properties. It forms cointegrates between two separate DNA molecules using two mechanisms. The well-known copy-in (replicative) route generates an additional IS copy and duplicates the target site. The recently discovered and more efficient and targeted conservative mechanism requires an IS in both participating molecules and does not generate any new sequence. The unit of movement for PCTs, known as a translocatable unit or TU, includes only one IS26. TU formed by homologous recombination between the bounding IS26s can be reincorporated via either cointegration route. However, the targeted conservative reaction is key to generation of arrays of overlapping PCTs seen in resistant pathogens. Using the copy-in route, IS26 can also act on a site in the same DNA molecule, either inverting adjacent DNA or generating an adjacent deletion plus a circular molecule carrying the DNA segment lost and an IS copy. If reincorporated, these circular molecules create a new PCT. IS26 is the best characterized IS in the IS26 family, which includes IS257/IS431, ISSau10, IS1216, IS1006, and IS1008 that are also implicated in spreading resistance genes in Gram-positive and Gram-negative pathogens.

Variation in the plasmid backbone and dif module content of R3-T33 Acinetobacter plasmids

The predominant type of plasmids found in Acinetobacter species encode a Rep_3 initiation protein and many of these carry their accessory genes in dif modules. Here, available sequences of the 14 members of the group of Rep_3 plasmids typed as R3-T33, using a threshold of 95% identity in the repA gene, were compiled and compared. These plasmids were from various Acinetobacter species. The pdif sites were identified allowing the backbone and dif modules to be defined. As for other Rep_3 plasmids carrying dif modules, orfX encoding a protein of unknown function was found downstream of repA followed by a pdif site in the orientation XerC binding site-spacer-XerD binding site. Most backbones (n = 12) also included mobA and mobC genes but the two plasmids with the most diverged repA and orfX genes had different backbone contents. Although the gene content of the plasmid backbone was largely conserved, extensive recombinational exchange was detected and only two small groups carried identical or nearly identical backbones. Individual plasmids were associated with 1 to 13 dif modules. Many different dif modules were identified, including ones containing antibiotic or chromate resistance genes and several toxin/antitoxin gene pairs. In some cases, modules carrying the same genes were significantly diverged. Generally, the orientation of the pdif sites alternated such that C modules (XerC binding sites internal) alternated with D modules (XerD binding sites internal). However, fusions of two dif modules via mutational inactivation or loss of a pdif site were also detected.

The extensively antibiotic resistant ST111 Acinetobacter baumannii isolate RBH2 carries an extensive mobile element complement of plasmids, transposons and insertion sequences

The complete genome of RBH2, a sporadic, carbapenem resistant ST111 Acinetobacter baumannii isolate from Brisbane, Australia was determined and analysed. RBH2 is extensively resistant and the chromosome includes two transposons carrying antibiotic resistance genes, AbaR4 (oxa23 in Tn2006) and Tn7::Tn2006 (dfrA1, sat2, aadA1, oxa23). The chromosome also includes two copies of Tn6175, a transposon carrying putative copper resistance genes, and 1-17 copies of six different insertion sequences. RBH2 has six plasmids ranging in size from 6 kb - 141 kb, four carrying antibiotic resistance genes. Plasmids pRBH2-1 (aadB) and pRBH2-2 (aphA6 in TnaphA6) were found to be essentially identical to known plasmids pRAY*-v1 and pS21-1, respectively. The largest plasmids, pRBH2-5 (oxa23 in AbaR4) and pRBH2-6 (oxa23 in AbaR4::ISAba11 and sul2, tet(B), strA and strB in Tn6172) have known transfer-proficient relatives. pRBH2-5, an RP-T1 (RepAci6) plasmid, also carries a different putative copper resistance transposon related to Tn6177 found in pS21-2. The backbone of pRBH2-5 is related to those of previously described RepAci6 plasmids pAb-G7-2 and pA85-3 but has some distinctive features. Three different RepAci6 backbone types were distinguished, Type 1 (pAb-G7-2), Type 2 (pA85-3) and Type 3 (pRBH2-5 and pS21-2). pRBH2-6 is closely related to pAB3 and their backbones differ by only 5 SNPs. Plasmids pRBH2-3 and pRBH2-4 do not carry antibiotic resistance genes. pRBH2-3 does not include an identifiable rep gene and is a novel plasmid type. pRBH2-4 is of the R3-T3 type and includes segments of the larger pABTJ2 that heads this group. Other ST111 genomes carry different plasmids.

Acinetobacter baumannii GC2 Sublineage Carrying the aac(6')-Im Amikacin, Netilmicin, and Tobramycin Resistance Gene Cassette

The aminoglycoside antibiotics amikacin, gentamicin, and tobramycin are important therapeutic options for Acinetobacter iinfections. Several genes that confer resistance to one or more of these antibiotics are prevalent in the globally distributed resistant clones of Acinetobacter baumannii, but the aac(6')-Im (aacA16) gene (amikacin, netilmicin, and tobramycin resistance), first reported in isolates from South Korea, has rarely been reported since. In this study, GC2 isolates (1999 to 2002) from Brisbane, Australia, carrying aac(6')-Im and belonging to the ST2:ST423:KL6:OCL1 type were identified and sequenced. The aac(6')-Im gene and surrounds have been incorporated into one end of the IS26-bounded AbGRI2 antibiotic resistance island and are accompanied by a characteristic 70.3-kbp deletion of adjacent chromosome. The compete genome of the 1999 isolate F46 (RBH46) includes only two copies of ISAba1 (in AbGRI1-3 and upstream of ampC) but later isolates, which differ from one another by <10 single nucleotide differences (SND), carry two to seven additional shared copies. Several complete GC2 genomes with aac(6')-Im in an AbGRI2 island (2004 to 2017; several countries) found in GenBank and two additional Australian A. baumannii isolates (2006) carry different gene sets, KL2, KL9, KL40, or KL52, at the capsule locus. These genomes include ISAba1 copies in a different set of shared locations. The distribution of SND between F46 and AYP-A2, a 2013 ST2:ST208:KL2:OCL1 isolate from Victoria, Australia, revealed that a 640-kbp segment that includes KL2 and the AbGRI1 resistance island replaces the corresponding region in F46. Over 1,000 A. baumannii draft genomes also include aac(6')-Im, indicating that it is currently globally disseminated and significantly underreported. IMPORTANCE Aminoglycosides are important therapeutic options for treatment of Acinetobacter infections. Here, we show that a little-known aminoglycoside resistance gene, aac(6')-Im (aacA16), that confers amikacin, netilmicin, and tobramycin resistance has been circulating undetected for many years in a sublineage of A. baumannii global clone 2 (GC2), generally with a second aminoglycoside resistance gene, aacC1, which confers resistance to gentamicin. These two genes are commonly found together in GC2 complete and draft genomes and globally distributed. One isolate appears to be ancestral, as its genome contains few ISAba1 copies, providing insight into the original source of this insertion sequence (IS), which is abundant in most GC2 isolates. Tracking ISAba1 spread can provide a simple means to track the development and ongoing evolution as well as the dissemination of specific lineages and detect the formation of many sublineages. The complete ancestral genome will provide an essential base point for tracking this process.

Insertion sequences related to ISAjo2 target pdif and dif sites and belong to a new IS family, the IS1202 family

Several insertion sequences (IS) found in various Acinetobacter species exhibit target specificity. They are found, in the same orientation, 5 bp from the XerC binding site of the pdif sites associated with dif modules in Acinetobacter plasmids, and searches revealed they are also found near chromosomal dif sites of Acinetobacter species. These IS are 1.5 kb long, bounded by 24-26 bp imperfect terminal inverted repeats (TIRs) and encode a large transposase of 441-457 aa. They generate 5 bp target site duplications (TSDs). Structural predictions of the ISAjo2 transposase, TnpAjo2, modelled on TnsB of Tn7 revealed two N-terminal HTH domains followed by an RNaseH fold (DDE domain), a β barrel and a C-terminal domain. Similar to Tn7, the outer IS ends are 5'-TGT and ACA-3', and an additional Tnp binding site, corresponding to the internal portion of the IR, is found near each end. However, the Acinetobacter IS do not encode further proteins related to those required by Tn7 for targeted transposition, and the transposase may interact directly with XerC bound to a dif-like site. We propose that these IS, currently in the IS1202 group in the not characterized yet (NCY) category in ISFinder, are part of a distinct IS1202 family. Other IS listed as in the IS1202 group encode transposases related to TnpAjo2 (25-56 % amino acid identity) and have similar TIRs but fall into three groups based on the TSD length (3-5, >15, 0 bp). Those with 3-5 bp TSDs may also target dif-like sites but targets were not found for the other groups.

Dynamic state of plasmid genomic architectures resulting from XerC/D-mediated site-specific recombination in Acinetobacter baumannii Rep_3 superfamily resistance plasmids carrying blaOXA-58 - and TnaphA6-resistance modules

The acquisition of bla _OXA genes encoding different carbapenem-hydrolyzing class-D β-lactamases (CHDL) represents a main determinant of carbapenem resistance in the nosocomial pathogen Acinetobacter baumannii. The bla_OXA-58 gene, in particular, is generally embedded in similar resistance modules (RM) carried by plasmids unique to the Acinetobacter genus lacking self-transferability. The ample variations in the immediate genomic contexts in which bla_OXA-58 -containing RMs are inserted among these plasmids, and the almost invariable presence at their borders of non-identical 28-bp sequences potentially recognized by the host XerC and XerD tyrosine recombinases (pXerC/D-like sites), suggested an involvement of these sites in the lateral mobilization of the gene structures they encircle. However, whether and how these pXerC/D sites participate in this process is only beginning to be understood. Here, we used a series of experimental approaches to analyze the contribution of pXerC/D-mediated site-specific recombination to the generation of structural diversity between resistance plasmids carrying pXerC/D-bounded bla _OXA-58- and TnaphA6-containing RM harbored by two phylogenetically- and epidemiologically-closely related A. baumannii strains of our collection, Ab242 and Ab825, during adaptation to the hospital environment. Our analysis disclosed the existence of different bona fide pairs of recombinationally-active pXerC/D sites in these plasmids, some mediating reversible intramolecular inversions and others reversible plasmid fusions/resolutions. All of the identified recombinationally-active pairs shared identical GGTGTA sequences at the cr spacer separating the XerC- and XerD-binding regions. The fusion of two Ab825 plasmids mediated by a pair of recombinationally-active pXerC/D sites displaying sequence differences at the cr spacer could be inferred on the basis of sequence comparison analysis, but no evidence of reversibility could be obtained in this case. The reversible plasmid genome rearrangements mediated by recombinationally-active pairs of pXerC/D sites reported here probably represents an ancient mechanism of generating structural diversity in the Acinetobacter plasmid pool. This recursive process could facilitate a rapid adaptation of an eventual bacterial host to changing environments, and has certainly contributed to the evolution of Acinetobacter plasmids and the capture and dissemination of bla _OXA-58 genes among Acinetobacter and non-Acinetobacter populations co-residing in the hospital niche.

Pdif-mediated antibiotic resistance genes transfer in bacteria identified by pdifFinder

Modules consisting of antibiotic resistance genes (ARGs) flanked by inverted repeat Xer-specific recombination sites were thought to be mobile genetic elements that promote horizontal transmission. Less frequently, the presence of mobile modules in plasmids, which facilitate a pdif-mediated ARGs transfer, has been reported. Here, numerous ARGs and toxin-antitoxin genes have been found in pdif site pairs. However, the mechanisms underlying this apparent genetic mobility is currently not understood, and the studies relating to pdif-mediated ARGs transfer onto most bacterial genera are lacking. We developed the web server pdifFinder based on an algorithm called PdifSM that allows the prediction of diverse pdif-ARGs modules in bacterial genomes. Using test set consisting of almost 32 thousand plasmids from 717 species, PdifSM identified 481 plasmids from various bacteria containing pdif sites with ARGs. We found 28-bp-long elements from different genera with clear base preferences. The data we obtained indicate that XerCD-dif site-specific recombination mechanism may have evolutionary adapted to facilitate the pdif-mediated ARGs transfer. Through multiple sequence alignment and evolutionary analyses of duplicated pdif-ARGs modules, we discovered that pdif sites allow an interspecies transfer of ARGs but also across different genera. Mutations in pdif sites generate diverse arrays of modules which mediate multidrug-resistance, as these contain variable numbers of diverse ARGs, insertion sequences and other functional genes. The identification of pdif-ARGs modules and studies focused on the mechanism of ARGs co-transfer will help us to understand and possibly allow controlling the spread of MDR bacteria in clinical settings. The pdifFinder code, standalone software package and description with tutorials are available at https://github.com/mjshao06/pdifFinder.

Plasmids as Key Players in Acinetobacter Adaptation

This review briefly summarizes the data on the mechanisms of development of the adaptability of Acinetobacters to various living conditions in the environment and in the clinic. A comparative analysis of the genomes of free-living and clinical strains of A. lwoffii, as well as the genomes of A. lwoffii and A. baumannii, has been carried out. It has been shown that plasmids, both large and small, play a key role in the formation of the adaptability of Acinetobacter to their living conditions. In particular, it has been demonstrated that the plasmids of various strains of Acinetobacter differ from each other in their structure and gene composition depending on the lifestyle of their host bacteria. Plasmids of modern strains are enriched with antibiotic-resistant genes, while the content of genes involved in resistance to heavy metals and arsenic is comparable to plasmids from modern and ancient strains. It is concluded that Acinetobacter plasmids may ensure the survival of host bacteria under conditions of various types of environmental and clinical stresses. A brief overview of the main mechanisms of horizontal gene transfer on plasmids inherent in Acinetobacter strains is also given.

GR13-type plasmids in Acinetobacter potentiate the accumulation and horizontal transfer of diverse accessory genes

Carbapenem and other antibiotic resistance genes (ARGs) can be found in plasmids in Acinetobacter, but many plasmid types in this genus have not been well-characterized. Here we describe the distribution, diversity and evolutionary capacity of rep group 13 (GR13) plasmids that are found in Acinetobacter species from diverse environments. Our investigation was prompted by the discovery of two GR13 plasmids in A. baumannii isolated in an intensive care unit (ICU). The plasmids harbour distinct accessory genes: pDETAB5 contains bla_NDM-1 and genes that confer resistance to four further antibiotic classes, while pDETAB13 carries putative alcohol tolerance determinants. Both plasmids contain multiple dif modules, which are flanked by pdif sites recognized by XerC/XerD tyrosine recombinases. The ARG-containing dif modules in pDETAB5 are almost identical to those found in pDETAB2, a GR34 plasmid from an unrelated A. baumannii isolated in the same ICU a month prior. Examination of a further 41 complete, publicly available plasmid sequences revealed that the GR13 pangenome consists of just four core but 1186 accessory genes, 123 in the shell and 1063 in the cloud, reflecting substantial capacity for diversification. The GR13 core genome includes genes for replication and partitioning, and for a putative tyrosine recombinase. Accessory segments encode proteins with diverse putative functions, including for metabolism, antibiotic/heavy metal/alcohol tolerance, restriction-modification, an anti-phage system and multiple toxin–antitoxin systems. The movement of dif modules and actions of insertion sequences play an important role in generating diversity in GR13 plasmids. Discrete GR13 plasmid lineages are internationally disseminated and found in multiple Acinetobacter species, which suggests they are important platforms for the accumulation, horizontal transmission and persistence of accessory genes in this genus.

Phylogenomics of two ST1 antibiotic-susceptible non-clinical Acinetobacter baumannii strains reveals multiple lineages and complex evolutionary history in global clone 1

Acinetobacter baumannii is an opportunistic pathogen that is difficult to treat due to its resistance to extreme conditions, including desiccation and antibiotics. Most strains causing outbreaks around the world belong to two main global lineages, namely global clones 1 and 2 (GC1 and GC2). Here, we used a combination of Illumina short read and MinION (Oxford Nanopore) long-read sequence data with a hybrid assembly approach to complete the genome sequence of two antibiotic-sensitive GC1 strains, Ex003 and Ax270, recovered in Lebanon from water and a rectal swab of a cat, respectively. Phylogenetic analysis of Ax270 and Ex003 with 186 publicly available GC1 genomes revealed two major clades, including five main lineages (L1–L5), and four single-isolate lineages outside of the two clades. Ax270 and Ex003, along with AB307-0294 and MRSN7213 (both predicted antibiotic-susceptible isolates) represent these individual lineages. Antibiotic resistance islands and transposons interrupting the comM gene remain important features in L1–L5, with L1 associated with the AbaR-type resistance islands, L2 with AbaR4, L3 strains containing either AbaR4 or its variants as well as Tn6022::ISAba42, and L4 and L5 associated with Tn6022 or its variants. Analysis of the capsule (KL) and outer core (OCL) polysaccharide loci further revealed a complex evolutionary history probably involving many recombination events. As more genomes become available, more GC1 lineages continue to emerge. However, genome sequence data from more diverse geographical regions are needed to draw a more accurate population structure of this globally distributed clone.

Two carbapenem-resistant ST1:ST231:KL1:OCL1 Acinetobacter baumannii strains recovered in Tehran, Iran, carry AbaR31 in the chromosome and AbaR4 and TnaphA6 in a RepAci6 plasmid

Objectives

To analyse the context of genes conferring antibiotic resistance in two carbapenem-resistant Acinetobacter baumannii isolates recovered in Tehran, Iran.

Methods

The antibiotic resistance phenotype for 28 antibiotics was determined using disc diffusion. The whole genome sequences of ABH008 and ABS200 were determined using the Illumina HiSeq X Ten platform. Resistance genes were identified using ResFinder and multilocus sequence types were determined using the Oxford and Institut Pasteur schemes.

Results

Isolates ABH008 and ABS200, recovered in 2012 and 2013, respectively, in two different Tehran hospitals, belong to the common global clone 1 lineage, ST1_IP and ST231_OX. They are resistant to sulfamethoxazole, tetracycline, gentamicin, amikacin, third-generation cephalosporins and carbapenems. Despite being isolated in different hospitals, phylogenetic analysis indicated they are closely related. Consistent with this, both isolates carry catA1, sul1, aacC1 and aadA1 in a novel variant of the AbaR3-type resistance island, named AbaR31. Both isolates are resistant to amikacin and carbapenems owing to aphA6 and oxa23, respectively. The oxa23 gene is located in the AbaR4 resistance island, and aphA6 in TnaphA6, and both mobile elements are in an ∼90 kbp plasmid encoding the putative RepAci6 replication initiation protein. Resistance to third-generation cephalosporins is due to the acquisition by homologous recombination of a 5 kb DNA segment that contains ISAba1-ampC from a ST623 strain.

Conclusions

The resistance gene complements of ABH008 and ABS200 were found in AbaR31 and a plasmid that encodes RepAci6. The close genetic relationship of ABH008 and ABS200, despite each being recovered from different hospitals, indicates transmission between the two hospitals.

Recovery of small plasmid sequences via Oxford Nanopore sequencing

Oxford Nanopore Technologies (ONT) sequencing platforms currently offer two approaches to whole-genome native-DNA library preparation: ligation and rapid. In this study, we compared these two approaches for bacterial whole-genome sequencing, with a specific aim of assessing their ability to recover small plasmid sequences. To do so, we sequenced DNA from seven plasmid-rich bacterial isolates in three different ways: ONT ligation, ONT rapid and Illumina. Using the Illumina read depths to approximate true plasmid abundance, we found that small plasmids (<20 kbp) were underrepresented in ONT ligation read sets (by a mean factor of ~4) but were not underrepresented in ONT rapid read sets. This effect correlated with plasmid size, with the smallest plasmids being the most underrepresented in ONT ligation read sets. We also found lower rates of chimaeric reads in the rapid read sets relative to ligation read sets. These results show that when small plasmid recovery is important, ONT rapid library preparations are preferable to ligation-based protocols.