APERO: a genome-wide approach for identifying bacterial small RNAs from RNA-Seq data.

Small non-coding RNAs (sRNAs) regulate numerous cellular processes in all domains of life. Several approaches have been developed to identify them from RNA-seq data, which are efficient for eukaryotic sRNAs but remain inaccurate for the longer and highly structured bacterial sRNAs. We present APERO, a new algorithm to detect small transcripts from paired-end bacterial RNA-seq data. In contrast to previous approaches that start from the read coverage distribution, APERO analyzes boundaries of individual sequenced fragments to infer the 5' and 3' ends of all transcripts. Since sRNAs are about the same size as individual fragments (50-350 nucleotides), this algorithm provides a significantly higher accuracy and robustness, e.g., with respect to spontaneous internal breaking sites. To demonstrate this improvement, we develop a comparative assessment on datasets from Escherichia coli and Salmonella enterica, based on experimentally validated sRNAs. We also identify the small transcript repertoire of Dickeya dadantii including putative intergenic RNAs, 5' UTR or 3' UTR-derived RNA products and antisense RNAs. Comparisons to annotations as well as RACE-PCR experimental data confirm the precision of the detected transcripts. Altogether, APERO outperforms all existing methods in terms of sRNA detection and boundary precision, which is crucial for comprehensive genome annotations. It is freely available as an open source R package on https://github.com/Simon-Leonard/APERO.

RhoTermPredict: an algorithm for predicting Rho-dependent transcription terminators based on Escherichia coli, Bacillus subtilis and Salmonella enterica databases.

BACKGROUND: In bacterial genomes, there are two mechanisms to terminate the DNA transcription: the "intrinsic" or Rho-independent termination and the Rho-dependent termination. Intrinsic terminators are characterized by a RNA hairpin followed by a run of 6-8 U residues relatively easy to identify using one of the numerous available prediction programs. In contrast, Rho-dependent termination is mediated by the Rho protein factor that, firstly, binds to ribosome-free mRNA in a site characterized by a C > G content and then reaches the RNA polymerase to induce its release. Conversely on intrinsic terminators, the computational prediction of Rho-dependent terminators in prokaryotes is a very difficult problem because the sequence features required for the function of Rho are complex and poorly defined. This is the reason why it still does not exist an exhaustive Rho-dependent terminators prediction program. RESULTS: In this study we introduce RhoTermPredict, the first published algorithm for an exhaustive Rho-dependent terminators prediction in bacterial genomes. RhoTermPredict identifies these elements based on a previously proposed consensus motif common to all Rho-dependent transcription terminators. It essentially searches for a 78 nt long RUT site characterized by a C > G content and with regularly spaced C residues, followed by a putative pause site for the RNA polymerase. We tested RhoTermPredict performances by using available genomic and transcriptomic data of the microorganism Escherichia coli K-12, both in limited-length sequences and in the whole-genome, and available genomic sequences from Bacillus subtilis 168 and Salmonella enterica LT2 genomes. We also estimated the overlap between the predictions of RhoTermPredict and those obtained by the predictor of intrinsic terminators ARNold webtool. Our results demonstrated that RhoTermPredict is a very performing algorithm both for limited-length sequences (F1-score obtained about 0.7) and for a genome-wide analysis. Furthermore the degree of overlap with ARNold predictions was very low. CONCLUSIONS: Our analysis shows that RhoTermPredict is a powerful tool for Rho-dependent terminators search in the three analyzed genomes and could fill this gap in computational genomics. We conclude that RhoTermPredict could be used in combination with an intrinsic terminators predictor in order to predict all the transcription terminators in bacterial genomes.

A synthetic growth switch based on controlled expression of RNA polymerase.

The ability to control growth is essential for fundamental studies of bacterial physiology and biotechnological applications. We have engineered an Escherichia coli strain in which the transcription of a key component of the gene expression machinery, RNA polymerase, is under the control of an inducible promoter. By changing the inducer concentration in the medium, we can adjust the RNA polymerase concentration and thereby switch bacterial growth between zero and the maximal growth rate supported by the medium. We show that our synthetic growth switch functions in a medium-independent and reversible way, and we provide evidence that the switching phenotype arises from the ultrasensitive response of the growth rate to the concentration of RNA polymerase. We present an application of the growth switch in which both the wild-type E. coli strain and our modified strain are endowed with the capacity to produce glycerol when growing on glucose. Cells in which growth has been switched off continue to be metabolically active and harness the energy gain to produce glycerol at a twofold higher yield than in cells with natural control of RNA polymerase expression. Remarkably, without any further optimization, the improved yield is close to the theoretical maximum computed from a flux balance model of E. coli metabolism. The proposed synthetic growth switch is a promising tool for gaining a better understanding of bacterial physiology and for applications in synthetic biology and biotechnology.

Characterization of the Escherichia coli σ(S) core regulon by Chromatin Immunoprecipitation-sequencing (ChIP-seq) analysis.

