# Whole Transcriptome Analysis of Acinetobacter baumannii Assessed by RNA-Sequencing Reveals Different mRNA Expression Profiles in Biofilm Compared to Planktonic Cells

Soraya Rumbo-Feal, Manuel J Gómez, Carmen Gayoso, Laura Alvarez-Fraga, María P Cabral, Ana M Aransay, Naiara Rodríguez-Ezpeleta, Ane Fullaondo, Jaione Valle, María Tomás, et al.

► To cite this version:

Soraya Rumbo-Feal, Manuel J Gómez, Carmen Gayoso, Laura Alvarez-Fraga, María P Cabral, et al. Whole Transcriptome Analysis of Acinetobacter baumannii Assessed by RNA-Sequencing Reveals Different mRNA Expression Profiles in Biofilm Compared to Planktonic Cells. PLoS ONE, 2013, 8 (8), pp.e72968. 10.1371/journal.pone.0072968. [https://doi.org/10.1371/journal.pone.0072968](https://doi.org/10.1371/journal.pone.0072968).

HAL Id: hal-04939843

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

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>

Distributed under a Creative Commons Attribution 4.0 International License
# Whole Transcriptome Analysis of *Acinetobacter baumannii* Assessed by RNA-Sequencing Reveals Different mRNA Expression Profiles in Biofilm Compared to Planktonic Cells

Soraya Rumbo-Feal¹*, Manuel J. Gómez²*, Carmen Gayoso¹*, Laura Álvarez-Fraga¹, María P. Cabral¹, Ana M. Aransay³, Naiara Rodríguez-Ezpeleta³,⁴, Ane Fullaondo³, Jaione Valle⁵, María Tomás¹, Germán Bou¹*, Margarita Poza¹*

1 Department of Microbiology, Biomedical Research Institute, University Hospital, A Coruña, Spain, 2 Department of Molecular Evolution, Center for Astrobiology, INTA-CSIC, Madrid, Spain, 3 Genome Analysis Platform, CIC bioGUNE & CIBERehd, Derio, Spain, 4 Marine Research Division, AZTI, Tecnalia, Sukarrieta, Spain, 5 Department of Microbial biofilms, Agrobiotechnology Institute, Navarra, Spain

## Abstract

*Acinetobacter baumannii* has emerged as a dangerous opportunistic pathogen, with many strains able to form biofilms and thus cause persistent infections. The aim of the present study was to use high-throughput sequencing techniques to establish complete transcriptome profiles of planktonic (free-living) and sessile (biofilm) forms of *A. baumannii* ATCC 17978 and thereby identify differences in their gene expression patterns. Collections of mRNA from planktonic (both exponential and stationary phase cultures) and sessile (biofilm) cells were sequenced. Six mRNA libraries were prepared following the mRNA-Seq protocols from Illumina. Reads were obtained in a HiScanSQ platform and mapped against the complete genome to describe the complete mRNA transcriptomes of planktonic and sessile cells. The results showed that the gene expression pattern of *A. baumannii* biofilm cells was distinct from that of planktonic cells, including 1621 genes over-expressed in biofilms relative to stationary phase cells and 55 genes expressed only in biofilms. These differences suggested important changes in amino acid and fatty acid metabolism, motility, active transport, DNA-methylation, iron acquisition, transcriptional regulation, and quorum sensing, among other processes. Disruption or deletion of five of these genes caused a significant decrease in biofilm formation ability in the corresponding mutant strains. Among the genes over-expressed in biofilm cells were those in an operon involved in quorum sensing. One of them, encoding an acyl carrier protein, was shown to be involved in biofilm formation as demonstrated by the significant decrease in biofilm formation by the corresponding knockout strain. The present work serves as a basis for future studies examining the complex network systems that regulate bacterial biofilm formation and maintenance.

**Citation:** Rumbo-Feal S, Gómez MJ, Gayoso C, Álvarez-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(8): e72968. doi: 10.1371/journal.pone.0072968

**Editor:** Gunnar F Kaufmann, The Scripps Research Institute and Sorrento Therapeutics, Inc., United States of America

**Received:** March 13, 2013; **Accepted:** July 15, 2013; **Published:** August 30, 2013

**Copyright:** © 2013 Rumbo-Feal et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work was supported by grants from Instituto de Salud Carlos III - Ministerio Economía y Competitividad Unión Europea - Fondo Europeo de Desarrollo Regional (FEDER) PI11/01034 to MP and PI10/00056 to MT, by Ministerio de Economía y Competitividad, Instituto de Salud Carlos III - co-financed by European Development Regional Fund "A way to achieve Europe" ERDF, Spanish Network for the Research in Infectious Diseases (REIPI RD12/0015) and Fondo de Investigaciones Sanitarias (PI081638, PI12/00552) to GB and by the project Consolider INGENIO CSD2007-0005 and INTA-CSIC (Madrid, Spain) to MJG. CIC bioGUNE support was provided from The Department of Industry, Tourism and Trade of the Government of the Autonomous Community of the Basque Country (Etortek Research Programs 2010/ 2012) and from the Innovation Technology Department of the Bizkaia County. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

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

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

© These authors contributed equally to this work.

## Introduction

*Acinetobacter baumannii* are non-fermentative, oxidase negative, non-flagellated gram-negative bacilli. Although this species is a normal inhabitant of the human skin flora, intestinal

tract, and respiratory system, it has been shown to cause severe disease, including bacteremia and pneumonia, especially in patients hospitalized in intensive care units and reanimation wards [1–5]. Consequently, *A. baumannii* was recently listed as one of the six most dangerous opportunistic
pathogens [6,7]. Its high genetic plasticity allows it to rapidly adapt to stressful or otherwise unfavorable conditions by acquiring mutations, plasmids, or transposable elements. Moreover, *A. baumannii* species exhibit a remarkable ability to develop antibiotic resistance, which may quickly evolve into a multiresistant pattern following the acquisition of different resistance mechanisms, including β-lactamases, efflux pumps, porins, penicillin-binding proteins (PBPs), and methylase enzymes [1,2,8–17].

Bacteria often adopt sessile lifestyles in the form of matrix-enclosed habitats referred to as biofilms [18]. These are dynamic structures in which transitions between planktonic and sessile modes of growth occur in response to different environmental signals. The bacterial species inhabiting biofilms differ physiologically and behaviorally from their free-living counterparts [19]. Importantly, the structural characteristics of biofilms make them resistant to most antibiotics and host defenses [19–24]. This persistence provides a source of recurrent infections. In the case of *A. baumannii*, the infection of mucous surfaces and bacterial contamination of medical devices, such as intravascular catheters or endotracheal intubation, may result in biofilm formation, increasing the risk of bloodstream and respiratory infections [9].

An understanding of the ability of *A. baumannii* to form biofilms that adhere to and persist on a broad range of surfaces may offer the key to revealing its pathogenic mechanisms. Biofilm formation in *A. baumannii* has been shown to involve several regulatory processes, including those based on the sensing of bacterial cell density, the presence of different nutrients, and the concentration of free cations available to bacterial cells. Some of these extracellular signals may be sensed by two-component regulatory systems such as BfmRS. This transcriptional regulatory system activates expression of the usher-chaperone assembly apparatus responsible for the production of pili, which are needed for cell attachment and subsequent biofilm formation on polystyrene surfaces [25–27]. Recently, three new two-component sensor/regulator systems involved in biofilm formation in *P. aeruginosa* were identified [28]. In *A. baumannii*, a homolog of the biofilm-associated protein (Bap) of *Staphylococcus aureus* has been described [29] and the involvement of the membrane protein OmpA in the development of solid biofilms on abiotic surfaces and in virulence was demonstrated [30,31].

Several studies have reported changes in amino acid
metabolism during biofilm development [32–35]. In their
analyses of planktonic and sessile cells from biofilms of *A.
baumannii* ATCC 17978, Cabral et al. [34] found differences in
their proteomic profiles. In general, processes involved in
bacterial adhesion and the formation and development of
biofilms engage complex regulatory networks that coordinate
the temporal expression of genes related to adhesion, motility,
and the synthesis of matrix components. Although the ability of
*A. baumannii* to form biofilms on abiotic surfaces contributes to
the unique survival pattern of this pathogen in hospital settings,
little is known about the mechanisms that promote and support
biofilm formation. However, this knowledge is essential to the
identification of new therapeutic targets and thus to the design

of drugs effective against persistent diseases caused by multi-
resistant biofilm-forming clones of *A. baumannii*.

Microarray technology has been used to obtain the complete transcriptional profiles of different microorganisms and offers an approach to studying biofilm formation [36,37]. For example, Whiteley et al. [38] found significant differences in gene expression between sessile and planktonic cells in *Pseudomonas aeruginosa*. Moreno-Paz et al. [19] demonstrated different profiles in cells of the iron-oxidizing bacteria *Leptospirillium* grown in biofilm vs. planktonic modes. In *A. baumannii*, Hood et al. described the distinct transcriptional profile of the bacterium in response to NaCl [39] while Eijkelkamp et al. [40] were able to analyze its transcriptome in cultures grown under iron-limiting conditions (which prevents biofilm formation), reporting major transcriptional changes mostly related to iron acquisition but also to motility processes.

Among the more recent techniques used to analyze the
genome-wide RNA profiles of a number of organisms is deep
sequencing, using the platforms 454 GS FLX (Roche),
Genome Analyzer or HiSeq (Illumina Inc.), and ABI SOLID (Life
Technologies). These are open platforms not limited to the
study of previously known genes, and they are sensitive as well
as fast [41–48]. RNA sequencing using the Illumina system has
developed as an extremely informative technique for the study
of transcriptional profiles of microbes [45,49,50].

The aim of the present study was to use bacterial mRNA and the Illumina RNA-sequencing technologies to gain insight into the mechanisms behind the remarkable ability of *A. baumannii* to form biofilms. We therefore obtained whole transcriptomes from planktonic and biofilm cells of *A. baumannii* strain ATCC 17978 and then compared them for differences in their gene expression profiles.

**Materials and Methods**

Strains and culture conditions

A. baumannii ATCC 17978 was routinely grown in Mueller-Hinton (MH) broth. *E. coli* TG1, used for cloning procedures, was grown in Luria-Bertani (LB) broth. Agar was added to a final concentration of 2% when necessary. All strains were grown at 37 °C with shaking (180 rpm), and stored at -80 °C in LB broth containing 10% glycerol. Kanamycin (50 µg/mL) and rifampicin (50 µg/mL) were from Sigma-Aldrich (St. Louis, MO) and were added to select transformant strains. Cultures of planktonic cells originated from a single colony of *A. baumannii* strain ATCC 17978 isolated in MH agar and then grown in 5 mL of MH broth overnight as described above. The resulting culture was diluted 100-fold in 500 mL of MH broth in 1-L flasks and again grown as described above, measuring the optical density at 600 nm (OD600nm) every 30 min. Cells were harvested during the exponential (OD600nm = 0.4) and late stationary phases (OD600nm = 2.0) of growth, 48 h after inoculation. Planktonic and sessile cells (obtained as described below) were resuspended in RNA Later reagent (Sigma-Aldrich), frozen using liquid nitrogen, and stored at -80 °C.
## Biofilm generation in Pyrex plates

A. baumannii ATCC 17978 biofilms were obtained in the Fermentation Laboratory of the Agrobiotechnology Institute (Navarra, Spain). A sample from an overnight culture of A. baumannii grown in MH broth was used to inoculate 60-mL microfermentors (Institute Pasteur, Paris, France), which were then maintained at 37 °C for 24 h. The bacterium was grown in MH broth medium under a continuous-flow culture system and continuous aeration consisting of 40 mL of compressed, sterile air/h. Submerged Pyrex slides served as the growth substratum. Biofilms that formed on the Pyrex slides were removed with a cell scraper and frozen in liquid nitrogen at -80 °C.

## Isolation of mRNA

Three samples, corresponding to exponential and stationary phase cells and sessile cells from biofilms, were reduced to powder under liquid nitrogen while grinding using a mortar and pestle. Total RNA was then isolated using the mirVana miRNA isolation kit (Ambion) following the manufacturer's protocols. Ten µg of each total RNA was further processed by removing 23S and 16S rRNAs using the MICROBExpress bacterial mRNA enrichment kit (Ambion). The rRNA-depleted samples (free of 16S and 23S rRNA) of exponentially growing, stationary phase, and biofilm cells were treated with DNAse I (Invitrogen), purified using phenol-chloroform, and concentrated by ethanol precipitation. Final concentrations and purity grades of the samples were determined using a NanoDrop ND-1000 (Thermo Scientific) and a BIOANALYZER 2100 (Agilent Technologies Inc., Germany).

## Transcription assays

A cDNA synthesis kit (Roche) was used to obtain double-stranded cDNA (ds-cDNA), following the manufacturer's instructions. The rRNA-depleted samples together with 5'-phosphorylated degenerated hexamers and the AMV reverse transcriptase (both from Roche) were used to obtain the first cDNA strand. The second cDNA strand was then generated and treated with RNase. The ds-cDNA products were purified using the High Pure PCR purification kit (Roche). The samples were further quantified using a Nanodrop ND-1000 (Thermo Scientific) and a BIOANALYZER 2100 (Agilent Technologies Inc., Germany). The ds-cDNAs were then used in subsequent steps of the study.

## Deep-sequencing procedures

To characterize the complete transcriptomes of the studied samples, mRNA libraries from three cellular conditions (exponential and stationary phase planktonic cells and sessile cells from biofilms) were prepared following the Truseq RNA sample preparation protocols from Illumina Inc. at CIC bioGUNE's genome analysis platform (Derio, Spain). Two biological replicates were studied for each sample.

## Read processing and comparisons of gene expression profiles

Fifty nucleotide reads from each mRNA library were obtained using HiScanSQ (Illumina Inc., CIC bioGUNE, Bilbao, Spain). Short reads were aligned against the complete genome of A. baumannii ATCC 17978 and plasmids pAB1 and pAB2 (GenBank accession codes: NC_009085.1, NC_009083.1 and NC_009084.1, respectively) using Bowtie [51], allowing a maximum of three mismatches within the first 50 bases. Reads were annotated with the R Bioconductor Genominator package and differences in expression levels estimated with the R DESeq package [52]. DESeq performs a count normalization to control the variation in the number of reads sequenced across samples. After normalization, fold changes and their significance (p values), indicating differential expression, were determined after a negative binomial distribution. Those mRNAs with p values (adjusted for a false discovery rate of 0.1%) < 0.001 were considered to be differentially expressed. Raw sequences were deposited at the NCBI Sequence Read Archive, under Bioproject accession number PRJNA191863 (experiment accessions codes SRX263965, SRX263966, SRX263968 and SRX263969 to SRX263977). Blast2GO [53] was used for the functional re-annotation of genes, the mapping of gene ontology terms, and the description of biological processes, molecular functions, cellular components, and metabolic pathways associated with the biofilm expression profiles.

## Quantitative biofilm assay

Biofilm formation was quantified using the procedure described by Tendolkar et al. [54], with slight modifications. A colony of A. baumannii was grown on MH agar medium for 18 h at 37 °C and used to inoculate 25 mL of liquid MH medium, supplemented with 50 µg kanamycin /mL when necessary. The culture was maintained overnight at 37 °C and 180 rpm. Cells were harvested by centrifugation (3500 g, 10 min), washed three times with 0.9% NaCl, and resuspended in fresh liquid medium without antibiotic. From this suspension, 100 µL (containing $10^8$ CFU) were dispensed into each well of a 96-well flat-bottom polypropylene microtiter plate containing MH medium and then incubated at 37 °C for 24–48 h. Next, the cells were stained with 25 µL of a 1% w/v crystal violet solution for 15 min at room temperature, washed twice with sterile 0.9% w/v NaCl, solubilized with 200 µL of a 4: 1 v/v mixture of ethanol and acetone, and finally quantified at 570 nm. All biofilm assays were performed with at least six replicates for each strain. ANOVA tests were used to evaluate the statistical significance of the measured differences.

## Gene disruption

Plasmids were inserted into the target genes as previously described [55], with slight modifications. Briefly, kanamycin- and zeocin-resistant plasmid pCR-BluntII-TOPO (Invitrogen), unable to replicate in A. baumannii, was used as a suicide vector. An internal fragment (~ 500 bp) of the target gene was PCR-amplified using the primers listed in Table 1 and genomic DNA from A. baumannii ATCC 17978 as template. The PCR products were cloned into the pCR-BluntII-TOPO vector and
the recombinant plasmids (0.1 µg) were introduced into kanamycin- and zeocin-susceptible *A. baumannii* ATCC 17978 by electroporation. Mutants were selected on kanamycin-containing plates. Inactivation of the target gene by insertion of the plasmid via single-crossover recombination was confirmed by sequencing the PCR-amplified products using the primers listed in Table 1.

## Construction of knockout strains

Knockout strains were constructed using the plasmid pMo130 (Genbank accession code EU862243), containing the genes *xylE*, *sacB* and a kanamycin resistance marker, as a suicide vector [56]. Briefly, 851–932 bp upstream and downstream of the *A. baumannii* ATCC 17978 (Genbank accession code NC_009085.1) gene of interest were cloned into the pMo130 vector using the primers listed in Table 1. The resulting plasmid was used to transform *A. baumannii* 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 cathecol-based method described by Hamad et al. [56]. Bright yellow and kanamycin-resistant colonies were grown overnight in LB supplemented with 15% sucrose and then plated on the same agar medium. The second crossover event was confirmed by PCR using the primers listed in Table 1. Quantitative biofilm assays were used to determine the phenotype of the mutants.

## Complementation of stable knockout mutant

To complement the stable knockout mutant, the target gene was amplified from *A. baumannii* ATCC 17978 genomic DNA using the primers listed in Table 1 and then cloned into the XbaI restriction site of the pET-RA plasmid under the control of the β-lactamase CXT-M-14 gene promoter, as described by Aranda et al. [57]. The new construction was used to transform the mutant strain. Transformants were selected on rifampicin- and kanamycin-containing plates and confirmed by PCR using the primers listed in Table 1. The mutant strain containing the pET-RA plasmid was used as the control.

## Real-time RT-PCR

Real-time reverse transcription-PCR (RT-PCR) was carried out to determine the expression levels of a collection of genes using Taqman probes (TIB Mol Biol) listed in Table 1. In all cases, the expression levels were standardized relative to the transcription levels of the housekeeping gene *gyrB*. The primers used were those listed in Table 1. Total RNA was isolated from exponentially growing and stationary phase cultures and from the biofilms using the High Pure RNA isolation kit (Roche, Germany) and then treated with RNase-free DNase I (Invitrogen Corporation, CA). The samples were further purified using the RNeasy MinElute Cleanup kit (Qiagen, Germany). For qRT-PCR, the LightCycler 480 RNA Master hydrolysis probes kit and a LightCycler 480 RNA instrument (both from Roche, Germany) were used together with the following protocol: initial incubation of 65 °C, 3 min, followed by a denaturation step at 95 °C for 30 s, 45 cycles at 95 °C, 15 s and 60 °C, 45 s, and a final elongation step at 40 °C,

<table><caption>**Table 1.** Oligonucleotides and probes used in the present work.</caption><thead><tr><th>Primer/Probe name</th><th>Sequence</th><th>Use in the present study</th></tr></thead><tbody><tr><td>0114intF</td><td>actggagcgcaatcattcgt</td><td>Disruption of gene A1S_0114</td></tr><tr><td>0114intR</td><td>atgaagcaactccctgctgc</td><td>Disruption of gene A1S_0114</td></tr><tr><td>0114extF</td><td>caaggagttgaaacgat</td><td>Confirm disruption of gene A1S_0114</td></tr><tr><td>0114extR</td><td>ctcgcagcaatagacca</td><td>Confirm disruption of gene A1S_0114</td></tr><tr><td>0302intF</td><td>cggaagcagtggtaaacttgc</td><td>Disruption of gene A1S_0302</td></tr><tr><td>0302intR</td><td>tggtgaaaacacgcgagagc</td><td>Disruption of gene A1S_0302</td></tr><tr><td>0302extF</td><td>acaccaactatttccgtg</td><td>Confirm disruption of gene A1S_0302</td></tr><tr><td>0302extR</td><td>cccaaaatcagtcacctc</td><td>Confirm disruption of gene A1S_0302</td></tr><tr><td>1507intF</td><td>ccacaccaactccgtttgct</td><td>Disruption of gene A1S_1507</td></tr><tr><td>1507intR</td><td>acttgcaaccgtgccaatga</td><td>Disruption of gene A1S_1507</td></tr><tr><td>1507extF</td><td>tgtgtgtgatcattttgac</td><td>Confirm disruption of gene A1S_1507</td></tr><tr><td>1507extR</td><td>aagagcggtttactcatc</td><td>Confirm disruption of gene A1S_1507</td></tr><tr><td>3168intF</td><td>atctcgagcagcttgtgcag</td><td>Disruption of gene A1S_3168</td></tr><tr><td>3168intR</td><td>attaagccgtggtgcaggtg</td><td>Disruption of gene A1S_3168</td></tr><tr><td>3168extF</td><td>actcttattgccaaacc</td><td>Confirm disruption of gene A1S_3168</td></tr><tr><td>3168extR</td><td>cttgcttaatgatggagg</td><td>Confirm disruption of gene A1S_3168</td></tr><tr><td>2042intF</td><td>tgactggatttacacagaaga</td><td>Disruption of gene A1S_2042</td></tr><tr><td>2042intR</td><td>tgttccatcattaataactcc</td><td>Disruption of gene A1S_2042</td></tr><tr><td>2042extF</td><td>ccagagcactagccttaa</td><td>Confirm disruption of gene A1S_2042</td></tr><tr><td>2042extR</td><td>ttgagtgagtgcagctaa</td><td>Confirm disruption of gene A1S_2042</td></tr><tr><td>0114UpFNotI</td><td>cccgcggccgcgggttggtagctgagcaactc</td><td>Construction of stable knockout strain ΔA1S_0114</td></tr><tr><td>0114UpRBamHI</td><td>gggggatccccgcgggtaatctcctttttaacc</td><td>Construction of stable knockout strain ΔA1S_0114</td></tr><tr><td>0114DownFBamHI</td><td>cccggtatccgggacaaccttgtcacgactagaa</td><td>Construction of stable knockout strain ΔA1S_0114</td></tr><tr><td>0114DownRXbal</td><td>gggtctagacccttcaagtcgacctgctacg</td><td>Construction of stable knockout strain ΔA1S_0114</td></tr></tbody></table>
<table><caption>Table 1 (continued).</caption><thead><tr><th>Primer/Probe name</th><th>Sequence</th><th>Use in the present study</th></tr></thead><tbody><tr><td>pMo130 site2 F</td><td>attcatgaccgtgctgac</td><td>Confirm construction of stable knockout strain ΔA1S_0114</td></tr><tr><td>pMo130 site2 R</td><td>cttgtctgtaagcggatg</td><td>Confirm construction of stable knockout strain ΔA1S_0114</td></tr><tr><td>0114XbaIF</td><td>ccctctagaggggtattcgctcgtattgctg</td><td>Cloning of the gene A1S_0114 into the pET-RA plasmid for complementation of the stable knockout strain ΔA1S_0114</td></tr><tr><td>0114XbaIR</td><td>ccctctagaggggactgggtgaccttcacatc</td><td>Cloning of the gene A1S_0114 into the pET-RA plasmid for complementation of the stable knockout strain ΔA1S_0114</td></tr><tr><td>pETRAF</td><td>ttcttcgtgaaatagtgtattttt</td><td>Confirm complementation of stable knockout strain ΔA1S_0114</td></tr><tr><td>pETRAR</td><td>ctgtttcatatgatctgggtatc</td><td>Confirm complementation of stable knockout strain ΔA1S_0114</td></tr><tr><td>A1S_0109F</td><td>caaacatcgaaatatccatcaatcgtc</td><td>qRT-PCR</td></tr><tr><td>A1S_0109R</td><td>cagccgtagatttttcaaatccg</td><td>qRT-PCR</td></tr><tr><td>A1S_0109 Taqman probe</td><td>cctctagcagtcaggctgtgtcatcacc</td><td>qRT-PCR</td></tr><tr><td>A1S_0112F</td><td>accagaagatgttggcctga</td><td>qRT-PCR</td></tr><tr><td>A1S_0112R</td><td>gagccgatcaaccccatca</td><td>qRT-PCR</td></tr><tr><td>A1S_0112 Taqman Probe</td><td>gctgcctg</td><td>qRT-PCR</td></tr><tr><td>A1S_0113F</td><td>tggctttaacaacgctgaaa</td><td>qRT-PCR</td></tr><tr><td>A1S_0113R</td><td>aacccctgaccttcttcacc</td><td>qRT-PCR</td></tr><tr><td>A1S_0113 Taqman Probe</td><td>tgccctga</td><td>qRT-PCR</td></tr><tr><td>A1S_0114F</td><td>gtagagcctgagacgattgatcca</td><td>qRT-PCR</td></tr><tr><td>A1S_0114R</td><td>gttggctcaagtcttaattttcgta</td><td>qRT-PCR</td></tr><tr><td>A1S_0114 Taqman Probe</td><td>ttctaaatccccagacacagacaaagcaa</td><td>qRT-PCR</td></tr><tr><td>A1S_0302F</td><td>gcaggtaaagcaataatatcgaaag</td><td>qRT-PCR</td></tr><tr><td>A1S_0302R</td><td>ttatcaactaaggagaagctagcaagt</td><td>qRT-PCR</td></tr><tr><td>A1S_0302 Taqman Probe</td><td>ggaagcag</td><td>qRT-PCR</td></tr><tr><td>A1S_1507F</td><td>acaccaactccgtttgcttt</td><td>qRT-PCR</td></tr><tr><td>A1S_1507R</td><td>ctgacacttcaaatagccaggtt</td><td>qRT-PCR</td></tr><tr><td>A1S_1507 Taqman Probe</td><td>tcagcagc</td><td>qRT-PCR</td></tr><tr><td>A1S_3168F</td><td>tcgcatctcgagcagctt</td><td>qRT-PCR</td></tr><tr><td>A1S_3168R</td><td>cgcagctggaatttttgctt</td><td>qRT-PCR</td></tr><tr><td>A1S_3168 Taqman Probe</td><td>cagccacc</td><td>qRT-PCR</td></tr></tbody></table>

<table><caption>Table 1 (continued).</caption><thead><tr><th>Primer/Probe name</th><th>Sequence</th><th>Use in the present study</th></tr></thead><tbody><tr><td>A1S_2042F</td><td>tggtatattgactggatttacacaga</td><td>qRT-PCR</td></tr><tr><td>A1S_2042R</td><td>catcattaataactccatcgagg</td><td>qRT-PCR</td></tr><tr><td>A1S_2042 Taqman Probe</td><td>tggctctatgagcttgtttttctatttt</td><td>qRT-PCR</td></tr><tr><td>gyrBF</td><td>tctctagtcatgagaagtgggtacatt</td><td>qRT-PCR</td></tr><tr><td>gyrBR</td><td>ggttatattcttcacggccaat</td><td>qRT-PCR</td></tr><tr><td>gyrB Taqman Probe</td><td>tggctgtg</td><td>qRT-PCR</td></tr></tbody></table>

doi: 10.1371/journal.pone.0072968.t001

30 s. All assays were performed in triplicate. The statistical significance of the determined differences was confirmed by ANOVA tests.

## Results

### Determination of the complete transcriptomes of planktonic and biofilm cells

The mRNA fractions purified from exponentially growing (Exp) and stationary-phase (Sta) cultures and from biofilms (Bio) of A. baumannii ATCC 17978 were analyzed to determine the respective gene expression level profiles and to identify differentially expressed genes. Six libraries, including two biological replicates *per* sample, were constructed (Exp 1, Exp 2, Sta 1, Sta 2, Bio 1, Bio 2) and paired-end sequenced using Illumina technology (50 bp x 2). Insert average sizes in the above mentioned libraries were 208, 240, 230, 239, 209 and 253 bp, respectively. Reads were aligned against the chromosome and plasmids of A. baumannii ATCC 17978. The number of reads that mapped against the genome is detailed in Table 2. Gene level read counts are shown in Figure S1, and MD plots and correlation between samples in Figure S2. The complete mRNA transcriptomic profiles of exponentially growing and stationary-phases cultures and from biofilm cells were obtained by Illumina procedures. Gene expression values are provided in the Supporting Information (Tables S1, S2, and S3). Table S1 shows the gene expression profile of cells obtained in exponential growth phase vs. stationary phase cultures. Table S2 and Table S3 show the gene expression profile of biofilm cells vs. exponentially growing and stationary phase cultures, respectively. Overall, the data confirmed the complete description of the whole transcriptome of each stage of growth. Approximately 97% of the genes described in the A. baumannii ATCC 17978 genome database (NC_009085.1) were transcribed using the method described herein.

### Different mRNA expression patterns of cells grown in exponential phase, stationary phase and in biofilms

The expression patterns of exponentially growing vs. stationary phase cells, stationary phase cells vs. biofilm cells, and exponentially growing cells vs. biofilm cells were compared to identify differentially expressed transcripts. Up-regulated and down-regulated genes were determined based on differences
<table><caption>**Table 2.** Total number of reads aligning with the regions of interest (coverage) of the six libraries constructed from the mRNA samples.</caption><thead><tr><th>Exp 1</th><th>Exp 2</th><th>Sta 1</th><th>Sta 1</th><th>Bio 1</th><th>Bio 2</th></tr></thead><tbody><tr><td>8189422</td><td>9497363</td><td>7707791</td><td>6317554</td><td>9084229</td><td>3111192</td></tr></tbody></table>

doi: 10.1371/journal.pone.0072968.t002

for which the *p* values were below 0.001. The results are shown in Supporting Information (Tables S4, S5, S6, S7, S8, S9). Although, as expected, many genes were constitutively expressed, the comparisons indicated the association of each cellular condition with a specific expression profile, with significant differences in the expression level of a large number of genes. Thus, in biofilm vs. stationary-phase cells, 31 genes were down-regulated and 35 up-regulated; in biofilm vs. exponentially growing cells 15 genes were down-regulated and 116 up-regulated, and in stationary-phase vs. exponentially growing cells 130 genes were up-regulated and 33 down-regulated in (*p* < 0.001 in all cases).

## The gene expression profile of biofilm cells

A comparison of gene expression levels in biofilms vs. stationary phase cells without applying any *p* value filter indicated that among the 1621 genes over-expressed in biofilms there were 408 genes whose expression was at least four-fold higher in sessile cells or completely inhibited in planktonic cells but with an expression level value of at least 2 in biofilm. With the aim of describing gene expression profile differences in terms of gene ontology, the complete proteome of A. baumannii strain ATCC1 7978 was re-annotated using Blast2GO. Biological processes, molecular functions, and cellular components associated with the set of 1621 up-regulated genes in biofilms, as determined using Blast2GO, are shown in Figure 1. The results showed that the largest group, made up of 129 genes, was involved in transcriptional regulation. Many genes were those involved in acyl carrier protein biosynthetic processes, amino acid metabolism, fatty acid metabolism, ion transport, carbohydrate biosynthesis, translation, transmembrane transport, and the stress response, among other biological processes. The cellular location of the majority of the proteins encoded by these 1621 genes was in most cases consistent with the proteins being integral to the inner membrane or members of a transcription factor complex. Fewer proteins were located in the outer membrane periplasmic space, in the cell outer membrane, or in a transcriptional repressor complex. Moreover, there were small groups of genes that encoded proteins associated with the peptidoglycan-based cell wall, the type II protein secretion system complex, the fatty acid synthase complex, or cell projection, among other cellular components. According to the molecular function ontology, most of the genes over-expressed in biofilms were related to transferase, hydrolase, and oxidoreductase activities and, to a lesser extent, to metal ion, ATP, coenzyme, or DNA binding activities, among other molecular functions.

When the comparison of gene expression levels in biofilms vs. exponentially growing and stationary phase planktonic cells was filtered by *p* value (< 0.001), similar results were obtained (Figure S3). A list of genes differentially expressed (*p* < 0.001) in the biofilm vs. exponential and stationary cultures is presented in Table 3. Among them, genes coding an acyl carrier protein, an allophanate hydrolase and the RND efflux pump AdeT (A1S_0114, A1S_1278 and A1S_1755, respectively) were only expressed in biofilms and were totally inhibited in planktonic cells. Genes corresponding to hypothetical proteins (A1S_0302, A1S_0644, A1S_1293, and A1S_2893), a transmembrane arsenate pump protein (A1S_1454), the CsuD, CsuC, and CsuA/B proteins (A1S_2214, A1S_2215, and A1S_2218), the BasD protein (A1S_2382), a ferric acinetobactin binding protein (A1S_2386), a sulfate transport protein (A1S_2534), and maleylacetoacetate isomerase (A1S_3415) were also highly expressed in biofilms but totally inhibited in stationary cells. Many genes involved in amino acid metabolism and transport (such as A1S_0115, A1S_0429, A1S_1357, A1S_3134, A1S_3185, A1S_3402, A1S_3404, A1S_3405, A1S_3406, A1S_3407, or A1S_3413), or related to iron acquisition and transport (A1S_0653, A1S_0742, A1S_0980, A1S_1631, A1S_1657, A1S_2385, or A1S_2390, encoding a ferrous iron transport protein, an iron-regulated protein, a ferric enterobactin receptor, an iron-binding protein, a siderophore biosynthesis protein, a ferric acinetobactin receptor, and an acinetobactin biosynthesis protein, respectively), transcriptional regulators (A1S_1377 and A1S_1687), or encoding efflux pumps (A1S_0009 and A1S_0538) were also up-regulated in biofilm vs. planktonic cells. A gene coding for a fimbrial protein (A1S_1507) was highly up-regulated and the outer membrane protein A (A1S_2840) was down-regulated in biofilm cells vs. either exponential or stationary phase planktonic cells. An operon containing a group of genes related to phenylacetate metabolism (with identifiers A1S_1335 to A1S_1340) was up-regulated in biofilm cells vs. exponential cells but down-regulated vs. stationary cells. A homoserine lactone synthase (A1S_0109) was over-expressed in biofilms vs. either growth form of planktonic cells, as was a group of seven genes (from A1S_0112 to A1S_0118). Among the latter, A1S_0114 was an extreme case because of its very high level of expression in the biofilm (ca. 127) and lack of detectable expression in planktonic cells.

## Genes only expressed in biofilm associated cells

Fifty-five genes were exclusively expressed in sessile cells (Table 4), including 12 genes assigned to uncharacterized proteins and nine encoding transcriptional regulators (A1S_0547, A1S_1256, A1S_1430, A1S_1763, A1S_1958, A1S_2042, A1S_2151, A1S_2208 and A1S_3255). Other genes in this group belonged to the Csu operon (A1S_2216 and A1S_2217), encoded a membrane protein (A1S_0595), or were related to iron acquisition systems (A1S_0945, A1S_1719, A1S_2380, and A1S_2388). Genes coding for a DNA polymerase (A1S_2015), a DNA helicase (A1S_1585), an extracellular nuclease (A1S_1198), and an endonuclease (A1S_2408) as well as two genes involved in DNA methylation
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Sequence distribution of the 1621 genes identified in the present work as up-regulated in biofilm vs. stationary phase cells. Genes involved in: A) biological processes, B) cellular components, and C) molecular functions. The results were filtered by the number of sequences (cutoff = 40, 5, and 80, respectively).</figcaption></figure>

