Comparative Profiling of Pseudomonas aeruginosa Strains Reveals Differential Expression of Novel Unique and Conserved Small RNAs

Abstract

Pseudomonas aeruginosa is a highly adaptable bacterium that thrives in a broad range of ecological niches and can infect multiple hosts as diverse as plants, nematodes and mammals. In humans, it is an important opportunistic pathogen. This wide adaptability correlates with its broad genetic diversity. In this study, we used a deep-sequencing approach to explore the complement of small RNAs (sRNAs) in P. aeruginosa as the number of such regulatory molecules previously identified in this organism is relatively low, considering its genome size, phenotypic diversity and adaptability. We have performed a comparative analysis of PAO1 and PA14 strains which share the same host range but differ in virulence, PA14 being considerably more virulent in several model organisms. Altogether, we have identified more than 150 novel candidate sRNAs and validated a third of them by Northern blotting. Interestingly, a number of these novel sRNAs are strain-specific or showed strain-specific expression, strongly suggesting that they could be involved in determining specific phenotypic traits.

About the Authors:

Silvia Ferrara

Affiliation: Dipartimento di Scienze Biomolecolari e Biotecnologie, Università degli Studi di Milano, Milan, Italy

Margherita Brugnoli

Affiliation: Dipartimento di Scienze Biomolecolari e Biotecnologie, Università degli Studi di Milano, Milan, Italy

Angela De Bonis

Affiliation: Dipartimento di Scienze Biomolecolari e Biotecnologie, Università degli Studi di Milano, Milan, Italy

Francesco Righetti

Affiliation: Dipartimento di Scienze Biomolecolari e Biotecnologie, Università degli Studi di Milano, Milan, Italy

Francesco Delvillani

Affiliation: Dipartimento di Scienze Biomolecolari e Biotecnologie, Università degli Studi di Milano, Milan, Italy

Gianni Dehò

Affiliation: Dipartimento di Scienze Biomolecolari e Biotecnologie, Università degli Studi di Milano, Milan, Italy

David Horner

Affiliation: Dipartimento di Scienze Biomolecolari e Biotecnologie, Università degli Studi di Milano, Milan, Italy

Federica Briani

* E-mail: [email protected] (GB); [email protected] (FB)

Affiliation: Dipartimento di Scienze Biomolecolari e Biotecnologie, Università degli Studi di Milano, Milan, Italy

Giovanni Bertoni

* E-mail: [email protected] (GB); [email protected] (FB)

Affiliation: Dipartimento di Scienze Biomolecolari e Biotecnologie, Università degli Studi di Milano, Milan, Italy

Introduction

Small RNAs (sRNAs) are widespread in bacteria and play critical regulatory roles in several cellular processes [1]–[4]. The prototype of a bacterial sRNA is a non-coding RNA 50-300 nucleotides long that acts by imperfect base pairing with trans-encoded RNA target(s). sRNA-target interaction may lead to modulation of mRNA translation and/or stability [2], [4]. Variations on this theme are also known. For instance, some sRNAs modulate the activity of target proteins or act as mRNAs coding for short proteins. Moreover, there is growing evidence that many sRNAs are cis-encoded and transcribed antisense to their target RNA [5]. The target genes of sRNAs-mediated regulation belong to several different functional groups. The prevalent view is that sRNAs might target almost all bacterial cell processes [6]. In pathogenic microbes, several sRNAs have been shown to be involved in host-microbe interactions and in the adaptation to the host environment [6]. In recent years, genome-scale searches have led to a remarkable increase in the number of identified sRNAs in bacteria [2]. In this context, our knowledge of the sRNA complement of Pseudomonas aeruginosa seemed limited.

P. aeruginosa is a highly adaptable bacterium which thrives in a broad range of ecological niches. In addition, it can infect multiple hosts as diverse as plants, nematodes and mammals. In humans, it is an important opportunistic pathogen in compromised individuals, such as patients with cystic fibrosis, severe burns and impaired immunity [7], [8]. The broad habitat and host ranges of P. aeruginosa reflect the large variety of structural, metabolic and virulence functions found in its pangenome (being 6.2–6.9 Mbp the size range of the sequenced strain genomes) [9]–[12] composed of a high proportion (approximately 90%) of conserved core genes and a rather small accessory genome, found in some strains but not in others, which includes genetic elements supposed to be acquired by horizontal transfer. Accessory genetic elements can confer specific phenotypes that are advantageous under the selective pressure of given habitat or host conditions [10]. Interestingly, a study on the highly virulent strain PA14 has suggested that pathogenicity requires not only virulence factors encoded in the two pathogenicity islands of the accessory genome, but also several core genes [13]. Thus, there seems to be some combinatorial effects between accessory and core functions. In addition, it seems likely that the coordination of the expression of such a panoply of functions is accomplished by regulatory networks based on a large number of regulators. Strikingly, the genome of the archetypal strain PAO1 was found to contain among the highest proportions (9–10%) of regulatory genes as compared to other sequenced bacterial genomes, there being more than 500 genes predicted to encode either transcriptional regulators or two-component regulatory system proteins [11], [12], [14]. In contrast, only a small number (about 40) of regulatory sRNAs have been reported in P. aeruginosa [15] whereas, for example, more than 100 sRNAs have been described in Escherichia coli and Salmonella [1], [3], [16], whose genomes are considerably smaller than P. aeruginosa.

The apparent low proportion of sRNAs in P. aeruginosa could reflect either a real paucity of regulatory sRNAs or the limited number of genome-wide searches that have been performed in this species [17]–[19]. In addition, only few of the sRNAs experimentally validated in P. aeruginosa have been functionally characterized to date; they have been implicated in carbon catabolite repression (CrcZ) [20], in virulence genes expression control (RsmY,Z) [21]–[23], or in other functions that can be important for survival in the infected host, such as iron uptake and storage (PrrF1) [24] and quorum sensing (PhrS) [25]. Finally, despite the variable degree of virulence shown by different P. aeruginosa isolates [13], experimental sRNAs screening has been performed only on PAO1. The identification of genes differentially expressed in virulent vs. attenuated strains, irrespective of whether they belong to core or accessory genome, can be a valuable approach for dissecting pathogenicity in this bacterium. This would be particularly true for genes encoding regulatory factors, such as sRNAs, whose expression level may in turn influence the expression of multiple target genes.

In this work we aimed at the systematic identification of sRNAs of P. aeruginosa by means of the recently developed “sRNA-Seq” approach, an unbiased high-throughput method for the screening of the entire sRNA complement of any organism based on “next-generation” sequencing technologies [26]. We applied the sRNA-Seq method both to PAO1 and to the highly virulent strain PA14, which differ for the presence of about 112 strain-specific gene clusters (54 PAO1-specific and 58 PA14-specific, including the two PA14 pathogenicity islands PAPI-1 and PAPI-2) [13].

By using this approach, we have identified more than 150 novel candidate sRNAs in P. aeruginosa. Interestingly, a relevant number of sRNA hits were strain-specific or showed strain-specific expression, strongly suggesting that they could be involved in determining strain-characteristic phenotypic traits. We probed by Northern blotting 71 candidates and confirmed the expression of 52 new sRNAs, with a validation rate above 73%. Our results expand the panel of P. aeruginosa sRNAs resulting from previous surveys and strongly indicate that the degree of sRNAs utilization as regulators is consistent with other bacterial species.

Methods

RNA Isolation and Generation of sRNA-Seq Amplicon Libraries