In bacteria, selective promoter recognition by RNA polymerase is achieved by its association with σ factors, accessory subunits able to direct RNA polymerase "core enzyme" (E) to different promoter sequences. Using Chromatin Immunoprecipitation-sequencing (ChIP-seq), we searched for promoters bound by the σ(S)-associated RNA polymerase form (Eσ(S)) during transition from exponential to stationary phase. We identified 63 binding sites for Eσ(S) overlapping known or putative promoters, often located upstream of genes (encoding either ORFs or non-coding RNAs) showing at least some degree of dependence on the σ(S)-encoding rpoS gene. Eσ(S) binding did not always correlate with an increase in transcription level, suggesting that, at some σ(S)-dependent promoters, Eσ(S) might remain poised in a pre-initiation state upon binding. A large fraction of Eσ(S)-binding sites corresponded to promoters recognized by RNA polymerase associated with σ(70) or other σ factors, suggesting a considerable overlap in promoter recognition between different forms of RNA polymerase. In particular, Eσ(S) appears to contribute significantly to transcription of genes encoding proteins involved in LPS biosynthesis and in cell surface composition. Finally, our results highlight a direct role of Eσ(S) in the regulation of non coding RNAs, such as OmrA/B, RyeA/B and SibC.

sigmaS, a major player in the response to environmental stresses in Escherichia coli: role, regulation and mechanisms of promoter recognition.

Bacterial cells often face hostile environmental conditions, to which they adapt by activation of stress responses. In Escherichia coli, environmental stresses resulting in significant reduction in growth rate stimulate the expression of the rpoS gene, encoding the alternative σ factor σ(S). The σ(S) protein associates with RNA polymerase, and through transcription of genes belonging to the rpoS regulon allows the activation of a 'general stress response', which protects the bacterial cell from harmful environmental conditions. Each step of this process is finely tuned in order to cater to the needs of the bacterial cell: in particular, selective promoter recognition by σ(S) is achieved through small deviations from a common consensus DNA sequence for both σ(S) and the housekeeping σ(70). Recognition of specific DNA elements by σ(S) is integrated with the effects of environmental signals and the interaction with regulatory proteins, in what represents a fascinating example of multifactorial regulation of gene expression. In this report, we discuss the function of the rpoS gene in the general stress response, and review the current knowledge on regulation of rpoS expression and on promoter recognition by σ(S).

Repression of flagellar genes in exponential phase by CsgD and CpxR, two crucial modulators of Escherichia coli biofilm formation.

Escherichia coli adapts its lifestyle to the variations of environmental growth conditions, swapping between swimming motility or biofilm formation. The stationary-phase sigma factor RpoS is an important regulator of this switch, since it stimulates adhesion and represses flagellar biosynthesis. By measuring the dynamics of gene expression, we show that RpoS inhibits the transcription of the flagellar sigma factor, FliA, in exponential growth phase. RpoS also partially controls the expression of CsgD and CpxR, two transcription factors important for bacterial adhesion. We demonstrate that these two regulators repress the transcription of fliA, flgM, and tar and that this regulation is dependent on the growth medium. CsgD binds to the flgM and fliA promoters around their -10 promoter element, strongly suggesting direct repression. We show that CsgD and CpxR also affect the expression of other known modulators of cell motility. We propose an updated structure of the regulatory network controlling the choice between adhesion and motility.

Expression dynamics of RpoS/Crl-dependent genes in Escherichia coli.

The alternative sigma factor RpoS is a central regulator of the stress response in many Proteobacteria, acting both during exponential growth and in stationary phase. The small protein Crl increases the interaction between RpoS and RNA polymerase and thereby activates certain RpoS-dependent promoters. However, the growth-phase dependence of the interaction of Crl with different forms of polymerase remains unknown. We use 41 GFP transcriptional fusions to study the dynamics of gene regulation by RpoS and Crl during growth transition from exponential to stationary phase in Escherichia coli. We confirm that RpoS can regulate gene expression in exponential phase, both positively and negatively. Crl slightly stimulates transcription by RpoS in exponential phase and controls a subset of RpoS-dependent genes in stationary phase. Growth temperature strongly affects induction of specific promoters by RpoS, whereas its impact on gene regulation by Crl is much less significant. In addition, we identify five new genes regulated by Crl (ada, cbpA, glgS, sodC and flgM) and demonstrate that Crl improves promoter binding and opening by RpoS-containing RNA polymerase at the hdeA promoter. Our study also shows that Crl is a cognate enhancer of RpoS activity under different growth conditions, since its deletion has no effect on genes transcribed by other sigma factors.

Substitutions in region 2.4 of sigma70 allow recognition of the sigmaS-dependent aidB promoter.