doi: 10.1371/journal.pone.0072968.g001
<table><caption>**Table 3.** Differentially expressed genes in biofilm-associated cells vs. both exponentially growing and stationary-phase cells.</caption><thead><tr><th>Gene Id*</th><th>Gene description</th><th>Fold change Biofilm vs. exponential phase cells</th><th>Fold change Biofilm vs. stationary phase cells</th></tr></thead><tbody><tr><td>A1S_0004</td><td>DNA gyrase</td><td>0.44</td><td>0.19</td></tr><tr><td>A1S_0009</td><td>RND type efflux pump</td><td>2.57</td><td>4.84</td></tr><tr><td>A1S_0032</td><td>signal peptide</td><td>32.18</td><td>4.82</td></tr><tr><td>A1S_0073</td><td>2-methylisocitrate lyase</td><td>6.60</td><td>1.61</td></tr><tr><td>A1S_0087</td><td>short-chain dehydrogenase</td><td>2.57</td><td>6.37</td></tr><tr><td>A1S_0103</td><td>3-hydroxyisobutyrate dehydrogenase</td><td>61.82</td><td>3.81</td></tr><tr><td>A1S_0107</td><td>enoyl-CoA hydratase</td><td>5.51</td><td>5.97</td></tr><tr><td>A1S_0109</td><td>homoserine lactone synthase</td><td>60.22</td><td>16.74</td></tr><tr><td>A1S_0112</td><td>acyl-CoA synthetase/AMP-acid ligases II</td><td>75.17</td><td>53.14</td></tr><tr><td>A1S_0113</td><td>acyl-CoA dehydrogenase</td><td>135.15</td><td>42.91</td></tr><tr><td>A1S_0114</td><td>acyl carrier protein</td><td>from zero to 127.96**</td><td>from zero to 127.96**</td></tr><tr><td>A1S_0115</td><td>amino acid adenylation</td><td>151.37</td><td>32.08</td></tr><tr><td>A1S_0116</td><td>RND superfamily exporter</td><td>56.18</td><td>79.67</td></tr><tr><td>A1S_0117</td><td>hypothetical protein</td><td>23.97</td><td>8.73</td></tr><tr><td>A1S_0118</td><td>hypothetical protein</td><td>9.31</td><td>5.26</td></tr><tr><td>A1S_0151</td><td>F0F1 ATP synthase subunit B</td><td>1.90</td><td>4.77</td></tr><tr><td>A1S_0153</td><td>F0F1 ATP synthase subunit alpha</td><td>1.13</td><td>3.75</td></tr><tr><td>A1S_0154</td><td>F0F1 ATP synthase subunit gamma</td><td>1.65</td><td>4.77</td></tr><tr><td>A1S_0155</td><td>F0F1 ATP synthase subunit beta</td><td>1.06</td><td>4.55</td></tr><tr><td>A1S_0156</td><td>F0F1 ATP synthase subunit epsilon</td><td>1.00</td><td>4.54</td></tr><tr><td>A1S_0179</td><td>NADPH-dependent FMN reductase</td><td>0.01</td><td>0.71</td></tr><tr><td>A1S_0279</td><td>elongation factor Tu</td><td>1.18</td><td>3.08</td></tr><tr><td>A1S_0283</td><td>50S ribosomal protein L11</td><td>1.99</td><td>1.94</td></tr><tr><td>A1S_0285</td><td>50S ribosomal protein</td><td>2.31</td><td>5.91</td></tr><tr><td>A1S_0292</td><td>outer membrane protein W</td><td>0.53</td><td>0.08</td></tr><tr><td>A1S_0302</td><td>hypothetical protein</td><td>23.30</td><td>from zero to 27.09**</td></tr><tr><td>A1S_0360</td><td>30S ribosomal protein S15</td><td>0.74</td><td>0.14</td></tr><tr><td>A1S_0429</td><td>DAACS family glutamate:aspartate symporter</td><td>3.04</td><td>1.84</td></tr><tr><td>A1S_0445</td><td>hypothetical protein</td><td>0.53</td><td>0.14</td></tr><tr><td>A1S_0449</td><td>coniferyl aldehyde dehydrogenase</td><td>0.16</td><td>0.12</td></tr><tr><td>A1S_0481</td><td>phosphate acetyltransferase</td><td>3.92</td><td>2.49</td></tr><tr><td>A1S_0482</td><td>acetate kinase</td><td>3.31</td><td>0.75</td></tr><tr><td>A1S_0496</td><td>phosphatidylglycerophosphatase B</td><td>10.13</td><td>1.05</td></tr><tr><td>A1S_0528</td><td>preprotein translocase subunit SecB</td><td>0.67</td><td>0.14</td></tr><tr><td>A1S_0538</td><td>RND efflux transporter</td><td>6.59</td><td>11.03</td></tr><tr><td>A1S_0570</td><td>hypothetical protein</td><td>0.70</td><td>0.12</td></tr><tr><td>A1S_0591</td><td>acyl-CoA synthetase</td><td>6.17</td><td>0.74</td></tr><tr><td>A1S_0628</td><td>transposase</td><td>4.32</td><td>2.80</td></tr><tr><td>1S_0644</td><td>hypothetical protein</td><td>18.27</td><td>from zero to 33.44**</td></tr><tr><td>A1S_0653</td><td>ferrous iron transport protein B</td><td>4.46</td><td>1.58</td></tr><tr><td>A1S_0661</td><td>phage integrase family protein</td><td>0.59</td><td>0.07</td></tr><tr><td>A1S_0670</td><td>protein tyrosine phosphatase</td><td>0.65</td><td>0.06</td></tr><tr><td>A1S_0671</td><td>protein tyrosine phosphatase</td><td>0.53</td><td>0.04</td></tr></tbody></table>

