Phenotypic and genomic comparison of dominant and nondominant sequence-type of Acinetobacter baumannii isolated in China

1 Introduction

A. baumannii is a Gram-negative bacterium that causes a number of clinically important, life-threatening infections including ventilator-associated pneumonia, bloodstream infections, and intracranial infections (Joly-Guillou, 2005; Maragakis and Perl, 2008; Freire et al., 2016). Furthermore, because of its inherent antibiotic resistance and the ease with which it acquires resistance elements from elsewhere to develop multidrug resistance, carbapenem-resistant A. baumannii has become a common clinical pathogen and is listed by WHO as a pathogenic bacterium requiring priority development of new antibacterial drugs (World Health Organization; Lee et al., 2017).

ST is often used to investigate the prevalence, transmission, and outbreaks of A. baumannii clones in different regions or hospitals, and thus to assist in formulating appropriate hospital infection control measures (Schultz et al., 2016). Several studies, in different regions of China, have shown that ST191, ST195, and ST208 (Oxford scheme) are the most prevalent A. baumannii ST types isolated (Deng et al., 2014; Xiao et al., 2016; Ning et al., 2017; Zhou et al., 2018). ST208 and ST191 accounted for 58.7% and 10.9%, respectively, of strains isolated in Shanghai, whereas 41.7% and 13.1% of the strains isolated in the south of China were ST195 and ST208, respectively. ST191 was also shown to be the most prevalent strain in Korea (Son et al., 2020). Consequently,we define ST191 ST195 and ST208 these top three dominant ST types as the DST, and the other ST types including ST462, ST547 and STn(new types of numbers not yet assigned) as non-dominant sequence type(NST).

Studies have shown that the outer membrane protein OmpA, phospholipase, capsule capsular polysaccharide(CPS), iron acquisition system, efflux pump, Csu chaperone usher-type pilus, and secretion system are all important virulence factors of A. baumannii (Harding et al., 2018). Moreover, the ability to survive in unfavorable environments has been suggested as an important virulence strategy. Indeed, it has been reported that specific ST strains are highly virulent and correlate with poor clinical prognosis in infected patients (Yoon et al., 2019).

Although desiccation resistance has been reported to contribute to the spread and persistence of type-specific A. baumannii in the hospital setting (Giannouli et al., 2013). The factors that facilitate the generation of DST, and the phenotypic and genotypic differences between DST and NST strains of A. baumannii, have not been fully elucidated. We combined biological experiments with genomic and transcriptomic analysis, therefore, to comprehensively elucidate the biological properties, especially those relating to virulence and resistance, that differentiate DST and NST strains.

2 Materials and methods

2.1 Isolates, culture conditions, and susceptibility tests

All A. baumannii strains were isolated from a tertiary, first-class teaching hospital in East China, and the strain number, ST type, etc. are shown in Table S1. 10 strains of ST191, 12 strains of ST195, 11 strains of ST208, 10 strains of ST462, and 9 strains of ST547, total 52 strains were used for each biological experiments. Two strains from each ST type above, plus two STn strains, a total of 12 strains were selected for subsequent genomic and transcriptomic analysis based on their ability to resist complement killing. The antimicrobial susceptibility of all A. baumannii isolates was determined using agar dilution according to the Clinical and Laboratory Standards Institute (CLSI) guidelines from 2020. Resistance breakpoints determined by the European Committee on Antimicrobial Susceptibility Testing (EUCAST) were used for polymyxin B.

2.2 Growth curves

Growth curves were performed as described previously with slight modifications (Hall et al., 2014) Briefly, A. baumannii in the log-phase of growth were adjusted to 0.5 absorbance units (OD600) and diluted 50-fold with LB medium. Aliquots (200 μl) were added to a 96-well plate, placed in a microplate reader at 37°C, and shaken for 5 s every 30 min, and the OD600 was measured for 24 h.

2.3 Desiccation resistance assays

Desiccation resistance assays were completed as previously described (Boll et al., 2015) Log-phase A. baumannii were harvested and washed twice with an equal volume of LB medium. Each sample was adjusted to 1×108 CFU/ml and then serially diluted and plated to determine the input CFU. Bacterial suspension (10 μl) was spotted onto 96-well polystyrene plates and desiccated at 30°C at a humidity of 40%. PBS (200 μl) was added to each well to resuspend the bacteria after 48 h and 96 h, respectively, and the resuspended bacteria were serially diluted and plated to calculate the output CFU. The percent of survival was defined as the ratio of output to input CFUs.

2.4 Oxidative killing assay

A. baumannii in the log-phase of growth were adjusted to OD₆₀₀ = 0.5 and spread evenly on an MH agar plate. A sterile filter paper disk with a diameter of 6 mm was placed in the middle of the inoculated agar, and 10 μl 20% (v/v) hydrogen peroxide was pipetted onto the filter paper. The plates were incubated at 37°C for 18 h, after which the growth inhibition zone around the disk was measured using Vernier calipers.

2.5 Antiserum complement killing

Blood collected from healthy volunteers was centrifuged and filtered with a 0.22 μm pore size syringe filter to obtain sterilized serum. Half of the serum was inactivated by heating in a 56°C water bath for 30 min. A. baumannii in the log-phase of growth were adjusted to 2 × 10⁶ CFU/ml and mixed with normal or inactivated serum at a ratio of 1:9, respectively. After incubation for 1 h at 37°C, samples were serially diluted and spread on MH agar plates. After overnight incubation at 37°C, bacteria colonies were counted, and the bacterial survival rate was calculated by the following formula: Bacterial survival rate = (number of colonies in normal serum/number of colonies in inactivated serum) × 100%.

2.6 Biofilm formation

