# Global assessment of small RNAs reveals a non-coding transcript involved in biofilm formation and attachment in Acinetobacter baumannii ATCC 17978

Laura Alvarez-Fraga, Soraya Rumbo-Feal, Astrid Pérez, Manuel J Gómez, Carmen Gayoso, Juan A Vallejo, Emily J Ohneck, Jaione Valle, Luis A Actis, Alejandro Beceiro, et al.

► To cite this version:

Laura Alvarez-Fraga, Soraya Rumbo-Feal, Astrid Pérez, Manuel J Gómez, Carmen Gayoso, et al. Global assessment of small RNAs reveals a non-coding transcript involved in biofilm formation and attachment in Acinetobacter baumannii ATCC 17978. PLoS ONE, 2017, 12 (8), pp.e0182084. 10.1371/journal.pone.0182084. hal-04939934

HAL Id: hal-04939934

https://hal.science/hal-04939934v1

Submitted on 11 Feb 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
RESEARCH ARTICLE

# Global assessment of small RNAs reveals a non-coding transcript involved in biofilm formation and attachment in *Acinetobacter baumannii* ATCC 17978

Laura Álvarez-Fraga¹*, Soraya Rumbo-Feal¹*, Astrid Pérez¹,², Manuel J. Gómez³,⁴, Carmen Gayoso¹, Juan A. Vallejo¹, Emily J. Ohneck², Jaione Valle⁵, Luis A. Actis², Alejandro Beceiro¹, Germán Bou¹*, Margarita Poza¹*

¹ Departamento de Microbiología, Instituto de Investigación Biomédica (INIBIC), Complejo Hospitalario Universitario (CHUAC), A Coruña, Spain, ² Department of Microbiology, Miami University, Oxford, Ohio, United States of America, ³ Department of Molecular Evolution, Center for Astrobiology, INTA-CSIC, Torrejón de Ardoz, Madrid, Spain, ⁴ Centro Nacional de Investigaciones Cardiovasculares Carlos III, Madrid, Spain, ⁵ Departamento de Biofilms Microbianos, Instituto de Agrobiotecnología, Navarra, Spain

* These authors contributed equally to this work.

* german.bou.arevalo@sergas.es (GB); margarita.poza.dominguez@sergas.es (MP)

<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

## OPEN ACCESS

**Citation:** Álvarez-Fraga L, Rumbo-Feal S, Pérez A, Gómez MJ, Gayoso C, Vallejo JA, 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(8): e0182084. https://doi.org/10.1371/journal.pone.0182084

**Editor:** Indranil Biswas, University of Kansas Medical Center, UNITED STATES

**Received:** April 5, 2017

**Accepted:** July 12, 2017

**Published:** August 1, 2017

**Copyright:** © 2017 Álvarez-Fraga et al. This is an open access article distributed under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** Raw data obtained from pyrosequencing have been deposited in the Short Read Archive database (SRA, NCBI) under accession numbers SRX591862, SRX591861 and SRX591774.

**Funding:** This work has been funded by Projects PI12/00552 and PI15/00860 to GB, PI11/01034 to MP, CP13/00226 to AB, and P14/000059 to MP and AB, all integrated in the National Plan for Scientific Research, Development and

## Abstract

Many strains of *Acinetobacter baumannii* have been described as being able to form biofilm. Small non-coding RNAs (sRNAs) control gene expression in many regulatory circuits in bacteria. The aim of the present work was to provide a global description of the sRNAs produced both by planktonic and biofilm-associated (sessile) cells of *A. baumannii* ATCC 17978, and to compare the corresponding gene expression profiles to identify sRNAs molecules associated to biofilm formation and virulence. sRNA was extracted from both planktonic and sessile cells and reverse transcribed. cDNA was subjected to 454-pyrosequencing using the GS-FLX Titanium chemistry. The global analysis of the small RNA transcriptome revealed different sRNA expression patterns in planktonic and biofilm associated cells, with some of the transcripts only expressed or repressed in sessile bacteria. A total of 255 sRNAs were detected, with 185 of them differentially expressed in the different types of cells. A total of 9 sRNAs were expressed only in biofilm cells, while the expression of other 21 coding regions were repressed only in biofilm cells. Strikingly, the expression level of the sRNA 13573 was 120 times higher in biofilms than in planktonic cells, an observation that prompted us to further investigate the biological role of this non-coding transcript. Analyses of an isogenic mutant and over-expressing strains revealed that the sRNA 13573 gene is involved in biofilm formation and attachment to A549 human alveolar epithelial cells. The present work serves as a basis for future studies examining the complex regulatory network that regulate biofilm biogenesis and attachment to eukaryotic cells in *A. baumannii* ATCC 17978.
Technological Innovation 2013-2016 and funded by the ISCIII - General Subdirection of Assessment and Promotion of the Research - European Regional Development Fund (FEDER) "A way of making Europe". Miami University Research Funds from the department of Microbiology supported this work. We also want to thank the Spanish Network for Research in Infectious Diseases (REIPI RD12/0015/0014 to GB), co-financed by the European Development Regional Fund (EDRF) "A Way to Achieve Europe, Instituto de Salud Carlos III, Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Economía y Competitividad as well as Sociedad Española de Enfermedades Infecciosas y Microbiología Clínica (SEIMC). A.P. was financially supported by the Galician Plan for Research, Innovation and Growth (I2C Plan 2012-2016). J.A. Vallejo was financially supported by the Sara Borrell Programme (ISCIII, Spain CD13/00373). S. Rumbo-Feal was financially supported by the Agustí Pumarola Grant (Societat Catalana de Malalties Infeccioses i Microbiologia Clínica, SCMIMC) and Sociedad Española de Enfermedades Infecciosas y Microbiología Clínica (SEIMC).

**Competing interests:** The authors have declared that no competing interests exist.

## Introduction

Acinetobacter baumannii is a non-fermentative, oxidase negative and non-flagellated Gram-negative bacillus. Although it is a normal inhabitant of the human skin flora, the intestinal tract, and the respiratory system, it has been recently described as a dangerous opportunistic pathogen [1]. A. baumannii exhibits a remarkable ability to develop antibiotic resistance that may lead to multiresistant patterns in a short time [2]. Its high genetic versatility facilitates the rapid adaptation to stressing or unfavourable conditions. In the last years, outbreaks caused by resistant strains of A. baumannii have emerged, causing serious health problems [3–11]. Many bacterial species can grow as biofilms: sessile lifestyle forms that increase protection against antimicrobial agents and host defences [12–17]. The growth of A. baumannii on mucous surfaces and medical devices, such as intravascular catheters or endotracheal intubation devices, may result in biofilm formation, increasing the risk of bloodstream and respiratory infections [16]. Biofilms are difficult to eliminate and may cause persistent and recurrent infections. The ability of A. baumannii to adhere and persist as biofilm structures onto surfaces may be the key to understanding its pathogenic mechanisms. In general, processes involved in adhesion and the subsequent formation and development of biofilms are controlled by complex regulatory networks that coordinate the expression of multiple genes [18–24]. Since biofilm formation on abiotic surfaces contributes to the unique survival pattern of A. baumannii in hospital settings, several studies have focused on the description of the proteins and genes involved in the adoption of that lifestyle. However, little is known about the involvement of small non-coding RNA transcripts (sRNAs) in controlling biofilm formation. sRNA molecules serve a wide range of regulatory functions in bacteria and modulate almost every aspect of cell metabolism [25–32]. RNA regulators have several advantages over protein regulators. They are less costly and can be produced faster than proteins, as they are shorter than most mRNAs and do not require the extra step of translation. Moreover, RNA regulators usually act at a level complementary to protein regulators, most often functioning at the post-transcriptional level as opposed to transcription factors that act before sRNAs or enzymes such as kinases or proteases that act after sRNAs. Different combinations of these proteins and RNA regulators can provide a variety of regulatory outcomes, such as extremely tight repression, an expansion in the genes regulated in response to a single signal or, conversely, an increase in the number of signals sensed by a given gene [33]. sRNAs, usually ranging from 50 to 500 nt in length, are typically encoded within intergenic regions and are independently transcribed, using their own promoter and terminator regions [30]. Although regulatory RNAs may control gene expression by a variety of mechanisms, including the modulation of transcription, translation, mRNA stability and DNA maintenance or silencing, most of the studied small RNA regulators act through base pairing with particular RNA targets, usually affecting the translation and stability of mRNAs. Moreover, many sRNAs have been described as involved in the regulation of important bacterial survival responses [30]. Although sRNAs were first discovered in bacteria many years ago, recent technical advances such as multilayered computational searches, deep sequencing or tiled microarrays, have facilitated the discovery of hundreds of potential regulatory sRNA genes [25, 34–40]. In fact, new high throughput sequencing approaches allowed the characterization of complete bacterial transcriptomes and the description of sRNAs in different microorganisms [29, 32, 41–46]. Although studies of sRNAs involved in bacterial infection are still in the early stages, certain sRNAs are related to virulence and pathogenicity in *Salmonella*, *Erwinia*, *Yersinia*, *Vibrio*, *Listeria*, *Pseudomonas*, *Shigella* and *Staphylococcus* [47–54]. Schilling et al. [28] detected a small RNA gene (Aar) in *A. baylyi* involved in the regulation of amino acid metabolism. Iron-regulated sRNAs have been associated to biofilm formation in *Aggregatibacter actinomycetemcomitans* [55]. Recently, Sharma et al. [56] identified novel
regulatory sRNAs in *A. baumannii* by bioinformatic approaches and described three of them as differentially expressed in different phases of the bacterial growth curve. One of them, sRNA AbsR25, was suggested to be involved in regulating the expression of a transporter.

The aim of the present study was to use deep sequencing technologies to characterize the small RNA transcriptome of *A. baumannii* and to gain insight into the mechanisms behind the ability of this organism to form biofilms. The small noncoding RNA transcriptomes of both planktonic and biofilm-associated cells of *A. baumannii* strain ATCC 17978 were then compared to identify differences in sRNA gene expression profiles. One of these sRNA molecules was described as involved in biofilm formation and in attachment to eukaryotic cells.

# Materials and methods

## Bacterial strains and culture conditions

*A. baumannii* ATCC 17978 was routinely grown in Mueller-Hinton (MH) or Luria-Bertani (LB) broth or agar. *Escherichia coli* TG1 was grown in LB broth and used for cloning procedures. *E. coli* OP50 (*Caenorhabditis Genetics Center*) was grown in LB and used for *Caenorhabditis briggsae* virulence assays. All strains were grown at 37°C and stored at -80°C in LB containing 10% glycerol. The concentration of antibiotic used for selection of transformants was 50 µg/mL of kanamycin (Sigma-Aldrich, St. Louis, MO). For obtaining planktonic cells, a single colony of *A. baumannii* ATCC 17978 was isolated on LB agar and grown in 5 mL of LB broth overnight at 37°C in an orbital shaker at 180 rpm. The resulting culture was diluted 100-fold in 500 mL of LB broth in 1L-flasks and incubated under the same conditions. Optical density (OD) was evaluated at 600 nm each 30 min. Cells were harvested at exponential (OD<sub>600nm</sub> = 0.4) and late stationary (OD<sub>600nm</sub> = 2.0) phases of growth after inoculation. Finally, exponential and stationary cells were resuspended in RNA later reagent (Sigma-Aldrich, St. Louis, MO), frozen using liquid nitrogen and stored at -80°C until RNA extraction.

