Novel small RNA-encoding genes in the intergenic regions of Bacillus subtilis

Type I toxin-antitoxin systems in bacteria: from regulation to biological functions

Toxin-antitoxin systems are ubiquitous in the prokaryotic world and widely distributed among chromosomes and mobile genetic elements. Several different toxin-antitoxin system types exist, but what they all have in common is that toxin activity is prevented by the cognate antitoxin. In type I toxin-antitoxin systems, toxin production is controlled by an RNA antitoxin and by structural features inherent to the toxin messenger RNA. Most type I toxins are small membrane proteins that display a variety of cellular effects. While originally discovered as modules that stabilize plasmids, chromosomal type I toxin-antitoxin systems may also stabilize prophages, or serve important functions upon certain stress conditions and contribute to population-wide survival strategies. Here, we will describe the intricate RNA-based regulation of type I toxin-antitoxin systems and discuss their potential biological functions.

Identification and initial characterization of Hfq-associated sRNAs in Histophilus somni strain 2336

Small RNAs (sRNA), in association with the global chaperone regulator Hfq, positively or negatively regulate gene expression in bacteria. For this study, Histophilus somni sRNAs that bind to Hfq were identified and then partially characterized. The Hfq-associated sRNAs in H. somni were isolated and identified by co-immunoprecipitation using anti-Hfq antibody, followed by sRNA sequencing. Sequence analysis of the sRNA samples identified 100 putative sRNAs, out of which 16 were present in pathogenic strain 2336, but not in non-pathogenic strain 129Pt. Bioinformatic analyses suggested that the sRNAs HS9, HS79, and HS97 could bind to many genes putatively involved in virulence/biofilm formation. Furthermore, multi-sequence alignment of the sRNA regions in the genome revealed that HS9 and HS97 could interact with sigma 54, which is a transcription factor linked to important bacterial traits, including motility, virulence, and biofilm formation. Northern blotting was used to determine the approximate size, abundance and any processing events attributed to the sRNAs. Selected sRNA candidates were confirmed to bind Hfq, as determined by electrophoretic mobility shift assays using sRNAs synthesized by in vitro transcription and recombinant Hfq. The exact transcriptional start site of the sRNA candidates was determined by RNA ligase-mediated rapid amplification of cDNA ends, followed by cloning and sequencing. This is the first investigation of H. somni sRNAs that show they may have important regulatory roles in virulence and biofilm formation.

The RNA-RNA interactome between a phage and its satellite virus reveals a small RNA that differentially regulates gene expression across both genomes

Satellite viruses are present across all domains of life, defined as subviral parasites that require infection by another virus for satellite progeny production. Phage satellites exhibit various regulatory mechanisms to manipulate phage gene expression to the benefit of the satellite, redirecting resources from the phage to the satellite, and often inhibiting phage progeny production. While small RNAs (sRNAs) are well documented as regulators of prokaryotic gene expression, they have not been shown to play a regulatory role in satellite-phage conflicts. Vibrio cholerae encodes the phage inducible chromosomal island-like element (PLE), a phage satellite, to defend itself against the lytic phage ICP1. Here, we use Hi-GRIL-seq to identify a complex RNA-RNA interactome between PLE and ICP1. Both inter- and intragenome RNA interactions were detected, headlined by the PLE sRNA, SviR. SviR is involved in regulating both PLE and ICP1 gene expression uniquely, decreasing ICP1 target translation and affecting PLE transcripts. The striking conservation of SviR across all known PLEs suggests the sRNA is deeply rooted in the PLE-ICP1 conflict and implicates sRNAs as unidentified regulators of gene expression in phage-satellite interactions.

Cis- and Trans-Encoded Small Regulatory RNAs in Bacillus subtilis

Small regulatory RNAs (sRNAs) that act by base-pairing are the most abundant posttranscriptional regulators in all three kingdoms of life. Over the past 20 years, a variety of approaches have been employed to discover chromosome-encoded sRNAs in a multitude of bacterial species. However, although largely improved bioinformatics tools are available to predict potential targets of base-pairing sRNAs, it is still challenging to confirm these targets experimentally and to elucidate the mechanisms as well as the physiological role of their sRNA-mediated regulation. Here, we provide an overview of currently known cis- and trans-encoded sRNAs from B. subtilis with known targets and defined regulatory mechanisms and on the potential role of RNA chaperones that are or might be required to facilitate sRNA regulation in this important Gram-positive model organism.

RNAxplorer: Harnessing the Power of Guiding Potentials to Sample RNA Landscapes

Motivation

Predicting the folding dynamics of RNAs is a computationally difficult problem, first and foremost due to the combinatorial explosion of alternative structures in the folding space. Abstractions are therefore needed to simplify downstream analyses, and thus make them computationally tractable. This can be achieved by various structure sampling algorithms. However, current sampling methods are still time consuming and frequently fail to represent key elements of the folding space.

Method

We introduce RNAxplorer, a novel adaptive sampling method to efficiently explore the structure space of RNAs. RNAxplorer uses dynamic programming to perform an efficient Boltzmann sampling in the presence of guiding potentials, which are accumulated into pseudo-energy terms and reflect similarity to already well-sampled structures. This way, we effectively steer sampling toward underrepresented or unexplored regions of the structure space.

Results

We developed and applied different measures to benchmark our sampling methods against its competitors. Most of the measures show that RNAxplorer produces more diverse structure samples, yields rare conformations that may be inaccessible to other sampling methods and is better at finding the most relevant kinetic traps in the landscape. Thus, it produces a more representative coarse graining of the landscape, which is well suited to subsequently compute better approximations of RNA folding kinetics.

Availabilityand implementation

https://github.com/ViennaRNA/RNAxplorer/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Movement of small RNAs in and between plants and fungi