A. baumannii in the log-phase of growth were adjusted to 0.5 absorbance units (OD₆₀₀), before 200 μl aliquots were pipetted into the wells of a 96-well plate. Bacteria were cultivated in an incubator for 24 h at 37°C and 5% CO₂, to allow biofilm formation, before nonattached bacteria were removed from the wells by aspiration. Biofilms were gently washed with ddH₂O, after which the 96-well plate was inverted on absorbent paper for 15 min at room temperature. Biofilm cells were fixed by adding 4% (v/v) paraformaldehyde to each well for 20 min and were then stained with 1% (w/v) crystal violet for 15 min. Biofilms were washed thoroughly with ddH₂O, before the bound crystal violet was eluted with 200 μl 95% ethanol, and were then incubated at 37°C for 30 min. The crystal violet eluted from each biofilm was measured colorimetrically at OD₅₇₀.

2.7 G. mellonella larvae infection

A. baumannii in the log-phase were adjusted to 1× 10⁷ CFU/ml. G. mellonella larvae weighing ~300 mg was divided randomly into groups, with each group containing 15 larvae. Each larva was infected with 10 μl of the adjusted bacterial suspension, incubated at 37°C, and observed once every 12 h for 3 days. Larval survival status was assessed using the acupuncture method, and no response was considered as death.

2.8 Whole genome sequencing and genome analysis