## Biofilm generation in microfermentors

Biofilm formation under flow conditions was performed as previously described by Tormo *et al*. with some modifications [57]. Briefly, *A. baumannii* ATCC 17978 biofilms were obtained in the Fermentation Laboratory of the Agrobiotechnology Institute (Navarra, Spain) under a continuous flow culture system in LB broth medium, consisting of 60-mL microfermentors (Institute Pasteur, Paris, France) with a continuous aeration flow of 40 mL/h of compressed sterile air. Submerged glass slides served as growing substratum. A sample from an overnight culture of *A. baumannii* ATCC 17978 grown in LB broth was used to inoculate the microfermentors, which were then maintained at 37°C for 24 h. Biofilms formed on the slides were removed with a cell scraper and immediately frozen in liquid nitrogen and stored at -80°C.

## Isolation of small RNA molecules

Three independent samples of planktonic bacteria collected at the exponential and stationary growth phases and the biofilm-associated bacteria generated in the microfermentor, all grown in LB media, were reduced to powder using mortar and pestle in the presence liquid nitrogen. RNA isolation was achieved using the *mirVana miRNA Isolation kit* (Ambion) following the manufacturer's protocols. This kit allowed the separation of a fraction enriched in small RNA molecules (sRNA). Ribosomal RNA species were removed using the Microbexpress kit (Ambion). sRNA samples were treated with DNAse I (Invitrogen), purified using phenol-chloroform and concentrated using standard ethanol precipitation. The concentration and the
purity grade of the RNA samples were evaluated using a NanoDrop ND-1000 (Thermo Fisher Scientific) and the integrity and size using a Bioanalyzer 2100 and a Small RNA Analysis Kit (Agilent Technologies Inc., Germany).

## Construction of cDNA libraries

The SOLID small RNA expression kit (SREK, Ambion) was employed for obtaining double stranded cDNA from sRNA following the manufacturer's instructions and using the primers listed in S1 Table. Briefly, sRNA molecules were hybridized and ligated to adaptors, reverse transcribed and treated with RNase H. Then, the second strand was obtained through amplification. Quantification integrity and size determination of samples were done using a Bioanalyzer 2100 as described above.

## Deep sequencing procedures

Three double-stranded cDNA libraries, derived from sRNA samples obtained from liquid cultures in the exponential and stationary phases of growth, as well as biofilm cells were pyrosequenced in the Roche 454 Sequencing Center (Connecticut, USA) using the GS FLX Titanium chemistry and following the manufacturer's protocols. Raw data obtained from pyrosequencing have been deposited in the Short Read Archive database (SRA, NCBI) under accession numbers SRX591862, SRX591861 and SRX591774.

## 454 read processing and gene expression profile comparison

DNA sequences were pre-processed with cross-match (Phred-Pharp-Consed package) and TrimSeq (Emboss package), to mask and trim the linkers added by the SREK kit (S1 Table). The sequences were then aligned with cross-match against the *A. baumannii* ATCC 17978 chromosome and plasmids [Genbank: NC_009085, NC_009083 and NC_009084]. Alignments were combined and processed with *ad-hoc* Perl scripts to identify read-covered regions and to exclude those overlapping with known genes. The average coverage for the remaining read-covered regions was calculated for each of the growing conditions, and normalized according to the number of mapped reads for each sample, to obtain relative abundance values (RA). RA for known sRNAs (tRNAs and 5S rRNAs) were also obtained, and the lowest (that corresponding to an arginine tRNA expressed in the exponential growth sample) was used as the upper threshold to filter out candidate transcribed regions. Remaining chromosomal regions were considered to be truly expressed. To identify those that could be differentially expressed in the different growth conditions, a 2-fold change was used as threshold. Additional details about 454 read data processing can be found in Supplementary material (S1 Text).

## Real time PCR assays

To analyse the expression levels of a set of differentially expressed sRNA molecules *Taqman Micro RNA Reverse Transcription Kit* and *Taqman Universal PCR Master Mix* as well as TaqMan probes and primers were ordered from Applied Biosystems (Life Technologies). Total RNA from planktonic and biofilm cells obtained as described above was isolated using the *RNeasy mini kit* (Qiagen). 200 ng of this RNA were reverse transcribed and 2 µL of the RT mixture were used for real time PCR assays. The expression of *gyrB* was used as a constitutively expressed control. The program followed for reverse transcription consisted in 30 min at 16°C, 30 min at 42°C and 5 min at 85°C. The program followed for real-time PCR was 1 cycle of 10 min at 95°C, 45 cycles of 95°C for 15 s, 60°C for 1 min and final cycle at 40°C for 20 s.
Three biological replicates were tested for each sample. Student's t-test was performed to evaluate the statistical significance of observed differences.

## Cloning of sRNA gene in an expression vector

The gene coding for sRNA 13573, which is over-expressed in biofilm cells, was PCR-amplified using ATCC 17978 genomic DNA as template using specific primers that were designed using the complete genome sequence [Genbank: NC_009085] and included XbaI and/or NcoI restriction sites (S1 Table). The amplicon was then ligated into the pETRA plasmid containing the CTXM promoter [58] and a kanamycin resistance cassette previously introduced into the PstI restriction site using primers listed in S1 Table. Cloning procedures were performed in E. coli TG1. Finally, an insert-containing pETRA derivative was used to transform A. baumannii ATCC 17978 as described by Aranda et al. [58].

## Constructions of knock-out strain

Plasmid pMO130 [Genbank: EU862243], a suicide vector containing *xylE*, *sacB* and a kanamycin resistance marker, was used as described by Hamad et al. [59]. The sRNA 13573 gene was targeted for deletion. Briefly, 900–1,000 bp upstream and downstream fragments flanking the sRNA 13573 gene were cloned into the pMO130 vector using primers listed in S1 Table. The resulting plasmid was used to transform A. baumannii ATCC 17978 by electroporation. Recombinant colonies representing the first crossover event were obtained using a combination of kanamycin selection and visual detection of XylE activity following the catechol-based method. Bright yellow kanamycin resistant colonies were grown overnight in LB supplemented with 15% sucrose and then plated on the same agar medium. Second crossover events were then confirmed by PCR using primers listed in S1 Table.

## Complementation of stable knockout mutant

To complement the isogenic Δ13573 mutant derivative, the cognate wild type allele was amplified from A. baumannii ATCC 17978 genomic DNA using the primers listed in S1 Table and then cloned into the pETRA plasmid as described above [58]. The construction was used to transform the mutant strain. Transformants were selected on kanamycin-containing plates and confirmed by PCR using primers listed in S1 Table. A strain containing the empty pETRA vector was used as control.

## Adhesion to and invasion of A549 human alveolar epithelial cells

Adhesion and invasion abilities were determined following the procedure described by Gaddy et al. with some modifications [22]. Briefly, A549 human alveolar epithelial cells were grown in 5% CO₂ at 37°C in Dulbecco's Modified Eagle Medium (DMEM) (Sigma-Aldrich) supplemented with 10% heat-inactivated fetal bovine serum (FBS) and 1% of penicillin-streptomycin (Gibco). Confluent monolayers were washed twice with saline solution and once with modified Hank's balanced salt solution (mHBSS, same as HBSS but without glucose). Then, A549 cells were infected with 10⁵ bacteria *per* well and incubated for 24 h in mHBSS at 37°C. To determine bacterial adhesion, the infected monolayers were washed three times with saline solution and then lysed in 500 µL of 0.5% sodium deoxycholate. To determine bacterial invasion A549 cells were infected as described above for 24 h and each well was treated for 2 h with gentamicin (256 µg/mL) before washing. Dilutions of the cell lysates were plated onto LB agar and incubated at 37°C for 24 h. Colony forming units were counted to determine the % of bacteria that had attached to or invaded A549 cells at 24 h compared to the growth control, this
being static conditions and same medium without cells as previously described by Álvarez-Fraga *et al.* [60]. Four independent replicates were done. Student's t-test was performed to evaluate the statistical significance of the observed differences.

## Scanning electron microscopy (SEM) of bacterial biofilms on plastic coverslips

Overnight cultures of *A. baumannii* were used to inoculate 5 ml of LB in 50-ml conical tubes at a 1:100 dilution. Sterile polystyrene coverslips were placed in the inoculated 50-mL conical tubes and the tubes were incubated for 48 h at 37°C without shaking as previously described [22]. Coverslips were washed, dehydrated in ethanol, processed with a critical point drier, and sputter coated as described previously [18]. Biofilms formed above, at and below the liquid-air interface were viewed using a Zeiss Supra Gemini Series 35V scanning electron microscope as described previously [18].

## Analysis of biofilms formed on polarized A549 human alveolar cells

A549 human alveolar epithelial cells were routinely maintained as previously described [61, 62]. Costar transwell permeable support polycarbonate membrane 24 well plates (Costar Transwell Polyester Supports, Corning Inc, Corning, New York) were preconditioned 24 h prior to seeding with DMEM on both sides of the membrane and incubated at 37°C and 5% CO₂. DMEM was removed from the conditioned transwell plates and the membranes were seeded with 10⁵ A549 cells per membrane. A549 cells were maintained submerged (DMEM on top and bottom) on the transwell membranes for one week. Following the initial week of submerged growth, DMEM was removed from the top of the membrane to allow the A549 cells to polarize and begin secreting surfactant. Cells were polarized for 2 weeks. One day prior to and for the duration of infection, A549 cells were fed DMEM supplemented with 10% heat inactivated FBS without penicillin-streptomycin (DMEM-). Bacteria, previously grown in LB at 37°C for 24 h in a shaking incubator at 180 rpm, were washed and resuspended in Hank's Buffered Salt Solution (HBSS) (Hyclone Laboratories, Inc, Logan, Utah). A concentration of 10² bacteria was applied to the apical surface of A549 cells by pipetting 1 µL of suspension onto the center of each membrane. The transwell plate was then incubated and maintained for 72 h at 37°C and 5% CO₂. After 72 h, the membranes were washed with HBSS to remove secretions and unattached bacterial cells. The membranes were then fixed for 24 h in 4% formaldehyde-HBSS at 4°C. The membranes were then prepared for SEM using the previously described [18].

## Quantitative biofilm assay

Biofilm formation was quantified following the procedure described by Tomaras *et al.* [18], with some modifications. *A. baumannii* was grown on LB agar for 18 h at 37°C and used to inoculate 5 mL of LB broth. Cultures were grown at 37°C with shaking. Overnight cultures were pelleted, washed and resuspended in 5 mL of LB. A 1:100 dilution of each strain was incubated at 37°C for 12, 24 and 48 h in 15 mL polyethylene tubes. Growth culture was measured at OD₆₀₀ to estimate total cell biomass. Biofilm formation was quantified by staining with crystal violet and solubilisation with ethanol-acetone. The OD₅₈₀/OD₆₀₀ ratio was used to normalize the amount of biofilm formed to the total cell content of each sample tested to avoid variations due to differences in bacterial growth under different experimental conditions. Eight independent replicates were performed. Student's t-test was performed to evaluate the statistical significance of observed differences.
# Results and discussion