RNA interference is a biological process whereby small RNAs inhibit gene expression through neutralizing targeted mRNA molecules. This process is conserved in eukaryotes. Here, recent work regarding the mechanisms of how small RNAs move within and between organisms is examined. Small RNAs can move locally and systemically in plants through plasmodesmata and phloem, respectively. In fungi, transportation of small RNAs may also be achieved by septal pores and vesicles. Recent evidence also supports bidirectional cross-kingdom communication of small RNAs between host plants and adapted fungal pathogens to affect the outcome of infection. We discuss several mechanisms for small RNA trafficking and describe evidence for transport through naked form, combined with RNA-binding proteins or enclosed by vesicles.

sRNA EsrE Is Transcriptionally Regulated by the Ferric Uptake Regulator Fur in Escherichia coli

Small RNAs (sRNAs) are widespread and play major roles in regulation circuits in bacteria. Previously, we have demonstrated that transcription of esrE is under the control of its own promoter. However, the regulatory elements involved in EsrE sRNA expression are still unknown. In this study, we found that different cis-regulatory elements exist in the promoter region of esrE. We then screened and analyzed seven potential corresponding trans-regulatory elements by using pull-down assays based on DNA affinity chromatography. Among these candidate regulators, we investigated the relationship between the ferric uptake regulator (Fur) and the EsrE sRNA. Electrophoresis mobility shift assays (EMSAs) and β-galactosidase activity assays demonstrated that Fur can bind to the promoter region of esrE, and positively regulate EsrE sRNA expression in the presence of Fe²⁺.

A new sRNA-mediated posttranscriptional regulation system for Bacillus subtilis

Many genetic tools for gene regulation have been developed during the past decades. Some of them edit genomic DNA, such as nucleotides deletions and insertions, while the others interfere with the gene transcriptions or messenger RNA translation. Here, we report a posttranscriptional regulation tool which is termed "Modulation via the small RNA (sRNA)-dependent operation system: MS-DOS" by engineering the type I toxin-antitoxin system in Bacillus subtilis. MS-DOS depends simply on insertion of an operation region (OPR; partial toxin-encoding region) downstream of a genomic open reading frame of interest and overexpression of the coupling antitoxin sRNA from a plasmid. Pairing between the OPR and the sRNA will trigger the RNAse degradation of the transcripts of selected genes. MS-DOS allows for the quantitative, specific, and reversible knockdown of single or multiple genomic genes in B. subtilis. We also showed that the truncated antitoxin SR4 with 53 nt length is sufficient to repress gene expression. Superior to other existing RNA based interfering systems, MS-DOS allows simultaneous knockdown of multiple genes with effortless expression of a single antitoxin RNA. This sRNA-guided repression system will further enrich the gene regulation tools and expand the gene regulation flexibility.

Dual Expression of the Salmonella Effector SrfJ in Mammalian Cells and Plants

SrfJ is an effector of the Salmonella pathogenicity island 2-encoded type III secretion system. Salmonella enterica serovar Typhimurium expresses srfJ under two disparate sets of conditions: media with low Mg2+ and low pH, imitating intravacuolar conditions, and media with myo-inositol (MI), a carbohydrate that can be used by Salmonella as sole carbon source. We investigated the molecular basis for this dual regulation. Here, we provide evidence for the existence of two distinct promoters that control the expression of srfJ. A proximal promoter, PsrfJ, responds to intravacuolar signals and is positively regulated by SsrB and PhoP and negatively regulated by RcsB. A second distant promoter, PiolE, is negatively regulated by the MI island repressor IolR. We also explored the in vivo activity of these promoters in different hosts. Interestingly, our results indicate that the proximal promoter is specifically active inside mammalian cells whereas the distant one is expressed upon Salmonella colonization of plants. Importantly, we also found that inappropriate expression of srfJ leads to reduced proliferation inside macrophages whereas lack of srfJ expression increases survival and decreases activation of defense responses in plants. These observations suggest that SrfJ is a relevant factor in the interplay between Salmonella and hosts of different kingdoms.

Iron-regulated small RNA expression as Neisseria gonorrhoeae FA 1090 transitions into stationary phase growth

BACKGROUND

For most pathogens, iron (Fe) homeostasis is crucial for maintenance within the host and the ability to cause disease. The primary transcriptional regulator that controls intracellular Fe levels is the Fur (ferric uptake regulator) protein, which exerts its action on transcription by binding to a promoter-proximal sequence termed the Fur box. Fur-regulated transcriptional responses are often fine-tuned at the post-transcriptional level through the action of small regulatory RNAs (sRNAs). Consequently, identifying sRNAs contributing to the control of Fe homeostasis is important for understanding the Fur-controlled bacterial Fe-response network.

RESULTS

In this study, we sequenced size-selected directional libraries representing sRNA samples from Neisseria gonorrhoeae strain FA 1090, and examined the Fe- and temporal regulation of these sRNAs. RNA-seq data for all time points identified a pool of at least 340 potential sRNAs. Differential analysis demonstrated that expression appeared to be regulated by Fe availability for at least fifteen of these sRNAs. Fourteen sRNAs were induced in high Fe conditions, consisting of both cis and trans sRNAs, some of which are predicted to control expression of a known virulence factor, and one SAM riboswitch. An additional putative cis-acting sRNA was repressed by Fe availability. In the pathogenic Neisseria species, one sRNA that contributes to Fe-regulated post-transcriptional control is the Fur-repressible sRNA NrrF. The expression of five Fe-induced sRNAs appeared to be at least partially controlled by NrrF, while the remainder was expressed independently of NrrF. The expression of the 14 Fe-induced sRNAs also exhibited temporal control, as their expression levels increased dramatically as the bacteria entered stationary phase.

CONCLUSIONS

Here we report the temporal expression of Fe-regulated sRNAs in N. gonorrhoeae FA 1090 with several appearing to be controlled by the Fe-repressible sRNA NrrF. Temporal regulation of these sRNAs suggests a regulatory role in controlling functions necessary for survival, and may be important for phenotypes often associated with altered growth rates, such as biofilm formation or intracellular survival. Future functional studies will be needed to understand how these regulatory sRNAs contribute to gonococcal biology and pathogenesis.

