Genome-wide detection of predicted non-coding RNAs in Rhizobium etli expressed during free-living and host-associated growth using a high-resolution tiling array
© Vercruysse et al; licensee BioMed Central Ltd. 2010
Received: 29 July 2009
Accepted: 20 January 2010
Published: 20 January 2010
Non-coding RNAs (ncRNAs) play a crucial role in the intricate regulation of bacterial gene expression, allowing bacteria to quickly adapt to changing environments. In the past few years, a growing number of regulatory RNA elements have been predicted by computational methods, mostly in well-studied γ-proteobacteria but lately in several α-proteobacteria as well. Here, we have compared an extensive compilation of these non-coding RNA predictions to intergenic expression data of a whole-genome high-resolution tiling array in the soil-dwelling α-proteobacterium Rhizobium etli.
Expression of 89 candidate ncRNAs was detected, both on the chromosome and on the six megaplasmids encompassing the R. etli genome. Of these, 11 correspond to functionally well characterized ncRNAs, 12 were previously identified in other α-proteobacteria but are as yet uncharacterized and 66 were computationally predicted earlier but had not been experimentally identified and were therefore classified as novel ncRNAs. The latter comprise 17 putative sRNAs and 49 putative cis-regulatory ncRNAs. A selection of these candidate ncRNAs was validated by RT-qPCR, Northern blotting and 5' RACE, confirming the existence of 4 ncRNAs. Interestingly, individual transcript levels of numerous ncRNAs varied during free-living growth and during interaction with the eukaryotic host plant, pointing to possible ncRNA-dependent regulation of these specialized processes.
Our data support the practical value of previous ncRNA prediction algorithms and significantly expand the list of candidate ncRNAs encoded in the intergenic regions of R. etli and, by extension, of α-proteobacteria. Moreover, we show high-resolution tiling arrays to be suitable tools for studying intergenic ncRNA transcription profiles across the genome. The differential expression levels of some of these ncRNAs may indicate a role in adaptation to changing environmental conditions.
The first bacterial non-coding RNAs (ncRNAs) were discovered over 25 years ago [1, 2]. Still, only in the past decade have we begun appreciating their crucial role in bacterial gene regulation in response to environmental changes. By controlling metabolic pathways or stress responses, these ncRNAs play a role in diverse biological processes, including regulation of outer membrane proteins or transporters, iron metabolism, pathogenesis, quorum sensing and plasmid copy number [3–7].
The regulatory mode of action of ncRNAs is diverse as well. The best characterized group of RNA regulators are short transcripts termed small RNAs (sRNAs) that regulate gene expression through base pairing with mRNA and are either cis- or trans-encoded . Other ncRNAs can bind to proteins in order to modulate protein activity . Regulatory RNAs also include mRNA leader sequences that control expression of the downstream genes. These cis-regulatory RNA elements can be antisense RNA controlling mRNA transcription or 5' untranslated regions (UTRs) modulating expression through conformation changes by temperature shift or binding of specific metabolites [10, 11]. This kind of regulation can lead to premature transcription termination of the 5' UTR, concomitantly producing a short transcript. Recently, a new group of RNA regulators was discovered, called CRISPR (clustered regulatory interspaced short palindromic repeats) RNAs. These RNAs provide resistance to bacteriophage infection and prevent plasmid conjugation .
Although just a small number of ncRNAs were known in E. coli initially, the use of computational predictions changed this dramatically [13–16]. Today, over 80 sRNAs are known in E. coli. In recent years the search was extended to many more bacterial species such as Bacillus subtilis, Vibrio cholerae, Pseudomonas aeruginosa, Staphylococcus aureus, Streptomyces coelicolor, Salmonella enterica, Mycobacterium tuberculosis and Listeria monocytogenes [17–23]. The experimental approaches for identification include computational predictions, direct detection by dedicated microarrays or Northern blotting, direct isolation (RNomics), co-purification with RNA-binding proteins and high-throughput pyrosequencing [24, 25]. In addition, advances in array technology and the growing list of sequenced microbial replicons make custom-design high-density arrays increasingly affordable and attractive for a multitude of organisms, and expression-based ncRNA discovery and transcription profiling on a genome-wide scale feasible. Still, this was not put into practice until very recently [26, 27].
In this study, we used a high-resolution tiling array representing the entire genome of Rhizobium etli, the nitrogen-fixing endosymbiont of the common bean plant Phaseolus vulgaris [28, 29], to perform a focused study of transcriptionally active intergenic regions (IGR). Loci showing significant expression were compared to an extensive compilation of recently published ncRNA predictions in R. etli and related α-proteobacteria. 89 candidate ncRNAs similar to one or more predicted or previously detected ncRNAs were detected, and a selection of these was confirmed by Northern analysis and 5' RACE. Numerous ncRNAs are differentially expressed in R. etli during free-living growth and symbiosis with the eukaryotic host. Our results therefore significantly expand the known repertoire of ncRNAs in α-proteobacteria and provide a wealth of information for future studies to build on.
Bacterial strains and growth conditions
RNA isolation and detection by tiling microarray
Total RNA was isolated without sRNA enrichment by adapting previously published protocols [33–35]. The RNA content of 20 ml and 40 ml bacterial culture in respectively exponential and stationary phase was stabilized by adding 1/5 volume of ice-cold phenol:ethanol (5:95). Cells were harvested by centrifugation, pellets were immediately frozen in liquid nitrogen and stored at -80°C. Liquid nitrogen was used to flash freeze nodules and to crush them with a mortar and pestle. Each sample of 10 g powder, obtained from 25 plants, was suspended in 1 volume of RNAprotect (Qiagen) and 2 volumes of 10 mM MgSO4 to further stabilize the bacteroid RNA. Bacteroid and plant material was separated by differential centrifugation. The bacteroid pellet was immediately frozen in liquid nitrogen and stored at -80°C.
Total RNA was extracted using the TRIzol Plus RNA Purification kit (Invitrogen). The cell or bacteroid pellets were resuspended in 1 ml of TRIzol and shaken twice by a Precellys 24 (Bertin Technologies) at 6500 rpm for 45 seconds with 0.25 ml of 0.1 mm glass beads before following the manufacturer's instructions. Phase Lock Gel tubes (heavy type) were used to efficiently separate the organic and aqueous phase. DNA contamination was removed by two treatments of 2 μl TURBO DNase (Ambion) and afterwards checked by PCR (45 cycles). To increase RNA yields and account for experimental variation, RNA from 6 different cultures or 4 batches of bacteroids was pooled. RNA was precipitated in 3 volumes of isopropanol and 1/10 volume of sodium acetate, washed twice in ethanol and dissolved in nuclease-free ultrapure water. RNA integrity was analyzed using Experion RNA StdSens Chips (Biorad, Hercules, CA, before and after precipitation) . All samples had an RNA Quality Indicator value of 10. The ncRNA peak could be detected in each sample (Figure 1B). RNA quantity and purity was assessed using the NanoDrop ND-1000. The A260/A280 ratio and A260/A230 ratio of all samples were ≥ 2.
cDNA was synthesized using random decamers (Ambion) and the SuperScript Double-Stranded cDNA Synthesis Kit (Invitrogen) according to the manufacturer's protocol.
A whole-genome tiling array covering the entire R. etli genome sequence (6.5 Mbp in total) was designed by NimbleGen Systems, Inc. (Madison, WI), with ~385,000 60mer probes having an average start-to-start spacing of 13 base pairs and consequently an average overlap of 47 base pairs. Samples were hybridized and scanned by NimbleGen Systems [37, 38].
Microarray data preprocessing
Data preprocessing was done by performing a nonlinear intensity-dependent rescaling on non-background corrected data. To this end, a loess fit normalization  with a span of 25% was performed for each array compared against an artificial reference array, consisting of the median intensity values of each probe across all arrays. To ensure that the artificial reference itself was not altered by this rescaling, the artificial reference expression levels were chosen for the average log intensity in the loess fit (instead of the mean expression levels of the respective array and the artificial reference). A robust median-polish procedure was used to combine measurements from multiple probes into a single value .
The data were deposited in the NCBI Gene Expression Omnibus (GEO) and can be accessed through accession numbers GSE18580, GSM462173, GSM462178, GSM462180, GSM462184 and GSM462187.
A list of 1814 computationally predicted ncRNAs in R. etli and other α-proteobacteria was compiled from literature [26, 41–45]. This list was used as a query in a BLASTN search against all R. etli intergenic regions (IGR) ≥60 bp with an E-value threshold of 10-5 and with otherwise default parameters, resulting in 447 non-redundant candidate ncRNA regions. The IGRs were extracted from the NCBI Genome database and defined as the regions separating annotated genes. To determine the significance of the obtained expression in each condition, a robust estimation of the noise in the expression data was carried out for each experiment (red curve, see Additional file 1, Figure S1). The cut off (red vertical line, see Additional file 1, Figure S1) was based on the normal inverse cumulative distribution function at 99.9%. 89 regions corresponding to putative ncRNAs were found to be expressed above this significance threshold (see Additional file 2, Table S1). Differential expression across the five conditions was determined by applying the same procedure to the distribution of the standard deviations of the normalized intensities over all conditions (see Additional file 3, Figure S2). Every significantly expressed candidate ncRNA was used to query the Rfam database and RibEx (Riboswitch Explorer) web server . RibEx compares query sequences to known riboswitches as well as other predicted, but highly conserved, bacterial regulatory elements. Detailed information on the analysis of the candidate ncRNAs is summarized in Additional file 2, Table S1.
Expression levels were verified by reverse transcription quantitative real-time PCR (RT-qPCR) using the StepOnePlus System and Power SYBR Green PCR Master Mix containing AmpliTaq Gold DNA Polymerase (Applied Biosystems). Primers were designed using Primer Express 3.0 (optimal primer length of 20 bases, GC-content of 40-60% and Tm of 58-60°C, see Additional file 4, Table S2) and purchased from MacroGen. Secondary structures and dimer formation were checked with Oligoanalyzer 3.1. 2 μg of pooled total RNA of each growth condition (early/late exponential phase, stationary phase, symbiotic state 2 and 3 weeks after inoculation) was reverse transcribed to single stranded cDNA using the SuperScript VILO cDNA Synthesis Kit, including SuperScript III Reverse Transcriptase and random primers, according to the manufacturer's instructions (Invitrogen). DNA contamination of the RNA samples was checked by PCR (45 cycles) before RT and a negative control without cDNA template was included during qPCR. cDNA was stored at -80°C and stock solutions were prepared to minimize freeze-thaw cycles. After dilution of cDNA, 2 μl of cDNA (20 ng/μl), 2 μl of each specific primer (200 nM) and 4 μl of nuclease-free water were mixed with 10 μl of Power SYBR Green PCR Master Mix. In order to confirm that there was no background contamination, a negative control was included in each run. PCR conditions were: a holding stage of 10 min at 95°C, a cycling stage of 45 cycles of 15 s at 95°C and 1 min at 60°C and a melting curve stage of 15 s at 95°C, 1 min at 60°C increased to 95°C with steps of 0.3°C. The last stage was used to verify the specificity of each PCR reaction. All reactions were performed in triplicate and carried out in fast optical 96-well reaction plates (MicroAmp using optical adhesive film (MicroAmp)) with heat bonding as the method of sealing. The absence of inhibitors and the efficiency of each primer were determined by standard curves with dilution series of cDNA (5 log10 concentrations) in each run for the reference gene and the ncRNA of interest. The calibration curve's linear interval included the interval for each ncRNA of interest being quantified. The raw data was analyzed using StepOne Software v2.1. 16S RNA (RHE_CH00059) showed an invariant expression under the experimental conditions and was used as reference gene. The early exponential phase was used as calibrator condition. Relative gene expression was calculated using the Pfaffl method that corrects for differences in amplification efficiency .
Northern hybridization was performed using 1 μg of single strand DNA probe 5' end-labeled with digoxigenin (see Additional file 4, Table S2). RiboRuler RNA Ladder, Low Range (Fermentas) was used to estimate the sizes of the RNA bands. Total RNA (10-15 μg) was separated on 6% PAGE-urea gels and transferred to Hybond-N nylon membranes (Amersham) by electroblotting. The membranes were hybridized overnight in ULTRAhyb-Oligo Buffer (Ambion) at 42°C. After hybridization, membranes were first washed with buffer 1 (100 mM maleic acid, 150 mM NaCl, 7 g NaOH, pH 7.5) followed by a wash with buffer 2 (5 ml blocking stock (10% blocking reagent in buffer 1) and 45 ml buffer 1). 4 μl of anti-dioxigenin AP-Fab fragments (Roche) were added to 20 ml of buffer 2 and incubated at room temperature for 30 min. Unbound antibodies were removed by two washes with buffer 1. The membrane was equilibrated for 2 min with buffer 3 (100 mM Tris-HCl, 100 mM NaCl, 50 mM MgCl2, pH 9.5). Finally, 10 μl of chloro-5-substituted adamantyl-1,2-dioxetane phosphate was added to the membrane and incubated for 10 min at 37°C. All membranes were exposed for 10 min to an X-ray film.
Rapid amplification of 5' complementary DNA ends (5' RACE) was performed using the FirstChoice RLM-RACE kit (Ambion) according to the manufacturer's instructions, except that the CIP treatment at the start was omitted because prokaryotic RNA was used. A control without tobacco acid pyrophosphatase (TAP) treatment was included each time . Sequences of gene specific inner and outer primers are listed in Additional file 4, Table S2. 40 PCR cycles were performed at 58°C using 2 μl reverse transcription reaction (SuperScript VILO cDNA Synthesis Kit, Invitrogen) and Taq polymerase (Westburg). 5' RACE products were analyzed using 2% agarose gels and the specific TAP treated products were cloned into pCRII-TOPO (Invitrogen). Between 5 and 10 clones carrying inserts of the expected size were sequenced.
Results and Discussion
Overview of the predicted and detected ncRNAs.
del Val et al. 2007
Weinberg et al. 2007
Livny et al. 2008
Ulve et al. 2008
Valverde et al. 2008
Landt et al. 2008
Classification of the detected candidate ncRNAs(a)
Functionally characterized ncRNAs
Expression of 4 sRNAs that are highly conserved among bacteria, including 6S RNA, the signal recognition particle RNA 4.5S (SRP), bacterial RNase P class A and tmRNA, was observed. Furthermore, several known riboswitches and replication incompatibility factors were also detected, thus providing a first validation of our approach to identify ncRNAs.
6S RNA is known to associate with RNA polymerase holoenzymes containing σ70 [48, 49], blocking σ70-dependent transcription during stationary phase when 6S is abundant. Transcription from many σ70-dependent promoters will be inhibited while transcription from σS-promoters will increase. This is one mechanism that allows the stationary phase sigma factor σS (rpoS) to be an effective regulator in E. coli . However, ε- and α-proteobacteria, including R. etli, do not have an rpoS homologue. Although an alternative mechanism functionally equivalent to rpoS is not known, 6S RNA is highly expressed in R. etli during the stationary phase and might therefore play a similar role as it does in E. coli.
SRP is a universally conserved ribonucleoprotein implicated in the translation and targeting of proteins to cell membranes. The SRP of most bacteria is composed of the Ffh protein and the 4.5S RNA molecule [50, 51]. The expression of 4.5S RNA in R. etli is highest during stationary phase.
RNase P is an omnipresent endoribonuclease, found in bacteria, Archaea and Eukarya including mitochondria and chloroplasts. The processing of precursor-tRNAs into tRNAs with mature 5'-ends is its best characterized function . RNase P of R. etli is highly expressed in both free-living and symbiotic conditions, but especially in the stationary phase.
tmRNA or SsrA directly affects gene expression in general. It rescues stalled ribosomes and tags incomplete polypeptides for degradation . Even though the sensitivity to tmRNA defects varies with species and growth conditions, tmRNA seems to play a role in the ability of cells to adapt to and survive in diverse environments . R. etli tmRNA is expressed in all conditions, but primarily in the stationary phase.
12 ncRNAs are homologous to one or more previously identified but functionally uncharacterized ncRNAs in other α-proteobacteria. Remarkably, there is a sizeable overlap for this class of ncRNA between the results of the various predictive studies used. For example, ReC06 is similar to a prediction of five out of the six studies used, while ReC25, ReC58 and ReC59 were predicted by four out of six. This is partly due to three studies focusing on Sinorhizobium meliloti yielding redundant predictions using different computational approaches [41, 42, 45]. ReC11 (Figure 2B), for example, is similar to a S. meliloti sRNA highly expressed during symbiosis. This is also the case in R. etli, based on array data and as verified by RT-qPCR.
Three candidate ncRNAs share similarity with elements identified previously in Agrobacterium tumefaciens and other α-proteobacteria . ReC58 and ReC59 are both similar to the regulatory RNA element suhB, an sRNA that is probably involved in antisense gene regulation. Neither ReC58 nor ReC59 is flanked by suhB homologues in R. etli. Also in accordance with previous detection in A. tumefaciens, ReC69 is located upstream of the R. etli serCA-operon.
Our tiling array expression data provide experimental evidence for 66 ncRNAs that are homologous to previously predicted ncRNAs for which no experimental evidence was reported prior to this study. These candidate ncRNAs comprise 17 putative sRNAs and 49 putative cis-regulatory ncRNA elements. The majority of these novel ncRNAs, 48 in total, correspond to predictions of Livny et al. .
ReC43 and ReC49 are possible transcriptionally independent trans-regulatory sRNAs (Figure 2C, D). Both are highly expressed under specific conditions, respectively stationary phase and symbiosis. The expression levels of both sRNAs were confirmed by RT-qPCR. 7 of the 49 cis-regulatory ncRNA elements, including ReC27 (Figure 2E), have one or more regulatory elements that are known to depend on structured RNA, called riboswitch-like elements . No expression was detected downstream of these putative RS, ReC33 and ReC41 excepted. This was also the case for the known RS (see above), suggesting that the RS-mediated regulation occurs mainly via transcription termination. This is somewhat unexpected as most described RS in Gram-negative bacteria function by inhibiting translation initiation .
The genes downstream of the remaining 42 novel cis-regulatory ncRNAs often showed a lower expression or no expression at all, indicating the presence of possible 5' UTR fragments giving rise to short transcripts. These sRNAs may have independent functions or, alternatively, be byproducts of (post-) transcriptional regulation. Similar findings were reported previously in E. coli [59–61] and L. monocytogenes . ReC62 is an example of a ncRNA that overlaps the 3' UTR (Figure 2F).
A well-studied regulator in rhizobia is the RpoN sigma factor (σ54) that is required for nitrogen assimilation and nitrogen fixation during symbiosis in particular. R. etli contains two RpoN paralogs . One is needed during free-living growth, the other during symbiosis. No significant matches to the RpoN binding site consensus sequence could be detected upstream of the detected candidate ncRNAs. Similarly, Livny et al. (2008) reported the absence of LexA, σ54 and Fur binding sites upstream of predicted R. etli ncRNAs .
Condition-specific expression and validation of array data
Differential expression of several ncRNAs (TPPb RS, ReC11, ReC27, ReC43, ReC49, ReC62) was independently confirmed using RT-qPCR, the data obtained using this complementary technique correlating well with the expression levels estimated from the array data (Figure 2).
In addition to the primary transcripts, a smaller band was detected for ReC12, ReC14 and ReC56 (Figure 4). These bands indicate endonucleolytic processing of the primary transcripts into smaller fragments. Processing is a common feature of ncRNAs. Stable sRNAs like 6S, 4.5S and tmRNA are observed to mature by 5'- and 3'-endonucleolytic cleavage, as is the case for other sRNAs such as RprA, SraC and SraG in E. coli . For ReC64, an additional transcript, larger than expected was observed. It is unclear where this fragment might have originated from and why it was not detected on the array. Evidence for processing was also observed during 5' RACE experiments, as a less specific signal was detected in the RNA samples that were not treated with TAP. Treatment with and without TAP allows the distinction between 5'-triphosphate ends (transcriptional initiation sites) and 5'-monophosphates generated by processing of primary transcripts. Therefore, the processing likely occurs upstream of the inner RACE primer. The processed RACE products were not analyzed in this study. A more detailed analysis in the future could clarify the precise processing sites.
Genomic distribution and conservation in other α-proteobacteria
In recent years, a number of studies were published predicting ncRNAs in many α-proteobacteria including R. etli. We therefore decided to put these predictions to the test by combining R. etli tiling array expression data and a comprehensive analysis of a large number of predicted ncRNAs. This allowed us to detect 89 ncRNAs out of 447 candidate ncRNA regions. Undoubtedly, there are still many more ncRNAs to be discovered. Improvements in computational analyses and the inclusion of more experimental conditions will surely contribute to this number, as will de novo ncRNA discovery starting from the expression data.
We were able to discern well-characterized ncRNAs like 6S RNA, tmRNA and a TPP RS, several previously reported but uncharacterized ncRNAs, as well as a large number of novel ncRNAs similar to earlier predictions that had not been detected experimentally before. Additional experimental evidence for the detected ncRNAs was obtained by Northern analysis for 4 novel ncRNAs and two of them were also verified by 5' RACE. Our results show an endonucleolytic processing of the selected ncRNAs. Conservation analysis showed that a significant number of ncRNAs is conserved beyond closely related species. With condition-specific expression patterns providing a first clue to the role that some of these ncRNAs may play, a further functional analysis will help to better understand the intricate details of ncRNA-mediated gene regulation allowing bacteria to adapt to different and alternating environmental conditions.
MV is indebted to the Institute for the Promotion of Innovation through Science and Technology in Flanders (IWT-Vlaanderen). This work was supported by grants from the Research Council of the K.U. Leuven (GOA/2003/09 and GOA/011/2008) and from the Fund for Scientific Research-Flanders (G.0287.04 and G.0637.06). We thank David De Coster and Sigrid De Keersmaecker (CMPG-S&P) for their generous advice and help in isolating RNA and Northern blotting. Ann Jans, Cyrielle Kint and Serge Beullens (CMPG-SPI) are acknowledged for their assistance in 5' RACE. Authors' website: http://www.biw.kuleuven.be/dtp/cmpg/spi/index.aspx.
- Brownlee GG: Sequence of 6S RNA of E. coli. Nat New Biol. 1971, 229: 147-149.PubMedView ArticleGoogle Scholar
- Stougaard P, Molin S, Nordstrom K: RNAs involved in copy-number control and incompatibility of plasmid R1. Proc Natl Acad Sci USA. 1981, 78: 6008-6012.PubMed CentralPubMedView ArticleGoogle Scholar
- Gottesman S: Micros for microbes: non-coding regulatory RNAs in bacteria. Trends Genet. 2005, 21: 399-404.PubMedView ArticleGoogle Scholar
- Storz G, Altuvia S, Wassarman KM: An abundance of RNA regulators. Annu Rev Biochem. 2005, 74: 199-217.PubMedView ArticleGoogle Scholar
- Vogel J, Sharma CM: How to find small non-coding RNAs in bacteria. Biol Chem. 2005, 386: 1219-1238.PubMedGoogle Scholar
- Wassarman KM: Small RNAs in bacteria: diverse regulators of gene expression in response to environmental changes. Cell. 2002, 109: 141-144.PubMedView ArticleGoogle Scholar
- Waters LS, Storz G: Regulatory RNAs in bacteria. Cell. 2009, 136: 615-628.PubMed CentralPubMedView ArticleGoogle Scholar
- Brantl S: Regulatory mechanisms employed by cis-encoded antisense RNAs. Curr Opin Microbiol. 2007, 10: 102-109.PubMedView ArticleGoogle Scholar
- Babitzke P, Romeo T: CsrB sRNA family: sequestration of RNA-binding regulatory proteins. Curr Opin Microbiol. 2007, 10: 156-163.PubMedView ArticleGoogle Scholar
- Winkler WC: Riboswitches and the role of noncoding RNAs in bacterial metabolic control. Curr Opin Chem Biol. 2005, 9: 594-602.PubMedView ArticleGoogle Scholar
- Johansson J: RNA thermosensors in bacterial pathogens. Contrib Microbiol. 2009, 16: 150-160.PubMedView ArticleGoogle Scholar
- Sorek R, Kunin V, Hugenholtz P: CRISPR--a widespread system that provides acquired resistance against phages in bacteria and archaea. Nat Rev Microbiol. 2008, 6: 181-186.PubMedView ArticleGoogle Scholar
- Argaman L, Hershberg R, Vogel J, Bejerano G, Wagner EG, Margalit H, Altuvia S: Novel small RNA-encoding genes in the intergenic regions of Escherichia coli. Curr Biol. 2001, 11: 941-950.PubMedView ArticleGoogle Scholar
- Chen S, Lesnik EA, Hall TA, Sampath R, Griffey RH, Ecker DJ, Blyn LB: A bioinformatics based approach to discover small RNA genes in the Escherichia coli genome. Biosystems. 2002, 65: 157-177.PubMedView ArticleGoogle Scholar
- Rivas E, Klein RJ, Jones TA, Eddy SR: Computational identification of noncoding RNAs in E. coli by comparative genomics. Curr Biol. 2001, 11: 1369-1373.PubMedView ArticleGoogle Scholar
- Wassarman KM, Repoila F, Rosenow C, Storz G, Gottesman S: Identification of novel small RNAs using comparative genomics and microarrays. Genes Dev. 2001, 15: 1637-1651.PubMed CentralPubMedView ArticleGoogle Scholar
- Saito S, Kakeshita H, Nakamura K: Novel small RNA-encoding genes in the intergenic regions of Bacillus subtilis. Gene. 2009, 428: 2-8.PubMedView ArticleGoogle Scholar
- Livny J, Brencic A, Lory S, Waldor MK: Identification of 17 Pseudomonas aeruginosa sRNAs and prediction of sRNA-encoding genes in 10 diverse pathogens using the bioinformatic tool sRNAPredict2. Nucleic Acids Res. 2006, 34: 3484-3493.PubMed CentralPubMedView ArticleGoogle Scholar
- Pichon C, Felden B: Small RNA genes expressed from Staphylococcus aureus genomic and pathogenicity islands with specific expression among pathogenic strains. Proc Natl Acad Sci USA. 2005, 102: 14249-14254.PubMed CentralPubMedView ArticleGoogle Scholar
- Swiercz JP, Hindra , Bobek J, Haiser HJ, Di Berardo C, Tjaden B, Elliot MA: Small non-coding RNAs in Streptomyces coelicolor. Nucleic Acids Res. 2008, 36: 7240-7251.PubMed CentralPubMedView ArticleGoogle Scholar
- Sittka A, Lucchini S, Papenfort K, Sharma CM, Rolle K, Binnewies TT, Hinton JC, Vogel J: Deep sequencing analysis of small noncoding RNA and mRNA targets of the global post-transcriptional regulator, Hfq. PLoS Genet. 2008, 4: e1000163-PubMed CentralPubMedView ArticleGoogle Scholar
- Mandin P, Repoila F, Vergassola M, Geissmann T, Cossart P: Identification of new noncoding RNAs in Listeria monocytogenes and prediction of mRNA targets. Nucleic Acids Res. 2007, 35: 962-974.PubMed CentralPubMedView ArticleGoogle Scholar
- Arnvig KB, Young DB: Identification of small RNAs in Mycobacterium tuberculosis. Mol Microbiol. 2009, 73: 397-408.PubMed CentralPubMedView ArticleGoogle Scholar
- Livny J, Waldor MK: Identification of small RNAs in diverse bacterial species. Curr Opin Microbiol. 2007, 10: 96-101.PubMedView ArticleGoogle Scholar
- Altuvia S: Identification of bacterial small non-coding RNAs: experimental approaches. Curr Opin Microbiol. 2007, 10: 257-261.PubMedView ArticleGoogle Scholar
- Landt SG, Abeliuk E, McGrath PT, Lesley JA, McAdams HH, Shapiro L: Small non-coding RNAs in Caulobacter crescentus. Mol Microbiol. 2008, 68: 600-614.PubMedView ArticleGoogle Scholar
- Toledo-Arana A, Dussurget O, Nikitas G, Sesto N, Guet-Revillet H, Balestrino D, Loh E, Gripenland J, Tiensuu T, Vaitkevicius K, Barthelemy M, Vergassola M, Nahori MA, Soubigou G, Regnault B, Coppee JY, Lecuit M, Johansson J, Cossart P: The Listeria transcriptional landscape from saprophytism to virulence. Nature. 2009, 459: 950-956.PubMedView ArticleGoogle Scholar
- Braeken K, Fauvart M, Vercruysse M, Beullens S, Lambrichts I, Michiels J: Pleiotropic effects of a rel mutation on stress survival of Rhizobium etli CNPAF512. BMC Microbiol. 2008, 8: 219-PubMed CentralPubMedView ArticleGoogle Scholar
- Fauvart M, Michiels J: Rhizobial secreted proteins as determinants of host specificity in the rhizobium-legume symbiosis. FEMS Microbiol Lett. 2008, 285: 1-9.PubMedView ArticleGoogle Scholar
- Moris M, Braeken K, Schoeters E, Verreth C, Beullens S, Vanderleyden J, Michiels J: Effective symbiosis between Rhizobium etli and Phaseolus vulgaris requires the alarmone ppGpp. J Bacteriol. 2005, 187: 5460-5469.PubMed CentralPubMedView ArticleGoogle Scholar
- D'Hooghe I, Michiels J, Vlassak K, Verreth C, Waelkens F, Vanderleyden J: Structural and functional analysis of the fixLJ genes of Rhizobium leguminosarum biovar phaseoli CNPAF512. Mol Gen Genet. 1995, 249: 117-126.PubMedView ArticleGoogle Scholar
- D'Hooghe I, Wauven Vander C, Michiels J, Tricot C, de Wilde P, Vanderleyden J, Stalon V: The arginine deiminase pathway in Rhizobium etli: DNA sequence analysis and functional study of the arcABC genes. J Bacteriol. 1997, 179: 7403-7409.PubMed CentralPubMedGoogle Scholar
- Bhagwat AA, Phadke RP, Wheeler D, Kalantre S, Gudipati M, Bhagwat M: Computational methods and evaluation of RNA stabilization reagents for genome-wide expression studies. J Microbiol Methods. 2003, 55: 399-409.PubMedView ArticleGoogle Scholar
- Becker A, Berges H, Krol E, Bruand C, Ruberg S, Capela D, Lauber E, Meilhoc E, Ampe F, de Bruijn FJ, Fourment J, Francez-Charlot A, Kahn D, Kuster H, Liebe C, Puhler A, Weidner S, Batut J: Global changes in gene expression in Sinorhizobium meliloti 1021 under microoxic and symbiotic conditions. Mol Plant Microbe Interact. 2004, 17: 292-303.PubMedView ArticleGoogle Scholar
- Capela D, Filipe C, Bobik C, Batut J, Bruand C: Sinorhizobium meliloti differentiation during symbiosis with Alfalfa: a transcriptomic dissection. Mol Plant Microbe Interact. 2006, 19: 363-372.PubMedView ArticleGoogle Scholar
- Spaniolas S, Bazakos C, Awad M, Kalaitzis P: Exploitation of the chloroplast trnL (UAA) intron polymorphisms for the authentication of plant oils by means of a lab-on-a-chip capillary electrophoresis system. J Agric Food Chem. 2008, 56: 6886-6891.PubMedView ArticleGoogle Scholar
- Nuwaysir EF, Huang W, Albert TJ, Singh J, Nuwaysir K, Pitas A, Richmond T, Gorski T, Berg JP, Ballin J, McCormick M, Norton J, Pollock T, Sumwalt T, Butcher L, Porter D, Molla M, Hall C, Blattner F, Sussman MR, Wallace RL, Cerrina F, Green RD: Gene expression analysis using oligonucleotide arrays produced by maskless photolithography. Genome Res. 2002, 12: 1749-1755.PubMed CentralPubMedView ArticleGoogle Scholar
- Nelson CM, Herron MJ, Felsheim RF, Schloeder BR, Grindle SM, Chavez AO, Kurtti TJ, Munderloh UG: Whole genome transcription profiling of Anaplasma phagocytophilum in human and tick host cells by tiling array analysis. BMC Genomics. 2008, 9: 364-PubMed CentralPubMedView ArticleGoogle Scholar
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002, 30: e15-PubMed CentralPubMedView ArticleGoogle Scholar
- Tukey JW: Exploratory Data Analysis. 1977, Addison-Wesley, Menlo Park, FirstGoogle Scholar
- del Val C, Rivas E, Torres-Quesada O, Toro N, Jimenez-Zurdo JI: Identification of differentially expressed small non-coding RNAs in the legume endosymbiont Sinorhizobium meliloti by comparative genomics. Mol Microbiol. 2007, 66: 1080-1091.PubMed CentralPubMedView ArticleGoogle Scholar
- Ulve VM, Sevin EW, Cheron A, Barloy-Hubler F: Identification of chromosomal alpha-proteobacterial small RNAs by comparative genome analysis and detection in Sinorhizobium meliloti strain 1021. BMC Genomics. 2007, 8: 467-PubMed CentralPubMedView ArticleGoogle Scholar
- Weinberg Z, Barrick JE, Yao Z, Roth A, Kim JN, Gore J, Wang JX, Lee ER, Block KF, Sudarsan N, Neph S, Tompa M, Ruzzo WL, Breaker RR: Identification of 22 candidate structured RNAs in bacteria using the CMfinder comparative genomics pipeline. Nucleic Acids Res. 2007, 35: 4809-4819.PubMed CentralPubMedView ArticleGoogle Scholar
- Livny J, Teonadi H, Livny M, Waldor MK: High-throughput, kingdom-wide prediction and annotation of bacterial non-coding RNAs. PLoS ONE. 2008, 3: e3197-PubMed CentralPubMedView ArticleGoogle Scholar
- Valverde C, Livny J, Schluter JP, Reinkensmeier J, Becker A, Parisi G: Prediction of Sinorhizobium meliloti sRNA genes and experimental detection in strain 2011. BMC Genomics. 2008, 9: 416-PubMed CentralPubMedView ArticleGoogle Scholar
- Abreu-Goodger C, Merino E: RibEx: a web server for locating riboswitches and other conserved bacterial regulatory elements. Nucleic Acids Res. 2005, 33: W690-692.PubMed CentralPubMedView ArticleGoogle Scholar
- Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-PubMed CentralPubMedView ArticleGoogle Scholar
- Wassarman KM, Storz G: 6S RNA regulates E. coli RNA polymerase activity. Cell. 2000, 101: 613-623.PubMedView ArticleGoogle Scholar
- Wassarman KM: 6S RNA: a regulator of transcription. Mol Microbiol. 2007, 65: 1425-1431.PubMedView ArticleGoogle Scholar
- Regalia M, Rosenblad MA, Samuelsson T: Prediction of signal recognition particle RNA genes. Nucleic Acids Res. 2002, 30: 3368-3377.PubMed CentralPubMedView ArticleGoogle Scholar
- Rosenblad MA, Gorodkin J, Knudsen B, Zwieb C, Samuelsson T: SRPDB: Signal Recognition Particle Database. Nucleic Acids Res. 2003, 31: 363-364.PubMed CentralPubMedView ArticleGoogle Scholar
- Evans D, Marquez SM, Pace NR: RNase P: interface of the RNA and protein worlds. Trends Biochem Sci. 2006, 31: 333-341.PubMedView ArticleGoogle Scholar
- Gur E, Sauer RT: Evolution of the ssrA degradation tag in Mycoplasma: specificity switch to a different protease. Proc Natl Acad Sci USA. 2008, 105: 16113-16118.PubMed CentralPubMedView ArticleGoogle Scholar
- Nudler E, Mironov AS: The riboswitch control of bacterial metabolism. Trends Biochem Sci. 2004, 29: 11-17.PubMedView ArticleGoogle Scholar
- Miranda-Rios J, Morera C, Taboada H, Davalos A, Encarnacion S, Mora J, Soberon M: Expression of thiamin biosynthetic genes (thiCOGE) and production of symbiotic terminal oxidase cbb3 in Rhizobium etli. J Bacteriol. 1997, 179: 6887-6893.PubMed CentralPubMedGoogle Scholar
- Miranda-Rios J, Navarro M, Soberon M: A conserved RNA structure (thi box) is involved in regulation of thiamin biosynthetic gene expression in bacteria. Proc Natl Acad Sci USA. 2001, 98: 9736-9741.PubMed CentralPubMedView ArticleGoogle Scholar
- MacLellan SR, Smallbone LA, Sibley CD, Finan TM: The expression of a novel antisense gene mediates incompatibility within the large repABC family of alpha-proteobacterial plasmids. Mol Microbiol. 2005, 55: 611-623.PubMedView ArticleGoogle Scholar
- Izquierdo J, Venkova-Canova T, Ramirez-Romero MA, Tellez-Sosa J, Hernandez-Lucas I, Sanjuan J, Cevallos MA: An antisense RNA plays a central role in the replication control of a repC plasmid. Plasmid. 2005, 54: 259-277.PubMedView ArticleGoogle Scholar
- Vogel J, Bartels V, Tang TH, Churakov G, Slagter-Jager JG, Huttenhofer A, Wagner EG: RNomics in Escherichia coli detects new sRNA species and indicates parallel transcriptional output in bacteria. Nucleic Acids Res. 2003, 31: 6435-6443.PubMed CentralPubMedView ArticleGoogle Scholar
- Kawano M, Reynolds AA, Miranda-Rios J, Storz G: Detection of 5'- and 3'-UTR-derived small RNAs and cis-encoded antisense RNAs in Escherichia coli. Nucleic Acids Res. 2005, 33: 1040-1050.PubMed CentralPubMedView ArticleGoogle Scholar
- Guillier M, Gottesman S: The 5' end of two redundant sRNAs is involved in the regulation of multiple targets, including their own regulator. Nucleic Acids Res. 2008, 36: 6781-6794.PubMed CentralPubMedView ArticleGoogle Scholar
- Michiels J, Van Soom T, D'Hooghe I, Dombrecht B, Benhassine T, de Wilde P, Vanderleyden J: The Rhizobium etli rpoN locus: DNA sequence analysis and phenotypical characterization of rpoN, ptsN, and ptsA mutants. J Bacteriol. 1998, 180: 1729-1740.PubMed CentralPubMedGoogle Scholar
- Kulkarni PR, Cui X, Williams JW, Stevens AM, Kulkarni RV: Prediction of CsrA-regulating small RNAs in bacteria and their experimental verification in Vibrio fischeri. Nucleic Acids Res. 2006, 34: 3361-3369.PubMed CentralPubMedView ArticleGoogle Scholar
- Wilderman PJ, Sowa NA, FitzGerald DJ, FitzGerald PC, Gottesman S, Ochsner UA, Vasil ML: Identification of tandem duplicate regulatory small RNAs in Pseudomonas aeruginosa involved in iron homeostasis. Proc Natl Acad Sci USA. 2004, 101: 9792-9797.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.