## Processing of sRNA 454 sequencing reads

sRNA fractions, obtained from *A. baumannii* ATCC 17978 planktonic cultures in the exponential and stationary phases of growth, and biofilms were pyrosequenced. The total number of reads was 689,097, 502,152 and 627,209 for the exponential, stationary and biofilm samples, respectively (S2 Table). Between 32% and 43% of pre-processed reads could be aligned for each sample (S2 Table). The total number of filtered alignments was around twice the number of aligned reads, indicating that a significant number of reads could align to more than one location with the same score. Additional information can be found in Supplementary material (S1 Text, S1 and S2 Figs, S2, S3 and S4 Tables and S1 Dataset).

## Quantification of the expression level for known sRNA genes

Filtered alignments were processed to estimate the expression level for known sRNA genes described for *A. baumannii* ATCC 17978. Relative abundance (RA) values were calculated for known protein coding genes and for 16S and 23S rRNA genes (S3 Table). tRNA and 5S rRNA genes had RA average values of 1,092 to 1,338 for the three samples, which were more than 1,000-fold and about 20-fold higher than RA values calculated for protein coding genes and 16S and 23S rRNA genes, respectively (S4 Table and S1 Fig), suggesting that sRNA fractions isolated as described in Methods were significantly depleted of mRNA and rRNA. A1S_2909, coding for Leu tRNA, was the known sRNA gene expressed at the highest level, with a RA value of 18,693 in the stationary phase sample. In contrast, A1S_2764, coding for Arg tRNA, was the known sRNA gene expressed at the lowest level, with a RA value of 7.6 in the exponential phase sample. For more details see S1 Text in Supplementary files.

## Identification of new sRNA gene candidates

A total of 26,956 read-covered regions (exprRegs) were detected. To identify read-covered regions overlapping with known genes, the *ad-hoc* script FindOlappingFeatures was used to compare their coordinates with the coordinates of 3,451 known protein, tRNA and rRNA coding genes. The minimal overlap required to consider two features as overlapping was set to 0.001% of the length of any of the features. A total of 21,392 exprRegs were found to overlap with already described genes. The remaining 5,564 exprRegs (S1 Dataset) were considered as new potential sRNA genes. New sRNA gene candidates were identified both in chromosome and plasmids (5,505 and 59 genes, respectively). For more details see S1 Text in Supplementary files.

## Quantification of the expression level of sRNA gene candidates

The set of 5,564 exprRegs that did not overlap with known protein, tRNA and rRNA coding genes had been defined by processing the combined filtered alignments of the reads from the three samples (exponential phase, stationary phase and biofilm) against the *A. baumannii* ATCC 17978 chromosome and plasmids. They provided a common set of candidate sRNA genes whose expression could now be quantified for each of the growing conditions. To quantify the expression of the 5,564 exprRegs, the script mapAlignHits was used again in "expression mode", using the filtered alignments obtained for each sample and the coordinate set for the 5,564 exprRegs as input. RA values were, in average, 500- to 1,000-fold lower than those calculated for known sRNAs (S5 Table). In fact, the distribution of RA values for the set of 5,564 expressed regions indicated that there were about 3,000 expRegs with RA values equal to zero for the biofilm and stationary phase samples and about 1,500 for the exponential phase sample
(S2 Fig). The significant abundance of exprRegs that were not expressed in some of the three growing conditions suggested that many of them could be differentially expressed under the conditions used in this study. Since the aim of this study was to identify new sRNA genes in *A. baumannii* and compare their expression levels in planktonic exponential and stationary phase cells as well as in sessile bacteria, we decided to focus on the subset of sRNA gene candidates with RA values higher than that calculated for A1S_2764 (coding for tRNA-Arg), which was the lowest RA value previously observed for a known sRNA gene as described above. Accordingly, a total of 255 exprRegs were selected because they had RA values equal or higher than 7.6 in some of the three growing conditions (S2 Dataset). Fig 1 shows the genomic location of these 255 exprRegs. Out of them, 108 exprRegs fulfilled the 7.6 requirement in the RNA samples isolated from sessile bacteria, while this requirement was fulfilled in 140 and 121 exprRegs identified in the exponential and stationary phase planktonic samples, respectively (S6 Table and S3 Fig). With the exception of exprReg_29, which was located on plasmid pAB2 (NC_009084), all of them were located in the chromosome.

This number of 255 exprReg selected varies significantly from those observed in a study done by Weiss *et al.*, where 78 sRNAs were detected in *A. baumannii* AB5075 [63]. However, the authors evidenced that only a small number of sRNAs showed high levels of conservation among different *A. baumannii* strains, including ATCC 17978, and that many transcripts are only present either in a small subgroup of strains or exclusively in AB5075 [63]. This issue has been observed in other species such as *Pseudomonas aeruginosa*, where numerous sRNAs exhibit a strain-specific expression pattern [64, 65]. Furthermore, the set of sRNAs identified could depend on the library preparation strategy used, as it was demonstrated by Gómez-Lozano *et al.* using three different approaches to detect novel transcripts of *P. aeruginosa* PAO1 strain [66].

## Identification of differentially expressed sRNA genes

To identify differentially expressed regions among the set of 255 exprRegs, fold changes were calculated and an induction/repression ratio of 4 was chosen arbitrarily as threshold. A total of 185 exprRegs were classified as differentially expressed in biofilm or stationary phase conditions, relative to exponential phase (Table 1). Twenty-eight and 32 exprRegs were up or down-regulated in biofilm cells with respect to both types of planktonic cells, respectively (Fig 1). It is important to remark the high over-expression level of sRNA 13573; about 120-fold higher in biofilm with respect to planktonic cells. Nine sRNA gene candidates (expression regions 1004, 8014, 9803, 11065, 14414, 15677, 18098, 21533, and 24430) were only expressed in biofilm cells. In contrast, 21 putative sRNA genes (expression regions 29, 424, 788, 2345, 2892, 4493, 4750, 5958, 9947, 10022, 12757, 13281, 15861, 17548, 19511, 19890, 23002, 25126, 25193, 25499, and 25840) were only expressed in planktonic cells and repressed in biofilm (Table 1 and Fig 1). It is important to note that the expression region 29 was located in the plasmid pAB2 as mentioned above.

In addition, qRT-PCR analysis of sRNA 13573 and a group of six sRNA only expressed in biofilm cells using proper Taqman probes (S4 Fig) validated the differential production of these non-coding transcripts detected by pyrosequencing as described above.

Sharma *et al.* [56] described 31 predicted sRNAs coding regions in the genome of *A. baumannii* ATCC 17978 using computational approaches. A group of 18 of these 31 predicted sRNAs were also detected as candidate sRNAs in the present work as showed in S1 Dataset. Furthermore, these authors analyzed 10 of the 31 predicted sRNAs by northern blotting of total RNA isolated from *A. baumannii* MTCC1425 (ATCC 15308) cells at different stages of the growth cycle. This analysis showed that, although sRNAs AbsR11 and AbsR28 had almost
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Representation of the 255 putative sRNA expression regions found by deep sequencing in the genome of the *A. baumannii* strain ATCC 17978. Bars indicate their location in the genome. *no DE* black bars indicate no expression differences of the sRNAs analyzed in the three growing conditions. Blue bars indicate sRNAs differentially expressed in stationary phase conditions related to exponential phase. Yellow bars indicate down regulated sRNAs in biofilm related to both types of planktonic cells. Light green bars indicate up regulated sRNAs in biofilm related to both types of planktonic cells. Red bars indicate sRNAs only expressed in planktonic cells and repressed in biofilm. Dark green bars indicate 9 sRNAs only expressed in biofilm and repressed in planktonic cells.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0182084.g001