The strict dependence of transcription from the aidB promoter (PaidB) on the Esigma(S) form of RNA polymerase is because of the presence of a C nucleotide as the first residue of the -10 promoter sequence (-12C), which does not allow an open complex formation by Esigma(70). In this report, sigma(70) mutants carrying either the Q437H or the T440I single amino acid substitutions, which allow -12C recognition by sigma(70), were tested for their ability to carry out transcription from PaidB. The Gln-437 and Thr-440 residues are located in region 2.4 of sigma(70) and correspond to Gln-152 and Glu-155 in sigma(S). Interestingly, the Q437H mutant of sigma(70), but not T440I, was able to promote an open complex formation and to initiate transcription at PaidB. In contrast to T440I, a T440E mutant was proficient in carrying out transcription from PaidB. No sigma(70) mutant displayed significantly increased interaction with a PaidB mutant in which the -12C was substituted by a T (PaidB((C12T))), which is also efficiently recognized by wild type sigma(70). The effect of the T440E mutation suggests that the corresponding Glu-155 residue in sigma(S) might be involved in -12C recognition. However, substitution to alanine of the Glu-155 residue, as well as of Gln-152, in the sigma(S) protein did not significantly affect Esigma(S) interaction with PaidB. Our results reiterate the importance of the -12C residue for sigma(S)-specific promoter recognition and strongly suggest that interaction with the -10 sequence and open complex formation are carried out by different determinants in the two sigma factors.

SigmaS-dependent gene expression at the onset of stationary phase in Escherichia coli: function of sigmaS-dependent genes and identification of their promoter sequences.

The sigma(S) subunit of RNA polymerase, the product of the rpoS gene, controls the expression of genes responding to starvation and cellular stresses. Using gene array technology, we investigated rpoS-dependent expression at the onset of stationary phase in Escherichia coli grown in rich medium. Forty-one genes were expressed at significantly lower levels in an rpoS mutant derived from the MG1655 strain; for 10 of these, we also confirmed rpoS and stationary-phase dependence by reverse transcription-PCR. Only seven genes (dps, osmE, osmY, sodC, rpsV, wrbA, and yahO) had previously been recognized as rpoS dependent. Several newly identified rpoS-dependent genes are involved in the uptake and metabolism of amino acids, sugars, and iron. Indeed, the rpoS mutant strain shows severely impaired growth on some sugars such as fructose and N-acetylglucosamine. The rpoS gene controls the production of indole, which acts as a signal molecule in stationary-phase cells, via regulation of the tnaA-encoded tryptophanase enzyme. Genes involved in protein biosynthesis, encoding the ribosome-associated protein RpsV (sra) and the initiation factor IF-1 (infA), were also induced in an rpoS-dependent fashion. Using primer extension, we determined the promoter sequences of a selection of rpoS-regulated genes representative of different functional classes. Significant fractions of these promoters carry sequence features specific for Esigma(S) recognition of the -10 region, such as cytosines at positions -13 (70%) and -12 (30%) as well as a TG motif located upstream of the -10 region (50%), thus supporting the TGN(0-2)C(C/T)ATA(C/A)T consensus sequence recently proposed for sigma(S).

Nucleotides from -16 to -12 determine specific promoter recognition by bacterial sigmaS-RNA polymerase.

The alternative sigma factor sigmaS, mainly active in stationary phase of growth, recognizes in vitro a -10 promoter sequence almost identical to the one for the main sigma factor, sigma70, thus raising the problem of how specific promoter recognition by sigmaS-RNA polymerase (EsigmaS) is achieved in vivo. We investigated the promoter features involved in selective recognition by EsigmaS at the strictly sigmaS-dependent aidB promoter. We show that the presence of a C nucleotide as first residue of the aidB -10 sequence (-12C), instead of the T nucleotide canonical for sigma70-dependent promoters, is the major determinant for selective recognition by EsigmaS. The presence of the -12C does not allow formation of an open complex fully proficient in transcription initiation by Esigma70. The role of -12C as specific determinant for promoter recognition by EsigmaS was confirmed by sequence analysis of known EsigmaS-dependent promoters as well as site-directed mutagenesis at the promoters of the csgB and sprE genes. We propose that EsigmaS, unlike Esigma70, can recognize both C and T as the first nucleotide in the -10 sequence. Additional promoter features such as the presence of a C nucleotide at position -13, contributing to open complex formation by EsigmaS, and a TG motif found at the unusual -16/-15 location, possibly contributing to initial binding to the promoter, also represent important factors for sigmaS-dependent transcription. We propose a new sequence, TG(N)0-2CCATA(c/a)T, as consensus -10 sequence for promoters exclusively recognized by EsigmaS.

Mechanism of specific recognition of the aidB promoter by sigma(S)-RNA polymerase.

Transcription of the Escherichia coli aidB gene is controlled by an Esigma(S)-dependent promoter (PaidB) and is poorly transcribed by the Esigma(70) form of RNA polymerase in the absence of additional factors. In this report, we investigate the interaction between PaidB and either the Esigma(70) or the Esigma(S) forms of RNA polymerase in vitro. We show that although Esigma(70) can bind the aidB promoter, its interaction with the promoter results in the formation of an open complex inefficient in transcription initiation and sensitive to heparin challenge. Deletion of the C residue at position -13 of PaidB (Delta-13C) slightly impaired transcription initiation by Esigma(S), consistent with the role of -13C as a specific feature of Esigma(S)-dependent promoters. However, Esigma(S) could still bind and initiate transcription from the Delta-13C mutant aidB promoter more efficiently than Esigma(70), suggesting that sequence elements other than the -13C play an important role in specific promoter recognition by Esigma(S).