<table><caption>**Table 3 (continued).**</caption><thead><tr><th>Gene Id*</th><th>Gene description</th><th>Fold change Biofilm vs. exponential phase cells</th><th>Fold change Biofilm vs. stationary phase cells</th></tr></thead><tbody><tr><td>A1S_0675</td><td>dihydropteroate synthase</td><td>0.88</td><td>0.09</td></tr><tr><td>A1S_0736</td><td>hypothetical protein</td><td>79.49</td><td>2.76</td></tr><tr><td>A1S_0737</td><td>methyltetrahydropteroyltriglutamate/homocysteine S-methyltransferase</td><td>25.32</td><td>2.26</td></tr><tr><td>A1S_0742</td><td>iron-regulated protein</td><td>2.50</td><td>1.74</td></tr><tr><td>A1S_0745</td><td>hypothetical protein</td><td>30.81</td><td>4.28</td></tr><tr><td>A1S_0869</td><td>elongation factor Tu</td><td>1.22</td><td>2.52</td></tr><tr><td>A1S_0884</td><td>outer membrane protein</td><td>0.50</td><td>4.35</td></tr><tr><td>A1S_0971</td><td>B12-dependent methionine synthase</td><td>0.05</td><td>0.07</td></tr><tr><td>A1S_0980</td><td>ferric enterobactin receptor precursor</td><td>4.38</td><td>4.01</td></tr><tr><td>A1S_1032</td><td>hypothetical protein</td><td>4.41</td><td>2.63</td></tr><tr><td>A1S_1077</td><td>hypothetical protein</td><td>9.41</td><td>2.43</td></tr><tr><td>A1S_1104</td><td>chlorogenate esterase</td><td>0.35</td><td>0.01</td></tr><tr><td>A1S_1266</td><td>hypothetical protein</td><td>1.09</td><td>11.88</td></tr><tr><td>A1S_1278</td><td>allophanate hydrolase subunit 2</td><td>from zero to 20.35**</td><td>from zero to 20.35**</td></tr><tr><td>A1S_1293</td><td>hypothetical protein</td><td>15.25</td><td>from zero to 30.54**</td></tr><tr><td>A1S_1316</td><td>major facilitator superfamily transporter cyanate permease</td><td>14.40</td><td>4.83</td></tr><tr><td>A1S_1319</td><td>hypothetical protein</td><td>22.56</td><td>50.37</td></tr><tr><td>A1S_1335</td><td>bifunctional aldehyde dehydrogenase/enoyl-CoA hydratase</td><td>21.33</td><td>0.22</td></tr><tr><td>A1S_1336</td><td>phenylacetate-CoA oxygenase subunit PaaA</td><td>93.43</td><td>0.28</td></tr><tr><td>A1S_1337</td><td>phenylacetate-CoA oxygenase subunit PaaB</td><td>22.63</td><td>0.45</td></tr><tr><td>A1S_1338</td><td>hypothetical protein</td><td>34.73</td><td>0.41</td></tr><tr><td>A1S_1339</td><td>phenylacetate-CoA oxygenase PaaJ subunit</td><td>196.37</td><td>0.78</td></tr><tr><td>A1S_1340</td><td>phenylacetate-CoA oxygenase/reductase PaaK subunit</td><td>161.34</td><td>0.77</td></tr><tr><td>A1S_1341</td><td>enoyl-CoA hydratase/carnithine racemase</td><td>28.43</td><td>0.56</td></tr><tr><td>A1S_1344</td><td>thiolase</td><td>14.31</td><td>0.40</td></tr><tr><td>A1S_1357</td><td>alanine racemase</td><td>4.59</td><td>1.23</td></tr><tr><td>A1S_1370</td><td>oxidoreductase</td><td>2.67</td><td>0.80</td></tr><tr><td>A1S_1376</td><td>acyl-CoA dehydrogenase</td><td>11.34</td><td>3.14</td></tr><tr><td>A1S_1377</td><td>acrR family transcriptional regulator</td><td>4.28</td><td>0.56</td></tr><tr><td>A1S_1385</td><td>hypothetical protein</td><td>9.16</td><td>8.68</td></tr><tr><td>A1S_1454</td><td>transmembrane arsenate pump protein</td><td>54.33</td><td>from zero to 27.07**</td></tr><tr><td>A1S_1507</td><td>fimbrial protein</td><td>17.73</td><td>18.49</td></tr><tr><td>A1S_1530</td><td>SSS family major sodium/proline symporter</td><td>0.29</td><td>1.03</td></tr><tr><td>A1S_1541</td><td>hypothetical protein</td><td>8.27</td><td>9.19</td></tr><tr><td>A1S_1572</td><td>30S ribosomal protein S1</td><td>1.74</td><td>0.88</td></tr><tr><td>A1S_1617</td><td>30S ribosomal protein S20</td><td>4.08</td><td>1.15</td></tr><tr><td>A1S_1631</td><td>iron-binding protein</td><td>0.71</td><td>0.13</td></tr><tr><td>A1S_1637</td><td>DNA-binding protein HU-beta</td><td>1.13</td><td>5.12</td></tr><tr><td>A1S_1657</td><td>siderophore biosynthesis protein</td><td>13.61</td><td>2.89</td></tr><tr><td>A1S_1687</td><td>transcriptional regulator</td><td>1.84</td><td>0.00</td></tr><tr><td>A1S_1726</td><td>aspartate ammonia-lyase</td><td>0.33</td><td>0.22</td></tr></tbody></table>
<table><caption>Table 3 (continued).</caption><thead><tr><th>Gene Id*</th><th>Gene description</th><th>Fold change Biofilm vs. exponential phase cells</th><th>Fold change Biofilm vs. stationary phase cells</th></tr></thead><tbody><tr><td>A1S_1731</td><td>acetoacetyl-CoA transferase subunit beta</td><td>29.63</td><td>3.12</td></tr><tr><td>A1S_1732</td><td>acetoacetyl-CoA transferase subunit alpha</td><td>78.74</td><td>5.03</td></tr><tr><td>A1S_1736</td><td>hypothetical protein</td><td>6.21</td><td>1.47</td></tr><tr><td>A1S_1755</td><td>RND efflux pump AdeT</td><td>from zero to 17.27**</td><td>from zero to 17.27**</td></tr><tr><td>A1S_1924</td><td>cytochrome d terminal oxidase polypeptide subunit I</td><td>0.22</td><td>0.21</td></tr><tr><td>A1S_1925</td><td>cytochrome d terminal oxidase polypeptide subunit II</td><td>0.27</td><td>0.22</td></tr><tr><td>A1S_1926</td><td>hypothetical protein</td><td>0.17</td><td>0.03</td></tr><tr><td>A1S_1932</td><td>hypothetical protein</td><td>1.88</td><td>0.11</td></tr><tr><td>A1S_1965</td><td>UDP-N-acetylglucosamine acyltransferase</td><td>3.57</td><td>1.48</td></tr><tr><td>A1S_2072</td><td>universal stress family protein</td><td>0.52</td><td>0.21</td></tr><tr><td>A1S_2091</td><td>hypothetical protein</td><td>24.78</td><td>3.28</td></tr><tr><td>A1S_2093</td><td>hypothetical protein</td><td>1.13</td><td>0.01</td></tr><tr><td>A1S_2098</td><td>alcohol dehydrogenase</td><td>13.14</td><td>130.77</td></tr><tr><td>A1S_2102</td><td>aldehyde dehydrogenase 1</td><td>2.59</td><td>8.55</td></tr><tr><td>A1S_2148</td><td>acetyl-CoA synthetase/AMP-(fatty) acid ligase</td><td>12.90</td><td>0.51</td></tr><tr><td>A1S_2149</td><td>acyl CoA dehydrogenase oxidoreductase protein</td><td>8.68</td><td>3.11</td></tr><tr><td>A1S_2150</td><td>oxidoreductase</td><td>5.52</td><td>1.13</td></tr><tr><td>A1S_2164</td><td>phosphoenolpyruvate synthase</td><td>1.04</td><td>0.28</td></tr><tr><td>A1S_2183</td><td>signal peptide</td><td>0.59</td><td>0.03</td></tr><tr><td>A1S_2214</td><td>protein CsuD</td><td>180.04</td><td>from zero to 89.72**</td></tr><tr><td>A1S_2215</td><td>protein CsuC</td><td>201.23</td><td>from zero to 33.43**</td></tr><tr><td>A1S_2218</td><td>protein CsuA/B</td><td>164.40</td><td>from zero to 1122.03**</td></tr><tr><td>A1S_2261</td><td>cold shock protein</td><td>5.09</td><td>1.12</td></tr><tr><td>A1S_2289</td><td>signal peptide</td><td>20.61</td><td>2.63</td></tr><tr><td>A1S_2322</td><td>elongation factor Ts</td><td>1.50</td><td>4.89</td></tr><tr><td>A1S_2382</td><td>BasD (iron acquisition systems)</td><td>72.89</td><td>from zero to 24.22**</td></tr><tr><td>A1S_2385</td><td>ferric acinetobactin receptor</td><td>6.48</td><td>6.48</td></tr><tr><td>A1S_2386</td><td>ferric acinetobactin binding protein</td><td>9.10</td><td>from zero to 48.57**</td></tr><tr><td>A1S_2390</td><td>acinetobactin biosynthesis protein</td><td>34.11</td><td>18.97</td></tr><tr><td>A1S_2447</td><td>EsvD ABC transporter</td><td>7.56</td><td>14.72</td></tr><tr><td>A1S_2449</td><td>aromatic amino acid APC transporter</td><td>16.58</td><td>1.15</td></tr><tr><td>A1S_2450</td><td>pyruvate decarboxylase</td><td>8.17</td><td>0.22</td></tr><tr><td>A1S_2452</td><td>NAD-dependent aldehyde dehydrogenases</td><td>1.71</td><td>0.15</td></tr><tr><td>A1S_2458</td><td>fatty acid desaturase</td><td>0.24</td><td>0.15</td></tr><tr><td>A1S_2496</td><td>phosphoserine phosphatase</td><td>0.3</td><td>0.01</td></tr><tr><td>A1S_2534</td><td>sulfate transport protein</td><td>21.12</td><td>from zero to 24.66**</td></tr><tr><td>A1S_2696</td><td>hypothetical protein</td><td>1.35</td><td>0.20</td></tr><tr><td>A1S_2705</td><td>hypothetical protein</td><td>0.21</td><td>0.09</td></tr><tr><td>A1S_2718</td><td>succinyl-CoA synthetase subunit beta</td><td>1.14</td><td>7.65</td></tr><tr><td>PLD-ONE</td><td>superoxide dismutase subunit alpha</td><td>1.07</td><td>5.76</td></tr></tbody></table>