the same expression level during all phases of planktonic growth, sRNA AbsR25 had no expression in the lag phase and a maximum expression in the exponential phase [56]. The AbsR25 and AbsR28 small transcripts were described in the present work as expression regions, as shown in S1 Dataset, which correspond to exprRegs 10452 and 22397. However, we did not find them as differentially expressed in the conditions analyzed. Moreover, the length of the sequences of these two regions in the ATCC 17978 strain described in the present work and in the MTCC1425 strain described by Sharma *et al.* [56] does not match completely. These variations could be due to genomic differences between these two strains or even to the different bioinformatic approaches used for prediction and analyses of these coding regions. This fact could partly explain the gene expression differences found between Sharma *et al.'s* work and the present study.
<table><caption>Table 1. List of the 185 differentially expressed coding regions.</caption><thead><tr><th>Expression region No.</th><th>Biofilm</th><th>Exponential</th><th>Stationary</th></tr></thead><tbody><tr><td>29**</td><td>1.00</td><td>14.20</td><td>3.39</td></tr><tr><td>165*</td><td>9.3</td><td>1.71</td><td>1.58</td></tr><tr><td>424**</td><td>1.00</td><td>1.97</td><td>20.67</td></tr><tr><td>497</td><td>1.00</td><td>11.91</td><td>1.00</td></tr><tr><td>598</td><td>4.24</td><td>16.09</td><td>1.00</td></tr><tr><td>601</td><td>4.25</td><td>12.83</td><td>1.00</td></tr><tr><td>602</td><td>6.37</td><td>12.42</td><td>1.00</td></tr><tr><td>611</td><td>3.86</td><td>16.89</td><td>1.00</td></tr><tr><td>786</td><td>13.85</td><td>1.42</td><td>25.22</td></tr><tr><td>788**</td><td>1.00</td><td>1.33</td><td>22.54</td></tr><tr><td>795</td><td>25.84</td><td>2.59</td><td>36.24</td></tr><tr><td>801</td><td>14.07</td><td>1.92</td><td>20.31</td></tr><tr><td>833</td><td>1.00</td><td>1.00</td><td>9.28</td></tr><tr><td>1004*</td><td>7.75</td><td>1.00</td><td>1.00</td></tr><tr><td>1136</td><td>7.52</td><td>7.68</td><td>1.69</td></tr><tr><td>1221*</td><td>7.78</td><td>6.59</td><td>1.00</td></tr><tr><td>1368</td><td>34.5</td><td>141.16</td><td>1.11</td></tr><tr><td>1435</td><td>2.53</td><td>1.00</td><td>12.49</td></tr><tr><td>1571*</td><td>13.75</td><td>1.54</td><td>1.73</td></tr><tr><td>2005</td><td>3.67</td><td>13.14</td><td>1.00</td></tr><tr><td>2006</td><td>1.56</td><td>9.58</td><td>1.00</td></tr><tr><td>2247</td><td>3.89</td><td>9.87</td><td>1.00</td></tr><tr><td>2345**</td><td>1.00</td><td>1.69</td><td>7.72</td></tr><tr><td>2350</td><td>36.03</td><td>60.87</td><td>5.13</td></tr><tr><td>2359</td><td>7.99</td><td>13.43</td><td>1.16</td></tr><tr><td>2492</td><td>1.00</td><td>7.71</td><td>1.00</td></tr><tr><td>2675</td><td>37.66</td><td>8.34</td><td>143.60</td></tr><tr><td>2711</td><td>1.22</td><td>8.06</td><td>1.00</td></tr><tr><td>2811*</td><td>8.3</td><td>1.00</td><td>2.58</td></tr><tr><td>2892**</td><td>1.00</td><td>8.52</td><td>1.69</td></tr><tr><td>2988**</td><td>1.65</td><td>8.89</td><td>3.41</td></tr><tr><td>3084</td><td>1.00</td><td>15.32</td><td>1.00</td></tr><tr><td>3275</td><td>1.00</td><td>1.00</td><td>21.31</td></tr><tr><td>3326</td><td>6.37</td><td>27.10</td><td>2.72</td></tr><tr><td>3333</td><td>4.00</td><td>7.90</td><td>1.00</td></tr><tr><td>3380**</td><td>1.74</td><td>8.66</td><td>29.98</td></tr><tr><td>3440*</td><td>15.43</td><td>2.94</td><td>1.00</td></tr><tr><td>3477</td><td>4.92</td><td>10.92</td><td>1.88</td></tr><tr><td>3641</td><td>4.65</td><td>10.11</td><td>1.00</td></tr><tr><td>3667</td><td>13.5</td><td>65.31</td><td>5.71</td></tr><tr><td>4493**</td><td>1.00</td><td>2.26</td><td>10.89</td></tr><tr><td>4505</td><td>1.00</td><td>8.20</td><td>1.00</td></tr><tr><td>4750**</td><td>1.00</td><td>1.50</td><td>9.92</td></tr><tr><td>4776</td><td>1.27</td><td>17.11</td><td>1.00</td></tr><tr><td>4846</td><td>1.00</td><td>1.00</td><td>15.24</td></tr><tr><td>4946*</td><td>9.91</td><td>1.38</td><td>1.00</td></tr><tr><td>5077</td><td>4.26</td><td>18.87</td><td>3.07</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1. (Continued)**</caption><thead><tr><th>Expression region No.</th><th>Biofilm</th><th>Exponential</th><th>Stationary</th></tr></thead><tbody><tr><td>5098**</td><td>1.13</td><td>29.95</td><td>1.54</td></tr><tr><td>5148</td><td>1.00</td><td>1.00</td><td>7.84</td></tr><tr><td>5168</td><td>1.00</td><td>8.00</td><td>1.00</td></tr><tr><td>5286</td><td>1.00</td><td>20.55</td><td>1.00</td></tr><tr><td>5428</td><td>1.00</td><td>1.00</td><td>8.47</td></tr><tr><td>5881</td><td>298.42</td><td>484.14</td><td>15.87</td></tr><tr><td>5954</td><td>1.00</td><td>9.56</td><td>1.00</td></tr><tr><td><b>5958**</b></td><td>1.00</td><td>34.00</td><td>7.70</td></tr><tr><td>5985</td><td>8.67</td><td>1.00</td><td>10.16</td></tr><tr><td>6714</td><td>7.74</td><td>15.43</td><td>1.34</td></tr><tr><td>6904*</td><td>9.90</td><td>5.09</td><td>1.00</td></tr><tr><td>7203**</td><td>1.76</td><td>8.58</td><td>3.96</td></tr><tr><td>7297</td><td>1.00</td><td>1.00</td><td>7.93</td></tr><tr><td>7468</td><td>1.00</td><td>1.00</td><td>18.63</td></tr><tr><td>7562</td><td>3.82</td><td>10.67</td><td>1.00</td></tr><tr><td>7565</td><td>1.00</td><td>1.00</td><td>8.40</td></tr><tr><td>7923</td><td>1.00</td><td>1.00</td><td>31.41</td></tr><tr><td><b>8014*</b></td><td>8.72</td><td>1.00</td><td>1.00</td></tr><tr><td>8136</td><td>453.88</td><td>68.62</td><td>745.54</td></tr><tr><td>8382*</td><td>16.36</td><td>1.71</td><td>1.00</td></tr><tr><td>8776</td><td>1.00</td><td>7.65</td><td>1.00</td></tr><tr><td>9139**</td><td>1.32</td><td>8.64</td><td>1.39</td></tr><tr><td>9255</td><td>49.00</td><td>67.51</td><td>1.63</td></tr><tr><td>9276</td><td>2.07</td><td>1.00</td><td>8.95</td></tr><tr><td>9507</td><td>7.86</td><td>45.51</td><td>4.12</td></tr><tr><td>9656</td><td>11.25</td><td>9.21</td><td>43.47</td></tr><tr><td>9720</td><td>1.00</td><td>1.00</td><td>8.15</td></tr><tr><td><b>9803*</b></td><td>8.74</td><td>1.00</td><td>1.00</td></tr><tr><td>9878**</td><td>1.54</td><td>12.62</td><td>4.23</td></tr><tr><td><b>9947**</b></td><td>1.00</td><td>1.95</td><td>9.52</td></tr><tr><td>10021</td><td>1.37</td><td>1.00</td><td>17.98</td></tr><tr><td><b>10022**</b></td><td>1.00</td><td>1.79</td><td>9.67</td></tr><tr><td>10085</td><td>1.00</td><td>11.00</td><td>1.00</td></tr><tr><td>10453</td><td>5.99</td><td>16.68</td><td>3.05</td></tr><tr><td>10493</td><td>1.00</td><td>1.00</td><td>17.87</td></tr><tr><td>10859</td><td>1.00</td><td>1.00</td><td>18.08</td></tr><tr><td>11025</td><td>1.00</td><td>1.00</td><td>10.35</td></tr><tr><td><b>11065*</b></td><td>8.65</td><td>1.00</td><td>1.00</td></tr><tr><td>11473*</td><td>12.61</td><td>1.80</td><td>1.00</td></tr><tr><td>11798</td><td>1.00</td><td>12.00</td><td>1.00</td></tr><tr><td>11841</td><td>1.00</td><td>1.00</td><td>10.16</td></tr><tr><td>11854</td><td>3.78</td><td>7.81</td><td>1.16</td></tr><tr><td>12046</td><td>2.44</td><td>8.45</td><td>1.00</td></tr><tr><td>12048**</td><td>1.46</td><td>12.61</td><td>1.55</td></tr><tr><td>12150</td><td>3.07</td><td>15.33</td><td>1.00</td></tr><tr><td>12158</td><td>1.00</td><td>1.00</td><td>7.71</td></tr><tr><td>12361</td><td>1.00</td><td>1.00</td><td>13.55</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1. (Continued)**</caption><thead><tr><th>Expression region No.</th><th>Biofilm</th><th>Exponential</th><th>Stationary</th></tr></thead><tbody><tr><td>12370</td><td>6.98</td><td>7.96</td><td>1.00</td></tr><tr><td>12371</td><td>4.40</td><td>10.40</td><td>1.34</td></tr><tr><td>12407</td><td>1.00</td><td>1.00</td><td>8.82</td></tr><tr><td>12460</td><td>1.00</td><td>1.00</td><td>22.09</td></tr><tr><td>12684*</td><td>10.17</td><td>4.62</td><td>1.00</td></tr><tr><td>12706</td><td>13.95</td><td>40.36</td><td>3.98</td></tr><tr><td>12757**</td><td>1.00</td><td>8.31</td><td>1.21</td></tr><tr><td>13053</td><td>5.16</td><td>16.31</td><td>1.00</td></tr><tr><td>13107</td><td>2.55</td><td>12.67</td><td>1.00</td></tr><tr><td>13281**</td><td>1.00</td><td>4.22</td><td>11.56</td></tr><tr><td>13573*</td><td>120.98</td><td>1.55</td><td>1.02</td></tr><tr><td>13621**</td><td>1.46</td><td>2.78</td><td>26.39</td></tr><tr><td>13625</td><td>1.00</td><td>1.00</td><td>7.90</td></tr><tr><td>13631</td><td>1.00</td><td>8.00</td><td>1.00</td></tr><tr><td>13953</td><td>1.00</td><td>1.00</td><td>8.47</td></tr><tr><td>14242</td><td>1.00</td><td>7.60</td><td>1.00</td></tr><tr><td>14286</td><td>1.00</td><td>9.00</td><td>1.00</td></tr><tr><td>14414*</td><td>17.52</td><td>1.00</td><td>1.00</td></tr><tr><td>14887</td><td>1.00</td><td>1.00</td><td>11.99</td></tr><tr><td>15147</td><td>1.00</td><td>1.00</td><td>10.16</td></tr><tr><td>15214</td><td>6.92</td><td>4.50</td><td>22.38</td></tr><tr><td>15677*</td><td>17.48</td><td>1.00</td><td>1.00</td></tr><tr><td>15861**</td><td>1.00</td><td>1.57</td><td>15.83</td></tr><tr><td>15955**</td><td>2.05</td><td>12.55</td><td>3.98</td></tr><tr><td>16106</td><td>2.52</td><td>11.71</td><td>1.00</td></tr><tr><td>16107*</td><td>10.66</td><td>2.50</td><td>1.00</td></tr><tr><td>16117</td><td>2.77</td><td>1.00</td><td>23.09</td></tr><tr><td>16257</td><td>1.00</td><td>1.00</td><td>13.55</td></tr><tr><td>16412</td><td>1.21</td><td>1.00</td><td>10.74</td></tr><tr><td>16574</td><td>1.00</td><td>11.68</td><td>1.00</td></tr><tr><td>16705</td><td>12.10</td><td>16.30</td><td>2.27</td></tr><tr><td>16723</td><td>12.16</td><td>19.32</td><td>3.89</td></tr><tr><td>16724</td><td>8.81</td><td>22.76</td><td>1.00</td></tr><tr><td>16912</td><td>1.00</td><td>1.00</td><td>10.16</td></tr><tr><td>17299**</td><td>1.13</td><td>9.73</td><td>1.31</td></tr><tr><td>17531</td><td>1.00</td><td>1.00</td><td>36.10</td></tr><tr><td>17548**</td><td>1.00</td><td>3.46</td><td>15.52</td></tr><tr><td>17642</td><td>1.00</td><td>1.00</td><td>37.26</td></tr><tr><td>18098*</td><td>9.58</td><td>1.00</td><td>1.00</td></tr><tr><td>18150</td><td>2.12</td><td>1.62</td><td>9.05</td></tr><tr><td>18205</td><td>19.44</td><td>35.13</td><td>1.00</td></tr><tr><td>18260</td><td>5.74</td><td>8.49</td><td>1.19</td></tr><tr><td>18451</td><td>6.23</td><td>11.35</td><td>1.54</td></tr><tr><td>18718</td><td>1.00</td><td>1.00</td><td>8.47</td></tr><tr><td>18763</td><td>1.00</td><td>1.00</td><td>8.31</td></tr><tr><td>18900</td><td>1.54</td><td>8.78</td><td>1.00</td></tr><tr><td>19276</td><td>1.00</td><td>1.00</td><td>11.62</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1. (Continued)**</caption><thead><tr><th>Expression region No.</th><th>Biofilm</th><th>Exponential</th><th>Stationary</th></tr></thead><tbody><tr><td>19347</td><td>1.00</td><td>7.72</td><td>1.00</td></tr><tr><td>19511**</td><td>1.00</td><td>10.26</td><td>1.12</td></tr><tr><td>19604</td><td>6.55</td><td>17.70</td><td>1.00</td></tr><tr><td>19647</td><td>1.00</td><td>1.00</td><td>22.14</td></tr><tr><td>19890**</td><td>1.00</td><td>10.06</td><td>4.29</td></tr><tr><td>19898</td><td>6.21</td><td>22.87</td><td>1.00</td></tr><tr><td>19931</td><td>1.00</td><td>1.00</td><td>10.65</td></tr><tr><td>20262</td><td>1.00</td><td>11.00</td><td>1.00</td></tr><tr><td>20685</td><td>2.89</td><td>8.11</td><td>1.00</td></tr><tr><td>20762*</td><td>11.33</td><td>1.78</td><td>1.00</td></tr><tr><td>20837</td><td>1.00</td><td>1.00</td><td>10.49</td></tr><tr><td>21011</td><td>34.08</td><td>52.04</td><td>5.05</td></tr><tr><td>21223*</td><td>8.51</td><td>1.54</td><td>1.00</td></tr><tr><td>21275</td><td>1.00</td><td>1.00</td><td>8.47</td></tr><tr><td>21392</td><td>6.13</td><td>7.95</td><td>1.31</td></tr><tr><td>21533*</td><td>11.65</td><td>1.00</td><td>1.00</td></tr><tr><td>21637</td><td>4.54</td><td>10.68</td><td>1.00</td></tr><tr><td>21762*</td><td>8.08</td><td>7.45</td><td>1.00</td></tr><tr><td>21850**</td><td>3.14</td><td>20.00</td><td>3.20</td></tr><tr><td>22036</td><td>8.30</td><td>11.04</td><td>1.00</td></tr><tr><td>22275</td><td>1.00</td><td>1.00</td><td>15.24</td></tr><tr><td>22490*</td><td>7.75</td><td>1.88</td><td>1.00</td></tr><tr><td>22686</td><td>2.59</td><td>10.49</td><td>2.30</td></tr><tr><td>22981</td><td>1.00</td><td>1.00</td><td>11.08</td></tr><tr><td>23002**</td><td>1.00</td><td>1.14</td><td>8.00</td></tr><tr><td>23283*</td><td>12.19</td><td>9.52</td><td>1.00</td></tr><tr><td>23620</td><td>1.00</td><td>1.00</td><td>13.55</td></tr><tr><td>23869</td><td>1.00</td><td>11.00</td><td>1.00</td></tr><tr><td>24430*</td><td>8.30</td><td>1.00</td><td>1.00</td></tr><tr><td>25126**</td><td>1.00</td><td>1.21</td><td>11.90</td></tr><tr><td>25193**</td><td>1.00</td><td>9.98</td><td>1.18</td></tr><tr><td>25478</td><td>4.26</td><td>16.55</td><td>2.78</td></tr><tr><td>25499**</td><td>1.00</td><td>1.88</td><td>13.81</td></tr><tr><td>25551</td><td>18.69</td><td>3.92</td><td>46.59</td></tr><tr><td>25571</td><td>14.09</td><td>6.32</td><td>70.99</td></tr><tr><td>25811</td><td>1.96</td><td>1.00</td><td>10.16</td></tr><tr><td>25840**</td><td>1.00</td><td>6.00</td><td>8.36</td></tr><tr><td>25942</td><td>1.00</td><td>1.00</td><td>8.16</td></tr><tr><td>26235</td><td>1.00</td><td>1.00</td><td>8.47</td></tr><tr><td>26247*</td><td>17.86</td><td>16.66</td><td>3.39</td></tr><tr><td>26276</td><td>2.62</td><td>12.32</td><td>1.22</td></tr><tr><td>26388</td><td>1.88</td><td>9.77</td><td>1.38</td></tr><tr><td>26511</td><td>1.00</td><td>8.00</td><td>1.00</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1. (Continued)**</caption><thead><tr><th>Expression region No.</th><th>Biofilm</th><th>Exponential</th><th>Stationary</th></tr></thead><tbody><tr><td>26627*</td><td>15.83</td><td>2.06</td><td>3.82</td></tr></tbody></table>