Total RNA was prepared from 25 ml samples of early stationary phase (A600 of about 2.6) cultures of P. aeruginosa strains PAO1 [14] and PA14 [27] grown at 37°C in 100 ml of Brain Heart Infusion (BHI) rich medium in 500-ml flasks vigorously shaken (120 rpm). The cells were recovered by centrifugation, resuspended in RNAprotect Cell Reagent (Qiagen) and incubated for 5 min at room temperature, pelleted by centrifugation and stored at −80°C until use. Cells were resuspended in TE–lysozyme (10 mM Tris HCl, 1 mM EDTA, 1 mg/ml lysozyme, pH 7.5), incubated at room temperature for 5 min, and lysed by QIAzol Lysis Reagent (Qiagen). Total RNA was then extracted by the RNeasy Mini Kits (Qiagen) according to the manufacturer’s instructions, including RNase-free DNase I in-column treatment and modifications to enrich for small RNAs (<200 nt). The quality of the RNA was assessed by denaturing (8 M urea) 6% polyacrylamide gel electrophoresis (dPAGE).

Size selection of RNA ranging from 20 to 500 nt was performed by fractionating 160 µg of total RNA on preparative dPAGE and cutting the gel slice containing 20 to 500 nt long transcripts. RNA from 20 to 500 nt was electroeluted from gel slices in a Model 422 Electro-Eluter (Biorad). For the preparation of amplicon libraries, the purified 20–500 nt RNA fraction of each strain was first tagged at the 3′-end with linker L1 (Table S1), a 5′-monophosphate oligonucleotide starting with three ribonucleotides followed by a sequence of 20 deoxyribonucleotides and terminally protected with an inverted dT (IDT, Integrated DNA Technologies). The sequence of this hybrid oligonucloeotide does not match any sequence in the P. aeruginosa genome and is predicted not to form complex secondary structures. 60 µg of RNA was ligated with 78 µg of L1 in T4 RNA ligase Buffer, 10% DMSO, at 16°C with 90 U of T4 RNA ligase (New England Biolabs). After 16 hrs, additional 90 U of T4 RNA ligase was added and incubation prolonged for 8 hrs. To check ligation efficiency, 0.5–1 µg of RNA from the ligation reaction was probed by Northern blotting with [32P]-labelled -oligos AL1 and PA5SRNA02 (Table S1), which probe L1 and 5S rRNA, respectively. To remove non-ligated L1, the ligation mixtures were run on preparative dPAGE and RNA ranging from 40 to 520 nt was electroeluted from gel slices as described above.

[Figure omitted. See PDF.]

Table 1. nstSGR distribution in PAO1 and PA14.

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

[Figure omitted. See PDF.]

Figure 1. Classes of candidate sRNAs and their distribution within the nstSGR groups resulting from sRNA-Seq.

The histogram summarizes the data of Table S2. Candidate sRNAs identified by sRNA-Seq were categorized into five structural/functional classes (I, sRNAs; II, 5′-UTRs; III, asRNAs; IV, CRISPRs; V, sRNAs overlapping annotated ORFs) according to the criteria depicted in Figure S3 and distributed within each nstSGR group (A and B, unique in PAO1 or PA14, respectively; C and E, conserved in both strains but expressed in either PAO1 or PA14, respectively; D, conserved and expressed in both strains).

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

RNase H depletion of tRNA and 5S rRNA was performed as previously described [26] with some modifications. 30 µg of L1-sRNA20–500 was annealed to 9 nmol of Oligo Mix, an equimolar mixture of 47 oligonucleotides (Table S1) complementary to the 3′-ends of P. aeruginosa tRNAs and 5S rRNA. The RNA-DNA hybrids were then digested with RNase H so as to remove the 3′-L1 tail from the small stable RNAs. Depletion efficiency was checked by Northern blotting with [32P]-labelled oligos AL1 and PA5SRNA02. RNA-L1 was then separated from Oligo Mix, tRNA and 5S rRNA degradation products by preparative dPAGE and electroelution from gel slices as described above.

[Figure omitted. See PDF.]

Figure 2. Validation of candidate sRNAs expressed from unique nstSGR in either PAO1 or PA14.

A selection of nstSGRs listed in Table S2, unique in either PAO1 or PA14, were inspected by Northern blot for the expression of sRNAs. Total RNA was extracted from both PAO1 (○) and PA14 (•) grown in the same conditions as for sRNA-Seq. Equal amounts of RNA (8 µg) from both strains were blotted and probed with radiolabelled riboprobes (0002 and 0021) or oligos (Table S1) complementary to nstSGR regions with the highest read coverage, as detailed in Materials and Methods. Validated unique sRNAs in PAO1 or PA14 are shown in (A) and (B), respectively. For SPA0014, 0015, 0018, 0019, signals detected in both strains (dots on the left of PAO1 lanes) can be due to aspecific probe hybridization. The ladder of molecular weight markers is shared by (A) and (B). (nt): nucleotides.

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

[Figure omitted. See PDF.]

Figure 3. Validation of candidate sRNAs differentially expressed from conserved nstSGR.

A selection of conserved nstSGRs listed in Table S2 that were supposed to be differentially expressed between the two strains according to sRNA-Seq data, were inspected by Northern blot. Total RNA was extracted from both PAO1 (○) and PA14 (•) grown in the same conditions as for sRNA-Seq. 8 µg of RNA from both strains were blotted and probed with radiolabelled oligos (Table S1) complementary to nstSGR regions with the highest read coverage. Validated sRNAs which showed higher levels of expression in PAO1 or PA14 are shown in (A) and (B), respectively. The ladder of molecular weight markers is shared by (A) and (B). (nt): nucleotides.

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

[Figure omitted. See PDF.]

Figure 4. Validation of candidate sRNAs similarly expressed from conserved nstSGR.

A selection of conserved nstSGRs listed in Table S2 that were supposed to be similarly expressed between the two strains according to sRNA-Seq data, were inspected by Northern blot. Total RNA was extracted from both PAO1 (○) and PA14 (•) grown in the same conditions as for sRNA-Seq. Equal amounts of RNA (8 µg) from both strains were blotted and probed with radiolabelled oligos or riboprobes (0104, 0112, 0118, 0131, 0143, 0150 and 0157) (Table S1) complementary to nstSGR regions with the highest read coverage. nstSGRs SPA0072, 0085, 0092, and 0122, corresponding to PAO1 loci for the known sRNAs tmRNA, RsmZ, CrcZ, and Spot42, respectively, were included in this analysis as positive controls. 5S RNA served as loading control and molecular weight marker. The ladder of molecular weight markers is indicated on the left of each panel. (nt): nucleotides.

https://doi.org/10.1371/journal.pone.0036553.g004

cDNA20–500 was generated from 1 µg of 5S/tRNA-depleted L1-sRNA20–500 using the SMARTer™ PCR cDNA Synthesis Kit (Clontech), which combines RNA reverse transcription with cDNA 3′-end SMART tailing activity, according to manufacturer’s instruction, except that AL1 oligo (Table S1) was used for reaction priming. cDNA was checked by Southern blotting with [32P]- labelled oligo SmarterII A (Table S1). RNA template was removed by RNase A digestion. To remove free AL1 oligo, the cDNA preparations were run on preparative dPAGE and AL1-cDNA-SMART ranging from 60 to 540 nt (AL1-cDNA20–500-SMART) was electroeluted from gel slices as described above.

[Figure omitted. See PDF.]

Table 2. Candidate sRNAs validated by Northern blot.

https://doi.org/10.1371/journal.pone.0036553.t002