<table><caption>Table 3 (continued).</caption><thead><tr><th>Gene Id*</th><th>Gene description</th><th>Fold change Biofilm vs. exponential phase cells</th><th>Fold change Biofilm vs. stationary phase cells</th></tr></thead><tbody><tr><td>A1S_2753</td><td>hypothetical protein</td><td>1.66</td><td>3.36</td></tr><tr><td>A1S_2840</td><td>outer membrane protein A</td><td>0.60</td><td>0.74</td></tr><tr><td>A1S_2889</td><td>signal peptide</td><td>46.50</td><td>25.85</td></tr><tr><td>A1S_2893</td><td>hypothetical protein</td><td>64.83</td><td>from zero to 32.31**</td></tr><tr><td>A1S_3043</td><td>hypothetical protein</td><td>3.91</td><td>1.5</td></tr><tr><td>A1S_3055</td><td>50S ribosomal protein L17</td><td>2.54</td><td>3.88</td></tr><tr><td>A1S_3056</td><td>DNA-directed RNA polymerase subunit alpha</td><td>1.89</td><td>3.04</td></tr><tr><td>A1S_3057</td><td>30S ribosomal protein S4</td><td>1.89</td><td>2.99</td></tr><tr><td>A1S_3058</td><td>30S ribosomal protein S11</td><td>2.00</td><td>2.58</td></tr><tr><td>A1S_3061</td><td>preprotein translocase subunit SecY</td><td>2.64</td><td>1.71</td></tr><tr><td>A1S_3062</td><td>50S ribosomal protein L15</td><td>2.79</td><td>2.07</td></tr><tr><td>A1S_3063</td><td>50S ribosomal protein L30</td><td>2.54</td><td>3.31</td></tr><tr><td>A1S_3064</td><td>30S ribosomal protein S5</td><td>2.97</td><td>3.54</td></tr><tr><td>A1S_3065</td><td>50S ribosomal protein L18</td><td>3.37</td><td>3.34</td></tr><tr><td>A1S_3066</td><td>50S ribosomal protein L6</td><td>2.55</td><td>2.29</td></tr><tr><td>A1S_3068</td><td>30S ribosomal protein S14</td><td>2.90</td><td>2.27</td></tr><tr><td>A1S_3069</td><td>50S ribosomal protein L5</td><td>2.13</td><td>1.66</td></tr><tr><td>A1S_3070</td><td>50S ribosomal protein L24</td><td>2.26</td><td>2.05</td></tr><tr><td>A1S_3073</td><td>50S ribosomal protein L29</td><td>2.01</td><td>3.00</td></tr><tr><td>A1S_3074</td><td>50S ribosomal protein L16</td><td>2.05</td><td>3.56</td></tr><tr><td>A1S_3075</td><td>30S ribosomal protein S3</td><td>1.73</td><td>3.06</td></tr><tr><td>A1S_3077</td><td>50S ribosomal protein L2</td><td>1.75</td><td>1.77</td></tr><tr><td>A1S_3079</td><td>50S ribosomal protein L4</td><td>1.87</td><td>1.59</td></tr><tr><td>A1S_3080</td><td>50S ribosomal protein L3</td><td>2.13</td><td>1.29</td></tr><tr><td>A1S_3104</td><td>ATP-dependent RNA helicase</td><td>1.64</td><td>0.23</td></tr><tr><td>A1S_3108</td><td>coproporphyrinogen III oxidase</td><td>0.28</td><td>0.33</td></tr><tr><td>A1S_3113</td><td>hypothetical protein</td><td>0.90</td><td>0.04</td></tr><tr><td>A1S_3134</td><td>glutamate dehydrogenase</td><td>1.26</td><td>3.21</td></tr><tr><td>A1S_3161</td><td>50S ribosomal protein L19</td><td>2.57</td><td>2.14</td></tr><tr><td>A1S_3185</td><td>glutamate synthase subunit alpha</td><td>0.40</td><td>0.48</td></tr><tr><td>A1S_3231</td><td>acetyl-CoA hydrolase/transferase</td><td>3.42</td><td>0.96</td></tr><tr><td>A1S_3297</td><td>outer membrane protein</td><td>1.17</td><td>3.80</td></tr><tr><td>A1S_3300</td><td>acetate permease</td><td>17.44</td><td>1.23</td></tr><tr><td>A1S_3301</td><td>hypothetical protein</td><td>5.77</td><td>0.52</td></tr><tr><td>A1S_3303</td><td>hypothetical protein</td><td>5.78</td><td>0.47</td></tr><tr><td>A1S_3309</td><td>acetyl-CoA synthetase</td><td>4.17</td><td>1.96</td></tr><tr><td>A1S_3328</td><td>pyruvate dehydrogenase subunit E1</td><td>0.50</td><td>0.92</td></tr><tr><td>A1S_3350</td><td>hypothetical protein</td><td>0.38</td><td>0.89</td></tr><tr><td>A1S_3402</td><td>arginase/agmatinase/ formimionoglutamate hydrolase</td><td>3.72</td><td>7.33</td></tr><tr><td>A1S_3404</td><td>amino acid APC transporter</td><td>3.88</td><td>4.18</td></tr><tr><td>A1S_3405</td><td>histidine ammonia-lyase</td><td>3.13</td><td>3.94</td></tr><tr><td>A1S_3406</td><td>urocanate hydratase</td><td>3.92</td><td>0.85</td></tr><tr><td>A1S_3407</td><td>urocanase</td><td>4.52</td><td>2.63</td></tr><tr><td>A1S_3413</td><td>APC family aromatic amino acid transporter</td><td>66.30</td><td>18.38</td></tr><tr><td>A1S_3414</td><td>fumarylacetoacetase</td><td>60.10</td><td>20.87</td></tr><tr><td>A1S_3415</td><td>maleylacetoacetate isomerase</td><td>24.49</td><td>from zero to 77.30**</td></tr></tbody></table>
<table><caption>Table 3 (continued).</caption><thead><tr><th>Gene Id*</th><th>Gene description</th><th>Fold change Biofilm vs. exponential phase cells</th><th>Fold change Biofilm vs. stationary phase cells</th></tr></thead><tbody><tr><td>A1S_3416</td><td>glyoxalase/bleomycin resistance protein/dioxygenas</td><td>24.26</td><td>1.02</td></tr><tr><td>A1S_3418</td><td>4-hydroxyphenylpyruvate dioxygenase</td><td>78.62</td><td>12.69</td></tr><tr><td>A1S_3463</td><td>diaminopimelate decarboxylase</td><td>0.41</td><td>0.08</td></tr><tr><td>A1S_3473</td><td>hypothetical protein</td><td>0.67</td><td>0.20</td></tr><tr><td>A1S_3475</td><td>hypothetical protein</td><td>1.15</td><td>0.19</td></tr></tbody></table>