Regulatory RNAs in Bacillus subtilis: a Gram-Positive Perspective on Bacterial RNA-Mediated Regulation of Gene Expression

Bacteria can employ widely diverse RNA molecules to regulate their gene expression. Such molecules include trans-acting small regulatory RNAs, antisense RNAs, and a variety of transcriptional attenuation mechanisms in the 5' untranslated region. Thus far, most regulatory RNA research has focused on Gram-negative bacteria, such as Escherichia coli and Salmonella. Hence, there is uncertainty about whether the resulting insights can be extrapolated directly to other bacteria, such as the Gram-positive soil bacterium Bacillus subtilis. A recent study identified 1,583 putative regulatory RNAs in B. subtilis, whose expression was assessed across 104 conditions. Here, we review the current understanding of RNA-based regulation in B. subtilis, and we categorize the newly identified putative regulatory RNAs on the basis of their conservation in other bacilli and the stability of their predicted secondary structures. Our present evaluation of the publicly available data indicates that RNA-mediated gene regulation in B. subtilis mostly involves elements at the 5' ends of mRNA molecules. These can include 5' secondary structure elements and metabolite-, tRNA-, or protein-binding sites. Importantly, sense-independent segments are identified as the most conserved and structured potential regulatory RNAs in B. subtilis. Altogether, the present survey provides many leads for the identification of new regulatory RNA functions in B. subtilis.

A Transcriptome Map of Actinobacillus pleuropneumoniae at Single-Nucleotide Resolution Using Deep RNA-Seq

Actinobacillus pleuropneumoniae is the pathogen of porcine contagious pleuropneumoniae, a highly contagious respiratory disease of swine. Although the genome of A. pleuropneumoniae was sequenced several years ago, limited information is available on the genome-wide transcriptional analysis to accurately annotate the gene structures and regulatory elements. High-throughput RNA sequencing (RNA-seq) has been applied to study the transcriptional landscape of bacteria, which can efficiently and accurately identify gene expression regions and unknown transcriptional units, especially small non-coding RNAs (sRNAs), UTRs and regulatory regions. The aim of this study is to comprehensively analyze the transcriptome of A. pleuropneumoniae by RNA-seq in order to improve the existing genome annotation and promote our understanding of A. pleuropneumoniae gene structures and RNA-based regulation. In this study, we utilized RNA-seq to construct a single nucleotide resolution transcriptome map of A. pleuropneumoniae. More than 3.8 million high-quality reads (average length ~90 bp) from a cDNA library were generated and aligned to the reference genome. We identified 32 open reading frames encoding novel proteins that were mis-annotated in the previous genome annotations. The start sites for 35 genes based on the current genome annotation were corrected. Furthermore, 51 sRNAs in the A. pleuropneumoniae genome were discovered, of which 40 sRNAs were never reported in previous studies. The transcriptome map also enabled visualization of 5'- and 3'-UTR regions, in which contained 11 sRNAs. In addition, 351 operons covering 1230 genes throughout the whole genome were identified. The RNA-Seq based transcriptome map validated annotated genes and corrected annotations of open reading frames in the genome, and led to the identification of many functional elements (e.g. regions encoding novel proteins, non-coding sRNAs and operon structures). The transcriptional units described in this study provide a foundation for future studies concerning the gene functions and the transcriptional regulatory architectures of this pathogen.

dRNA-Seq Reveals Genomewide TSSs and Noncoding RNAs of Plant Beneficial Rhizobacterium Bacillus amyloliquefaciens FZB42

Bacillus amyloliquefaciens subsp. plantarum FZB42 is a representative of Gram-positive plant-growth-promoting rhizobacteria (PGPR) that inhabit plant root environments. In order to better understand the molecular mechanisms of bacteria-plant symbiosis, we have systematically analyzed the primary transcriptome of strain FZB42 grown under rhizosphere-mimicking conditions using differential RNA sequencing (dRNA-seq). Our analysis revealed 4,877 transcription start sites for protein-coding genes, identified genes differentially expressed under different growth conditions, and corrected many previously mis-annotated genes. We also identified a large number of riboswitches and cis-encoded antisense RNAs, as well as trans-encoded small noncoding RNAs that may play important roles in the gene regulation of Bacillus. Overall, our analyses provided a landscape of Bacillus primary transcriptome and improved the knowledge of rhizobacteria-host interactions.

Transcriptional analysis of Deinococcus radiodurans reveals novel small RNAs that are differentially expressed under ionizing radiation

Small noncoding RNAs (sRNAs) are posttranscriptional regulators that have been identified in multiple species and shown to play essential roles in responsive mechanisms to environmental stresses. The natural ability of specific bacteria to resist high levels of radiation has been of high interest to mechanistic studies of DNA repair and biomolecular protection. Deinococcus radiodurans is a model extremophile for radiation studies that can survive doses of ionizing radiation of >12,000 Gy, 3,000 times higher than for most vertebrates. Few studies have investigated posttranscriptional regulatory mechanisms of this organism that could be relevant in its general gene regulatory patterns. In this study, we identified 199 potential sRNA candidates in D. radiodurans by whole-transcriptome deep sequencing analysis and confirmed the expression of 41 sRNAs by Northern blotting and reverse transcriptase PCR (RT-PCR). A total of 8 confirmed sRNAs showed differential expression during recovery after acute ionizing radiation (15 kGy). We have also found and confirmed 7 sRNAs in Deinococcus geothermalis, a closely related radioresistant species. The identification of several novel sRNAs in Deinococcus bacteria raises important questions about the evolution and nature of global gene regulation in radioresistance.

Genome-wide analyses in bacteria show small-RNA enrichment for long and conserved intergenic regions

