Abundance of Colistin-Resistant, OXA-23- and ArmA-Producing Acinetobacter baumannii Belonging to International Clone 2 in Greece

Carbapenem resistant Acinetobacter baumannii (CRAB) represents one of the most challenging pathogens in clinical settings. Colistin is routinely used for treatment of infections by this pathogen, but increasing colistin resistance has been reported. We obtained 122 CRAB isolates from nine Greek hospitals between 2015 and 2017, and those colistin resistant (ColR; N = 40, 32.8%) were whole genome sequenced, also by including two colistin susceptible (ColS) isolates for comparison. All ColR isolates were characterized by a previously described mutation, PmrB^A226V, which was associated with low-level colistin resistance. Some isolates were characterized by additional mutations in PmrB (E140V or L178F) or PmrA (K172I or D10N), first described here, and higher colistin minimum inhibitory concentrations (MICs), up to 64 mg/L. Mass spectrometry analysis of lipid A showed the presence of a phosphoethanolamine (pEtN) moiety on lipid A, likely resulting from the PmrA/B-induced pmrC overexpression. Interestingly, also the two ColS isolates had the same lipid A modification, suggesting that not all lipid A modifications lead to colistin resistance or that other factors could contribute to the resistance phenotype. Most of the isolates (N = 37, 92.5%) belonged to the globally distributed international clone (IC) 2 and comprised four different sequence types (STs) as defined by using the Oxford scheme (ST 425, 208, 451, and 436). Three isolates belonged to IC1 and ST1567. All the genomes harbored an intrinsic bla_OXA–51 group carbapenemase gene, where bla_OXA–66 and bla_OXA–69 were associated with IC2 and IC1, respectively. Carbapenem resistance was due to the most commonly reported acquired carbapenemase gene bla_OXA–23, with ISAba1 located upstream of the gene and likely increasing its expression. The armA gene, associated with high-level resistance to aminoglycosides, was detected in 87.5% of isolates. Collectively, these results revealed a convergent evolution of different clonal lineages toward the same colistin resistance mechanism, thus limiting the effective therapeutic options for the treatment of CRAB infections.

Introduction

Acinetobacter baumannii is now recognized as a major hospital pathogen by its ability to resist major antimicrobials and to survive in the healthcare environment (Peleg et al., 2008). Currently, carbapenem resistant A. baumannii (CRAB) is widespread, with rates reaching or exceeding 90% in some clinical settings in Southern and Eastern European countries (ECDC, 2018) and elsewhere¹, and mortality rates for the most common CRAB infections such as bloodstream infections and hospital acquired pneumoniae approaching 60% (Wong et al., 2017). OXA-type carbapenemases constitute the most prevalent mechanism of carbapenem resistance in this species, with OXA-23, OXA-24, and OXA-58 being the most prevalent enzymes (Poirel and Nordmann, 2006). Molecular epidemiological studies usually revealed an oligoclonal distribution of CRAB, with outbreak strains mostly belonging to international clones (IC) 1 and 2 (Zarrilli et al., 2013).

Since their first emergence in 2000, CRAB have become endemic, and the percentage of carbapenem resistance reached 94% in 2017 (Tsakris et al., 2003; ECDC, 2018). Regarding the CRAB clonal nature and carbapenemase gene content, a study conducted from 2000 to 2009 in Greece showed that CRAB were harboring only the OXA-58 carbapenemase gene; while IC1 was prevalent until 2004, IC2 became dominant during 2005–2009 (Gogou et al., 2011). Between 2009 and 2011, OXA-23 producers emerged and replaced the previously predominant OXA-58 producing A. baumannii strains (Liakopoulos et al., 2012). Recently, a molecular epidemiological study on contemporary CRAB clinical isolates derived from hospitals throughout Greece demonstrated the predominance of OXA-23 producers belonging to IC2 (Pournaras et al., 2017).

Colistin-based treatment often represents the only therapeutic option for CRAB infections (Viehman et al., 2014). However, CRAB isolates that are also colistin resistant (ColR) are being reported more frequently. Data from the EARS-Net study in 2016 collected from 30 European countries showed that 4.0% of the tested isolates were resistant to colistin, with the vast majority (70.7%) of the resistant isolates reported from Greece and Italy (ECDC, 2017). A study from Greece reported an increase in colistin resistance from 1% in 2012 to 21.1% in 2014 (Oikonomou et al., 2015), while Pournaras et al. (2017) reported a resistance rate of 27.3% in 2015. More alarmingly, the colistin resistance rate was 56.8% in isolates collected from patients with ventilator-associated pneumonia in Greece during 2015 (Nowak et al., 2017). Colistin resistance has been linked to mutations in the two-component transcriptional regulator genes pmrA/B and consequent pmrC overexpression in most instances. The phosphoethanolamine (pEtN) phosphotransferase PmrC adds a pEtN group to the lipid A of the lipopolysaccharide, lowering the net negative charge of the cell membrane, thus impacting the binding of colistin and preventing the cell membrane leakage (Poirel et al., 2017). Colistin resistance may also result from the overexpression of etpA, a pmrC homolog. This is mediated by insertional inactivation of a gene encoding an H-NS family transcriptional regulator (Lucas et al., 2018) or by integration of insertion sequence elements upstream of the eptA gene itself (Gerson et al., 2019; Potron et al., 2019; Trebosc et al., 2019).

In this study, 40 ColR and two colistin susceptible (ColS) CRAB isolates collected from nine Greek hospitals between 2015 and 2017 were studied. Whole genome sequencing was performed to investigate the mechanisms of antibiotic resistance as well as the genomic relatedness between the strains.

Materials and Methods

Bacterial Strains and Antimicrobial Susceptibility Testing