The data were filtered based on a *p* value < 0.001.

\* Genes that significantly differed in their expression values with a *p* value below 0.001 in at least one of the two profile comparisons are listed in this table.

\*\* In these cases the expression values are absolute and no expression was detected under planktonic conditions.
doi: 10.1371/journal.pone.0072968.t003

(A1S_1146 and A1S_1147) were likewise only expressed in biofilms. Other groups of genes comprised those involved in efflux systems (A1S_1117, A1S_1751, and A1S_1755), or amino acid metabolism and transport (A1S_0956 and A1S_2302), or encoded an acyl carrier protein (A1S_0114). The complete list is shown in Table 4.

Further analysis of these 55 genes using Blast2GO revealed that most were involved in regulation of transcription and fewer to processes such as electron transport, acyl carrier biosynthesis, transmembrane transport, DNA replication, and siderophore biosynthesis (Figure 2a). The main molecular functions ascribed to the 55 genes are shown in Figure 2b, with oxidoreductase, transporters, DNA binding, and transcription factors activities predominating. The cellular location of the proteins encoded by the 55 genes is illustrated in Figure 2c, which shows that most of the proteins were located in a transcription factor complex or in the cell membrane.

Decrease in biofilm formation ability by gene disruption and knockout mutants

Five genes over-expressed in the biofilm vs. planktonic cells, as previously confirmed by qRT-PCR (Table 5), were selected for gene disruption by insertion of the plasmid pCR-Blunt-II-TOPO via single crossover recombination, as described in Materials and methods. These genes were A1S_0114 (encoding an acyl carrier protein expressed only in biofilms and inhibited in planktonic cells), A1S_0302 (encoding a hypothetical protein whose expression was ca. 27-fold higher in biofilms than in stationary-phase cells), A1S_1507 (encoding a fimbrial protein with ca. 18-fold higher expression in biofilms than in planktonic cells), A1S_3168 (encoding a pilus assembly protein PilW expressed in biofilms and repressed in stationary-phase cells, see Table S3), and A1S_2042 (a transcriptional regulator of the TetR family expressed in biofilms but inhibited in planktonic cells). The resulting mutant strains were used to evaluate their ability to form biofilms compared to the wild-type strain. As shown in Figure 3, biofilm formation ability was

<table><caption>Table 4. List of genes expressed only in biofilm cells and inhibited in planktonic cells.</caption><thead><tr><th>Gene Id</th><th>Expression value in biofilm cells</th><th>Gene description</th></tr></thead><tbody><tr><td>A1S_0079</td><td>0.47</td><td>N-acetyltransferase GNAT family (98% Ab SDF)*</td></tr><tr><td>A1S_0114</td><td>127.96</td><td>acyl carrier protein</td></tr><tr><td>A1S_0547</td><td>1.22</td><td>transcriptional regulator</td></tr><tr><td>A1S_0595</td><td>0.56</td><td>membrane protein (100% Ab MDR-TJ)*</td></tr><tr><td>A1S_0648</td><td>1.41</td><td>hypothetical protein</td></tr><tr><td>A1S_0741</td><td>0.19</td><td>hypothetical protein</td></tr><tr><td>A1S_0945</td><td>0.75</td><td>ferredoxin</td></tr><tr><td>A1S_0946</td><td>1.13</td><td>hypothetical protein</td></tr><tr><td>A1S_0956</td><td>1.13</td><td>L-aspartate dehydrogenase</td></tr><tr><td>A1S_0969</td><td>0.19</td><td>transketolase</td></tr><tr><td>A1S_1116</td><td>1.03</td><td>vanillate O-demethylase oxygenase subunit</td></tr><tr><td>A1S_1117</td><td>2.16</td><td>MFS superfamily vanillate transporter</td></tr><tr><td>A1S_1121</td><td>0.19</td><td>lipase/esterase</td></tr><tr><td>A1S_1125</td><td>0.38</td><td>transferase</td></tr><tr><td>A1S_1133</td><td>0.75</td><td>flavin-binding monooxygenase</td></tr><tr><td>A1S_1146</td><td>1.41</td><td>site-specific DNA-methyltransferase</td></tr><tr><td>A1S_1147</td><td>1.88</td><td>DNA methylase-like protein</td></tr><tr><td>A1S_1198</td><td>0.19</td><td>extracellular nuclease</td></tr><tr><td>A1S_1256</td><td>0.38</td><td>transcriptional regulator</td></tr><tr><td>A1S_1276</td><td>0.28</td><td>hypothetical protein</td></tr><tr><td>A1S_1278</td><td>20.35</td><td>allophanate hydrolase subunit 2</td></tr><tr><td>A1S_1349</td><td>0.47</td><td>thioesterase</td></tr><tr><td>A1S_1366</td><td>1.03</td><td>transporter LysE family</td></tr><tr><td>A1S_1430</td><td>0.28</td><td>LysR family malonate utilization transcriptional regulator</td></tr><tr><td>A1S_1452</td><td>0.94</td><td>arsenate reductase</td></tr><tr><td>A1S_1583</td><td>0.84</td><td>hypothetical protein</td></tr><tr><td>A1S_1585</td><td>0.56</td><td>replicative DNA helicase</td></tr><tr><td>A1S_1590</td><td>0.94</td><td>peptidase U35 phage prohead HK97</td></tr><tr><td>A1S_1622</td><td>1.13</td><td>hypothetical protein</td></tr><tr><td>A1S_1699</td><td>3.28</td><td>pyruvate/2-oxoglutarate dehydrogenase complex</td></tr><tr><td>A1S_1719</td><td>0.38</td><td>4Fe-4S ferredoxin iron-sulfur binding</td></tr><tr><td>A1S_1751</td><td>2.34</td><td>AdeA membrane fusion protein</td></tr><tr><td>A1S_1755</td><td>17.27</td><td>AdeT</td></tr><tr><td>A1S_1763</td><td>1.50</td><td>transcriptional regulator</td></tr><tr><td>A1S_1853</td><td>0.38</td><td>hypothetical protein</td></tr><tr><td>A1S_1887</td><td>0.28</td><td>major facilitator superfamily permease</td></tr><tr><td>A1S_1958</td><td>0.38</td><td>transcriptional regulator</td></tr><tr><td>A1S_2015</td><td>0.38</td><td>DNA-directed DNA polymerase</td></tr><tr><td>A1S_2028</td><td>0.19</td><td>phage putative head morphogenesis protein</td></tr><tr><td>A1S_2029</td><td>11.30</td><td>hypothetical protein</td></tr><tr><td>A1S_2033</td><td>0.47</td><td>hypothetical protein</td></tr><tr><td>A1S_2035</td><td>0.28</td><td>hypothetical protein</td></tr><tr><td>A1S_2042</td><td>2.72</td><td>transcriptional regulator (TetR family)</td></tr><tr><td>A1S_2151</td><td>0.75</td><td>transcriptional regulator (AraC family)</td></tr><tr><td>A1S_2208</td><td>0.38</td><td>transcriptional regulator</td></tr><tr><td>A1S_2216</td><td>11.96</td><td>CsuB</td></tr><tr><td>A1S_2217</td><td>3.84</td><td>CsuA</td></tr><tr><td>A1S_2302</td><td>0.75</td><td>ABC lysine-arginine-ornithine transporter</td></tr><tr><td>A1S_2380</td><td>8.81</td><td>acinetobactin biosynthesis protein</td></tr></tbody></table>
<table><caption>Table 4 (continued).</caption><thead><tr><th rowspan="2">Gene Id</th><th colspan="2">Expression value in</th></tr><tr><th>biofilm cells</th><th>Gene description</th></tr></thead><tbody><tr><td>A1S_2388</td><td>1.69</td><td>putative ferric acinetobactin transport system</td></tr><tr><td>A1S_2408</td><td>0.09</td><td>HNH endonuclease (93% Ab MDR-TJ)*</td></tr><tr><td>A1S_2580</td><td>1.40</td><td>23-dihydro-2,3-dihydroxybenzoate synthetase, isochorismatase</td></tr><tr><td>A1S_3120</td><td>0.09</td><td>hypothetical protein</td></tr><tr><td>A1S_3255</td><td>0.09</td><td>transcriptional regulator AraC/XylS family protein</td></tr><tr><td>A1S_3260</td><td>1.59</td><td>hypothetical protein</td></tr></tbody></table>