A cDNA20–500-derived amplicon library for 454 pyrosequencing (Roche) was obtained by PCR amplification of AL1-cDNA20–500-SMART using Advantage® 2 PCR polymerase (Clontech) with primers (Table S1) tailored for 454-sequencing with Roche Multiplex Identifiers (MID) for “barcoding”. In particular, MID42 (TCGATCACGT) and MID47 (TGTGAGTAGT) were used to tag amplicons from PAO1 and PA14, respectively. To remove free primers, the PCR reactions were run on preparative dPAGE and amplicons ranging from 130 to 610 nt were electroeluted from gel slices as described above. Amplicons quality and length distribution was checked by Southern blot with [32P]-labelled SmarterII A oligo and by capillary electrophoresis in Agilent Bioanalyzer DNA 7500 (Figure S2A). To enrich RNA transcripts ranging from 130 to 500 nt, a second sRNA-derived amplicon library (Figure S2B) was generated from 5S/tRNA-depleted L1-sRNA ranging in size from 150 to 520 nt (selected by preparative dPAGE and electroelution from gel slices) as described above.

Northern and Southern Blot Analyses

The following procedure was used for both Northern and Southern blot analyses. RNA or DNA samples were heated at 95°C for 5 minutes in loading buffer (5 mM EDTA, 0.025% xylene cyanol, 0.025% bromophenol blue dissolved in formamide) and resolved by dPAGE. Nucleic acids were transferred onto Hybond N+ nylon membranes (GeHealthcare) using a semi-dry electroblotter apparatus (Fastblot B33, Biometra) set at 25 V, 400 mA for 1 hour. The blots were UV-crosslinked and hybridized with [32P]-labelled oligos or riboprobes (Table S1) as described previously [28]. Visualization of radioactive bands was performed by Typhoon™ 8600 variable mode Imager scanner (GE Healthcare BioSciences). All DNA oligonucleotide probes were 5′-end labeled with [γ-32P]ATP and T4 polynucleotide kinase. Riboprobes were prepared as described previously [28] by T7 RNA polymerase transcription of DNA templates obtained by PCR using oligos listed in Table S1 and P. aeruginosa genomic DNA as template. For each validation of candidate sRNA, the probe was complementary to the cognate genomic region with the highest read coverage.

454-pyrosequencing and Data Analysis

Equal amounts of the PAO1 and PA14 amplicon cDNA libraries were combined and submitted to deep-sequencing by a 454 Roche Titanium sequencer using 2/8 of PicoTiterPlate, which should assure at least 140,000 reads. The MID-containing reads were trimmed to eliminate both terminal adaptors, i.e. MIDs, SMART and 454 pyrosequencing primer A-B sequences. Reads were then mapped and clustered throughout the corresponding genome sequence (Genbank accession numbers NC_002516 and NC_008463 for PAO1 and PA14, respectively) as follows. The mapping step was performed using the software SEGEMEHL [29] with default settings but reporting all equal best hits. Mapping positions were considered reliable only if 90% of the read was aligned with ≥90% identity with the genome sequence. Then, the search for read clusters on genome sequences was performed by a sliding window of 200 bp shifted 100 bp at a time along the genome sequence. Significance of clustering of mapped reads was estimated under a null hypothesis of random distribution of reads along the genome using a cumulative Poisson probability. Significant Genomic Regions (SGRs) were defined as consecutive windows where at least one window showed a significant clustering of reads under the Poisson test described above (P≤0.1). SGRs were divided into “structural” (stSGRs, if the genomic annotation reported the keyword “ribosomal” or “tRNA”) and “non-structural” SGRs (nstSGRs, in the other cases). nstSGRs orthology between PAO1 and PA14 strains was determined by reciprocal BLAST. The clustered reads were visualized by GBrowse interface at www.pseudomonas.com database. Sequencing data are accessible at GEO (accession number, GSE36340).

Results

Deep-sequencing of the Low Molecular Weight RNA Fraction of Pseudomonas aeruginosa PAO1 and PA14 Strains

We aimed at sRNA profiling in the P. aeruginosa PAO1 and PA14 strains by sRNA-Seq [26], a massive sequencing approach tailored for unbiased identification of low molecular weight RNA (see Figure S1 for an overview of the procedure). To this end, total RNA was purified from late-exponential cultures of both PAO1 and PA14 strains, respectively, and transcripts ranging from 20 to 500 nt (sRNA20–500) were isolated by gel electrophoresis. The 3′ ends of PAO1 and PA14 sRNA20–500 were tagged by ligation with linker L1, a mixed ribo-deoxyribo-oligonucleotide with its 3′-end protected by an inverted dT (Table S1), obtaining L1-sRNA20–500. sRNA preparations are expected to contain a high proportion of the stable and very abundant tRNAs and 5S rRNA that may interfere with the efficiency of sRNA profiling. We thus selectively degraded the stable RNA component as described previously [26]. Briefly, L1-sRNA20–500 was mixed with a pool of DNA oligos (Table S1) complementary to P. aeruginosa 5S rRNA and tRNAs, and digested with RNase H. Using the SMARTer™ PCR cDNA Synthesis Kit (Clontech), cDNA20–500 was then generated from 5S/tRNA-depleted L1-sRNA20–500 by reverse transcription with an oligonucleotide primer complementary to L1 (AL1; Table S1) and the cDNA 3′-end was tailed with a specific sequence (see Materials and Methods for details). An amplicon library for 454 pyrosequencing (Roche) was then generated by PCR amplification of cDNA20–500 with modular primers complementary to cDNA ends and carrying sequences tailored for 454 sequencing priming and multiplex identification (MID). The PAO1 and PA14 cDNA20–500 amplicons described above were combined in a 1∶1 ratio (amplicon library 1), and submitted to pyrosequencing. This resulted in a raw pool of 101,019 reads (Figure S2A) among which, 0.3% did not show any identifiable linker sequence. The 100,680 linkers-containing reads were examined for MID sequences. 32,156 and 41,514 reads included MID42 (PAO1) and MID47 (PA14) identifiers, respectively, and were at least 17 bases long. After trimming both terminal linker sequences, the reads showed an average length of 34 and 31 nt for PAO1 and PA14, respectively.

As shown in Figure S2A, sRNA molecules longer than 130 nt were poorly represented in this amplicon library. To increase the abundance of longer RNA molecules (corresponding to a read length of about 230 nt in Figure S2A), additional sRNA-derived amplicons were generated for each strain from 150 to 520 nt long RNA fractionated by gel electrophoresis and processed as described above. These PAO1 and PA14 L1-sRNA150–500 amplicons were then combined in a 1∶1 ratio, thus producing a second library (amplicon library 2). The pyrosequencing of the latter resulted in a raw pool of 61,490 reads (Figure S2B), among which 59,132 contained identifiable linker sequences. MID analysis showed that 23,608 and 29,107 reads derived from PAO1 and PA14, respectively. Following terminal trimming, the average read length was about 100 nt for PAO1 and 80 nt for PA14.

Identification of Candidate sRNA Loci and Comparative Analysis between PAO1 and PA14 Strains