In the period 2015–2017, a total of 122 consecutive non-duplicate clinical CRAB isolates were obtained from routine microbiological cultures of clinical samples (e.g., urine, blood, skin, bronchial aspirate) from different patients admitted to nine Greek hospitals involved in this study (Figure 1). Bacteria were not isolated by the authors but provided by the respective medical centers. Therefore, an ethics approval was not required as per institutional and national guidelines and regulations. Antimicrobial susceptibility testing was performed using the Vitek2 instrument (bioMérieux, Marcy l’Étoile, France) and the results were interpreted following the EUCAST breakpoints (EUCAST, 2019). Since EUCAST doesn’t provide breakpoints for cephalosporins and Acinetobacter spp., CLSI breakpoints were used for those antibiotics (CLSI, 2019). Colistin minimum inhibitory concentrations (MICs) were obtained by broth microdilution following the CLSI guidelines (CLSI, 2019), and the results were interpreted following the EUCAST susceptibility breakpoint of 2 mg/L (EUCAST, 2019). Only the ColR CRAB isolates plus two randomly selected ColS CRAB isolates were retained for further experiments.

FIGURE 1

Figure 1. Phylogenetic tree of the A. baumannii clinical isolates belonging to IC2. Cps: capsular polysaccharides.

Genome Sequencing and Assembly

Whole DNA of the selected CRAB isolates was extracted using the QIAGEN UltraClean Microbial kit and sequenced with a NovaSeq sequencer (Illumina, United States), generating paired end reads of 100 bp. Raw reads were assembled using SPAdes v.3.11.1 (Bankevich et al., 2012) and annotated with Prokka (Seemann, 2014). Whole genome sequencing data have been deposited under BioProject PRJNA578598.

Bioinformatics Analysis

Sequence types (STs) were assigned by the mlst tool² by using the Oxford (gltA, gyrB, gdhB, recA, cpn60, gpi, and rpoD genes) and the Pasteur (cpn60, fusA, gltA, pyrG, recA, rplB, and rpoB genes) schemes available on pubMLST.org. The ABRicate tool³ was used for the detection of antimicrobial resistance genes, by using the ResFinder (Zankari et al., 2012), CARD (Jia et al., 2017), BLDB (Naas et al., 2017), and ARG-ANNOT (Gupta et al., 2014) databases. The minimum percentage of coverage and identity used were 60 and 90%, respectively. The Kaptive tool was used to detect the KL and OC locus (Wyres et al., 2019). BLAST+ (2.7.1) was used to detect mutations in genes previously demonstrated to be potentially involved in colistin resistance (i.e., pmrCAB, eptA), and only those leading to amino acid variations were considered. The pmrA/B/C and eptA genes were compared to the reference genome ACICU (accession no. CP031380.1). The presence of insertion sequence elements in the 500 bp region upstream of the bla_ADC, bla_OXA–23, armA, eptA, and pmrC genes was determined using the ISfinder tool (Siguier et al., 2006). Core genes were defined by Roary (v3.12.0) (Page et al., 2015) by using the annotated genomes, and genomes belonging to different ICs were treated separately. The alignment of these genes was screened for further recombination using Gubbins (v2.3.4) (Croucher et al., 2015), while an ML phylogeny was obtained by using RAxML (v8.2.12) (Stamatakis, 2014) with the GTRGAMMA model and 100 bootstrap replicates. The phylogenetic tree was visualized together with associated metadata using Microreact (v7.0.0) (Argimón et al., 2016). Single nucleotide polymorphisms (SNPs) were obtained with the snp-dists tool⁴ by using the Roary core genes alignment as input.

Analysis of Lipid A

Lipid A was extracted using an acetic acid-based procedure as previously described (Kocsis et al., 2017). Once extracted, 0.7 μL of the concentrate was spotted on a matrix-assisted laser desorption/ionization–time of flight mass spectrometry (MALDI-TOF MS) plate followed by 0.7 μL of norharmane matrix (Sigma-Aldrich, St Louis, MO, United States) and then air-dried. The samples were analyzed on a Vitek MS instrument (bioMérieux, Marcy l’Étoile, France) in the negative-ion mode. The resulting spectra were compared to that obtained for the ColS reference strain A. baumannii ATCC 19606.

Results

Bacterial Strains and Antimicrobial Susceptibilities

Of the 122 CRAB isolates, 40 (32.8%) were also ColR, with colistin MICs ranging from 4 to 64 mg/L. All following data concern only the ColR isolates. Antimicrobial susceptibility testing revealed that all isolates were resistant to cephalosporins (ceftazidime and cefepime), carbapenems (imipenem and meropenem) fluoroquinolones (ciprofloxacin and levofloxacin) and tobramycin. Resistance rates for gentamicin and trimethoprim/sulfamethoxazole were 87.5% (N = 35) and 92.5% (N = 37), respectively. The two ColS CRAB isolates, included in this study for comparative purposes, had a colistin MIC of 0.5 mg/L (Supplementary Table S1).

Genomic Epidemiology

The majority of the ColR CRAB isolates (N = 37, 92.5%) were ST2, belonging to the previously described IC2 as defined by the Pasteur MLST scheme (Diancourt et al., 2010; Figure 1). The Oxford MLST scheme allowed to further differentiate the IC2 isolates in 4 different STs, all belonging to the clonal complex (CC) 208: the majority of isolates (N = 29, 78.4%) belonged to ST425, while ST208, ST451 and ST436 represented the 8.1% (N = 3), 8.1% (N = 3) and 5.4% (N = 2), respectively. These four STs shared six out of seven alleles and differed only by the gpi gene. The gpi gene is one of the capsular polysaccharide synthesis genes; therefore, the Oxford MLST scheme suffers from limitations, as the gpi gene is prone to homologous recombination (Gaiarsa et al., 2019). From a total of 4,612 different genes detected in all the isolates, 3,031 (65.7%) were core genes. Core gene SNPs among IC2 genomes varied between 2 and 1,652 (mean: 569, median: 531).