severely hindered (~8-fold reduction) in all of the mutant strains.

To obtain a stable mutant free of antibiotic resistance markers or potential polar effects, the A1S_0114 gene was deleted from the genome of A. baumannii ATCC 17978 using the pMo130 vector, as described in Material and methods. As shown in Figure 4, the biofilm formation ability of the stable A1S_0114 knock-out (KO) mutant was significantly reduced (< 3-fold) compared to the wild-type strain.

The relationship of the gene A1S_0114 to genes related to homoserine lactone synthesis (A1S_0109, A1S_0112 and A1S_0113) was examined in qRT-PCR assays. As shown in Figure 5, genes A1S_0109, A1S_0112, A1S_0113 and A1S_0114 were over-expressed in the late stationary phase of growth compared to the exponential phase in the wild-type strain. When gene A1S_0114 was deleted from the chromosome (yielding the stable A1S_0114 KO mutant strain), the expression levels of genes related to homoserine lactone synthesis (A1S_0109, A1S_0112 and A1S_0113) were considerably reduced in the late stationary phase of growth (83, 68 and 73%, respectively) of the resulting mutant compared with the wild-type strain.

Discussion

In the present work, we successfully used Illumina RNA-sequencing to establish the complete transcriptional profile of A. baumannii strain ATCC 17978 grown in planktonic and sessile (biofilm) modes. To obtain an overview of the temporal regulation of gene expression, planktonic cells were harvested during the exponential and late stationary phases of growth. A similar strategy was previously used in a proteomic study demonstrating the growth-dependent regulation of many proteins [58]. In another proteomic study of A. baumannii ATCC 17978 [34], planktonic and sessile cells were shown to exhibit distinct proteomic profiles, indicating that biofilms are not simply surface-attached stationary-phase cells.

Our data revealed that although many genes were constitutively expressed in both biofilm and planktonic cells, others differed in their growth-dependent expression, with clearly distinct and specific expression patterns between

sessile biofilm cells and cells in either phase of planktonic growth. Among the 1621 genes over-expressed in biofilms, 55 genes were only expressed in sessile cells and were totally inhibited in planktonic cells. The majority of the 55 genes encoded proteins involved in functions and mechanisms already known to be related with biofilm formation and maintenance whereas others were detected in this study for the first time. The presence of 12 genes encoding uncharacterized proteins further highlights the deficits in our knowledge of the specific genes associated with biofilm in A. baumannii. One of these genes (A1S_0302) was selected for gene disruption procedures because of its high level of expression in biofilm cells; indeed, the corresponding mutant strain was significantly deficient in biofilm formation. In addition, nine transcriptional regulators were found to be expressed only by biofilm cells, suggesting that biofilm formation and maintenance is controlled by specific molecules that are either not expressed, silenced, or not operative in planktonic cells. Gaddy and Actis [26] suggested that the regulatory process associated with biofilm formation includes the sensing of bacterial density, the presence of nutrients, and the concentration of free cations. Some of these extracellular signals are controlled by two-component regulatory systems such as BfmR/S. The *bfmS* gene encodes a sensor kinase that receives extracellular signals and phosphorylates the product of the *bfmR* gene, a response regulator. Tomaras et al. [27] studied the BfmR/S system in A. baumannii strain 19606, where this two-component regulator is required for the activation of the usher-chaperone assembly system involved in pili formation, a feature of biofilms. Based on their study of *P. aeruginosa*, Petrova et al. [59] proposed a role for BfmR in biofilm development, by limiting bacteriophage-mediated lysis and subsequent DNA release. According to our data, expression of the *bfmR* gene, identified as A1S_0748 by Liou et al. [60], was *ca.* five-fold higher in biofilm cells than in stationary cells. However, BfmR cannot be claimed as a biofilm-specific molecule of the strain 17978 since it was also expressed in the planktonic cells. No significant similarities were found in the databases for the nine transcriptional regulators described herein as biofilm specific. The mutant strain generated by the disruption of one of these genes (A1S_2042) showed an important decrease in biofilm formation ability relative to the wild-type strain. A1S_2042 appears to be a transcriptional regulator of the TetR family and could play an important role in biofilm regulation. Together with the other uncharacterized biofilm-specific transcriptional regulators, all of which were expressed at low but significant levels in biofilms, A1S_2042 merits further study to gain insight into the complex regulatory networks involved in biofilm formation and maintenance.

The CsuA/BABCDE chaperone-usher pili assembly system is involved in the adherence of A. baumannii strain 19606 biofilm to abiotic surfaces [26,27]. In the present work, Csu A/B, C, D and E were highly over-expressed in biofilms vs. planktonic cells. Of particular interest is CsuA/B, which is predicted to form part of the type I pili rod [40]. While the gene encoding CsuA/B was not expressed in stationary cells, its expression was greatly enhanced (value of 1122) in biofilms, although it was also expressed in exponentially growing cells

<footnote>* These genes are annotated in strain ATCC 17978. Similarities to sequences in the databases are indicated.</footnote>

<footnote>doi: 10.1371/journal.pone.0072968.t004</footnote>
<figure><img src="image_3.png" /></figure>

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

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Sequence distribution of genes expressed only in biofilm-associated cells and inhibited in planktonic cells. Genes involved in A) biological processes, B) molecular functions, and C) cellular components. The results were filtered by the number of sequences (cutoff = 1, 4, and 1, respectively).</figcaption></figure>

doi: 10.1371/journal.pone.0072968.g002
<figure><img src="image_6.png" /><figcaption>Figure 3. Quantification of biofilm formation by the wild-type strain (ATCC 17978) and strains with chromosomal disruptions in the genes A1S_0114, A1S_0302, A1S_1507, A1S_3168 and A1S_2042.</figcaption></figure>

doi: 10.1371/journal.pone.0072968.g003

<table><caption>**Table 5.** Expression levels of genes A1S_0114, A1S_0302, A1S_1507, A1S_2042, and A1S_3168 in biofilm and planktonic cells as measured by qRT-PCR.</caption><thead><tr><th>Gene Id</th><th>Expression level in exponential phase cells*</th><th>Expression level at the stationary phase cells</th><th>Expression level in biofilm cells</th></tr></thead><tbody><tr><td>A1S_0114</td><td>1 ± 0.319</td><td>0.219 ± 0.234</td><td>6.023 ± 1.493</td></tr><tr><td>A1S_0302</td><td>1 ± 0.339</td><td>3.099 ± 0.847</td><td>4.069 ± 0.599</td></tr><tr><td>A1S_1507</td><td>1 ± 0.073</td><td>1.535 ± 0.215</td><td>7.761 ± 0.719</td></tr><tr><td>A1S_2042</td><td>1 ± 0.488</td><td>4.589 ± 2.152</td><td>39.35 ± 12.670</td></tr><tr><td>A1S_3168</td><td>1 ± 0.158</td><td>1.237 ± 0.076</td><td>2.087 ± 0.522</td></tr></tbody></table>

doi: 10.1371/journal.pone.0072968.t005

(value of ca 164). Moreover, we detected CsuA (A1S_2216) and CsuB (A1S_2217) transcripts only in biofilms and not in planktonic cells. These results indicate that the complete Csu operon is highly active in the biofilms analyzed in this study. Nevertheless, it should be noted that MacQueary and Actis [61], found strong variations in the CsuA/BABCDE chaperone-usher pili assembly system and other motility factors among different strains of *A. baumannii* attached to abiotic surfaces.

This finding may pose a challenge in the treatment of the infections caused by this bacterium, if biofilm formation on abiotic surfaces is chosen as a target for the development of new antimicrobial agents.

Two genes coding for a fimbrial protein (A1S_1507) and a pilus assembly protein PilW (A1S_3168), different from the CsuA/BABCDE chaperone-usher pili assembly system, were over-expressed in biofilm vs. planktonic cells. The disruption of these two genes in the genome of *A. baumannii* revealed their involvement in biofilm formation and suggested that the biofilm analyzed here could require multiple pili systems to maintain its cohesive structure. Pilus and fimbriae are important for the initial step of bacterial adhesion, which is followed by the production of exopolysaccharides, an important constituent of mature biofilms that suppresses neutrophil activity and contributes to resistance. Variation in the expression of factors involved in these pathways may account for the different capacity of bacterial strains to form biofilms and therefore to colonize or infect the host environment [62].

*A. baumannii* secretes a variety of molecules involved in iron acquisition including siderophores such as acinetobactin. The iron concentration in the medium acts as an important environmental signal that induces the expression of adhesion factors, thus playing a critical role in biofilm formation [63]. However, there is wide variability in the expression of iron uptake molecules, even between strains isolated during the

<footnote>* The expression levels of each of the five genes were determined with respect to the exponential growth phase value, defined as 1.</footnote>
<figure><img src="image_7.png" /><figcaption>Figure 4. Quantification of biofilm formation by the wild-type strain (ATCC 17978), a stable knockout mutant strain lacking the gene A1S_0114 (ATCC Δ0114), the same mutant strain containing the pET-RA plasmid (ATCC Δ0114 + PETRA), and a mutant strain containing the pET-RA plasmid harboring the A1S_0114 gene (ATCC Δ0114 + PETRA + 0114).</figcaption></figure>

doi: 10.1371/journal.pone.0072968.g004

<figure><img src="image_8.png" /><figcaption>Figure 5. Comparison of the expression levels of genes related to homoserine lactone synthesis (A1S_0109, A1S_0112 and A1S_0113) in the wild-type strain A. baumannii 17978 and in the A1S_0114 knock-out (KO) strain as determined by real-time qRT-PCR assays.</figcaption></figure>

doi: 10.1371/journal.pone.0072968.g005
same outbreak [61]. In our experimental model, several genes involved in iron acquisition were over-expressed in biofilm vs. planktonic cells, while some genes related to acinetobactin (A1S_2380 and A1S_2388) and ferredoxin (A1S_0945 or A1S_1719) were expressed only in biofilms and totally inhibited in planktonic cells. The exclusively expression of acinetobactin genes in biofilm cells could be explained in terms of iron starved conditions in the sessile cells compared with an iron-rich medium used for growing planktonic cells. Eijkelkamp et al. [40] found transcriptional changes in genes involved in motility when *A. baumannii* was grown under iron-limiting conditions. As shown by our data, the biofilm is a resistance mode where cells clearly over-express many genes related to iron acquisition systems. It is known that the ability of *A. baumannii* to obtain and utilize resources such as iron is an important factor for bacterial survival but it seems to also be essential for biofilm formation and maintenance, given that bacteria able to form biofilms actively search for iron [63]. This scenario was reflected in our study by the over-expression of many genes involved in iron acquisition and transport.

The detection of two genes involved in DNA methylation and expressed exclusively in biofilms suggests a role for DNA methylation in the regulation of biofilm-associated processes. In addition, several genes encoding efflux system components were activated in the biofilm cells, including the gene encoding resistance-nodulation-cell division type efflux pump (RND pump), involved in bacterial resistance to a number of antibiotics. Our results indicate that the up-regulation of efflux pumps is a mechanism of antibiotic resistance that operates in the mature biofilm [34].

Another factor previously described as involved in biofilm formation is the homolog of the staphylococcal protein Bap, studied in *A. baumannii* 307-0294. The protein is a surface adhesin that mediates primary attachment to both biotic and abiotic surfaces and is involved in intercellular adhesion within the mature biofilm [29]. The *A. baumannii* 17978 genome (NC_009085.1) contains two loci homologous to the 5' and 3' ends of the *bap* locus defined in *A. baumannii* strain 307-0294 [29]. These two regions correspond to genes A1S_2724 and A1S_2696 (annotated as a hemaglutinin/hemolysin like protein and a hypothetical protein, respectively) [46] that were over-expressed in biofilms vs. exponential cells, suggesting that *A. baumannii* Bap-related proteins in the strain 17978 could also enhance the cell to cell interactions that support biofilm maturation.

Amino acid metabolism also clearly differed in our biofilm experimental design with respect to planktonic cells, as several genes involved in the metabolism and transport of amino acids were differentially expressed. Our results not only corroborate the hypothesis formulated by Cabral et al. [34] regarding the importance of histidine metabolism in biofilm formation but also extend it, based on our detection of genes involved in amino acid metabolism that were differentially expressed in biofilm cells and were not previously detected by proteomic analysis.

Cell surface membrane proteins may be essential to biofilm formation. Some of these proteins were differentially expressed in our biofilm cells, such as CarO (A1S_2538) and OprD-like (A1S_0201), which were up-regulated, while OmpA