Interest in finding small RNAs (sRNAs) in bacteria has significantly increased in recent years due to their regulatory functions. Development of high-throughput methods and more sophisticated computational algorithms has allowed rapid identification of sRNA candidates in different species. However, given their various sizes (50 to 500 nucleotides [nt]) and their potential genomic locations in the 5' and 3' untranslated regions as well as in intergenic regions, identification and validation of true sRNAs have been challenging. In addition, the evolution of bacterial sRNAs across different species continues to be puzzling, given that they can exert similar functions with various sequences and structures. In this study, we analyzed the enrichment patterns of sRNAs in 13 well-annotated bacterial species using existing transcriptome and experimental data. All intergenic regions were analyzed by WU-BLAST to examine conservation levels relative to species within or outside their genus. In total, more than 900 validated bacterial sRNAs and 23,000 intergenic regions were analyzed. The results indicate that sRNAs are enriched in intergenic regions, which are longer and more conserved than the average intergenic regions in the corresponding bacterial genome. We also found that sRNA-coding regions have different conservation levels relative to their flanking regions. This work provides a way to analyze how noncoding RNAs are distributed in bacterial genomes and also shows conserved features of intergenic regions that encode sRNAs. These results also provide insight into the functions of regions surrounding sRNAs and into optimization of RNA search algorithms.

sRNA antitoxins: more than one way to repress a toxin

Bacterial toxin-antitoxin loci consist of two genes: one encodes a potentially toxic protein, and the second, an antitoxin to repress its function or expression. The antitoxin can either be an RNA or a protein. For type I and type III loci, the antitoxins are RNAs; however, they have very different modes of action. Type I antitoxins repress toxin protein expression through interacting with the toxin mRNA, thereby targeting the mRNA for degradation or preventing its translation or both; type III antitoxins directly bind to the toxin protein, sequestering it. Along with these two very different modes of action for the antitoxin, there are differences in the functions of the toxin proteins and the mobility of these loci between species. Within this review, we discuss the major differences as to how the RNAs repress toxin activity, the potential consequences for utilizing different regulatory strategies, as well as the confirmed and potential biological roles for these loci across bacterial species.

Small regulatory RNAs from low-GC Gram-positive bacteria

Small regulatory RNAs (sRNAs) that act by base-pairing were first discovered in so-called accessory DNA elements—plasmids, phages, and transposons—where they control replication, maintenance, and transposition. Since 2001, a huge body of work has been performed to predict and identify sRNAs in a multitude of bacterial genomes. The majority of chromosome-encoded sRNAs have been investigated in E. coli and other Gram-negative bacteria. However, during the past five years an increasing number of sRNAs were found in Gram-positive bacteria. Here, we outline our current knowledge on chromosome-encoded sRNAs from low-GC Gram-positive species that act by base-pairing, i.e., an antisense mechanism. We will focus on sRNAs with known targets and defined regulatory mechanisms with special emphasis on Bacillus subtilis.

RNASurface: fast and accurate detection of locally optimal potentially structured RNA segments

MOTIVATION

During the past decade, new classes of non-coding RNAs (ncRNAs) and their unexpected functions were discovered. Stable secondary structure is the key feature of many non-coding RNAs. Taking into account huge amounts of genomic data, development of computational methods to survey genomes for structured RNAs remains an actual problem, especially when homologous sequences are not available for comparative analysis. Existing programs scan genomes with a fixed window by efficiently constructing a matrix of RNA minimum free energies. A wide range of lengths of structured RNAs necessitates the use of many different window lengths that substantially increases the output size and computational efforts.

RESULTS

In this article, we present an algorithm RNASurface to efficiently scan genomes by constructing a matrix of significance of RNA secondary structures and to identify all locally optimal structured RNA segments up to a predefined size. RNASurface significantly improves precision of identification of known ncRNA in Bacillus subtilis.

AVAILABILITY AND IMPLEMENTATION

RNASurface C source code is available from http://bioinf.fbb.msu.ru/RNASurface/downloads.html.

RNA-Seq of Bacillus licheniformis: active regulatory RNA features expressed within a productive fermentation

BACKGROUND

The production of enzymes by an industrial strain requires a complex adaption of the bacterial metabolism to the conditions within the fermenter. Regulatory events within the process result in a dynamic change of the transcriptional activity of the genome. This complex network of genes is orchestrated by proteins as well as regulatory RNA elements. Here we present an RNA-Seq based study considering selected phases of an industry-oriented fermentation of Bacillus licheniformis.

RESULTS

A detailed analysis of 20 strand-specific RNA-Seq datasets revealed a multitude of transcriptionally active genomic regions. 3314 RNA features encoded by such active loci have been identified and sorted into ten functional classes. The identified sequences include the expected RNA features like housekeeping sRNAs, metabolic riboswitches and RNA switches well known from studies on Bacillus subtilis as well as a multitude of completely new candidates for regulatory RNAs. An unexpectedly high number of 855 RNA features are encoded antisense to annotated protein and RNA genes, in addition to 461 independently transcribed small RNAs. These antisense transcripts contain molecules with a remarkable size range variation from 38 to 6348 base pairs in length. The genome of the type strain B. licheniformis DSM13 was completely reannotated using data obtained from RNA-Seq analyses and from public databases.

CONCLUSION

The hereby generated data-sets represent a solid amount of knowledge on the dynamic transcriptional activities during the investigated fermentation stages. The identified regulatory elements enable research on the understanding and the optimization of crucial metabolic activities during a productive fermentation of Bacillus licheniformis strains.

Association of RNAs with Bacillus subtilis Hfq