The phylogenetic analysis of IC2 isolates shows two major clusters of ST425, well differentiated in the tree. These two clusters were characterized by two different capsular polysaccharides, KL4 and KL40.

Different capsular polysaccharides were observed in the other IC2 isolates (Figure 1), while all the IC2 isolates were characterized by the lipooligosaccharide outer core (OC) locus 1 (OCL1).

Isolates of ST425:KL4 were only observed in Athens, within two hospitals in 2015 (Aglaia Kyriakoy and Agia Olga) and one isolate in the Thriassio General hospital in 2017. The 13 isolates obtained from the Aglaia Kyriakoy hospital had an average of 12 core SNPs, suggesting cross-transmission of isolates between different patients. Isolates of ST425:KL40 were retrieved between 2015 and 2016 from four hospitals in Athens and one isolate from the University Hospital in Patras (200 km west of Athens). This underscores the endemicity at the local level of this clone, moreover, suggesting inter-hospital cross infections, given the absence of a clear clustering in the tree of isolates obtained from different hospitals.

The remaining three isolates belonged to ST1 (IC1) and ST1567 according to the Pasteur and Oxford MLST schemes, respectively, and harbored a capsule and lipooligosaccharide of type KL40 and OCL2. Core gene SNPs varied between 29 and 305.

The two ColS CRAB isolates belonged to IC2, or ST208 (isolate PU_2016_41) and ST195 (GE_2017_62) by using the Oxford MLST scheme, and had a median of 706 (min: 8, max: 1,568) and 1175 SNPs (min: 725, max: 1,652) compared to the ColR isolates, respectively.

Colistin Resistance Mechanisms

Several chromosomal mutations in genes potentially involved in colistin resistance were detected, in comparison with the ACICU ColS reference genome. The mutation PmrB^A138T was detected in all ColR and ColS isolates, indicating that it may not contribute significantly to the resistance phenotype, as previously reported (Oikonomou et al., 2015). Conversely, the mutation A226V in the histidine kinase A (phosphoacceptor) domain of PmrB was observed in all ColR isolates, and not in the ColS ones (Figure 1). This mutation has been described in several prior studies, associated with ColR strains (Arroyo et al., 2011; Mavroidi et al., 2015, 2017; Dortet et al., 2018; Trebosc et al., 2019).

Isolates with PmrB^A226V without other alterations had colistin MICs ranging from 4 to 8 mg/L. When an additional mutation in PmrB occurred (PmrB^E140F in AK_2015_33 and PmrB^L178F in SI_2017_69), strains showed a colistin MIC of 16 mg/L. Two strains belonging to IC1 (FK_2016_46 and FK_2016_47) had an additional K172I mutation in the transcriptional regulatory protein C-terminal domain of PmrA, and showed colistin MICs of 32 mg/L. Finally, the strain AO_2015_54 had an additional D10N mutation in the CheY-homologous receiver PmrA domain and was associated with colistin MIC of 64 mg/L. All these additional mutations are, to the best of our knowledge, first described here.

The susceptible strain GE_2017_62 had no additional mutations in pmrA/B genes. However, it had an ISAba1 positioned 110 bp upstream of the pmrC gene, in reverse orientation. This is, to the best of our knowledge, the first report of an insertion sequence transposition upstream of the pmrC gene. However, this transposition event doesn’t seem to alter the colistin susceptibility in this isolate. The second susceptible strain PU_2016_41 had PmrA^M12V and PmrB^R181H+Y388N. These mutations are firstly described here, and in this strain they don’t seem to impact the colistin susceptibility.

The pmrC homolog eptA was detected in all the isolates of the IC2 except the susceptible isolate GE_2017_62, while it was absent in the IC1 isolates. The obtained eptA gene sequences were identical to that of the susceptible reference ACICU, and did not present insertion elements in the upstream region.

The mcr genes, encoding for acquired colistin resistance, have not been described in A. baumannii yet, and were not detected in our strain collection.

Lipid A Modifications

An increased expression of pmrC or eptA results to the addition of pEtN to lipid A. The lipid A of the ColR and ColS CRAB isolates was extracted and analyzed by MALDI-TOF MS, and the resulting spectra were compared to that of the ColS reference strain A. baumannii ATCC 19606.

Several lipid A species were detected in the reference strain ATCC 19606 and in all clinical isolates: hepta-acylated lipid A (m/z 1,910), hexa-acylated lipid A (m/z 1,728) and tetra-acylated lipid A (m/z 1,404). The addition of pEtN (m/z 124) to lipid A was shown by the mass at m/z 2,034, and unexpectedly it was observed in all the clinical strains, including the ColS ones (Figure 2).

FIGURE 2

Figure 2. Mass spectrometry analysis of Lipid A. From the bottom, isolate (colistin MIC, resistant/susceptible): ATCC-19606 (0.5, S), TR_2016_35 (4, R), SI_2017_69 (16, R), AO_2015_54 (64, R) and GE_2017_62 (0.5, S).

Isolates with colistin MICs higher than 8 mg/L also showed the peak at m/z 1,954, representing the pEtN-modified hepta-acylated lipid A (m/z 2,034) minus one phosphate group (m/z 80), as previously reported (Kim et al., 2014).

The addition of galactosamine to lipid A, which is indicated by a mass at m/z 2,071 (Pelletier et al., 2013), was not observed in any isolate.

Antimicrobial Resistance Mechanisms and Phenotype Correlation