(A1S_2840) was down-regulated. These results conflict somewhat with those of Cabral et al. [34], who found that OmpA was up-regulated in biofilm cells. Gaddy et al. [30] also described the importance of OmpA in biofilm formation in *A. baumannii* strain 19606. The discrepancy in the results can be explained by strain-dependent variations or different adhesion phenomena in response to diverse biotic or abiotic surface materials, as previously described [64]. However, Marti et al. [63], analyzed the proteome of *A. baumannii* strain 77 and found three mass isoforms identified as OmpA. In accordance with our results, OmpA was down-regulated in the biofilm, leading the authors to suggest that this porin participates in the initiation step of biofilm formation and that the subsequent iron starvation conditions encountered during biofilm maturation trigger a decrease in its expression. This may have been the case in our experimental model. The extracellular matrix that surrounds the biofilm protects the resident bacterial cells against a number of agents but it also limits bacterial access to fresh nutrients. Accordingly, an increase in the expression of transmembrane channels may be essential for the entrance of important nutrients. In the present work, the under-expression of OmpA was complemented by an over-expression of the porins OprD-like and CarO, which may have helped to maintain the permeability of the cells in the biofilm.

Although little is known about the factors involved in biofilm regulation, cell to cell signaling mediated by N-acyl-homoserine lactones has been implicated in gram-negative bacteria [65-67]. Indeed, we identified a group of genes (identifiers A1S_0112 to A1S_0118) over-expressed in biofilms vs. planktonic cells. This group of genes has been described as an operon related to quorum sensing and may be involved in the expression of the protein encoded by A1S_0109, the only homoserine lactone synthase described thus far in *A. baumannii* [46,68-70]. In our experimental model, this homoserine lactone synthase (A1S_0109) was over-expressed in biofilm vs. planktonic cells. Among the genes contained in the above-mentioned operon, A1S_0114 was exclusively expressed at high levels in biofilms but totally inhibited in planktonic cells. This gene encodes a small acidic acyl carrier protein (ACP) that is very abundant in bacteria, where it serves as an important acyl donor. ACP is first synthesized in its inactive form (apo-ACP) and then activated by an acyl carrier protein synthase [71]. In its activated form, ACP is essential for the synthesis of N-acyl-homoserine lactone, which is a substrate for the homoserine lactone synthase [72]. In this work, proteins encoded by the genes A1S_0112 and A1S_0113 were identified as an acyl-CoA synthetase and an acyl-CoA dehydrogenase, respectively, both of which are necessary for ACP activation. In the gene disruption and in the stable knock out A1S_0114 (ACP) mutants there was a notable decrease in biofilm formation ability compared to the wild-type strain, demonstrating the importance of this gene in biofilm formation. Our qRT-PCR results indicated reduced expression of the genes A1S_0112, A1S_0113 as well as the N-acyl-homoserine lactone synthase gene A1S_0109 in the stable A1S_0114 KO mutant, which presumably could affect quorum sensing and biofilm formation. Moreover, since our results were consistent with alterations in fatty acid metabolism in biofilms
vs. planktonic cells an alternative explanation for the decrease
in biofilm formation ability of the A1S_0114 mutant is that the
encoded ACP acts as an acyl donor associated with general
fatty acid metabolism.

Concluding remarks

The main goal of this study was to provide insight into the
molecular mechanisms underlying the ability of *A. baumannii* to
form biofilms. The expression profiles described herein allow
the definition of many genetic elements involved in the sessile
lifestyle of *A. baumannii*, including 55 genes exclusively
expressed in biofilm. Five genes were disrupted in the
chromosome and the corresponding mutant strains were
significantly hindered in their biofilm formation ability,
demonstrating their involvement in biofilm development. An
ACP-encoding gene that belongs to an operon involved in
quorum sensing mediated by a homoserine lactone was highly
over-expressed in our biofilm experimental model and its
inactivation significantly limited biofilm formation by cells of the
corresponding mutant strain.

The results described in this work constitute a basis for the
identification of new therapeutic targets and the design of new
drugs able to prevent infectious diseases related to biofilm
production by *A. baumannii*. It also serves as a starting point
for future studies of the complex network systems involved in
biofilm formation and maintenance, as well as the regulation of
these processes.

**Supporting Information**

**Figure S1. Gene level counts.** Left: boxplot (median, first and third quartiles and standard deviation) of the number of reads per gene. Right: density functions of the number of reads *per gene*.

(TIF)

**Figure S2. MD plots and correlation between samples.**
Upper right: MD plots showing (countsA+countsB)/2 against (countsA-countsB), with A and B being the samples shown on the diagonal.

(TIF)

**Figure S3. Sequence distribution of genes up-regulated in biofilm-associated cells.** The data were filtered based on *p* < 0.001 and with respect to biological processes. A) Exponentially growing cells, filtered by the number of sequences (cutoff 6). B) Stationary phase cells, filtered by the number of sequences (cutoff 1).

(TIF)

**Table S1.** Gene expression data from the complete transcriptome analysis of *Acinetobacter baumannii* ATCC 17978, showing gene expression levels in exponentially growing vs. stationary phase cells. Id: name or code of the region of interest; baseMean: mean of the two next columns; baseMeanA: normalized number of counts for sample A; baseMeanB: normalized number of counts for sample B; Fold-

change: baseMeanB/baseMeanA, log2Fold-change: log₂
baseMeanB/baseMeanA, pval: *p* value, padj: *p* value adjusted
for multiple testing, resVarA: variance of A, resVarB: variance
of B.A: stationary phase cells. B: exponential phase cells. NA,
non-applicable because of zero expression.

(XLSX)

**Table S2.** Gene expression data from the complete transcriptome analysis of *Acinetobacter baumannii* ATCC 17978, showing gene expression levels in biofilm-associated cells vs. exponentially growing cells. Id: name or code of the region of interest; baseMean: mean of the two next columns; baseMeanA: normalized number of counts for sample A; baseMeanB: normalized number of counts for sample B; Fold-change: baseMeanB/baseMeanA; log2Fold-change: log₂ baseMeanB/baseMeanA; pval: *p* value; padj: *p* value adjusted for multiple testing; resVarA: variance of A; resVarB: variance of B.A: exponential phase cells. B: biofilm-associated cells. NA, non-applicable because of zero expression.

(XLSX)

**Table S3.** Gene expression data from the complete transcriptome analysis of *Acinetobacter baumannii* ATCC 17978, showing gene expression levels of biofilm-associated vs. stationary phase cells. Id: name or code of the region of interest; baseMean: mean of the two next columns; baseMeanA: normalized number of counts for sample A; baseMeanB: normalized number of counts for sample B; Fold-change: baseMeanB/baseMeanA; log2Fold-change: log₂ baseMeanB/baseMeanA; pval: *p* value; padj: *p* value adjusted for multiple testing; resVarA: variance of A; resVarB: variance of B.A: stationary phase cells. B: biofilm-associated cells. NA: non-applicable because of zero expression.

(XLSX)

**Table S4.** The expression levels of genes down-regulated in biofilm-associated vs. stationary phase cells. The data were filtered based on *p* < 0.001. Id: name or code of the region of interest; baseMean: mean of the two next columns; baseMeanA: normalized number of counts for sample A; baseMeanB: normalized number of counts for sample B; Fold-change: baseMeanB/baseMeanA; log2Fold-change: log₂ baseMeanB/baseMeanA; pval: *p* value; padj: *p* value adjusted for multiple testing; resVarA: variance of A; resVarB: variance of B.A: stationary phase cells. B: biofilm-associated cells.

(XLSX)

**Table S5.** The expression levels of genes up-regulated in biofilm-associated vs. stationary phase cells. The data

were filtered based on *p* < 0.001. Id: name or code of the
region of interest; baseMean: mean of the two next columns;
baseMeanA: normalized number of counts for sample A;
baseMeanB: normalized number of counts for sample B; Fold-
change: baseMeanB/baseMeanA; log2Fold-change: log₂
baseMeanB/baseMeanA; pval: *p* value; padj: *p* value adjusted
for multiple testing; resVarA: variance of A; resVarB: variance
of B.A: stationary phase cells. B: biofilm-associated cells.
(XLSX)

**Table S6. The expression levels of genes down-regulated in exponentially growing vs. stationary phase cells.** The data were filtered based on *p* < 0.001. Id: name or code of the region of interest; baseMean: mean of the two next columns; baseMeanA: normalized number of counts for sample A; baseMeanB: normalized number of counts for sample B; Fold-change: baseMeanB/baseMeanA; log2Fold-change: log₂ baseMeanB/baseMeanA; pval: *p* value; padj: *p* value adjusted for multiple testing; resVarA: variance of A; resVarB: variance of B.A: stationary phase cells. B: exponential phase cells.
(XLSX)

**Table S7. The expression levels of genes up-regulated in exponentially growing vs. stationary phase cells, filtered based on *p* < 0.001.** Id: name or code of the region of interest; baseMean: mean of the two next columns; baseMeanA: normalized number of counts for sample A; baseMeanB: normalized number of counts for sample B; Fold-change: baseMeanB/baseMeanA; log2Fold-change: log₂ baseMeanB/baseMeanA; pval: *p* value; padj: *p* value adjusted for multiple testing; resVarA: variance of A; resVarB: variance of B.A: stationary phase cells. B: exponential phase cells.
(XLSX)

**Table S8. The expression levels of genes down-regulated in biofilm-associated vs. exponentially growing cells.** The data were filtered based on *p* < 0.001. Id: name or code of the region of interest; baseMean: mean of the two next columns; baseMeanA: normalized number of counts for sample A; baseMeanB: normalized number of counts for sample B; Fold-change: baseMeanB/baseMeanA; log2Fold-change: log₂ baseMeanB/baseMeanA; pval: *p* value; padj: *p* value adjusted for multiple testing; resVarA: variance of A; resVarB: variance of B.A: stationary phase cells. B: exponential phase cells.
(XLSX)

## References

1. Corbella X, Montero A, Pujol M, Domínguez 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. PubMed: 11060073.

2. 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 colonisation and infection. Clin Microbiol Infect 11: 540-546. doi: 10.1111/j.1469-0691.2005.01184.x. PubMed: 15966971.

3. Domenech de Cellès 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. doi:10.1371/journal.pone.0045758. PubMed: 23029226.

4. Montero A, Ariza J, Corbella X, Doménech A, Cabellos C et al. (2002) Efficacy of colistin versus beta-lactams, aminoglycosides, and rifampin as monotherapy in a mouse model of pneumonia caused by multiresistant *Acinetobacter baumannii*. Antimicrob Agents Chemother 46: 1946-1952. doi:10.1128/AAC.46.6.1946-1952.2002. PubMed: 12019113.

5. 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. doi:10.1016/j.biomaterials.2010.12.015. PubMed: 21295343.

6. 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. doi: 10.1186/1476-0711-3-14. PubMed: 15283864.

7. Ong CW, Lye DC, Khoo KL, Chua GS, Yeoh SF et al. (2009) Severe community-acquired *Acinetobacter baumannii* pneumonia: an emerging

change: baseMeanB/baseMeanA; log2Fold-change: log₂ baseMeanB/baseMeanA; pval: *p* value; padj: *p* value adjusted for multiple testing; resVarA: variance of A; resVarB: variance of B.A: exponential phase cells. B: biofilm-associated cells.
(XLSX)

**Table S9. The expression levels of genes up-regulated in biofilm-associated vs. exponentially growing cells.** The data were filtered based on *p* < 0.001. Id: name or code of the region of interest; baseMean: mean of the two next columns; baseMeanA: normalized number of counts for sample A; baseMeanB: normalized number of counts for sample B; Fold-change: baseMeanB/baseMeanA; log2Fold-change: log₂ baseMeanB/baseMeanA; pval: *p* value, padj: *p* value adjusted for multiple testing; resVarA: variance of A; resVarB: variance of B.A: exponential phase cells. B: biofilm-associated cells.
(XLSX)

## Acknowledgements

The authors thank Fiona Veira McTiernan for kindly revising the manuscript, Juan Andrés Vallejo for his technical support, and M.I. Voskuil for pMo130.

## Author Contributions

Conceived and designed the experiments: MP CG GB. Performed the experiments: SRF CG MP LAF MPC JV. Analyzed the data: AMA NRE AF MJG. Contributed reagents/ materials/analysis tools: AMA NRE AF MJG MT. Wrote the manuscript: MP MJG.

highly lethal infectious disease in the Asia-Pacific. Respirology 14: 1200-1205. doi:10.1111/j.1440-1843.2009.01630.x. PubMed: 19909464.

8. Bou G, Cerveró G, Domínguez MA, Quereda C, Martínez-Beltrán 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. doi: 10.1046/j.1469-0691.2000.00181.x. PubMed: 11284921.

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

10. Santillana E, Beceiro A, Bou G, Romero A (2007) Crystal structure of the carbapenemase OXA-24 reveals insights into the mechanism of carbapenem hydrolysis. Proc Natl Acad Sci U S A 104: 5354-5359. doi: 10.1073/pnas.0607557104. PubMed: 17374723.

11. Higgins PG, Dammhayn C, Hackel M, Seifert H (2010) Global spread of carbapenem-resistant *Acinetobacter baumannii*. J Antimicrob Chemother 65: 233-238. doi:10.1093/jac/dkp428. PubMed: 19996144.

12. De la Cruz MA, Calva E (2010) The complexities of porin genetic regulation. J Mol Microbiol Biotechnol 18: 24-36. doi: 10.1159/000274309. PubMed: 20068355.