The prevalence and characteristics of small regulatory RNAs (sRNAs) have not been well characterized for Bacillus subtilis, an important model system for Gram-positive bacteria. However, B. subtilis was recently found to synthesize many candidate sRNAs during stationary phase. In the current study, we performed deep sequencing on Hfq-associated RNAs and found that a small subset of sRNAs associates with Hfq, an enigmatic RNA-binding protein that stabilizes sRNAs in Gram-negatives, but whose role is largely unknown in Gram-positive bacteria. We also found that Hfq associated with antisense RNAs, antitoxin transcripts, and many mRNA leaders. Several new candidate sRNAs and mRNA leader regions were also discovered by this analysis. Additionally, mRNA fragments overlapping with start or stop codons associated with Hfq, while, in contrast, relatively few full-length mRNAs were recovered. Deletion of hfq reduced the intracellular abundance of several representative sRNAs, suggesting that B. subtilis Hfq-sRNA interactions may be functionally significant in vivo. In general, we anticipate this catalog of Hfq-associated RNAs to serve as a resource in the functional characterization of Hfq in B. subtilis.

Computational discovery and RT-PCR validation of novel Burkholderia conserved and Burkholderia pseudomallei unique sRNAs

BACKGROUND

The sRNAs of bacterial pathogens are known to be involved in various cellular roles including environmental adaptation as well as regulation of virulence and pathogenicity. It is expected that sRNAs may also have similar functions for Burkholderia pseudomallei, a soil bacterium that can adapt to diverse environmental conditions, which causes the disease melioidosis and is also able to infect a wide variety of hosts.

RESULTS

By integrating several proven sRNA prediction programs into a computational pipeline, available Burkholderia spp. genomes were screened to identify sRNA gene candidates. Orthologous sRNA candidates were then identified via comparative analysis. From the total prediction, 21 candidates were found to have Rfam homologs. RT-PCR and sequencing of candidate sRNA genes of unknown functions revealed six putative sRNAs which were highly conserved in Burkholderia spp. and two that were unique to B. pseudomallei present in a normal culture conditions transcriptome. The validated sRNAs include potential cis-acting elements associated with the modulation of methionine metabolism and one B. pseudomallei-specific sRNA that is expected to bind to the Hfq protein.

CONCLUSIONS

The use of the pipeline developed in this study and subsequent comparative analysis have successfully aided in the discovery and shortlisting of sRNA gene candidates for validation. This integrated approach identified 29 B. pseudomallei sRNA genes - of which 21 have Rfam homologs and 8 are novel.

Acting antisense: plasmid- and chromosome-encoded sRNAs from Gram-positive bacteria

sRNAs that act by base pairing were first discovered in plasmids, phages and transposons, where they control replication, maintenance and transposition. Since 2001, however, computational searches were applied that led to the discovery of a plethora of sRNAs in bacterial chromosomes. Whereas the majority of these chromsome-encoded sRNAs have been investigated in Escherichia coli, Salmonella and other Gram-negative bacteria, only a few well-studied examples are known from Gram-positive bacteria. Here, the author summarizes our current knowledge on plasmid- and chromosome-encoded sRNAs from Gram-positive species, thereby focusing on regulatory mechanisms used by these RNAs and their biological role in complex networks. Furthermore, regulatory factors that control the expression of these RNAs will be discussed and differences between sRNAs from Gram-positive and Gram-negative bacteria highlighted. The main emphasis of this review is on sRNAs that act by base pairing (i.e., by an antisense mechanism). Thereby, both plasmid-encoded and chromosome-encoded sRNAs will be considered.

Regulation of ykrL (htpX) by Rok and YkrK, a novel type of regulator in Bacillus subtilis

Expression of ykrL of Bacillus subtilis, encoding a close homologue of the Escherichia coli membrane protein quality control protease HtpX, was shown to be upregulated under membrane protein overproduction stress. Using DNA affinity chromatography, two proteins were found to bind to the promoter region of ykrL: Rok, known as a repressor of competence and genes for extracytoplasmic functions, and YkrK, a novel type of regulator encoded by the gene adjacent to ykrL but divergently transcribed. Electrophoretic mobility shift assays showed Rok and YkrK binding to the ykrL promoter region as well as YkrK binding to the ykrK promoter region. Comparative bioinformatic analysis of the ykrL promoter regions in related Bacillus species revealed a consensus motif, which was demonstrated to be the binding site of YkrK. Deletion of rok and ykrK in a PykrL-gfp reporter strain showed that both proteins are repressors of ykrL expression. In addition, conditions which activated PykrL (membrane protein overproduction, dissipation of the membrane potential, and salt and phenol stress) point to the involvement of YkrL in membrane protein quality control.

Condition-dependent transcriptome reveals high-level regulatory architecture in Bacillus subtilis

Bacteria adapt to environmental stimuli by adjusting their transcriptomes in a complex manner, the full potential of which has yet to be established for any individual bacterial species. Here, we report the transcriptomes of Bacillus subtilis exposed to a wide range of environmental and nutritional conditions that the organism might encounter in nature. We comprehensively mapped transcription units (TUs) and grouped 2935 promoters into regulons controlled by various RNA polymerase sigma factors, accounting for ~66% of the observed variance in transcriptional activity. This global classification of promoters and detailed description of TUs revealed that a large proportion of the detected antisense RNAs arose from potentially spurious transcription initiation by alternative sigma factors and from imperfect control of transcription termination.

BsrG/SR4 from Bacillus subtilis--the first temperature-dependent type I toxin-antitoxin system

Here, we describe bsrG/SR4, a novel type I toxin-antitoxin system from the SPβ prophage region of the Bacillus subtilis chromosome. The 294-nucleotide bsrG RNA encodes a 38-amino-acid toxin, whereas SR4 is a 180-nucleotide antisense RNA that acts as the antitoxin. Both genes overlap by 123 nucleotides. BsrG expression increases at the onset of stationary phase. The sr4 promoter is 6- to 10-fold stronger than the bsrG promoter. Deletion of sr4 stabilizes bsrG mRNA and causes cell lysis on agar plates, which is due to the BsrG peptide and not the bsrG mRNA. SR4 overexpression could compensate cell lysis caused by overexpression of bsrG. SR4 interacts with the 3' UTR of bsrG RNA, thereby promoting its degradation. RNase III cleaves the bsrG RNA/SR4 duplex at position 185 of bsrG RNA, but is not essential for the function of the toxin-antitoxin system. Endoribonuclease Y and 3'-5' exoribonuclease R participate in the degradation of both bsrG RNA and SR4, whereas PnpA processes three SR4 precursors to the mature RNA. A heat shock at 48°C results in faster degradation and, therefore, significantly decreased amounts of bsrG RNA.