All CRAB genomes harbored a chromosomal bla_ADC cephalosporinase, an intrinsic bla_OXA–51 group carbapenemase and an acquired bla_OXA–23. IC2 genomes contained the bla_ADC–73 (accession no. KP881233), a variant of bla_ADC with a sequence identity of 1,151/1,152 nucleotides compared to that of bla_ADC–30, and previously observed in IC2 isolates (Karah et al., 2016). An ISAba1 element was present 9 bp upstream of the bla_ADC–73 gene in reverse orientation in all IC2 genomes, and it is responsible to increase the cephalosporinase gene expression (Héritier et al., 2006). Conversely, IC1 genomes contained bla_ADC–175 (MH594297) with an ISAba125 element positioned 66 bp upstream the gene in reverse orientation, as also previously reported (Lopes and Amyes, 2012). ISAba125 was shown to increase the cephalosporinase expression 6 times more than ISAba1 (Lopes and Amyes, 2012). The allelic variants of the intrinsic bla_OXA–51-like carbapenemase genes were bla_OXA–66 and bla_OXA–69, associated with IC2 and IC1, respectively, as previously observed (Zander et al., 2012). All the bla_OXA–23 genes were characterized by the presence of an ISAba1 located upstream of the gene, which has been previously demonstrated to increase its expression (Turton et al., 2006). In particular, the bla_OXA–23 gene was part of a Tn2006 transposon in the IC2 genomes. Conversely, a Tn2008 embedded within a TnaphA6 was found in the three IC1 genomes, matching 100% with the sequence of plasmid pABKp1 (KP074966.1) obtained from A. baumannii isolates from Romania (Gheorghe et al., 2014). Consistently with the mentioned genes and their genetic environments, all isolates were resistant to cephalosporins, including ceftazidime (3rd generation) and cefepime (4th generation), and carbapenems (imipenem and meropenem).

Aminoglycoside resistance genes were observed among the isolates, namely aac(3)-I, aac(3)-Ia, ant(3″)-1a, aph(3′)-Ia, aph(3′)-Via, aph(6)-Id, armA, and strA (Supplementary Table S1). ArmA is a 16S ribosomal RNA methyltransferase, which protects the 30S ribosomal subunit from aminoglycoside binding and conferring high aminoglycosides MICs. Consistently, all the strains carrying armA (35/40, 87.5%) were resistant to gentamicin and tobramycin. The armA gene was located in the chromosome aboard on the widely disseminated Tn1548, and it was found downstream of a cluster of genes encoding proteins annotated as paraquat-inducible protein A and protein B, as previously described for ST195 strain AC29 (Lean et al., 2016).

All strains contained substitutions within the QRDR, namely GyrA^S83L and ParC^S80L, previously associated to quinolone resistance (Vila et al., 1995, 1997). As expected, all strains were non-susceptible to ciprofloxacin and levofloxacin.

Discussion

Carbapenems represent first-line agents for the treatment of A. baumannii infections, consequently the rise of infections due to carbapenem-resistant strains is of major concern. The carbapenem resistance in the isolates described here was associated to the ISAba1-mediated overexpression of bla_OXA–23 located either in Tn2006 (IC2 isolates) or Tn2008 (IC1) transposons. Previous studies reported that OXA-23 producers emerged and replaced the previously predominant OXA-58 A. baumannii isolates (Liakopoulos et al., 2012), and this phenomenon could be linked to the stronger hydrolytic activity of OXA-23 compared to OXA-58 (Peleg et al., 2008).

Additionally, most CRAB isolates are susceptible to only 1 or 2 agents, making them extensively drug-resistant (XDR) pathogens (Viehman et al., 2014).

Because of the increasing use of colistin, resistance to this antibiotic has rapidly increased, especially in CRAB isolates (Giamarellou, 2016; Jeannot et al., 2017), and now reached critical levels in some countries (Nowak et al., 2017). From the nine hospitals involved in this study, the 32.8% of the CRAB isolates were also ColR. These results indicate that colistin resistance rates among CRAB isolates from Greece is on the rise, as a previous study reported a resistance rate of 27.3% in 2015 (Pournaras et al., 2017). While the mcr genes, encoding for acquired colistin resistance, were absent among our isolates, we found several mutations in the pmrCAB operon associated with the colistin resistance phenotype. Interestingly, the previously described PmrB^A226V mutation, previously associated to low-level colistin resistance, was detected in all the ColR isolates but no in the ColS ones. In a recent study, Trebosc et al. (2019) investigated the colistin resistance mechanisms of 12 clinical A. baumannii strains. The authors concluded that colistin resistance was conferred, in most cases, by mutations in the PmrB sensor kinase that led to PmrC overexpression. Two of those strains were isolated in Greece in 2012, belonged to either IC1 or IC2 and had the mutation PmrB^A226V. Such findings support the important role of the mentioned PmrB mutation in the colistin resistance phenotype. Moreover, a similar substitution of the alanine in position 226 of PmrB was reported to confer stable colistin resistance in clinical A. baumannii isolates (Charretier et al., 2018). Some of our isolates had additional mutations in either PmrB or PmrA, and were associated with higher colistin MICs, up to 64 mg/L. Multiple mutations may result in an increased expression of pmrC, as recently shown by RNA-Seq experiments (Wright et al., 2017) and by qRT-PCR (Gerson et al., 2020). However, the same studies reported clinical isolates characterized by pmrC overexpression due to PmrA/B mutations, but with an unexpected ColS phenotype. Similarly, the two ColS isolates from our study had pmrCAB alterations and a modified lipid A, as observed with mass spectrometry. All these observations support the hypothesis that additional and still unknown factors are involved in colistin resistance of clinical A. baumannii isolates (Jeannot et al., 2017; Gerson et al., 2019, 2020). Determination of the cell-envelope charge could be useful in the elucidation process of the complex mechanism of colistin resistance in A. baumannii (Cafiso et al., 2019).