Differentially expressed coding regions assessed by 454 pyrosequencing corresponding to sRNA molecules and their normalized expression values in sessile (biofilm) and planktonic cells (in exponential and stationary phase of growth).

*sRNA molecules up-regulated in biofilm-related to both planktonic cell types.

**sRNA molecules down-regulated in biofilm related to both planktonic cell types.

Coding regions only expressed or repressed in biofilm are shown in bold.

https://doi.org/10.1371/journal.pone.0182084.t001

Taken together, all these observations strongly indicate that the production of a relatively large number of sRNA is significantly affected when *A. baumannii* cells grow as planktonic or sessile cells. Furthermore, the type and amount of particular sRNA are significantly affected by the growth stage of planktonic bacteria. Thus, sRNA transcripts could play a critical role in the physiology of *A. baumannii* by mechanisms that are poorly understood.

## Involvement of sRNA 13573 in biofilm formation

Special attention was paid to sRNA 13573 due to its high level of expression, *ca.* 120-fold in biofilm cells compared to planktonic cells as assessed by pyrosequencing. Thus, sRNA 13573 was first cloned into the pETRA vector under the control of the CTXM promoter using primers listed in S1 Table. This approach resulted in the over-expression of this coding region as confirmed by qRT-PCR (Table 2). The transformant colonies were grown and their ability to form biofilm was quantified using the wild type strain harbouring the empty pETRA vector as a control. Since the highest rate of biofilm formation was assessed at 48 h, the comparison between the different isogenic strains was performed at 48 h (Fig 2A, S5 Fig). Cells over-expressing sRNA 13573 (13573) showed the highest value of biofilm formation; about 2-fold compared to the strain containing the empty vector (*P* value 0.0008).

Also, the sRNA 13573 region was deleted from the ATCC 17978 genome and the biofilm formation ability was investigated for the corresponding knockout strain showing a no significant decrease (*P* value 0.22) in biofilm formation compared to the wild type strain as shown in Fig 2A. When a pETRA derivative containing sRNA 13573 was introduced into the knockout strain (Δ13573) the biofilm formation ability was higher (*P* value 0.0014) than that shown by the parental strain (17978). Expression level of the sRNA 13573 gene in the parental, the

<table><caption>**Table 2. Expression level of the sRNA 13573 coding region.**</caption><thead><tr><th>Strain</th><th>Expression level of sRNA 13573</th></tr></thead><tbody><tr><td>17978</td><td>0.86 ± 0.19</td></tr><tr><td>Δ13573</td><td>0</td></tr><tr><td>17978 with empty pETRA</td><td>0.47 ± 0.23</td></tr><tr><td>13573</td><td>151.45 ± 50.84</td></tr><tr><td>Δ13573 complemented</td><td>151.70 ± 15.98</td></tr></tbody></table>

Expression level of the sRNA 13573 was determined by qRT-PCR using Taqman probes of the strains: wild type *A. baumannii* (17978), *A. baumannii* ATCC 17978 without sRNA 13573 (Δ13573), *A. baumannii* ATCC 17978 harbouring pETRA (ATCC with empty pETRA), *A. baumannii* ATCC 17978 harbouring pETRA with over-expressed sRNA 13575 (13573) and *A. baumannii* ATCC 17978 without sRNA 13573 harbouring pETRA with sRNA 13575 (Δ13573 complemented).

https://doi.org/10.1371/journal.pone.0182084.t002
<figure><img src="image_5.png" /><figcaption>**Fig 2. Biofilm formation assays.** A) Quantification of biofilm formation by *A. baumannii* ATCC 17978 (17978), *A. baumannii* ATCC 17978 Δ13573 (Δ13573), *A. baumannii* ATCC 17978 harbouring pETRA (17978 with empty pETRA), *A. baumannii* ATCC 17978 harbouring pETRA with sRNA 13575 (13573), and *A. baumannii* ATCC Δ13573 harbouring pETRA with sRNA 13575 (Δ13573 complemented). B) SEM analysis of *A. baumannii* ATCC 17978 (a), *A. baumannii* sRNA 13573 over-producing strain (b) and *A. baumannii* ATCC 17978 Δ13573 (c). Micrographs were taken at 5,000x and bars indicate the scale marks (2 µm).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0182084.g002

Δ13573, the 13573 and the control strains was confirmed by qRT-PCR (Table 2). Growth curves performed in parallel (S6 Fig) showed the same growing profile in all cases. SEM was performed for elucidating biofilm structure in *A. baumannii* ATCC 17978 and the strain over-expressing the sRNA 13573 or the strain lacking the 13573 sRNA (Fig 2B). SEM analysis revealed that cells over-expressing sRNA 13573 were able to develop tridimensional and organized biofilm structures on the surface of plastic coverslips while cells from both the parental and Δ13573 strains showed unorganized cells. However, cells from the Δ13573 strain showed a more disorganized pattern and more dispersion than cells from the parental strain (Fig 2B).
Little is known about the involvement of sRNA molecules in the regulation of biofilm formation and maintenance. To date, some sRNA molecules have been described as regulators involved in biofilm formation in *Pseudomonas*, *Vibrio*, *Salmonella* or *Erwinia* [67, 68]. Most of the sRNAs involved in biofilm have been described in *E. coli* [69] including RprA [70], McaS [71], GcvB [72] and OmrA-B [73], being OmrA-B sRNAs involved in the biosynthesis of the biofilm matrix by controlling the expression of $\sigma^s$ or CsgD, which is in turn related to biofilm formation. Thomason *et al.* [71] described that under specific conditions *E. coli* over-expressing sRNA McaS increased its biofilm formation ability while it was reduced in the corresponding McaS knockout strain. Attending to the SEM analysis described in the present work, differences in biofilm formation can be appreciated between the wild type and the Δ13573 knockout strain although the crystal violet assay did not showed significant differences between both strains. This could be due to the existence of a complex regulation network where sRNA 13573 seems to play a role in biofilm formation although its deletion causes no remarkable effects. A similar effect has been shown for sRNAs RseX, CsrC or SgrS from *E. coli* whose over-expression increased biofilm formation while their deletion did not affect biofilm formation. Bak *et al.* [74] reported that after over-expressing and deleting 99 sRNAs in *E. coli*, in only a few cases the deletion of sRNA genes had effects on biofilm formation. This observation suggests that some sRNA may not be expressed highly enough to affect this cellular process under the experimental conditions used in this study. Alternatively, it is possible that there might be redundancy of sRNAs acting on the process. Furthermore, multiple sRNAs have been identified that modulate the activity of transcriptional regulators important for biofilm formation. CsrA is a RNA-binding protein that represses biofilm formation and its activity is repressed by sRNAs CsrB and CsrC [75, 76]. Also, the transcriptional regulator CsgG, which has been shown to be required for attachment and biofilm formation in *E. coli* is regulated at the mRNA level by several sRNAs, including McaS, RprA or OmrA/OmrB, in response to different environmental cues [71, 73]. Thus, the effect caused by the deletion of a single sRNA could be counteracted by the action of other sRNAs.

## Involvement of the sRNA 13573 in attachment to human cells