Under stringent mapping criteria (>90% read coverage aligned at >90% identity to reference genome), 13,438 and 22,691 reads gave at least one satisfactory match with the genome sequences of PAO1 and PA14, respectively (GenBank accession numbers NC_002516 and NC_008463). A non uniform distribution of reads across the genomes was observed. In fact, more than 99% of genomic positions showed zero coverage, while a limited proportion of sites showed high levels of coverage. To map candidate sRNA loci, genomic regions showing significant reads clustering, hereafter referred to as significant genomic regions (SGRs), were identified as detailed in Materials and Methods. For each strain, about half of the mapped reads clustered in SGRs overlapping stable RNA genes (i.e. tRNAs, 5S rRNA); these were classified as structural SGRs (stSGRs) and not included in further analysis. Around 90% of the remaining mapped reads fell in other significant clusters (non-structural SGRs, nstSGRs), whereas about 10% were not clustered. As stable RNAs are expected to be much more abundant than other RNAs, the observed 1∶1 ratio between the number of reads mapping in stSGRs over nstSGR reads indicates the high efficiency of tRNAs and 5S rRNA depletion achieved in amplicon library preparation.

As a whole, we defined 128 and 168 nstSGRs in PAO1 and PA14 genomes, respectively (Table 1 and Table S2) mapping within different loci: i) genes for housekeeping RNAs (tmRNA, 6S, 4.5S and RNase P RNAs); ii) genes for sRNAs previously identified in PAO1 (14) (Table S3), and iii) both intergenic and intragenic loci not previously known to express sRNAs (201).

By reciprocal BLAST, we determined whether the identified loci were conserved or not in the two strains. Both unique (in either strain) and conserved loci were found. Therefore, the corresponding nstSGRs were classified in 5 groups as shown in Table 1. Group A and B nstSGRs map within loci unique to PAO1 and PA14, respectively; groups C, D and E include nstSGRs mapping in conserved loci. Group D nstSGRs were found in both strains, whereas group C and E nstSGRs were identified only in PAO1 and PA14, respectively. Thus, the comparative profiling of sRNAs from PAO1 and PA14 suggested the existence of both strain-specific (groups A and B) and conserved candidate sRNA loci; the latter, in a number of cases (groups C and E), appeared differentially expressed in the two strains.

Within each group described above, we classified the candidate sRNAs according to functional/structural categories established for regulatory RNAs in bacteria [1] as follows (Figure S3). Class I groups nstSGRs located in intergenic regions (>30 nt from flanking ORFs). Trans-encoded sRNAs (sRNA) would belong to this class; class II groups nstSGRs with read clustering spanning 5′-untranslated regions (5′-UTRs) in sense orientation. This class would encompass mRNA riboswitches and sRNAs generated by mRNA transcription attenuation or processing; class III includes nstSGRs with intragenic (<30 nt from flanking ORFs) reads clustering in antisense orientation. Cis-encoded antisense sRNAs (asRNAs) would cluster in this class; class IV groups intergenic nstSGRs containing CRISPR-like arrays [30]; finally, nstSGRs with read clustering within ORFs and/or 3′-UTRs in sense orientation belong to class V.

The results of this analysis are summarized in Figure 1 and details of each nstSGRs are listed in Table S2. Since class V nstSGRs may correspond to stable mRNA degradation fragments, whose regulatory role is uncertain, they were excluded from further analysis and not reported in Table S2, with the exception of nstSGRs encompassing small putative ORFs.

Remarkably 19 hits of Table S2 corresponded to members of the panel of about 40 sRNAs previously identified in PAO1 [15] including sRNAs annotated in the Pseudomonas genome databasev2 (www.pseudomonas.com) such as the housekeeping tmRNA, 6S, 4.5S and RNase P RNAs, and sRNAs already characterized such as CrcZ, RsmY, RsmZ, PhrS and AmiL [15] and a putative Spot42 sRNA (SPA0122) which is located in a conserved genomic context in E. coli, Salmonella and pseudomonads [31] (Table S3). We show here that this panel of known sRNAs previously detected in PAO1 is comparably expressed in PA14. Many of the previously identified sRNAs that escaped our analysis have been reported to be expressed at low level or in response to environmental stimuli (e.g. iron limitation for Prrf1 and 2) [24]. However, in a recent deeper transcriptomic survey of PA14 [32] all known P. aeruginosa sRNAs were detected. Therefore, it is possible that in our sRNA-seq approach we missed scarcely expressed sRNAs.

Taken together, the data described above, subtracted from those sRNAs already known in P. aeruginosa, represent a panel of 163 novel sRNA candidates.

Validation by Northern-blot Analysis of sRNAs Expression from nstSGRs

We tested by Northern blotting the expression of a sample of 71 novel candidates covering all groups and classes (Table S2). Our sample for validation was not random, as we gave priority to strain-specific candidates for validation, but disregarding those having features typical of antisense sRNA regulating a transposase genes (see below). Moreover, we favored class I and III candidates, i.e. trans-encoded sRNAs and asRNAs, respectively. In particular, we analyzed 29 class I, 12 class II, 29 class III and 1 class V nstSGRs throughout the A-E groups. The previously identified sRNAs RsmZ (SPA0085) [22], CrcZ (SPA0092) [20], Spot42 (SPA0122) [31] and tmRNA (SPA0072) [19] were used as positive controls. Moreover, a class IV nstSGR, corresponding to CRISPR-2 [33] was included in this validation panel.

Out of 71 novel candidates tested, 52 showed signals in Northern blot experiments (Figures 2, 3 and 4); for 19 we could detect only very faint signals, barely above the background, or no signal (data not shown). Thus, the validation rate was above 73%. Among the validated sRNAs (Table 2), 22 belonged to class I (sRNA), 19 to class III (asRNA), 10 to class II (5′-UTR), 1 to class V. The expression of CRISPR-2 (class IV) was validated and a major band, corresponding to processed crRNA [34], was observed.

The majority of sRNAs tested, which were expected to be equally expressed in both strains (group D) (Table S2), showed signals whose intensity in the two strains was consistent with sRNA-Seq data (i.e. the read number of the corresponding nstSGRs). One exception was SPA0101 (Figure 3B) which showed a comparable read number in both strains, but gave a sharp signal corresponding to a ∼70 nt long transcript only in PA14. On the contrary, many group E sRNAs tested (Table S2), whose corresponding nstSGRs displayed expression in PA14 only in sRNA-Seq, showed comparable expression in the two strains in Northern blot analysis. However, all the corresponding nstSGRs of these sRNAs were identified in PA14 by a read number at best slightly above the significance threshold (from 3 to 6). Thus, stochastic fluctuations in amplicon library preparation may have been sufficient to keep the read number below the threshold in PAO1. On the whole, we validated 13 novel unique sRNAs (Figures 2A and B), 30 conserved sRNA with comparable expression in both strains (Figure 4), and 9 conserved sRNAs showing differential expression (Figure 3A and B).

In most cases, we found that transcript size predicted by sRNASeq (Table 2) corresponded to the strongest Northern blotting signal (Figures 2, 3, and 4). The sRNA-Seq reads were scattered within two class I (SPA0015 and SPA0074) and several class II and class III nstSGRs (SPA0023, 0112, 0118, 0150 and 0156). Accordingly, these nstSGRs showed complex patterns with multiple signals of comparable intensity. Degradation of unstable primary transcripts by cellular nucleases may explain these results. However, the presence of scattered reads within an nstSGR was observed also for other sRNAs such as CRISPR2, tmRNA, CrcZ and SPA102, for which a major signal was clearly visible by Northern hybridization. Moreover, in some cases (i.e. SPA0011, 0013, 0018, 0157, 0167, 0168), the regions covered by the reads were smaller than the observed transcripts. As samples preparation for sRNA-Seq included an RNase H digestion step, unspecific RNA degradation by this enzyme may account for these results.

Discussion

Increasing Complexity of P. aeruginosa RNA World by sRNA-Seq

