- 1Department of Biomedical and Biotechnological Sciences, University of Catania, Catania, Italy
- 2Department of Translational Research and New Technology in Medicine and Surgery, Azienda Ospedaliero Universitaria Pisana, University of Pisa, Pisa, Italy
- 3Department of Clinical and Experimental Medicine, University of Catania, Catania, Italy
- 4Departament de Genètica i Microbiologia, Facultat de Biociènces, Universitat Autònoma de Barcelona, Barcelona, Spain
Multidrug-Resistant (MDR) and Extensively Drug Resistant (XDR) Acinetobacter baumannii (Ab) represent a serious cause of healthcare-associated infections worldwide. Currently, the available treatment options are very restricted and colistin-based therapies are last-line treatments of these infections, even though colistin resistant (COLR) Ab have rarely been isolated yet. In bacteria, small non-coding RNAs (sRNAs) have been implicated in regulatory pathways of different biological functions, however, no knowledge exists about the sRNA role on the biological adaptation in COLR Ab. Our study investigated two Italian XDR isogenic colistin-susceptible/resistant (COLS/R) Ab strain-pairs to discover new sRNA signatures. Comparative sRNA transcriptome (sRNAome) analyses were carried out by Illumina RNA-seq using both a Tru-Seq and a Short Insert library, whilst Ab ATCC 17978 and ACICU Reference Genome assembly, mapping, annotation and statistically significant differential expression (q-value ≤ 0.01) of the raw reads were performed by the Rockhopper tool. A computational filtering, sorting only similarly statistically significant differentially expressed (DE) sRNAs mapping on the same gene in both COLR Ab isolates was conducted. COLR vs. COLS sRNAome, analyzed integrating the DE sRNAs obtained from the two different libraries, revealed some statistically significant DE sRNAs in COLR Ab. In detail, we found: (i) two different under-expressed cis-acting sRNAs (AbsRNA1 and AbsRNA2) mapping in antisense orientation the 16S rRNA gene A1S_r01, (ii) one under-expressed cis-acting sRNA (AbsRNA3) targeting the A1S_2505 gene (hypothetical protein), (iii) one under-expressed microRNA-size small RNA fragment (AbsRNA4) and its pre-microAbsRNA4 targeting the A1S_0501 gene (hypothetical protein), (iv) as well as an over-expressed microRNA-size small RNA fragment (AbsRNA5) and its pre-microAbsRNA5 targeting the A1S_3097 gene (signal peptide). Custom TaqMan® probe-based real-time qPCRs validated the expression pattern of the selected sRNA candidates shown by RNA-seq. Furthermore, analysis on sRNA ΔA1S_r01, ΔA1S_2505 as well as the over-expressed A1S_3097 mutants revealed no effects on colistin resistance. Our study, for the first time, found the sRNAome signatures of clinical COLR Ab with a computational prediction of their targets related to protein synthesis, host-microbe interaction and other different biological functions, including biofilm production, cell-cycle control, virulence, and antibiotic-resistance.
Introduction
The Multidrug-Resistant (MDR) Gram-negative pathogens included recently in the WHO black list (Tacconelli et al., 2018), i.e., Acinetobacter, Pseudomonas, and various Enterobacterales (including Klebsiella, E. coli, Serratia, and Proteus) represent a serious health problem worldwide. In particular, focusing on the Acinetobacter baumannii infections, the therapeutic options for their treatment are very limited and mainly draw on colistin-based therapies (Hancock and Chapple, 1999; Zavascki et al., 2007; Vila and Pachón, 2012). Therefore, the increasing use of colistin (COL) and the spread of colistin resistant (COLR) can consequently determine an increase in polymyxin resistance onset (Ko et al., 2007; Gales et al., 2011).
Small non-coding RNAs (sRNAs) [∼30–300 nucleotides (nt) length] have been recognized as a major class of regulatory molecules in bacteria (Barquist and Vogel, 2015). They unroll a regulatory function affecting gene expression – by base pairing to their related target mRNA – modulating transcription, translation, mRNA stability, DNA maintenance, and silencing. Functionally, sRNAs are involved in the regulation of a wide range of physiological responses, reacting to environmental signals, such as pH or temperature shifts (Wassarman, 2002). They can help the modulation of changes in cellular metabolism to optimize use of available nutrients and improve the probability for survival, as well as contributing to virulence (Chabelskaya et al., 2010; Gottesman and Storz, 2011; Sayed et al., 2012; Álvarez-Fraga et al., 2017; Kröger et al., 2018).
sRNAs act in trans or cis depending on the transcription start position where the sRNA is transcribed with respect to their regulated gene. Trans-encoded sRNAs are transcribed at a genetic locus separated from the gene that they regulate and they often work via an imperfect base pairing with their target mRNAs. Conversely, cis-acting sRNAs are transcribed from the same genetic locus that they regulate even thought in an antisense orientation to the target gene. Since transcribed proximally to their target, cis-antisense RNAs share a perfect match with their targets allowing duplex formation and a stringent regulation (Thomason and Storz, 2010; Chang et al., 2015; Georg and Hess, 2018).
A third class of small RNA, the microRNA-size small RNA fragments (∼15–26 nt) have been identified in a few different species of bacteria by next generation sequencing (NGS). These molecules could originate from a cut of a longer precursor (pre-sRNA) in analogy with maturation of the eukaryotic microRNAs (Bloch et al., 2017).
Few mechanisms have been proposed to the underlying colistin resistance in Ab (Li et al., 2006; Adams et al., 2009; Falagas et al., 2010; Moffatt et al., 2010; Arroyo et al., 2011; Cai et al., 2012; Hood et al., 2013; Parra-Millán et al., 2018; Cafiso et al., 2019), however, no investigations have been conducted so far on the role that sRNAs exert on the biological adaptations of colistin resistance acquisition. Few studies have analyzed the sRNA contribution to the A. baumannii biology and antimicrobial resistance. Weiss et al. (2016) offer a detailed view of the sRNA content of Ab and provide new insights into the evolution and role of these regulatory molecules. In detail, using RNA-seq, 78 Ab sRNAs were identified in the AB5075 background, grouped in six classes of similar sRNAs, with one particularly abundant and homologous to regulatory C4 antisense RNAs found in bacteriophages. Sharma et al. (2014) found three new sRNAs, namely AbsR11, AbsR25, and AbsR28, hypothesizing an sRNA involvement in the regulation of antibiotic resistance in bacteria, specifically in cryptic A. baumannii.
Our study aimed to gain new insight on the sRNA signature and biological target of clinical COLR Ab by high-throughput technology RNA-seq.
Our data define the distinctive signatures of COLR Ab sRNAs revealing computational predicted targets involved in the protein synthesis machinery, host-microbe interactions, pathways involved in biofilm production, cell cycle control, virulence, and antibiotic-resistance.
Materials and Methods
Bacterial Strains
Two Italian Extensively Drug Resistant (XDR) isogenic colistin-susceptible/resistant (COLS/R) clinical Ab strain-pairs (1-S/R, 2-S/R) were investigated. Ab strain-pair source, antimicrobial susceptibility, molecular characterization, genomic epidemiology, and phylogeny were previously characterized (Cafiso et al., 2019).
RNA-Seq
To optimize data, RNA-seq was carried out using the Illumina Mi-seq Standard pipeline on two biological replicates consisting of two different libraries of each strain i.e., a Single-end library with 50 bp reads (Short-Insert library) and a Paired-end Read library with 150 bp reads (Tru-Seq library).
For RNA extraction, a single colony of each Ab strain was grown in 10 ml of Cation-adjusted Mueller-Hinton broth (Ca-MHB) (Becton Dickinson) and incubated at 37°C overnight. The overnight cultures were then diluted 1:50 in 50 ml of Ca-MHB in a sterile 250 ml flask and incubated with shaking at 250 rpm at 37°C. Bacterial pellets were harvested at mid-log growth phase (OD600 0.6 ∼1.5 × 107 CFU/ml, 18 h) according to the strain growth-curves (Supplementary Figure S1).
RNA extraction was performed using the NucleoSpin RNA kit (Macherey-Nagel, Düren, Germany) following the manufacturer’s protocol with minor modifications according to the previously published protocols (Cafiso et al., 2019).
The total RNA quality was checked by the 2200 TapeStation RNA Screen Tape device (Agilent, Santa Clara, CA, United States) and the RNA concentrations were determined using an ND-1000 spectrophotometer (NanoDrop, Wilmington, DE, United States). RNA Integrity Number (RIN) values, ranging from 1 to 10 with 10 being the highest quality, was determined by the Agilent TapeStation 2200 system. Only RNAs with preserved 16S and 23S peaks and with RIN values >8 were used for library construction. The RIN values >8 indicated intact and high quality RNA samples usable for downstream applications as previously published (Fleige and Pfaffl, 2006).
Library Preparation and Sequencing
The Tru-Seq library (TS) was an A Paired-end library with reads of 150 bp and average insert size of 350/400 bp. After sequencing, raw reads were processed using FastQC v0.11.2 to check data quality and, then, reads were trimmed by Trimmomatic v.0.33.2 to remove the adapters for Paired-end reads. A minimum base quality of 15 over a 4-bases sliding-window was required. Only trimmed reads with a length above 36 nucleotides were included in the downstream analysis (Cafiso et al., 2019).
The Short-Insert library (SI) was processed with an A Single-end stranded library with reads of 50 bp. After sequencing, raw reads were processed using FastQC v0.11.2 to evaluate data quality. Reads were then trimmed using Trimmomatic v.0.33.2 to remove sequencing adapters for Single-end reads, requiring a minimum base quality of 15 (Phred scale) and a minimum read length of 15 nucleotides. Only trimmed reads were included in the downstream analysis (Cafiso et al., 2019).
Tru-Seq and Short-Insert Library Analysis
TS and SI RNA-seq reads were aligned on Ab ATCC 17978 (CP000521.1) and on Ab ACICU (CP000863.1) Reference Genomes, assembled and quantified using Rockhopper v.2.03 opportunely developed to investigate the bacterial gene structures and transcriptomes as well as validated to identify novel sRNAs (McClure et al., 2013; Tjaden, 2015; Cafiso et al., 2019).
To obtain the data, analyses were carried out with default parameters and verbose output. Rockhopper normalizes read counts for each sample using the upper quartile gene expression level. Starting from the p-values calculated according to the Anders and Huber approach, differentially expressed genes (DEGs) were selected as statistically significant by computing q-values ≤ 0.01 based on the Benjamini–Hochberg correction with a false discovery rate of <1%. In addition, Rockhopper is a tool using biological replicates when available and surrogates when biological replicates for two different conditions are unavailable, considering the two conditions under investigation surrogate replicates for each other (McClure et al., 2013).
Comparative sRNA Prediction
A double computational filtering was carried out on the library analysis data outputs, first for the DE sRNAs in the COLR strains versus their COLS counterparts and, thus, to identify only those mapping on the same gene in both COLR Ab isolates with a statistically significant q-value ≤ 0.01.
Determination of Small RNA Functional Categories
Functional categories of small RNA target genes were investigated by different bioinformatic tools including BLAST, PANTHER (Protein ANalysis THrough Evolutionary Relationships) Classification System, Gene Ontology (GO) Consortium, ExPASy, and KEGG.
RNA-Seq Data Accession Number
RNA-seq data of the two different libraries have been deposited in the NCBI GEO database under study accession no. GSE109951.
I-Tasser ab initio Structure Modeling
For the I-Tasser (Iterative Threading asseMBLY Refinement) analysis (Zhang, 2009; Roy et al., 2010, 2012; Yang and Zhang, 2015), the first N-terminal 120 AA of the Ab ATCC 17978 AbsRNA2 target (A1S_2505) was selected as a target to obtain more significant predictions, as the C-terminal part of the A1S_2505 protein was excluded because of its low coverage with threading templates identified from the PDB-library: MTYQYHDESIVTELPEDTVFVFGSNMAGQHGSGAARVASQ HFGAVEGVGRGWAGQSFAIPTLNEHIQQMPLSQIEHYVEDF KVYAKNHPKMKYFVTALGCGIAGYKVSEIAPLFKGIHHN.
For the analysis of the microAbsRNA4 target (A1S_0501) the whole FASTA sequence of the protein was used as a template.
Validation of sRNA-Seq Expression
AbsRNA3 and microAbsRNA5 were selected to be validated since they are representative of under-expression and over-expression, respectively. Particularly, their RNA-seq expression levels were validated by Custom TaqMan® Small RNA Assays as follow: a single colony of each strain was grown in 10 ml of Ca-MHB (Becton Dickinson) and incubated at 37°C overnight. The overnight cultures were diluted 1:50 in 50 ml of Ca-MHB in a sterile 250 ml flask and incubated with shaking at 250 rpm under normal atmospheric conditions at 37°C. Bacterial pellets were harvested in mid-log growth phase, lysed by adding lysozyme (10 mg/ml) and incubated for 1 h at 37°C. Small RNA extraction was performed by using the mirVanaTM miRNA Isolation Kit (Ambion, Austin, TX, United States) according to the manufacture’s protocol following the enrichment procedure for small RNAs. Extracted sRNAs were quantified using Eppendorf BioPhotometer D30 to assess their quality and to properly dilute them to the amount suggested by Custom TaqMan® Small RNA Assays protocol (Applied BiosystemsTM), then a stem-loop qRT-PCR was performed, one of the most commonly used real-time PCR approaches to quantify sRNAs. The quantification assay was divided into two steps: i) RNA was reverse-transcribed into cDNA using a stem-loop specific primer, and ii) the quantified RT product was consequently used as a template for real-time qPCR with TaqMan® Fast Advanced Master Mix (Applied Biosystems). The stem-loop specific primer used for RT and the specific TaqMan® probes and primers used for real-time qPCR were provided by Custom TaqMan® Small RNA Assay Design Tool on the Applied Biosystems websiteTM. All real-time qPCRs were performed in triplicate, using Agilent AriaMx Real-Time PCR System, with three different biological replicates using one of the thermal profiles suggested by the Custom TaqMan® Small RNA Assays protocol, which provided a first enzymatic activation step of 95°C for 10 min, followed by 40 cycles of 95°C for 15 s and 60°C for 60 s (Salone and Rederstorff, 2015). The expression levels of AbsRNA3 and microAbsRNA5 are shown as the increment/decrement fold-change (FC) in COLR (1-R, 2-R) vs. COLS strains (1-S, 2-S) in RNA-seq and real-time qPCR.
Statistics
AbsRNA3 and microAbsRNA5 expression levels found by Custom TaqMan® probe-based real-time qPCRs were expressed as means ± standard deviations and analyzed by the one-way analysis of variance (ANOVA) using the on-line Free Statistics Calculators-DanielSoper (Soper, 2019) considering a p-value ≤ 0.01 as statistically significant.
Ab ATCC 17978 Mutant Construction
Due to the difficult to manipulate the clinical A. baumannii strain-pairs related to the lack of antibiotic markers for the mutants selection, A. baumannii ATCC 17978 strain was used to generate ΔA1S_r01 and ΔA1S_2505 mutants as well as to overexpress the A1S_3097 gene.
Gene inactivation was carried out as previously described (Aranda et al., 2010). Briefly, an internal fragment of the target gene was PCR-amplified from the A. baumannii ATCC 17978 genome, using the appropriated primers (listed in Table 1). The internal fragment was cloned into the pCR-BluntII-TOPO plasmid (Invitrogen), introduced by electroporation in Escherichia coli DH5α (Clontech) and selected in kanamycin-containing LB plates. Purified plasmids were then introduced by electroporation in A. baumannii ATCC 17978 strain and selected on kanamycin-containing plates. Recombinant clones were confirmed by sequencing (Macrogen) of the PCR products obtained by using the appropriate primers (Table 1). For overexpression, the A1S_3097 gene was cloned, using the indicated primers (Table 1), into XbaI-NcoI sites of the pET-RA vector (Aranda et al., 2010). The recombinant plasmid was introduced in E. coli DH5α and once correct construction was verified by both PCR and sequencing (Macrogen), in ATCC 17978. Finally, A. baumannii transformants overexpressing the A1S_3097 gene were selected on rifampicin- and kanamycin-containing plates and confirmed by PCR with the pETRAFW and pETRARV primers (Table 1).
Results
Comparative Transcriptome Analysis
To define the characterizing sRNA traits of the two COLR vs. COLS clinical Ab strains, a comparative analysis of the DE sRNAs was conducted by a computational double cross-filtering.
RNA-sequencing generated 1,307,792 – 1,175,327 – 1,173,332 – 1,270,020 total reads in 1-S, 1-R, 2-S, and 2-R, respectively, with 96, 96, 72, and 93% mapped reads on Ab ATCC 17978; as well as 97, 97, 76, and 95% reads aligned on Ab ACICU for the TS library, whilst 2,353,045 – 2,041,858 – 1,804,167 – 1,819,349 total reads in 1-S, 1-R, 2-S, and 2-R, respectively, with 57, 56, 53, and 56% mapped reads on Ab ATCC 17978 and 59, 58, 54, and 57% on Ab ACICU for the SI library. RNome structures were previously published (Cafiso et al., 2019).
As shown in Tables 2, 3 and Supplementary Table S1, the comparative statistically significant filtering-analysis of the sRNAome sorting for DE sRNAs of the COLR vs. COLS Ab parental strains returned two different under-expressed cis-acting sRNAs (AbsRNA1 and AbsRNA2) mapping in antisense orientation the 16S rRNA gene A1S_r01, one under-expressed cis-acting sRNA (AbsRNA3) targeting the A1S_ 2505/ACICU_02783 gene (hypothetical protein), one under-expressed microRNA-size small RNA fragment (AbsRNA4) and its pre-microAbsRNA4 targeting the A1S_0501 gene (hypothetical protein), as well as an over-expressed microRNA-size small RNA fragment (AbsRNA5) and its pre-microAbsRNA5 targeting the A1S_3097 gene coding a signal peptide involved in the cytokinin biosynthesis. In detail, the two different cis-acting sRNAs – mapping on the Ab ATCC 17978 Reference Genome – targeted two different positions (Table 2) of the A1S_r01 gene with a size of 35 bp in strain-pair 1 (AbsRNA1) and 39 bp in strain-pair 2 (AbsRNA2), no mapped regions were found on the Ab ACICU Reference Genome. The 70–75 nt cis-acting AbsRNA3 mapped in an analog position both on Ab ATCC 17978 and Ab ACICU Reference Genomes targeting the A1S_2505/ACICU_02783 gene in both the Ab strain-pairs. The 107 nt pre-microAbsRNA4 covering the A1S_0501 mapping on Ab ATCC 17978 Reference Genome and its smaller fragment of 21 bp (microAbsRNA4) were found in Ab strain-pair 2 and in Ab strain-pair 1, respectively. Similarly, a 228 pre-microAbsRNA5 targeting the A1S_3097 gene in Ab strain-pair 2 and its inner fragment of 20 nt in Ab strain-pair 1 were found. Furthermore, both Ab strain-pairs presented the same aforementioned sRNAs with similar expression profiles (over- or under-expression), though the q-value did not allow, in some cases, to return the same fragments in both strain-pairs, considering the two different genomic annotations. Moreover, the only sRNA with a statistically significant expression in Ab ACICU Reference Genome was ACICU_02783 (AbsRNA3) in both strain-pairs (Supplementary Table S1).
The sRNA nucleotide positions (transcription start and stop) reported by Rockhopper tool were shown in Table 2. In addition, none of these sRNAs targeted the 5′ or 3′ untranslated regions (UTR) of their target genes.
I-Tasser
To predict the putative role of AbsRNA2, an A1S_2505/ACICU_02783 conserved domain (CD) BLAST search and I-Tasser ab initio protein structure prediction were computationally investigated. The three-dimensional structure of a protein can be very informative and useful to understand functional characteristics of proteins with unknown functions. This is because the structure of a protein provides the precise molecular details that often facilitate experimental characterization of an expected function. In a case in which there is no expected function, the structure of a protein can be used to facilitate its functional predictions by using the structure as a search template for better-characterized proteins that share regions of structural similarity (Kemege et al., 2011). The CD-BLAST search provided a match with the PHA00684 super family (cl10259) domain related to a protein of unknown function. On the contrary, analyzing the concordances of the highest significant prediction of the I-Tasser TM-align structural alignment and the COACH Predicted biological function, we resolved the structure and biological function as similar to the Orphan Macrodomain Protein (human C6orf130) with O-Acyl-ADP-ribose deacylase activity. In particular, the closest structural similarity of the targeted A1S_2505 was the PDB-Hit 2lgrA (TM-score 0.925) matching the human protein C6orf130, previously published as an Orphan Macrodomain Protein (human C6orf130) with an O-Acyl-ADP-ribose deacylase activity, which catalyzes the deacylation of O-acetyl-ADP-ribose, O-propionyl-ADP-ribose, and O-butyryl-ADP-ribose to produce ADP-ribose (ADP-r) with acetate, propionate, and butyrate, respectively. Due to the structural similarity, we can speculate that A1S_2505/ACICU_02783 could have a function similar to O-Acyl-ADP-ribose deacylase. This structural prediction was also supported by the COACH Predicted biological function, 2l8rA, defining the human protein C6orf130 in complex with ADP-ribose (C-score 0.59), matching the O-acetyl-ADP-ribose deacetylase receptor binding the ADP-ribose in the AA residues – G19 D20 L21 F22 H32 C33 I34 S35 R39 A42 I44 A45 L47 A87 P118 R119 I120 G121 C122 G123 L124 D125 Y150 L152- representing the binding sites. This ligand binding site (BL0101984) showed GO Molecular Functions of purine nucleoside binding (GO:0001883), hydrolase activity (GO:0016787) and deacetylase activity (GO:0019213) as well as the GO Biological Process of the purine nucleoside metabolic process (GO:0042278). Regarding the computational prediction of the putative role of AbsRNA3, the A1S_0501 hypothetical protein referred to an integral membrane component (GO:0016021), whilst the CD BLAST search provided a CD of the cytoskeletal protein RodZ containing Xre-like HTH and DUF4115 domains related to the cell-cycle control, cell division, and chromosome partitioning (cl34261 Superfamily). This result was also supported by the I-Tasser predictions. By LOMETS, the A1S_0501 protein showed homology with the highest Norm Z-score (1.71) and 0.33 coverage with 2wus hit referred as a bacterial structural protein actin MreB that can be complexed with the cell shape protein RodZ. As regards the A1S_0501 GO and the consensus prediction of GO term, obtained from I-Tasser, the A1S_0501 protein showed a molecular function of DNA polymerase activity (GO:0034061) (GO-score 0.38) and DNA binding (GO:0003677) (GO-score 0.34) and the biological process of nucleic acid metabolism (GO:0090304) (GO-score 0.38).
Validation of sRNA Expression Profile
Custom TaqMan® probe-based real-time qPCRs, dedicated for the analysis of bacterial sRNAs, validated and confirmed the RNA-seq expression profiles of two selected sRNA candidates: AbsRNA3 and microAbsRNA5 in both Ab strain-pairs, as shown in Figure 1. In detail, AbsRNA3 had a statistically significant under-expression (p-value ≤ 0.01), whilst microAbsRNA5 showed a statistically significant over-expression (p-value ≤ 0.01) in both COLR Ab strains compared with their COLS parents.
Figure 1. Custom TaqMan probe-based real-time qPCR validation of RNA-seq expression data on colistin resistant (COLR) characterizing sRNAs.a value of 5-Fold Changes (FC) was indicated for incalculable FC due to the presence of a 0 value in one of the strains. ∗p-value ≤ 0.01 obtained by a Student’s T-test and †q-value ≤ 0.01 according to the Rockhopper guidelines were considered statistically significant.
sRNA Target Mutants
No COL MIC changes with respect to the wild-type Ab ATCC 17978 (COL-MIC 1 mg/L) were observed in the ΔA1S_r01 Ab ATCC 17978 and ΔA1S_2505 Ab ATCC 17978 mutants as well as in the WT + A1S_3097 (Table 4).
Discussion
The comparative (sRNAome) integrated to bioinformatics, computational cross-double filtering and experimental validations of COLR versus COLS Ab strain-pair revealed distinctive small RNA signatures in COLR Ab.
Small non-coding RNAs have been identified so far as crucial regulatory elements in bacteria showing high structural diversity and molecular action mechanisms. The most intensively studied prokaryotic sRNA regulators are cis-acting sRNAs or trans-encoded sRNAs, however, the more recent discovery of microRNAs in prokaryotes represents a challenging field of investigation regarding bacterial regulatory mechanisms.
Our clinical COLR Ab were distinguished by 3 cis-acting sRNAs and 2 microRNA-size small RNA fragments involved in different biological networks. These distinctive features span a different area of bacterial biology involving protein synthesis apparatus, host-microbe interactions, biofilm production, cell-cycle control, virulence and antibiotic resistance that emerge as sRNA targets. Notably, they do not seem apparently related to colistin resistance mechanism as demonstrated by our preliminary data on ΔA1S_r01, ΔA1S_2505 and high expressed A1S_3097 Ab ATCC 17978 mutants showing no COL MIC variations of mutants compared to the WT Ab ATCC 17978 – reflecting, however, the wide biological adaptations that the co-existence of colistin resistance and the XDR profile implies.
Regarding the novel AbsRNA1, AbsRNA2 and AbsRNA3, we can speculatively assume that via a cis-antisense regulation mechanism they could post-transcriptionally regulate the target translation. On the contrary, no regulation appears at the transcriptional level as demonstrated by the lack of statistically significant differential expression in these AbsRNA targets according to previously published transcriptomic data (Cafiso et al., 2019). On top, in COLR strains, AbsRNA1 and AbsRNA2 under-expression could indicate that they are characterized by sRNAs modulating the protein synthesis machinery and the amount of active ribosomes in agreement with other previous findings (Fernández-Reyes et al., 2009). The occurrence of the under-expressed AbsRNA3 – targeting the O-acetyl-ADP-ribose – provided evidence that COLR Ab are distinguished by an AbsRNA3 involved in the RNase III inhibition previously associated with different biological functions, including biofilm production, virulence, and antibiotic resistance in Gram-negative bacteria (Chen et al., 2011; Kim et al., 2013; Song et al., 2014). In fact, the O-acetyl-ADP-ribose is a substrate for several related macrodomain proteins such as the human MacroD1, human MacroD2, Escherichia coli YmdB. In E. coli, YmdB is an RNase III inhibitor that modulates many different functions including biofilm formation (Kim et al., 2013), and E. coli adaptive resistance to aminoglycosides via an antibiotic stress-induced sequential modulation of the endoribonucleolytic activity of RNase III and RNase G (Song et al., 2014).
Regarding the microRNA-size small RNA fragments (microAbsRNA4 and its pre-microAbsRNA4 as well as the microAbsRNA5 and its pre-microAbsRNA5), our data suggest that these sRNAs exist in a premature form (pre-microAbsRNA) that act as precursors of their mature form, microAbsRNAs, likely obtained cutting the premature form. Moreover, the under-expressed microAbsRNA4 and its pre-microAbsRNA4 targeting a gene coding a structural membrane protein with a CD similar to cytoskeletal protein RodZ related to the cell-cycle control, cell division and chromosome partitioning, could speculatively address its regulation of these functions. In fact, the GO-term prediction highlighted two possible molecular functions, a DNA polymerase activity and a DNA-binding function. Furthermore, the microAbsRNA4 target (A1S_0501) was previously listed as a transcript that decreased significantly upon exposure to NaCl in A. baumannii, but no relationship with the mechanisms of antimicrobial tolerance in response to monovalent cations was previously found (Hood et al., 2010).
For the over-expressed microAbsRNA5 and its pre-microAbsRNA5, targeting a signal peptide involved in the cytokinin biosynthetic process, we have to keep in mind that A. baumannii is a versatile pathogen that can adhere and invade numerous cell types displaying varying degrees of susceptibility to invasion, stimulating the pro-inflammatory immune response (Choi et al., 2008). The stimuli and signaling pathways implicated in cell death are not yet established; however, they involve imbalanced calcium homeostasis, pro-inflammatory cytokines, and oxidative stress, all traits related to strain virulence (Smani et al., 2011; Mortensen and Skaar, 2012). MicroAbsRNA5 can regulate host-microbe A. baumannii interactions that shape the pathogenesis of Ab infection mediated by the host immune response.
In this study, we experimentally and computationally discovered five statistically significant DE sRNAs characterizing COLR Ab strains speculatively implicated, as cis-antisense sRNAs, in the regulation of protein synthesis machinery via under-expressed AbsRNA1 and AbsRNA2 and in different biological functions, including biofilm production, virulence and aminoglycoside-resistance via an under-expressed AbsRNA3. Likewise, we found two microAbsRNAs that may be involved in cell-cycle control via an under-expressed microAbsRNA4 and in the host-microbe interaction via an over-expressed microAbsRNA5.
Colistin resistance onset in A. baumannii entails dissimilar biological adaptations – not exclusively related to colistin-resistance – supporting the extremely complex and dynamic nature of this life-threatening microorganism and the urgent need to elucidate the role of small RNAs, whose only the tip of the iceberg is known.
This work offers a model for the identification of sRNA signatures and the prediction of their targets in A. baumannii. Although we do not have clear information on their functions yet, our bioinformatic analysis may provide indications regarding the cellular roles of these new sRNAs.
Data Availability Statement
The datasets generated for this study can be found in the NCBI GEO Database: GSE109951.
Author Contributions
VC and SSte conceived and designed the study. VC, SStr, FL, VD, and AZ performed the genomics, transcriptomics, real time qPCR, and bioinformatics. GP contributed to the bioinformatics analysis. JA carried out the mutant construction. All authors analyzed the data and contributed and approved the manuscript.
Funding
This study was supported by a research grant PRIN 2017SFBFER from MIUR, Italy. This study was also supported by grant BIO2016-77011-R from the Ministerio de Economía y Competitividad. JA is a Serra Húnter Fellow, Generalitat de Catalunya, Barcelona, Spain.
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.
Acknowledgments
We wish to thank the Scientific Bureau of the University of Catania (Italy) for language support.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.03075/full#supplementary-material
References
Adams, M. D., Nickel, G. C., Bajaksouzian, S., Lavender, H., Murthy, A. R., Jacobs, M. R., et al. (2009). Resistance to colistin in Acinetobacter baumannii associated with mutations in the PmrAB two-component system. Antimicrob. Agents Chemother. 53, 3628–3634. doi: 10.1128/AAC.00284-09
Álvarez-Fraga, L., Rumbo-Feal, S., Pérez, A., Gómez, M. J., Gayoso, C., Vallejo, J. A., et al. (2017). Global assessment of small RNAs reveals a non-coding transcript involved in biofilm formation and attachment in Acinetobacter baumannii ATCC 17978. PLoS One 12:e0182084. doi: 10.1371/journal.pone.0182084
Aranda, J., Poza, M., Pardo, B. G., Rumbo, S., Rumbo, C., Parreira, J. R., et al. (2010). A rapid and simple method for constructing stable mutants of Acinetobacter baumannii. BMC Microbiol. 10:279. doi: 10.1186/1471-2180-10-279
Arroyo, L. A., Herrera, C. M., Fernandez, L., Hankins, J. V., Trent, M. S., and Hancock, R. E. (2011). The pmrCAB operon mediates polymyxin resistance in Acinetobacter baumannii ATCC 17978 and clinical isolates through phosphoethanolamine modification of lipid A. Antimicrob. Agents Chemother. 55, 3743–3751. doi: 10.1128/AAC.00256-11
Barquist, L., and Vogel, J. (2015). Accelerating discovery and functional analysis of small RNAs with new technologies. Annu. Rev. Genet. 49, 367–394. doi: 10.1146/annurev-genet-112414-054804
Bloch, S., Węgrzyn, A., Węgrzyn, G., and Nejman-Faleńczyk, B. (2017). Small and smaller-sRNAs and MicroRNAs in the regulation of toxin gene expression in prokaryotic cells: a mini-review. Toxins 30:9. doi: 10.3390/toxins9060181
Cafiso, V., Stracquadanio, S., Lo Verde, F., Gabriele, G., Mezzatesta, M. L., Caio, C., et al. (2019). Colistin resistant A. baumannii: genomic and transcriptomic traits acquired under colistin therapy. Front. Microbiol. 9:3195. doi: 10.3389/fmicb.2018.03195
Cai, Y., Chai, D., Wang, R., Liang, B., and Bai, N. (2012). Colistin resistance of Acinetobacter baumannii: clinical reports, mechanisms and antimicrobial strategies. J. Antimicrob. Chemother. 67, 1607–1615. doi: 10.1093/jac/dks084
Chabelskaya, S., Gaillot, O., and Felden, B. (2010). A Staphylococcus aureus small RNA is required for bacterial virulence and regulates the expression of an immune-evasion molecule. PLoS Pathog. 6:e1000927. doi: 10.1371/journal.ppat.1000927
Chang, H., Replogle, J. M., Vather, N., Tsao-Wu, M., Mistry, R., and Liu, J. M. (2015). A cis-regulatory antisense RNA represses translation in Vibrio cholerae through extensive complementarity and proximity to the target locus. RNA Biol. 12, 136–148. doi: 10.1080/15476286.2015.1017203
Chen, D., Vollmar, M., Rossi, M. N., Phillips, C., Kraehenbuehl, R., Slade, D., et al. (2011). Identification of macrodomain proteins as novel O-acetyl-ADP-ribose deacetylases. J. Biol. Chem. 286, 13261–13271. doi: 10.1074/jbc.M110.206771
Choi, C. H., Lee, J. S., Lee, Y. C., Park, T. I., and Lee, J. C. (2008). Acinetobacter baumannii invades epithelial cells and outer membrane protein A mediates interactions with epithelial cells. BMC Microbiol. 10:216. doi: 10.1186/1471-2180-8-216
Falagas, M. E., Rafailidis, P. I., and Matthaiou, D. K. (2010). Resistance to polymyxins: mechanisms, frequency and treatment options. Drug Resist Updat. 13, 132–138. doi: 10.1016/j.drup.2010.05.002
Fernández-Reyes, M., Rodríguez-Falcón, M., Chiva, C., Pachón, J., Andreu, D., and Rivas, L. (2009). The cost of resistance to colistin in Acinetobacter baumannii: a proteomic perspective. Proteomics 9, 1632–1645. doi: 10.1002/pmic.200800434
Fleige, S., and Pfaffl, M. W. (2006). RNA integrity and the effect on the real-time qRT-PCR performance. Mol. Aspects Med. 27, 126–139. doi: 10.1016/j.mam.2005.12.003
Gales, A. C., Jones, R. N., and Sader, H. S. (2011). Contemporary activity of colistin and polymyxin B against a worldwide collection of gram-negative pathogens: results from the SENTRY antimicrobial surveillance program (2006–09). J. Antimicrob. Chemother. 66, 2070–2074. doi: 10.1093/jac/dkr239
Georg, J., and Hess, W. R. (2018). Widespread antisense transcription in prokaryotes. Microbiol. Spectr. 6. doi: 10.1128/microbiolspec.RWR-0029-2018
Gottesman, S., and Storz, G. (2011). Bacterial small RNA regulators: versatile roles and rapidly evolving variations. Cold Spring Harb. Perspect. Biol. 3:a003798. doi: 10.1101/cshperspect.a003798
Hancock, R. E., and Chapple, D. S. (1999). Peptide antibiotics. Antimicrob. Agents Chemother. 43, 1317–1323.
Hood, M. I., Becker, K. W., Roux, C. M., Dunman, P. M., and Skaar, E. P. (2013). Genetic determinants of intrinsic colistin tolerance in Acinetobacter baumannii. Infect Immun. 81, 542–551. doi: 10.1128/IAI.00704-12
Hood, M. I., Jacobs, A. C., Sayood, K., Dunman, P. M., and Skaar, E. P. (2010). Acinetobacter baumannii increases tolerance to antibiotics in response to monovalent cations. Antimicrob. Agents Chemother. 54, 1029–1041. doi: 10.1128/AAC.00963-09
Kemege, K. E., Hickey, J. M., Lovell, S., Battaile, K. P., Zhang, Y., and Hefty, P. S. (2011). Ab initio structural modeling of and experimental validation for chlamydia trachomatis protein CT296 reveal structural similarity to Fe(II) 2-oxoglutarate-dependent enzymes. J. Bacteriol. 193, 6517–6528. doi: 10.1128/JB.05488-11
Kim, T., Lee, J., and Kim, K. S. (2013). Escherichia coli YmdB regulates biofilm formation independently of its role as an RNase III modulator. BMC Microbiol. 24:266. doi: 10.1186/1471-2180-13-266
Ko, K. S., Suh, J. Y., Kwon, K. T., Jung, S. I, Park, K. H., and Kang, C. I. (2007). High rates of resistance to colistin and polymyxin B in subgroups of Acinetobacter baumannii isolates from Korea. J. Antimicrob. Chemother. 60, 1163–1167. doi: 10.1093/jac/dkm305
Kröger, C., MacKenzie, K. D., Alshabib, E. Y., Kirzinger, M. W. B., Suchan, D. M., Chao, T. C., et al. (2018). The primary transcriptome, small RNAs and regulation of antimicrobial resistance in Acinetobacter baumannii ATCC 17978. Nucleic Acids Res. 12, 9684–9698. doi: 10.1093/nar/gky603
Li, J., Rayner, C. R., Nation, R. L., Owen, R. J., Spelman, D., Tan, K. E., et al. (2006). Heteroresistance to colistin in multidrug-resistant Acinetobacter baumannii. Antimicrob. Agents Chemother. 50, 2946–2950. doi: 10.1128/AAC.00103-06
McClure, R., Balasubramanian, D., Sun, Y., Bobrovskyy, M., Sumby, P., Genco, C. A., et al. (2013). Computational analysis of bacterial RNA-seq data. Nucleic Acids Res. 41:e140. doi: 10.1093/nar/gkt444
Moffatt, J. H., Harper, M., Harrison, P., Hale, J. D., Vinogradov, E., Seemann, T., et al. (2010). Colistin resistance in Acinetobacter baumannii is mediated by complete loss of lipopolysaccharide production. Antimicrob Agents Chemother. 54, 4971–4977. doi: 10.1128/AAC.00834-10
Mortensen, B. L., and Skaar, E. P. (2012). Host-microbe interactions that shape the pathogenesis of Acinetobacter baumannii infection. Cell Microbiol. 14, 1336–1344. doi: 10.1111/j.1462-5822.2012.01817.x
Parra-Millán, R., Vila-Farrés, X., Ayerbe-Algaba, R., Varese, M., Sánchez-Encinales, V., Bayó, N., et al. (2018). Synergistic activity of an OmpA inhibitor and colistin against colistin-resistant Acinetobacter baumannii: mechanistic analysis and in vivo efficacy. J. Antimicrob. Chemother. 73, 3405–3412. doi: 10.1093/jac/dky343
Roy, A., Kucukural, A., and Zhang, Y. (2010). I-TASSER: a unified platform for automated protein structure and function prediction. Nat. Protoc. 5, 725–738. doi: 10.1038/nprot.2010.5
Roy, A., Yang, J., and Zhang, Y. (2012). COFACTOR: an accurate comparative algorithm for structure-based protein function annotation. Nucleic Acids Res. 40, W471–W477. doi: 10.1093/nar/gks372
Salone, V., and Rederstorff, M. (2015). Stem-loop RT-PCR based quantification of small non-coding RNAs. Methods Mol. Biol. 2015, 103–108. doi: 10.1007/978-1-4939-2547-6-1
Sayed, N., Jousselin, A., and Felden, B. (2012). A cis-antisense RNA acts in trans in Staphylococcus aureus to control translation of a human cytolytic peptide. Nat. Struct. Mol. Biol. 19, 105–112. doi: 10.1038/nsmb.2193
Sharma, R., Arya, S., Patil, S. D., Sharma, A., Jain, P. K., Navani, N. K., et al. (2014). Identification of novel regulatory small RNAs in Acinetobacter baumannii. PLoS One 9:e93833. doi: 10.1371/journal.pone.0093833
Smani, Y., Docobo-Perez, F., McConnell, M. J., and Pachon, J. (2011). Acinetobacter baumannii-induced lung cell death: role of inflammation, oxidative stress and cytosolic calcium. Microb. Pathog. 50, 224–232. doi: 10.1016/j.micpath.2011.01.008
Song, W., Kim, Y. H., Sim, S. H., Hwang, S., Lee, J. H., Lee, Y., et al. (2014). Antibiotic stress-induced modulation of the endoribonucleolytic activity of RNase III and RNase G confers resistance to aminoglycoside antibiotics in Escherichia coli. Nucleic Acids Res. 42, 4669–4681. doi: 10.1093/nar/gku093
Soper, D. S. (2019). Analysis of Variance (ANOVA) Calculator - One-Way ANOVA from Summary Data [Software]. Available from http://www.danielsoper.com/statcalc (accessed September 18, 2019).
Tacconelli, E., Magrini, N., Carmeli, Y., Harbarth, S., Kahlmeter, G., Kluytmans, J., et al. (2018). Discovery, research, and development of new antibiotics: the WHO priority list of antibiotic-resistant bacteria and tuberculosis. Lancet Infect Dis. 18, 318–327. doi: 10.1016/S1473-3099(17)30753-3
Thomason, M. K., and Storz, G. (2010). Bacterial antisense RNAs: how many are there, and what are they doing? Annu. Rev. Genet. 44, 167–188. doi: 10.1146/annurev-genet-102209-163523
Tjaden, B. (2015). De novo assembly of bacterial transcriptomes from RNA-seq data. Genome Biol. 13:1. doi: 10.1186/s13059-014-0572-2
Vila, J., and Pachón, J. (2012). Therapeutic options for Acinetobacter baumannii infections: an update. Expert. Opin Pharmacother. 13, 2319–2336. doi: 10.1517/14656566.2012.729820
Wassarman, K. M. (2002). Small RNAs in bacteria: diverse regulators of gene expression in response to environmental changes. Cell. 19. 109, 141–144. doi: 10.1016/s0092-8674(02)00717-1
Weiss, A., Broach, W. H., Lee, M. C., and Shaw, L. N. (2016). Towards the complete small RNome of Acinetobacter baumannii. Microb. Genom. 2:e000045. doi: 10.1099/mgen.0.000045
Yang, J., and Zhang, Y. (2015). I-TASSER server: new development for protein structure and function predictions. Nucleic Acids Res. 43, W174–W181. doi: 10.1093/nar/gkv342
Zavascki, A. P., Goldani, L. Z., Li, J., and Nation, R. L. (2007). Polymyxin B for the treatment of multidrug-resistant pathogens: a critical review. J. Antimicrob. Chemother. 60, 1206–1215. doi: 10.1093/jac/dkm357
Keywords: XDR Acinetobacter baumannii, colistin resistance, small RNAs, Illumina RNA-seq, bioinformatics
Citation: Cafiso V, Stracquadanio S, Lo Verde F, Dovere V, Zega A, Pigola G, Aranda J and Stefani S (2020) COLR Acinetobacter baumannii sRNA Signatures: Computational Comparative Identification and Biological Targets. Front. Microbiol. 10:3075. doi: 10.3389/fmicb.2019.03075
Received: 29 October 2019; Accepted: 19 December 2019;
Published: 17 January 2020.
Edited by:
Benjamin Andrew Evans, University of East Anglia, United KingdomReviewed by:
Younes Smani, Institute of Biomedicine of Seville (IBIS), SpainMargarita Poza, Biomedical Research Institute, Spain
Copyright © 2020 Cafiso, Stracquadanio, Lo Verde, Dovere, Zega, Pigola, Aranda and Stefani. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Viviana Cafiso, v.cafiso@unict.it