Bacterial adhesion to A549 human alveolar cells was measured for the sRNA Δ13573 strain compared with the strain over-expressing sRNA 13573 from pETRA and using i) the 17978 strain, ii) the 17978 strain containing the empty pETRA and iii) the Δ13573 strain harbouring pETRA containing sRNA 13573, as controls. As shown in Fig 3A and considering the 17978 as the sample representing 100% of attached cells, the 17978 strain over-expressing 13573 (13573) attached more (30-fold) to A549 alveolar epithelial cells than the wild type strain harbouring the empty pETRA plasmid (P value < 0.0001). Once more, the 13573 strain revealed a significant phenotypic change when over-expressed, showing an increase in attachment to human cells ability. Besides, an adhesion ability decrease in the Δ13573 strain was detected when compared to the wild type strain being this difference significant (P value 0.036). The Δ13573 complemented derivative reached a higher (13-fold) value than the wild type strain (P value 0.0001). Furthermore, no colonies were detected when invasiveness was checked for all the strains, indicating that all the bacterial counts obtained after a 24 h-incubation were due to attachment (data not shown). Overall data suggest that the deletion of the sRNA 13573 produces a slight but significant change in the adhesion to A549 human epithelial cells and that its over-expression clearly increases attachment abilities.

The SEM analysis of polarized A549 alveolar epithelial cells co-incubated with bacteria revealed that bacteria over-expressing sRNA 13573 were able to form a thick and tridimensional biofilm (micrographs b and e, Fig 3B) over the A549 cell layer while the wild type
<figure><img src="image_6.png" /><figcaption>**Fig 3. Adhesion assays.** A) Attachment to A549 epithelial cells of *A. baumannii* wild type (17978), mutant derivative (Δ13573), mutant derivative over-expressing sRNA 13575 (Δ13573 complemented), harbouring empty pETRA vector (17978 with empty pETRA) and over-expressing sRNA 13575 (13573). The strain 17978 represents 100% of attached bacteria. B) SEM analysis of cells attached to A549 human alveolar cells of *A. baumannii* ATCC 17978 (a and d), *A. baumannii* ATCC 17978 over-expressing sRNA 13575 (b and e) and *A. baumannii* Δ13573 (c and f). Micrographs were taken at 5,000x (a, b and c) and 10,000x magnification (d, e and f). Bars indicate the scale marks (2 and 1 µm).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0182084.g003

showed poor biofilm formation ability (micrographs a and d, Fig 3B). Accordingly to previous results, the capacity of the Δ13573 strain to attach to eukaryotic cells was remarkably lower than the wild type strain (micrographs c and f, Fig 3B).

Little is known about the role that sRNAs can play in the adhesion to eukaryotic cells. For example, sRNA FasX from *Streptococcus* inhibited the expression of a cell surface pilus, therefore reducing the ability of the bacteria to adhere to host cells [77]. Eijkelkamp et al. [78] compared the ability to form biofilms and to adhere to eukaryotic cells of different strains and reported that there is not a clear relationship between these two capacities. Furthermore, previous data revealed that while the truncation of a gene as *csuE* caused a decrease in biofilm formation, the resulting mutant strain was able to adhere more to bronchial epithelial cells [79].
Taken together, these results suggest that there is no direct correlation between biofilm formation on abiotic and biotic surfaces, and that there is wide variation in the cell-surface and cell-cell interactions that result in adherence and biofilm formation by different *A. baumannii* clinical isolates [80]. However, the 13573 sRNA here described showed to be involved in both adhesion and biofilm abilities, being both functions widely related with the pathogenesis of *A. baumannii* [60].

## Conclusions

The global analysis of the small RNA transcriptome of *A. baumannii* ATCC 17978 revealed different sRNA expression patterns in planktonic and biofilm associated cells, with some of the transcripts only expressed or repressed in sessile bacteria. *In vitro* and *in vivo* results demonstrated that the sRNA 13573 is involved in the control of the sessile lifestyle adoption and attachment to eukaryotic cells. The present work serves as a basis for future studies examining the complex regulatory network that regulate biofilm biogenesis and adhesion properties in *A. baumannii* ATCC 17978.

## Supporting information

S1 Text. Additional details about 454 read data processing.

(DOCX)

S1 Table. Oligonucleotides and probes used in the present study.

(DOCX)

S2 Table. Number of initial reads, mapped reads and alignment hits for each sample.

(DOCX)

S3 Table. *A. baumannii* known sRNA and non-sRNA genes used as reference.

(DOCX)

S4 Table. Average, standard deviation, maximal and minimal values for the normalized expression scores calculated for protein coding genes, 16S and 23S rRNA genes and known sRNA genes, in each of the growing conditions.

(DOCX)

S5 Table. Average, standard deviation, maximal and minimal values for the normalized expression scores calculated for 5564 expressed regions not overlapping with known genes.

(DOCX)

S6 Table. Sets of expression regions not overlapping with known genes and having a normalized expression score equal or higher than 7.6 in some of the three growing conditions.

(DOCX)

S1 Dataset. Expressed regions. List of 5,564 expressed regions, their length and locations in the genome (NC_009085.1) or plasmids (NC_009083.1 and NC_009084.1) of *A. baumannii* ATCC 17978.

(DOCX)

S2 Dataset. Coding regions. List of normalized expression of 255 coding regions (expression score should be equal or above 7.6 in some growing condition: Bio, Exp or Sta), their length and locations in the genome (NC_009085.1) or plasmid (NC_009084.1) of *A. baumannii* ATCC 17978.

(DOCX)
**S1 Fig. Distribution of normalized expression scores for known sRNA genes.** Blue: biofilm samples. Orange: exponential phase samples. Yellow: stationary phase samples.
(DOCX)

**S2 Fig. Distribution of normalized expression scores for the 5564 expressed regions not overlapping with known genes.** Blue: biofilm samples. Orange: exponential phase samples. Yellow: stationary phase samples.
(DOCX)

**S3 Fig. Distribution of normalized expression scores for the 255 expressed regions not overlapping with known genes and having a normalized expression score equal or higher than 7.6 in some of the growing conditions.** Blue: biofilm samples. Orange: exponential phase samples. Yellow: stationary phase samples.
(DOCX)

**S4 Fig. qRT-PCR assays.** Expression levels of 8 sRNA regions in planktonic and sessile cells determined by qRT-PCR using Taqman probes. Y axis represents the relative expression of the genes taking the housekeeping gene *gyrB* as value 1.
(DOCX)

**S5 Fig. Biofilm formation at different times.** Staining at 12, 24 and 48 h of biofilm formation of *A. baumannii* ATCC 17978 (17978), *A. baumannii* ATCC 17978 Δ13573 (Δ13573), *A. baumannii* ATCC Δ13573 harbouring pETRA with sRNA 13575 (Δ13573 complemented), *A. baumannii* ATCC 17978 harbouring pETRA with sRNA 13575 (13573), and *A. baumannii* ATCC 17978 harbouring pETRA (17978 with empty pETRA).
(DOCX)

**S6 Fig. Growth curves.** Growth curves of wild type *A. baumannii* (17978), *A. baumannii* ATCC 17978 lacking the sRNA 13573 (Δ13573), *A. baumannii* Δ13573 over-expressing sRNA 13575 (Δ13573 complemented), *A. baumannii* ATCC 17978 harbouring pETRA over-expressing sRNA 13575 (13573) and *A. baumannii* ATCC 17978 harbouring the empty pETRA vector (17978 with empty pETRA).
(DOCX)

## Acknowledgments

This work has been funded by Projects PI12/00552 and PI15/00860 to GB, PI11/01034 to MP, CP13/00226 to AB, and P14/000059 to MP and AB, all integrated in the National Plan for Scientific Research, Development and Technological Innovation 2013–2016 and funded by the ISCIII—General Subdirection of Assessment and Promotion of the Research—European Regional Development Fund (FEDER) “A way of making Europe”. Miami University Research Funds from the department of Microbiology supported this work. We also want to thank the Spanish Network for Research in Infectious Diseases (REIPI RD12/0015/0014 to GB), co-financed by the European Development Regional Fund (EDRF) “A Way to Achieve Europe, Instituto de Salud Carlos III, Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Economía y Competitividad as well as Sociedad Española de Enfermedades Infecciosas y Microbiología Clínica (SEIMC). A.P. was financially supported by the Galician Plan for Research, Innovation and Growth (I2C Plan 2012–2016). J.A. Vallejo was financially supported by the Sara Borrell Programme (ISCIII, Spain CD13/00373). S. Rumbo-Feal was financially supported by the Agustí Pumarola Grant (Societat Catalana de Malalties Infeccioses i Microbiologia Clínica, SCMIMC) and Sociedad Española de Enfermedades
Infecciosas y Microbiología Clínica (SEIMC). Finally, we thank Chua Kim Lee (University of Singapore) for providing pMo130-TelR.

**Author Contributions**

**Conceptualization:** Astrid Pérez, Carmen Gayoso, Emily J. Ohneck, Margarita Poza.

**Data curation:** Laura Álvarez-Fraga, Soraya Rumbo-Feal, Astrid Pérez, Manuel J. Gómez, Carmen Gayoso, Juan A. Vallejo, Emily J. Ohneck, Jaione Valle.

**Formal analysis:** Laura Álvarez-Fraga, Soraya Rumbo-Feal, Astrid Pérez, Manuel J. Gómez, Carmen Gayoso, Juan A. Vallejo, Emily J. Ohneck, Jaione Valle.

**Funding acquisition:** Luis A. Actis, Alejandro Beceiro, Germán Bou, Margarita Poza.

**Investigation:** Laura Álvarez-Fraga, Soraya Rumbo-Feal, Manuel J. Gómez, Carmen Gayoso, Juan A. Vallejo, Emily J. Ohneck, Jaione Valle.

**Methodology:** Laura Álvarez-Fraga, Soraya Rumbo-Feal, Astrid Pérez, Manuel J. Gómez, Carmen Gayoso, Juan A. Vallejo, Emily J. Ohneck, Jaione Valle.

**Resources:** Luis A. Actis, Germán Bou.

**Software:** Manuel J. Gómez.

**Supervision:** Luis A. Actis, Alejandro Beceiro, Germán Bou, Margarita Poza.

**Validation:** Luis A. Actis, Alejandro Beceiro, Margarita Poza.

**Writing – original draft:** Margarita Poza.

**Writing – review & editing:** Germán Bou, Margarita Poza.

**References**

1. Peleg AY, Seifert H, Paterson DL (2008) Acinetobacter baumannii: emergence of a successful pathogen. Clin Microbiol Rev 21: 538–582. https://doi.org/10.1128/CMR.00058-07 PMID: 18625687

2. Villar M, Cano ME, Gato E, Garnacho-Montero J, Miguel Cisneros J, et al. (2014) Epidemiologic and clinical impact of Acinetobacter baumannii colonization and infection: a reappraisal. Medicine (Baltimore) 93: 202–210.

3. Corbella X, Montero A, Pujol M, Dominguez MA, Ayats J, et al. (2000) Emergence and rapid spread of carbapenem resistance during a large and sustained hospital outbreak of multiresistant Acinetobacter baumannii. J Clin Microbiol 38: 4086–4095. PMID: 11060073

4. Bou G, Cervero G, Dominguez MA, Quereda C, Martinez-Beltran J (2000) PCR-based DNA fingerprinting (REP-PCR, AP-PCR) and pulsed-field gel electrophoresis characterization of a nosocomial outbreak caused by imipenem- and meropenem-resistant Acinetobacter baumannii. Clin Microbiol Infect 6: 635–643. PMID: 11284921