CsfG, a sporulation-specific, small non-coding RNA highly conserved in endospore formers

Endospore formation is a characteristic shared by some Bacilli and Clostridia that involves the creation of two cell types, the forespore and the mother cell. Hundreds of protein-encoding genes have been shown to be transcribed in a cell-specific fashion during this developmental process in Bacillus subtilis. We have used a phylogenetic profiling procedure to identify clusters of B. subtilis coding and non-coding sequences that co-occur in other endospore formers. One such cluster shows a strong bias for sporulation-related genes (42 % among 156 genes) and is enriched in potential non-coding RNAs. We have studied one RNA candidate, encoded in the ylbG-ylbH interval. In vivo analysis using a transcriptional fusion to the Escherichia coli lacZ gene demonstrates that this region of the chromosome contains a gene, csfG, encoding a 147-nucleotide RNA that is transcribed only during sporulation, specifically in the forespore. csfG is present in many endospore formers, mostly Bacilli and some Clostridia, whereas it is absent from bacteria that do not produce endospores. All CsfG RNAs contain a strongly conserved, pyrimidine-rich, central motif that overlaps a potential stem-loop structure. The remarkable conservation of this sequence in widely divergent bacteria suggests that it plays a conserved physiological role, presumably by interacting with an unidentified target in the forespore, where it contributes to the acquisition of the spore properties.

An integrated approach for finding overlooked genes in Shigella

BACKGROUND

The completion of numerous genome sequences introduced an era of whole-genome study. However, many genes are missed during genome annotation, including small RNAs (sRNAs) and small open reading frames (sORFs). In order to improve genome annotation, we aimed to identify novel sRNAs and sORFs in Shigella, the principal etiologic agents of bacillary dysentery.

METHODOLOGY/PRINCIPAL FINDINGS

We identified 64 sRNAs in Shigella, which were experimentally validated in other bacteria based on sequence conservation. We employed computer-based and tiling array-based methods to search for sRNAs, followed by RT-PCR and northern blots, to identify nine sRNAs in Shigella flexneri strain 301 (Sf301) and 256 regions containing possible sRNA genes. We found 29 candidate sORFs using bioinformatic prediction, array hybridization and RT-PCR verification. We experimentally validated 557 (57.9%) DOOR operon predictions in the chromosomes of Sf301 and 46 (76.7%) in virulence plasmid.We found 40 additional co-expressed gene pairs that were not predicted by DOOR.

CONCLUSIONS/SIGNIFICANCE

We provide an updated and comprehensive annotation of the Shigella genome. Our study increased the expected numbers of sORFs and sRNAs, which will impact on future functional genomics and proteomics studies. Our method can be used for large scale reannotation of sRNAs and sORFs in any microbe with a known genome sequence.

Stress response regulators identified through genome-wide transcriptome analysis of the (p)ppGpp-dependent response in Rhizobium etli

BACKGROUND

The alarmone (p)ppGpp mediates a global reprogramming of gene expression upon nutrient limitation and other stresses to cope with these unfavorable conditions. Synthesis of (p)ppGpp is, in most bacteria, controlled by RelA/SpoT (Rsh) proteins. The role of (p)ppGpp has been characterized primarily in Escherichia coli and several Gram-positive bacteria. Here, we report the first in-depth analysis of the (p)ppGpp-regulon in an α-proteobacterium using a high-resolution tiling array to better understand the pleiotropic stress phenotype of a relA/rsh mutant.

RESULTS

We compared gene expression of the Rhizobium etli wild type and rsh (previously rel) mutant during exponential and stationary phase, identifying numerous (p)ppGpp targets, including small non-coding RNAs. The majority of the 834 (p)ppGpp-dependent genes were detected during stationary phase. Unexpectedly, 223 genes were expressed (p)ppGpp-dependently during early exponential phase, indicating the hitherto unrecognized importance of (p)ppGpp during active growth. Furthermore, we identified two (p)ppGpp-dependent key regulators for survival during heat and oxidative stress and one regulator putatively involved in metabolic adaptation, namely extracytoplasmic function sigma factor EcfG2/PF00052, transcription factor CH00371, and serine protein kinase PrkA.

CONCLUSIONS

The regulatory role of (p)ppGpp in R. etli stress adaptation is far-reaching in redirecting gene expression during all growth phases. Genome-wide transcriptome analysis of a strain deficient in a global regulator, and exhibiting a pleiotropic phenotype, enables the identification of more specific regulators that control genes associated with a subset of stress phenotypes. This work is an important step toward a full understanding of the regulatory network underlying stress responses in α-proteobacteria.

Identification of regulatory RNAs in Bacillus subtilis

Post-transcriptional regulatory mechanisms are widespread in bacteria. Interestingly, current published data hint that some of these mechanisms may be non-random with respect to their phylogenetic distribution. Although small, trans-acting regulatory RNAs commonly occur in bacterial genomes, they have been better characterized in Gram-negative bacteria, leaving the impression that they may be less important for Firmicutes. It has been presumed that Gram-positive bacteria, in particular the Firmicutes, are likely to utilize cis-acting regulatory RNAs located within the 5' mRNA leader region more often than trans-acting regulatory RNAs. In this analysis we catalog, by a deep sequencing-based approach, both classes of regulatory RNA candidates for Bacillus subtilis, the model microorganism for Firmicutes. We successfully recover most of the known small RNA regulators while also identifying a greater number of new candidate RNAs. We anticipate these data to be a broadly useful resource for analysis of post-transcriptional regulatory strategies in B. subtilis and other Firmicutes.