In this work we have performed a parallel sRNAs search in P. aeruginosa by sRNA-Seq, a powerful unbiased method that allows the analysis by deep sequencing of the whole small transcriptome (i.e. both primary and processed transcripts) [26]. Unlike previous surveys performed in PAO1, our search for sRNA loci was not biased by a priori assumptions about sRNA-based regulation mechanisms, such as binding by Hfq [35], whose role in sRNA-mediated regulation system is not clearly established in P. aeruginosa, or genetic features putatively associated with sRNA-coding loci (e.g. mapping within intergenic regions with predicted promoters and terminators), which were employed in previous bioinformatics-based analyses [17]–[19]. The first goal of our analysis was to expand the P. aeruginosa sRNA panel resulting from previous surveys in terms of both amplitude and sRNA typologies (potential antisense RNAs for example were completely disregarded by previous analyses) [15], [17]–[19]. Moreover, we did not overlook putative “bifunctional” sRNAs, such as short transcripts encompassing 5′-UTRs or encoding small peptides [3], [36].

Our approach resulted in the definition of 163 loci expressing new candidate sRNAs. We found a comparable number of class I (sRNA) and III (asRNA) sRNAs, which altogether accounted for more than 75% of our sRNAs panel (Figure 1). In addition, several (34/181) class II sRNAs (mapping within 5′-UTR) were found. These short transcripts, also identified in previous genome-wide searches for sRNAs [37], [38], can be generated by premature transcription termination, or 5′-UTR processing as by-products of post-transcriptional gene regulation or mRNA degradation. However, in some cases they can also act as trans-encoded sRNAs. In fact, it has recently been reported that two S-adenosylmethionine riboswitches of Listeria monocytogenes, SreA and SreB, can base pair with the mRNA of prfA, the master regulator of Listeria virulence, and repress its expression [39].

Finally, both the annotated CRISPR of PA14 [33] and ten class V loci for potential peptide-coding sRNAs were detected (Figure 1). These latter sRNAs may have the dual status of short mRNAs (encoding low molecular weight proteins) and trans-acting sRNAs, as it has been established for the E. coli SgrS and the Staphylococcus aureus RNAIII sRNAs [40]–[42].

We assayed the expression of a large sample of candidate sRNAs by Northern blotting. Remarkably, the expression of many ntsSGRs defined by a read number only slightly above the significance threshold (e.g. SPA0102, 0110, 0112 and 0156) could be demonstrated by Northern blotting. Furthermore, the majority of validated sRNAs showed expression levels in Northern assays that were consistent with sRNA-Seq analysis. Thus sRNA-Seq not only appears a sensitive approach to sRNA identification but could also represent a reliable method for estimating their expression levels in comparative analyses.

On the whole, we could validate the expression of 52 novel sRNAs, more than doubling the number of P. aeruginosa sRNAs annotated so far. Interestingly, several validated 5′-UTR nstSGRs (e.g. SPA0101-0104) showed one or two sharp signals in Northern blotting experiments (Figs. 2, 3, and 4) corresponding to discrete RNA species and may thus be good candidates for trans-acting sRNAs, as mentioned above. Overall, our data significantly increase the complexity of sRNA complement in P. aeruginosa and suggest that RNA-mediated regulation in this organism may be as common and multifaceted as it is in other bacteria [1], [3].

sRNA-mediated Regulation May Contribute to Pseudomonas Strain-specific Phenotypic Traits

Another purpose of our work was to get hints on sRNA-mediated regulatory mechanisms possibly involved in strain-specific phenotypic traits such as pathogenicity and virulence. To this aim, we performed a comparative analysis of PAO1 and PA14 strains that, although sharing the same host range, differ in virulence, being PA14 considerably more virulent in several model organisms [43].