5. Jones ME, Draghi DC, Thornsberry C, Karlowsky JA, Sahm DF, et al. (2004) Emerging resistance among bacterial pathogens in the intensive care unit—a European and North American Surveillance study (2000–2002). Ann Clin Microbiol Antimicrob 3: 14. https://doi.org/10.1186/1476-0711-3-14 PMID: 15283864

6. del Mar Tomas M, Cartelle M, Pertega S, Beceiro A, Linares P, et al. (2005) Hospital outbreak caused by a carbapenem-resistant strain of Acinetobacter baumannii: patient prognosis and risk-factors for colonization and infection. Clin Microbiol Infect 11: 540–546. https://doi.org/10.1111/j.1469-0691.2005.01184.x PMID: 15966971

7. Dijkshoorn L, Nemec A, Seifert H (2007) An increasing threat in hospitals: multidrug-resistant Acinetobacter baumannii. Nat Rev Microbiol 5: 939–951. https://doi.org/10.1038/nrmicro1789 PMID: 18007677

8. Ong CW, Lye DC, Khoo KL, Chua GS, Yeoh SF, et al. (2009) Severe community-acquired Acinetobacter baumannii pneumonia: an emerging highly lethal infectious disease in the Asia-Pacific. Respirology 14: 1200–1205. https://doi.org/10.1111/j.1440-1843.2009.01630.x PMID: 19909464
9. Higgins PG, Dammhayn C, Hackel M, Seifert H (2010) Global spread of carbapenem-resistant Acinetobacter baumannii. J Antimicrob Chemother 65: 233-238. https://doi.org/10.1093/jac/dkp428 PMID: 19996144

10. Acosta J, Merino M, Viedma E, Poza M, Sanz F, et al. (2011) Multidrug-resistant Acinetobacter baumannii Harboring OXA-24 carbapenemase, Spain. Emerg Infect Dis 17: 1064–1067. https://doi.org/10.3201/eid/1706.091866 PMID: 21749771

11. Domenech de Celles M, Salomon J, Marinier A, Lawrence C, Gaillard JL, et al. (2012) Identifying more epidemic clones during a hospital outbreak of multidrug-resistant Acinetobacter baumannii. PLoS One 7: e45758. https://doi.org/10.1371/journal.pone.0045758 PMID: 23029226

12. Raad II, Mohamed JA, Reitzel RA, Jiang Y, Dvorak TL, et al. (2011) The prevention of biofilm colonization by multidrug-resistant pathogens that cause ventilator-associated pneumonia with antimicrobial-coated endotracheal tubes. Biomaterials 32: 2689-2694. https://doi.org/10.1016/j.biomaterials.2010.12.015 PMID: 21295343

13. Costerton JW, Stewart PS, Greenberg EP (1999) Bacterial biofilms: a common cause of persistent infections. Science 284: 1318-1322. PMID: 10334980

14. Vidal R, Dominguez M, Urrutia H, Bello H, Gonzalez G, et al. (1996) Biofilm formation by Acinetobacter baumannii. Microbials 86: 49-58. PMID: 8771775

15. Lee HW, Koh YM, Kim J, Lee JC, Lee YC, et al. (2008) Capacity of multidrug-resistant clinical isolates of Acinetobacter baumannii to form biofilm and adhere to epithelial cell surfaces. Clin Microbiol Infect 14: 49-54. https://doi.org/10.1111/j.1469-0691.2007.01842.x PMID: 18005176

16. Rodriguez-Bano J, Marti S, Soto S, Fernandez-Cuenca F, Cisneros JM, et al. (2008) Biofilm formation in Acinetobacter baumannii: associated features and clinical implications. Clin Microbiol Infect 14: 276-278. https://doi.org/10.1111/j.1469-0691.2007.01916.x PMID: 18190568

17. Marti S, Nait Chabane Y, Alexandre S, Coquet L, Vila J, et al. (2011) Growth of Acinetobacter baumannii in pellicle enhanced the expression of potential virulence factors. PLoS One 6: e26030. https://doi.org/10.1371/journal.pone.0026030 PMID: 22046254

18. Tomaras AP, Dorsey CW, Edelmann RE, Actis LA (2003) Attachment to and biofilm formation on abiotic surfaces by Acinetobacter baumannii: involvement of a novel chaperone-usher pili assembly system. Microbiology 149: 3473-3484. https://doi.org/10.1099/mic.0.26541-0 PMID: 14663080

19. Gaddy JA, Actis LA (2009) Regulation of Acinetobacter baumannii biofilm formation. Future Microbiol 4: 273-278. https://doi.org/10.2217/fmb.09.5 PMID: 19327114

20. Tomaras AP, Flagler MJ, Dorsey CW, Gaddy JA, Actis LA (2008) Characterization of a two-component regulatory system from Acinetobacter baumannii that controls biofilm formation and cellular morphology. Microbiology 154: 3398-3409. https://doi.org/10.1099/mic.0.2008/019471-0 PMID: 18957593

21. Loehfelm TW, Luke NR, Campagnari AA (2008) Identification and characterization of an Acinetobacter baumannii biofilm-associated protein. J Bacteriol 190: 1036-1044. https://doi.org/10.1128/JB.01416-07 PMID: 18024522

22. Gaddy JA, Tomaras AP, Actis LA (2009) The Acinetobacter baumannii 19606 OmpA protein plays a role in biofilm formation on abiotic surfaces and in the interaction of this pathogen with eukaryotic cells. Infect Immun 77: 3150-3160. https://doi.org/10.1128/IAI.00096-09 PMID: 19470746

23. Cabral MP, Soares NC, Aranda J, Parreira JR, Rumbo C, et al. (2011) Proteomic and functional analyses reveal a unique lifestyle for Acinetobacter baumannii biofilms and a key role for histidine metabolism. J Proteome Res 10: 3399-3417. https://doi.org/10.1021/pr101299j PMID: 21612302

24. Rumbo-Feal S, Gomez MJ, Gayoso C, Alvarez-Fraga L, Cabral MP, et al. (2013) Whole transcriptome analysis of Acinetobacter baumannii assessed by RNA-sequencing reveals different mRNA expression profiles in biofilm compared to planktonic cells. PLoS One 8: e72968. https://doi.org/10.1371/journal.pone.0072968 PMID: 24023660

25. Sittka A, Lucchini S, Papenfort K, Sharma CM, Rolle K, et al. (2008) Deep sequencing analysis of small noncoding RNA and mRNA targets of the global post-transcriptional regulator, Hfq. PLoS Genet 4: e1000163. https://doi.org/10.1371/journal.pgen.1000163 PMID: 18725932

26. Liu JM, Livny J, Lawrence MS, Kimball MD, Waldor MK, et al. (2009) Experimental discovery of sRNAs in Vibrio cholerae by direct cloning, 5S/tRNA depletion and parallel sequencing. Nucleic Acids Res 37: e46. https://doi.org/10.1093/nar/gkp080 PMID: 19223322

27. Guillier M, Gottesman S, Storz G (2006) Modulating the outer membrane with small RNAs. Genes Dev 20: 2338-2348. https://doi.org/10.1101/gad.1457506 PMID: 16951250

28. Schilling D, Findeiss S, Richter AS, Taylor JA, Gerischer U (2010) The small RNA Aar in Acinetobacter baylyi: a putative regulator of amino acid metabolism. Arch Microbiol 192: 691-702. https://doi.org/10.1007/s00203-010-0592-6 PMID: 20559624
29. Koo JT, Alleyne TM, Schiano CA, Jafari N, Lathem WW (2011) Global discovery of small RNAs in Yersinia pseudotuberculosis identifies Yersinia-specific small, noncoding RNAs required for virulence. Proc Natl Acad Sci U S A 108: E709–E717. https://doi.org/10.1073/pnas.1101655108 PMID: 21876162

30. Waters LS, Storz G (2009) Regulatory RNAs in bacteria. Cell 136: 615–628. https://doi.org/10.1016/j.cell.2009.01.043 PMID: 19239884

31. Swiercz JP, Hindra, Bobek J, Haiser HJ, Di Berardo C, et al. (2008) Small non-coding RNAs in Streptomyces coelicolor. Nucleic Acids Res 36: 7240–7251. https://doi.org/10.1093/nar/gkn898 PMID: 19008244

32. Filiatrault MJ, Stodghill PV, Bronstein PA, Moll S, Lindeberg M, et al. (2010) Transcriptome analysis of Pseudomonas syringae identifies new genes, noncoding RNAs, and antisense activity. J Bacteriol 192: 2359–2372. https://doi.org/10.1128/JB.01445-09 PMID: 20190049

33. Shimoni Y, Friedlander G, Hetzroni G, Niv G, Altuvia S, et al. (2007) Regulation of gene expression by small non-coding RNAs: a quantitative view. Mol Syst Biol 3: 138. https://doi.org/10.1038/msb4100181 PMID: 17893699

34. Kwenda S, Gorshkov V, Ramesh AM, Naidoo S, Rubagotti E, et al. (2016) Discovery and profiling of small RNAs responsive to stress conditions in the plant pathogen Pectobacterium atrosepticum. BMC Genomics 17: 47. https://doi.org/10.1186/s12864-016-2376-0 PMID: 26753530

35. Khoo JS, Chai SF, Mohamed R, Nathan S, Firdaus-Raih M (2012) Computational discovery and RT-PCR validation of novel Burkholderia conserved and Burkholderia pseudomallei unique sRNAs. BMC Genomics 13 Suppl 7: S13.

36. Livny J, Teonadi H, Livny M, Waldor MK (2008) High-throughput, kingdom-wide prediction and annotation of bacterial non-coding RNAs. PLoS One 3: e3197. https://doi.org/10.1371/journal.pone.0003197 PMID: 18787707

37. Yao Z, Barrick J, Weinberg Z, Neph S, Breaker R, et al. (2007) A computational pipeline for high-throughput discovery of cis-regulatory noncoding RNA in prokaryotes. PLoS Comput Biol 3: e126. https://doi.org/10.1371/journal.pcbi.0030126 PMID: 17616982

38. Pitman S, Cho KH (2015) The Mechanisms of Virulence Regulation by Small Noncoding RNAs in Low GC Gram-Positive Pathogens. Int J Mol Sci 16: 29797–29814. https://doi.org/10.3390/ijms161226194 PMID: 26694351

39. Tjaden B (2012) Computational identification of sRNA targets. Methods Mol Biol 905: 227–234. https://doi.org/10.1007/978-1-61779-949-5_14 PMID: 22736007

40. Vogel J, Sharma CM (2005) How to find small non-coding RNAs in bacteria. Biol Chem 386: 1219–1238. https://doi.org/10.1515/BC.2005.140 PMID: 16336117

41. Khandige S, Kronborg T, Uhlin BE, Moller-Jensen J (2015) sRNA-Mediated Regulation of P-Fimbriae Phase Variation in Uropathogenic Escherichia coli. PLoS Pathog 11: e1005109. https://doi.org/10.1371/journal.ppat.1005109 PMID: 26291711