A. baumannii genomes were extracted using QIAamp DNA Mini Kit (QIAGEN, Valencia, CA), whole genome sequencing was performed by Illumina HiSeq 2500 (Illumina, San Diego, CA), and the raw sequencing results were quality-controlled using FastQC v.0.11.5. Trimmomatic v.0.40 was used to trim splice regions. The trimmed reads were assembled and annotated by SPAdes(http://cab.spbu.ru/software/spades/) v.3.6 and RAST(https://rast.nmpdr.org/), respectively. Antibiotic resistance genes and virulence genes were identified using the ResFinder (Bortolaia et al., 2020) and VFDB databases (Chen et al., 2016), respectively, with 80% identity and 80% query coverage as the cutoff values.

2.9 RNA sequencing and quality control

A. baumannii RNA library preparation, construction, sequencing, and processing of reads were performed at the Novogene Co., Ltd. (Beijing, China). The raw data were first processed by a Perl script to remove reads with adapters, those containing poly-N and low-quality reads. The sequencing error rate for a single base position should be less than 1% and no more than 6% (Table S2). All the analyses below were based on high-quality clean data.

2.10 Quantification of gene expression levels and analysis of differentially expressed genes (DEGs)

Clean reads were mapped to A. baumannii MDR-ZJ06 (accession No. CP001937.2) using Bowtie2-2.2.3. Counting the reads numbers mapped to each gene was achieved using HTSeq v0.6.1. Differential expression analysis for the DST and NST groups was performed using the DESeq R package (1.18.0), and p-values were corrected using the Benjamini–Hochberg method to control for false discovery rate (Glickman et al., 2014). Genes for which corrected p-values <0.05 were obtained, after DESeq processing, were defined as differentially expressed.

2.11 GO, COG, and KEGG enrichment analysis

Gene ontology (GO) enrichment analysis of DEGs was performed using the GOseq R package. GO terms with p-values less than 0.05 were considered significantly enriched for DEGs. Subsequently, the COG database was used to identify the functions of the proteins encoded by the DEGs. KOBAS software was used to test the statistical enrichment of DEGs in the KEGG pathway.

2.12 Quantitative real-time PCR

For analysis of expression of specific genes, A. baumannii RNA was first stabilized with RNA protection solution (Qiagen, 74124) and then extracted and purified according to the instructions of the RNeasy Mini Kit (Qiagen, 74104). One microliter RNA was added to nine microliter reaction buffer from the PrimeScript™ RT reagent Kit (Takara, RR037A) in the appropriate ratio, and then reverse transcribed to generate cDNA as follows: 37°C for 15 min, 85°C for 5 sec, and 4°C. cDNA was added to the reaction system of the TB Green^® Premix Ex Taq™ (Takara, 420A), and PCR was carried out as follows: 95°C, 30 s for one cycle (stage 1); 95°C, 5 s, 60°C, 34 s for 40 cycles (stage 2); and 95°C, 15 s, 60°C,1 min, 95°C 15 s (stage 3). The genes and primers used are shown in Table S3, and recA was selected as the internal reference. The relative expression of different genes was calculated by the 2^-ΔΔCt method.

2.13 Statistics

Statistics were performed using prism9. Two groups were compared using the Mann-Whitney t-test, and the five groups were compared using a one-way ANOVA, with horizontal lines representing means ± standard error of means; *p < 0.05, **p < 0.01, ***p < 0.001, **** p < 0.0001.The biological experiments were repeated at least three times, except for the growth curve experiment and anti-desiccation assay, where two biological replicates were done, and the biofilm experiment, where four biological replicates were done.

3 Results

3.1 Comparison of the biological functions of DST and NST groups

We first compared their growth and there was no difference in the growth rates of strains in the two groups (Figure S1, p = 0.42).

As resistance to desiccation and oxidative stress are essential for the persistence of A. baumannii in the medical environment, we compared their abilities in these two aspects. After 96 h of desiccation, strains in the DST group had a mean survival rate of 10.38%, while strains in the NST group had a mean survival rate of 5.13% (p < 0.0001; Figure 1A). Moreover, the oxidative stress analysis showed that the hydrogen-peroxide-induced growth inhibition zone for the DST group was lower than that for the NST group (zone diameter 36.29 mm vs. 37.16 mm; p =0.027), indicating that the DST strains were more resistant to oxidative stress than the NST strains (Figures 1B, C).

FIGURE 1

Figure 1 Comparison of (A) baumannii DST and NST group phenotypes. (A) Antidesiccation in vitro. (B) Antioxidation in vitro. (C) Representative pictures of the antioxidation in vitro. The upper figure is the inhibition circle of hydrogen peroxide in the DST group, while the lower is the inhibition circle in the NST group. (D) Antiserum killing in vitro. (E) Biofilm formation ability of bacteria in vitro. (F) The survival rate of G. mellonella after being infected with A. baumannii at 72 h. Each point represents data from one biological experiment for one bacterium. Each experiment included 10 strains of ST191, 12 strains of ST195, 11 strains of ST208, 10 strains of ST462, and 9 strains of ST547, and bar plots illustrate means ± SEM. Meaning of symbols “*, ***, **** and ns” indicate p<0.05, p<0.001, p<0.0001 and no statistical difference (p>0.05), respectively.

The DST group had higher resistance to complement-mediated killing (mean survival 11.00% vs. 5.31%; p < 0.0001; Figure 1D), but lower biofilm formation ability (OD₅₇₀ 0.31 vs. 0.42; p < 0.0001), than the NST group (Figure 1E).

When G. mellonella was used as an in vivo bacterial infection model, larvae infected with the DST group had a higher mean survival rate than those infected with the NST group in the early (12 h), middle (36 h), and late (72 h) stages of infection, however, all this difference was not statistically significant (Figures 1F, S2).

3.2 Antimicrobial susceptibility

All strains belonging to the DST group were resistant to all three cephalosporins, while the strains of the NST group displayed resistance rates of 33.34%, 66.66%, and 33.34% to ceftazidime, cefepime, and ceftriaxone, respectively. Similarly, the DST strains were all resistant to carbapenems, whereas 66.66% of the NST group strains showed resistance. However, for aminoglycosides, the resistance rate was 66.66% for the DST group and 16.66–33.34% for the NST group. All the DST group strains were resistant to quinolones, while the NST group strains had a resistance rate of 66.66%. The overall MICs of the DST group were also higher than those of the ST group for minocycline and combined antibiotics (Table 1).

TABLE 1

Table 1 Antibiotic susceptibility of A. baumannii DST and NST strains.

3.3 Virulence and resistance genes analysis based on genomes of DST group and NST group

The general features of the twelve genomes were shown in Table S4. The phylogenetic tree built using roary software (default parameters) indicates that strains of the same ST type were closely related, but ST208 of the DST group was evolutionarily distant from ST191 and ST195. ST547 of the NST group was also distant from ST462 and the STn type strains (Figure S3).

The results appeared that the two groups shared significant virulence genes. The encoded products of which included adhesion-related outer membrane protein OmpA, serum resistance-associated penicillin-binding protein PbpG, cleavage cell membrane related phospholipase Plc and PlcD, iron acquisition protein BarAB, BauABCDE, AbaI and AbaR related to quorum sensing, PgaABCD involved in biofilm formation, the CPS synthesis component, Wza, Wzb and Wzc related to immune escape, etc (Figure 2). This is consistent with the similar virulence seen for strains in the two groups in the G. mellonella infection model.

FIGURE 2

Figure 2 Comparison of virulence genes present in A. baumannii DST and NST strains. The first row represents the virulence gene name, and the second row represents gene functional classification. The columns from left to right represent the name and ST type. Both identity and coverage are greater than 80% as the threshold for the presence of genes.

Twenty-four virulence genes associated with CPS were found only in strains of the DST group. These include glycosyltransferase-related genes itrB2, gtr3/4/5/8/20/21/22, kpsS, CMP-glycan pathway-related pseB/C/F/H/I, UDP-glycan pathway-related fnlA/fnlB, dehydrogenase pgt1, and isomerase wecB. However, 17 virulence-associated genes were found exclusively in strains from the NST group. These include 15 genes related to the CPS (CMP-glycan pathway-related legC, lgaA/F/G gene, glycosyltransferase gtr14/15/52, M3Q_295/296, and weeH) and two genes related to iron acquisition.

Overall, strains of the DST group harbored a greater number of CPS synthesis-related genes than those of the NST group, consistent with the increased resistance to complement killing and desiccation previously reported for DST group strains (Russo et al., 2010; Tipton et al., 2018). Many different genes related to the K locus, which is associated with CPS synthesis in A. baumannii, were unique to members of each of the two groups, consistent with the high variability previously reported for the middle region at the K locus (Kenyon and Hall, 2013).

To explore why the DST group strains had higher antibiotic resistance rates than those of the NST group, the profile of resistance genes in each strain was analyzed using the ResFinder database. Although strains in the two groups had several resistance genes in common, including genes encoding aminoglycoside, β-lactam, and sulfonamide resistance, many other resistance genes were found only in the DST group, correlating with the differences in drug susceptibility observed for the two groups. These unique genes included ant(3”)-Ia_1 and armA_1 (aminoglycoside resistance), blaTEM-1D_1 (β-lactam resistance), mph(E)_1 (macrolide resistance), and sul1_5 (sulfonamide resistance; Figure 3). Interestingly, the only resistance genes unique to the NST group were β-lactam-related genes, such as blaOXA-120_1, blaOXA-51_1, and blaOXA-531_1 (Figure 3).

FIGURE 3

Figure 3 Antibiotic resistance gene profiles of A. baumannii DST and NST strains. The first row represents the resistance gene designation, and the second row represents the gene functional classification. The columns from left to right represent the strain name and ST type. Both identity and coverage are greater than 80% as the threshold for the presence of genes.

3.4 Transcriptome analysis

Transcriptome analysis showed that a total of 620 protein-coding genes were differentially expressed between the two groups (padj < 0.05), with 492 genes upregulated and 128 genes downregulated, in DST compared with NST, accounting for approximately 15% of all protein-coding genes (Figure 4). After screening by log₂ FoldChange (DST/NST) ≥ 2, 128 genes were identified as being significantly upregulated in the DST group, including nine CPS synthesis-related genes, one efflux pump gene, and three T6SS secretion system-related genes. Loci corresponding to gene_ids for all these genes were identified in the genomes of the DST strains being analyzed (Table S5). After screening by log₂ FoldChange (DST/NST) ≤-2, 53 genes were identified as being significantly downregulated in the DST group, including three genes related to CPS biosynthesis and six genes related to iron uptake (Table S6). Meanwhile, we randomly selected two genes for Q-PCR validation, respectively. The results showed that the relative trend in the expression of these genes, in strains from the different groups, was consistent with the transcriptomic data (Figures 4B, C).

FIGURE 4

Figure 4 Relative expression of genes in DST and NST strains. (A) The y-axis shows -log₁₀padj values for each gene, while the x-axis shows log₂ fold change of DST vs NST. Red represents upregulation genes, and green represents downregulation genes; FDR ≤ 0.05 was the criterion for significant differences. (B), (C) RNA-seq data of several highly expressed genes were validated using Q-PCR.

GO analysis showed that the most dramatically enriched gene set related to “lipid biosynthesis process,” which contains the wecB, capD, and galE genes involved in CPS synthesis, as well as the fabB, fabF, and psd genes involved in lipid syntheses. Other functional groups enriched in the DST upregulated genes group included “lipid metabolism process” and “carbohydrate transmembrane transport activity” (Figure 5A). The most significantly enriched, downregulated genes in the DST group were “potassium-transporting ATPase activity” and “ATPase activity, coupled to transmembrane movement of ions, phosphorylative mechanism,” with both of these functional sets containing the potassium ion transport genes kdpA and kdpC. Other functional sets with higher levels of enrichment in the DST downregulated genes group included “chemotaxis” and “taxis,” both of which include the chpA and pilI genes (Figure 5B).

FIGURE 5

Figure 5 GO term functional analysis. Vertical coordinates are GO term functional modules; horizontal coordinates represent the number of enriched genes, where green represents biological processes, red represents cellular components, and gray represents molecular functions. We selected the 30 most significantly enriched GO terms to be shown in the figure, or all of them if there are less than 30. (A) GO enrichment of up-regulated DEGs in DST strains; (B) GO enrichment of down-regulated DEGs in DST strains. (B) Enrichment of genes downregulated in DST strains.

Similar to the GO analysis, the COG analysis revealed that the “lipid transport and metabolism” category was enriched in the upregulated genes in the DST group (Figure 6). Additionally, the functional categories with the next largest differences in the number of genes expressed between the two groups were “Transcription” and “Replication,”. Although upregulated genes outnumbered downregulated genes in many functional categories, the “intracellular trafficking, secretion, and vesicular transport” category only contained downregulated genes. This suggests that strong secretion and vesicular transport abilities may be less important to DST strains.

FIGURE 6

Figure 6 COG functional classification of up- and downregulated genes in the DST group relative to the NST group. Y-axis represents the number of up- or downregulated genes, and X-axis represents the COG functional classification.

The KEGG pathway enrichment analysis indicated that the most significant pathway enriched for upregulated genes was “starch and sucrose metabolism” (p = 0.0054), which contains three genes, galU, ugd, and pgi, related to the CPS biosynthesis. This is in agreement with the fact that the most significant terms in the GO analysis also contain three CPS synthesis-related genes. Other enriched pathways were “Pentose phosphate pathway,” “Biotin metabolism,” and “Fatty acid biosynthesis” (Figure 7). The most significant of the downregulated gene-enriched pathways were the “two-component systems” (p=1.20E-07), which include proteins associated with potassium uptake (kdpA, kdpB, and kdpC) and motility (pilJ, pilH, pilI, pilG, PilR).

FIGURE 7

Figure 7 Statistics of pathway enrichment. We selected the 20 most significantly enriched pathway entries and displayed them in the graph. If the enriched pathway entries were less than 20, we displayed all of them. The ordinate represents the name of the pathway, the abscissa represents the Rich factor, the size of the point represents the number of DEGs in this pathway, and the color of the point corresponds to a different Q value. (A) KEGG enrichment of upregulated DEGs in DST strains; (B) KEGG enrichment of downregulated DEGs in DST strains.

3.5 Virulence and resistance genes in DEGs based on transcriptomic analysis

Differences in DEGs associated with virulence and antibiotic resistance were next investigated since these two phenotypes are closely linked to the clinical significance of bacteria and to their survival in the environment.

The VFDB database and the one-way blast approach for SwissProt annotation were combined to identify 42 virulence-related genes among the upregulated genes in strains of the DST group. The genes identified were similar to those identified by genomic analysis, the functions of which related mainly to CPS synthesis and transport (wza/wzb/wzc, etc.), iron uptake (Heme utilization genes, TonB-dependent iron carrier receptors, and iron regulatory proteins), type VI secretion systems (vgrg, vgrgA, contractile small subunits, and PAAR domain-containing proteins), and biofilm formation (pilin/pgaB, ompW, oprM, abaR). In addition, the toxin-antitoxin cdiA/cdiB, RNA chaperone hfq, sensor protein pilS, and the redox-sensitive transcriptional activator soxR were also highly expressed in strains of the DST group (Table S7). The functions of 36 virulence genes identified as downregulated in the DST groups mainly involved type IV pilus synthesis, and included genes encoding pilus assembly proteins (pilC, pilB, pilJ, pilG, pilH, pilI, pilT) and the pilus regulatory protein gene (PilR). In addition, genes downregulated in the DST group contained genes related to iron uptake, CPS synthesis, potassium transport, chemotactic and competition functions, type II secretion system, and phospholipase synthesis (Table S8).

Alignment of DEGs with the antibiotic resistance database ResFinder identified that genes encoding resistance to several antibiotic classes, including aminoglycosides, tetracyclines, β-lactam, chloramphenicol, and cationic antimicrobial peptide, were upregulated in the DST group (Table S9), while chlorhexidine, chloramphenicol, and aminoglycoside resistance genes were downregulated in the DST group (Table S10).

4 Discussion

A. baumannii is a significant cause of pneumonia and bloodstream infection in critically ill patients, and its high rate of multidrug resistance poses a great challenge to clinical treatment. Consequently, research on A. baumannii has focused mainly on revealing its antibiotic resistance mechanism. However, as the understanding of this pathogen deepens, it has become clear that virulence mechanisms of A. baumannii also play a very important role in human infections and disease prognosis. Sequence typing has attracted the attention of public health clinicians because it helps to study the transmission, outbreak, and prevalence of A. baumannii strains. Although sporadic studies have investigated the A. baumannii characteristics that contribute to specific STs becoming DSTs (Ali et al., 2017; Zahra et al., 2018; de Azevedo et al., 2019), research in this area is still in its infancy. Thus, a more integrated approach, which explores the relationship between the virulence, antibiotic resistance, and ST of specific dominant strains, is required.

Based on previously published studies, we identified ST191, ST195, and ST208 as the DSTs and ST 547 and ST462 as the NSTs of A. baumannii, and selected around ten strains of each ST for experimental research in this study. As serum complement killing assay is an important indicator of virulence of A. baumannii, we selected two strains from each of the three DST with high complement killing resistance, and two strains from each of the ST462, ST547 plus two STn that were sensitive to complement killing, for subsequent genomic and transcriptomic analyses.

As we expected, the DST group was more tolerant to various adverse environments, including oxidation, desiccation, and complement-mediated killing. This also correlated with the demonstration that the DST group possessed more genes related to CPS biosynthesis, and that the expression of these genes was higher than related genes in the NST strains. This is in agreement with previous studies, in which CPS was shown to play a vital role in desiccation resistance and anticomplement killing in Acinetobacter (Ophir and Gutnick, 1994; Tipton et al., 2018). Norton et al. (2013) revealed that desiccation stress can induce resistance to antibiotics, corresponding to higher resistance rates in the A. baumannii DST group (Norton et al., 2013).

Previous studies have shown that A. baumannii can form robust biofilms on wounds and medical devices such as endotracheal tubes (Thompson et al., 2014; Greene et al., 2016), as well as increase the tolerance to various extracellular stress (Espinal et al., 2012; Greene et al., 2016). To our surprise, the biofilm-forming ability of the DST group was not as strong as that of the NST group. Moreover, the biofilm formation-related genes adeFGH, bap, csu fimbriae, pga/b/c/d, abaR, and abaI all existed in the genomes of both groups, and there was no significant difference in the transcription of these genes by either group. The evolution of ST dominant strains may partly be driven by their sustained exposure to nutrient-deficient environments, and as such, 24 h biofilm formation experiments in a nutrient-rich environment do not accurately mimic these conditions. Alternatively, other, as yet unidentified, genes of A. baumannii may be important in the ability bacteria of this organism to form biofilms.

The evolutionary process of bacteria, on the one hand, desires to be resistant to as many antibiotics as possible, and on the other hand to retain their ability to survive in unfavorable environments and infect hosts, however, no bacteria can be highly pathogenic, pan-resistant, permanently prevalent and existence in all (Imperi et al., 2011; Da Silva and Domingues, 2017). This may explain why the DST group, which has a higher drug resistance rate and is more resistant to complement killing, has a weaker biofilm formation ability than the NST group.

Among the highly expressed virulence factors genes in the DST group (log₂Foldchange(DST/NST) ≥2) were CPS biosynthesis and T6SS component-associated genes, whereas the downregulated virulence factors genes (log₂Foldchange(DST/NST) ≤-2) were iron uptake related. This may be due to the fact that the capsule promotes the survival of A. baumannii both in the hospital setting and in patients. T6SS also contributes to the dominance of A. baumannii owing to its competition with other bacteria (Repizo et al., 2015; Weber et al., 2016). However, although the iron acquisition system helps bacteria to acquire iron in hosts with iron-deficient environments, it does not necessarily play an important role in promoting dominant clone formation (Penwell et al., 2015; Runci et al., 2019).

Interestingly, many type 4 pili genes related to adherence and motility, including pilR, pilB, pilC, pilJ, pilG, pilH, pilI, and pilT, were downregulated in the DST group, which showed higher resistance to serum complement. By contrast, a previous study showed that type IV pili of A. baumannii were upregulated while growing in serum (Jacobs et al., 2012). There are several conceivable explanations for this apparent anomaly. First, the media used for the growth of A. baumannii differed between the two studies. In our study, transcriptome analysis was carried out following growth of the bacterium in LB broth, whereas in the other study, the organism was grown in the presence of serum. Second, the A. baumannii strains used were different. Third, since glycosylation of pili has been shown to contribute to the survival of P. aeruginosa in the pulmonary environment (Smedley et al., 2005), post-transcriptional modifications may also play important roles in A. baumannii phenotypes; such modifications were not investigated in this study.

The transposon-sequence (Tn-seq) technique was used to identify 50 genes required for survival, in human serum, of a strain of A. baumannii isolated from a case of osteomyelitis (Sanchez-Larrayoz et al., 2017). Several of the genes identified in that study, such as CPS-related genes wza, and wzb, correlate with those we found to be upregulated in the DST group.

Lipid biosynthesis plays an important role in the synthesis of phospholipids in bacteria, and asymmetric bacterial cell membranes are an important permeability barrier that contributes to bacterial resistance to many antibiotics (Zhang and Rock, 2008; Simpson and Trent, 2019). This is consistent with one of our most important findings, namely, that genes associated with lipid biosynthesis and transport pathways are enriched in those genes upregulated in the DST group. Moreover, the upregulated fabB and fabF involved in this pathway exhibit a prolonged effect on fatty acid synthesis, and the knockdown of fabB in E. coli leads to unsaturated fatty acid nutrient deficiency and failure to grow properly (Wang and Cronan, 2004). The knockdown of fabF also reduces motility in P. aeruginosa (Overhage et al., 2007). To date, however, there is little information on the roles in the virulence and adaptation of Acinetobacter baumannii of many of the genes in this pathway.

Naturally, our study has some shortcomings; for example, the number of the strains used in the genome and transcriptome analysis is small and may not be particularly representative. What’s more, it would be impossible to collect and analyze all possible environmental conditions (e.g., temperature, humidity, rainfall, and other climatic factors) that may influence the evolution of the strains.

In conclusion, tolerance to desiccation, oxidation, and complement killing all play an important role in the evolution of DST A. baumannii strains, while CPS biosynthesis, T6SS, and especially lipid synthesis pathway-related genes also play an indispensable role in this process. Our study further revealed the relationship between antibiotic resistance, virulence, and specific dominant ST of A. baumannii, laying a foundation for better prevention, control, and treatment of infections by this organism in the future.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/,PRJNA906176.

Author contributions

XK completed most of the experiments and wrote the first draft of the manuscript. TC and PL contributed in toxicity testing of bacteria analysis. LG and YZ performed the bioinformatics analysis. YX contributed to conception and design of the study and revised the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by the National Key Research and Development Program of China (2021YFC2300300), the Key research and development program of Zhejiang province (No. 2021C03068), the Research Project of Jinan Microecological Biomedicine Shandong Laboratory (JNL-2022006B), and the Fundamental Research Funds for the Central Universities (2022ZFJH003).

Acknowledgments

We thank Jinru Ji, Chaoqun Ying, and Zhiying Liu for their help with sample collection in the laboratory, Wangxiao Zhou for providing assistance with Genome Analysis. Qixia Luo, and Beiwen Zheng for providing assistance with manuscript revision.

Conflict of interest

The 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.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

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

References

Ali, H. M., Salem, M. Z. M., El-Shikh, M. S., Megeed, A. A., Alogaibi, Y. A., Talea, I. A. (2017). Investigation of the virulence factors and molecular characterization of the clonal relations of multidrug-resistant acinetobacter baumannii isolates. J. AOAC Int. 100, 152–158. doi: 10.5740/jaoacint.16-0139

PubMed Abstract | CrossRef Full Text | Google Scholar

Boll, J. M., Tucker, A. T., Klein, D. R., Beltran, A. M., Brodbelt, J. S., Davies, B. W., et al. (2015). Reinforcing lipid a acylation on the cell surface of acinetobacter baumannii promotes cationic antimicrobial peptide resistance and desiccation survival. mBio 6, e00478–e00415. doi: 10.1128/mBio.00478-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Bortolaia, V., Kaas, R. S., Ruppe, E., Roberts, M. C., Schwarz, S., Cattoir, V., et al. (2020). ResFinder 4.0 for predictions of phenotypes from genotypes. J. Antimicrob. Chemother. 75, 3491–3500. doi: 10.1093/jac/dkaa345

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., Zheng, D., Liu, B., Yang, J., Jin, Q. (2016). VFDB 2016: hierarchical and refined dataset for big data analysis–10 years on. Nucleic Acids Res. 44, D694–D697. doi: 10.1093/nar/gkv1239

PubMed Abstract | CrossRef Full Text | Google Scholar

Da Silva, G. J., Domingues, S. (2017). Interplay between colistin resistance, virulence and fitness in acinetobacter baumannii. Antibiotics (Basel) 6(4), 28–38. doi: 10.3390/antibiotics6040028

PubMed Abstract | CrossRef Full Text | Google Scholar

de Azevedo, F., Dutra, V., Nakazato, L., Mello, C. M., Pepato, M. A., de Sousa, A., et al. (2019). Molecular epidemiology of multidrug-resistant acinetobacter baumannii infection in two hospitals in central Brazil: the role of ST730 and ST162 in clinical outcomes. J. Med. Microbiol. 68, 31–40. doi: 10.1099/jmm.0.000853

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, M., Zhu, M. H., Li, J. J., Bi, S., Sheng, Z. K., Hu, F. S., et al. (2014). Molecular epidemiology and mechanisms of tigecycline resistance in clinical isolates of acinetobacter baumannii from a Chinese university hospital. Antimicrob. Agents Chemother. 58, 297–303. doi: 10.1128/AAC.01727-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Espinal, P., Marti, S., Vila, J. (2012). Effect of biofilm formation on the survival of acinetobacter baumannii on dry surfaces. J. Hosp Infect. 80, 56–60. doi: 10.1016/j.jhin.2011.08.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Freire, M. P., de Oliveira Garcia, D., Garcia, C. P., Campagnari Bueno, M. F., Camargo, C. H., Kono Magri, A. S. G., et al. (2016). Bloodstream infection caused by extensively drug-resistant acinetobacter baumannii in cancer patients: high mortality associated with delayed treatment rather than with the degree of neutropenia. Clin. Microbiol. Infect. 22, 352–358. doi: 10.1016/j.cmi.2015.12.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Giannouli, M., Antunes, L. C., Marchetti, V., Triassi, M., Visca, P., Zarrilli, R. (2013). Virulence-related traits of epidemic acinetobacter baumannii strains belonging to the international clonal lineages I-III and to the emerging genotypes ST25 and ST78. BMC Infect. Dis. 13, 282. doi: 10.1186/1471-2334-13-282

PubMed Abstract | CrossRef Full Text | Google Scholar

Glickman, M. E., Rao, S. R., Schultz, M. R. (2014). False discovery rate control is a recommended alternative to bonferroni-type adjustments in health studies. J. Clin. Epidemiol. 67, 850–857. doi: 10.1016/j.jclinepi.2014.03.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Greene, C., Vadlamudi, G., Newton, D., Foxman, B., Xi, C. (2016). The influence of biofilm formation and multidrug resistance on environmental survival of clinical and environmental isolates of acinetobacter baumannii. Am. J. Infect. Control 44, e65–e71. doi: 10.1016/j.ajic.2015.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Greene, C., Wu, J., Rickard, A. H., Xi, C. (2016). Evaluation of the ability of acinetobacter baumannii to form biofilms on six different biomedical relevant surfaces. Lett. Appl. Microbiol. 63, 233–239. doi: 10.1111/lam.12627

PubMed Abstract | CrossRef Full Text | Google Scholar

Hall, B. G., Acar, H., Nandipati, A., Barlow, M. (2014). Growth rates made easy. Mol. Biol. Evol. 31, 232–238. doi: 10.1093/molbev/mst187

PubMed Abstract | CrossRef Full Text | Google Scholar

Harding, C. M., Hennon, S. W., Feldman, M. F. (2018). Uncovering the mechanisms of acinetobacter baumannii virulence. Nat. Rev. Microbiol. 16, 91–102. doi: 10.1038/nrmicro.2017.148

PubMed Abstract | CrossRef Full Text | Google Scholar

Imperi, F., Antunes, L. C., Blom, J., Villa, L., Iacono, M., Visca, P., et al. (2011). The genomics of acinetobacter baumannii: insights into genome plasticity, antimicrobial resistance and pathogenicity. IUBMB Life 63, 1068–1074. doi: 10.1002/iub.531

PubMed Abstract | CrossRef Full Text | Google Scholar

Jacobs, A. C., Sayood, K., Olmsted, S. B., Blanchard, C. E., Hinrichs, S., Russell, D., et al. (2012). Characterization of the acinetobacter baumannii growth phase-dependent and serum responsive transcriptomes. FEMS Immunol. Med. Microbiol. 64, 403–412. doi: 10.1111/j.1574-695X.2011.00926.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Joly-Guillou, M. L. (2005). Clinical impact and pathogenicity of acinetobacter. Clin. Microbiol. Infect. 11, 868–873. doi: 10.1111/j.1469-0691.2005.01227.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kenyon, J. J., 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

Lee, C. R., Lee, J. H., Park, M., Park, K. S., Bae, I. K., Kim, Y. B., et al. (2017). Biology of acinetobacter baumannii: Pathogenesis, antibiotic resistance mechanisms, and prospective treatment options. Front. Cell Infect. Microbiol. 7, 55. doi: 10.3389/fcimb.2017.00055

PubMed Abstract | CrossRef Full Text | Google Scholar

Maragakis, L. L., Perl, T. M. (2008). Acinetobacter baumannii: epidemiology, antimicrobial resistance, and treatment options. Clin. Infect. Dis. 46, 1254–1263. doi: 10.1086/529198

PubMed Abstract | CrossRef Full Text | Google Scholar

Ning, N. Z., Liu, X., Bao, C. M., Chen, S. M., Cui, E. B., Zhang, J. L., et al. (2017). Molecular epidemiology of bla OXA-23 -producing carbapenem-resistant acinetobacter baumannii in a single institution over a 65-month period in north China. BMC Infect. Dis. 17, 14. doi: 10.1186/s12879-016-2110-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Norton, M. D., Spilkia, A. J., Godoy, V. G. (2013). Antibiotic resistance acquired through a DNA damage-inducible response in acinetobacter baumannii. J. Bacteriol 195, 1335–1345. doi: 10.1128/JB.02176-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Ophir, T., Gutnick, D. L. (1994). A role for exopolysaccharides in the protection of microorganisms from desiccation. Appl. Environ. Microbiol. 60, 740–745. doi: 10.1128/aem.60.2.740-745.1994

PubMed Abstract | CrossRef Full Text | Google Scholar

Overhage, J., Lewenza, S., Marr, A. K., Hancock, R. E. (2007). Identification of genes involved in swarming motility using a pseudomonas aeruginosa PAO1 mini-Tn5-lux mutant library. J. Bacteriol 189, 2164–2169. doi: 10.1128/JB.01623-06

PubMed Abstract | CrossRef Full Text | Google Scholar

Penwell, W. F., DeGrace, N., Tentarelli, S., Gauthier, L., Gilbert, C. M., Arivett, B. A., et al. (2015). Discovery and characterization of new hydroxamate siderophores, baumannoferrin a and b, produced by acinetobacter baumannii. Chembiochem 16, 1896–1904. doi: 10.1002/cbic.201500147

PubMed Abstract | CrossRef Full Text | Google Scholar

Repizo, G. D., Gagne, S., Foucault-Grunenwald, M. L., Borges, V., Charpentier, X., Limansky, A. S., et al. (2015). Differential role of the T6SS in acinetobacter baumannii virulence. PLoS One 10, e0138265. doi: 10.1371/journal.pone.0138265

PubMed Abstract | CrossRef Full Text | Google Scholar

Runci, F., Gentile, V., Frangipani, E., Rampioni, G., Leoni, L., Lucidi, M., et al. (2019). Contribution of active iron uptake to acinetobacter baumannii pathogenicity. Infect. Immun. 87, e00755–18. doi: 10.1128/IAI.00755-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Russo, T. A., Luke, N. R., Beanan, J. M., Olson, R., Sauberan, S. L., MacDonald, U., et al. (2010). The K1 capsular polysaccharide of acinetobacter baumannii strain 307-0294 is a major virulence factor. Infect. Immun. 78, 3993–4000. doi: 10.1128/IAI.00366-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanchez-Larrayoz, A. F., Elhosseiny, N. M., Chevrette, M. G., Fu, Y., Giunta, P., Spallanzani, R. G., et al. (2017). Complexity of complement resistance factors expressed by acinetobacter baumannii needed for survival in human serum. J. Immunol. 199, 2803–2814. doi: 10.4049/jimmunol.1700877

PubMed Abstract | CrossRef Full Text | Google Scholar

Schultz, M. B., Pham Thanh, D., Tran Do Hoan, N., Wick, R. R., Ingle, D. J., Hawkey, J., et al. (2016). Repeated local emergence of carbapenem-resistant acinetobacter baumannii in a single hospital ward. Microb. Genom 2, e000050. doi: 10.1099/mgen.0.000050

PubMed Abstract | CrossRef Full Text | Google Scholar

Simpson, B. W., Trent, M. S. (2019). Pushing the envelope: LPS modifications and their consequences. Nat. Rev. Microbiol. 17, 403–416. doi: 10.1038/s41579-019-0201-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Smedley, J. G., 3rd, Jewell, E., Roguskie, J., Horzempa, J., Syboldt, A., Stolz, D. B., et al. (2005). Influence of pilin glycosylation on pseudomonas aeruginosa 1244 pilus function. Infect. Immun. 73, 7922–7931. doi: 10.1128/IAI.73.12.7922-7931.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Son, H. J., Cho, E. B., Bae, M., Lee, S. C., Sung, H., Kim, M. N., et al. (2020). Clinical and microbiological analysis of risk factors for mortality in patients with carbapenem-resistant acinetobacter baumannii bacteremia. Open Forum Infect. Dis. 7, ofaa378. doi: 10.1093/ofid/ofaa378

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, M. G., Black, C. C., Pavlicek, R. L., Honnold, C. L., Wise, M. C., Alamneh, Y. A., et al. (2014). Validation of a novel murine wound model of acinetobacter baumannii infection. Antimicrob. Agents Chemother. 58, 1332–1342. doi: 10.1128/AAC.01944-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Tipton, K. A., Chin, C. Y., Farokhyfar, M., Weiss, D. S., Rather, P. N. (2018). Role of capsule in resistance to disinfectants, host antimicrobials, and desiccation in acinetobacter baumannii. Antimicrob. Agents Chemother. 62, e00188-18. doi: 10.1128/AAC.01188-18

CrossRef Full Text | Google Scholar

Wang, H., Cronan, J. E. (2004). Functional replacement of the FabA and FabB proteins of escherichia coli fatty acid synthesis by enterococcus faecalis FabZ and FabF homologues. J. Biol. Chem. 279, 34489–34495. doi: 10.1074/jbc.M403874200

PubMed Abstract | CrossRef Full Text | Google Scholar

Weber, B. S., Hennon, S. W., Wright, M. S., Scott, N. E., de Berardinis, V., Foster, L. J., et al. (2016). Genetic dissection of the type VI secretion system in acinetobacter and identification of a novel peptidoglycan hydrolase, TagX, required for its biogenesis. mBio 7, e01253. doi: 10.1128/mBio.01253-16

PubMed Abstract | CrossRef Full Text | Google Scholar

World Health Organization. WHO priority pathogens list for research and development of new antibiotics. World Health Organization (WHO). Availbale at: https://www.who.int/news/item/27-02-2017-who-publishes-list-of-bacteriafor-which-new-antibiotics-are-urgently-needed(2June2018).

Google Scholar

Xiao, S. Z., Chu, H. Q., Han, L. Z., Zhang, Z. M., Li, B., Zhao, L., et al. (2016). Resistant mechanisms and molecular epidemiology of imipenem-resistant acinetobacter baumannii. Mol. Med. Rep. 14, 2483–2488. doi: 10.3892/mmr.2016.5538

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoon, E. J., Kim, D., Lee, H., Lee, H. S., Shin, J. H., Uh, Y., et al. (2019). Counter clinical prognoses of patients with bloodstream infections between causative acinetobacter baumannii clones ST191 and ST451 belonging to the international clonal lineage II. Front. Public Health 7, 233. doi: 10.3389/fpubh.2019.00233

PubMed Abstract | CrossRef Full Text | Google Scholar

Zahra, R., Javeed, S., Malala, B., Babenko, D., Toleman, M. A. (2018). Analysis of escherichia coli STs and resistance mechanisms in sewage from Islamabad, Pakistan indicates a difference in e. coli carriage types between south Asia and Europe. J. Antimicrob. Chemother. 73, 1781–1785. doi: 10.1093/jac/dky109

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y. M., Rock, C. O. (2008). Membrane lipid homeostasis in bacteria. Nat. Rev. Microbiol. 6, 222–233. doi: 10.1038/nrmicro1839

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, K., Tang, X., Wang, L., Guo, Z., Xiao, S., Wang, Q., et al. (2018). An emerging clone (ST457) of acinetobacter baumannii clonal complex 92 with enhanced virulence and increasing endemicity in south China. Clin. Infect. Dis. 67, S179–S188. doi: 10.1093/cid/ciy691

PubMed Abstract | CrossRef Full Text | Google Scholar