The architecture and ppGpp-dependent expression of the primary transcriptome of Salmonella Typhimurium during invasion gene expression
- Vinoy K Ramachandran†1, 3,
- Neil Shearer†1, 3,
- Jobin J Jacob1, 3,
- Cynthia M Sharma2 and
- Arthur Thompson1, 3Email author
© Ramachandran et al; licensee BioMed Central Ltd. 2012
Received: 24 October 2011
Accepted: 17 January 2012
Published: 17 January 2012
Invasion of intestinal epithelial cells by Salmonella enterica serovar Typhimurium (S. Typhimurium) requires expression of the extracellular virulence gene expression programme (STEX), activation of which is dependent on the signalling molecule guanosine tetraphosphate (ppGpp). Recently, next-generation transcriptomics (RNA-seq) has revealed the unexpected complexity of bacterial transcriptomes and in this report we use differential RNA sequencing (dRNA-seq) to define the high-resolution transcriptomic architecture of wild-type S. Typhimurium and a ppGpp null strain under growth conditions which model STEX. In doing so we show that ppGpp plays a much wider role in regulating the S. Typhimurium STEX primary transcriptome than previously recognised.
Here we report the precise mapping of transcriptional start sites (TSSs) for 78% of the S. Typhimurium open reading frames (ORFs). The TSS mapping enabled a genome-wide promoter analysis resulting in the prediction of 169 alternative sigma factor binding sites, and the prediction of the structure of 625 operons. We also report the discovery of 55 new candidate small RNAs (sRNAs) and 302 candidate antisense RNAs (asRNAs). We discovered 32 ppGpp-dependent alternative TSSs and determined the extent and level of ppGpp-dependent coding and non-coding transcription. We found that 34% and 20% of coding and non-coding RNA transcription respectively was ppGpp-dependent under these growth conditions, adding a further dimension to the role of this remarkable small regulatory molecule in enabling rapid adaptation to the infective environment.
The transcriptional architecture of S. Typhimurium and finer definition of the key role ppGpp plays in regulating Salmonella coding and non-coding transcription should promote the understanding of gene regulation in this important food borne pathogen and act as a resource for future research.
Pathogenic strains of Salmonella continue to pose an unacceptable worldwide threat to the health of humans and livestock. Infection of humans with S. Typhimurium results in a debilitating case of severe gastroenteritis that may result in death in immunocompromised individuals. There are about 1.3 billion cases of non-typhoidal salmonellosis worldwide each year and it is estimated that there are 17 million cases and over 500,000 deaths each year caused by typhoid fever . In the current study we focus on S. Typhimurium, which once ingested via contaminated food or water, invades human gut epithelial cells resulting in bloody diarrhoea. S. Typhimurium is able to invade intestinal epithelial cells due to the expression of a horizontally acquired set of virulence genes (Salmonella Pathogenicity Island 1; SPI1), which encode a type 3 secretion system (T3SS) . In the case of murine infection, S. Typhimurium can become systemic and cause a typhoid-like fever due to its ability to replicate and survive within macrophages; this is achieved by the expression of a second T3SS encoded by genes within SPI2 . The complex expression patterns of SPI1 and SPI2 during infection led us and others to develop the concept of the Salmonella extracellular (STEX) and intracellular (STIN) virulence gene expression programmes, [4, 5] Successful host invasion and colonisation requires expression of the STEX virulence gene programme followed by expression of the STIN programme (characterised by SPI1 and SPI2 expression respectively) .
The environmentally-dependent expression of nearly all of the STEX and STIN genes in S. Typhimurium is mediated by the bacterial alarmone, guanosine tetraphosphate (ppGpp) . In Salmonella and all beta- and gammaproteobacteria, ppGpp is produced by the activity of two enzymes, RelA and SpoT [for review see ]. Whilst RelA is only able to synthesise ppGpp, SpoT contains both synthetase and hydrolase activities. In most other bacteria RelA and SpoT are combined into a single enzyme referred to as Rel or RSH (R elA S poT h omologue) . Previous work implicates SpoT rather than RelA in Salmonella pathogenicity since an S. Typhimurium ΔrelA strain is almost fully virulent in BALB/c mouse infection studies, whereas a ΔrelA ΔspoT strain is severely attenuated . It has also been shown that ppGpp plays a key role in coupling virulence to metabolic status in several other pathogenic bacteria including Mycobacterium tuberculosis[10, 11], Listeria monocytogenes, Legionella pneumophilia[13, 14], Vibrio cholera and Pseudomonas aeruginosa. A complete understanding of the pathways and mechanisms by which ppGpp mediates bacterial virulence may suggest targets for antimicrobial therapies .
Guanosine tetraphosphate appears to exert most of its physiological effects by direct or indirect transcriptional control of target genes and binds near the active centre of RNA polymerase (RNAP) to modulate its activity, resulting in the direct repression of stable RNA operons (for review see . This is suggested to increase the availability of RNAP for activation of genes required for survival under various stressful conditions . One mechanism by which this occurs is via sigma factor competition, whereby ppGpp reduces the affinity of core RNAP for σ70 resulting in an increase in the availability of RNAP to bind alternative stress-response sigma factors . Although this model suggests an indirect mechanism for ppGpp activation of gene expression, direct activation has been observed at some promoters [20, 21]. The effect of ppGpp on transcription can also be potentiated by the RNAP accessory protein DksA which may help to stabilise the binding of ppGpp to RNAP .
Recently global transcriptome analysis using high-density tiling arrays and high throughput RNA sequencing (RNA-seq) has revealed an unexpected complexity of bacterial and archaeal transcriptomes [23–26]. A major advance in this area has been the development of differential RNA sequencing (dRNA-seq) which allows global and unambiguous mapping of transcription start sites (TSSs) [24, 27]. In this study we utilise dRNA-seq technology to define the primary transcriptomes of wild-type S. Typhimurium and an isogenic ΔrelA ΔspoT mutant, in order to define the extent of ppGpp-dependent expression. We identified primary TSSs for 78% of the annotated S. Typhimurium genes as well as ppGpp-dependent and independent alternative TSSs. We confirm the expression of known and predicted sRNAs , identify new candidate sRNAs, and report the discovery of 302 candidate antisense transcripts for the entire S. Typhimurium genome. Our data provides further insights into the regulatory roles of ppGpp, confirming and extending a previously reported link to global regulation of non-coding RNAs . The high resolution transcriptomic datasets presented here should facilitate future research on transcriptional and post-transcriptional regulation of virulence and other adaptive mechanisms within Salmonella.
Identification of transcriptional start sites
The nucleotide position of TSSs were identified from a dRNA-seq analysis of RNA samples isolated from the S. Typhimurium wild type strain (SL1344) and an isogenic ΔrelA ΔspoT strain grown to early stationary phase (SPI1 inducing conditions). The dRNA-seq analysis was performed according to Sharma et al. For each strain two cDNA libraries were prepared from the same total RNA sample. One library, referred to as (+), was enriched for primary transcripts by treating with terminator exonuclease (see Materials and Methods) and the second library, referred to as (-) or non-enriched, was untreated and contained both primary and processed transcripts. Following sequencing of the cDNA libraries on Roche-454 and Illumina-Solexa platforms the reads were mapped onto the SL1344 genome (including the endogenous SLP1-3 plasmids) and the number of reads mapping to each nucleotide position were visualised using the integrated genome browser (IGB; http://bioviz.org/igb/). Elevated read numbers at the 5' end of transcripts in the (+) library relative to the (-) library were identified as an increased presence of transcripts with 5' PPP end sequences compared to 5' P end sequences as described previously . Although 454 sequencing provided longer read lengths, the read numbers of the Solexa dataset were considerably higher and were used primarily for identification of TSSs (see additional file 1: Table S1 for sequencing and mapping statistics).
Summary of TSS types and categorisations.
Known & predicted sRNAs
New candidate sRNAs
Other candidate ncRNAs
The direct visualisation of transcribed genomic loci and unambiguous mapping of primary TSSs enabled optimisation of the SL1344 genome annotation. Where transcription was observed in regions where no gene was previously annotated the Artemis genome browser and annotation tool (http://www.sanger.ac.uk/resources/software/artemis/) was used to search for potential ORFs possessing upstream Shine-Dalgarno sequences and putative homologues were identified using BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi). We applied a similar procedure to re-annotate start codons where the TSS was found to be downstream of the previously annotated start. These procedures resulted in the re-annotation of 60 start codons (additional file 2: Table S9), the identification of 23 potential new ORFs (additional file 2: Table S7), and the re-designation of 2 ORFs previously annotated in different reading frames. Five of the new ORFs (ibs123 and ldrAB) were predicted to be small toxic peptides of the Type 1 toxin-antitoxin systems found in E. coli[31, 32].
Promoter analysis of transcriptional start sites
The dRNA-seq derived promoter database was also exploited to perform a MEME analysis of all of the ORF promoters to identify candidate and known targets for the alternative sigma factors σ24, σ28, σ32, σ38 & σ54. Firstly, a position specific probability matrix (PSPM) was derived by MEME analysis for each of E. coli alternative sigma factors from a promoter dataset from regulonDB http://regulondb.ccg.unam.mx/. The PSPM was then interrogated with FIMO (Find Individual Motif Occurrences) with a default p-value cut-off of 0.00001 to identify sigma factor specific promoters from the dRNA-seq promoter database. FIMO identified candidate binding sites for σ24, σ28, σ32, σ38 & σ54 factors at 20, 21, 34, 92 and 2 promoters respectively (additional file 2: Table S1). Our analysis was able to identify the majority of the genes previously shown to be dependent on sigma factors (6, 8, 13 and 9 genes for σ24, σ28, σ32 and σ38 respectively from S. Typhimurium or E. coli)[33–36]. The published genes shown to be dependent on alternative sigma factors are indicated in additional file 2: Table S1. A MEME analysis was performed on the remaining 592 promoters that did not contain identifiable potential sigma factor motifs. The analysis identified a conserved region for 264 promoters which was similar to the -10 motif, but contained a longer (5 nt) region between the first conserved T of the -10 motif (at -6) and the first nucleotide of the TSS (Figure 3C). A functional categorisation of the 264 promoters revealed that a strikingly high percentage (35%) belonged to genes involved in nucleoside and nucleotide interconversions (see additional file 3: Figure S4).
5' leader regions and leaderless mRNAs
Canonical bacterial mRNAs contain a 5' untranslated region (5' UTR) upstream of the initiation codon. As a minimal requirement, this region contains the Shine-Dalgarno (SD) ribosome binding site (RBS), and may contain additional sequence motifs required for efficient ribosome binding [37, 38]. In addition, many mRNAs include longer 5' leader regions which may possess regulatory functions. The structure and sequence of 5' leader regions can affect gene expression by modulating the synthesis of full length mRNAs or via regulation of post-transcriptional processing (e.g. synthesis of leader peptides, formation of secondary structures including riboswitches, or binding of proteins or regulatory sRNAs) . In addition to the majority of mRNAs, a few bacterial transcripts lack upstream UTRs and are termed "leaderless transcripts". In such cases transcription starts at, or up to 6 nucleotides upstream of the "A" residue of an AUG start codon [40, 41] and the resulting transcripts lack a SD RBS. Published examples include the cI repressor of bacteriophage λ , and the tetR gene of transposon Tn1721. Both the cl and tetR genes encode relatively low abundance regulatory proteins which is consistent with the reduced translation of genes lacking ribosomal binding sites . Interestingly Sullivan et al recently showed that the leaderless mRNA transcript of a regulatory gene, acuR, (involved in regulation of an operon encoding products involved in dimethylsulfoniopropionate catabolism in Rhodobacter sphaeroides), is transcribed at least as efficiently as downstream genes, but is translated at far lower levels, thus providing an elegant mechanism for differential control of operon-encoded protein levels. Although previously thought to be rare, leaderless genes are now known to be fairly common in prokaryotes and in the archaea, where as many as 69% of protein coding transcripts are leaderless [24, 38, 46]. Our dRNA-seq analysis identified 16 completely leaderless mRNAs where transcription started precisely at the A residue of the AUG start codon. A further 17 genes contained a leader of between 1 and 6 nt in length but lacked a SD sequence (additional file 2: Table S2). Functional analysis of these 33 leaderless genes showed that 6 encoded transcriptional regulators (including a TetR homologue), expected to be expressed at relatively low levels, 7 encoded membrane proteins which are also often required at low levels, and 4 were related to pathogenicity functions (additional file 2: Table S2).
Operons were identified during manual inspection of the strand-specific dRNA-seq data (additional file 2: Table S1). We found no difference in the operon structures determined from the wild-type and ΔrelA ΔspoT genomes (data not shown). The majority of 5' ends of operons were assigned from the primary TSS at the start of the first gene or from TSS's located within the first gene of the operon (where present). A further criterion for identification was the presence of 3' UTRs located at the ends of operons. Our inspection identified 1538 monocistronic transcripts and 625 polycistronic operons resulting in a mean of 1.65 genes per operon. We compared our operon map to operons predicted using DOOR (Database of prOkaryotic OpeRons; [49, 50]; see additional file 2: Table S1). DOOR predicts operons based on a comparison of 675 prokaryotic organisms and accuracy can reach 90.2% and 93.7% for the B. subtilis and E. coli genomes respectively (http://csbl1.bmb.uga.edu/OperonDB/DOOR.php). DOOR analysis of the S. Typhimurium SL1344 genome predicted 955 operons. Our comparison of dRNA-seq to DOOR predicted operons identified 60% (372) with an exact match and 24 (4%) new operons, not predicted by DOOR. We found 36% (229) of the operons identified from our dRNA-seq data were either extended or shortened (by one or two genes) compared to the DOOR predicted operons. In these cases we found that the dRNA-seq data identified TSSs that were located within operons predicted by DOOR (e.g. hypA and hypD contain internal TSSs within the hyp operon which encodes hydrogenase maturation factors; additional file 2: Table S1 ). Since the DOOR algorithm does not take into consideration TSS information, we suggest that our dRNA-seq identified operons are likely to be more accurate than the DOOR predicted operon structures.
dRNA-seq identification of sRNA expression in S. Typhimurium
Manual inspection of the dRNA-seq transcriptome of the wild-type and ΔrelA ΔspoT strains identified a total of 83 predicted and known sRNAs and we discovered a further 55 new candidate sRNAs (Table 1). We validated expression of 3 known (RprA, InvR, GcvB), 1 predicted (STnc1020) and 6 new candidate sRNAs (SLnc0011, SLnc0027, SLP1_ncRNA3, SLP1_ncRNA6, SLP2_ncRNA12, SLP2_ncRNA1) using Northern blotting (see additional file 2: Figure S5, additional file 1: Tables S3 and S4), and verified the presence of conserved -10 and -35 regions for the known and predicted sRNAs (Figure 3D). In order to further validate the S. Typhimurium new candidate and predicted sRNAs we determined whether they were conserved within the recently published S. Typhi ncRNA transcriptome . Of the 25 newly identified S. Typhi sRNAs, 20 sRNA homologues were found within the SL1344 genome. Of the 20 homologues, 2 new candidate sRNAs, 1 predicted sRNA and 1 asRNA (see following section) were expressed in SL1344 under the growth conditions used in this study (SLnc1039, SLnc1005, STnc560 and SLaRNA0247 respectively; additional file 2: Tables S3, S4 and S5). Interestingly 19% of the known, predicted and new candidate sRNAs have secondary TSSs, indicating that they may be subject to differential regulation (Table 1). Finally, we predicted intrinsic transcriptional terminators for 32 new candidate sRNAs and using TargetRNA software (http://snowwhite.wellesley.edu/targetRNA/) we were also able to predict potential targets (additional file 2: Table S4) .
Extensive antisense transcription in S. Typhimurium
Antisense RNAs (asRNAs) have been shown to be particularly abundant in eukaryotes, and recently a large proportion of the primary TSSs have been shown to be antisense to ORFs in E. coli and H. pylori, suggesting that asRNAs have a widespread regulatory function in bacteria [24, 25, 54, 55]. The dRNA-seq analysis detected 302 potential asRNAs in S. Typhimurium which were located directly opposite to coding regions of chromosomal genes (Table 1, additional file 2: Table S5). We also annotated ncRNAs which were within or close to the 3' and 5' UTRs of genes but which could not unambiguously be identified as asRNAs according to our strict definition (Table 1, additional file 2: Table S5). Finally, we found 94 ncRNAs which were located in intergenic regions (i.e. greater than 250 nt from the 3' or 5' ends of a gene) where one or both of the flanking genes were located on the opposing or same strand of the ncRNA (additional file 2: Table S5). For validating the presence of selected candidate antisense and ncRNAs in RNA samples two methods were used, Northern blotting and the more sensitive adapter assisted PCR. We verified the presence of 2 candidate asRNAs, SLasRNA0330 and SLasRNA0183 (additional file 3: Figs. S5 and S6). SLasRNA0330 is opposite to the sipA ORF, which encodes a SPI1 effector protein and SLasRNA0183 is opposite to the ycfQ ORF which encodes a putative transcriptional repressor (additional file 2: Table S1). The presence of 5 candidate ncRNAs were also verified (the short read lengths precluded classification as asRNAs or sRNAs). The 5 candidate ncRNAs were chosen to be representative of the various locations of ncRNAs on the genome with respect to adjacent genes and included putative transcripts found opposite to either the 5' or 3' ends of genes or classified as opposite intergenic (see additional file 3: Figs S5 and S6). The detection of antisense and ncRNAs using these techniques suggests that the observed antisense transcriptional initiation from the dRNA-seq data is not an artefact of library construction. A functional analysis of the genes opposite to asRNAs revealed that the most highly represented categories contained genes related to pathogenic or chemotaxis and motility functions (additional file 3: Figure S7). Similar to the TSSs located directly upstream of protein coding ORFs, we found that 79% of the candidate asRNA transcripts started with a purine residue (A - 51.40%, G - 27.70%). This provides further evidence that the asRNAs are primary transcripts rather than processing fragments or artefacts of the sequencing protocol, and reflects a recent study in E. coli where it was shown that 74% of the asRNA transcripts began with a purine .
A MEME analysis of the 302 asRNA TSSs revealed a strongly conserved -10 binding site in 280 promoters (TATAAT), however the -35 site was only weakly conserved (Figure 3E). Since the -35 region has been shown to enhance stability of the RNAP-promoter complex, this could suggest that the majority of asRNAs in S. Typhimurium are a consequence of promiscuous transcription initiation, as has been suggested for E. coli and eukaryotes . One of the mechanisms by which asRNAs opposing 5'UTRs may inhibit translation is by obscuring the ribosome binding site (RBS) . None of the ncRNAs we identified which were opposite to 5'UTRs appeared to obscure the RBS, however, the short Illumina read lengths may preclude this possibility. Alternately, the ncRNA may prevent transcription via transcriptional interference or attenuation . Indeed, some of the genes that were opposite to asRNAs were transcriptionally silent, suggesting a possible role for asRNAs in their regulation (e.g.SLaRNA310 which is antisense to the 3' end of SL1344_2729). Interestingly we found candidate asRNAs and ncRNAs to 41% of the rRNA genes (see additional file 2: Table S5); similarly in H. pylori, ~28% of the tRNA and rRNA genes were found to have antisense TSSs . In addition to stable RNA genes, we discovered putative candidate asRNAs to 18 virulence genes, 9 and 7 of which are located on the opposite strand to SPI1 and SPI2 encoded genes (see additional file 2: Table S5). Their potential role in regulating the expression of these virulence genes is currently being investigated.
Defining ppGpp-dependent gene expression using dRNA-seq
ppGpp-dependent transcriptional start sites within the S.Typhimurium genome
% ppGpp activated
ORF P, I
ORF S, I
ORF S, I, Alt
Known & predicted sRNAs
New candidate sRNAS
Of the ppGpp-activated genes, by far the largest functional category contained pathogenicity-related genes (22 genes), which supports our previous microarray based analysis of ppGpp-dependent virulence gene regulation in S. Typhimurium (Figure 8). We also found 8 ppGpp-activated genes to possess regulatory functions. These include rtsA and flhD which encode major transcriptional activators of SPI1 and flagella biosynthesis respectively [64, 65]. Previous work has shown that ppGpp-activated genes such as amino acid biosynthetic genes tend to contain AT-rich discriminator regions which allow optimal binding with the σ-subunit of RNAP [62, 66]. In confirmation, a MEME comparison of the ppGpp-activated genes revealed a tendency towards AT-rich discriminator regions (an average of 68% AT content for ppGpp activated promoters compared to 57% for all promoters), however no conserved motifs could be identified using MEME.
ppGpp-dependent expression of non-coding RNAs
Of the total known and predicted sRNAs we found that 18% of the TSSs (18) were ppGpp-dependent out of a total of 100 start sites, and 25% (15) of the new candidate sRNAs were ppGpp-dependent out of a total of 65 start sites (Table 2, additional file 2: Tables S3 and S4). This is less than the proportion of ppGpp-dependent TSSs identified for SL1344 chromosomal ORFs (34%). Similar ppGpp dependent control of small non-coding RNA abundance has been observed in other bacteria including Rhizobium etli and Staphylococcus aureus[29, 67]. Of the total number of ppGpp-dependent known, predicted and new candidate sRNAs, 18 were elevated and 15 repressed by ppGpp (Table 2). As noted previously, a characteristic of ppGpp-repressed genes is the presence of a GC rich discriminator region located between the TSS and the -10 region; however, there was no clear correlation between the GC content of the sRNA discriminator region and fold-repression or activation by ppGpp (data not shown). This suggests that the majority of ppGpp repressed sRNAs may be indirectly regulated, or the size of the dataset was too small to identify a conserved motif. We confirmed the dRNA-seq defined ppGpp-activation of of 2 sRNAs (STnc1020 and InvR) by Northern blotting, one of which (InvR) was previously found to be ppGpp-dependently elevated  (additional file 3: Figure S5, additional file 2: Table S3). We also confirmed the dRNA-seq defined ppGpp-repression of RprA (additional file 3: Figure S5, additional file 2: Table S3).
Of the 302 asRNAs that were directly opposite ORFs, we note that 32 were ppGpp-dependently elevated and 32 repressed (Table 2, additional file 2: Table S5). This represents a total of 21% of the candidate asRNAs and is similar to the percentage of ppGpp-dependent sRNAs. Interestingly we note that antisense transcripts to the sipA and invH genes were ppGpp-activated by 4.3 and 17-fold respectively. The remaining 190 start sites assigned to ncRNAs, (which we could not unambiguously identify as asRNAs) showed a similar level of ppGpp-dependency to the antisense and sRNA TSSs (19%).
We have determined the TSSs for 78% of the S. Typhimurium ORFs during growth conditions in which model the extracellular virulence gene expression programme (STEX). To date this is the most extensive and accurate map of the TSSs for this bacterium. Our analysis also identified secondary TSSs for many genes and operon structures. Our MEME based promoter analysis of the first genes of operons identified conserved regions in the promoters which were found to closely resemble consensus binding sites for σ24, σ28, σ32, σ38 & σ54 factors; many of the predicted sigma factor-dependent genes had previously been experimentally verified in either E. coli or Salmonella. We verified the expression of 38 out of 87 predicted sRNAs and 45 out of 62 known sRNAs [28, 68](and J. Vogel; pers. comm.) and also extended the repertoire of sRNAs encoded within the S. Typhimurium genome by 55. Of the predicted sRNAs we were unable to verify, it is possible that they were not expressed under the growth condition studied here. We also observed that the location of the TSSs of a subset of the predicted sRNAs did not correspond to the predicted start sites; from this we infer the bioinformatic approach used to identify the TSSs may require experimental-based refinements to enhance accuracy. We identified 302 candidate antisense transcripts for the S. Typhimurium genome for which we defined a conserved -10 hexamer upstream of the TSS. Although from this study, we cannot rule out the possibility that the expression of asRNAs are a result of promiscuous transcription, other work suggests this is not the case, at least for H. pylori.
Our dRNA-seq approach to identifying ppGpp-dependent transcription was validated by comparison with a microarray-based determination of the ppGpp-dependent transcriptome performed under identical growth conditions. The GC rich discriminator region located between the TSS and -10 region of ppGpp-repressed genes has been shown to play a role in destabilising the RNAP-ppGpp complex of rRNA promoters . We were able to correlate decreased transcript levels of ppGpp-repressed genes with the abundance of GC residues within the discriminator region . Our data showed no correlation between AT content of the discriminator region and the level of ppGpp-activation. However in agreement with Da Costa et al, we did find that in general, ppGpp-activated genes contained a higher overall discriminator AT content. Interestingly we note that SPI1 and SPI2 encoded genes contain AT-rich discriminator regions and the only sigma factor known to contribute to SPI1 expression is σ70; this suggests the possibility of a direct activation of SPI1 regulatory genes by ppGpp, rather than via sigma factor competition, as has already been suggested [9, 71].
Many regulons controlled by alternative sigma factors, including σ38 and σ32 are poorly induced in cells lacking ppGpp . In order to determine whether this was also the case for S. Typhimurium, we analysed our alternative sigma factor promoter database for ppGpp-dependency. We found that almost all the genes belonging to the σ28 and σ32 regulons and more than half of the σ24 -dependent genes were ppGpp-repressed. In contrast, the σ38 regulon showed no tendency towards ppGpp-activation or repression (additional file 2: Table S1). Previous work has also shown that, in contrast to E. coli, ppGpp does not control RpoS levels in S. Typhimurium during late-log and stationary phase growth . We conjecture that the ppGpp-repression of some of the alternative sigma factor regulons may represent an adaptation to favour σ70 dependent virulence gene expression under the STEX growth conditions studied here.
It is generally accepted that elevated levels of ppGpp during amino acid starvation (stringent response) result in repression of stable RNAs (rRNA and tRNA). Consistent with this we observed repression of the rRNA operons in the wild-type compared to the ΔrelA ΔspoT strain (additional file 2: Table S6). However, all but one of the tRNA mono- and polycistronic operons showed elevated transcript levels in the wild-type compared to the ΔrelA ΔspoT strain; a similar ppGpp-dependent elevation of tRNA levels was found in stationary phase Rhizobium etli relative to early exponential phase . It is possible that elevation of tRNA levels could be a consequence of ppGpp-dependent differential processing or stability rather than direct ppGpp-dependent regulation. Indeed tRNA has been reported to remain stable under starvation conditions that induce rRNA degradation in E. coli. In support of the possibility of ppGpp-dependent differential processing or stability of tRNAs we observed that expression of RNaseP, a ribozyme responsible for 5' end processing of tRNAs, was ppGpp-activated in S. Typhimurium (additional file 2: Table S1). Similarly, we note that the R. etli RNaseP was also ppGpp-activated (29). We hypothesise that the ppGpp-dependent activation of RNaseP may result in reduced tRNA processing in the ΔrelA ΔspoT strain and subsequent removal of incorrectly processed tRNAs via RNA quality control mechanisms .
For the known and predicted sRNAs described in this study, a MEME analysis was able to identify conserved -10 (TATTNT) and -35 (TTGaCA) regions upstream of the predicted TSSs (Figure 3D). A manual inspection of the smaller new candidate sRNA dataset identified AT rich -10 hexamers in 69% of the promoters (data not shown). A manual inspection of all of the sRNA promoters described in this study failed to find any of the well-defined alternative sigma factor binding motifs and in fact only four sRNAs, (MicA, RybB, GlmZ and GlmY) have so far been shown to be positively controlled by σE and σ54 in E. coli. This suggests that, at least for the sRNAs transcribed under the growth conditions studied here, their expression is mostly σ70 dependent and perhaps reflects the major role sRNAs play in maintaining house-keeping functions and regulating virulence determinants. In contrast to the discriminator regions of ppGpp-repressed genes, we were unable to identify a conserved GC rich region in the set of ppGpp-repressed sRNAs. Several of the ppGpp-repressed sRNAs (OmrA, OmrB, MicA, MicF and CyaR) have been shown to act as repressors of genes encoding porins and outer membrane proteins (OMPs) suggesting that ppGpp may indirectly activate these target genes. OMPs are important virulence factors and play a significant role in the bacterial adaptation to environmental conditions. Other highly ppGpp-repressed sRNAs shown to play a role in Salmonella virulence include IsrI, IsrP and CsrB. In addition, the sRNA chaperone, Hfq was ppGpp-repressed by 5.6-fold thus expanding the role of ppGpp in the regulation of Salmonella virulence gene expression (Figure 5) . The IsrI and IsrP sRNAs are expressed during infection of J774 macrophages . IsrI is also expressed during stationary phase, and under low oxygen or magnesium levels; IsrP is expressed under low magnesium and extreme acid conditions of pH2.5 . CsrB is part of the csr system shown to play a role in the regulation of invasion gene expression in S. Typhimurium . The SPI1 encoded sRNA, InvR, has previously been reported to be ppGpp-activated . Our data confirms InvR as the most highly ppGpp-activated sRNA we detected under these growth conditions (12.5-fold; additional file 2: Table S1). Although Hfq has been shown to reduce the stability of InvR , we note that despite a 5.6-fold ppGpp-dependent repression of Hfq trasncript levels, InvR remains highly ppGpp-activated. This suggests that ppGpp is able to modulate InvR transcript levels via a Hfq independent mechanism. InvR represses synthesis of the major outer membrane protein, OmpD . It is suggested sRNAs such as InvR have evolved to modulate OMP levels, which can be deleterious to the cell . OmpD has also been shown to facilitate Salmonella adherence to human macrophages and intestinal epithelial cell lines [78, 79]. Potential targets for the new ppGpp-dependent sRNAs include fabH, involved in the initiation of fatty acid biosynthesis, 6 genes involved in transport of sugars, nitrite, peptides and branched chain amino acids, and 3 transcriptional regulators, nadR, rob, and STM2275 (see additional file 2: Table S4).
We discovered extensive antisense transcription within the S. Typhimurium genome under the growth conditions studied here. Similarly, a considerable abundance of asRNA transcription was also discovered in E. coli and H. pylori[24, 54]. Interestingly, we observed candidate asRNAs to several virulence genes from SPI1, 2 and 6 and identified 4 putative ncRNAs which were classified as opposite intergenic between genes encoding several major SPI1 regulators including hilA and hilD (additional file 2: Table S5). One of the two candidate ncRNAs between hilA and hilD (Sla0508) was highly ppGpp-activated by a factor of at least 34-fold. HilD has been implicated in cross-talk between SPI1 and SPI2 expression . Under the growth conditions used in this study, SPI2 is not highly expressed compared to SPI1. It is therefore tantalising to suggest that ncRNAs may play a role in modulating expression of the STEX and STIN virulence gene expression programmes.
Here we used dRNA-seq to define the transcriptomic architecture for an S. Typhimurium wild-type and ppGpp0 strain during growth conditions where the invasion (SPI1) genes are expressed. We identify the precise location of the TSSs for 78% of the S. Typhimurium genome, the reannotation of 60 start codons and the identification of 23 potential new ORFs. The nucleotide position of the TSSs enabled us to perform a promoter analysis, which resulted in the prediction of binding sites for 6 sigma factors, the analysis of 5' leader lengths and the prediction of 625 operons. The definition of the ncRNA transcriptome resulted in confirmation of the expression of 83 predicted and known sRNAs under these growth conditions and the prediction of 55 new candidate sRNAs, for which potential targets were inferred. Extensive asRNA transcription was also discovered for 302 candidate asRNAs, 18 and 11 of which were opposite virulence genes and candidate sRNAs respectively. The dRNA-seq predicted ppGpp-dependent TSSs and ppGpp-dependent expression for the SL1344 genome and we showed that ppGpp is involved in regulating an average of 20% of S. Typhimurium ncRNA expression.
Bacterial Strains and Growth Conditions
The wild-type, virulent S. Typhimurium strain SL1344 (31) was provided by F. Norel (Institut Pasteur) and re-isolated from the spleens of infected BALB/c mice. The SL1344 ΔrelA ΔspoT strain was a kind gift from Dr. Karsten Tedin, Freie Universität, Berlin, and was constructed by lambda red mutagenesis . Deletion of the relA and spoT genes and intactness of flanking regions were determined by sequencing; antibiotic resistance markers were removed after construction (pers. comm. Dr. K. Tedin). Bacterial cultures were grown overnight in Luria-Bertani broth (LB) at 37°C, 250 rpm, from -70°C glycerol stocks and used to inoculate into 50 ml of fresh LB in 250 ml conical flasks. The cultures were grown aerobically at 37°C with shaking at 250 rpm to an OD600 of 2.1, (Early Stationary Phase; ESP), conditions previously shown to induce SPI1 gene expression . We confirmed that the wild-type and ΔrelA ΔspoT strains have almost identical growth rates in LB under these growth conditions (additional file 3: Table S8; 9) and that the wild-type is highly invasive compared to the ΔrelA ΔspoT strain in a HeLa cell infection assay (additional file 3: Table S9; 82).
RNA Extraction and Purification
Bacterial cultures were harvested at ESP, added to one-fifth volume of stop solution (5% (v/v) phenol in ethanol), and incubated on ice for 30 min to stabilize total RNA . Bacterial cells were harvested at 10,000 rpm at 4°C for 5 min and re-suspended in 500 μl of re-suspension buffer (25 mM Tris-HCl pH7.4, 1 mM EDTA). An equal volume of lysis solution (0.6 M Sodium acetate pH5.2, 4 mM EDTA, 3% SDS) was added and the mixture was boiled for 30 - 60 sec or until the suspension cleared. The solution was then incubated at 20 to 25°C for 5 min before extracting twice with phenol and chloroform. Total RNA was precipitated overnight at -20°C with 2.5 volumes of ice-cold ethanol and pelleted at 13,500 rpm for 60 min at 4°C. The pellet was washed once with 70% ethanol and vacuum dried for 5 min. Finally, the pellet was resuspended in 100 μl of nuclease free water. Chromosomal DNA was removed by digestion with 50-100 units of Turbo DNA-free DNase (Ambion). The removal of contaminating DNA was verified by performing PCR using primers targeted to bacterial housekeeping genes. The quantity and quality of the total RNA was determined using a 2100 Bioanalyzer™(Agilent) before and after DNase treatment.
Library preparation and sequencing
To differentiate primary transcripts with a 5' triphosphate end (5' PPP) from processed transcripts with a 5' monophosphate end (5' P), total RNA from each strain was divided into equal amounts and one half was treated with Terminator™ 5'-phosphate-dependent exonuclease (TEX, Epicentre Biotechnologies) as previously described . TEX specifically degrades RNAs with a 5' P end but does not degrade transcripts with a 5' PPP end and therefore enriches for primary transcripts. Both the untreated (Non-enriched (NE)) and treated (enriched (EN)) samples were then treated with Tobacco acid pyrophosphatase (TAP; Epicentre Biotechnologies) to generate 5'-mono-phosphate ends for linker ligation. After 5' linker ligation and poly(A) tailing, strand-specific cDNA libraries were constructed by Vertis Biotechnology, AG, Germany (http://www.vertis-biotech.com) as described . For each Roche-454 and Illumina-Solexa sequencing run, four strand specific cDNA libraries were prepared: Wild-type non-enriched (Wt-NE), Wild-type enriched (Wt-EN), ΔrelA ΔspoT non-enriched (DM-NE) and ΔrelA ΔspoT enriched (DM-EN). Each library was tagged with a different barcode at the 5' end to enable multiplexing during sequencing. Four of the Roche-454 cDNA libraries were pooled and sequenced (Liverpool University, UK), yielding a total of ~400,000 reads. Similarly, four of the Illumina-Solexa libraries were pooled and sequenced in a single lane using 36 single-read cycles on an Illumina Genome Analyzer II sequencing machine (GATC Biotech, Germany), yielding a total of ~107 reads.
Analysis of sequences and statistics
The mapping statistics for both the Roche-454 and Illumina-Solexa cDNA libraries are shown in additional file 1: Table S1. Following sequencing, custom PERL scripts were used to separate the cDNA libraries based on their barcodes and to remove 5' linker regions. Any Roche-454 sequencing reads less than 12 nt in length were removed to avoid mapping errors. Both Roche-454 and Illumina-Solexa generated reads were mapped onto the S. Typhimurium SL1344 genome (Genbank ID. FQ312003.1) including the three native virulence plasmids (SLP1, SLP2 and SLP3) using the segemehl program . Mapped reads were converted to a graph file and visualized on an Integrated Genome Browser (IGB) .
5' RACE determination of selected TSSs
The transcriptional start sites of selected genes were determined using 5' RACE as described previously . Briefly, 12 μg of total RNA was treated with TAP and the RNA oligonucleotide adaptor A3 was ligated to the 5' end of the treated RNA. TAP cleaves the 5'-triphosphate of primary transcripts to a monophosphate, thus making them available for ligation of the RNA adaptor. This results in an enrichment of 5'-RACE products for primary transcripts in TAP treated RNA, compared to an untreated control. First strand cDNA synthesis was performed using either random hexamers oligonucleotide primers (for SL1344_1204 and SL1344_1122) or gene-specific primers followed by PCR amplification with nested gene-specific primers and 5' Adaptor-specific DNA primer B6. Resulting PCR products were cloned into the pGEM®-T Easy vector (Promega) and sequenced using standard protocols. All primers used are detailed in additional file 1: Table S2.
Adaptor assisted RT-PCR to detect asRNAs
In general we found that Northern blotting failed to detect non-coding RNAs (ncRNAs) which had Solexa read numbers of less than 200. We therefore employed an RT-PCR technique to selectively amplify predicted ncRNAs and thus confirm their existence. Primers complementary to the 3' end of putative ncRNA but including a universal 5' adaptor sequence were used to prime cDNA synthesis from total RNA samples. RNA (0.5 μg) was mixed with 2 pmol of the RT primer and reverse transcribed at 55°C for 1 hour using AffinityScript (Agilent) according to the manufacturer's instructions. Following heat denaturation of the enzyme, 2 μl of the reaction was used as a template in a PCR reaction using a primer matching the 5' end of the putative ncRNA and a second universal primer (U5), matching the 5' adaptor sequence of the RT primer. For each target, PCR reactions were carried out with cDNA from both wild-type and ΔrelA ΔspoT strains as well as a genomic DNA negative control. Primers used are listed in additional file 1: Table S4.
Microarray analysis was performed as described previously . RNA was extracted from S. Typhimurium SL1344 wild type and isogenic ΔrelA ΔspoT strains as described above under identical growth conditions and grown to the same OD. The RNA was labelled and hybridised to IFR SALSA2 whole ORF microarrays (http://www.ifr.ac.uk/Safety/Microarrays/default.html#protocols), and data processed and analysed using GeneSpring™ (Agilent). The data was from 3 biological replicates, statistically filtered (P = 0.05) and a 2-fold cut off applied.
Microarray accession number
The microarray data discussed in this publication are MIAME compliant and have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE34269 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=gse34269).
RNA sequencing - transcriptomic studies utilising high-throughput deep sequencing of cDNA libraries
differential RNA sequencing - RNA-seq based technique that differentiates between primary and processed transcripts
transcription start site(s)
open reading frame(s)
Salmonella pathogenicity island
Salmonella plasmids 1-3
integrated genome browser
basic local alignment search tool
Shine Dalgarno sequence
ribosome binding site
type 3 secretion system
- STEX STIN:
Salmonella extracellular and intracellular virulence gene expression programs respectively
tobacco acid pyrophosphatase.
We thank Joerg Vogel for helpful discussions and Arnoud van Vliet for providing the adapter assisted PCR method. NS and AT were supported by a core funded programme from the Biotechnology Biological Science Research Council (BBSRC), UK. VKR was supported by BBSRC grant number BB/F00978X/1.
- Chimalizeni Y, Kawaza K, Molyneux E: The epidemiology and management of non typhoidal Salmonella infections. Advances in experimental medicine and biology. 659: 33-46.
- Hensel M: Evolution of pathogenicity islands of Salmonella enterica. Int J Med Microbiol. 2004, 294 (2-3): 95-102. 10.1016/j.ijmm.2004.06.025.View ArticlePubMedGoogle Scholar
- Hensel M: Salmonella pathogenicity island 2. Molecular microbiology. 2000, 36 (5): 1015-1023. 10.1046/j.1365-2958.2000.01935.x.View ArticlePubMedGoogle Scholar
- Cummings LA, Barrett SL, Wilkerson WD, Fellnerova I, Cookson BT: FliC-specific CD4+ T cell responses are restricted by bacterial regulation of antigen expression. J Immunol. 2005, 174 (12): 7929-7938.View ArticlePubMedGoogle Scholar
- Thompson A, Rolfe MD, Lucchini S, Schwerk P, Hinton JC, Tedin K: The bacterial signal molecule, ppGpp, mediates the environmental regulation of both the invasion and intracellular virulence gene programs of Salmonella. J Biol Chem. 2006, 281 (40): 30112-30121. 10.1074/jbc.M605616200.View ArticlePubMedGoogle Scholar
- Deiwick J, Nikolaus T, Erdogan S, Hensel M: Environmental regulation of Salmonella pathogenicity island 2 gene expression. Molecular microbiology. 1999, 31 (6): 1759-1773. 10.1046/j.1365-2958.1999.01312.x.View ArticlePubMedGoogle Scholar
- Potrykus K, Cashel M: (p)ppGpp: still magical?. Annu Rev Microbiol. 2008, 62: 35-51. 10.1146/annurev.micro.62.081307.162903.View ArticlePubMedGoogle Scholar
- Mittenhuber G: Comparative genomics and evolution of genes encoding bacterial (p)ppGpp synthetases/hydrolases (the Rel, RelA and SpoT proteins). J Mol Microbiol Biotechnol. 2001, 3 (4): 585-600.PubMedGoogle Scholar
- Pizarro-Cerda J, Tedin K: The bacterial signal molecule, ppGpp, regulates Salmonella virulence gene expression. Molecular microbiology. 2004, 52 (6): 1827-1844. 10.1111/j.1365-2958.2004.04122.x.View ArticlePubMedGoogle Scholar
- Klinkenberg LG, Lee JH, Bishai WR, Karakousis PC: The stringent response is required for full virulence of Mycobacterium tuberculosis in guinea pigs. J Infect Dis. 202 (9): 1397-1404.
- Primm TP, Andersen SJ, Mizrahi V, Avarbock D, Rubin H, Barry CE: The stringent response of Mycobacterium tuberculosis is required for long-term survival. Journal of bacteriology. 2000, 182 (17): 4889-4898. 10.1128/JB.182.17.4889-4898.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Taylor CM, Beresford M, Epton HA, Sigee DC, Shama G, Andrew PW, Roberts IS: Listeria monocytogenes relA and hpt mutants are impaired in surface-attached growth and virulence. Journal of bacteriology. 2002, 184 (3): 621-628. 10.1128/JB.184.3.621-628.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Hammer BK, Tateda ES, Swanson MS: A two-component regulator induces the transmission phenotype of stationary-phase Legionella pneumophila. Molecular microbiology. 2002, 44 (1): 107-118. 10.1046/j.1365-2958.2002.02884.x.View ArticlePubMedGoogle Scholar
- Zusman T, Gal-Mor O, Segal G: Characterization of a Legionella pneumophila relA insertion mutant and roles of RelA and RpoS in virulence gene expression. Journal of bacteriology. 2002, 184 (1): 67-75. 10.1128/JB.184.1.67-75.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Haralalka S, Nandi S, Bhadra RK: Mutation in the relA gene of Vibrio cholerae affects in vitro and in vivo expression of virulence factors. Journal of bacteriology. 2003, 185 (16): 4672-4682. 10.1128/JB.185.16.4672-4682.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Erickson DL, Lines JL, Pesci EC, Venturi V, Storey DG: Pseudomonas aeruginosa relA contributes to virulence in Drosophila melanogaster. Infect Immun. 2004, 72 (10): 5638-5645. 10.1128/IAI.72.10.5638-5645.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Na HS, Kim HJ, Lee HC, Hong Y, Rhee JH, Choy HE: Immune response induced by Salmonella typhimurium defective in ppGpp synthesis. Vaccine. 2006, 24 (12): 2027-2034. 10.1016/j.vaccine.2005.11.031.View ArticlePubMedGoogle Scholar
- Cashel M, Gentry DM, Hernandez VJ, Vinella D: The stringent response. 1996, Washington DC: ASM Press, 1:Google Scholar
- Jishage M, Kvint K, Shingler V, Nystrom T: Regulation of sigma factor competition by the alarmone ppGpp. Genes Dev. 2002, 16 (10): 1260-1270. 10.1101/gad.227902.PubMed CentralView ArticlePubMedGoogle Scholar
- Paul BJ, Berkmen MB, Gourse RL: DksA potentiates direct activation of amino acid promoters by ppGpp. Proc Natl Acad Sci USA. 2005, 102 (22): 7823-7828. 10.1073/pnas.0501170102.PubMed CentralView ArticlePubMedGoogle Scholar
- Potrykus K, Wegrzyn G, Hernandez VJ: Direct stimulation of the lambda paQ promoter by the transcription effector guanosine-3',5'-(bis)pyrophosphate in a defined in vitro system. J Biol Chem. 2004, 279 (19): 19860-19866. 10.1074/jbc.M313378200.View ArticlePubMedGoogle Scholar
- Perederina A, Svetlov V, Vassylyeva MN, Tahirov TH, Yokoyama S, Artsimovitch I, Vassylyev DG: Regulation through the secondary channel--structural framework for ppGpp-DksA synergism during transcription. Cell. 2004, 118 (3): 297-309. 10.1016/j.cell.2004.06.030.View ArticlePubMedGoogle Scholar
- Sorek R, Cossart P: Prokaryotic transcriptomics: a new view on regulation, physiology and pathogenicity. Nat Rev Genet. 2010, 11 (1): 9-16.View ArticlePubMedGoogle Scholar
- Sharma CM, Hoffmann S, Darfeuille F, Reignier J, Findeiss S, Sittka A, Chabas S, Reiche K, Hackermuller J, Reinhardt R: The primary transcriptome of the major human pathogen Helicobacter pylori. Nature. 464 (7286): 250-255.
- Thomason MK, Storz G: Bacterial antisense RNAs: how many are there, and what are they doing?. Annu Rev Genet. 2010, 44: 167-188. 10.1146/annurev-genet-102209-163523.PubMed CentralView ArticlePubMedGoogle Scholar
- van Vliet AH: Next generation sequencing of microbial transcriptomes: challenges and opportunities. FEMS microbiology letters. 2010, 302 (1): 1-7. 10.1111/j.1574-6968.2009.01767.x.View ArticlePubMedGoogle Scholar
- Irnov I, Sharma CM, Vogel J, Winkler WC: Identification of regulatory RNAs in Bacillus subtilis. Nucleic acids research. 2010, 38 (19): 6637-6651. 10.1093/nar/gkq454.PubMed CentralView ArticlePubMedGoogle Scholar
- Pfeiffer V, Sittka A, Tomer R, Tedin K, Brinkmann V, Vogel J: A small non-coding RNA of the invasion gene island (SPI-1) represses outer membrane protein synthesis from the Salmonella core genome. Molecular microbiology. 2007, 66 (5): 1174-1191. 10.1111/j.1365-2958.2007.05991.x.View ArticlePubMedGoogle Scholar
- Vercruysse M, Fauvart M, Jans A, Beullens S, Braeken K, Cloots L, Engelen K, Marchal K, Michiels J: Stress response regulators identified through genome-wide transcriptome analysis of the (p)ppGpp-dependent response in Rhizobium etli. Genome Biol. 2011, 12 (2): R17-10.1186/gb-2011-12-2-r17.PubMed CentralView ArticlePubMedGoogle Scholar
- Hawley DK, McClure WR: Compilation and analysis of Escherichia coli promoter DNA sequences. Nucleic Acids Res. 1983, 11 (8): 2237-2255. 10.1093/nar/11.8.2237.PubMed CentralView ArticlePubMedGoogle Scholar
- Fozo EM, Kawano M, Fontaine F, Kaya Y, Mendieta KS, Jones KL, Ocampo A, Rudd KE, Storz G: Repression of small toxic protein synthesis by the Sib and OhsC small RNAs. Molecular microbiology. 2008, 70 (5): 1076-1093. 10.1111/j.1365-2958.2008.06394.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Kawano M, Oshima T, Kasai H, Mori H: Molecular characterization of long direct repeat (LDR) sequences expressing a stable mRNA encoding for a 35-amino-acid cell-killing peptide and a cis-encoded small antisense RNA in Escherichia coli. Molecular microbiology. 2002, 45 (2): 333-349. 10.1046/j.1365-2958.2002.03042.x.View ArticlePubMedGoogle Scholar
- Ide N, Ikebe T, Kutsukake K: Reevaluation of the promoter structure of the class 3 flagellar operons of Escherichia coli and Salmonella. Genes & Genetic Systems. 1999, 74 (3): 113-116. 10.1266/ggs.74.113.View ArticleGoogle Scholar
- Ibanez-Ruiz M, Robbe-Saule V, Hermant D, Labrude S, Norel F: Identification of RpoS (sigma(S))-regulated genes in Salmonella enterica serovar typhimurium. Journal of bacteriology. 2000, 182 (20): 5749-5756. 10.1128/JB.182.20.5749-5756.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Loewen PCH B, Strutinsky J, Sparling R: Regulation of rpoS regulon of Escherichia coli. Canadian Journal of Microbiology. 1998, 44: 707-717.View ArticleGoogle Scholar
- Zhao K, Liu M, Burgess RR: The Global Transcriptional Response of Escherichia coli to Induced σ32 Protein Involves σ32 Regulon Activation Followed by Inactivation and Degradation of σ32 in Vivo. Journal of Biological Chemistry. 2005, 280 (18): 17758-17768.View ArticlePubMedGoogle Scholar
- Laursen BS, Sorensen HP, Mortensen KK, Sperling-Petersen HU: Initiation of protein synthesis in bacteria. Microbiol Mol Biol Rev. 2005, 69 (1): 101-123. 10.1128/MMBR.69.1.101-123.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Shine J, Dalgarno L: The 3'-terminal sequence of Escherichia coli 16S ribosomal RNA: complementarity to nonsense triplets and ribosome binding sites. Proc Natl Acad Sci USA. 1974, 71 (4): 1342-1346. 10.1073/pnas.71.4.1342.PubMed CentralView ArticlePubMedGoogle Scholar
- Grundy FJ, Henkin TM: From ribosome to riboswitch: control of gene expression in bacteria by RNA structural rearrangements. Crit Rev Biochem Mol Biol. 2006, 41 (6): 329-338. 10.1080/10409230600914294.View ArticlePubMedGoogle Scholar
- Krishnan KM, Van Etten WJ, Janssen GR: Proximity of the start codon to a leaderless mRNA's 5' terminus is a strong positive determinant of ribosome binding and expression in Escherichia coli. Journal of bacteriology. 192 (24): 6482-6485.
- Brock JE, Pourshahian S, Giliberti J, Limbach PA, Janssen GR: Ribosomes bind leaderless mRNA in Escherichia coli through recognition of their 5'-terminal AUG. RNA. 2008, 14 (10): 2159-2169. 10.1261/rna.1089208.PubMed CentralView ArticlePubMedGoogle Scholar
- Walz A, Pirrotta V, Ineichen K: Lambda repressor regulates the switch between PR and Prm promoters. Nature. 1976, 262 (5570): 665-669. 10.1038/262665a0.View ArticlePubMedGoogle Scholar
- Klock G, Hillen W: Expression, purification and operator binding of the transposon Tn1721-encoded Tet repressor. J Mol Biol. 1986, 189 (4): 633-641. 10.1016/0022-2836(86)90493-6.View ArticlePubMedGoogle Scholar
- O'Donnell SM, Janssen GR: The initiation codon affects ribosome binding and translational efficiency in Escherichia coli of cI mRNA with or without the 5' untranslated leader. Journal of bacteriology. 2001, 183 (4): 1277-1283. 10.1128/JB.183.4.1277-1283.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Sullivan MJ, Curson AR, Shearer N, Todd JD, Green RT, Johnston AW: Unusual regulation of a leaderless operon involved in the catabolism of dimethylsulfoniopropionate in Rhodobacter sphaeroides. PLoS One. 2011, 6 (1): e15972-10.1371/journal.pone.0015972.PubMed CentralView ArticlePubMedGoogle Scholar
- Wurtzel O, Sapra R, Chen F, Zhu Y, Simmons BA, Sorek R: A single-base resolution map of an archaeal transcriptome. Genome Res. 2010, 20 (1): 133-141. 10.1101/gr.100396.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Cho BK, Zengler K, Qiu Y, Park YS, Knight EM, Barrett CL, Gao Y, Palsson BO: The transcription unit architecture of the Escherichia coli genome. Nat Biotechnol. 2009, 27 (11): 1043-1049. 10.1038/nbt.1582.View ArticlePubMedGoogle Scholar
- Tsui HC, Feng G, Winkler ME: Transcription of the mutL repair, miaA tRNA modification, hfq pleiotropic regulator, and hflA region protease genes of Escherichia coli K-12 from clustered Esigma32-specific promoters during heat shock. J Bacteriol. 1996, 178 (19): 5719-5731.PubMed CentralPubMedGoogle Scholar
- Dam P, Olman V, Harris K, Su Z, Xu Y: Operon prediction using both genome-specific and general genomic information. Nucleic acids research. 2007, 35 (1): 288-298.PubMed CentralView ArticlePubMedGoogle Scholar
- Mao F, Dam P, Chou J, Olman V, Xu Y: DOOR: a database for prokaryotic operons. Nucleic acids research. 2009, D459-463. 37 Database
- Messenger SL, Green J: FNR-mediated regulation of hyp expression in Escherichia coli. FEMS microbiology letters. 2003, 228 (1): 81-86. 10.1016/S0378-1097(03)00726-2.View ArticlePubMedGoogle Scholar
- Perkins TT, Kingsley RA, Fookes MC, Gardner PP, James KD, Yu L, Assefa SA, He M, Croucher NJ, Pickard DJ: A strand-specific RNA-Seq analysis of the transcriptome of the typhoid bacillus Salmonella typhi. PLoS Genet. 2009, 5 (7): e1000569-10.1371/journal.pgen.1000569.PubMed CentralView ArticlePubMedGoogle Scholar
- Tjaden B: TargetRNA: a tool for predicting targets of small RNA action in bacteria. Nucleic acids research. 2008, W109-113. 36 Web Server
- Dornenburg JE, Devita AM, Palumbo MJ, Wade JT: Widespread antisense transcription in Escherichia coli. MBio. 1 (1):
- Thomason MK, Storz G: Bacterial antisense RNAs: how many are there, and what are they doing?. Annual review of genetics. 44: 167-188.
- Muirhead H: Isoenzymes of pyruvate kinase. Biochemical Society transactions. 1990, 18 (2): 193-196.View ArticlePubMedGoogle Scholar
- Croucher NJ, Thomson NR: Studying bacterial transcriptomes using RNA-seq. Curr Opin Microbiol. 13 (5): 619-624.
- Oliver HF, Orsi RH, Ponnala L, Keich U, Wang W, Sun Q, Cartinhour SW, Filiatrault MJ, Wiedmann M, Boor KJ: Deep RNA sequencing of L. monocytogenes reveals overlapping and extensive stationary phase and sigma B-dependent transcriptomes, including multiple highly transcribed noncoding RNAs. BMC Genomics. 2009, 10: 641-10.1186/1471-2164-10-641.PubMed CentralView ArticlePubMedGoogle Scholar
- Bernardo LM, Johansson LU, Solera D, Skarfstad E, Shingler V: The guanosine tetraphosphate (ppGpp) alarmone, DksA and promoter affinity for RNA polymerase in regulation of sigma-dependent transcription. Molecular microbiology. 2006, 60 (3): 749-764. 10.1111/j.1365-2958.2006.05129.x.View ArticlePubMedGoogle Scholar
- Potrykus K, Murphy H, Philippe N, Cashel M: ppGpp is the major source of growth rate control in E. coli. Environ Microbiol. 2011, 13 (3): 563-575. 10.1111/j.1462-2920.2010.02357.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Nagalakshmi U, Wang Z, Waern K, Shou C, Raha D, Gerstein M, Snyder M: The transcriptional landscape of the yeast genome defined by RNA sequencing. Science. 2008, 320 (5881): 1344-1349. 10.1126/science.1158441.PubMed CentralView ArticlePubMedGoogle Scholar
- Haugen SP, Berkmen MB, Ross W, Gaal T, Ward C, Gourse RL: rRNA promoter regulation by nonoptimal binding of sigma region 1.2: an additional recognition element for RNA polymerase. Cell. 2006, 125 (6): 1069-1082. 10.1016/j.cell.2006.04.034.View ArticlePubMedGoogle Scholar
- Zhou YN, Jin DJ: The rpoB mutants destabilizing initiation complexes at stringently controlled promoters behave like "stringent" RNA polymerases in Escherichia coli. Proc Natl Acad Sci USA. 1998, 95 (6): 2908-2913. 10.1073/pnas.95.6.2908.PubMed CentralView ArticlePubMedGoogle Scholar
- Ellermeier CD, Slauch JM: RtsA and RtsB coordinately regulate expression of the invasion and flagellar genes in Salmonella enterica serovar Typhimurium. Journal of bacteriology. 2003, 185 (17): 5096-5108. 10.1128/JB.185.17.5096-5108.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Kutsukake K, Iino T, Komeda Y, Yamaguchi S: Functional homology of fla genes between Salmonella typhimurium and Escherichia coli. Mol Gen Genet. 1980, 178 (1): 59-67. 10.1007/BF00267213.View ArticlePubMedGoogle Scholar
- Da Costa XJ, Artz SW: Mutations that render the promoter of the histidine operon of Salmonella typhimurium insensitive to nutrient-rich medium repression and amino acid downshift. Journal of bacteriology. 1997, 179 (16): 5211-5217.PubMed CentralPubMedGoogle Scholar
- Anderson KL, Roberts C, Disz T, Vonstein V, Hwang K, Overbeek R, Olson PD, Projan SJ, Dunman PM: Characterization of the Staphylococcus aureus heat shock, cold shock, stringent, and SOS responses and their effects on log-phase mRNA turnover. J Bacteriol. 2006, 188 (19): 6739-6756. 10.1128/JB.00609-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Sittka A, Lucchini S, Papenfort K, Sharma CM, Rolle K, Binnewies TT, Hinton JC, Vogel J: Deep sequencing analysis of small noncoding RNA and mRNA targets of the global post-transcriptional regulator, Hfq. PLoS Genet. 2008, 4 (8): e1000163-10.1371/journal.pgen.1000163.PubMed CentralView ArticlePubMedGoogle Scholar
- Zacharias M, Goringer HU, Wagner R: Influence of the GCGC discriminator motif introduced into the ribosomal RNA P2- and tac promoter on growth-rate control and stringent sensitivity. The EMBO journal. 1989, 8 (11): 3357-3363.PubMed CentralPubMedGoogle Scholar
- Travers AA: Promoter sequence for stringent control of bacterial ribonucleic acid synthesis. Journal of bacteriology. 1980, 141 (2): 973-976.PubMed CentralPubMedGoogle Scholar
- Dalebroux ZD, Svensson SL, Gaynor EC, Swanson MS: ppGpp conjures bacterial virulence. Microbiol Mol Biol Rev. 74 (2): 171-199.
- Davis BD, Luger SM, Tai PC: Role of ribosome degradation in the death of starved Escherichia coli cells. Journal of bacteriology. 1986, 166 (2): 439-445.PubMed CentralPubMedGoogle Scholar
- Li Z, Reimers S, Pandit S, Deutscher MP: RNA quality control: degradation of defective transfer RNA. The EMBO journal. 2002, 21 (5): 1132-1138. 10.1093/emboj/21.5.1132.PubMed CentralView ArticlePubMedGoogle Scholar
- Johansen J, Rasmussen AA, Overgaard M, Valentin-Hansen P: Conserved small non-coding RNAs that belong to the sigmaE regulon: role in down-regulation of outer membrane proteins. J Mol Biol. 2006, 364 (1): 1-8. 10.1016/j.jmb.2006.09.004.View ArticlePubMedGoogle Scholar
- Padalon-Brauch G, Hershberg R, Elgrably-Weiss M, Baruch K, Rosenshine I, Margalit H, Altuvia S: Small RNAs encoded within genetic islands of Salmonella typhimurium show host-induced expression and role in virulence. Nucleic acids research. 2008, 36 (6): 1913-1927. 10.1093/nar/gkn050.PubMed CentralView ArticlePubMedGoogle Scholar
- Altier C, Suyemoto M, Lawhon SD: Regulation of Salmonella enterica serovar typhimurium invasion genes by csrA. Infect Immun. 2000, 68 (12): 6790-6797. 10.1128/IAI.68.12.6790-6797.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Douchin V, Bohn C, Bouloc P: Down-regulation of porins by a small RNA bypasses the essentiality of the regulated intramembrane proteolysis protease RseP in Escherichia coli. J Biol Chem. 2006, 281 (18): 12253-12259. 10.1074/jbc.M600819200.View ArticlePubMedGoogle Scholar
- Hara-Kaonga B, Pistole TG: OmpD but not OmpC is involved in adherence of Salmonella enterica serovar typhimurium to human cells. Can J Microbiol. 2004, 50 (9): 719-727. 10.1139/w04-056.View ArticlePubMedGoogle Scholar
- Negm RS, Pistole TG: Macrophages recognize and adhere to an OmpD-like protein of Salmonella typhimurium. FEMS Immunol Med Microbiol. 1998, 20 (3): 191-199. 10.1111/j.1574-695X.1998.tb01127.x.View ArticlePubMedGoogle Scholar
- Bustamante VH, Martinez LC, Santana FJ, Knodler LA, Steele-Mortimer O, Puente JL: HilD-mediated transcriptional cross-talk between SPI-1 and SPI-2. Proc Natl Acad Sci USA. 2008, 105 (38): 14591-14596. 10.1073/pnas.0801205105.PubMed CentralView ArticlePubMedGoogle Scholar
- Datsenko KA, Wanner BL: One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci USA. 2000, 97 (12): 6640-6645. 10.1073/pnas.120163297.PubMed CentralView ArticlePubMedGoogle Scholar
- Song M, Kim HJ, Kim EY, Shin M, Lee HC, Hong Y, Rhee JH, Yoon H, Ryu S, Lim S: ppGpp-dependent stationary phase induction of genes on Salmonella pathogenicity island 1. J Biol Chem. 2004, 279 (33): 34183-34190. 10.1074/jbc.M313491200.View ArticlePubMedGoogle Scholar
- Hoffmann S, Otto C, Kurtz S, Sharma CM, Khaitovich P, Vogel J, Stadler PF, Hackermuller J: Fast mapping of short sequences with mismatches, insertions and deletions using index structures. PLoS Comput Biol. 2009, 5 (9): e1000502-10.1371/journal.pcbi.1000502.PubMed CentralView ArticlePubMedGoogle Scholar
- Nicol JW, Helt GA, Blanchard SG, Raja A, Loraine AE: The Integrated Genome Browser: free software for distribution and exploration of genome-scale datasets. Bioinformatics. 2009, 25 (20): 2730-2731. 10.1093/bioinformatics/btp472.PubMed CentralView ArticlePubMedGoogle Scholar
- Wagner EGH, Vogel J: Approaches to Identify Novel Non-messenger RNAs in Bacteria and to Investigate their Biological Functions: Functional Analysis of Identified Non-mRNAs. Handbook of RNA Biochemistry. Edited by: Hartmann RK, Bindereif A, Schön A, Westhof E. 2005, Weinheim; [Great Britain]: Wiley-VCH, 2:Google Scholar