- Open Access
Association of virulence gene expression with colistin-resistance in Acinetobacter baumannii: analysis of genotype, antimicrobial susceptibility, and biofilm formation
Annals of Clinical Microbiology and Antimicrobials volume 17, Article number: 24 (2018)
Acinetobacter baumannii causes difficult-to-treat nosocomial infections, which often lead to morbidity due to the development of antimicrobial drug resistance and expression of virulence genes. Data regarding the association of resistance to colistin, a last treatment option, and the virulence gene expression of A. baumannii is scarce.
We evaluated the MLVA genotype, antimicrobial resistance, and biofilm formation of 100 A. baumannii isolates from burn patients, and further compared the in vitro and in vivo expression of four virulence genes among five colistin-resistant A. baumannii (Cst-R-AB) isolates. Five Cst-R-AB isolates were tested; one from the present study, and four isolated previously.
Our results showed that reduced expression of recA, along with increased in vivo expression of lpsB, dnaK, and blsA; are associated with colistin resistance among Cst-R-AB isolates. Differences in virulence gene expressions among Cst-R-AB isolates, may in part explain common discrepant in vitro vs. in vivo susceptibility data during treatment of infections caused by Cst-R-AB.
Our findings highlight the intricate relationship between colistin-resistance and virulence among A. baumannii isolates, and underscore the importance of examining the interactions between virulence and antimicrobial resistance toward efforts to control the spread of multidrug-resistant A. baumannii (MDR-AB) isolates, and also to reduce disease severity in burn patients with MDR-AB infection.
Acinetobacter baumannii is an opportunistic pathogen that can cause formidable infections among patients with burn wounds worldwide . Also, the World Health Organization (WHO) has announced A. baumannii as “Critical” priority pathogenic bacteria that pose the greatest threat to human health. Treatment of A. baumannii wound infection is often difficult; primarily due to the pathogenic factors, which enable the establishment of persistent infections within burn patients leading to high morbidity and mortality [2,3,4]. These pathogenic factors mainly include the development of multi-drug resistance, and production of virulence determinants, such as biofilm formation [5,6,7]. Reports from various parts of the world have indicated a worrisome growing trend of isolating multi-, extensive- and pan-drug resistant (MDR, XDR, and PDR) strains of A. baumannii, some of which are even resistant to colistin, a last resort drug [8,9,10,11,12,13]. Whether the expression of virulence determinants, such as genes involved in biofilm formation, play an important role in the high rate of treatment failure among wound infections caused by highly resistant A. baumannii isolates remains to be explored .
To date, the expression of several A. baumannii virulence genes has been linked to the persistence and enhanced survival of A. baumannii within the host, including quorum sensing genes , and genes that confer biofilm production, which ultimately lead to increased antimicrobial resistance [16, 17]. Furthermore, the activation of A. baumannii virulence determinant genes, such as recA and dnaK, can render robust strains, which are less vulnerable to host response stresses . Moreover, the expression of recA by A. baumannii isolates increases resistance to stress , while recA inactivation increases susceptibility to a variety of antimicrobial agents [15, 20].
Then again, mounting evidence suggests that increased virulence gives rise to isolates that are less “fit” to survive in their host, and renders them more susceptible to antimicrobial agents . For instance, the loss of glycosyl transferase gene (lpsB)], essential for A. baumannii lipopolysaccharide (LPS) core biosynthesis, corresponds to the phenotypes of “reduced survival rate”, attenuated biofilm formation, and increased antimicrobial susceptibility [22,23,24]. Similarly, the expression of blue light sensing (blsA) gene and production of BlsA protein has been demonstrated to inhibit biofilm formation, which also leads to increased susceptibility to antimicrobial agents [25, 26]. In addition, among A. baumannii clinical isolates, the development of colistin resistance (CstR) has been shown to enhance biofilm formation, but reduce their invasiveness [27, 28]. Although these reports present evidence of likely interactions between virulence and the antimicrobial susceptibility profile of A. baumannii, the association of the expression of specific virulence genes with antimicrobial susceptibility has not been thoroughly investigated. Likewise, data comparing the expression of virulence genes among colistin-resistant A. baumannii (Cst-R-AB) strains is scarce.
In the present study, we aim to examine whether resistance to colistin in A. baumannii isolates is associated with the expression of specific virulence genes. We also analyze the genotypes of the A. baumannii isolates, their biofilm formation ability, and antimicrobial susceptibility profiles. By revealing the interaction of virulence genes with CstR among A. baumannii isolates, we attempt to explain, at least in part, the apparent discrepancies between in vitro antimicrobial resistance data, and in vivo clinical outcomes in A. baumannii infected burn patients. Our findings may ultimately help devise strategies toward effective treatments of A. baumannii burn wound infections, as well as make prudent decisions regarding other infection control measures to thwart nosocomial outbreaks of MDR A. baumannii.
Specimens, bacterial isolates and cultures
One hundred non-repetitive clinical specimens were collected from the burn unit of the Burn Medical Center Complex in Tehran, Iran; during April–December 2014. After admission, written informed consent forms were obtained from either the patients, or their authorized representatives. This study was approved by Tehran University of Medical Sciences (TUMS) Research Ethics Committee (Application No. TREC-89. 01-30-10430). All 100 burn wound isolates were initially identified as A. baumannii using the API20NE system (bioMérieux, Marcy-l’Etoile, France), and later confirmed using the gyrB multiplex PCR, as described previously . Mueller–Hinton Broth (MHB) and brain heart infusion (BHI) agar plates (both from Merck, Darmstadt, Germany) were used to culture the bacterial isolates under aerobic conditions for 24 h at 37 °C. In total, five Cst-R-AB isolates were analyzed including, a single Cst-R-AB isolates from this study, and four additional Cst-R-AB isolates from previous cohort studies.
Antimicrobial susceptibility testing
To assess antimicrobial susceptibility patterns of A. baumannii clinical isolates, we carried out the disk agar diffusion (DAD) method according to the Clinical and Laboratory Standards Institute (CLSI) procedures  and breakpoint interpretations, using antimicrobial disks containing the 17 antimicrobial agents: ampicillin–sulbactam, cefepime, ceftazidime, ciprofloxacin, colistin, gentamicin, imipenem, levofloxacin, meropenem, minocycline, piperacillin, piperacillin-tazobactam, rifampicin, tetracycline, tigecycline, tobramycin, and trimethoprim-sulfamethoxazole (Mast Diagnostics, Bootle, UK). The CLSI guideline for Broth microdilution test for minimum inhibitory concentration (MIC) was used to assess isolate susceptibility to colistin, rifampicin, and tigecycline. Colistin MICs were interpreted according to the CLSI breakpoints . For tigecycline susceptibility tests, the criteria of the European Committee on Antimicrobial Susceptibility Testing (EUCAST)  for Enterobacteriaceae were used, in which an MIC of < 1 µg/mL was defined as susceptible and > 2 µg/mL was considered resistant . Rifampicin susceptibility was interpreted according to CLSI criteria using breakpoint values suggested for Staphylococcus aureus, in which susceptible and resistant were defined as ≤ 1 and ≥ 4 µg/mL, respectively . A. baumannii isolates were defined as multi-, extensive- and pan-drug resistant (MDR, XDR, and PDR) according to the definitions described previously . Antimicrobial agents were categorized into three groups; namely, Group A that were agents deemed appropriate for primary testing panel; Group B agents, for which primary testing may be warranted, but they may be reported selectively (such as when the organism is resistant to agents of the same class, as in group A). Also, Group O (other) that comprised agents that had a clinical indication for A. baumannii, but were not candidates for primary testing panel in the USA . Colistin resistance (CstR) among A. baumannii isolates was confirmed by the genetic analysis of the pmr operon showing mutations in pmrA/pmrB signal transduction system, which confers CstR by modifying lipid A phosphoethanolamine moiety of the lipopolysaccharide (LPS) layer in the outer membrane.
International clones, MLVA genotype analysis
International clone (IC) types of the isolates were determined based on the presence of ompA, csuE, and blaOXA-51-like allele amplicons using two complementary multiplex PCR assays, as previously described . Table 1 shows the primers used in the multiplex PCR assays, which selectively amplified the outer membrane protein A (ompA), chaperone-subunit usher E (csuE), and blaOXA-51-like intrinsic carbapenemase genes of the A. baumannii isolates. Isolates not assigned as either IC type I, II, or type III were reported as the variant clonal type (IC-V).
To examine the genotypic diversity of A. baumannii clinical isolates, multi-loci variable-number tandem repeat analysis (MLVA) was carried out, as previously described . Briefly, DNA from the A. baumannii isolates was extracted using a GeneJET DNA purification kit and its concentration was assessed using a NanoDrop 1000 spectrophotometer (Thermo Scientific, USA). Genomic DNA (25 ng/µL) was assayed by PCR using a 30 μL final volume containing 5 ng of DNA, 10× reaction buffer (SinaClon BioScience, Iran), 1.5 mM MgCl [1, 2] U of Taq DNA polymerase, 200 μM each deoxynucleotides (dNTP), and 0.3 μM of each primer (all provided by SinaClon BioScience, Iran). Amplification cycles included an initial 94 °C denaturation for 5 min, followed by 35 cycles of 94 °C for 30, 30 s of annealing at 50 or 55 °C depend on melting temperature of the primers, elongation at 72 °C for 30 s, and a final elongation at 72 °C for 10 min. Then, 5 μL of each PCR amplicon solution were analyzed on a 25 cm 3% agarose gel (SinaClon BioScience, Iran) by electrophoresis (3 h at 6 V/cm) in 1× Tris–borate-EDTA (TBE) buffer. Size markers included a 100-bp (SinaClon Bioscience), and 50-bp DNA standards (Thermo Scientific) for L- and S-variable number tandem repeats (VNTRs), respectively. The amplicon bands were visualized with UV illumination after ethidium bromide (0.5 μg/mL) staining. Band sizes were determined using GeneTools v.3.08 automatic image analysis software (Syngene, Cambridge, United Kingdom). For cluster analysis, allele strings were analyzed by BioNumerics software v.7.0 as character values (Applied Maths, Sint-Martens-Latem, Belgium). Clusters were defined using cut-off values of 95% similarity, and MLVA type was determined using a 100% similarity cutoff, as previously described .
Semi-quantitative biofilm formation assay
The semi-quantitative assessment of biofilm formation was performed in triplicate using crystal violet staining, as previously described . Briefly, A. baumannii isolates were cultured on BHI agar overnight, and an isolated colony was suspended in LB Broth (Himedia, India) for 4 h. Bacterial suspensions (200 µL) in logarithmic growth phase were adjusted to 0.5 McFarland’s standard [1.5 × 108 colony forming units (CFU)/mL] and incubated at 35 °C in flat-bottomed 96-well plates. After 48 h, media was removed, and wells were washed three times with dH2O. Then, 200 µL of 1% crystal violet dye was added and plates were incubated at room temperature (RT) for 20 min. The dye was then decanted, wells were washed twice with dH2O and dried. Finally, 200 µL of 95% ethanol (200 µL/well) was added, and the optical absorbance (A) was measured at 570 nm (Thermo Scientific GmbH, Driesch, Germany). Reference strain (ATCC DH5α) was used as negative control (NC) to assign scores for biofilm formation, according to following formula: negative isolate (N) = AI ≤ ANC; weakly positive isolate (W) = ANC < AI ≤ 2ANC; moderately positive isolate (M) = 2ANC < AI ≤ 4ANC; and strongly positive isolate (S) = AI > 4ANC.
Virulence genes detection
To detect virulence factor genes dnaK, recA, lpsB and blsA, a series of PCR assays were utilized using primers (Table 1; Primer 3 software v.4.0). Briefly, 1 μL of genomic DNA (25 ng) from Cst-R-AB strains was assessed by PCR using 12.5 μL of 2× PCR Master Mix, and 1 μL (2.5 pM) of each primer (both from SinaClon BioScience, Iran). Amplifications were performed using a Mastercycler Personal (Eppendorf®, Germany), with initial denaturation for 5 min at 94 °C, and 35 cycles at 94 °C for 30 s, annealing for 30 s at 52 °C (Table 1), followed by elongation at 72 °C for 30 s, and a final step at 72 °C for 5 min. The presence of gene specific amplicons was verified visually by 2% agarose gel electrophoresis using ethidium bromide staining.
Animal study design
All animal experiments were carried out in accordance to the protocols approved by the Animal Ethics Committee of Tehran University of Medical Sciences (Application No. TUMS-AEC-89-0130-10430). Male C57BL/6 mice (6–8 week old; 18–23 g; Pasteur Institute, Tehran, Iran) were housed one mouse per cage, under sanitary conditions at 22–25 °C and at 12 h light/dark cycles, with access to sanitized pellet food and water. Mice were acclimated to room conditions for 1 week prior to each experiment. To increase the accuracy of microbiological assessments, the cages were disinfected with povidone iodine solution 10%, and the bedding materials were autoclaved, and replaced every day.
Based on mortality rates of burn wound infection from previous studies, the effect size (No. of mice in experimental groups) was estimated as five mice per each isolate, using power analysis with power arbitrarily set at 80%. Initially, two groups (five mice per group) were infected with the only Cst-R-AB (MT10) isolated in this study. One group was infected with colistin-treated Cst-R-AB isolate, and the other with Cst-R-AB isolate with no colistin treatment. Mock-infected burned mice (C-Burn) served as controls, and received physiological saline instead of A. baumannii infection (third group).
Similarly, groups of mice were later infected with additional four Cst-R-AB isolates (MIC 32–256 µg/mL), recovered from burn wounds in a previous study . In order to expose the Cst-R-AB isolates to colistin, and evaluate the putative colistin-induced changes in virulence gene expression, the group of colistin-treated Cst-R-AB infected mice were administered with colistin (20 µg/g weight); which yielded sub-MIC serum levels of approximately 25 µg/µL . For analysis of all five (i.e. 1 + 4) Cst-R-AB isolates, all in vivo data were combined and presented as Mean + SD.
Murine burn wound infection and bacterial quantification
Burn wounds were generated, as described previously . Briefly, mice were anesthetized by an intraperitoneal (i.p.) injection of a ketamine–xylazine cocktail, after which the dorsal/posterior surface was shaved, and each mouse was then injected with 0.5 mL (i.p.) of sterile saline to prevent dehydration. A brass 10 mm × 10 mm block preheated to 95 °C was applied to the shaved area for a 10 s period. This procedure consistently generated a nonlethal, full-thickness, third-degree burn. Burn wound infections were established as previously described . Briefly, 5 min after inducing the burn wound, using a sterile inoculating loop, a 50 µL suspension of each Cst-R-AB isolate containing 1.1–2.1 × 105 CFU/mL was smeared onto the burn wound site. After 72 h, bacterial concentrations in wound sites of mice infected with un-treated Cst-R-AB were compared to mice infected with Cst-treated Cst-R-AB isolate (treated with sub-optimal i.e. < MIC). Viable bacteria quantitated by colony enumeration of bacteria harvests from the wound bed, and reported as CFU/mL by using biopsy punch method as previously described . All mice were euthanized via an overdose of ketamine (250 mg/kg) and xylazine (25 mg/kg) prior to bacterial harvest from the wound site.
RNA extraction and cDNA synthesis and in vivo virulence gene expression
To measure in vivo virulence genes expression within the wound bed skin, tissue specimens were collected from the infected tissue from each sacrificed mouse were processed separately. Briefly, wound bed skin tissue specimens were aseptically excised (2 mm deep), placed in PBS, and frozen in liquid N2 for immediate RNA extraction.
To compare the relative quantities of gene-specific mRNA of Cst-R-AB isolates under in vivo vs. in vitro conditions RNA extraction was carried out. For the “in vivo samples” mRNA was extracted from Cst-R-AB isolates in their respective burn wound sites, 72 h after inoculation. The in vitro mRNA samples were, however, prepared using Cst-R-AB isolates grown in culture media, as described above. In all comparative analyses, 16S rRNA gene mRNA served as control to normalize mRNA quantities among the samples. Total bacterial RNA of planktonic mid-log phase cultures of Cst-R-AB isolates were extracted using the Total RNA Extraction Kit (iNtRON Biotech Inc., Seoul, South Korea). To extract bacterial RNA from tissue, wound specimens were ground and homogenized in liquid nitrogen, and the cellular lysates were run through All Prep DNA/RNA columns. Total RNA extraction was performed using RNeasy Plus Mini Kit (QIAGEN, Hilden, Germany) following the manufacturer’s recommendations. Bacterial mRNA enrichment was achieved by step-wise treatment of total RNA samples to eliminate mouse mRNA [NEBNext Poly(A) mRNA Magnetic Isolation Module, New England Biolabs, USA], followed by the removal of mouse and bacterial rRNA (Ribo-Zero Gold rRNA Removal Kit, Epicentre Biotechnologies, USA). RNA purity was determined by 260/280 and 260/230 nm absorbance ratio using a NanoDrop 1000 spectrophotometer (Thermo Scientific, USA), and size-verified by 2% agarose gel electrophoresis prior to RT-PCR and qRT-PCR assays. Before generating cDNA, bacterial mRNA samples were treated with DNase I (Thermo Scientific, USA) and RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific, USA). Samples of 1 μg total RNA were used to generate cDNA libraries, for the reverse transcription (RT)-PCR, and quantitative real time RT-PCR (qRT-PCR) assays.
Reverse transcriptase (RT-PCR) and quantitative (qRT)-PCR assays
RT-PCR assay was used to assess the expression of dnaK, blsA, recA, and lpsB genes. Briefly, first-strand cDNA served as a template for the amplification of virulence genes, and qRT-PCR was used to measure the relative changes in mRNA levels as indicators of virulence gene expression. Samples of cDNA from Cst-R-AB isolates were analyzed by qRT-PCR assay using the primer sets shown in Table 1. The Line-Gene K Real Time PCR System (BIOER Technology Co., China) was used to perform the qRT-PCR according to the minimum information for publication of quantitative real-time PCR experiments (MIQE) guidelines  under the following conditions: 95 °C for 15 min, followed by 45 cycles of 95 °C for 15 s, annealing for 30 s at 60, and 72 °C for 30 s, using 10 µL of SYBR® Premix Ex Taq™ II (Tli RNaseH Plus) (TaKaRa Bio Inc., Japan), 2 µL of gene-specific forward and reverse primers (2.5 pmol), cDNA template (50 ng), and 7 µL of ddH2O (total 20 µL). The quantity of gene-specific amplicons was determined by qPCR relative to an internal standard (16S rRNA) used as a calibrator. Target amplicons were checked for size and specificity by agarose gel electrophoresis and melting curve analysis, respectively. The changes in the expression level of each virulence gene was calculated by the 2−∆∆Ct method using the Relative Expression Software Tool (REST) 2009 software (version 2.0.13; Qiagen, Valencia, CA, USA) . Difference of > 4-fold relative to basal gene expression levels was considered significant, and referred to as n-fold difference.
Data were analyzed using the Student’s t-test, and the Chi square test using the SPSS software package (version 22). Results were considered significant if P < 0.05.
Antimicrobial susceptibility profiles
As shown in Fig. 1, the majority (95%) of A. baumannii isolates had a MDR (31%) or XDR (64%) profile of antimicrobial resistance, while the frequency of non-MDR, or PDR isolates was 4%, and 1%, respectively. While most isolates (96%) were resistant to at least three classes of antimicrobial agents, nearly 90% were resistant to at least 10 tested antimicrobial agents including; cephalosporins, trimethoprim/sulfamethoxazole, tetracycline, piperacillin/tazobactam, and fluoroquinolones. Figure 1 also shows the susceptibility profile of A. baumannii isolates against seven antimicrobial agents representing the A, B, and O groups of antimicrobial agents . The highest resistance rates were among group A (45–96%), followed by group B (19–97%), and group O antimicrobial agents (1%). The lowest rates of resistance were against colistin (1%), tigecycline (19%), and minocycline (42%). Moreover, resistance rates to ciprofloxacin, tobramycin, and imipenem were 89, 64 and, 60% respectively. Surprisingly, the rate of resistance to tobramycin was higher among MDR isolates (82%) than the XDR (87%) isolates (P > 0.05).
MLVA genotyping and international clone analysis
International clone (IC) type analysis showed that the majority of A. baumannii isolates were IC-V variants (43%), followed by IC-II (37%), IC-I (19%), and IC-III types (1%). Also, MLVA genotype analysis (95% allelic similarity) revealed eight different strain clusters, including strain clusters A, B, and C in one group, as well as clusters D through H in another group (Fig. 1).
Overall, we identified 32 MLVA types (MT), which consisted of 1–24 members each (at 100% similarity cutoff). While most MTs had a single member, the three main MTs were MT29 (n = 24), MT13 (n = 17), and MT14 (n = 10), which together comprised over 50% of all isolates (Fig. 1). The majority of MT29 (18/24) and MT14 (8/10) members were IC-V variants, which were predominantly XDR isolates. Specifically, 83% (15/18) of MT29 IC-V variants, and 50% (4/8) of MT14 variants showed XDR profiles. The small MT12 cluster (n = 4) was remarkable since all four were IC-V and XDR isolates. Approximately 75% of the MT29 isolates and 40% of the MT14 isolates were categorized as IC-V. Among the MT29, MT13, and MT14 isolates, XDR isolates predominated and showed similar frequencies. Most MT13 isolates (76%) were IC-II types with nearly 70% of these being XDR strains (Fig. 1).
Association of MDR profile with IC types and biofilm formation index
We further examined the association of the MDR profile of antimicrobial resistance of each isolate with the ability to form biofilm as a virulence determinant, as well as the IC type of each isolate (Fig. 2a, b). Our analysis focused on comparing the biofilm formation index of each isolate with the isolate MDR and XDR profile, since 95% of all isolates had either an MDR or XDR profile. Figure 2a shows that the majority of XDR isolates were IC variants, whereas MDR isolates were predominantly IC-II type. The frequency of IC-II isolates with MDR profile (55%) was almost two folds higher than that of those with XDR profiles (28%; P = 0.031). In contrast, the frequency of IC-V isolates with MDR profiles was nearly half the frequency of those with XDR profiles (23% vs. 53%; P = 0.021). The frequencies of MDR vs. XDR profiles among IC-I and IC-III isolates was similar (Fig. 2a).
Comparison of the antimicrobial resistance profile of the isolates’ with their biofilm strength index revealed that 77% of the W- and M-biofilm forming isolates had either an MDR, or XDR profile. While, the S-biofilm formers comprised only 23% of MDR isolates, the single PDR isolate was a W-biofilm former. All MT12 members (n = 4) were IC-V type and S-biofilm formers that showed an XDR antimicrobial susceptibility profile. Among XDR isolates, the frequency of S- and M-biofilm formers were 48 and 43%, respectively; while among the MDR isolates, W- and M-biofilm formers were isolated at an identical 39% rate (Fig. 2b). In contrast, the frequency of W-biofilm formers among MDR isolates was 39%, whereas only 12% of XDR isolates were W-biofilm formers (P = 0.009). Similarly, the frequency of S-biofilm forming isolates with the XDR, or MDR antimicrobial susceptibility profile was 47 and 23%, respectively (P = 0.004). Finally, all four non-MDR isolates were either W- or S-biofilm formers (50% frequency of isolation); i.e. they included no M-biofilm formers.
Association of IC type with biofilm formation
Additionally, we assessed whether the strength of biofilm formation was associated with the IC type of A. baumannii isolates. Figure 3 shows the frequency of W-, M- and S-biofilm formers among IC-I, IC-II, and IC-V isolates, which comprised 99% of all isolates. Among IC-V isolates, the frequency of S-biofilm formers (63%) was six folds higher than the W-biofilm formers (12%; P = 0.001); while, among IC-I and IC-II isolates, the frequencies of W-, M-, and S-biofilm formers were comparable, showing only 4–7% variation: i.e. 26% vs. 19% for W-biofilm formers (P = 0.382), and 26% vs. 30% for S-biofilm formers (P = 0.521), and finally 47% vs. 51% for M-biofilm formers (P = 0.437), in that order. As shown in Fig. 3, M-biofilm-forming isolates were almost twice as often seen among the IC-I and IC-II isolates compared to the IC-V isolates [47% vs. 26% for IC-I (P = 0.036), and 51% vs. 26% for IC-II (P = 0.028), respectively.
Quantification of CFU concentration in wound bed tissues
There was no significant difference between the log CFU concentrations in the wound sites infected with un-treated Cst-R-AB isolate (7.78 ± 6.0), versus the log CFU concentrations (6.99 ± 5.82; P > 0.05) in wound beds that were infected with Cst-R-AB, which was treated with sub-MIC colistin. Thus, Cst-R-AB isolate concentrations (i.e. growth) were unaffected by the Cst-treatment of wound bed infection site.
Virulence genes detection and expression in vitro and in vivo
Virulence factor genes dnaK, recA, lpsB and blsA, were detected in the Cst-R-AB strain by PCR. In order to examine whether the expression of virulence genes dnaK, blsA, recA, and lpsB is associated with Cst resistance among isolates, we compared the relative gene-specific mRNA quantities of Cst-R-AB isolates under in vivo, vs. in vitro conditions (Fig. 4). The final data analysis included five Cst-R-AB isolates; including a single Cst-R-AB isolate from the present study (Additional file 1: Figure S1), and four additional Cst-R-AB isolates from previous studies.
Overall, as shown in Fig. 4, the expression levels of dnaK, recA, and lpsB genes in Cst-R-AB isolates were markedly different between in vivo vs. in vitro conditions. The relative change over basal expression levels of the four virulence genes, in vivo vs. in vitro, among Cst-R-AB isolates was reported as the Mean ± SD (Fig. 4).
Among the Cst-R-AB isolates, the in vivo expression of dnaK and lpsB was 23 and 14 folds higher than basal; respectively (Fig. 4); while recA mRNA levels showed a slight 3-fold increase (P < 0.05). Conversely, Cst-R-AB in vitro cultures revealed no significant increase in the expression of dnaK, blsA, and lpsB, with only a modest 4-fold rise in recA expression. Interestingly, under either condition, blsA gene mRNA levels of the Cst-R-AB isolates was up-regulated by 3–4-fold.
The emergence of multidrug-resistant A. baumannii (MDR-AB), as a major nosocomial pathogen among patients with burn wound infections, has drawn clinical attention to an old antimicrobial agent that is active against MDR-AB, namely colistin (Cst) [43, 44]. Unfortunately, reports from several regions of the world, including Iran, indicate that the widespread use of Cst, as a last resort to control MDR-AB infections, has led to a worrisome growing trend of Cst-resistance (CstR) among MDR-AB strains [8, 45, 46]. Recurrently, the eradication of MDR-AB burn wound infections prove very challenging because, not only do clinical strains often develop antimicrobial resistance, but these strains also produce virulence factors (e.g. biofilm) that transform these MDR-AB into even more formidable wound pathogens [10, 12].
While, some studies have associated the presence of specific virulence genes with an increase in the pathogenic potential of clinical isolates [47, 48], the association between the virulence gene expression of A. baumannii and the susceptibility to specific antimicrobial agents has not been assessed. In fact, despite several reports regarding the increase in the prevalence of Cst-resistance among MDR-AB isolates worldwide [8,9,10,11,12,13,14, 45], few studies have focused on the possible link between CstR, and the expression of virulence genes that can potentially affect therapy outcome. We have used a murine burn infection model to demonstrate that mRNA levels of virulence genes dnaK, blsA, recA, and lpsB vary in Cst-R-AB isolates depending on whether they are cultured in vitro or in vivo. Moreover, to examine the genotypic diversity among A. baumannii isolates, we have characterized their MLVA genotypes, IC lineage. We further evaluated any associations between the isolates’ IC lineage, susceptibility patterns, and their biofilm formation profiles.
Here, we also present evidence suggesting that the XDR or MDR profile of A. baumannii is associated with the biofilm strength. As shown in Fig. 2b, the frequency of W-biofilm formers among MDR isolates was about three times higher than among XDR isolates, which suggests that the strength of biofilm formation is related to the extended spectrum of antimicrobial resistance among isolates. Conversely, W- and M-biofilm formers were equally common (39%) among MDR isolates. That the XDR isolates had more M- and S-biofilm formers than MDR isolates suggests that the antimicrobial susceptibility profile is associated with biofilm strength among these isolates. However, since the rates of M-biofilm formation were similar among XDR and MDR isolates, the M-biofilm formation phenotype does not seem to follow this association pattern.
Interestingly, the frequency of S- and M-biofilm formers among XDR isolates were similar (48% vs. 43%, respectively); whereas among the MDR isolates, the W- and M-biofilm formers were isolated at an identical 39% rate (Fig. 2b). Conversely, the frequency of W-biofilm formers among MDR isolates were about three times higher than XDR isolates (P = 0.009), which suggests that strength of biofilm formation is associated with the isolate antimicrobial profile. Likewise, isolates with XDR antimicrobial susceptibility profiles were more than twice as likely to be S-biofilm formers, as compared to the frequency of S-formers among MDR isolates (P = 0.004). This suggests that stronger biofilm formation may be associated with broader antimicrobial resistance. The remarkable finding that all MT12 members were S-biofilm formers, and shared the same IC-V type and antimicrobial susceptibility profile is suggestive of a common source of isolation.
Interestingly, while IC-V variants had the highest frequency of S-biofilm formers, they included lowest frequency of W-biofilm formers. In addition, our data revealed that most of S-biofilm forming IC-V isolates also had XDR profiles, which is contrary to the reports that indicate biofilm formation is a common attribute of clinical A. baumannii isolates, regardless of their IC minor or major clonal variation [49,50,51,52]. This apparent discrepancy might be due to the unique features of our local IC-V isolates that have neither been isolated, nor characterized elsewhere. Consequently, the conclusion of studies that have reported no correlation between A. baumannii biofilm formation and IC lineage may be limited to the IC-I, -II, and -III categories, which apparently do not encompass a global representation of all A. baumannii isolates.
Here, we present evidence suggesting that the development of CstR may also lead to changes in A. baumannii virulence determinant genes expression, which in turn might enhance the ability to establish, and persist in burn wound infections. In order to further strengthen the scope of our initial findings, which were based on a single Cst-R-AB isolate, we expanded our in vivo experimental data by using four additional Cst-R-AB isolates in a similar murine model (Fig. 4). Since CstR has been shown to cause increased outer membrane permeability, higher division rate, and structural changes in the cell wall of A. baumannii isolates, it could be proposed that during prolonged infections, persistent CstR can lead to reduced fitness, and lower A. baumannii virulence over time. However, our findings suggest that the opposite may occur, where enhanced A. baumannii virulence can be of greater benefit to some highly resistant local variants of MDR-AB. While the present study focused on the expression profile of four virulence genes among five Cst-R-AB isolates, the variation in the levels of virulence gene mRNA in vitro vs. in vivo conditions among Cst-R-AB isolates could be due to a myriad of additional factors, such as differences between the cellular environment of the wound site as compared to the cell-free culture conditions. Moreover, the role of tissue specific factors that dictate gene regulation, non-specific immune responses, remains to be investigated.
To our knowledge, this is the first report that attempts to establish an association between CstR among A. baumannii isolates and changes in the virulence genes expression. Our findings show the merits of a comprehensive investigation using a larger number of Cst-R-AB isolates from various regions of the world. Future detailed studies will shed light on the causal relationship between expression of a wider array of A. baumannii virulence genes, and genes that confer resistance to colistin, as well as other antimicrobial agents.
Finally, our findings highlight the intricate relationship between resistance to colistin, and the expression of A. baumannii virulence genes. It is commonly supposed that an increase in antimicrobial resistance leads, directly or indirectly, to reduced virulence of A. baumannii isolates. However, our data suggest that the opposite may occur, and it provides evidence that thorough investigations are warranted to directly measure the correlation between CstR, reduced virulence, and “fitness” of A. baumannii isolates. Virulence gene expression variations among Cst-R-AB isolates might gain further clinical importance in health care settings, where isolates often transfer from patient to patient. Moreover, that the virulence of Cst-R-AB isolates may depend on their in vitro or in vivo conditions can help explain, at least in part, commonly encountered discrepancies between in vivo and in vitro laboratory test results. For instance, the recurrent failure of therapeutic antimicrobial agents in eradicating A. baumannii in patients’ burn wounds, despite the in vitro susceptibility of isolated pathogen to the same agents.
This study concludes that the development of resistance to colistin among A. baumannii may be associated with virulence gene expression, which fluctuates under in vivo and in vitro conditions. We also provide evidence to speculate that enhanced virulence may be of greater benefit to highly-resistant variants of MDR-AB, and may give rise to formidable pathogens that increase morbidity among burn patients, and threaten the public health systems worldwide. Our findings warrant detailed examination of the interactions between virulence and antimicrobial resistance toward efforts to control the spread of multidrug-resistant A. baumannii (MDR-AB) isolates, and also to reduce disease severity in burn patients with MDR-AB infection.
Acinetobacter baumannii or A. baumannii
multidrug-resistant A. baumannii
colistin-resistant A. baumannii
- recA :
gene encoding the RecA protein
- lpsB :
gene encoding the lipopolysaccharide core biosynthesis mannosyltransferase protein LpsB
- dnaK :
gene encoding the chaperone protein DnaK
- blsA :
gene encoding the BlsA (blue-light-sensing) photoreceptor protein
Gesellschaft mit beschränkter Haftung, which means “company with limited liability”
- CstR :
Tehran University of Medical Sciences
- gyrB :
gene encoding the GyrB subunit of DNA gyrase protein
polymerase chain reaction
disk agar diffusion
Clinical and Laboratory Standards Institute
minimum inhibitory concentration
European Committee on Antimicrobial Susceptibility Testing
Mock-infected burned mice that serve as controls
- ompA :
gene encoding outer membrane protein A
- bla OXA-51 :
the gene encoding the carbapenemase OXA-51
outer membrane protein A
chaperone-subunit usher E
multi-loci variable-number tandem repeat analysis
- MgCl2 :
minimum information for publication of quantitative real-time PCR experiments
- et al.:
Et alii meaning “and other people”
variable number tandem repeats
large variable number tandem repeats
small variable number tandem repeats
colony forming unit(s)
American Type Culture Collection
moderately positive isolate
strongly positive isolate
- AI :
- ANC :
absorbance of negative control
- ® :
- N2 :
svedberg unit(s) (e.g. 16S)
ribosomal ribonucleic acid
mitochondrial ribonucleic acid
United States of America
reverse transcription PCR
quantitative real-time RT-PCR
- DNase I:
double distilled water
Relative Expression Software Tool
- n-fold difference:
Relative fold difference as compared to the basal gene expression
- W-biofilm formers:
weak biofilm formers
- M-biofilm formers:
moderate biofilm formers
- S-biofilm formers:
strong biofilm formers
Girerd-Genessay I, Bénet T, Vanhems P. Multidrug-resistant bacterial outbreaks in burn units: a synthesis of the literature according to the ORION statement. J Burn Care Res. 2016;37(3):172–80.
Leseva M, Arguirova M, Nashev D, Zamfirova E, Hadzhyiski O. Nosocomial infections in burn patients: etiology, antimicrobial resistance, means to control. Ann Burns Fire Disasters. 2013;26(1):5–11.
Golberg A, Broelsch GF, Vecchio D, et al. Eradication of multidrug-resistant A. baumannii in burn wounds by antiseptic pulsed electric field. Technology (Singap World Sci). 2014;2(2):153–60.
Atilla A, Tomak L, Katrancı AO, Ceylan A, Kılıç SS. Mortality risk factors in burn care units considering the clinical significance of acinetobacter infections. Ulus Travma Acil Cerrahi Derg. 2015;21(1):34–8.
Mortensen BL, Skaar EP. Host-microbe interactions that shape the pathogenesis of Acinetobacter baumannii infection. Cell Microbiol. 2012;14(9):1336–44.
Longo F, Vuotto C, Donelli G. Biofilm formation in Acinetobacter baumannii. New Microbiol. 2014;37(2):119–27.
Weber BS, Harding CM, Feldman MF. Pathogenic acinetobacter: from the cell surface to infinity and beyond. J Bacteriol. 2015;198(6):880–7.
Bahador A, Taheri M, Pourakbari B, et al. Emergence of rifampicin, tigecycline, and colistin-resistant Acinetobacter baumannii in Iran; spreading of MDR strains of novel international clone variants. Microb Drug Resist. 2013;19(5):397–406.
Moradi J, Hashemi FB, Bahador A. Antibiotic resistance of Acinetobacter baumannii in Iran: a systemic review of the published literature. Osong Public Health Res Perspect. 2015;6(2):79–86.
Doi Y, Murray GL, Peleg AY. Acinetobacter baumannii: evolution of antimicrobial resistance-treatment options. Semin Respir Crit Care Med. 2015;36(1):85–98.
Inchai J, Liwsrisakun C, Theerakittikul T, Chaiwarith R, Khositsakulchai W, Pothirat C. Risk factors of multidrug-resistant, extensively drug-resistant and pandrug-resistant Acinetobacter baumannii ventilator-associated pneumonia in a Medical Intensive Care Unit of University Hospital in Thailand. J Infect Chemother. 2015;21(8):570–4.
Potron A, Poirel L, Nordmann P. Emerging broad-spectrum resistance in Pseudomonas aeruginosa and Acinetobacter baumannii: mechanisms and epidemiology. Int J Antimicrob Agents. 2015;45(6):568–85.
Gupta M, Lakhina K, Kamath A, et al. Colistin-resistant Acinetobacter baumannii ventilator-associated pneumonia in a tertiary care hospital: an evolving threat. J Hosp Infect. 2016;94(1):72–3.
Liu H, Wu YQ, Chen LP, et al. Biofilm-related genes: analyses in multi-antibiotic resistant Acinetobacter Baumannii isolates from Mainland China. Med Sci Monit. 2016;28(22):1801–7.
Zarrilli R. Acinetobacter baumannii virulence determinants involved in biofilm growth and adherence to host epithelial cells. Virulence. 2016;7(4):367–8.
Lee J, Page R, García-Contreras R, et al. Structure and function of the Escherichia coli protein YmgB: a protein critical for biofilm formation and acid-resistance. J Mol Biol. 2007;373(1):11–26.
Wroblewska MM, Sawicka-Grzelak A, Marchel H, Luczak M, Sivan A. Biofilm production by clinical strains of Acinetobacter baumannii isolated from patients hospitalized in two tertiary care hospitals. FEMS Immunol Med Microbiol. 2008;53(1):140–4.
Fiester SE, Actis LA. Stress responses in the opportunistic pathogen Acinetobacter baumannii. Future Microbiol. 2013;8(3):353–65.
Schröder W, Goerke C, Wolz C. Opposing effects of aminocoumarins and fluoroquinolones on the SOS response and adaptability in Staphylococcus aureus. J Antimicrob Chemother. 2013;68(3):529–38.
Aranda J, Bardina C, Beceiro A, et al. Acinetobacter baumannii RecA protein in repair of DNA damage, antimicrobial resistance, general stress response, and virulence. J Bacteriol. 2011;193(15):3740–7.
Beceiro A, Tomás M, Bou G. Antimicrobial resistance and virulence: a successful or deleterious association in the bacterial world? Clin Microbiol Rev. 2013;26(2):185–230.
Luke NR, Sauberan SL, Russo TA, et al. Identification and characterization of a glycosyltransferase involved in Acinetobacter baumannii lipopolysaccharide core biosynthesis. Infect Immun. 2010;78(5):2017–23.
McConnell MJ, Actis L, Pachón J. Acinetobacter baumannii: human infections, factors contributing to pathogenesis and animal models. FEMS Microbiol Rev. 2013;37(2):130–55.
Lean SS, Yeo CC, Suhaili Z, Thong KL. Comparative genomics of Two ST 195 carbapenem-resistant Acinetobacter baumannii with different susceptibility to polymyxin revealed underlying resistance mechanism. Front Microbiol. 2016;6:1445. https://doi.org/10.3389/fmicb.2015.01445.
Mussi MA, Gaddy JA, Cabruja M, et al. The opportunistic human pathogen Acinetobacter baumannii senses and responds to light. J Bacteriol. 2010;192(24):6336–45.
Golic A, Vaneechoutte M, Nemec A, Viale AM, Actis LA, Mussi MA. Staring at the cold sun: blue light regulation is distributed within the genus Acinetobacter. PLoS ONE. 2013;8(1):e55059. https://doi.org/10.1371/journal.pone.0055059.
Pournaras S, Poulou A, Dafopoulou K, et al. Growth retardation, reduced invasiveness, and impaired colistin-mediated cell death associated with colistin resistance development in Acinetobacter baumannii. Antimicrob Agents Chemother. 2014;58(2):828–32.
Dafopoulou K, Xavier BB, Hotterbeekx A, et al. Colistin-resistant Acinetobacter baumannii clinical strains with deficient biofilm formation. Antimicrob Agents Chemother. 2015;60(3):1892–5.
Higgins PG, Schneiders T, Hamprecht A, Seifert H. In vivo selection of a missense mutation in adeR and conversion of the novel blaOXA-164 gene into blaOXA-58 in carbapenem-resistant Acinetobacter baumannii isolates from a hospitalized patient. Antimicrob Agents Chemother. 2010;54(12):5021–7.
Clinical and Laboratory Standards Institute. Performance standards for antimicrobial susceptibility testing: twenty fourth informational supplement M100-S28. Wayne: CLSI; 2018.
The European Committee on Antimicrobial Susceptibility Testing. Breakpoint tables for interpretation of MICs and zone diameters Version 5.0. 2015. http://www.eucastorg/clinical_breakpoints. Accessed 01 Jan 2015.
Magiorakos AP, Srinivasan A, Carey RB, et al. Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin Microbiol Infect. 2012;18(3):268–81.
Turton JF, Gabriel SN, Valderrey C, Kaufmann ME, Pitt TL. Use of sequence-based typing and multiplex PCR to identify clonal lineages of outbreak strains of Acinetobacter baumannii. Clin Microbiol Infect. 2007;13(8):807–15.
Pourcel C, Minandri F, Hauck Y, et al. Identification of variable-number tandem-repeat (VNTR) sequences in Acinetobacter baumannii and interlaboratory validation of an optimized multiple-locus VNTR analysis typing scheme. J Clin Microbiol. 2011;49(2):539–48.
Nhu NT, Lan NP, Campbell JI, et al. Emergence of carbapenem-resistant Acinetobacter baumannii as the major cause of ventilator-associated pneumonia in intensive care unit patients at an infectious disease hospital in southern Vietnam. J Med Microbiol. 2014;63(Pt10):1386–94.
Zhang D, Xia J, Xu Y, et al. Biological features of biofilm-forming ability of Acinetobacter baumannii strains derived from 121 elderly patients with hospital-acquired pneumonia. Clin Exp Med. 2016;16(1):73–80.
Bahador A, Raoo An R, Farshadzadeh Z, et al. The prevalence of IS Aba 1 and IS Aba 4 in Acinetobacter baumannii species of different international clone lineages among patients with burning in Tehran, Iran. Jundishapur J Microbiol. 2015;8(7):e17167. https://doi.org/10.5812/jjm.17167v2.
Pachón-Ibáñez ME, Docobo-Pérez F, López-Rojas R, et al. Efficacy of rifampin and its combinations with imipenem, sulbactam, and colistin in experimental models of infection caused by imipenem-resistant Acinetobacter baumannii. Antimicrob Agents Chemother. 2010;54(3):1165–72.
Dai T, Tegos GP, Lu Z, et al. Photodynamic therapy for Acinetobacter baumannii burn infections in mice. Antimicrob Agents Chemother. 2009;53(9):3929–34.
Ha U, Jin S. Expression of the soxR gene of Pseudomonas aeruginosa is inducible during infection of burn wounds in mice and is required to cause efficient bacteremia. Infect Immun. 1999;67(10):5324–31.
Thompson MG, Black CC, Pavlicek RL, et al. Validation of a novel murine wound model of Acinetobacter baumannii infection. Antimicrob Agents Chemother. 2014;58(3):1332–42.
Bustin SA, Beaulieu JF, Huggett J, Jaggi R, Kibenge FSB, Olsvik PA, Penning LC, Toegel S. MIQE précis: practical implementation of minimum standard guidelines for fluorescence-based quantitative realtime PCR experiments. BMC Mol Biol. 2010;11:74.
Saun TJ, Rogers AD, Leis JA, Cartotto R. The use of intravenous and inhaled colistin therapy during a burn center outbreak of multidrug-resistant Acinetobacter baumannii. J Burn Care Res. 2017. https://doi.org/10.1097/BCR.
Fishbain J, Peleg AY. Treatment of acinetobacter infections. Clin Infect Dis. 2010;51(1):79–84.
Sepahvand S, Doudi M, Davarpanah MA, Bahador A, Ahmadi M. Analyzing pmrA and pmrB genes in Acinetobacter baumannii resistant to colistin in Shahid Rajai Shiraz, Iran Hospital by PCR: first report in Iran. Pak J Pharm Sci. 2016;29(4 Suppl):1401–6.
Garnacho-Montero J, Dimopoulos G, Poulakou G, et al. Task force on management and prevention of Acinetobacter baumannii infections in the ICU. Intensive Care Med. 2015;41(12):2057–75.
Eijkelkamp BA, Stroeher UH, Hassan KA, Paulsen IT, Brown MH. Comparative analysis of surface-exposed virulence factors of Acinetobacter baumannii. BMC Genomics. 2014;25(15):1020. https://doi.org/10.1186/1471-2164-15-1020.
Yoon EJ, Balloy V, Fiette L, Chignard M, Courvalin P, Grillot-Courvalin C. Contribution of the Ade resistance-nodulation-cell division-type efflux pumps to fitness and pathogenesis of Acinetobacter baumannii. MBio. 2016;7(3):pii: e00697–716. https://doi.org/10.1128/mbio.00697-16.
Lee NY, Chang TC, Wu CJ, et al. Clinical manifestations, antimicrobial therapy, and prognostic factors of monomicrobial Acinetobacter baumannii complex bacteremia. J Infect. 2010;61(3):219–27.
Azim A, Dwivedi M, Rao PB, et al. Epidemiology of bacterial colonization at intensive care unit admission with emphasis on extended-spectrum beta-lactamase- and metallo-beta-lactamase-producing Gram-negative bacteria–an Indian experience. J Med Microbiol. 2010;59(Pt 8):955–60.
Kang YS, Park W. Contribution of quorum-sensing system to hexadecane degradation and biofilm formation in Acinetobacter sp. strain DR1. J Appl Microbiol. 2010;109(5):1650–9.
Sanchez CJ, Kumar N, Lizcano A, et al. Streptococcus pneumoniae in biofilms are unable to cause invasive disease due to altered virulence determinant production. PLoS ONE. 2011;6(12):e28738. https://doi.org/10.1371/journal.pone.0028738.
FBH and AB analyzed and processed the data and drafted the manuscript. ZF and MM did the PCR and MLVA assays. MM did the specimen cultures, animal studies and biofilm assays. MP performed the antimicrobial assays. RR helped with the RNA extraction and virulence gene expression assays, BP helped with clinical coordination and obtaining the patient specimens from the hospitals. All authors read and approved of the final manuscript.
We would like to thank all the medical staff, who assisted with collection of the specimens used in this study. We would also like to acknowledge Ms. Jocelyn C. Gile for her invaluable editorial assistance.
The authors declare that they have no competing interests.
All documents and additional data are available from the corresponding author upon reasonable request.
Ethics approval and consent to participate
This study was approved by Tehran University of Medical Sciences (TUMS) Research Ethics Committee for Human Research (Application No. TRECHR-89. 01-30-10430). Written informed consent forms were obtained from either the patients, or their authorized representatives. All animal experiments were carried out in accordance to the protocols approved by the Animal Ethics Committee of Tehran University of Medical Sciences.
This research was supported by Tehran University of Medical Sciences & Health Services Grant No. 89. 01-30-10430.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1.
Comparison of relative change in the mRNA levels of virulence determinant genes dnaK, blsA, recA, and lpsB in cultures of a single Cst-R A. baumannii isolate, under in vitro and in vivo conditions. Bar indicates Mean + SD; five mice in each group. #P < 0.01; *P < 0.05; ***P < 0.001.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Bahador, A., Farshadzadeh, Z., Raoofian, R. et al. Association of virulence gene expression with colistin-resistance in Acinetobacter baumannii: analysis of genotype, antimicrobial susceptibility, and biofilm formation. Ann Clin Microbiol Antimicrob 17, 24 (2018). https://doi.org/10.1186/s12941-018-0277-6
- Acinetobacter baumannii
- Colistin resistance
- Virulence gene expression
- Murine burn model