In this study, we observed a clear predominance of IC2, which is globally distributed (Higgins et al., 2010) and which is gradually replacing IC1 (Gogou et al., 2011; Villalon et al., 2011). The major sequence type within IC2 was ST425, as defined by the Oxford MLST scheme. To the best of our knowledge, only one study reported such ST, with one clinical isolate collected in 2002 in Sydney, Australia (Nigro and Hall, 2016). However, WGS data were not provided. Both capsular polysaccharides reported within our ST425 isolates, KL4 and KL40, were rarely observed (0.2%) or completely absent, respectively, within IC2 genomes in a recent study where 3,416 publicly available A. baumannii genomes were analyzed (Wyres et al., 2019). Conversely, KL4 and KL40 represented the second (20.1%) and third (11.9%) most common capsular types observed within IC1 genomes. It is conceivable that ST425 resulted from homologous recombination between a CC208 and an IC1 genomes, where the IC1 capsular polysaccharides genes were acquired by the CC208 strain, as this region was previously shown to be a frequent subject of homologous recombination (Adams et al., 2008; Snitkin et al., 2011; Kenyon and Hall, 2013).

Conclusion

Genomic analysis of ColR CRAB isolates from different Greek hospitals revealed a convergent evolution of different clonal lineages toward the same colistin resistance mechanism, characterized by the mutation PmrB^A226V. The prevalence of ColR CRAB isolates belonging to IC2 and expressing OXA-23 and ArmA is increasing, and it represents a huge threat within clinical settings, given the very limited effective agents for the treatment of infections caused by such strains.

Data Availability Statement

The datasets generated for this study can be found in the PRJNA578598.

Ethics Statement

Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

AB, GR, and HG conceived and designed the study. MP and NP performed the phenotypic and genomics experiments. MP and MD’A performed the bioinformatics analysis. All authors analyzed the data and contributed to the manuscript.

Funding

The work performed by MP and AP was funded by the European Union’s Horizon 2020 Research and Innovation Program under the Marie Skłodowska-Curie Grant Agreement No 675412 [New Diagnostics for Infectious Diseases (ND4ID)].

Conflict of Interest

NL was employed by IASO Group Hospitals.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We would like to thank the members of the study group for epidemiology of carbapenem and colistin resistant Acinetobacter baumannii for providing the strain collection: S. Tsiplakou and V. Papaioannou (KAT University Hospital, Athens); E. Trikka-Graphakos, N. Charalampakis and C. Sereti (Thriassio General Hospital, Athens); E. Lebessi and E. Bozavoutoglou (Aglaia Kyriakou Children’s Hospital, Athens); A. Xanthaki (Ippokratio University Hospital, Athens); E. Peteinaki (General University Hospital, Larissa); S. Chatzipanagiotou and A. Ioannidis (Aeginitio University Hospital, Athens); E. Platsouka (Agia Olga General Hospital, Athens); M. Skarmoutsou (Sismanogleio General Hospital, Athens); I. Spiliopoulou and M. Christofidou (University Hospital, Patras); N. Legakis and T. Skalidis (Central Laboratories, IASO Group Hospitals, Athens, Greece).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.00668/full#supplementary-material

TABLE S1 | Metadata, genotypic and phenotypic results of the 42 CRAB from Greece.

Footnotes

1. ^ https://resistancemap.cddep.org
2. ^ https://github.com/tseemann/mlst
3. ^ https://github.com/tseemann/abricate
4. ^ https://github.com/tseemann/snp-dists

References

Adams, M. D., Goglin, K., Molyneaux, N., Hujer, K. M., Lavender, H., Jamison, J. J., et al. (2008). Comparative genome sequence analysis of multidrug-resistant Acinetobacter baumannii. J. Bacteriol. 190, 8053–8064. doi: 10.1128/JB.00834-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Argimón, S., Abudahab, K., Goater, R. J. E., Fedosejev, A., Bhai, J., Glasner, C., et al. (2016). Microreact: visualizing and sharing data for genomic epidemiology and phylogeography. Microb. Genom. 2:e000093. doi: 10.1099/mgen.0.000093

PubMed Abstract | CrossRef Full Text | Google Scholar

Arroyo, L. A., Herrera, C. M., Fernandez, L., Hankins, J. V., Trent, M. S., and Hancock, R. E. W. (2011). The pmrCAB operon mediates polymyxin resistance in Acinetobacter baumannii ATCC 17978 and clinical isolates through phosphoethanolamine modification of lipid A. Antimicrob. Agents Chemother. 55, 3743–3751. doi: 10.1128/AAC.00256-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021

PubMed Abstract | CrossRef Full Text | Google Scholar

Cafiso, V., Stracquadanio, S., Lo Verde, F., Gabriele, G., Mezzatesta, M. L., Caio, C., et al. (2019). Colistin resistant A. baumannii: genomic and transcriptomic traits acquired under colistin therapy. Front. Microbiol. 9:3195. doi: 10.3389/fmicb.2018.03195

CrossRef Full Text | Google Scholar

Charretier, Y., Diene, S. M., Baud, D., Chatellier, S., Santiago-Allexant, E., Van Belkum, A., et al. (2018). Colistin heteroresistance and involvement of the PmrAB regulatory system in Acinetobacter baumannii. Antimicrob. Agents Chemother. 62:e788-18. doi: 10.1128/AAC.00788-18

PubMed Abstract | CrossRef Full Text | Google Scholar

CLSI (2019). Performance Standards for Antimicrobial Susceptibility Testing. 29th ed. CLSI supplement M100. Wayne, PA: Clinical and Laboratory Standars Institute.

Google Scholar

Croucher, N. J., Page, A. J., Connor, T. R., Delaney, A. J., Keane, J. A., Bentley, S. D., et al. (2015). Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 43:e15. doi: 10.1093/nar/gku1196

PubMed Abstract | CrossRef Full Text | Google Scholar

Diancourt, L., Passet, V., Nemec, A., Dijkshoorn, L., and Brisse, S. (2010). The population structure of Acinetobacter baumannii: expanding multiresistant clones from an ancestral susceptible genetic pool. PLoS ONE 5:e10034. doi: 10.1371/journal.pone.0010034