42. De Lay NR, Garsin DA (2016) The unmasking of ‘junk’ RNA reveals novel sRNAs: from processed RNA fragments to marooned riboswitches. Curr Opin Microbiol 30: 16–21. https://doi.org/10.1016/j.mib.2015.12.006 PMID: 26771674

43. Barquist L, Vogel J (2015) Accelerating Discovery and Functional Analysis of Small RNAs with New Technologies. Annu Rev Genet 49: 367–394. https://doi.org/10.1146/annurev-genet-112414-054804 PMID: 26473381

44. Camarena L, Bruno V, Euskirchen G, Poggio S, Snyder M (2010) Molecular mechanisms of ethanol-induced pathogenesis revealed by RNA-sequencing. PLoS Pathog 6: e1000834. https://doi.org/10.1371/journal.ppat.1000834 PMID: 20368969

45. Lee HJ, Hong SH (2012) Analysis of microRNA-size, small RNAs in Streptococcus mutans by deep sequencing. FEMS Microbiol Lett 326: 131–136. https://doi.org/10.1111/j.1574-6968.2011.02441.x PMID: 22092283

46. Dotsch A, Eckweiler D, Schniederjans M, Zimmermann A, Jensen V, et al. (2012) The Pseudomonas aeruginosa transcriptome in planktonic cultures and static biofilms using RNA sequencing. PLoS One 7: e31092. https://doi.org/10.1371/journal.pone.0031092 PMID: 22319605

47. Suar M, Ryan D (2015) Small RNA in the acid tolerance response of Salmonella and their role in virulence. Virulence 6: 105–106. https://doi.org/10.4161/21505594.2014.988543 PMID: 25853732

48. Zeng Q, Sundin GW (2014) Genome-wide identification of Hfq-regulated small RNAs in the fire blight pathogen Erwinia amylovora discovered small RNAs with virulence regulatory function. BMC Genomics 15: 414. https://doi.org/10.1186/1471-2164-15-414 PMID: 24885615

49. Yan Y, Su S, Meng X, Ji X, Qu Y, et al. (2013) Determination of sRNA expressions by RNA-seq in Yersinia pestis grown in vitro and during infection. PLoS One 8: e74495. https://doi.org/10.1371/journal.pone.0074495 PMID: 24040259
50. Bardill JP, Hammer BK (2012) Non-coding sRNAs regulate virulence in the bacterial pathogen Vibrio cholerae. RNA Biol 9: 392–401. https://doi.org/10.4161/rna.19975 PMID: 22546941

51. Sievers S, Sternkopf Lillebaek EM, Jacobsen K, Lund A, Mollerup MS, et al. (2014) A multicopy sRNA of Listeria monocytogenes regulates expression of the virulence adhesin LapB. Nucleic Acids Res 42: 9383–9398. https://doi.org/10.1093/nar/gku630 PMID: 25034691

52. Wenner N, Maes A, Cotado-Sampayo M, Lapouge K (2014) NrsZ: a novel, processed, nitrogen-dependent, small non-coding RNA that regulates Pseudomonas aeruginosa PAO1 virulence. Environ Microbiol 16: 1053–1068. https://doi.org/10.1111/1462-2920.12272 PMID: 24308329

53. Wang L, Yang G, Qi L, Li X, Jia L, et al. (2016) A Novel Small RNA Regulates Tolerance and Virulence in Shigella flexneri by Responding to Acidic Environmental Changes. Front Cell Infect Microbiol 6: 24. https://doi.org/10.3389/fcimb.2016.00024 PMID: 27014636

54. Olson PD, Kuechenmeister LJ, Anderson KL, Daily S, Beenken KE, et al. (2011) Small molecule inhibitors of Staphylococcus aureus RnpA alter cellular mRNA turnover, exhibit antimicrobial activity, and attenuate pathogenesis. PLoS Pathog 7: e1001287. https://doi.org/10.1371/journal.ppat.1001287 PMID: 21347352

55. Amarasinghe JJ, Connell TD, Scannapieco FA, Haase EM (2012) Novel iron-regulated and Fur-regulated small regulatory RNAs in Aggregatibacter actinomycetemcomitans. Mol Oral Microbiol 27: 327–349. https://doi.org/10.1111/j.2041-1014.2012.00645.x PMID: 22958383

56. Sharma R, Arya S, Patil SD, Sharma A, Jain PK, et al. (2014) Identification of novel regulatory small RNAs in Acinetobacter baumannii. PLoS One 9: e93833. https://doi.org/10.1371/journal.pone.0093833 PMID: 24705412

57. Tormo MA, Marti M, Valle J, Manna AC, Cheung AL, et al. (2005) SarA is an essential positive regulator of Staphylococcus epidermidis biofilm development. J Bacteriol 187: 2348–2356. https://doi.org/10.1128/JB.187.7.2348-2356.2005 PMID: 15774878

58. Aranda J, Poza M, Pardo BG, Rumbo S, Rumbo C, et al. (2010) A rapid and simple method for constructing stable mutants of Acinetobacter baumannii. BMC Microbiol 10: 279. https://doi.org/10.1186/1471-2180-10-279 PMID: 21062436

59. Hamad MA, Zajdowicz SL, Holmes RK, Voskuil MI (2009) An allelic exchange system for compliant genetic manipulation of the select agents Burkholderia pseudomallei and Burkholderia mallei. Gene 430: 123–131. https://doi.org/10.1016/j.gene.2008.10.011 PMID: 19010402

60. Alvarez-Fraga L, Perez A, Rumbo-Feal S, Merino M, Vallejo JA, et al. (2016) Analysis of the role of the LH92_11085 gene of a biofilm hyper-producing Acinetobacter baumannii strain on biofilm formation and attachment to eukaryotic cells. Virulence: 0.

61. Rumbo-Feal S, Pérez A, Ramelot TA, Álvarez-Fraga L, Vallejo JA, et al. (2017) Contribution of the A. baumannii A1S_0114 Gene to the Interaction with Eukaryotic Cells and Virulence. Frontiers in Cellular and Infection Microbiology 7.

62. Pérez A, Merino M, Rumbo-Feal S, Álvarez-Fraga L, Vallejo JA, et al. (2016) The FhaB/FhaC two-partner secretion system is involved in adhesion of Acinetobacter baumannii AbH12O-A2 strain. Virulence: 1–16.

63. Weiss A, Broach WH, Lee MC, Shaw LN (2016) Towards the complete small RNome of Acinetobacter baumannii. Microb Genom 2: e000045. https://doi.org/10.1099/mgen.0.000045 PMID: 28348845

64. Gomez-Lozano M, Marvig RL, Molina-Santiago C, Tribelli PM, Ramos JL, et al. (2015) Diversity of small RNAs expressed in Pseudomonas species. Environ Microbiol Rep 7: 227–236. https://doi.org/10.1111/1758-2229.12233 PMID: 25394275

65. Ferrara S, Brugnoli M, De Bonis A, Righetti F, Delvillani F, et al. (2012) Comparative profiling of Pseudomonas aeruginosa strains reveals differential expression of novel unique and conserved small RNAs. PLoS One 7: e36553. https://doi.org/10.1371/journal.pone.0036553 PMID: 22590564

66. Gomez-Lozano M, Marvig RL, Molin S, Long KS (2012) Genome-wide identification of novel small RNAs in Pseudomonas aeruginosa. Environ Microbiol 14: 2006–2016. https://doi.org/10.1111/j.1462-2920.2012.02759.x PMID: 22533370

67. Martinez LC, Vadyvaloo V (2014) Mechanisms of post-transcriptional gene regulation in bacterial biofilms. Front Cell Infect Microbiol 4: 38. https://doi.org/10.3389/fcimb.2014.00038 PMID: 24724055

68. Ghaz-Jahanian MA, Khodaparastan F, Berenjian A, Jafarizadeh-Malmiri H (2013) Influence of small RNAs on biofilm formation process in bacteria. Mol Biotechnol 55: 288–297. https://doi.org/10.1007/s12033-013-9700-6 PMID: 24062263

69. Mika F, Hengge R (2013) Small Regulatory RNAs in the Control of Motility and Biofilm Formation in E. coli and Salmonella. Int J Mol Sci 14: 4560–4579. https://doi.org/10.3390/ijms14034560 PMID: 23443158
70. Majdalani N, Hernandez D, Gottesman S (2002) Regulation and mode of action of the second small RNA activator of RpoS translation, RprA. Mol Microbiol 46: 813–826. PMID: 12410838

71. Thomason MK, Fontaine F, De Lay N, Storz G (2012) A small RNA that regulates motility and biofilm formation in response to changes in nutrient availability in Escherichia coli. Mol Microbiol 84: 17–35. https://doi.org/10.1111/j.1365-2958.2012.07965.x PMID: 22289118

72. Jorgensen MG, Nielsen JS, Boysen A, Franch T, Moller-Jensen J, et al. (2012) Small regulatory RNAs control the multi-cellular adhesive lifestyle of Escherichia coli. Mol Microbiol 84: 36–50. https://doi.org/10.1111/j.1365-2958.2012.07976.x PMID: 22250746

73. Holmqvist E, Reimegard J, Sterk M, Grantcharova N, Romling U, et al. (2010) Two antisense RNAs target the transcriptional regulator CsgD to inhibit curli synthesis. Embo J 29: 1840–1850. https://doi.org/10.1038/emboj.2010.73 PMID: 20407422

74. Bak G, Lee J, Suk S, Kim D, Young Lee J, et al. (2015) Identification of novel sRNAs involved in biofilm formation, motility, and fimbriae formation in Escherichia coli. Sci Rep 5: 15287. https://doi.org/10.1038/srep15287 PMID: 26469694

75. Chambers JR, Sauer K (2013) Small RNAs and their role in biofilm formation. Trends Microbiol 21: 39–49. https://doi.org/10.1016/j.tim.2012.10.008 PMID: 23178000

76. Weilbacher T, Suzuki K, Dubey AK, Wang X, Gudapaty S, et al. (2003) A novel sRNA component of the carbon storage regulatory system of Escherichia coli. Mol Microbiol 48: 657–670. PMID: 12694612

77. Liu Z, Trevino J, Ramirez-Pena E, Sumby P (2012) The small regulatory RNA FasX controls pilus expression and adherence in the human bacterial pathogen group A Streptococcus. Mol Microbiol 86: 140–154. https://doi.org/10.1111/j.1365-2958.2012.08178.x PMID: 22882718

78. Eijkelkamp BA, Stroeher UH, Hassan KA, Papadimitrious MS, Paulsen IT, et al. (2011) Adherence and motility characteristics of clinical Acinetobacter baumannii isolates. FEMS Microbiol Lett 323: 44–51. https://doi.org/10.1111/j.1574-6968.2011.02362.x PMID: 22092679

79. de Breij A, Dijkshoorn L, Lagendijk E, van der Meer J, Koster A, et al. (2010) Do biofilm formation and interactions with human cells explain the clinical success of Acinetobacter baumannii? PLoS One 5: e10732. https://doi.org/10.1371/journal.pone.0010732 PMID: 20505779

80. McConnell MJ, Actis L, Pachon J (2013) Acinetobacter baumannii: human infections, factors contributing to pathogenesis and animal models. FEMS Microbiol Rev 37: 130–155. https://doi.org/10.1111/j.1574-6976.2012.00344.x PMID: 22568581