Identification of novel non-coding small RNAs from Streptococcus pneumoniae TIGR4 using high-resolution genome tiling arrays

Background

The identification of non-coding transcripts in human, mouse, and Escherichia coli has revealed their widespread occurrence and functional importance in both eukaryotic and prokaryotic life. In prokaryotes, studies have shown that non-coding transcripts participate in a broad range of cellular functions like gene regulation, stress and virulence. However, very little is known about non-coding transcripts in Streptococcus pneumoniae (pneumococcus), an obligate human respiratory pathogen responsible for significant worldwide morbidity and mortality. Tiling microarrays enable genome wide mRNA profiling as well as identification of novel transcripts at a high-resolution.

Results

Here, we describe a high-resolution transcription map of the S. pneumoniae clinical isolate TIGR4 using genomic tiling arrays. Our results indicate that approximately 66% of the genome is expressed under our experimental conditions. We identified a total of 50 non-coding small RNAs (sRNAs) from the intergenic regions, of which 36 had no predicted function. Half of the identified sRNA sequences were found to be unique to S. pneumoniae genome. We identified eight overrepresented sequence motifs among sRNA sequences that correspond to sRNAs in different functional categories. Tiling arrays also identified approximately 202 operon structures in the genome.

Conclusions

In summary, the pneumococcal operon structures and novel sRNAs identified in this study enhance our understanding of the complexity and extent of the pneumococcal 'expressed' genome. Furthermore, the results of this study open up new avenues of research for understanding the complex RNA regulatory network governing S. pneumoniae physiology and virulence.

A dual-function sRNA from B. subtilis: SR1 acts as a peptide encoding mRNA on the gapA operon

Small non-coding RNAs (sRNAs) have been found to regulate gene expression in all three kingdoms of life. So far, relatively little is known about sRNAs from Gram-positive bacteria. SR1 is a regulatory sRNA from the Bacillus subtilis chromosome that inhibits by base-pairing translation initiation of ahrC mRNA encoding a transcriptional activator of the arginine catabolic operons. Here we present a novel target of SR1, the glycolytic gapA operon. Both microarray and Northern blot analyses show that the amount of gapA operon mRNA is significantly higher in the presence of SR1 when cells were grown in complex medium until stationary phase. Translational lacZ fusions and toeprinting analyses demonstrate that SR1 does not promote translation of gapA mRNA. By contrast, the half-life of gapA operon mRNA is strongly reduced in the sr1 knockout strain. SR1 does not act as a base-pairing sRNA on gapA operon mRNA. Instead, we demonstrate that the 39 aa peptide encoded by SR1, SR1P, is responsible for the effect of SR1 on the gapA operon. We show that SR1P binds GapA, thereby stabilizing the gapA operon mRNA by a hitherto unknown mechanism. SR1 is the first dual-function sRNA found in B. subtilis.

Comparative genomics reveals 104 candidate structured RNAs from bacteria, archaea, and their metagenomes

BACKGROUND

Structured noncoding RNAs perform many functions that are essential for protein synthesis, RNA processing, and gene regulation. Structured RNAs can be detected by comparative genomics, in which homologous sequences are identified and inspected for mutations that conserve RNA secondary structure.

RESULTS

By applying a comparative genomics-based approach to genome and metagenome sequences from bacteria and archaea, we identified 104 candidate structured RNAs and inferred putative functions for many of these. Twelve candidate metabolite-binding RNAs were identified, three of which were validated, including one reported herein that binds the coenzyme S-adenosylmethionine. Newly identified cis-regulatory RNAs are implicated in photosynthesis or nitrogen regulation in cyanobacteria, purine and one-carbon metabolism, stomach infection by Helicobacter, and many other physiological processes. A candidate riboswitch termed crcB is represented in both bacteria and archaea. Another RNA motif may control gene expression from 3'-untranslated regions of mRNAs, which is unusual for bacteria. Many noncoding RNAs that likely act in trans are also revealed, and several of the noncoding RNA candidates are found mostly or exclusively in metagenome DNA sequences.

CONCLUSIONS

This work greatly expands the variety of highly structured noncoding RNAs known to exist in bacteria and archaea and provides a starting point for biochemical and genetic studies needed to validate their biologic functions. Given the sustained rate of RNA discovery over several similar projects, we expect that far more structured RNAs remain to be discovered from bacterial and archaeal organisms.

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

Background

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.

Results

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.

Conclusions

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.

An overview of RNAs with regulatory functions in gram-positive bacteria

During the last decade, RNA molecules with regulatory functions on gene expression have benefited from a renewed interest. In bacteria, recent high throughput computational and experimental approaches have led to the discovery that 10-20% of all genes code for RNAs with critical regulatory roles in metabolic, physiological and pathogenic processes. The trans-acting RNAs comprise the noncoding RNAs, RNAs with a short open reading frame and antisense RNAs. Many of these RNAs act through binding to their target mRNAs while others modulate protein activity or target DNA. The cis-acting RNAs include regulatory regions of mRNAs that can respond to various signals. These RNAs often provide the missing link between sensing changing conditions in the environment and fine-tuning the subsequent biological responses. Information on their various functions and modes of action has been well documented for gram-negative bacteria. Here, we summarize the current knowledge of regulatory RNAs in gram-positive bacteria.

Identification and characterization of noncoding small RNAs in Streptococcus pneumoniae serotype 2 strain D39