PubMed Abstract | CrossRef Full Text | Google Scholar

Dortet, L., Potron, A., Bonnin, R. A., Plesiat, P., Naas, T., Filloux, A., et al. (2018). Rapid detection of colistin resistance in Acinetobacter baumannii using MALDI-TOF-based lipidomics on intact bacteria. Sci. Rep. 8:16910. doi: 10.1038/s41598-018-35041-y

PubMed Abstract | CrossRef Full Text | Google Scholar

ECDC (2017). Surveillance of antimicrobial resistance in Europe 2016 Annual Report of the European Antimicrobial Resistance Surveillance Network (EARS-Net). Stockholm: ECDC, 2017.

Google Scholar

ECDC (2018). Surveillance of antimicrobial resistance in Europe – Annual report of the European Antimicrobial Resistance Surveillance Network (EARS-Net) 2017, Stockholm: ECDC2018.

Google Scholar

EUCAST (2019). The European Committee on Antimicrobial Susceptibility Testing. Breakpoint tables for interpretation of MICs and zone diameters. Version 9.0. Available online at:http://www.eucast.org.

Google Scholar

Gaiarsa, S., Batisti Biffignandi, G., Esposito, E. P., Castelli, M., Jolley, K. A., Brisse, S., et al. (2019). Comparative analysis of the two Acinetobacter baumannii multilocus sequence typing (MLST) schemes. Front. Microbiol. 10:930. doi: 10.3389/fmicb.2019.00930

CrossRef Full Text | Google Scholar

Gerson, S., Betts, J. W., Lucaßen, K., Nodari, C. S., Wille, J., Josten, M., et al. (2019). Investigation of novel pmrB and eptA mutations in isogenic Acinetobacter baumannii isolates associated with colistin resistance and increased virulence In Vivo. Antimicrob. Agents Chemother. 63, 1–15. doi: 10.1128/aac.01586-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerson, S., Lucaßen, K., Wille, J., Nodari, C. S., Stefanik, D., Nowak, J., et al. (2020). Diversity of amino acid substitutions in PmrCAB associated with colistin resistance in clinical isolates of Acinetobacter baumannii. Int. J. Antimicrob. Agents. 55:105862. doi: 10.1016/j.ijantimicag.2019.105862

PubMed Abstract | CrossRef Full Text | Google Scholar

Gheorghe, I., Novais, Â, Grosso, F., Rodrigues, C., Chifiriuc, M. C., Lazar, V., et al. (2014). Snapshot on carbapenemase-producing Pseudomonas aeruginosa and Acinetobacter baumannii in bucharest hospitals reveals unusual clones and novel genetic surroundings for blaOXA-23. J. Antimicrob. Chemother. 70, 1016–1020. doi: 10.1093/jac/dku527

PubMed Abstract | CrossRef Full Text | Google Scholar

Giamarellou, H. (2016). Epidemiology of infections caused by polymyxin-resistant pathogens. Int. J. Antimicrob. Agents 48, 614–621. doi: 10.1016/j.ijantimicag.2016.09.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Gogou, V., Pournaras, S., Giannouli, M., Voulgari, E., Piperaki, E.-T., Zarrilli, R., et al. (2011). Evolution of multidrug-resistant Acinetobacter baumannii clonal lineages: a 10 year study in Greece (2000–09). J. Antimicrob. Chemother. 66, 2767–2772. doi: 10.1093/jac/dkr390

PubMed Abstract | CrossRef Full Text | Google Scholar

Gupta, S. K., Padmanabhan, B. R., Diene, S. M., Lopez-Rojas, R., Kempf, M., Landraud, L., et al. (2014). ARG-annot, a new bioinformatic tool to discover antibiotic resistance genes in bacterial genomes. Antimicrob. Agents Chemother. 58, 212–220. doi: 10.1128/AAC.01310-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Héritier, C., Poirel, L., and Nordmann, P. (2006). Cephalosporinase over-expression resulting from insertion of ISAba1 in Acinetobacter baumannii. Clin. Microbiol. Infect. 12, 123–130. doi: 10.1111/j.1469-0691.2005.01320.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Higgins, P. G., Dammhayn, C., Hackel, M., and Seifert, H. (2010). Global spread of carbapenem-resistant Acinetobacter baumannii. J. Antimicrob. Chemother. 65, 233–238. doi: 10.1093/jac/dkp428

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeannot, K., Bolard, A., and Plésiat, P. (2017). Resistance to polymyxins in Gram-negative organisms. Int. J. Antimicrob. Agents 49, 526–535. doi: 10.1016/j.ijantimicag.2016.11.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Jia, B., Raphenya, A. R., Alcock, B., Waglechner, N., Guo, P., Tsang, K. K., et al. (2017). CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 45, 566–573. doi: 10.1093/nar/gkw1004

PubMed Abstract | CrossRef Full Text | Google Scholar

Karah, N., Dwibedi, C. K., Sjöström, K., Edquist, P., Johansson, A., Wai, S. N., et al. (2016). Novel aminoglycoside resistance transposons and transposon-derived circular forms detected in carbapenem-resistant Acinetobacter baumannii clinical isolates. Antimicrob. Agents Chemother. 60, 1801–1818. doi: 10.1128/AAC.02143-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Kenyon, J. J., and Hall, R. M. (2013). Variation in the complex carbohydrate biosynthesis loci of Acinetobacter baumannii genomes. PLoS ONE 8:e62160. doi: 10.1371/journal.pone.0062160

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, Y., Bae, I. K., Lee, H., Jeong, S. H., Yong, D., and Lee, K. (2014). In vivo emergence of colistin resistance in Acinetobacter baumannii clinical isolates of sequence type 357 during colistin treatment. Diagn. Microbiol. Infect. Dis. 79, 362–366. doi: 10.1016/j.diagmicrobio.2014.03.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Kocsis, B., Kilár, A., Péter, S., Dörnyei, Á, Sándor, V., and Kilár, F. (2017). “Mass spectrometry for profiling LOS and lipid A structures from whole-cell lysates: directly from a few bacterial colonies or from liquid broth cultures,” in Methods in Molecular Biology (Clifton, N.J.), ed. J. M. Walker (Berlin: Springer), 187–198. doi: 10.1007/978-1-4939-6958-6_17