13. Fernández-Cuenca F, Martínez-Martínez L, Conejo MC, Ayala JA, Perea EJ et al. (2003) Relationship between β-lactamase production, outer membrane protein and penicillin-binding protein profiles on the activity of carbapenems against clinical isolates of *Acinetobacter baumannii*. J Antimicrob Chemother 51: 565-574. doi:10.1093/jac/dkg097. PubMed: 12615856.

14. Srinivasan VB, Rajamohan G, Gebreyes WA (2009) Role of AbeS, a novel efflux pump of the SMR family of transporters, in resistance to antimicrobial agents in *Acinetobacter baumannii*. Antimicrob Agents
Chemother 53: 5312-5316. doi:10.1128/AAC.00748-09. PubMed: 19770280.

15. Roca I, Marti S, Espinal P, Martínez P, Gibert I et al. (2009) CraA, a major facilitator superfamily efflux pump associated with chloramphenicol resistance in *Acinetobacter baumannii*. Antimicrob Agents Chemother 53: 4013-4014. doi:10.1128/AAC.00584-09. PubMed: 19581458.

16. Coyne S, Guigon G, Courvalin P, Périchon B (2010) Screening and quantification of the expression of antibiotic resistance genes in *Acinetobacter baumannii* with a microarray. Antimicrob Agents Chemother 54: 333-340. doi:10.1128/AAC.01037-09. PubMed: 19884373.

17. Rajamohan G, Srinivasan VB, Gebreyes WA (2010) Novel role of *Acinetobacter baumannii* RND efflux transporters in mediating decreased susceptibility to biocides. J Antimicrob Chemother 65: 228-232. doi:10.1093/jac/dkp427. PubMed: 20008046.

18. Costerton JW, Stewart PS, Greenberg EP (1999) Bacterial biofilms: a common cause of persistent infections. Science 284: 1318-1322. doi: 10.1126/science.284.5418.1318. PubMed: 10334980.

19. Moreno-Paz M, Gómez MJ, Arcas A, Parro V (2010) Environmental transcriptome analysis reveals physiological differences between biofilm and planktonic modes of life of the iron oxidizing bacteria *Leptospirillum* spp. in their natural microbial community. BMC Genomics 11: 404. doi:10.1186/1471-2164-11-404. PubMed: 20576116.

20. Vidal R, Dominguez M, Urrutia H, Bello H, Gonzalez G et al. (1996) Biofilm formation by *Acinetobacter baumannii*. Microbios 86: 49-58. PubMed: 8771775.

21. 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. doi:10.1111/j.1469-0691.2007.01842.x. PubMed: 18005176.

22. Rodríguez-Baño J, Martí S, Soto S, Fernández-Cuenca F, Cisneros JM et al. (2008) Biofilm formation in *Acinetobacter baumannii*: associated features and clinical implications. Clin Microbiol Infect 14: 276-278. doi: 10.1111/j.1469-0691.2007.01916.x. PubMed: 18190568.

23. Martí S, Rodríguez-Baño J, Catel-Ferreira M, Jouenne T, Vila J et al. (2011) Biofilm formation at the solid-liquid and air-liquid interfaces by *Acinetobacter* species. BMC Res Notes 4: 5. doi: 10.1186/1756-0500-4-5. PubMed: 21223561.

24. Espinal P, Martí S, Vila J (2012) Effect of biofilm formation on the survival of *Acinetobacter baumannii* on dry surfaces. J Hosp Infect 80: 56-60. doi:10.1016/j.jhin.2011.08.013. PubMed: 21975219.

25. 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. doi:10.1099/mic.0.26541-0. PubMed: 14663080.

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

27. 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. doi:10.1099/mic. 0.2008/019471-0. PubMed: 18957593.

28. Petrova OE, Sauer K (2009) A novel signaling network essential for regulating *Pseudomonas aeruginosa* biofilm development. PLOS Pathog 5: e1000668.

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

30. 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. doi:10.1128/IAI.00096-09. PubMed: 19470746.

31. McConnell MJ, Actis L, Pachón J (2013) *Acinetobacter baumannii*: human infections, factors contributing to pathogenesis and animal models. FEMS Microbiol Rev 37: 130-155. PubMed: 22568581.

32. Hamilton S, Bongaerts RJ, Mulholland F, Cochrane B, Porter J et al. (2009) The transcriptional programme of *Salmonella enterica* serovar Typhimurium reveals a key role for tryptophan metabolism in biofilms. BMC Genomics 10: 599. doi:10.1186/1471-2164-10-599. PubMed: 20003355.

33. Valle J, Da Re S, Schmid S, Skurnik D, D'Ari R et al. (2008) The amino acid valine is secreted in continuous-flow bacterial biofilms. J Bacteriol 190: 264-274. doi:10.1128/JB.01405-07. PubMed: 17981982.

34. 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. doi:10.1021/pr101299j. PubMed: 21612302.

35. Kolodkin-Gal I, Romero D, Cao S, Clardy J, Kolter R et al. (2010) D-amino acids trigger biofilm disassembly. Science 328: 627-629. doi: 10.1126/science.1188628. PubMed: 20431016.

36. Vebø HC, Snipen L, Nes IF, Brede DA (2009) The transcriptome of the nosocomial pathogen *Enterococcus faecalis* V583 reveals adaptive responses to growth in blood. PLOS ONE 4: e7660. doi:10.1371/journal.pone.0007660. PubMed: 19888459.

37. Toledo-Arana A, Dussurget O, Nikitas G, Sesto N, Guet-Revillet H et al. (2009) The *Listeria* transcriptional landscape from saprophytism to virulence. Nature 459: 950-956. doi:10.1038/nature08080. PubMed: 19448609.

38. Whiteley M, Bangera MG, Bumgarner RE, Parsek MR, Teitzel GM et al. (2001) Gene expression in *Pseudomonas aeruginosa* biofilms. Nature 413: 860-864. doi:10.1038/35101627. PubMed: 11677611.

39. Hood MI, Jacobs AC, Sayood K, Dunman PM, Skaar EP (2010) *Acinetobacter baumannii* increases tolerance to antibiotics in response to monovalent cations. Antimicrob Agents Chemother 54: 1029-1041. doi:10.1128/AAC.00963-09. PubMed: 20028819.

40. Eijkelkamp BA, Hassan KA, Paulsen IT, Brown MH (2011) Investigation of the human pathogen *Acinetobacter baumannii* under iron limiting conditions. BMC Genomics 12: 126. doi:10.1186/1471-2164-12-126. PubMed: 21342532.

41. Mao C, Evans C, Jensen RV, Sobral BW (2008) Identification of new genes in *Sinorhizobium meliloti* using the Genome Sequencer FLX system. BMC Microbiol 8: 72. doi:10.1186/1471-2180-8-72. PubMed: 18454850.

42. Droege M, Hill B (2008) The Genome Sequencer FLX System--longer reads, more applications, straight forward bioinformatics and more complete data sets. J Biotechnol 136: 3-10. doi:10.1016/j.jbiotec. 2008.07.1845. PubMed: 18616967.

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

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

45. 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. PubMed: 20368969.

46. Smith MG, Gianoulis TA, Pukatzki S, Mekalanos JJ, Ornston LN et al. (2007) New insights into *Acinetobacter baumannii* pathogenesis revealed by high-density pyrosequencing and transposon mutagenesis. Genes Dev 21: 601-614. doi:10.1101/gad.1510307. PubMed: 17344419.

47. Jäger D, Sharma CM, Thomsen J, Ehlers C, Vogel J et al. (2009) Deep sequencing analysis of the *Methanosarcina mazei* Gö1 transcriptome in response to nitrogen availability. Proc Natl Acad Sci U S A 106: 21878-21882. doi:10.1073/pnas.0909051106. PubMed: 19996181.

48. Sharma CM, Hoffmann S, Darfeuille F, Reignier J, Findeiss S et al. (2010) The primary transcriptome of the major human pathogen *Helicobacter pylori*. Nature 464: 250-255. doi:10.1038/nature08756. PubMed: 20164839.

49. 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. doi:10.1128/JB.01445-09. PubMed: 20190049.

50. Dötsch 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. doi:10.1371/journal.pone.0031092. PubMed: 22319605.

51. Langmead B, Trapnell C, Pop M, Salzberg SL (2009) Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 10: R25. doi:10.1186/gb-2009-10-3-r25. PubMed: 19261174.

52. Anders S, Huber W (2010) Differential expression analysis for sequence count data. Genome Biol 11: R106. doi:10.1186/gb-2010-11-10-r106. PubMed: 20979621.

53. Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH et al. (2008) High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res 36: 3420-3435. doi:10.1093/nar/gkn176. PubMed: 18445632.
54. Tendolkar PM, Baghdayan AS, Gilmore MS, Shankar N (2004) Enterococcal surface protein, Esp, enhances biofilm formation by *Enterococcus faecalis*. *Infect Immun* 72: 6032-6039. doi:10.1128/IAI.72.10.6032-6039.2004. PubMed: 15385507.

55. Héritier C, Poirel L, Lambert T, Nordmann P (2005) Contribution of acquired carbapenem-hydrolyzing oxacillinases to carbapenem resistance in *Acinetobacter baumannii*. *Antimicrob Agents Chemother* 49: 3198-3202. doi:10.1128/AAC.49.8.3198-3202.2005. PubMed: 16048925.

56. 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. doi:10.1016/j.gene.2008.10.011. PubMed: 19010402.

57. 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. doi:10.1186/1471-2180-10-279. PubMed: 21062436.

58. Soares NC, Cabral MP, Gayoso C, Mallo S, Rodriguez-Velo P et al. (2010) Associating growth-phase-related changes in the proteome of *Acinetobacter baumannii* with increased resistance to oxidative stress. *J Proteome* Res 9: 1951-1964. doi:10.1021/pr901116r. PubMed: 20108952.

59. Petrova OE, Schurr JR, Schurr MJ, Sauer K (2011) The novel *Pseudomonas aeruginosa* two-component regulator BfmR controls bacteriophage-mediated lysis and DNA release during biofilm development through PhdA. *Mol Microbiol* 81: 767-783. doi:10.1111/j. 1365-2958.2011.07733.x. PubMed: 21696457.

60. Liou ML, Soo PC, Ling SR, Kuo HY, Tang CY et al. (2013) The sensor kinase BfmS mediates virulence in *Acinetobacter baumannii*. *J Microbiol Immunol Infect*. PubMed: 23453128. In press

61. McQueary CN, Actis LA (2011) *Acinetobacter baumannii* biofilms: variations among strains and correlations with other cell properties. *J Microbiol* 49: 243-250. doi:10.1007/s12275-011-0343-7. PubMed: 21538245.

62. Gordon NC, Wareham DW (2010) Multidrug-resistant *Acinetobacter baumannii*: mechanisms of virulence and resistance. *Int J Antimicrob Agents* 35: 219-226. doi:10.1016/j.ijantimicag.2009.10.024. PubMed: 20047818.

63. 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. doi:10.1371/journal.pone.0026030. PubMed: 22046254.

64. de Breij A, Gaddy J, van der Meer J, Koning R, Koster A et al. (2009) Csua/BABCDE-dependent pili are not involved in the adherence of *Acinetobacter baumannii* ATCC19606(T) to human airway epithelial cells and their inflammatory response. *Res Microbiol* 160: 213-218. doi: 10.1016/j.resmic.2009.01.002. PubMed: 19530313.

65. Domka J, Lee J, Bansal T, Wood TK (2007) Temporal gene-expression in *Escherichia coli* K-12 biofilms. *Environ Microbiol* 9: 332-346. doi: 10.1111/j.1462-2920.2006.01143.x. PubMed: 17222132.

66. Camilli A, Bassler BL (2006) Bacterial small-molecule signaling pathways. *Science* 311: 1113-1116. doi:10.1126/science.1121357. PubMed: 16497924.

67. Anbazhagan D, Mansor M, Yan GO, Yusof MY, Hassan H, et al. (2012) Detection of quorum sensing signal molecules and identification of an autoinducer synthase gene among biofilm forming clinical isolates of Acinetobacter spp. *PLOS ONE* 7: e36696.

68. Clemmer KM, Bonomo RA, Rather PN (2011) Genetic analysis of surface motility in *Acinetobacter baumannii*. *Microbiology* 157: 2534-2544. doi:10.1099/mic.0.049791-0. PubMed: 21700662.

69. Saroj SD, Rather PN (2013) Streptomycin inhibits quorum sensing in *Acinetobacter baumannii*. *Antimicrob Agents Chemother* 57: 1926-1929. doi:10.1128/AAC.02161-12. PubMed: 23318804.

70. Niu C, Clemmer KM, Bonomo RA, Rather PN (2008) Isolation and characterization of an autoinducer synthase from *Acinetobacter baumannii*. *J Bacteriol* 190: 3386-3392. doi:10.1128/JB.01929-07. PubMed: 18281398.

71. McAllister KA, Peery RB, Zhao G (2006) Acyl carrier protein synthases from gram-negative, gram-positive, and atypical bacterial species: Biochemical and structural properties and physiological implications. *J Bacteriol* 188: 4737-4748. doi:10.1128/JB.01917-05. PubMed: 16788183.

72. Moré MI, Finger LD, Stryker JL, Fuqua C, Eberhard A et al. (1996) Enzymatic synthesis of a quorum-sensing autoinducer through use of defined substrates. *Science* 272: 1655-1658. doi:10.1126/science. 272.5268.1655. PubMed: 8658141.