We report a search for small RNAs (sRNAs) in the low-GC, gram-positive human pathogen Streptococcus pneumoniae. Based on bioinformatic analyses by Livny et al. (J. Livny, A. Brencic, S. Lory, and M. K. Waldor, Nucleic Acids Res. 34:3484-3493, 2006), we tested 40 candidates by Northern blotting and confirmed the expression of nine new and one previously reported (CcnA) sRNAs in strain D39. CcnA is one of five redundant sRNAs reported by Halfmann et al. (A. Halfmann, M. Kovacs, R. Hakenbeck, and R. Bruckner, Mol. Microbiol. 66:110-126, 2007) that are positively controlled by the CiaR response regulator. We characterized 3 of these 14 sRNAs: Spd-sr17 (144 nucleotides [nt]; decreased in stationary phase), Spd-sr37 (80 nt; strongly expressed in all growth phases), and CcnA (93 nt; induced by competence stimulatory peptide). Spd-sr17 and CcnA likely fold into structures containing single-stranded regions between hairpin structures, whereas Spd-sr37 forms a base-paired structure. Primer extension mapping and ectopic expression in deletion/insertion mutants confirmed the independent expression of the three sRNAs. Microarray analyses indicated that insertion/deletion mutants in spd-sr37 and ccnA exerted strong cis-acting effects on the transcription of adjacent genes, indicating that these sRNA regions are also cotranscribed in operons. Deletion or overexpression of the three sRNAs did not cause changes in growth, certain stress responses, global transcription, or virulence. Constitutive ectopic expression of CcnA reversed some phenotypes of D39 Delta ciaR mutants, but attempts to link CcnA to -E to comC as a target were inconclusive in ciaR(+) strains. These results show that S. pneumoniae, which lacks known RNA chaperones, expresses numerous sRNAs, but three of these sRNAs do not strongly affect common phenotypes or transcription patterns.

A community-curated consensual annotation that is continuously updated: the Bacillus subtilis centred wiki SubtiWiki

Bacillus subtilis is the model organism for Gram-positive bacteria, with a large amount of publications on all aspects of its biology. To facilitate genome annotation and the collection of comprehensive information on B. subtilis, we created SubtiWiki as a community-oriented annotation tool for information retrieval and continuous maintenance. The wiki is focused on the needs and requirements of scientists doing experimental work. This has implications for the design of the interface and for the layout of the individual pages. The pages can be accessed primarily by the gene designations. All pages have a similar flexible structure and provide links to related gene pages in SubtiWiki or to information in the World Wide Web. Each page gives comprehensive information on the gene, the encoded protein or RNA as well as information related to the current investigation of the gene/protein. The wiki has been seeded with information from key publications and from the most relevant general and B. subtilis-specific databases. We think that SubtiWiki might serve as an example for other scientific wikis that are devoted to the genes and proteins of one organism.

Database URL: The wiki can be accessed at http://subtiwiki.uni-goettingen.de/

The transcriptionally active regions in the genome of Bacillus subtilis

The majority of all genes have so far been identified and annotated systematically through in silico gene finding. Here we report the finding of 3662 strand-specific transcriptionally active regions (TARs) in the genome of Bacillus subtilis by the use of tiling arrays. We have measured the genome-wide expression during mid-exponential growth on rich (LB) and minimal (M9) medium. The identified TARs account for 77.3% of the genes as they are currently annotated and additionally we find 84 putative non-coding RNAs (ncRNAs) and 127 antisense transcripts. One ncRNA, ncr22, is predicted to act as a translational control on cstA and an antisense transcript was observed opposite the housekeeping sigma factor sigA. Through this work we have discovered a long conserved 3' untranslated region (UTR) in a group of membrane-associated genes that is predicted to fold into a large and highly stable secondary structure. One of the genes having this tail is efeN, which encodes a target of the twin-arginine translocase (Tat) protein translocation system.

CodY activates transcription of a small RNA in Bacillus subtilis

Regulatory small RNAs (sRNAs) in bacterial genomes have become a focus of research over the past 8 years. Whereas more than 100 such sRNAs have been found in Escherichia coli, relatively little is known about sRNAs in gram-positive bacteria. Using a computational approach, we identified two sRNAs in intergenic regions of the Bacillus subtilis genome, SR1 and SR2 (renamed BsrF). Recently, we demonstrated that SR1 inhibits the translation initiation of the transcriptional activator AhrC. Here, we describe detection of BsrF, its expression profile, and its regulation by CodY. Furthermore, we mapped the secondary structure of BsrF. BsrF is expressed in complex and minimal media in all growth phases in B. subtilis and, with a similar expression profile, also in Bacillus amyloliquefaciens. Neither overexpression nor deletion of bsrF affected the growth of B. subtilis. BsrF was found to be long-lived in complex and minimal media. Analysis of 13 putative transcription factor binding sites upstream of bsrF revealed only an effect for CodY. Here, we showed by using Northern blotting, lacZ reporter gene fusions, in vitro transcription, and DNase I footprinting that the transcription of bsrF is activated by CodY in the presence of branched-chain amino acids and GTP. Furthermore, BsrF transcription was increased 1.5- to 2-fold by glucose in the presence of branched-chain amino acids, and this increase was independent of the known glucose-dependent regulators. BsrF is the second target for which transcriptional activation by CodY has been discovered.

From a consortium sequence to a unified sequence: the Bacillus subtilis 168 reference genome a decade later

Comparative genomics is the cornerstone of identification of gene functions. The immense number of living organisms precludes experimental identification of functions except in a handful of model organisms. The bacterial domain is split into large branches, among which the Firmicutes occupy a considerable space. Bacillus subtilis has been the model of Firmicutes for decades and its genome has been a reference for more than 10 years. Sequencing the genome involved more than 30 laboratories, with different expertises, in a attempt to make the most of the experimental information that could be associated with the sequence. This had the expected drawback that the sequencing expertise was quite varied among the groups involved, especially at a time when sequencing genomes was extremely hard work. The recent development of very efficient, fast and accurate sequencing techniques, in parallel with the development of high-level annotation platforms, motivated the present resequencing work. The updated sequence has been reannotated in agreement with the UniProt protein knowledge base, keeping in perspective the split between the paleome (genes necessary for sustaining and perpetuating life) and the cenome (genes required for occupation of a niche, suggesting here that B. subtilis is an epiphyte). This should permit investigators to make reliable inferences to prepare validation experiments in a variety of domains of bacterial growth and development as well as build up accurate phylogenies.