PubMed Abstract | CrossRef Full Text | Google Scholar

Lean, S. S., Yeo, C. C., Suhaili, Z., and Thong, K. L. (2016). Comparative genomics of two ST 195 carbapenem-resistant Acinetobacter baumannii with different susceptibility to polymyxin revealed underlying resistance mechanism. Front. Microbiol. 6:1445. doi: 10.3389/fmicb.2015.01445

CrossRef Full Text | Google Scholar

Liakopoulos, A., Miriagou, V., Katsifas, E. A., Karagouni, A. D., Daikos, G. L., Tzouvelekis, L. S., et al. (2012). Identification of OXA-23-producing Acinetobacter baumannii in Greece, 2010 to 2011. Euro Surveill. 17:20117.

Google Scholar

Lopes, B. S., and Amyes, S. G. B. (2012). Role of ISAba1 and ISAba125 in governing the expression of blaADC in clinically relevant Acinetobacter baumannii strains resistant to cephalosporins. J. Med. Microbiol. 61, 1103–1108. doi: 10.1099/jmm.0.044156-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Lucas, D. D., Crane, B., Wright, A., Han, M.-L., Moffatt, J., Bulach, D., et al. (2018). Emergence of high-level colistin resistance in an Acinetobacter baumannii clinical isolate mediated by inactivation of the global regulator H-NS. Antimicrob. Agents Chemother 30, 1–17. doi: 10.1128/AAC.02442-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Mavroidi, A., Katsiari, M., Palla, E., Likousi, S., Roussou, Z., Nikolaou, C., et al. (2017). Investigation of extensively drug-resistant blaOXA-23-producing Acinetobacter baumannii spread in a greek hospital. Microb. Drug Resist. 23, 488–493. doi: 10.1089/mdr.2016.0101

PubMed Abstract | CrossRef Full Text | Google Scholar

Mavroidi, A., Likousi, S., Palla, E., Katsiari, M., Roussou, Z., Maguina, A., et al. (2015). Molecular identification of tigecycline- and colistin-resistant carbapenemase-producing Acinetobacter baumannii from a Greek hospital from 2011 to 2013. J. Med. Microbiol. 64, 993–997. doi: 10.1099/jmm.0.000127

PubMed Abstract | CrossRef Full Text | Google Scholar

Naas, T., Oueslati, S., Bonnin, R. A., Dabos, M. L., Zavala, A., Dortet, L., et al. (2017). Beta-lactamase database (BLDB) – structure and function. J. Enzyme Inhib. Med. Chem. 32, 917–919. doi: 10.1080/14756366.2017.1344235

PubMed Abstract | CrossRef Full Text | Google Scholar

Nigro, S. J., and Hall, R. M. (2016). Loss and gain of aminoglycoside resistance in global clone 2 Acinetobacter baumannii in Australia via modification of genomic resistance islands and acquisition of plasmids. J. Antimicrob. Chemother. 71, 2432–2440. doi: 10.1093/jac/dkw176

PubMed Abstract | CrossRef Full Text | Google Scholar

Nowak, J., Zander, E., Stefanik, D., Higgins, P. G., Roca, I., Vila, J., et al. (2017). High incidence of pandrug-resistant Acinetobacter baumannii isolates collected from patients with ventilator-associated pneumonia in Greece, Italy and Spain as part of the MagicBullet clinical trial. J. Antimicrob. Chemother. 72, 3277–3282. doi: 10.1093/jac/dkx322

PubMed Abstract | CrossRef Full Text | Google Scholar

Oikonomou, O., Sarrou, S., Papagiannitsis, C. C., Georgiadou, S., Mantzarlis, K., Zakynthinos, E., et al. (2015). Rapid dissemination of colistin and carbapenem resistant Acinetobacter baumannii in Central Greece: mechanisms of resistance, molecular identification and epidemiological data. BMC Infect. Dis. 15:559. doi: 10.1186/s12879-015-1297-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Page, A. J., Cummins, C. A., Hunt, M., Wong, V. K., Reuter, S., Holden, M. T. G., et al. (2015). Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics 31, 3691–3693. doi: 10.1093/bioinformatics/btv421

PubMed Abstract | CrossRef Full Text | Google Scholar

Peleg, A. Y., Seifert, H., and Paterson, D. L. (2008). Acinetobacter baumannii: emergence of a successful pathogen. Clin. Microbiol. Rev. 21, 538–582. doi: 10.1128/CMR.00058-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Pelletier, M. R., Casella, L. G., Jones, J. W., Adams, M. D., Zurawski, D. V., Hazlett, K. R. O., et al. (2013). Unique structural modifications are present in the lipopolysaccharide from colistin-resistant strains of Acinetobacter baumannii. Antimicrob. Agents Chemother. 57, 4831–4840. doi: 10.1128/AAC.00865-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Poirel, L., Jayol, A., and Nordmann, P. (2017). Polymyxins: antibacterial activity, susceptibility testing, and resistance mechanisms encoded by plasmids or chromosomes. Clin. Microbiol. Rev. 30, 557–596.

Google Scholar

Poirel, L., and Nordmann, P. (2006). Carbapenem resistance in Acinetobacter baumannii: mechanisms and epidemiology. Clin. Microbiol. Infect. 12, 826–836. doi: 10.1111/j.1469-0691.2006.01456.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Potron, A., Vuillemenot, J.-B., Puja, H., Triponney, P., Bour, M., Valot, B., et al. (2019). ISAba1-dependent overexpression of eptA in clinical strains of Acinetobacter baumannii resistant to colistin. J. Antimicrob. Chemother. 74, 2544–2550. doi: 10.1093/jac/dkz241