26 nstSGRs identified by sRNA-Seq consisted of unique loci in either PAO1 or PA14 (groups A and B, respectively; Table S2). In PA14, these loci mostly mapped within regions of genome plasticity (RGPs, defined as polymorphic strain-specific segments encompassing at least 4 contiguous ORFs) [9], with SPA0016 representing the only exception (Table S2). As for the 9 nstSGRs unique to PAO1, 2 mapped in RGPs (SPA0001, which corresponded to the already known PhrD sRNA [15], [19] and SPA0003. Remarkably, 6 overlapped in antisense orientation the 5′-UTRs of a gene encoding a putative transposase of the IS116/IS110/IS902 family (SPA0004-8 and SPA0066). This gene is identically repeated six times in PAO1 genome and sRNA-Seq reads were randomly distributed by the mapping software among the six loci. Transposase translation regulation by antisense RNAs has been extensively studied in the IS10 system, where a short RNA (RNA-OUT), which is transcribed in antisense orientation from the 5′-end of the transposase locus, interacts with transposase mRNA to hinder ribosome binding site [44], [45]. We did not check by Northern blotting the expression of this PAO1 putative transposase antisense RNA; however, the high overall sRNA-Seq read number (899; Table S2) suggests that it can be actively transcribed and thus may play a role in transposase regulation. Another sRNA antisense to a transposase gene could be expressed by the SPA0022 locus, which maps within the PAPI-1 pathogenicity island unique to PA14 and encodes a polypeptide belonging to IS66 OrfC family [46]. However, in this case only few reads were detected by sRNA-Seq that mapped within the 3′-end region of the transposase.

We validated the expression of 13 novel strain-specific sRNAs by Northern blotting, 2 unique to PAO1 and 11 to PA14. Interestingly, three PA14 novel sRNAs, SPA0015, SPA0021 and SPA0023, fall within the pathogenicity island PAPI-1. SPA0015 locus maps in the intergenic region between genes RL003 (PA14_60130), encoding a homolog of DNA relaxase [47], and RL004 (PA14_60120; dcd2), encoding a putative deoxycytidine deaminase. RL003 mutants showed virulence-attenuated phenotype [48], and reduced efficiency in PAPI-1 horizontal transfer [47]. SPA0021 and SPA0023 asRNAs are cis-encoded antisense to RL076 (PA14_59370) and RL033 (PA14_59840) genes, respectively, both encoding hypothetical proteins. An insertion mutation in RL076 reduced the efficiency of PAPI-1 horizontal transfer [47], whereas an RL033 mutant showed attenuated virulence [48]. It will be interesting to assess whether sRNA-mediated regulation at these loci may be involved in PA14 virulence.

Most nstSGRs (155/181; Table S2) mapped in conserved loci and were identified by sRNA-Seq in both strains (group D) or in either one (groups C and E). At the pangenome level, these conserved loci mostly belong to the core genome, which constitutes approximately 90% of the total genome and is highly conserved in all strains analyzed so far [9], [10]. Only three conserved loci, SPA0097, SPA0169 and SPA182 belong to the accessory genome and mapped within RGPs (Table S2). Out of the 52 novel validated sRNAs loci, 30 belong to the core genome and exhibited comparable expression in PAO1 and PA14 (Figure 4). On the contrary, 8 sRNAs belonging to the core genome and 1 to RPG43 showed differential expression between the two strains, being more highly expressed either in PAO1 (Figure 3A) or in PA14 (Figure 3B). It has been pointed out that for sRNAs with multiple targets, a hierarchy in target binding due to sRNA-mRNA interaction strength may exist [49]. Thus, it is possible that differential expression of a sRNA may result not only in quantitative differences in the strength of target mRNA(s) regulation, but also ultimately change the number of mRNA species targeted by the same sRNA. It will be interesting to assess whether this is indeed the case for P. aeruginosa differentially expressed sRNAs.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Steps of the comparative analysis of the small transcriptome of PAO1 and PA14 strains. From RNA extraction to sRNAs verification, the sequence of steps followed for the comparative analysis of the small transcriptome of the strains PAO1 and PA14 is depicted as a flow chart. The whole approach was performed in three phases: (A) RNA preparation and 454 pyrosequencing; (B) Deep-sequencing data analysis; (C) Comparative PAO1 vs PA14 analysis.

https://doi.org/10.1371/journal.pone.0036553.s001

(TIF)

Figure S2.

Read length distribution of amplicon libraries. 454 pyrosequencing results in terms of length distribution of untrimmed reads are shown for amplicon library 1 (A) and 2 (B). Terminal adaptors of amplicons for pyrosequencing are altogether 100 nt long. Note the enrichment in (B) of reads with actual length longer than 130 nt, which are scarcely represented in (A).

https://doi.org/10.1371/journal.pone.0036553.s002

(TIF)

Figure S3.

Criteria for categorization of nstSGRs into classes. nstSGRs are represented by thick black arrows. By way of example, the upper nstSGR bears on top the read cluster by which nstSGR has been defined. Grey tip-ended segments represent annotated ORFs located farther than 30 nt from the nearest end of read cluster. Orange tip-ended segments represent annotated ORFs located within 30 nt from or overlapping with (ol) at least one end of the read cluster. Class I (intergenic), II (5′-UTRs), III (antisense) and V (intragenic) nstSGRs differ for their relative positioning to annotated ORFs. Class IV nstSGRs corresponding to CRISPR-like array are not depicted in this figure.

https://doi.org/10.1371/journal.pone.0036553.s003

(TIF)

Table S1.

Oligonucleotides.

https://doi.org/10.1371/journal.pone.0036553.s004

(PDF)

Table S2.

Compilation of nstSGRs identified by parallel sRNA-Seq approach in the P. aeruginosa strains PAO1 and PA14.

https://doi.org/10.1371/journal.pone.0036553.s005

(PDF)

Table S3.

Previously identified P. aeruginosa sRNAs found in this work.

https://doi.org/10.1371/journal.pone.0036553.s006

(PDF)

Acknowledgments

The authors acknowledge all members of the lab for helpful discussion and technical support, F. Dal Pero and A. Albiero for 454-pyrosequencing and data analysis, respectively, and G. Pavesi and M. Borsani for L1 oligonucleotide design.

Author Contributions

Conceived and designed the experiments: SF FB GB. Performed the experiments: SF MB ADB FR FD. Analyzed the data: SF MB ADB FR FD GD DH FB GB. Wrote the paper: SF GD DH FB GB. Designed the software used in this analysis: DH. Conceived the project: GB.

Citation: Ferrara S, Brugnoli M, De Bonis A, Righetti F, Delvillani F, Dehò G, et al. (2012) Comparative Profiling of Pseudomonas aeruginosa Strains Reveals Differential Expression of Novel Unique and Conserved Small RNAs. PLoS ONE7(5): e36553. https://doi.org/10.1371/journal.pone.0036553

References

1. Waters LS, Storz G (2009) Regulatory RNAs in bacteria. Cell 136: 615–628.LS WatersG. Storz2009Regulatory RNAs in bacteria.Cell136615628

2. Liu JM, Camilli A (2010) A broadening world of bacterial small RNAs. Curr Opin Microbiol 13: 18–23.JM LiuA. Camilli2010A broadening world of bacterial small RNAs.Curr Opin Microbiol131823

3. Gottesman S, Storz G (2010) Bacterial small RNA regulators: versatile roles and rapidly evolving variations. Cold Spring Harb Perspect Biol. S. GottesmanG. Storz2010Bacterial small RNA regulators: versatile roles and rapidly evolving variations.Cold Spring Harb Perspect Biol

4. Storz G, Vogel J, Wassarman KM (2011) Regulation by small RNAs in bacteria: expanding frontiers. Mol Cell 43: 880–891.G. StorzJ. VogelKM Wassarman2011Regulation by small RNAs in bacteria: expanding frontiers.Mol Cell43880891

5. Georg J, Hess WR (2011) cis-antisense RNA, another level of gene regulation in bacteria. Microbiol Mol Biol Rev 75: 286–300.J. GeorgWR Hess2011cis-antisense RNA, another level of gene regulation in bacteria.Microbiol Mol Biol Rev75286300

6. Papenfort K, Vogel J (2010) Regulatory RNA in bacterial pathogens. Cell Host Microbe 8: 116–127.K. PapenfortJ. Vogel2010Regulatory RNA in bacterial pathogens.Cell Host Microbe8116127

7. Pier GB, Ramphal R (2005) Pseudomonas aeruginosa. In: Mandell GL, Bennett JE, Dolin R, editors. Principles and practice of infectious diseases. Philadelphia, PA: Elsevier Churchill Livingstone. pp. 2587–2615.GB PierR. Ramphal2005Pseudomonas aeruginosa.GL MandellJE BennettR. DolinPrinciples and practice of infectious diseasesPhiladelphia, PAElsevier Churchill Livingstone25872615

8. Wagner VE, Filiatrault MJ, Picardo KF, Iglewski BH (2008) Pseudomonas aeruginosa virulence and pathogenesis issues. In: Cornelis P, editor. Pseudomonas genomics and molecular biology. Norfolk: Caister Academic Press. pp. 129–158.VE WagnerMJ FiliatraultKF PicardoBH Iglewski2008Pseudomonas aeruginosa virulence and pathogenesis issues.P. CornelisPseudomonas genomics and molecular biologyNorfolkCaister Academic Press129158

9. Mathee K, Narasimhan G, Valdes C, Qiu X, Matewish JM, et al. (2008) Dynamics of Pseudomonas aeruginosa genome evolution. Proc Natl Acad Sci U S A 105: 3100–3105.K. MatheeG. NarasimhanC. ValdesX. QiuJM Matewish2008Dynamics of Pseudomonas aeruginosa genome evolution.Proc Natl Acad Sci U S A10531003105

10. Kung VL, Ozer EA, Hauser AR (2010) The accessory genome of Pseudomonas aeruginosa. Microbiol Mol Biol Rev 74: 621–641.VL KungEA OzerAR Hauser2010The accessory genome of Pseudomonas aeruginosa.Microbiol Mol Biol Rev74621641

11. Klockgether J, Cramer N, Wiehlmann L, Davenport CF, Tummler B (2011) Pseudomonas aeruginosa genomic structure and diversity. Front Microbiol 2: 150.J. KlockgetherN. CramerL. WiehlmannCF DavenportB. Tummler2011Pseudomonas aeruginosa genomic structure and diversity.Front Microbiol2150

12. Silby MW, Winstanley C, Godfrey SA, Levy SB, Jackson RW (2011) Pseudomonas genomes: diverse and adaptable. FEMS Microbiol Rev 35: 652–680.MW SilbyC. WinstanleySA GodfreySB LevyRW Jackson2011Pseudomonas genomes: diverse and adaptable.FEMS Microbiol Rev35652680

13. Lee DG, Urbach JM, Wu G, Liberati NT, Feinbaum RL, et al. (2006) Genomic analysis reveals that Pseudomonas aeruginosa virulence is combinatorial. Genome Biol 7: R90.DG LeeJM UrbachG. WuNT LiberatiRL Feinbaum2006Genomic analysis reveals that Pseudomonas aeruginosa virulence is combinatorial.Genome Biol7R90

14. Stover CK, Pham XQ, Erwin AL, Mizoguchi SD, Warrener P, et al. (2000) Complete genome sequence of Pseudomonas aeruginosa PAO1, an opportunistic pathogen. Nature 406: 959–964.CK StoverXQ PhamAL ErwinSD MizoguchiP. Warrener2000Complete genome sequence of Pseudomonas aeruginosa PAO1, an opportunistic pathogen.Nature406959964

15. Sonnleitner E, Haas D (2011) Small RNAs as regulators of primary and secondary metabolism in Pseudomonas species. Appl Microbiol Biotechnol 91: 63–79.E. SonnleitnerD. Haas2011Small RNAs as regulators of primary and secondary metabolism in Pseudomonas species.Appl Microbiol Biotechnol916379

16. Vogel J (2009) A rough guide to the non-coding RNA world of Salmonella. Mol Microbiol 71: 1–11.J. Vogel2009A rough guide to the non-coding RNA world of Salmonella.Mol Microbiol71111

17. Livny J, Brencic A, Lory S, Waldor MK (2006) Identification of 17 Pseudomonas aeruginosa sRNAs and prediction of sRNA-encoding genes in 10 diverse pathogens using the bioinformatic tool sRNAPredict2. Nucleic Acids Res 34: 3484–3493.J. LivnyA. BrencicS. LoryMK Waldor2006Identification of 17 Pseudomonas aeruginosa sRNAs and prediction of sRNA-encoding genes in 10 diverse pathogens using the bioinformatic tool sRNAPredict2.Nucleic Acids Res3434843493

18. Gonzalez N, Heeb S, Valverde C, Kay E, Reimmann C, et al. (2008) Genome-wide search reveals a novel GacA-regulated small RNA in Pseudomonas species. BMC Genomics 9: 167.N. GonzalezS. HeebC. ValverdeE. KayC. Reimmann2008Genome-wide search reveals a novel GacA-regulated small RNA in Pseudomonas species.BMC Genomics9167

19. Sonnleitner E, Sorger-Domenigg T, Madej MJ, Findeiss S, Hackermuller J, et al. (2008) Detection of small RNAs in Pseudomonas aeruginosa by RNomics and structure-based bioinformatic tools. Microbiology 154: 3175–3187.E. SonnleitnerT. Sorger-DomeniggMJ MadejS. FindeissJ. Hackermuller2008Detection of small RNAs in Pseudomonas aeruginosa by RNomics and structure-based bioinformatic tools.Microbiology15431753187

20. Sonnleitner E, Abdou L, Haas D (2009) Small RNA as global regulator of carbon catabolite repression in Pseudomonas aeruginosa. Proc Natl Acad Sci U S A 106: 21866–21871.E. SonnleitnerL. AbdouD. Haas2009Small RNA as global regulator of carbon catabolite repression in Pseudomonas aeruginosa.Proc Natl Acad Sci U S A1062186621871

21. Kay E, Humair B, Denervaud V, Riedel K, Spahr S, et al. (2006) Two GacA-dependent small RNAs modulate the quorum-sensing response in Pseudomonas aeruginosa. J Bacteriol 188: 6026–6033.E. KayB. HumairV. DenervaudK. RiedelS. Spahr2006Two GacA-dependent small RNAs modulate the quorum-sensing response in Pseudomonas aeruginosa.J Bacteriol18860266033

22. Heurlier K, Williams F, Heeb S, Dormond C, Pessi G, et al. (2004) Positive control of swarming, rhamnolipid synthesis, and lipase production by the posttranscriptional RsmA/RsmZ system in Pseudomonas aeruginosa PAO1. J Bacteriol 186: 2936–2945.K. HeurlierF. WilliamsS. HeebC. DormondG. Pessi2004Positive control of swarming, rhamnolipid synthesis, and lipase production by the posttranscriptional RsmA/RsmZ system in Pseudomonas aeruginosa PAO1.J Bacteriol18629362945

23. Brencic A, McFarland KA, McManus HR, Castang S, Mogno I, et al. (2009) The GacS/GacA signal transduction system of Pseudomonas aeruginosa acts exclusively through its control over the transcription of the RsmY and RsmZ regulatory small RNAs. Mol Microbiol 73: 434–445.A. BrencicKA McFarlandHR McManusS. CastangI. Mogno2009The GacS/GacA signal transduction system of Pseudomonas aeruginosa acts exclusively through its control over the transcription of the RsmY and RsmZ regulatory small RNAs.Mol Microbiol73434445

24. Wilderman PJ, Sowa NA, FitzGerald DJ, FitzGerald PC, Gottesman S, et al. (2004) Identification of tandem duplicate regulatory small RNAs in Pseudomonas aeruginosa involved in iron homeostasis. Proc Natl Acad Sci U S A 101: 9792–9797.PJ WildermanNA SowaDJ FitzGeraldPC FitzGeraldS. Gottesman2004Identification of tandem duplicate regulatory small RNAs in Pseudomonas aeruginosa involved in iron homeostasis.Proc Natl Acad Sci U S A10197929797

25. Sonnleitner E, Gonzalez N, Sorger-Domenigg T, Heeb S, Richter AS, et al. (2011) The small RNA PhrS stimulates synthesis of the Pseudomonas aeruginosa quinolone signal. Mol Microbiol 80: 868–885.E. SonnleitnerN. GonzalezT. Sorger-DomeniggS. HeebAS Richter2011The small RNA PhrS stimulates synthesis of the Pseudomonas aeruginosa quinolone signal.Mol Microbiol80868885

26. Liu JM, Livny J, Lawrence MS, Kimball MD, Waldor MK, et al. (2009) Experimental discovery of sRNAs in Vibrio cholerae by direct cloning, 5S/tRNA depletion and parallel sequencing. Nucleic Acids Res 37: e46.JM LiuJ. LivnyMS LawrenceMD KimballMK Waldor2009Experimental discovery of sRNAs in Vibrio cholerae by direct cloning, 5S/tRNA depletion and parallel sequencing.Nucleic Acids Res37e46

27. Rahme LG, Stevens EJ, Wolfort SF, Shao J, Tompkins RG, et al. (1995) Common virulence factors for bacterial pathogenicity in plants and animals. Science 268: 1899–1902.LG RahmeEJ StevensSF WolfortJ. ShaoRG Tompkins1995Common virulence factors for bacterial pathogenicity in plants and animals.Science26818991902

28. Briani F, Zangrossi S, Ghisotti D, Deho G (1996) A Rho-dependent transcription termination site regulated by bacteriophage P4 RNA immunity factor. Virology 223: 57–67.F. BrianiS. ZangrossiD. GhisottiG. Deho1996A Rho-dependent transcription termination site regulated by bacteriophage P4 RNA immunity factor.Virology2235767

29. Hoffmann S, Otto C, Kurtz S, Sharma CM, Khaitovich P, et al. (2009) Fast mapping of short sequences with mismatches, insertions and deletions using index structures. PLoS Comput Biol 5: e1000502.S. HoffmannC. OttoS. KurtzCM SharmaP. Khaitovich2009Fast mapping of short sequences with mismatches, insertions and deletions using index structures.PLoS Comput Biol5e1000502

30. Sorek R, Kunin V, Hugenholtz P (2008) CRISPR–a widespread system that provides acquired resistance against phages in bacteria and archaea. Nat Rev Microbiol 6: 181–186.R. SorekV. KuninP. Hugenholtz2008CRISPR–a widespread system that provides acquired resistance against phages in bacteria and archaea.Nat Rev Microbiol6181186

31. Gottesman S, McCullen CA, Guillier M, Vanderpool CK, Majdalani N, et al. (2006) Small RNA regulators and the bacterial response to stress. Cold Spring Harb Symp Quant Biol 71: 1–11.S. GottesmanCA McCullenM. GuillierCK VanderpoolN. Majdalani2006Small RNA regulators and the bacterial response to stress.Cold Spring Harb Symp Quant Biol71111

32. Dotsch A, Eckweiler D, Schniederjans M, Zimmermann A, Jensen V, et al. (2012) The Pseudomonas aeruginosa transcriptome in planktonic cultures and static biofilms using RNA sequencing. PLoS One 7: e31092.A. DotschD. EckweilerM. SchniederjansA. ZimmermannV. Jensen2012The Pseudomonas aeruginosa transcriptome in planktonic cultures and static biofilms using RNA sequencing.PLoS One7e31092

33. Zegans ME, Wagner JC, Cady KC, Murphy DM, Hammond JH, et al. (2009) Interaction between bacteriophage DMS3 and host CRISPR region inhibits group behaviors of Pseudomonas aeruginosa. J Bacteriol 191: 210–219.ME ZegansJC WagnerKC CadyDM MurphyJH Hammond2009Interaction between bacteriophage DMS3 and host CRISPR region inhibits group behaviors of Pseudomonas aeruginosa.J Bacteriol191210219

34. Cady KC, O’Toole GA (2011) Non-identity-mediated CRISPR-bacteriophage interaction mediated via the Csy and Cas3 proteins. J Bacteriol 193: 3433–3445.KC CadyGA O’Toole2011Non-identity-mediated CRISPR-bacteriophage interaction mediated via the Csy and Cas3 proteins.J Bacteriol19334333445

35. Vogel J, Luisi BF (2011) Hfq and its constellation of RNA. Nat Rev Microbiol 9: 578–589.J. VogelBF Luisi2011Hfq and its constellation of RNA.Nat Rev Microbiol9578589

36. Vanderpool CK, Balasubramanian D, Lloyd CR (2011) Dual-function RNA regulators in bacteria. Biochimie 93: 1943–1949.CK VanderpoolD. BalasubramanianCR Lloyd2011Dual-function RNA regulators in bacteria.Biochimie9319431949

37. Vogel J, Bartels V, Tang TH, Churakov G, Slagter-Jager JG, et al. (2003) RNomics in Escherichia coli detects new sRNA species and indicates parallel transcriptional output in bacteria. Nucleic Acids Res 31: 6435–6443.J. VogelV. BartelsTH TangG. ChurakovJG Slagter-Jager2003RNomics in Escherichia coli detects new sRNA species and indicates parallel transcriptional output in bacteria.Nucleic Acids Res3164356443

38. Kawano M, Reynolds AA, Miranda-Rios J, Storz G (2005) Detection of 5′- and 3′-UTR-derived small RNAs and cis-encoded antisense RNAs in Escherichia coli. Nucleic Acids Res 33: 1040–1050.M. KawanoAA ReynoldsJ. Miranda-RiosG. Storz2005Detection of 5′- and 3′-UTR-derived small RNAs and cis-encoded antisense RNAs in Escherichia coli.Nucleic Acids Res3310401050

39. Loh E, Dussurget O, Gripenland J, Vaitkevicius K, Tiensuu T, et al. (2009) A trans-acting riboswitch controls expression of the virulence regulator PrfA in Listeria monocytogenes. Cell 139: 770–779.E. LohO. DussurgetJ. GripenlandK. VaitkeviciusT. Tiensuu2009A trans-acting riboswitch controls expression of the virulence regulator PrfA in Listeria monocytogenes.Cell139770779

40. Novick RP, Ross HF, Projan SJ, Kornblum J, Kreiswirth B, et al. (1993) Synthesis of staphylococcal virulence factors is controlled by a regulatory RNA molecule. EMBO J 12: 3967–3975.RP NovickHF RossSJ ProjanJ. KornblumB. Kreiswirth1993Synthesis of staphylococcal virulence factors is controlled by a regulatory RNA molecule.EMBO J1239673975

41. Wadler CS, Vanderpool CK (2007) A dual function for a bacterial small RNA: SgrS performs base pairing-dependent regulation and encodes a functional polypeptide. Proc Natl Acad Sci U S A 104: 20454–20459.CS WadlerCK Vanderpool2007A dual function for a bacterial small RNA: SgrS performs base pairing-dependent regulation and encodes a functional polypeptide.Proc Natl Acad Sci U S A1042045420459

42. Chevalier C, Boisset S, Romilly C, Masquida B, Fechter P, et al. (2010) Staphylococcus aureus RNAIII binds to two distant regions of coa mRNA to arrest translation and promote mRNA degradation. PLoS Pathog 6: e1000809.C. ChevalierS. BoissetC. RomillyB. MasquidaP. Fechter2010Staphylococcus aureus RNAIII binds to two distant regions of coa mRNA to arrest translation and promote mRNA degradation.PLoS Pathog6e1000809

43. Rahme LG, Ausubel FM, Cao H, Drenkard E, Goumnerov BC, et al. (2000) Plants and animals share functionally common bacterial virulence factors. Proc Natl Acad Sci U S A 97: 8815–8821.LG RahmeFM AusubelH. CaoE. DrenkardBC Goumnerov2000Plants and animals share functionally common bacterial virulence factors.Proc Natl Acad Sci U S A9788158821

44. Simons RW, Kleckner N (1983) Translational control of IS10 transposition. Cell 34: 683–691.RW SimonsN. Kleckner1983Translational control of IS10 transposition.Cell34683691

45. Ma C, Simons RW (1990) The IS10 antisense RNA blocks ribosome binding at the transposase translation initiation site. EMBO J 9: 1267–1274.C. MaRW Simons1990The IS10 antisense RNA blocks ribosome binding at the transposase translation initiation site.EMBO J912671274

46. Gourbeyre E, Siguier P, Chandler M (2010) Route 66: investigations into the organisation and distribution of the IS66 family of prokaryotic insertion sequences. Res Microbiol 161: 136–143.E. GourbeyreP. SiguierM. Chandler2010Route 66: investigations into the organisation and distribution of the IS66 family of prokaryotic insertion sequences.Res Microbiol161136143

47. Carter MQ, Chen J, Lory S (2010) The Pseudomonas aeruginosa pathogenicity island PAPI-1 is transferred via a novel type IV pilus. J Bacteriol 192: 3249–3258.MQ CarterJ. ChenS. Lory2010The Pseudomonas aeruginosa pathogenicity island PAPI-1 is transferred via a novel type IV pilus.J Bacteriol19232493258

48. He J, Baldini RL, Deziel E, Saucier M, Zhang Q, et al. (2004) The broad host range pathogen Pseudomonas aeruginosa strain PA14 carries two pathogenicity islands harboring plant and animal virulence genes. Proc Natl Acad Sci U S A 101: 2530–2535.J. HeRL BaldiniE. DezielM. SaucierQ. Zhang2004The broad host range pathogen Pseudomonas aeruginosa strain PA14 carries two pathogenicity islands harboring plant and animal virulence genes.Proc Natl Acad Sci U S A10125302535

49. Levine E, Zhang Z, Kuhlman T, Hwa T (2007) Quantitative characteristics of gene regulation by small RNA. PLoS Biol 5: e229.E. LevineZ. ZhangT. KuhlmanT. Hwa2007Quantitative characteristics of gene regulation by small RNA.PLoS Biol5e229