The transcriptionally active regions in the genome of Bacillus subtilis

Genome-wide annotation of transcript boundaries using bacterial Rend-seq datasets

Accurate annotation to single-nucleotide resolution of the transcribed regions in genomes is key to optimally analyse RNA-seq data, understand regulatory events and for the design of experiments. However, currently most genome annotations provided by GenBank generally lack information about untranslated regions. Additionally, information regarding genomic locations of non-coding RNAs, such as sRNAs, or anti-sense RNAs is frequently missing. To provide such information, diverse RNA-seq technologies, such as Rend-seq, have been developed and applied to many bacterial species. However, incorporating this vast amount of information into annotation files has been limited and is bioinformatically challenging, resulting in UTRs and other non-coding elements being overlooked or misrepresented. To overcome this problem, we present pyRAP (python Rend-seq Annotation Pipeline), a software package that analyses Rend-seq datasets to accurately resolve transcript boundaries genome-wide. We report the use of pyRAP to find novel transcripts, transcript isoforms, and RNase-dependent sRNA processing events. In Bacillus subtilis we uncovered 63 novel transcripts and provide genomic coordinates with single-nucleotide resolution for 2218 5'UTRs, 1864 3'UTRs and 161 non-coding RNAs. In Escherichia coli, we report 117 novel transcripts, 2429 5'UTRs, 1619 3'UTRs and 91 non-coding RNAs, and in Staphylococcus aureus, 16 novel transcripts, 664 5'UTRs, 696 3'UTRs, and 81 non-coding RNAs. Finally, we use pyRAP to produce updated annotation files for B. subtilis 168, E. coli K-12 MG1655, and S. aureus 8325 for use in the wider microbial genomics research community.

Phylogenomic characterization and pangenomic insights into the surfactin-producing bacteria Bacillus subtilis strain RI4914

Bacillus subtilis is a versatile bacterial species able to produce surfactin, a lipopeptide biosurfactant. We carried out the phylogenomic characterization and pangenomic analyses using available B. subtilis complete genomes. Also, we report the whole genome of the biosurfactant-producing B. subtilis strain RI4914 that was isolated from effluent water from an oil exploration field. We applied a hybrid sequencing approach using both long- and short-read sequencing technologies to generate a highly accurate, single-chromosome genome. The pangenomics analysis of 153 complete genomes classified as B. subtilis retrieved from the NCBI shows an open pangenome composed of 28,511 accessory genes, which agrees with the high genetic plasticity of the species. Also, this analysis suggests that surfactin production is a common trait shared by members of this species since the srfA operon is highly conserved among the B. subtilis strains found in most of the assemblies available. Finally, increased surfactin production corroborates the higher srfAA gene expression in B. subtilis strain RI4914.

Model-Based Design of Synthetic Antisense RNA for Predictable Gene Repression

Our enhanced understanding of RNA folding and function has increased the use of small RNA regulators. Among these RNA regulators, synthetic antisense RNA (asRNA) is designed to contain an RNA sequence complementary to the target mRNA sequence, and the formation of double-stranded RNA (dsRNA) facilitates gene repression due to dsRNA degradation or prevention of ribosome access to the mRNA. Despite the simple complementarity rule, however, predictably tunable repression has been challenging when synthetic asRNAs are used. Here, the protocol for model-based asRNA design is described. This model can predict synthetic asRNA-mediated repression efficiency using two parameters: the change in free energy of complex formation (ΔG_CF) and percent mismatch of the target binding region (TBR). The model has been experimentally validated in both Gram-positive and Gram-negative bacteria as well as for target genes in both plasmids and chromosomes. These asRNAs can be created by simply replacing the TBR sequence with one that is complementary to the target mRNA sequence of interest. In principle, this protocol can be applied to design and build asRNAs for predictable gene repression in various contexts, including multiple target genes and organisms, making asRNAs predictably tunable regulators for broad applications.

A High-Content Microscopy Screening Identifies New Genes Involved in Cell Width Control in Bacillus subtilis

How cells control their shape and size is a fundamental question of biology. In most bacteria, cell shape is imposed by the peptidoglycan (PG) polymeric meshwork that surrounds the cell. Thus, bacterial cell morphogenesis results from the coordinated action of the proteins assembling and degrading the PG shell. Remarkably, during steady-state growth, most bacteria maintain a defined shape along generations, suggesting that error-proof mechanisms tightly control the process. In the rod-shaped model for the Gram-positive bacterium Bacillus subtilis, the average cell length varies as a function of the growth rate, but the cell diameter remains constant throughout the cell cycle and across growth conditions. Here, in an attempt to shed light on the cellular circuits controlling bacterial cell width, we developed a screen to identify genetic determinants of cell width in B. subtilis. Using high-content screening (HCS) fluorescence microscopy and semiautomated measurement of single-cell dimensions, we screened a library of ∼4,000 single knockout mutants. We identified 13 mutations significantly altering cell diameter, in genes that belong to several functional groups. In particular, our results indicate that metabolism plays a major role in cell width control in B. subtilis.

IMPORTANCE Bacterial shape is primarily dictated by the external cell wall, a vital structure that, as such, is the target of countless antibiotics. Our understanding of how bacteria synthesize and maintain this structure is therefore a cardinal question for both basic and applied research. Bacteria usually multiply from generation to generation while maintaining their progenies with rigorously identical shapes. This implies that the bacterial cells constantly monitor and maintain a set of parameters to ensure this perpetuation. Here, our study uses a large-scale microscopy approach to identify at the whole-genome level, in a model bacterium, the genes involved in the control of one of the most tightly controlled cellular parameters, the cell width.

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.

Adaptive gene misregulation

Because gene expression is important for evolutionary adaptation, its misregulation is an important cause of maladaptation. A misregulated gene can be incorrectly silent ("off") when a transcription factor (TF) that is required for its activation does not binds its regulatory region. Conversely, a misregulated gene can be incorrectly active ("on") when a TF not normally involved in its activation binds its regulatory region, a phenomenon also known as regulatory crosstalk. DNA mutations that destroy or create TF binding sites on DNA are an important source of misregulation and crosstalk. Although misregulation reduces fitness in an environment to which an organism is well-adapted, it may become adaptive in a new environment. Here, I derive simple yet general mathematical expressions that delimit the conditions under which misregulation can be adaptive. These expressions depend on the strength of selection against misregulation, on the fraction of DNA sequence space filled with TF binding sites, and on the fraction of genes that must be expressed for optimal adaptation. I then use empirical data from RNA sequencing, protein-binding microarrays, and genome evolution, together with population genetic simulations to ask when these conditions are likely to be met. I show that they can be met under realistic circumstances, but these circumstances may vary among organisms and environments. My analysis provides a framework in which improved theory and data collection can help us demonstrate the role of misregulation in adaptation. It also shows that misregulation, like DNA mutation, is one of life's many imperfections that can help propel Darwinian evolution.

Quasi-essentiality of RNase Y in Bacillus subtilis is caused by its critical role in the control of mRNA homeostasis

RNA turnover is essential in all domains of life. The endonuclease RNase Y (rny) is one of the key components involved in RNA metabolism of the model organism Bacillus subtilis. Essentiality of RNase Y has been a matter of discussion, since deletion of the rny gene is possible, but leads to severe phenotypic effects. In this work, we demonstrate that the rny mutant strain rapidly evolves suppressor mutations to at least partially alleviate these defects. All suppressor mutants had acquired a duplication of an about 60 kb long genomic region encompassing genes for all three core subunits of the RNA polymerase—α, β, β′. When the duplication of the RNA polymerase genes was prevented by relocation of the rpoA gene in the B. subtilis genome, all suppressor mutants carried distinct single point mutations in evolutionary conserved regions of genes coding either for the β or β’ subunits of the RNA polymerase that were not tolerated by wild type bacteria. In vitro transcription assays with the mutated polymerase variants showed a severe decrease in transcription efficiency. Altogether, our results suggest a tight cooperation between RNase Y and the RNA polymerase to establish an optimal RNA homeostasis in B. subtilis cells.

Identification of an RNA sponge that controls the RoxS riboregulator of central metabolism in Bacillus subtilis

sRNAs are a taxonomically-restricted but transcriptomically-abundant class of post-transcriptional regulators. While of major importance for adaption to the environment, we currently lack global-scale methodology enabling target identification, especially in species without known RNA hub proteins (e.g. Hfq). Using psoralen RNA cross-linking and Illumina-sequencing we identify RNA-RNA interacting pairs in vivo in Bacillus subtilis, resolving previously well-described interactants. Although sRNA-sRNA pairings are rare (compared with sRNA-mRNA), we identify a robust example involving the conserved sRNA RoxS and an unstudied sRNA RosA (Regulator of sRNA A). We show RosA to be the first confirmed RNA sponge described in a Gram-positive bacterium. RosA interacts with at least two sRNAs, RoxS and FsrA. The RosA/RoxS interaction not only affects the levels of RoxS but also its processing and regulatory activity. We also found that the transcription of RosA is repressed by CcpA, the key regulator of carbon-metabolism in B. subtilis. Since RoxS is already known to be transcriptionally controlled by malate via the transcriptional repressor Rex, its post-transcriptional regulation by CcpA via RosA places RoxS in a key position to control central metabolism in response to varying carbon sources.

BSGatlas: a unified Bacillus subtilis genome and transcriptome annotation atlas with enhanced information access

A large part of our current understanding of gene regulation in Gram-positive bacteria is based on Bacillus subtilis, as it is one of the most well studied bacterial model systems. The rapid growth in data concerning its molecular and genomic biology is distributed across multiple annotation resources. Consequently, the interpretation of data from further B. subtilis experiments becomes increasingly challenging in both low- and large-scale analyses. Additionally, B. subtilis annotation of structured RNA and non-coding RNA (ncRNA), as well as the operon structure, is still lagging behind the annotation of the coding sequences. To address these challenges, we created the B. subtilis genome atlas, BSGatlas, which integrates and unifies multiple existing annotation resources. Compared to any of the individual resources, the BSGatlas contains twice as many ncRNAs, while improving the positional annotation for 70 % of the ncRNAs. Furthermore, we combined known transcription start and termination sites with lists of known co-transcribed gene sets to create a comprehensive transcript map. The combination with transcription start/termination site annotations resulted in 717 new sets of co-transcribed genes and 5335 untranslated regions (UTRs). In comparison to existing resources, the number of 5' and 3' UTRs increased nearly fivefold, and the number of internal UTRs doubled. The transcript map is organized in 2266 operons, which provides transcriptional annotation for 92 % of all genes in the genome compared to the at most 82 % by previous resources. We predicted an off-target-aware genome-wide library of CRISPR-Cas9 guide RNAs, which we also linked to polycistronic operons. We provide the BSGatlas in multiple forms: as a website (https://rth.dk/resources/bsgatlas/), an annotation hub for display in the UCSC genome browser, supplementary tables and standardized GFF3 format, which can be used in large scale -omics studies. By complementing existing resources, the BSGatlas supports analyses of the B. subtilis genome and its molecular biology with respect to not only non-coding genes but also genome-wide transcriptional relationships of all genes.

Bacterial 3'UTRs: A Useful Resource in Post-transcriptional Regulation

Bacterial messenger RNAs (mRNAs) are composed of 5′ and 3′ untranslated regions (UTRs) that flank the coding sequences (CDSs). In eukaryotes, 3′UTRs play key roles in post-transcriptional regulatory mechanisms. Shortening or deregulation of these regions is associated with diseases such as cancer and metabolic disorders. Comparatively, little is known about the functions of 3′UTRs in bacteria. Over the past few years, 3′UTRs have emerged as important players in the regulation of relevant bacterial processes such as virulence, iron metabolism, and biofilm formation. This MiniReview is an update for the different 3′UTR-mediated mechanisms that regulate gene expression in bacteria. Some of these include 3′UTRs that interact with the 5′UTR of the same transcript to modulate translation, 3′UTRs that are targeted by specific ribonucleases, RNA-binding proteins and small RNAs (sRNAs), and 3′UTRs that act as reservoirs of trans-acting sRNAs, among others. In addition, recent findings regarding a differential evolution of bacterial 3′UTRs and its impact in the species-specific expression of orthologous genes are also discussed.

Endogenous antisense walR RNA modulates biofilm organization and pathogenicity of Enterococcus faecalis

Enterococcus faecalis (E. faecalis) is regarded as the major pathogen for persistent periapical periodontitis. The aim of the present study was to investigate the role of antisense walR RNA in the regulation of adjacent downstream genes. Reverse transcription-PCR assays were performed to validate walR. Adjacent downstream genes walK, EF1195, EF1196, and EF1197 were co-transcribed and detect antisense walR RNA. Northern blotting and 5'-rapid amplification of cDNA ends (5'-RACE) assays were conducted to detect and confirm a novel walR antisense (ASwalR) RNA. ASwalR overexpression mutants were constructed, and the biofilm biomass was determined using a crystal violet microtiter assay. The present study detected and confirmed a 550-bp noncoding antisense RNA with the potential to attenuate the activities of the essential response regulator WalR. The levels of antisense walR RNA transcripts were inversely associated with the production of WalR protein. It was showed that overexpression of ASwalR leads to reduced biofilm formation and exopolysaccharide synthesis. Furthermore, the pathogenicity of E. faecalis was markedly decreased by ASwalR overexpression in an in vivo periapical periodontitis model. In summary, the present study detected a novel antisense walR RNA that leads to a reduction in biofilm formation and the pathogenicity of E. faecalis. Collectively, the data suggest a role for ASwalR as a post-transcriptional modulator of the WalR regulator in E. faecalis.

Intermolecular Communication in Bacillus subtilis: RNA-RNA, RNA-Protein and Small Protein-Protein Interactions

In bacterial cells we find a variety of interacting macromolecules, among them RNAs and proteins. Not only small regulatory RNAs (sRNAs), but also small proteins have been increasingly recognized as regulators of bacterial gene expression. An average bacterial genome encodes between 200 and 300 sRNAs, but an unknown number of small proteins. sRNAs can be cis- or trans-encoded. Whereas cis-encoded sRNAs interact only with their single completely complementary mRNA target transcribed from the opposite DNA strand, trans-encoded sRNAs are only partially complementary to their numerous mRNA targets, resulting in huge regulatory networks. In addition to sRNAs, uncharged tRNAs can interact with mRNAs in T-box attenuation mechanisms. For a number of sRNA-mRNA interactions, the stability of sRNAs or translatability of mRNAs, RNA chaperones are required. In Gram-negative bacteria, the well-studied abundant RNA-chaperone Hfq fulfils this role, and recently another chaperone, ProQ, has been discovered and analyzed in this respect. By contrast, evidence for RNA chaperones or their role in Gram-positive bacteria is still scarce, but CsrA might be such a candidate. Other RNA-protein interactions involve tmRNA/SmpB, 6S RNA/RNA polymerase, the dual-function aconitase and protein-bound transcriptional terminators and antiterminators. Furthermore, small proteins, often missed in genome annotations and long ignored as potential regulators, can interact with individual regulatory proteins, large protein complexes, RNA or the membrane. Here, we review recent advances on biological role and regulatory principles of the currently known sRNA-mRNA interactions, sRNA-protein interactions and small protein-protein interactions in the Gram-positive model organism Bacillus subtilis. We do not discuss RNases, ribosomal proteins, RNA helicases or riboswitches.

Recommendations for bacterial ribosome profiling experiments based on bioinformatic evaluation of published data

Ribosome profiling (RIBO-Seq) has improved our understanding of bacterial translation, including finding many unannotated genes. However, protocols for RIBO-Seq and corresponding data analysis are not yet standardized. Here, we analyzed 48 RIBO-Seq samples from nine studies of Escherichia coli K12 grown in lysogeny broth medium and particularly focused on the size-selection step. We show that for conventional expression analysis, a size range between 22 and 30 nucleotides is sufficient to obtain protein-coding fragments, which has the advantage of removing many unwanted rRNA and tRNA reads. More specific analyses may require longer reads and a corresponding improvement in rRNA/tRNA depletion. There is no consensus about the appropriate sequencing depth for RIBO-Seq experiments in prokaryotes, and studies vary significantly in total read number. Our analysis suggests that 20 million reads that are not mapping to rRNA/tRNA are required for global detection of translated annotated genes. We also highlight the influence of drug-induced ribosome stalling, which causes bias at translation start sites. The resulting accumulation of reads at the start site may be especially useful for detecting weakly expressed genes. As different methods suit different questions, it may not be possible to produce a "one-size-fits-all" ribosome profiling data set. Therefore, experiments should be carefully designed in light of the scientific questions of interest. We propose some basic characteristics that should be reported with any new RIBO-Seq data sets. Careful attention to the factors discussed should improve prokaryotic gene detection and the comparability of ribosome profiling data sets.

The Many Facets of the Small Non-coding RNA RsaE (RoxS) in Metabolic Niche Adaptation of Gram-Positive Bacteria

Small regulatory RNAs (sRNAs) are increasingly recognized as players in the complex regulatory networks governing bacterial gene expression. RsaE (synonym RoxS) is an sRNA that is highly conserved in bacteria of the Bacillales order. Recent analyses in Bacillus subtilis, Staphylococcus aureus and Staphylococcus epidermidis identified RsaE/RoxS as a potent riboregulator of central carbon metabolism and energy balance with many molecular RsaE/RoxS functions and targets being shared across species. Similarities and species-specific differences in cellular processes modulated by RsaE/RoxS suggest that this sRNA plays a prominent role in the adaptation of Gram-positive bacteria to niches with varying nutrient availabilities and environmental cues. This review summarizes recent findings on the molecular function of RsaE/RoxS and its interaction with mRNA targets. Special emphasis will be on the integration of RsaE/RoxS into metabolic regulatory circuits and, derived from this, the role of RsaE/RoxS as a putative driver to generate phenotypic heterogeneity in bacterial populations. In this respect, we will particularly discuss heterogeneous RsaE expression in S. epidermidis biofilms and its possible contribution to metabolic niche diversification, programmed bacterial lysis and biofilm matrix production.

Activity of Streptococcus mutans VicR Is Modulated by Antisense RNA

The VicRK 2-component system of Streptococcus mutans regulates genes associated with cell wall biogenesis and biofilm formation. A putative RNase III-encoding gene ( rnc) is located downstream from the vicRKX operon. The goals of this study were to investigate the potential role of VicR in the regulation of adjacent downstream genes and evaluate transcription levels of vicR during planktonic and biofilm growth. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to investigate whether vicRKX and adjacent downstream genes were cotranscribed. Binding of purified recombinant VicR protein to promoter regions of vicR, rnc, and syfA genes was confirmed by electrophoretic mobility shift assay and by chromatin immunoprecipitation analyses. VicR antisense (AS vicR) RNA was detected by Northern blotting and qRT-PCR assays. AS vicR overexpression mutants were constructed, and the biofilm biomass was determined by crystal violet microtiter assay. Adjacent downstream genes rnc, smc, syfA, smu.1511, and syfB were cotranscribed with vicRKX. The predicted promoter regions of vicR, rnc, and syfA genes were directly regulated by VicR. An AS vicR RNA transcript was detected upstream of the rnc gene. Expression of the AS vicR RNA transcript was elevated in planktonic cultures and repressed during biofilm growth. In addition, Western blot data showed that expression of the VicR protein decreased by 35% in planktonic as compared with biofilm cultures. Furthermore, we show that overexpression of AS vicR led to a reduction in biofilm formation. The downstream genes rnc, smc, syfA, smu.1511, and syfB are cotranscribed with vicRKX. VicR is autophosphorylated, and rnc and syfA are directly regulated by VicR. Expression of VicR protein correlated inversely with different levels of AS vicR RNA transcript and growth conditions. The biofilm biomass decreased in the AS vicR overexpression mutant. These data suggest a role for the AS vicR RNA transcript in posttranscriptional regulation of VicR protein production in S. mutans.

Antisense transcription in Pseudomonas aeruginosa

A large number of antisense transcripts have been detected in diverse microbial genomes and considerable effort has been devoted to elucidating the functional role of antisense transcription. In this study, we reanalysed extensive RNA sequencing data from the opportunistic pathogen Pseudomonas aeruginosa and found that the majority of genes have a propensity for antisense transcription. Although antisense transcripts were found in more than 80 % of the genes of the P. aeruginosa genome, the majority of sequencing reads were mapping sense and only a minority (<2 %) were mapping antisense to genes. Similarly to the sense expression levels, the antisense expression levels varied under different environmental conditions, with the sense and antisense expression levels often being inversely regulated and modulated by the activity of alternative sigma factors. Environment-modulated antisense transcription showed a bias towards being antisense to genes within regions of genomic plasticity and to those encoding small regulatory RNAs. In the future, the validation and functional characterization of antisense transcripts, and novel transcripts that are antisense to small regulatory RNAs in particular, have the potential to contribute to our understanding of the various levels of transcriptional regulation and its dynamics in the bacterial pathogen P. aeruginosa.

Untranslated Parts of Genes Interpreted: Making Heads or Tails of High-Throughput Transcriptomic Data via Computational Methods: Computational methods to discover and quantify isoforms with alternative untranslated regions

In this review we highlight the importance of defining the untranslated parts of transcripts, and present a number of computational approaches for the discovery and quantification of alternative transcription start and poly-adenylation events in high-throughput transcriptomic data. The fate of eukaryotic transcripts is closely linked to their untranslated regions, which are determined by the position at which transcription starts and ends at a genomic locus. Although the extent of alternative transcription starts and alternative poly-adenylation sites has been revealed by sequencing methods focused on the ends of transcripts, the application of these methods is not yet widely adopted by the community. We suggest that computational methods applied to standard high-throughput technologies are a useful, albeit less accurate, alternative to the expertise-demanding 5' and 3' sequencing and they are the only option for analysing legacy transcriptomic data. We review these methods here, focusing on technical challenges and arguing for the need to include better normalization of the data and more appropriate statistical models of the expected variation in the signal.

Genomic and transcriptomic landscape of Escherichia coli BL21(DE3)

Escherichia coli BL21(DE3) has long served as a model organism for scientific research, as well as a workhorse for biotechnology. Here we present the most current genome annotation of E. coli BL21(DE3) based on the transcriptome structure of the strain that was determined for the first time. The genome was annotated using multiple automated pipelines and compared to the current genome annotation of the closely related strain, E. coli K-12. High-resolution tiling array data of E. coli BL21(DE3) from several different stages of cell growth in rich and minimal media were analyzed to characterize the transcriptome structure and to provide supporting evidence for open reading frames. This new integrated analysis of the genomic and transcriptomic structure of E. coli BL21(DE3) has led to the correction of translation initiation sites for 88 coding DNA sequences and provided updated information for most genes. Additionally, 37 putative genes and 66 putative non-coding RNAs were also identified. The panoramic landscape of the genome and transcriptome of E. coli BL21(DE3) revealed here will allow us to better understand the fundamental biology of the strain and also advance biotechnological applications in industry.

Revisiting the in vivo GlnR-binding sites at the genome scale in Bacillus subtilis

BACKGROUND

In Bacillus subtilis, two major transcriptional factors, GlnR and TnrA, are involved in a sophisticated network of adaptive responses to nitrogen availability. GlnR was reported to repress the transcription of the glnRA, tnrA and ureABC operons under conditions of excess nitrogen. As GlnR and TnrA regulators share the same DNA binding motifs, a genome-wide mapping of in vivo GlnR-binding sites was still needed to clearly define the set of GlnR/TnrA motifs directly bound by GlnR.

METHODS

We used chromatin immunoprecipitation coupled with hybridization to DNA tiling arrays (ChIP-on-chip) to identify the GlnR DNA-binding sites, in vivo, at the genome scale.

RESULTS

We provide evidence that GlnR binds reproducibly to 61 regions on the chromosome. Among those, 20 regions overlap the previously defined in vivo TnrA-binding sites. In combination with real-time in vivo transcriptional profiling using firefly luciferase, we identified the alsT gene as a new member of the GlnR regulon. Additionally, we characterized the GlnR secondary regulon, which is composed of promoter regions harboring a GlnR/TnrA box and bound by GlnR in vivo. However, the growth conditions revealing a GlnR-dependent regulation for this second category of genes are still unknown.

CONCLUSIONS

Our findings show an extended overlap between the GlnR and TnrA in vivo binding sites. This could allow efficient and fine tuning of gene expression in response to nitrogen availability. GlnR appears to be part of complex transcriptional regulatory networks, which involves interactions between different regulatory proteins.

Common and phylogenetically widespread coding for peptides by bacterial small RNAs

Background

While eukaryotic noncoding RNAs have recently received intense scrutiny, it is becoming clear that bacterial transcription is at least as pervasive. Bacterial small RNAs and antisense RNAs (sRNAs) are often assumed to be noncoding, due to their lack of long open reading frames (ORFs). However, there are numerous examples of sRNAs encoding for small proteins, whether or not they also have a regulatory role at the RNA level.

Methods

Here, we apply flexible machine learning techniques based on sequence features and comparative genomics to quantify the prevalence of sRNA ORFs under natural selection to maintain protein-coding function in 14 phylogenetically diverse bacteria. Importantly, we quantify uncertainty in our predictions, and follow up on them using mass spectrometry proteomics and comparison to datasets including ribosome profiling.

Results

A majority of annotated sRNAs have at least one ORF between 10 and 50 amino acids long, and we conservatively predict that 409±191.7 unannotated sRNA ORFs are under selection to maintain coding (mean estimate and 95% confidence interval), an average of 29 per species considered here. This implies that overall at least 10.3±0.5% of sRNAs have a coding ORF, and in some species around 20% do. 165±69 of these novel coding ORFs have some antisense overlap to annotated ORFs. As experimental validation, many of our predictions are translated in published ribosome profiling data and are identified via mass spectrometry shotgun proteomics. B. subtilis sRNAs with coding ORFs are enriched for high expression in biofilms and confluent growth, and S. pneumoniae sRNAs with coding ORFs are involved in virulence. sRNA coding ORFs are enriched for transmembrane domains and many are predicted novel components of type I toxin/antitoxin systems.

Conclusions

We predict over two dozen new protein-coding genes per bacterial species, but crucially also quantified the uncertainty in this estimate. Our predictions for sRNA coding ORFs, along with predicted novel type I toxins and tools for sorting and visualizing genomic context, are freely available in a user-friendly format at http://disco-bac.web.pasteur.fr. We expect these easily-accessible predictions to be a valuable tool for the study not only of bacterial sRNAs and type I toxin-antitoxin systems, but also of bacterial genetics and genomics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3932-y) contains supplementary material, which is available to authorized users.

Small RNA mediated repression of subtilisin production in Bacillus licheniformis

The species Bacillus licheniformis includes important strains that are used in industrial production processes. Currently the physiological model used to adapt these processes is based on the closely related model organism B. subtilis. In this study we found that both organisms reveal significant differences in the regulation of subtilisin, their main natural protease and a product of industrial fermentation processes. We identified and characterized a novel antisense sRNA AprAs, which represents an RNA based repressor of apr, the gene encoding for the industrial relevant subtilisin protease. Reduction of the AprAs level leads to an enhanced proteolytic activity and an increase of Apr protein expression in the mutant strain. A vector based complementation of the AprAs deficient mutant confirmed this effect and demonstrated the necessity of cis transcription for full efficiency. A comparative analysis of the corresponding genome loci from B. licheniformis and B. subtilis revealed the absence of an aprAs promoter in B. subtilis and indicates that AprAs is a B. licheniformis species specific phenomenon. The discovery of AprAs is of great biotechnological interest since subtilisin Carlsberg is one of the main products of industrial fermentation by B. licheniformis.

Regulatory 3' Untranslated Regions of Bacterial mRNAs

The untranslated regions (UTRs) of mRNA contain important features that are relevant to the post-transcriptional and translational regulation of gene expression. Most studies of bacterial UTRs have focused on the 5′regions; however, 3′UTRs have recently emerged as a new class of post-transcriptional regulatory elements. 3′UTRs were found to regulate the decay and translation initiation in their own mRNAs. In addition, 3′UTRs constitute a rich reservoir of small regulatory RNAs, regulating target gene expression. In the current review, we describe several recently discovered examples of bacterial regulatory 3′UTRs, discuss their modes of action, and illustrate how they facilitate gene regulation in various environments.

Deep sequencing approaches for the analysis of prokaryotic transcriptional boundaries and dynamics

The identification of the protein-coding regions of a genome is straightforward due to the universality of start and stop codons. However, the boundaries of the transcribed regions, conditional operon structures, non-coding RNAs and the dynamics of transcription, such as pausing of elongation, are non-trivial to identify, even in the comparatively simple genomes of prokaryotes. Traditional methods for the study of these areas, such as tiling arrays, are noisy, labour-intensive and lack the resolution required for densely-packed bacterial genomes. Recently, deep sequencing has become increasingly popular for the study of the transcriptome due to its lower costs, higher accuracy and single nucleotide resolution. These methods have revolutionised our understanding of prokaryotic transcriptional dynamics. Here, we review the deep sequencing and data analysis techniques that are available for the study of transcription in prokaryotes, and discuss the bioinformatic considerations of these analyses.

Genome-wide identification and characterization of small RNAs in Rhodobacter capsulatus and identification of small RNAs affected by loss of the response regulator CtrA

Small non-coding RNAs (sRNAs) are involved in the control of numerous cellular processes through various regulatory mechanisms, and in the past decade many studies have identified sRNAs in a multitude of bacterial species using RNA sequencing (RNA-seq). Here, we present the first genome-wide analysis of sRNA sequencing data in Rhodobacter capsulatus, a purple nonsulfur photosynthetic alphaproteobacterium. Using a recently developed bioinformatics approach, sRNA-Detect, we detected 422 putative sRNAs from R. capsulatus RNA-seq data. Based on their sequence similarity to sRNAs in a sRNA collection, consisting of published putative sRNAs from 23 additional bacterial species, and RNA databases, the sequences of 124 putative sRNAs were conserved in at least one other bacterial species; and, 19 putative sRNAs were assigned a predicted function. We bioinformatically characterized all putative sRNAs and applied machine learning approaches to calculate the probability of a nucleotide sequence to be a bona fide sRNA. The resulting quantitative model was able to correctly classify 95.2% of sequences in a validation set. We found that putative cis-targets for antisense and partially overlapping sRNAs were enriched with protein-coding genes involved in primary metabolic processes, photosynthesis, compound binding, and with genes forming part of macromolecular complexes. We performed differential expression analysis to compare the wild type strain to a mutant lacking the response regulator CtrA, an important regulator of gene expression in R. capsulatus, and identified 18 putative sRNAs with differing levels in the two strains. Finally, we validated the existence and expression patterns of four novel sRNAs by Northern blot analysis.

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.

RodZ and PgsA Play Intertwined Roles in Membrane Homeostasis of Bacillus subtilis and Resistance to Weak Organic Acid Stress

Weak organic acids like sorbic and acetic acid are widely used to prevent growth of spoilage organisms such as Bacilli. To identify genes involved in weak acid stress tolerance we screened a transposon mutant library of Bacillus subtilis for sorbic acid sensitivity. Mutants of the rodZ (ymfM) gene were found to be hypersensitive to the lipophilic weak organic acid. RodZ is involved in determining the cell's rod-shape and believed to interact with the bacterial actin-like MreB cytoskeleton. Since rodZ lies upstream in the genome of the essential gene pgsA (phosphatidylglycerol phosphate synthase) we hypothesized that expression of the latter might also be affected in rodZ mutants and hence contribute to the phenotype observed. We show that both genes are co-transcribed and that both the rodZ::mini-Tn10 mutant and a conditional pgsA mutant, under conditions of minimal pgsA expression, were sensitive to sorbic and acetic acid. Both strains displayed a severely altered membrane composition. Compared to the wild-type strain, phosphatidylglycerol and cardiolipin levels were lowered and the average acyl chain length was elongated. Induction of rodZ expression from a plasmid in our transposon mutant led to no recovery of weak acid susceptibility comparable to wild-type levels. However, pgsA overexpression in the same mutant partly restored sorbic acid susceptibility and fully restored acetic acid sensitivity. A construct containing both rodZ and pgsA as on the genome led to some restored growth as well. We propose that RodZ and PgsA play intertwined roles in membrane homeostasis and tolerance to weak organic acid stress.

AhpA is a peroxidase expressed during biofilm formation in Bacillus subtilis

Organisms growing aerobically generate reactive oxygen species such as hydrogen peroxide. These reactive oxygen molecules damage enzymes and DNA, potentially causing cell death. In response, Bacillus subtilis produces at least nine potential peroxide-scavenging enzymes; two belong to the alkylhydroperoxide reductase (Ahp) class of peroxidases. Here, we explore the role of one of these Ahp homologs, AhpA. While previous studies demonstrated that AhpA can scavenge peroxides and thus defend cells against peroxides, they did not clarify when during growth the cell produces AhpA. The results presented here show that the expression of ahpA is regulated in a manner distinct from that of the other peroxide-scavenging enzymes in B. subtilis. While the primary Ahp, AhpC, is expressed during exponential growth and stationary phase, these studies demonstrate that the expression of ahpA is dependent on the transition-state regulator AbrB and the sporulation and biofilm formation transcription factor Spo0A. Furthermore, these results show that ahpA is specifically expressed during biofilm formation, and not during sporulation or stationary phase, suggesting that derepression of ahpA by AbrB requires a signal other than those present upon entry into stationary phase. Despite this expression pattern, ahpA mutant strains still form and maintain robust biofilms, even in the presence of peroxides. Thus, the role of AhpA with regard to protecting cells within biofilms from environmental stresses is still uncertain. These studies highlight the need to further study the Ahp homologs to better understand how they differ from one another and the unique roles they may play in oxidative stress resistance.

Precise Identification of Genome-Wide Transcription Start Sites in Bacteria by 5'-Rapid Amplification of cDNA Ends (5'-RACE)

Transcription start sites are commonly used to locate promoter elements in bacterial genomes. TSS were previously studied one gene at a time, often through 5'-rapid amplification of cDNA ends (5'-RACE). This technique has now been adapted for high-throughput sequencing and can be used to precisely identify TSS in a genome-wide fashion for practically any bacterium, which greatly contributes to our understanding of gene regulatory networks in microorganisms.

Staphylococcus aureus Transcriptome Architecture: From Laboratory to Infection-Mimicking Conditions

Staphylococcus aureus is a major pathogen that colonizes about 20% of the human population. Intriguingly, this Gram-positive bacterium can survive and thrive under a wide range of different conditions, both inside and outside the human body. Here, we investigated the transcriptional adaptation of S. aureus HG001, a derivative of strain NCTC 8325, across experimental conditions ranging from optimal growth in vitro to intracellular growth in host cells. These data establish an extensive repertoire of transcription units and non-coding RNAs, a classification of 1412 promoters according to their dependence on the RNA polymerase sigma factors SigA or SigB, and allow identification of new potential targets for several known transcription factors. In particular, this study revealed a relatively low abundance of antisense RNAs in S. aureus, where they overlap only 6% of the coding genes, and only 19 antisense RNAs not co-transcribed with other genes were found. Promoter analysis and comparison with Bacillus subtilis links the small number of antisense RNAs to a less profound impact of alternative sigma factors in S. aureus. Furthermore, we revealed that Rho-dependent transcription termination suppresses pervasive antisense transcription, presumably originating from abundant spurious transcription initiation in this A+T-rich genome, which would otherwise affect expression of the overlapped genes. In summary, our study provides genome-wide information on transcriptional regulation and non-coding RNAs in S. aureus as well as new insights into the biological function of Rho and the implications of spurious transcription in bacteria.

The major human pathogen Staphylococcus aureus can survive under a wide range of conditions, both inside and outside the human body. The goal of this study was to determine how S. aureus adapts to such different conditions and, additionally, we wanted to identify general factors governing the staphylococcal transcriptome architecture. Therefore, we performed a precise analysis of all RNA transcripts of S. aureus across experimental conditions ranging from in vitro growth in different media to internalization by eukaryotic host cells. We systematically mapped all transcription units, annotated non-coding RNAs, and assigned promoters controlled by particular RNA polymerase sigma factors and transcription factors. By a comparison with data available for the related Gram-positive bacterium Bacillus subtilis, we made key observations concerning the abundance and origin of antisense RNAs. Intriguingly, these findings support the view that many antisense RNAs in a bacterium like B. subtilis could be byproducts of spurious promoter recognition by condition-specific alternative sigma factors. We also report that the transcription termination factor Rho prevents widespread antisense transcription, presumably caused by pervasive transcription initiation in the A+T-rich genome of S. aureus. Altogether our study presents new perspectives on the biological significance of antisense and pervasive transcription in bacteria.

Long 5' untranslated regions regulate the RNA stability of the deep-sea filamentous phage SW1

Virus production in the deep-sea environment has been found to be high, and viruses have been suggested to play significant roles in the overall functioning of this ecosystem. Nevertheless, little is known about these viruses, including the mechanisms that control their production, which makes them one of the least understood biological entities on Earth. Previously, we isolated the filamentous phage SW1, whose virus production and gene transcription were found to be active at low temperatures, from a deep-sea bacterium, Shewanella piezotolerans WP3. In this study, the operon structure of phage SW1 is presented, which shows two operons with exceptionally long 5′ and 3′ untranslated regions (UTRs). In addition, the 5′UTR was confirmed to significantly influence the RNA stability of the SW1 transcripts. Our study revealed novel regulation of the operon and led us to propose a unique regulatory mechanism for Inoviruses. This type of RNA-based regulation may represent a mechanism for significant viral production in the cold deep biosphere.

The hmsT 3' untranslated region mediates c-di-GMP metabolism and biofilm formation in Yersinia pestis

Yersinia pestis, the cause of plague, forms a biofilm in the proventriculus of its flea vector to enhance transmission. Biofilm formation in Y. pestis is regulated by the intracellular levels of cyclic diguanylate (c-di-GMP). In this study, we investigated the role of the 3' untranslated region (3'UTR) in hmsT mRNA, a transcript that encodes a diguanylate cyclase that stimulates biofilm formation in Y. pestis by synthesizing the second messenger c-di-GMP. Deletion of the 3'UTR increased the half-life of hmsT mRNA, thereby upregulating c-di-GMP levels and biofilm formation. Our findings indicate that multiple regulatory sequences might be present in the hmsT 3'UTR that function together to mediate mRNA turnover. We also found that polynucleotide phosphorylase is partially responsible for hmsT 3'UTR-mediated mRNA decay. In addition, the hmsT 3'UTR strongly repressed gene expression at 37°C and 26°C, but affected gene expression only slightly at 21°C. Our findings suggest that the 3'UTR might be involved in precise and rapid regulation of hmsT expression, allowing Y. pestis to fine-tune c-di-GMP synthesis and consequently regulate biofilm production to adapt to the changing host environment.

Insights into the Function of a Second, Nonclassical Ahp Peroxidase, AhpA, in Oxidative Stress Resistance in Bacillus subtilis

Organisms growing aerobically generate reactive oxygen-containing molecules, such as hydrogen peroxide (H2O2). These reactive oxygen molecules damage enzymes and DNA and may even cause cell death. In response, Bacillus subtilis produces at least nine potential peroxide-scavenging enzymes, two of which appear to be the primary enzymes responsible for detoxifying peroxides during vegetative growth: a catalase (encoded by katA) and an alkylhydroperoxide reductase (Ahp, encoded by ahpC). AhpC uses two redox-active cysteine residues to reduce peroxides to nontoxic molecules. A specialized thioredoxin-like protein, AhpF, is then required to restore oxidized AhpC back to its reduced state. Curiously, B. subtilis has two genes encoding Ahp: ahpC and ahpA. Although AhpC is well characterized, very little is known about AhpA. In fact, numerous bacterial species have multiple ahp genes; however, these additional Ahp proteins are generally uncharacterized. We seek to understand the role of AhpA in the bacterium's defense against toxic peroxide molecules in relation to the roles previously assigned to AhpC and catalase. Our results demonstrate that AhpA has catalytic activity similar to that of the primary enzyme, AhpC. Furthermore, our results suggest that a unique thioredoxin redox protein, AhpT, may reduce AhpA upon its oxidation by peroxides. However, unlike AhpC, which is expressed well during vegetative growth, our results suggest that AhpA is expressed primarily during postexponential growth.

IMPORTANCE

B. subtilis appears to produce nine enzymes designed to protect cells against peroxides; two belong to the Ahp class of peroxidases. These studies provide an initial characterization of one of these Ahp homologs and demonstrate that the two Ahp enzymes are not simply replicates of each other, suggesting that they instead are expressed at different times during growth of the cells. These results highlight the need to further study the Ahp homologs to better understand how they differ from one another and to identify their function, if any, in protection against oxidative stress. Through these studies, we may better understand why bacteria have multiple enzymes designed to scavenge peroxides and thus have a more accurate understanding of oxidative stress resistance.

Constitutive Stringent Response Restores Viability of Bacillus subtilis Lacking Structural Maintenance of Chromosome Protein

Bacillus subtilis mutants lacking the SMC-ScpAB complex are severely impaired for chromosome condensation and partitioning, DNA repair, and cells are not viable under standard laboratory conditions. We isolated suppressor mutations that restored the capacity of a smc deletion mutant (Δsmc) to grow under standard conditions. These suppressor mutations reduced chromosome segregation defects and abrogated hypersensitivity to gyrase inhibitors of Δsmc. Three suppressor mutations were mapped in genes involved in tRNA aminoacylation and maturation pathways. A transcriptomic survey of isolated suppressor mutations pointed to a potential link between suppression of Δsmc and induction of the stringent response. This link was confirmed by (p)ppGpp quantification which indicated a constitutive induction of the stringent response in multiple suppressor strains. Furthermore, sublethal concentrations of arginine hydroxamate (RHX), a potent inducer of stringent response, restored growth of Δsmc under non permissive conditions. We showed that production of (p)ppGpp alone was sufficient to suppress the thermosensitivity exhibited by the Δsmc mutant. Our findings shed new light on the coordination between chromosome dynamics mediated by SMC-ScpAB and other cellular processes during rapid bacterial growth.

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.

γ-PGA Hydrolases of Phage Origin in Bacillus subtilis and Other Microbial Genomes

Poly-γ-glutamate (γ-PGA) is an industrially interesting polymer secreted mainly by members of the class Bacilli which forms a shield able to protect bacteria from phagocytosis and phages. Few enzymes are known to degrade γ-PGA; among them is a phage-encoded γ-PGA hydrolase, PghP. The supposed role of PghP in phages is to ensure access to the surface of bacterial cells by dismantling the γ-PGA barrier. We identified four unannotated B. subtilis genes through similarity of their encoded products to PghP; in fact these genes reside in prophage elements of B. subtilis genome. The recombinant products of two of them demonstrate efficient polymer degradation, confirming that sequence similarity reflects functional homology. Genes encoding similar γ-PGA hydrolases were identified in phages specific for the order Bacillales and in numerous microbial genomes, not only belonging to that order. The distribution of the γ-PGA biosynthesis operon was also investigated with a bioinformatics approach; it was found that the list of organisms endowed with γ-PGA biosynthetic functions is larger than expected and includes several pathogenic species. Moreover in non-Bacillales bacteria the predicted γ-PGA hydrolase genes are preferentially found in species that do not have the genetic asset for polymer production. Our findings suggest that γ-PGA hydrolase genes might have spread across microbial genomes via horizontal exchanges rather than via phage infection. We hypothesize that, in natural habitats rich in γ-PGA supplied by producer organisms, the availability of hydrolases that release glutamate oligomers from γ-PGA might be a beneficial trait under positive selection.

Whole-genome mapping of 5' RNA ends in bacteria by tagged sequencing: a comprehensive view in Enterococcus faecalis

Enterococcus faecalis is the third cause of nosocomial infections. To obtain the first snapshot of transcriptional organizations in this bacterium, we used a modified RNA-seq approach enabling to discriminate primary from processed 5' RNA ends. We also validated our approach by confirming known features in Escherichia coli. We mapped 559 transcription start sites (TSSs) and 352 processing sites (PSSs) in E. faecalis. A blind motif search retrieved canonical features of SigA- and SigN-dependent promoters preceding transcription start sites mapped. We discovered 85 novel putative regulatory RNAs, small- and antisense RNAs, and 72 transcriptional antisense organizations. Presented data constitute a significant insight into bacterial RNA landscapes and a step toward the inference of regulatory processes at transcriptional and post-transcriptional levels in a comprehensive manner.

Tracking the Elusive Function of Bacillus subtilis Hfq

RNA-binding protein Hfq is a key component of the adaptive responses of many proteobacterial species including Escherichia coli, Salmonella enterica and Vibrio cholera. In these organisms, the importance of Hfq largely stems from its participation to regulatory mechanisms involving small non-coding RNAs. In contrast, the function of Hfq in Gram-positive bacteria has remained elusive and somewhat controversial. In the present study, we have further addressed this point by comparing growth phenotypes and transcription profiles between wild-type and an hfq deletion mutant of the model Gram-positive bacterium, Bacillus subtilis. The absence of Hfq had no significant consequences on growth rates under nearly two thousand metabolic conditions and chemical treatments. The only phenotypic difference was a survival defect of B. subtilis hfq mutant in rich medium in stationary phase. Transcriptomic analysis correlated this phenotype with a change in the levels of nearly one hundred transcripts. Albeit a significant fraction of these RNAs (36%) encoded sporulation-related functions, analyses in a strain unable to sporulate ruled out sporulation per se as the basis of the hfq mutant’s stationary phase fitness defect. When expressed in Salmonella, B. subtilis hfq complemented the sharp loss of viability of a degP hfq double mutant, attenuating the chronic σ^E-activated phenotype of this strain. However, B. subtilis hfq did not complement other regulatory deficiencies resulting from loss of Hfq-dependent small RNA activity in Salmonella indicating a limited functional overlap between Salmonella and B. subtilis Hfqs. Overall, this study confirmed that, despite structural similarities with other Hfq proteins, B. subtilis Hfq does not play a central role in post-transcriptional regulation but might have a more specialized function connected with stationary phase physiology. This would account for the high degree of conservation of Hfq proteins in all 17 B. subtilis strains whose genomes have been sequenced.

Dual-specificity anti-sigma factor reinforces control of cell-type specific gene expression in Bacillus subtilis

Gene expression during spore development in Bacillus subtilis is controlled by cell type-specific RNA polymerase sigma factors. σ^Fand σ^E control early stages of development in the forespore and the mother cell, respectively. When, at an intermediate stage in development, the mother cell engulfs the forespore, σ^F is replaced by σ^G and σ^E is replaced by σ^K. The anti-sigma factor CsfB is produced under the control of σ^F and binds to and inhibits the auto-regulatory σ^G, but not σ^F. A position in region 2.1, occupied by an asparagine in σ^G and by a glutamate in ο^F, is sufficient for CsfB discrimination of the two sigmas, and allows it to delay the early to late switch in forespore gene expression. We now show that following engulfment completion, csfB is switched on in the mother cell under the control of σ^K and that CsfB binds to and inhibits σ^E but not σ^K, possibly to facilitate the switch from early to late gene expression. We show that a position in region 2.3 occupied by a conserved asparagine in σ^E and by a conserved glutamate in σ^K suffices for discrimination by CsfB. We also show that CsfB prevents activation of σ^G in the mother cell and the premature σ^G-dependent activation of σ^K. Thus, CsfB establishes negative feedback loops that curtail the activity of σ^E and prevent the ectopic activation of σ^G in the mother cell. The capacity of CsfB to directly block σ^E activity may also explain how CsfB plays a role as one of the several mechanisms that prevent σ^E activation in the forespore. Thus the capacity of CsfB to differentiate between the highly similar σ^F/σ^G and σ^E/σ^K pairs allows it to rinforce the cell-type specificity of these sigma factors and the transition from early to late development in B. subtilis, and possibly in all sporeformers that encode a CsfB orthologue.

Precise temporal and cell-type specific regulation of gene expression is required for development of differentiated cells even in simple organisms. Endospore development by the bacterium Bacillus subtilis involves only two types of differentiated cells, a forespore that develops into the endospore, and a mother cell that nurtures the developing endospore. During development temporal and cell-type specific regulation of gene expression is controlled by transcription factors called sigma factors (σ). An anti-sigma factor known as CsfB binds to σ^G to prevent its premature activity in the forespore. We found that CsfB is also expressed in the mother cell where it blocks ectopic activity of σ^G, and blocks the activity σ^E to allow σ^K to take over control of gene expression during the final stages of development. Our finding that CsfB directly blocks σ^E activity also explains how CsfB plays a role in preventing ectopic activity of σ^E in the forespore. Remarkably, each of the major roles of CsfB, (i.e., control of ectopic σ^G and σ^E activities, and the temporal limitation of σ^E activity) is also accomplished by redundant regulatory processes. This redundancy reinforces control of key regulatory steps to insure reliability and stability of the developmental process.

Cis-encoded non-coding antisense RNAs in streptococci and other low GC Gram (+) bacterial pathogens

Due to recent advances of bioinformatics and high throughput sequencing technology, discovery of regulatory non-coding RNAs in bacteria has been increased to a great extent. Based on this bandwagon, many studies searching for trans-acting small non-coding RNAs in streptococci have been performed intensively, especially in the important human pathogen, group A and B streptococci. However, studies for cis-encoded non-coding antisense RNAs in streptococci have been scarce. A recent study shows antisense RNAs are involved in virulence gene regulation in group B streptococcus, S. agalactiae. This suggests antisense RNAs could have important roles in the pathogenesis of streptococcal pathogens. In this review, we describe recent discoveries of chromosomal cis-encoded antisense RNAs in streptococcal pathogens and other low GC Gram (+) bacteria to provide a guide for future studies.

Multiple ways to regulate translation initiation in bacteria: Mechanisms, regulatory circuits, dynamics

To adapt their metabolism rapidly and constantly in response to environmental variations, bacteria often target the translation initiation process, during which the ribosome assembles on the mRNA. Here, we review different mechanisms of regulation mediated by cis-acting elements, sRNAs and proteins, showing, when possible, their intimate connection with the translational apparatus. Indeed the ribosome itself could play a direct role in several regulatory mechanisms. Different features of the regulatory signals (sequences, structures and their positions on the mRNA) are contributing to the large variety of regulatory mechanisms. Ribosome heterogeneity, variation of individual cells responses and the spatial and temporal organization of the translation process add more layers of complexity. This hampers to define manageable set of rules for bacterial translation initiation control.

Genome-wide mapping of TnrA-binding sites provides new insights into the TnrA regulon in Bacillus subtilis

Under nitrogen limitation conditions, Bacillus subtilis induces a sophisticated network of adaptation responses. More precisely, the B. subtilis TnrA regulator represses or activates directly or indirectly the expression of a hundred genes in response to nitrogen availability. The global TnrA regulon have already been identified among which some directly TnrA-regulated genes have been characterized. However, a genome-wide mapping of in vivo TnrA-binding sites was still needed to clearly define the set of genes directly regulated by TnrA. Using chromatin immunoprecipitation coupled with hybridization to DNA tiling arrays (ChIP-on-chip), we now provide in vivo evidence that TnrA reproducibly binds to 42 regions on the chromosome. Further analysis with real-time in vivo transcriptional profiling, combined with results from previous reports, allowed us to define the TnrA primary regulon. We identified 35 promoter regions fulfilling three criteria necessary to be part of this primary regulon: (i) TnrA binding in ChIP-on-chip experiments and/or in previous in vitro studies; (ii) the presence of a TnrA box; (iii) TnrA-dependent expression regulation. In addition, the TnrA primary regulon delimitation allowed us to improve the TnrA box consensus. Finally, our results reveal new interconnections between the nitrogen regulatory network and other cellular processes.

Genome-wide analysis of phosphorylated PhoP binding to chromosomal DNA reveals several novel features of the PhoPR-mediated phosphate limitation response in Bacillus subtilis

The PhoPR two-component signal transduction system controls one of three responses activated by Bacillus subtilis to adapt to phosphate-limiting conditions (PHO response). The response involves the production of enzymes and transporters that scavenge for phosphate in the environment and assimilate it into the cell. However, in B. subtilis and some other Firmicutes bacteria, cell wall metabolism is also part of the PHO response due to the high phosphate content of the teichoic acids attached either to peptidoglycan (wall teichoic acid) or to the cytoplasmic membrane (lipoteichoic acid). Prompted by our observation that the phosphorylated WalR (WalR∼P) response regulator binds to more chromosomal loci than are revealed by transcriptome analysis, we established the PhoP∼P bindome in phosphate-limited cells. Here, we show that PhoP∼P binds to the chromosome at 25 loci: 12 are within the promoters of previously identified PhoPR regulon genes, while 13 are newly identified. We extend the role of PhoPR in cell wall metabolism showing that PhoP∼P binds to the promoters of four cell wall-associated operons (ggaAB, yqgS, wapA, and dacA), although none show PhoPR-dependent expression under the conditions of this study. We also show that positive autoregulation of phoPR expression and full induction of the PHO response upon phosphate limitation require PhoP∼P binding to the 3' end of the phoPR operon.

IMPORTANCE

The PhoPR two-component system controls one of three responses mounted by B. subtilis to adapt to phosphate limitation (PHO response). Here, establishment of the phosphorylated PhoP (PhoP∼P) bindome enhances our understanding of the PHO response in two important ways. First, PhoPR plays a more extensive role in adaptation to phosphate-limiting conditions than was deduced from transcriptome analyses. Among 13 newly identified binding sites, 4 are cell wall associated (ggaAB, yqgS, wapA, and dacA), revealing that PhoPR has an extended involvement in cell wall metabolism. Second, amplification of the PHO response must occur by a novel mechanism since positive autoregulation of phoPR expression requires PhoP∼P binding to the 3' end of the operon.

Genome-wide identification of Bacillus subtilis Zur-binding sites associated with a Zur box expands its known regulatory network

Background

The Bacillus subtilis Zur transcription factor recognizes a specific DNA motif, the Zur box, to repress expression of genes in response to zinc availability. Although several Zur-regulated genes are well characterized, a genome-wide mapping of Zur-binding sites is needed to define further the set of genes directly regulated by this protein.

Results

Using chromatin immunoprecipitation coupled with hybridization to DNA tiling arrays (ChIP-on-chip), we reported the identification of 80 inter- and intragenic chromosomal sites bound by Zur. Seven Zur-binding regions constitute the Zur primary regulon while 35 newly identified targets were associated with a predicted Zur box. Using transcriptional fusions an intragenic Zur box was showed to promote a full Zur-mediated repression when placed within a promoter region. In addition, intragenic Zur boxes appeared to mediate a transcriptional cis-repressive effect (4- to 9-fold) but the function of Zur at these sites remains unclear. Zur binding to intragenic Zur boxes could prime an intricate mechanisms of regulation of the transcription elongation, possibly with other transcriptional factors. However, the disruption of zinc homeostasis in Δzur cells likely affects many cellular processes masking direct Zur-dependent effects. Finally, most Zur-binding sites were located near or within genes responsive to disulfide stress. These findings expand the potential Zur regulon and reveal unknown interconnections between zinc and redox homeostasis regulatory networks.

Conclusions

Our findings considerably expand the potential Zur regulon, and reveal a new level of complexity in Zur binding to its targets via a Zur box motif and via a yet unknown mechanism that remains to be characterized.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-015-0345-4) contains supplementary material, which is available to authorized users.

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.

Expression of antisense small RNAs in response to stress in Pseudomonas aeruginosa

BACKGROUND

RNA sequencing technologies reveal that bacteria express RNA molecules other than mRNA, rRNA or tRNA. During the last years genome-wide bacterial transcriptomes have been shown to comprise intergenic RNA, antisense RNA, and untranslated regions, all capable of performing diverse regulatory functions.

RESULTS

In this study we used RNA-seq to identify 232 antisense RNAs (asRNAs) in the opportunistic pathogen Pseudomonas aeruginosa grown under 13 different conditions. The conditions studied include exponential and stationary growth as well as osmotic, oxidative and antibiotic stress. We found a significant overrepresentation of asRNAs that are transcribed opposite to genes involved in cell division and in cell wall, lipopolysaccharide (LPS), and capsule biosynthesis, most likely reflecting the conditions used in this study. A substantial number of asRNAs significantly changed their expression under osmotic, oxidative and antibiotic stress, suggesting that asRNAs may play regulatory roles during these conditions. We also made a comparison between the asRNAs detected in this study in P. aeruginosa PAO1 with the asRNAs detected in two previous studies in P. aeruginosa PA14, and found that the extent of overlap between the studies is very limited.

CONCLUSIONS

RNA-seq experiments are revealing hundreds of novel transcripts in all bacterial genomes investigated. The comparison between independent studies that used RNA-seq to detect novel asRNAs in P. aeruginosa shows that the overlap between the results reported is very narrow. It is necessary to address how reproducibility of these kind of studies should be reported in order to avoid misleading conclusions when comparing data generated by non-identical methods.

Global analysis of mRNA decay intermediates in Bacillus subtilis wild-type and polynucleotide phosphorylase-deletion strains

Messenger RNA decay in Bacillus subtilis is accomplished by a combination of exoribonucleases and endoribonucleases. Intermediates in the decay process have not been readily detectable, and previous studies on mRNA decay have used a handful of highly expressed transcripts as models. Here, we use RNA-Seq analysis to probe mRNA turnover globally. A significant fraction of messages showed differential accumulation of RNA fragments that mapped near the 5' or 3' end of the coding sequence, consistent with initiation of decay from either the 5' end or from an internal cleavage site. Patterns of mRNA decay in the wild type were compared with patterns in a mutant strain lacking polynucleotide phosphorylase (PNPase), which is considered the major 3' exonuclease activity in mRNA decay and which is one of four known 3' exonucleases in B. subtilis. The results showed a striking dependence on PNPase for mRNA turnover in many cases, suggesting specificity in the ability of 3' exonucleases to degrade from 3'-hydroxyl termini. RNA-Seq data demonstrated a sharp decrease in expression of Sigma D in the PNPase-deletion strain. Reduction in sigD regulon expression explained the chain growth phenotype of the PNPase mutant and also predicted a defect in swarming motility.

Impact of Hfq on the Bacillus subtilis transcriptome

The RNA chaperone Hfq acts as a central player in post-transcriptional gene regulation in several Gram-negative Bacteria, whereas comparatively little is known about its role in Gram-positive Bacteria. Here, we studied the function of Hfq in Bacillus subtilis, and show that it confers a survival advantage. A comparative transcriptome analysis revealed mRNAs with a differential abundance that are governed by the ResD-ResE system required for aerobic and anaerobic respiration. Expression of resD was found to be up-regulated in the hfq- strain. Furthermore, several genes of the GerE and ComK regulons were de-regulated in the hfq- background. Surprisingly, only six out of >100 known and predicted small RNAs (sRNAs) showed altered abundance in the absence of Hfq. Moreover, Hfq positively affected the transcript abundance of genes encoding type I toxin-antitoxin systems. Taken the moderate effect on sRNA levels and mRNAs together, it seems rather unlikely that Hfq plays a central role in RNA transactions in Bacillus subtilis.