PubMed Abstract | CrossRef Full Text | Google Scholar

Pournaras, S., Dafopoulou, K., Del Franco, M., Zarkotou, O., Dimitroulia, E., Protonotariou, E., et al. (2017). Predominance of international clone 2 OXA-23-producing- Acinetobacter baumannii clinical isolates in Greece, 2015: results of a nationwide study. Int. J. Antimicrob. Agents 49, 749–753. doi: 10.1016/j.ijantimicag.2017.01.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069. doi: 10.1093/bioinformatics/btu153

PubMed Abstract | CrossRef Full Text | Google Scholar

Siguier, P., Perochon, J., Lestrade, L., Mahillon, J., and Chandler, M. (2006). ISfinder: the reference centre for bacterial insertion sequences. Nucleic Acids Res. 34, D32–D36. doi: 10.1093/nar/gkj014

PubMed Abstract | CrossRef Full Text | Google Scholar

Snitkin, E. S., Zelazny, A. M., Montero, C. I., Stock, F., Mijares, L., Mullikin, J., et al. (2011). Genome-wide recombination drives diversification of epidemic strains of Acinetobacter baumannii. Proc. Natl. Acad. Sci. U.S.A. 108, 13758–13763. doi: 10.1073/pnas.1104404108

PubMed Abstract | CrossRef Full Text | Google Scholar

Stamatakis, A. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313. doi: 10.1093/bioinformatics/btu033

PubMed Abstract | CrossRef Full Text | Google Scholar

Trebosc, V., Gartenmann, S., Tötzl, M., Lucchini, V., Schellhorn, B., Pieren, M., et al. (2019). Dissecting colistin resistance mechanisms in extensively drug-resistant Acinetobacter baumannii clinical isolates. MBio 10:e01083-19. doi: 10.1128/mBio.01083-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsakris, A., Tsioni, C., Pournaras, S., Polyzos, S., Maniatis, A. N., and Sofianou, D. (2003). Spread of low-level carbapenem-resistant Acinetobacter baumannii clones in a tertiary care Greek hospital. J. Antimicrob. Chemother. 52, 1046–1047. doi: 10.1093/jac/dkg470

PubMed Abstract | CrossRef Full Text | Google Scholar

Turton, J. F., Ward, M. E., Woodford, N., Kaufmann, M. E., Pike, R., Livermore, D. M., et al. (2006). The role of ISAba1 in expression of OXA carbapenemase genes in Acinetobacter baumannii. FEMS Microbiol. Lett. 258, 72–77. doi: 10.1111/j.1574-6968.2006.00195.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Viehman, J. A., Nguyen, M. H., and Doi, Y. (2014). Treatment options for carbapenem-resistant and extensively drug-resistant Acinetobacter baumannii infections. Drugs 74, 1315–1333. doi: 10.1007/s40265-014-0267-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Vila, J., Ruiz, J., Goñi, P., and Jimenez de Anta, T. (1997). Quinolone-resistance mutations in the topoisomerase IV parC gene of Acinetobacter baumannii. J. Antimicrob. Chemother. 39, 757–762.

Google Scholar

Vila, J., Ruiz, J., Goñi, P., Marcos, A., and Jimenez de Anta, T. (1995). Mutation in the gyrA gene of quinolone-resistant clinical isolates of Acinetobacter baumannii. Antimicrob. Agents Chemother. 39, 1201–1203.

Google Scholar

Villalon, P., Valdezate, S., Medina-Pascual, M. J., Rubio, V., Vindel, A., and Saez-Nieto, J. A. (2011). Clonal diversity of nosocomial epidemic Acinetobacter baumannii strains isolated in Spain. J. Clin. Microbiol. 49, 875–882. doi: 10.1128/JCM.01026-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Wong, D., Nielsen, T. B., Bonomo, R. A., Pantapalangkoor, P., Luna, B., and Spellberg, B. (2017). Clinical and pathophysiological overview of Acinetobacter infections: a century of challenges. Clin. Microbiol. Rev. 30, 409–447.

Google Scholar

Wright, M. S., Jacobs, M. R., Bonomo, R. A., and Adams, M. D. (2017). Transcriptome remodeling of Acinetobacter baumannii during infection and treatment. MBio 8:e2193-16. doi: 10.1128/mBio.02193-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Wyres, K. L., Cahill, S. M., Holt, K. E., Hall, R. M., and Kenyon, J. J. (2019). Identification of Acinetobacter baumannii loci for capsular polysaccharide (KL) and lipooligosaccharide outer core (OCL) synthesis in genome assemblies using curated reference databases compatible with Kaptive. Microb Genom. doi: 10.1101/869370 [Epub ahead of print].

CrossRef Full Text | Google Scholar

Zander, E., Nemec, A., Seifert, H., and Higgins, P. G. (2012). Association between β-lactamase-encoding blaOXA-51 variants and DiversiLab rep-PCR-based typing of Acinetobacter baumannii isolates. J. Clin. Microbiol. 50, 1900–1904. doi: 10.1128/JCM.06462-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Zankari, E., Hasman, H., Cosentino, S., Vestergaard, M., Rasmussen, S., Lund, O., et al. (2012). Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother. 67, 2640–2644. doi: 10.1093/jac/dks261

PubMed Abstract | CrossRef Full Text | Google Scholar

Zarrilli, R., Pournaras, S., Giannouli, M., and Tsakris, A. (2013). Global evolution of multidrug-resistant Acinetobacter baumannii clonal lineages. Int. J. Antimicrob. Agents 41, 11–19. doi: 10.1016/j.ijantimicag.2012.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar