Phenotypic and Molecular Characteristics of Carbapenem-Resistant Acinetobacter baumannii Isolates from Bulgarian Intensive Care Unit Patients

Carbapenem-resistant Acinetobacter baumannii (CRAB) is designated as an urgent public health threat, both due to its remarkable multidrug resistance and propensity for clonal spread. This study aimed to explore the phenotypic and molecular characteristics of antimicrobial resistance in CRAB isolates (n = 73) from intensive care unit (ICU) patients in two university hospitals in Bulgaria (2018–2019). The methodology included antimicrobial susceptibility testing, PCR, whole-genome sequencing (WGS), and phylogenomic analysis. The resistance rates were as follows: imipenem, 100%; meropenem, 100%; amikacin, 98.6%; gentamicin, 89%; tobramycin, 86.3%; levofloxacin, 100%; trimethoprim–sulfamethoxazole, 75.3%; tigecycline, 86.3%; colistin, 0%; and ampicillin–sulbactam, 13.7%. All isolates harbored blaOXA-51-like genes. The frequencies of distribution of other antimicrobial resistance genes (ARGs) were: blaOXA-23-like, 98.6%; blaOXA-24/40-like, 2.7%; armA, 86.3%; and sul1, 75.3%. The WGS of selected extensively drug-resistant A. baumannii (XDR-AB) isolates (n = 3) revealed the presence of OXA-23 and OXA-66 carbapenem-hydrolyzing class D β-lactamases in all isolates, and OXA-72 carbapenemase in one of them. Various insertion sequencies, such as ISAba24, ISAba31, ISAba125, ISVsa3, IS17, and IS6100, were also detected, providing increased ability for horizontal transfer of ARGs. The isolates belonged to the widespread high-risk sequence types ST2 (n = 2) and ST636 (n = 1) (Pasteur scheme). Our results show the presence of XDR-AB isolates, carrying a variety of ARGs, in Bulgarian ICU settings, which highlights the crucial need for nationwide surveillance, especially in the conditions of extensive antibiotic usage during COVID-19.


Introduction
Acinetobacter baumannii belongs to a group of antimicrobial-resistant nosocomial pathogens, designated by the acronym "ESKAPE" (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter spp.). These bacteria represent a substantial therapeutic problem worldwide [1]. They frequently "escape" from the most commonly used antimicrobial treatment via acquisition and/or development of various resistance mechanisms [2].
A. baumannii infections are usually hard to treat owing to innate and acquired resistance to multiple antimicrobial agents. Carbapenems, recognized as the therapy of choice for serious HAIs and for the care of patients affected by MDR-AB, are no longer effective in some cases [15,16]. Carbapenem-resistant A. baumannii (CRAB) is designated as an urgent public health threat by the Centers for Disease Control and Prevention [17]. Its increasing importance has been recognized by the World Health Organization (WHO), which classified CRAB amongst the critical-priority pathogens in the 'Global priority list of antibiotic-resistant bacteria to guide research, discovery, and development of new antibiotics, 2017' [18]. Recent trends reveal that many infections are caused by CRAB, MDR-AB, extensively drug-resistant (XDR-AB) or even pan drug-resistant A. baumannii (PDR-AB), for which effective treatment options are severely limited [16,19,20].
Carbapenem resistance in A. baumannii is predominantly related to the production of carbapenem-hydrolyzing β-lactamases (carbapenemases), and rarely to the loss or modification of outer membrane proteins and/or the overexpression of active efflux pump systems [16,21]. In addition to the naturally occurring class D β-lactamases, also called oxacillinases (OXAs), several acquired carbapenem-hydrolyzing class D β-lactamases (CHDLs) have been identified in A. baumannii, among which OXA-23-like enzymes are the most common [22]. Other commonly found CHDLs include OXA-24/40-like and OXA-58-like enzymes [23]. It is worth mentioning that most acquired CHDLs possess weak carbapenemase activity, but the presence of specific insertion sequences (ISs) upstream of the CHDL gene leads to their overexpression, conferring carbapenem resistance [24]. A variety of class B metallo-β-lactamases (MBLs) have been identified in A. baumannii, but these occur at a much lower frequency than CHDLs. Class A serine β-lactamases with carbapenemase activity have been reported occasionally [3,25]. Multiple reports have demonstrated higher resistance percentages of CRAB in Southern and Eastern Europe than Northern Europe, as well as very high carbapenem resistance rates in the Balkan countries [26].
In addition to its remarkable multidrug resistance, A. baumannii is characterized by a propensity for clonal spread. During the past five decades, the evolution of A. baumannii was mainly driven by two global clones or international clones (ICs) of high risk (IC1 and IC2) [27]. Molecular epidemiological studies have assigned nine ICs (IC1-9) [28,29], the most widespread of which is IC2, commonly harboring the acquired OXA-23 carbapenemase [28,30,31].
The aim of this study was to explore the phenotypic and molecular characteristics of antimicrobial resistance in CRAB isolates from ICU patients treated at two large hospitals in Bulgaria (2018-2019). A phylogenomic analysis was also performed.
A. baumannii ATCC 19,606 was used as a control strain for species identification, antimicrobial susceptibility testing and quinolone-resistance-determining regions (QRDR) analysis.
All procedures involving patients were performed in accordance with the ethical standards of the Medical University of Sofia, Bulgaria, and the Helsinki Declaration of 1964 and its later amendments. The current study was focused solely on bacterial isolates and no personal patient information was used; therefore, formal consent was not required.

Species Identification of the Isolates
Species identification was done using the BD Phoenix M50 automated system (BD, Franklin Lakes, NJ, USA). Confirmation of species identity was performed by a polymerase chain reaction (PCR)-based gyrB method, as described previously [32]. The method distinguished between the closely related species A. baumannii and Acinetobacter nosocomialis. The identification of three selected isolates was further confirmed by analyzing the assembled draft genome sequence using the Microbial Genomes Atlas (MiGA) Web server [33]. The included workflow for the National Center for Biotechnology Information (NCBI) Genome Database, prokaryotic section was followed with default settings.

DNA Isolation
Total DNA from all investigated strains was isolated by the DNeasy Blood and Tissue Kit (QIAGEN, Hilden, Germany), according to the manufacturer's instructions, from 3 mL of overnight cultures inoculated with a single colony.
DNA was amplified in a Gene ProThermal Cycler (Bioer Technology, Hangzhou, China) using the following protocol: initial denaturation at 95 • C for 5 min, followed by 30 cycles of denaturation at 95 • C for 45 s, annealing at 50-60 • C for 45 s and extension at 72 • C for 45 s, and a single final extension at 72 • C for 7 min. PCR products were separated in 1.5% agarose gel for 50 min at 130 V, stained with SimplySafe (0.05 µL/mL) (EURx, Gdansk, Poland) and detected by ultraviolet light (wavelength 312 nm). Amplified genes were identified on the basis of fragment size (Table 1).

Whole-Genome Sequencing (WGS)
Three XDR-AB isolates (Aba52, Aba176 and Aba190) were subjected to WGS in order to perform a resistome analysis. The inclusion criteria were: different place of isolation (Sofia and Varna), different year of isolation and origin (regarding the isolates from Sofia), different level of aminoglycoside resistance, different ARGs previously determined by PCR. Furthermore, these were the only three isolates in which ISAba1 was not identified using PCR. The WGS was carried out using DNA nanoball sequencing technology, as previously described [45]. Briefly, genomic DNA obtained from all selected strains was randomly fragmented using a Covaris g-TUBE device, and fragments were size-selected by magnetic beads to an average size of 200 to 400 bp. The purified fragments from each sample were end-repaired, 3 -adenylated, ligated to adapters, and, then, PCR-amplified. All libraries generated in this way were then uploaded onto an MGISEQ-2000 platform (BGI Group, Hong Kong, China). The following sequencing step was done generating 2 × 150 bp paired-end reads.

Draft Genome Assembly
All steps of quality control, raw reads preprocessing and draft genome assembly were carried out through the Galaxy online platform [46,47]. Default parameters were used for all following software tools, unless otherwise specified. The entire procedure was performed using the following tools: FastQC v0.11.9, Trimmomatic v0.38 and SPAdes v3.12.0.

Resistome and Multilocus Sequence Typing (MLST) Analyses
The draft genome contigs were screened for ARGs using the ABRicate tool (Galaxy Version 1.0.1) with the following settings: NCBI Bacterial Antimicrobial Resistance Reference Gene Database, Minimum DNA identity (70%) and Minimum DNA coverage (60%). The read-mapping against the pA105-2 plasmid was performed applying the Bowtie2 tool (Galaxy version 2.5.0) and then visualized by the BLAST Ring Image Generator (BRIG) program [48]. All mobile genetic elements were detected using the MobileElementFinder tool with default settings [49]. Afterwards, MLST was done on the assembled draft genome sequences using the MLST tool (Galaxy Version 2.19.0, https://usegalaxy.eu/, accessed on 18 January 2023). Classification of the isolates into sequence types (ST) was done according to the Pasteur scheme. This scheme includes the identification of seven internal housekeeping genes: cpn60 (60-kDa chaperonin), fusA (elongation factor EF-G), gltA (citrate synthase), pyrG (CTP synthase), recA (homologous recombination factor), rplB (50S ribosomal protein L2) and rpoB (RNA polymerase subunit B) [50].

Phylogenomic Analysis
A set of 66 genomes of colistin-susceptible A. baumannii isolates (Table S1) were uploaded to the KBase platform [51], together with our WGS-analyzed isolates (Aba52, Aba176 and Aba190). Then, all sequences were annotated with Prokka (v1.14.5) and SpeciesTree (v2.2.0) with default settings, used to build a phylogenetic tree of all isolates based on a collection of 49 core universal genes defined by COG (cluster of orthologous groups) gene families.

Statistical Analysis
The antimicrobial resistance rates of the studied A. baumannii isolates were compared with recent results reported by other authors, using the Student's t-test. For simple comparison tests, a p-value of <0.05 was considered statistically significant. To counteract the problem of multiple comparisons, when used, a Bonferroni correction was applied.
A total of 63 of the 73 studied CRAB isolates (86.3%) revealed high-level resistance to AMK, GEN and TOB (MICs > 256 mg/L), and were later identified as ArmA and OXA-23 co-producers. SXT-resistant A. baumannii isolates (75.3%) were shown to be sul1-positive.

Draft Genome Assemblies: Evaluation and Comparison
Three XDR-AB isolates (Aba52, Aba176 and Aba190) were subjected to WGS-based resistome analysis. The phenotypic and genotypic characteristics of the selected isolates are summarized in Table 2.  The three assembled draft genomes varied in size between 3.96 and 4.14 Mbp, and their GC content was 38.9% (Table 3).

Resistome Analysis
All identified ARGs found in the three sequenced genomes are listed in Table 4. Moreover, a wide variety of ISs were identified in all analyzed isolates, while only the Aba190 strain harbored a single Tn6018 transposon.
The 5 ends of the bla ADC-25 (encoding class C extended-spectrum β-lactamase ADC-25) and bla OXA-23 β-lactam resistance determinants were always found next to a contig end, suggesting the presence of ISs at the start of the corresponding promotor regions.
Isolate Aba190 was found to harbor a bla OXA-72 gene. The sequencing coverage of its corresponding contig was more than six times higher compared to the contig including the bla OXA-66 gene with a known chromosomal location. Such coverage bias suggested that bla OXA-72 was located on a plasmid. This hypothesis was also supported by the fact that just next to it, two genes encoding for proteins involved in the replication of DNA were found. Blast search with the 3785 bp long contig generated a perfect hit (100% query coverage and identity) against the bla OXA-72 -carrying plasmid pA105-2 (KR535993.1), isolated from the A. baumannii A105 strain. Mapping the sequencing reads against its sequence showed that Aba190 also possessed the remaining plasmid elements, as illustrated in Figure 1.   A variety of structural genes encoding aminoglycoside-modifying enzymes were identified in the WGS-subjected XDR-AB isolates. Moreover, armA genes were found in the isolates from Sofia (Aba52 and Aba176), which is in accordance with their high-level aminoglycoside resistance.
All isolates investigated were positive for sul1, qacEΔ1 and intI1 genes, suggesting the presence of class 1 integrons. The complete resistance gene cassettes embedded into the variable regions were only partially recovered due to the limitations of the used sequencing technology with a short reading frame. They contained aac(3)-Ia and ant(3")-Ia A variety of structural genes encoding aminoglycoside-modifying enzymes were identified in the WGS-subjected XDR-AB isolates. Moreover, armA genes were found in the isolates from Sofia (Aba52 and Aba176), which is in accordance with their high-level aminoglycoside resistance.
All isolates investigated were positive for sul1, qacE∆1 and intI1 genes, suggesting the presence of class 1 integrons. The complete resistance gene cassettes embedded into the variable regions were only partially recovered due to the limitations of the used sequencing technology with a short reading frame. They contained aac(3)-Ia and ant(3 )-Ia genes, encoding for aminoglycoside-modifying enzymes.
Comparison of the QRDRs of the three isolates with the corresponding region of A. baumannii ATCC 19,606 revealed that all of them possessed identical amino acid substitutions (pS83L in gyrA and pS80I in parC), with the previously reported role in quinolone resistance [52].

Phylogenomic Analysis
To analyze the local epidemiology of our isolates, a phylogenetic tree based on SNPs in 49 universal genes was built, as described above. The 66 genomes, included in the analysis were previously described [53]. They all originated from colistin-susceptible isolates. The resulting tree is presented in Figure 2.

Phylogenomic Analysis
To analyze the local epidemiology of our isolates, a phylogenetic tree based on SNPs in 49 universal genes was built, as described above. The 66 genomes, included in the analysis were previously described [53]. They all originated from colistin-susceptible isolates. The resulting tree is presented in Figure 2.

Discussion
The investigated clinical A. baumannii isolates were problematic ICU pathogens. Antimicrobial susceptibility testing revealed retained activity of only two of the antibiotics studied: COL (0% resistance) and SAM (13.7% resistance). These antibiotics have been involved in the recommended therapeutic regimens of severe infections caused by MDR-AB at the respective hospitals [54]. In line with global trends, to treat severe CRAB infections in critically ill patients, such as hospital-acquired pneumonia, VAP, BSIs with severe sepsis/septic shock, combined therapy of two agents active in vitro was increasingly being used, including COL, high-dose SAM, an aminoglycoside, high-dose TGC, and high-dose minocycline [15]. Treatment recommendations for PDR-AB infections included cefiderocol if available, or triple-drug combination of COL, high-dose MEM and high-dose SAM [15,55,56]. A recent study demonstrated the clinical efficacy of rescue therapy with cefiderocol for severe BSIs and VAP caused by CRAB among Italian critically ill patients [57]. Additionally, over the past five years, there has been an increasing number of studies reporting positive results in the treatment of CRAB/MDR-AB infections with the durlobactam-sulbactam combination [58]. Durlobactam is a new member of the diazabicyclooctane class of β-lactamase inhibitors, with broad-spectrum activity against Ambler class A, C, and D β-lactamases. Sulbactam is a first generation β-lactamase inhibitor, which also possesses antibacterial activity against A. baumannii, due to its selective binding to the penicillin-binding proteins 1, 2, and 3 [59]. To date, the frequency of A. baumannii spontaneous resistance to this dual β-lactamase inhibitor combination has been low [60]. It is worth mentioning that the global prevalence of MDR-AB and XDR-AB isolates has renewed interest in the therapy with non-antibiotic options, such as lytic bacteriophages [61,62] and bacteriophage-encoded endolysin [63].
It is estimated that resistance to carbapenems alone is sufficient for A. baumannii to be considered a highly resistant and therapeutically problematic pathogen [64]. In the current study, the prevalence of XDR-AB isolates (78.1%) was significantly higher (p < 0.001) than that observed among nosocomial CRAB isolated at four Bulgarian university hospitals during the period of 2014-2016 (12.4%), as well as than their proportion (21.7%) in CRAB isolates obtained from critically ill patients undergoing renal replacement therapy at a large university hospital in Sofia, Bulgaria (2016-2018) [65,66]. Analysis of the antimicrobial resistance of the investigated CRAB isolates compared to that among the isolates from the aforementioned studies showed a chronological increase in resistance to both TOB (86.3% resistant isolates vs. 55.8% (2014-2016) and 60.9% (2016-2018); p < 0.001 and p < 0.05, respectively) and TGC (86.3% vs. 58.4% non-susceptible isolates (2014-2016), p < 0.001, and 69.6% resistant isolates (2016-2018)). Moreover, it is worth noting a significant decrease in SAM resistance reported in the present study (13.7%) compared to the two earlier Bulgarian studies (41.6-43.5%; p < 0.01-0.001) [65,66], which confirms the role of this antimicrobial agent in treatment regimens for ICU patients [15]. Our finding contrasted with the large proportion of contemporary CRAB isolates demonstrating high-level resistance to sulbactam, mainly associated with the spread of OXA-23 producers in high-prevalence geographic areas, such as China [67].
In the present study, the frequency of genes encoding intrinsic OXAs (OXA-51-like) was absolute. Previous studies have shown that genes for naturally occurring OXA-51like enzymes of A. baumannii are localized on the chromosome and their expression is regulated by ISAba1 [23], found in almost all our isolates. Insertion of ISAba1 at the start of the promotor region, immediately upstream the structural gene, results in carbapenemase activity of these enzymes [24,71]. In addition, all CRAB isolates harbored genes for acquired CHDLs as follows: 71 were identified to be positive for bla OXA-23-like ; one isolate was bla OXA -  [72]. The bla OXA-23 gene, associated with the upstream-located ISAba1, was identified as the mechanism responsible for carbapenem resistance in all isolates. Thereafter, widespread dissemination of clonally-related OXA-23-producing CRAB isolates was also found at large hospitals in other cities of the country [65,[73][74][75][76]. bla OXA-23 is the predominant carbapenemase-coding gene, being present in two-thirds of the CRAB isolates from countries of all continents [26,71]. Further, this gene has been reported to be the most prevalent antimicrobial resistance determinants among CRAB isolates in several studies in the Balkan region, including Croatia [68], Greece [69], Romania [77], and Turkey [78], as well as in the EURECA study (67.7%) conducted on invasive CRAB isolates from 10 European countries, most of them, again, from the Balkans [26]. The bla OXA-23 gene is either located on the chromosome or on plasmids and is associated with different genetic structures [22,79].
In the present study, the frequency of distribution of bla OXA-24/40-like genes found only in SM-V isolates was significantly lower than that detected in earlies CRAB isolates (2014-2016) from four Bulgarian large hospitals (2.7% vs. 66.4%; p < 0.001) [65]. In accordance with our finding, in the cited study, OXA-24/40-like CHDLs were reported as the most common mechanism, conferring carbapenem resistance among CRAB isolates in SM-V. The isolates were clonally related and approximately 86.4% of them were identified as bla OXA-24/40-like -positive [65]. The first Bulgarian bla OXA-24-like -carrynig CRAB isolates were identified in 2012 [80]. Through sequencing techniques, OXA-72 carbapenemase, a variant of OXA-24/40-like CHDLs, was later found in the country [65,76]. One of our two bla OXA-24/40-like -postive isolates (Aba190) was also determined as an OXA-72 producer by WGS. The OXA-72 carbapenemase was reported in the Balkan region for the first time in Serbia in 2016 [81]. Moreover, 30% of the CRAB isolates in the EURECA study, predominantly from Serbia, harbored bla OXA-72 [26] . Additionally, the coexistence of bla OXA-72 and bla NDM-1 was found in four ST492 isolates from a Serbian hospital. Like OXA-23, OXA-72 CHDLs (chromosomally or plasmid-encoded) have been reported worldwide, including in Taiwan, China, Brazil, Colombia, the USA, France, Poland, Italy, Belgium, Lithuania, Sweden and Croatia [71,79].
The WGS of three selected A. baumannii isolates revealed that the resistance to carbapenems was mainly mediated by the presence of CHDLs (OXA-23, OXA-66 and OXA-72). Moreover, Aba52 and Aba176 also harbored ISAba24 and ISAba125 insertion sequences that were positioned in the promotor regions of the bla ADC-25 and bla OXA-23 genes, likely contributing to their transcriptional upregulation. Previous studies have highlighted ISAba125 as a particularly efficient element in this regard due to its better promoter sequence compared to ISAba1 [82]. The Aba190 isolate compensates for the absence of these elements with the presence of a bla OXA-72 gene, located on a small plasmid. This finding raised concerns since OXA-72 can be easily transmitted via horizontal gene transfer due to its localization. Moreover, the plasmid also harbored an ISAba31 element, which creates an additional possibility for gene transfer via recombinational events between insertion sequences.
Aminoglycosides have traditionally been included in therapeutic regimens for the treatment of infections caused by CRAB [15]. Unfortunately, aminoglycoside resistance has also increased in the recent years [70]. A total of 86.3% of the isolates tested in this study demonstrated high-level resistance to the three aminoglycosides (AMK, GEN and TOB), with MIC values above 256 mg/L. Methylation of the 16S rRNA in the 30S ribosomal subunit is recognized as the most common molecular mechanism, conferring high-level resistance to most clinically useful aminoglycosides [3,43]. All high-level-resistant A. baumannii isolates in this study were identified to be ArmA 16S rRNA methylase producers. A previous study found 32.9% ArmA-producing CRAB isolates among MDR-AB isolates obtained from ICU patients at a Bulgarian university hospital during the period of 2010-2011 [74]. The authors reported plasmid localization of the encoding gene, as well as coexistence of armA and bla OXA-23 . Moreover, the co-production of ArmA 16S rRNA methylase and OXA-23 carbapenemase was determined as the most prevalent mechanism conferring MDR in clinical A. baumannii isolates in many countries worldwide, such as Italy, Greece, China, South Korea, India, Yemen etc. [69,71,83].
Similar to other pathogenic bacteria, class 1, 2 and 3 integrons that harbored sul genes, responsible for sulfonamide resistance, have been found in nosocomial MDR-AB isolates [84,85]. As described above, all SXT-resistant CRAB isolates in the present study were sul1-positive and the three WGS-subjected ones were confirmed as class-1-carrying isolates.
Two of the three WGS analyzed XDR-AB isolates belonged to the high-risk ST2 (CC2, IC2) and one of them was classified into ST636 (a triple-locus variant of ST2). ST2 is a CRAB lineage that has been reported from various parts of the world [24,26,79,86]. Of note, ST2 was found to be the most common ST among the CRAB isolates (2016-2018), included in the EURECA multicenter study (67.7%) [26]. Additionally, a large-scale WGS-based study determined ST2 as the most prevalent ST in CRAB isolates (2013-2017) from three Mediterranean countries as follows: Greece (93.5%), Italy (82.9%) and Israel (40%) [86]. Recently, Cherubini et al. investigated A. baumannii isolates recovered from the blood cultures of 43 critically ill COVID-19 patients treated in the ICU ward of a hospital in Pescara, Central Italy. Interestingly, all isolates demonstrated a remarkable antimicrobial resistance, including MEM and COL, and belonged to ST2 high-risk international clone [87]. ST2 was predominant (58.1%), followed by ST636 (14%), in clinical A. baumannii isolates collected across Belgium 2014-2017. The collection included CRAB, MDR-AB, XDR-AB, and even PDR-AB isolates [79]. Recently, Gajic et al. reported an outbreak of BSI in preterm neonates caused by A. baumannii bla OXA-66 /bla OXA-72 /ST636 in a neonatal ICU in Serbia [88]. Moreover, OXA-72-producing CRAB isolates, belonging to ST636, were previously found in a nationwide study conducted in Serbia [89].
The phylogenomic analysis revealed that the Bulgarian isolates were closely related to those collected from Greece, Croatia and Romania. These results suggest that human migration and travel facilitate the spread of high-risk CRAB clones in the Balkans.

Conclusions
The studied nosocomial CRAB isolates also demonstrated cross-resistance to antimicrobial agents from other classes, which was the basis for categorizing them as MDR-AB (21.9%) or XDR-AB (78.1%). Antimicrobial treatment choices for these problematic isolates were severely limited. There were only a few effective options available, such as COL and SAM and, rarely, SXT.
Production of acquired CHDLs, predominantly OXA-23, was found in all isolates. The WGS revealed the presence of different ISs, such as ISAba24, ISAba31, ISAba125, ISVsa3, IS17, and IS6100, which, in addition to the most common ISAba1, allow A. baumannii the high ability to capture and mobilize ARGs. Moreover, the emergence of OXA-72, associated with a plasmid and transposable mobile genetic elements, such as ISAba31, increases the risk of horizontal transfer of ARGs among nosocomial A. baumannii isolates, as well as other Gram-negative pathogens in Bulgaria. A geographical difference between the isolates from Sofia and Varna was established, both in terms of the identified antimicrobial resistance determinants and mobile genetic elements, as well as their ST affiliations (high-risk ST2 and ST636 belonging to IC2).
To the best of our knowledge, this is the first WGS-based resistome analysis of A. baumannii conducted in Bulgaria. The implementation of a rapid technique, such as next-generation sequencing, to identify resistomes and clonal relatedness of problematic pathogens in hospitals could be useful for effective antimicrobial stewardship programs. Our results show a widespread dissemination of XDR-AB isolates, carrying a variety of ARGs in Bulgarian ICU settings, which highlights the crucial need for nationwide surveillance, especially in the conditions of extensive antibiotic usage during the COVID-19 pandemic.

Informed Consent Statement:
All patients signed the usual informed consent statement before hospital admission. In the present study, we worked only with bacterial isolates and did not use personal patient information.
Data Availability Statement: Whole-genome shotgun sequencing projects of the A. baumannii Aba52, Aba176 and Aba190 isolates have been deposited in GenBank under Accession numbers JARNMV000000000, JARNMW000000000 and JAR NMX000000000, respectively.

Conflicts of Interest:
The authors declare no conflict of interest.