Central role for RNase YbeY in Hfq-dependent and Hfq-independent small-RNA regulation in bacteria
© Pandey et al.; licensee BioMed Central Ltd. 2014
Received: 9 January 2014
Accepted: 5 February 2014
Published: 11 February 2014
Conceptual parallels exist between bacterial and eukaryotic small-RNA (sRNA) pathways, yet relatively little is known about which protein may recognize and recruit bacterial sRNAs to interact with targets. In eukaryotes, Argonaute (AGO) proteins discharge such functions. The highly conserved bacterial YbeY RNase has structural similarities to the MID domain of AGOs. A limited study had indicated that in Sinorhizobium meliloti the YbeY ortholog regulates the accumulation of sRNAs as well as the target mRNAs, raising the possibility that YbeY may play a previously unrecognized role in bacterial sRNA regulation.
We have applied a multipronged approach of loss-of-function studies, genome-wide mRNA and sRNA expression profiling, pathway analysis, target prediction, literature mining and network analysis to unravel YbeY-dependent molecular responses of E. coli exposed to hydroxyurea (HU). Loss of ybeY function, which results in a marked resistance to HU, had global affects on sRNA-mediated gene expression. Of 54 detectable E. coli sRNAs in our microarray analysis, 30 sRNAs showed a differential expression upon HU stress, of which 28 sRNAs displayed a YbeY-dependent change in expression. These included 12 Hfq-dependent and 16 Hfq-independent sRNAs. We successfully identified at least 57 experimentally inferred sRNA-mRNA relationships. Further applying a ‘context likelihood of relatedness’ algorithm, we reverse engineered the YbeY-dependent Hfq-dependent sRNA-mRNA network as well as YbeY-dependent Hfq-independent sRNA-mRNA network.
YbeY extensively modulates Hfq-dependent and independent sRNA-mRNA interactions. YbeY-dependent sRNAs have central roles in modulating cellular response to HU stress.
KeywordsYbeY Small-RNA Hfq Hydroxyurea Stress adaptation E. coli
Small regulatory RNAs (sRNAs) play key roles in modulating gene expression in both prokaryotes and eukaryotes. In bacteria, numerous sRNAs that range in size from ~50-300 nucleotides in length act on independently transcribed mRNA targets. In E. coli, approximately 80 such sRNAs have been validated . The most extensively characterized class of bacterial sRNAs are the trans-encoded sRNAs that are encoded distant from the genes for their mRNA targets and that typically have only limited complementarity (10–30 nt) . Trans- encoded sRNAs form imperfect base-pairing interactions with complementary sequences in their mRNA targets, which are often located at or near ribosome binding sites (RBS), but can also be located upstream of the translation start site as well as deep in the coding regions (CDS) . Such interactions generally result in a decrease in protein synthesis, either by stimulating mRNA degradation or by inhibiting translation. Other bacterial sRNAs, referred to as cis-encoded sRNAs, are present in close proximity to their targets, such as upstream, opposite of the 5′ UTR of the target, or between two genes in an operon . Such sRNAs have extensive (>75 nt) complementarity to their targets .
Gaining a deeper insight into how sRNAs recognize and interact with their targets is crucial to understanding the mechanism of sRNA action and function in bacteria at a molecular level. For certain sRNAs (e.g. RNAIII, RsaE, and SprD in Staphylococcus aureus as well as CyaR sRNAs), target recognition is structure-driven using C-rich stretches that are located within accessible loop regions . For many other trans-acting sRNAs, however, the interacting region is not located in a structured region, rather in a single-stranded region that is often located at the 5′ end of the sRNA. Like microRNAs (miRNAs) in eukaryotes, trans-acting bacterial sRNAs appear to recognize their targets by a seed-pairing mechanism using seeds as small as 6–7 nucleotides. Fusion studies have revealed that seed regions of RybB or MicC sRNAs are sufficient to guide the recognition of targets . Moreover, as in miRNAs, 3′ adenosine (A) residues have been reported recently to occur adjacent to the pairing region [3, 4].
In numerous bacteria, mainly Gram negatives, the RNA-binding protein Hfq is required for the action and stability of many trans-encoded sRNAs [5, 6]. The Hfq chaperone binds to A/U rich regions of sRNAs that are often located near the stem-loop structures as well as to the poly (U) regions at the 3′ end of the sRNAs . Structural and binding studies have revealed several RNA binding sites on the proximal as well as distal faces of Hfq hexameric ring, which may facilitate the interaction of sRNA and their target mRNAs . Hfq may assist duplex formation by enhancing local RNA concentrations, changing RNA structures and accelerating strand exchange and annealing. Although most of the E. coli trans-acting sRNAs that have been characterized require Hfq for base pairing, some in Vibrio cholerae do not require Hfq for pairing with target . Loss-of-function studies suggest that Hfq is essential for virulence of several pathogens as well as symbiosis of Sinorhizobium meliloti with plants [5, 9].
Hfq interacts with RNase E to serve to recruit the RNA degradation machinery once the sRNAs have base-paired with targets . Moreover, recent work has shown that a 5′-monophosphorylated sRNA seed both guides RNase E to its mRNA target and also stimulates the degradation. Hfq is needed for optimal RNase E activity in this sRNA-guided mRNA cleavage . RNase III also participates in sRNA-mediated modulation of mRNAs . Further, RNases also play an important role in generation of mature sRNAs. For example, 6S RNA maturation involves multi-layered pathway involving endonucleolytic digestion by RNase E or G and exonucleolytic trimming at 5′ and 3′ ends .
Although Hfq is known to bind to the C-terminus of RNase E and recruit it to sRNA-mediated interactions, much remains to be learned concerning the molecular mechanism and function of various RNases and other RNA binding proteins during bacterial RNA-interference. Recent progress has offered insights into bacterial sRNAs that are recognized and loaded on to the Hfq protein scaffold during their interaction with the target mRNAs [12–14] but the identity of the proteins that facilitate sRNA-mRNA interactions in bacteria lacking Hfq or facilitate the interactions of the sRNAs that are Hfq-independent remains largely unknown.
Striking conceptual parallels exist between the bacterial and the eukaryotic sRNA-pathways. In eukaryotes, numerous structural studies have revealed the molecular details of how the miRNAs and siRNAs are recognized by, and loaded onto, the Argonaute proteins and then guided to the targets by the RNA-induced silencing complex (RISC) . In contrast, considerably less is known about the molecular details of how bacterial sRNAs recognize their target mRNAs or about the roles of Hfq and other proteins in this process.
We have recently reported evidence that the highly conserved bacterial protein YbeY (SMc01113 in Sinorhizobium meliloti) may play a major, previously unrecognized role in bacterial sRNA regulation . ybeY, which is one of the 206 genes that comprise the postulated minimal bacterial genome set , is essential in some bacteria [18, 19]. In contrast, in certain other bacteria such as Escherichia coli and Sinorhizobium meliloti, ybeY is not essential but its loss sensitizes cells to a wide variety of physiologically diverse stresses and causes striking defects that affect ribosome activity, translational fidelity and ribosome assembly [20–22]. Several observations had suggested that YbeY might interact with RNA. Our mapping of 16S, 23S and 5S rRNA termini in an E. coli ΔybeY mutant showed that YbeY influences the maturation of all three rRNAs, with a particularly strong effect on maturation at both the 5′- and 3′-end of 16S rRNA as well as maturation of the 5′-termini of 23S and 5S rRNAs . Furthermore, we demonstrated that there are strong genetic interactions between ybeY and rnc (encoding RNase III), ybeY and rnr (encoding RNase R), and ybeY and pnp (encoding PNPase) .
We have recently shown that YbeY is a previously undiscovered single-strand RNase with a combination of biochemical properties that distinguish it from all previously reported RNases . Additionally, we have shown that YbeY plays in a key role in a previously unrecognized system of 70S ribosome quality control, in which YbeY and RNase R act together to degrade defective 70S ribosomes but not properly matured 70S ribosomes or individual subunits. In addition, we discovered that there is essentially no fully matured 16S rRNA in a ∆ybeY mutant at 45°C, making YbeY the first endoribonuclease to be implicated in the critically important processing of the 16S rRNA 3′ terminus.
Two key observations stimulated us to investigate the possibility of involvement of YbeY in sRNA regulation . First, we had observed that YbeY displays high structural similarities to the MID domains of Argonaute proteins. In silico modeling of substrate and protein binding suggested that YbeY has the potential to bind to sRNA seeds and we identified a possible phosphate-binding center in YbeY cavity. Second, we noted striking parallels between the phenotypes of S. meliloti hfq and smc01113 (ybeY) mutants, an observation that suggested there might be an underlying mechanistic connection. To test this hypothesis, we carried out a limited study in S. meliloti in which we evaluated the accumulation of 13 target genes and 9 sRNAs in S. meliloti compromised for their expression of Hfq or SMc01113 (YbeY). We showed that both mutants exhibited similar deregulation of sRNAs and targets . This study raised the possibility that YbeY might play a role in sRNA regulation in bacteria whose importance is comparable to that of Hfq.
To test our hypothesis that YbeY plays an important role in sRNA regulation in addition to its key roles in 70S ribosome quality control and rRNA processing, we have used an integrative biology approach to evaluate the YbeY-dependent molecular response of E. coli cells exposed to hydroxyurea (HU), a widely used inhibitor of E. coli’s class I ribonucleotide reductase. Our results indicate that YbeY plays an extremely important role in bacteria, modulating the levels of both Hfq-dependent and Hfq-independent sRNAs as well as their targets.
Design of experiment and overview of microarray analysis
Summary of microarray analysis: number of differentially regulated genes between MC4100 and Δ ybeY with or without (Unt) HU treatment
ΔybeY Unt vs MC4100 Unt
MC4100 HU vs MC4100 Unt
ΔybeY HU vs ΔybeY Unt
ΔybeY HU vs MC4100 HU
Loss of ybeY function results in changes in gene expression even in the absence of any stress, with 1151 genes significantly differentially expressed; 673 genes were up-regulated in the ΔybeY mutant relative to WT and 478 were down-regulated (Table 1, Additional file 1: Table S1). As we have previously reported , treatment of E. coli with HU causes major changes in gene expression with 1326 genes significantly differentially expressed upon HU treatment; 718 genes were significantly up-regulated upon HU treatment and 608 were significantly down-regulated (Table 1, Additional file 1: Table S1). As we have noted previously, up-regulated genes include those associated with three classes of survival responses: ribonucleotide reductases (e.g. nrdA, nrdB, nrdD, nrdE, nrdF), primosome components for replication restart (priA and priB), and the SOS response (e.g. recA, recN, sulA, umuC). We had also noted that numerous genes associated with iron import were strongly up-regulated upon treatment of WT with HU  (e.g. fepC, fepD, fepG, fhuA, fhuB, fhuC, fhuD, fhuE, fhuF, tonB, exbB, exbD). Treatment of the ΔybeY mutant with HU also resulted in major changes in gene expression with 1178 genes significantly differentially expressed: 506 genes were up-regulated upon HU treatment and 672 were down-regulated.
A comparison of the genes whose expression was altered when the WT strain was treated with HU with those whose expression was altered when the ΔybeY mutant was treated with HU revealed that 535 genes displayed a significant change in expression independent of the ybeY status of the cells (Figure 3; Additional file 2: Table S2). Genes in this category included many of those mentioned above that are potentially associated with survival when WT is treated with HU: ribonucleotide reductases (e.g. nrdA, nrdB, nrdD, nrdE, nrdF) and the SOS response (response (e.g. recA, recN, sulA, umuC). In addition, many of the genes associated with iron import were induced by HU regardless of the ybeY status of the cells (e.g. fepC, fepD, fepG, fhuE, fhuF, exbB, exbD). Thus, the presence or absence of YbeY had only subtle affects on expression of these genes. However, many of the genes whose expression levels are altered by HU treatment are strongly influenced by the ybeY status of the cells. In the case of the 1326 genes whose expression was affected when the WT strain was treated with HU, the remaining 470 up-regulated genes and 321 down-regulated genes only displayed a significant change when the cells were proficient for YbeY function. Interestingly, the primosome/replication restart genes (priA and priB) were in this category, as were certain of the iron import genes (fhuA, fhuB, fhuC, fhuD, tonB). Reciprocally, in the case of the 1178 genes whose expression was affected when the ΔybeY strain was treated with HU, the remaining 267 up-regulated genes and 376 down-regulated genes only displayed a significant change when the cells lacked YbeY function. Thus the presence or absence of YbeY function is a huge factor in determining how cells respond to HU.
YbeY dependent reprogramming of sRNA and target expression
Number of differentially regulated sRNAs between MC4100 and ΔybeY with or without (Unt) HU treatment
ΔybeY Unt vs MC4100 Unt
MC4100 HU vs MC4100 Unt
ΔybeY HU vs ΔybeY Unt
ΔybeY HU vs MC4100 HU
Of the 28 differentially regulated sRNAs (detailed in above comparisons), 12 sRNAs were Hfq-dependent and 16 sRNAs were Hfq-independent. Hfq-dependent sRNAs that also showed a YbeY dependence included OxyS, DsrA, CyaR, ArcZ, RybB, MicA, GlmZ, RyeA, RyeB, CydC, RyjB and SgrS; whereas Hfq-independent YbeY dependent sRNAs are Ffs, RygD, CsrC, RdlD, RdlA, GadY, PsrD, RygC, RyfD, RyfA, RybA, RyrC, RyeD, RydB, IsrB, and IsrC.
Insights into the biological pathways regulated by YbeY-dependent sRNAs during HU treatment in ΔybeY and WT genotypes were obtained by identifying bona-fide, experimentally inferred targets and genes that are associated with differentially regulated sRNAs using two databases, sRNATarBase and sRNAMap, and by primary literature mining ([1, 28]; Figures 2 and 4). 45 experimentally inferred genes associated with 11 Hfq-dependent (OxyS, DsrA, CyaR, ArcZ, RybB, MicA, GlmZ, RyeA, RyeB, CydC, and SgrS) and 6 Hfq-independent (Ffs, RygD, CsrC, RdlD, RdlA, and GadY) sRNAs were obtained. Again, cluster analysis revealed the differences in expression of individual genes between the four comparisons (Figure 4).
In E. coli, sRNAs can both up-regulate or down-regulate expression of their targets . Therefore, for negatively regulated targets, expression levels of sRNA-mRNA target pair should be inversely correlated i.e. for a given sRNA, if its expression was up-regulated, the level of expression of its target genes were down-regulated and vice-versa. Similarly, for positively regulated sRNA-target pairs, levels of sRNAs as well as their targets changed in the same direction. A total of 57 sRNA-mRNA interactions for 17 sRNAs were correctly correlated to change in expression of 45 experimentally inferred genes (Figure 4), indicating that these relationships are functional YbeY-dependent-sRNA-mRNA combinations during HU response.
Combining the information on known experimentally inferred targets, obtained from sRNAMap, sRNATarbase, and primary literature, as well as clustering of targets into functional groups revealed several important pathways that most likely changed in ∆ybeY cells when they were exposed to HU (Figure 4). For example pathways related to envelope stress, redox stress and oxidative stress, porin synthesis, translocation of membrane protein and their assembly, signal recognition particles (SRPs), anti-mutation response, metabolism, toxin-anti-toxin pairs, transporters, and cell survival (Figure 4) are altered in ΔybeY cells when they are exposed to HU. This analysis highlights the broad consequences of YbeY-dependent sRNA regulation on cellular physiology in response to HU stress (elaborated in the following sections).
Inferring YbeY dependent sRNA-mRNA interactions using CLR
Physiological responses of ∆ybeY cells to HU treatment and cell fate
Functional clustering of genes was performed to gain insight into the biological processes that potentially changed after HU-treatment in MC4100 and ∆ybeY genotypes (Additional file 4: Table S4, Additional file 5: Table S5, Additional file 6: Table S6, Additional file 7: Table S7, Additional file 8: Table S8 and Additional file 9: Table S9). Our analysis of transcriptional reprogramming in ∆ybeY against HU-stress suggested alteration of several pathways (Additional file 4: Table S4, Additional file 5: Table S5, Additional file 6: Table S6, Additional file 7: Table S7, Additional file 8: Table S8 and Additional file 9: Table S9) that might possibly contribute to the resistance of a ∆ybeY mutant to HU. For example, genes contributing to envelope stress as well as those responsible for synthesis and repair of cell wall, membrane, lipoproteins and polysaccharides were strongly altered in the ∆ybeY mutant. Similarly multiple two-component signal transduction systems that enable bacteria to sense, respond and adapt to environmental stresses were specifically regulated in the ∆ybeY genotype. ∆ybeY cells adjust their TCA cycle and components of electron transfer chain in response to HU stress, alterations that could in principle contribute to a reduction in the production of harmful oxidizing radicals so that the damage to genetic material may be reduced . Intriguingly, components of base excision DNA repair, which are employed by cells to repair DNA damage due to oxidizing agents, were up-regulated only in ∆ybeY cells (Additional file 8: Table S8). Furthermore, components of non-coding RNA biogenesis pathways were reprogrammed in ∆ybeY genotype (Additional file 8: Table S8). Moreover, components of drug resistance pathways were also evidently regulated only in cells with the ∆ybeY genotype (Additional file 8: Table S8). ∆ybeY cells displayed several molecular characteristics similar to those undergoing adaptation to antibiotic exposure . Exposure of WT to HU results in up-regulation of iron-uptake systems (Additional file 9: Table S9), which is highly detrimental to E. coli and since it could cause cell death during HU stress by promoting Fenton chemistry . Most of these genes of iron uptake system (e.g. tonB-exbB-exbD, fhu system genes) were expressed at lower level in HU-exposed ∆ybeY as compared to HU-exposed WT MC4100 (Additional file 7: Table S7). When the WT cells are exposed to HU treatment, genes like tonB, fhuA, B, C, D are strongly up-regulated, whereas HU-mediated up-regulation of these genes does not occur in ∆ybeY cells. These results are consistent with our earlier observation that loss of expression of tonB provides a protective effect in HU-exposed cells .
Taken together, our results indicate that the highly conserved bacterial protein YbeY plays a major role in bacterial sRNA regulation. Extending our more limited study in S. meliloti, our analyses suggest that YbeY participates in both Hfq-dependent and Hfq-independent sRNAs-mediated interactions in E. coli. This means that YbeY is a remarkably central protein in RNA metabolism in bacteria, as we have additionally shown that YbeY is also involved in 70S ribosome quality control and in rRNA processing . Thus YbeY can help a cell withstand stress both by modulating changes in gene expression through its role in sRNA regulation and by helping to maintain the fidelity of protein translation. YbeY’s three distinct RNA-related physiological roles offer an explanation for why it is one of 206 genes in the postulated minimal bacterial gene set .
We have constructed a whole-genome sRNA-target interaction network to explore organism-wide interactions for most of E. coli sRNAs and evaluated the subset of network that is reconfigured in YbeY-dependent manner during the response to HU. Our results suggest novel mechanistic insights into how cells respond to HU and reveals that YbeY and sRNAs play central roles. The complexity of E.coli’ s response to HU that is evident in our experiments is fully consistent with a recent genome-wide screening study with HU that revealed a link between non-essential ribosomal proteins and reactive oxygen species  and the subsequent demonstration that a tRNA thiolation pathway, which modulates the intracellular redox state, affects sensitivity to HU .
Several toxin-antitoxin pairs encoded by E. coli genome were also differentially reprogrammed in ∆ybeY, so was the SRP machinery (in which the Ffs sRNA is a central component) that regulates synthesis of membrane proteins and their targeting . We also observed differential expression of several sRNA modules that regulate porins and outer membrane proteins (OMPs) e.g. MicA, CyaR and RybB and their targets [(Figures 4 and 6); [33, 39]]. Reprogramming of SRP machinery along with sRNA-OMP modules of regulation could help to reduce the amount of mistranslated membrane proteins, repair damaged membrane proteins, and help circumventing envelope stress, thus enhancing survival (Figures 7 and 8). Our results suggest that oxidative and free radicals stresses pose a multilayered threat to cell and that YbeY dependent sRNAs (which include both Hfq- dependent as well as independent sRNAs) have hitherto unappreciated roles in helping cells to adapt to such stresses.
Although around 70 sRNAs have been estimated in E. coli, the targets and functions of many of them still remain unknown. We could obtain targets for only 17 differentially expressed sRNAs from literature and database sources (Figures 2 and 4). The deficiency in experimental information underscores the need to develop complementary approaches for identification of targets and functions of sRNAs as compared to standard, tedious biological and genetic protocols, and large-scale screens of mutant libraries. One such complementary approach has been recently developed in our laboratory that uses the ‘Context Likelihood of Relatedness’ (CLR) algorithm to infer the Hfq-dependent sRNA regulatory network in E. coli[26, 27]. In our current study, we have used an E. coli compendium of 759 microarrays maintained in one of our laboratories  to expand this algorithm to all the E. coli sRNAs (Hfq-dependent as well as independent). We have inferred potential targets of every sRNA with a very high q-value (FDR corrected p-value; <0.005). Our results suggest presence of a complex global regulatory network modulated by sRNAs (Figures 5 and 6) that would help the bacteria to program its response effectively to changes in their environment. It is not difficult to imagine that bacteria may elicit specific smaller sub-networks and sRNA-circuits in order to adapt to specific stresses. We indeed observed a similar situation as we were able to determine the YbeY-dependent sRNA sub-network that responded to HU (Figure 6). Our analysis is also significant as it demonstrates the relevance of integrative biology approaches to infer the global and specific regulatory circuits of sRNAs, as well as presents examples of different regulatory sRNA circuits. Several models of sRNA circuitry modules have been proposed e.g. single input module, dense overlapping regulon, positive- and negative feedback modules, and feed forward modules . All these modules were readily observed in our global network.
It is evident that sRNAs act post-transcriptionally and modulate gene expression through both extensive and limited base-pairing interactions with their targets. Several of these require the RNA chaperone Hfq for pairing with their targets, with Hfq assisting in the trans-annealing of the sRNAs to target mRNA in an antisense manner . Recent progress has offered insights into how bacterial sRNAs are recognized and loaded on to the Hfq protein scaffold during their interaction with the target mRNAs [12–14]. However, Hfq has not been identified in numerous sequenced bacteria, whereas YbeY is extremely highly conserved. Thus, our results raise the possibility that YbeY may play an especially important role in sRNA regulation in bacteria that do not encode Hfq.
The phenomenon of sRNA recognition and its guidance to target mRNA is quite well understood in higher organisms, where miRNAs and siRNAs are loaded onto Argonaute proteins and then guided and assembled on to the targets by the RISC. Structures of Argonaute proteins include a ‘MID domain’ that specifically recognize 5′-phosphate of the sRNAs and anchors them on to the Argonaute/RISC, and a ‘PIWI domain’ that have hydrolytic (RNase) capacity . Interestingly, it was our observation that YbeY has structural similarities to the MID domain of Argonaute proteins that stimulated us to consider the possibility that YbeY might play a role in bacterial sRNA regulation and to carry out modeling studies indicating that it had the potential to bind a small seed RNA.
In contrast, considerably less is known about how bacterial sRNAs are recognized to interact with their targets. Although important recent work has shown that a 5′-monophosphorylated sRNA seed both guides RNase E to its mRNA target and stimulates degradation, and that Hfq is needed for optimal RNase E activity in this sRNA-guided mRNA cleavage , many mechanistic questions remain. Furthermore, as noted above, many bacteria lack Hfq yet exhibit sRNA regulation, while our results suggest that YbeY-dependent, Hfq-independent sRNA regulation may be considerably more important in bacteria that possess Hfq than has hitherto been recognized.
It will be extremely interesting to determine how YbeY participates mechanistically in sRNA regulation. YbeY possesses a metal-dependent, single strand endoribonuclease activity that is relatively weak compared to many strictly degradative RNases . Furthermore, as its RNase activity on naked RNA substrates is not particularly specific (a preference for cleavage after U’s), yet its in vivo activities in rRNA processing are highly specific, its RNase activity within living cells must be highly controlled .
A particular biochemical characteristic of YbeY that is potentially of interest with respect to its possible roles in sRNA regulation is that YbeY can bind to and cut single stranded oligoribonucleotides as short as 10 nucleotides in length, but is not able to cut a 7 base oligoribonucleotide despite the presence of a site that is cleaved in the context of larger oligoribonucleotides . This raises the possibility that YbeY could play two possible mechanistic roles in sRNA regulation in bacteria. A possible non-catalytic role of YbeY could be that it binds a seed RNA and subsequently influences the interaction of that seed with its target mRNA or with other proteins. A possible catalytic role of YbeY could be that it participates in the degradation of the target mRNA and/or the sRNA.
A particular structural characteristic of YbeY is also of potential interest with respect to its possible roles in sRNA regulation. Both the MID domain of AGO and RNase E have a seed-binding site that is constrained at the 5′-end of the seed RNA and, in fact, both particularly recognize the 5′-phosphate of the RNA seed [11, 15]. In contrast, the RNA binding site of YbeY is an open channel . Although our modeling studies suggest that YbeY could bind a short seed RNA and even recognize a 5′-phosphate, the nature of the YbeY structure suggests that it could potentially interact with a seed sequence that is internal to the sRNA. Since some seed sequences of sRNA are internal and thus are not expected to end with a 5′-phosphate [3, 39], YbeY might possibly play a particularly important role in their regulatory action.
Taken together, our study places YbeY in the centre of sRNA-mediated gene regulation in bacterial genomes. Our study also offers mechanistic insights into regulatory basis of response of E. coli to HU stress. Along with demonstrating the role of YbeY, this study places sRNA pathway at the center of cellular response to oxidative stress caused by exposure of cells to HU.
Strains, growth conditions, treatment and microarrays
All strains were grown in Luria-Berani (LB) medium at 37°C with constant aeration. Strains were grown with or without 100 mM HU in liquid cultures for microarrays and with 10 mM HU on LB-Agar plates for survival assays . The effect of ΔybeY mutation on growth rate is very modest and does not affect the cell density that is achieved . For isolation of RNA, MC4100 and ∆ybeY were grown with or without HU as described previously . Three independent cultures were grown for isolating RNA as previously described; cDNA was prepared and microarray were performed for three independent biological replicates as described .
*.CEL files obtained from microarray hybridizations were combined with those in the ‘E. coli CDS compendium’ regularly maintained in J.J. Collins’ lab . The backbone of this compendium is the publically available M3D database at http://m3d.mssm.edu/. The raw intensities were background adjusted, log2-transformed and RMA-normalized with RMAexpress. The E. coli compendium used here comprised of a total of 759 RMA-normalized E. coli expression arrays, which are also publically available at the M3D database. Standard deviation (SD) of expression, σ, was calculated across the entire compendium for each gene.
We used the ‘z scale difference’ statistic described previously [25, 29] and defined: ΔZ = [(Xt – Xc)/σ], where Xt and Xc are the normalized gene expression values for a give gene in treatment and control arrays respectively. ΔZ values were calculated for all the four comparisons, WT MC4100 HU treated vs. untreated, ΔybeY HU treated vs. untreated, ΔybeY HU treated vs. WT MC4100 HU treated, and ΔybeY untreated vs. WT MC4100 untreated (Figure 2). Genes with a |ΔZ| score of >1 was considered significant [25, 29]. ΔZ allows the measurement of change in expression of each gene in any given comparison in ‘SD units’ in form of a z-test . The SD of the z-score standardization allows comparison of each observation from different normal distributions, and the average of zero avoids introducing aggregation distortions stemming from differences in the normal gene expression means. Thus, genes with extreme expression values will have intrinsically greater effects on the composite standard z-scores.
Functional clustering of genes, identification of sRNA and their targets, and network maps
Clustering of gene sets was performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID; http://david.abcc.ncifcrf.gov/; ) to identify potential gene pathways and key functional groups that may modulate response of bacteria to HU or knock out of ybeY (Figure 1). Curated gene sets from KEGG pathways, Swiss-Prot (SP) and Protein Information Resource (SP PIR keywords), Uniprot sequence features (UP Seq Feature), COG (clusters of orthologous groups) ontology, GO term analysis and SMART (simple modular architectural tool) that includes well-studied metabolic and signaling pathways, were used for annotation and clustering of genes into functional groups. Categories of function in the differentially expressed genes were determined for both up- and down-regulated genes.
For evaluating sRNAs, names of all the sRNAs were extracted from sRNAMap database . This list of sRNAs was mapped onto the microarray annotation file and all the information was extracted for the sRNAs that were found on E. coli microarray. The nomenclature of E. coli sRNAs is still not standardized and many sRNAs are known by alternate names e.g. 4.5S sRNA is also called Ffs sRNA. Therefore, we performed literature mining for all the remaining sRNAs that could not be found on the E. coli expression array in our first round of matching; sRNAs with alternate names were remapped to microarrays. Expression of a total of 54 sRNAs was detectable. Level of significance among four comparisons was determined as described above.
For the list of differentially regulated sRNAs in any one of the four comparisons, experimentally inferred information related to targets was extracted from sRNAMap and sRNATarbase databases [1, 28]. In parallel, we mined primary literature to extract such information about their targets and functions of all the differentially expressed sRNAs (Figure 2). Similarly, classification of sRNAs into Hfq- dependent and independent categories was based on primary literature and information in these databases. Functional categories for experimentally inferred targets as well as pathways in which these sRNAs may act were determined as described above.
In complementation to the above approach, we adapted a recently developed CLR based strategy to evaluate network of sRNA targets for all the sRNAs and extracted sub-network that responded to HU in YbeY dependent manner. CLR algorithm is based on relevance network theory, infers cellular regulatory interactions from a compendium of expression profiles (the algorithm is available at http://m3d.mssm.edu/network_inference.html). Although CLR was originally designed to identify regulatory interactions of transcription factors, in this work we adapt the algorithm to examine sRNA regulatory influences because of their role as posttranscriptional regulators of mRNA stability. An sRNA and a gene are predicted to interact if the mutual information between their expression levels is above a set threshold. Mutual information is a measure of the statistical dependence between two variables and, in contrast to correlation, does not assume linearity, continuity, or other specific properties of the dependence. CLR computes the significance of mutual information by assembling a background distribution from the mutual information scores of all other microarray probe sets in the compendium. This adaptive background correction allows the algorithm to eliminate false correlations and indirect functional influences. sRNA–gene interactions found to be significant using this procedure are represented in the network diagram as edges between nodes. Network analysis presented in this work focuses on the regulatory influences of Hfq-dependent sRNAs. At the time this work began, 27 Hfq-dependent sRNAs had been identified in E. coli. Because of the nature of Affymetrix annotation, our methodology is restricted to inferring relationships for genes associated with a Blattner ID, constraining our network to 24 Hfq-dependent sRNAs. We used a compendium of 759 Affymetrix E. coli Antisense2 arrays uniformly normalized with RMA to serve as input to the algorithm. This compendium includes arrays from the Many Microbe Microarray Database (E_coli _v3_Build_3), as well as other 235 arrays run in-house. Experiments did not include genetic or environmental perturbations specifically related to sRNAs but were generally focused on bacterial stress response. The reconstructed sRNA-mRNA regulatory network can help us to gain insights into the functional roles of these sRNAs. All computations were run in Matlab (Mathworks).
Hydroxyurea sensitivity assays, measurement of reactive oxygen species, and microscopy
For determining the effect of knocking out ybeY on the survival of E. coli, we conducted spot test assays on LB-Agar plates containing 10 mM hydroxyurea. Overnight grown WT MC4100 and ΔybeY cultures were diluted 1:1000 (OD600 of ~0.01) and grown to OD600 0.5-0.6; these were serially diluted and 5 μl of each dilution was spotted on LB-Agar-HU plates. Plates were incubated at 37°C overnight. For measurement of reactive oxygen species and microscopy, cells were grown to early exponential phase and then treated with 100 mM HU. 100 μL samples were collected hourly, centrifuged at 10,000 rpm, and resuspended in 100 μL PBS + 5 mM 3′-(p-hydroxyphenyl fluorescein (HPF). Cells were incubated in the dark for 15 minutes at room temperature, then centrifuged at 10,000 rpm and resuspended in 1X PBS for microscopy and FACS analysis. For FACS analysis, fluorescence data were collected using a Becton Dickinson FACSCalibur flow cytometer. For microscopy, images were obtained using the 100X oil-immersion objective lens.
This work was supported by NIH grants R01 CA021615 and GM31010 (to G.C.W.), and DP1 OD003644 (to J.C.C.), and P30 ES002019 (to the Massachusetts Institute of Technology Center for Environmental Sciences), as well as the Howard Hughes Medical Institute. S.P.P. acknowledges support by IISER-Kolkata startup grant. G.C.W. is an American Cancer Society Professor. S.P.P. is Max Planck-DST partner group leader. H.L. is supported by Mayo Clinic Center for Individualized Medicine. We thank Sheetal Modi for her comments on the manuscript.
- Huang HY, Chang HY, Chou CH, Tseng CP, Ho SY, Yang CD, Ju YW, Huang HD: sRNAMap: genomic maps for small non-coding RNAs, their regulators and their targets in microbial genomes. Nucleic Acids Res. 2009, 37 (Database issue): D150-D154.PubMed CentralPubMedView ArticleGoogle Scholar
- Waters LS, Storz G: Regulatory RNAs in bacteria. Cell. 2009, 136 (4): 615-628. 10.1016/j.cell.2009.01.043.PubMed CentralPubMedView ArticleGoogle Scholar
- Storz G, Vogel J, Wassarman KM: Regulation by small RNAs in bacteria: expanding frontiers. Mol Cell. 2011, 43 (6): 880-891. 10.1016/j.molcel.2011.08.022.PubMed CentralPubMedView ArticleGoogle Scholar
- Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120 (1): 15-20. 10.1016/j.cell.2004.12.035.PubMedView ArticleGoogle Scholar
- Barra-Bily L, Pandey SP, Trautwetter A, Blanco C, Walker GC: The Sinorhizobium meliloti RNA chaperone Hfq mediates symbiosis of S. meliloti and alfalfa. J Bacteriol. 2010, 192 (6): 1710-1718. 10.1128/JB.01427-09.PubMed CentralPubMedView ArticleGoogle Scholar
- Brennan RG, Link TM: Hfq structure, function and ligand binding. Curr Opin Microbiol. 2007, 10 (2): 125-133. 10.1016/j.mib.2007.03.015.PubMedView ArticleGoogle Scholar
- Link TM, Valentin-Hansen P, Brennan RG: Structure of Escherichia coli Hfq bound to polyriboadenylate RNA. Proc Natl Acad Sci U S A. 2009, 106 (46): 19292-19297. 10.1073/pnas.0908744106.PubMed CentralPubMedView ArticleGoogle Scholar
- Bardill JP, Hammer BK: Non-coding sRNAs regulate virulence in the bacterial pathogen Vibrio cholerae. RNA Biol. 2012, 9 (4): 392-401. 10.4161/rna.19975.PubMed CentralPubMedView ArticleGoogle Scholar
- Vogel J, Luisi BF: Hfq and its constellation of RNA. Nat Rev Microbiol. 2011, 9 (8): 578-589. 10.1038/nrmicro2615.PubMed CentralPubMedView ArticleGoogle Scholar
- Aiba H: Mechanism of RNA silencing by Hfq-binding small RNAs. Curr Opin Microbiol. 2007, 10 (2): 134-139. 10.1016/j.mib.2007.03.010.PubMedView ArticleGoogle Scholar
- Bandyra KJ, Said N, Pfeiffer V, Gorna MW, Vogel J, Luisi BF: The seed region of a small RNA drives the controlled destruction of the target mRNA by the endoribonuclease RNase E. Mol Cell. 2012, 47 (6): 943-953. 10.1016/j.molcel.2012.07.015.PubMed CentralPubMedView ArticleGoogle Scholar
- Sauer E, Schmidt S, Weichenrieder O: Small RNA binding to the lateral surface of Hfq hexamers and structural rearrangements upon mRNA target recognition. Proc Natl Acad Sci U S A. 2012, 109 (24): 9396-9401. 10.1073/pnas.1202521109.PubMed CentralPubMedView ArticleGoogle Scholar
- Sauer E, Weichenrieder O: Structural basis for RNA 3′-end recognition by Hfq. Proc Natl Acad Sci U S A. 2011, 108 (32): 13065-13070. 10.1073/pnas.1103420108.PubMed CentralPubMedView ArticleGoogle Scholar
- Papenfort K, Bouvier M, Mika F, Sharma CM, Vogel J: Evidence for an autonomous 5′ target recognition domain in an Hfq-associated small RNA. Proc Natl Acad Sci U S A. 2010, 107 (47): 20435-20440. 10.1073/pnas.1009784107.PubMed CentralPubMedView ArticleGoogle Scholar
- Hutvagner G, Simard MJ: Argonaute proteins: key players in RNA silencing. Nat Rev Mol Cell Biol. 2008, 9 (1): 22-32. 10.1038/nrm2321.PubMedView ArticleGoogle Scholar
- Pandey SP, Minesinger BK, Kumar J, Walker GC: A highly conserved protein of unknown function in Sinorhizobium meliloti affects sRNA regulation similar to Hfq. Nucleic Acids Res. 2011, 39 (11): 4691-4708. 10.1093/nar/gkr060.PubMed CentralPubMedView ArticleGoogle Scholar
- Gil R, Silva FJ, Pereto J, Moya A: Determination of the core of a minimal bacterial gene set. Microbiol Mol Biol Rev. 2004, 68 (3): 518-537. 10.1128/MMBR.68.3.518-537.2004.PubMed CentralPubMedView ArticleGoogle Scholar
- Kobayashi K, Ehrlich SD, Albertini A, Amati G, Andersen KK, Arnaud M, Asai K, Ashikaga S, Aymerich S, Bessieres P: Essential Bacillus subtilis genes. Proc Natl Acad Sci U S A. 2003, 100 (8): 4678-4683. 10.1073/pnas.0730515100.PubMed CentralPubMedView ArticleGoogle Scholar
- Akerley BJ, Rubin EJ, Novick VL, Amaya K, Judson N, Mekalanos JJ: A genome-scale analysis for identification of genes required for growth or survival of Haemophilus influenzae. Proc Natl Acad Sci U S A. 2002, 99 (2): 966-971. 10.1073/pnas.012602299.PubMed CentralPubMedView ArticleGoogle Scholar
- Davies BW, Walker GC: A highly conserved protein of unknown function is required by Sinorhizobium meliloti for symbiosis and environmental stress protection. J Bacteriol. 2008, 190 (3): 1118-1123. 10.1128/JB.01521-07.PubMed CentralPubMedView ArticleGoogle Scholar
- Rasouly A, Davidovich C, Ron EZ: The heat shock protein YbeY is required for optimal activity of the 30S ribosomal subunit. J Bacteriol. 2010, 192 (18): 4592-4596. 10.1128/JB.00448-10.PubMed CentralPubMedView ArticleGoogle Scholar
- Rasouly A, Schonbrun M, Shenhar Y, Ron EZ: YbeY - a heat shock protein involved in translation in Escherichia coli. J Bacteriol. 2009, 191: 2649-2655. 10.1128/JB.01663-08.PubMed CentralPubMedView ArticleGoogle Scholar
- Davies BW, Köhrer C, Jacob AI, Simmons LA, Zhu J, Aleman LM, Rajbhandhari UL, Walker GC: Role of Escherichia coli YbeY, a highly conserved protein, in rRNA processing. Mol Microbiol. 2010, 78 (2): 506-518. 10.1111/j.1365-2958.2010.07351.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Jacob AI, Kohrer C, Davies BW, Rajbhandary UL, Walker GC: Conserved bacterial RNase YbeY plays Key roles in 70S ribosome quality control and 16S rRNA maturation. Mol Cell. 2013, 49 (3): 427-438. 10.1016/j.molcel.2012.11.025.PubMed CentralPubMedView ArticleGoogle Scholar
- Davies BW, Kohanski MA, Simmons LA, Winkler JA, Collins JJ, Walker GC: Hydroxyurea induces hydroxyl radical-mediated cell death in Escherichia coli. Mol Cell. 2009, 36 (5): 845-860. 10.1016/j.molcel.2009.11.024.PubMed CentralPubMedView ArticleGoogle Scholar
- Faith JJ, Hayete B, Thaden JT, Mogno I, Wierzbowski J, Cottarel G, Kasif S, Collins JJ, Gardner TS: Large-scale mapping and validation of Escherichia coli transcriptional regulation from a compendium of expression profiles. PLoS Biol. 2007, 5 (1): e8-10.1371/journal.pbio.0050008.PubMed CentralPubMedView ArticleGoogle Scholar
- Modi SR, Camacho DM, Kohanski MA, Walker GC, Collins JJ: Functional characterization of bacterial sRNAs using a network biology approach. Proc Natl Acad Sci U S A. 2011, 108 (37): 15522-15527. 10.1073/pnas.1104318108.PubMed CentralPubMedView ArticleGoogle Scholar
- Cao Y, Wu J, Liu Q, Zhao Y, Ying X, Cha L, Wang L, Li W: sRNATarBase: a comprehensive database of bacterial sRNA targets verified by experiments. RNA. 2010, 16 (11): 2051-2057. 10.1261/rna.2193110.PubMed CentralPubMedView ArticleGoogle Scholar
- Kohanski MA, Dwyer DJ, Wierzbowski J, Cottarel G, Collins JJ: Mistranslation of membrane proteins and two-component system activation trigger antibiotic-mediated cell death. Cell. 2008, 135 (4): 679-690. 10.1016/j.cell.2008.09.038.PubMed CentralPubMedView ArticleGoogle Scholar
- Nakayashiki T, Mori H: Genome-wide screening with hydroxyurea reveals a link between nonessential ribosomal proteins and reactive oxygen species production. J Bacteriol. 2013, 195 (6): 1226-1235. 10.1128/JB.02145-12.PubMed CentralPubMedView ArticleGoogle Scholar
- Nakayashiki T, Saito N, Takeuchi R, Kadokura H, Nakahigashi K, Wanner BL, Mori H: The tRNA thiolaton pathway modulates the intracellular redox state in Escherchia coli. J Bacteriol. 2013, 195 (9): 2039-2049. 10.1128/JB.02180-12.PubMed CentralPubMedView ArticleGoogle Scholar
- Beisel CL, Storz G: Base pairing small RNAs and their roles in global regulatory networks. FEMS Microbiol Rev. 2010, 34 (5): 866-882.PubMed CentralPubMedView ArticleGoogle Scholar
- Papenfort K, Vogel J: Multiple target regulation by small noncoding RNAs rewires gene expression at the post-transcriptional level. Res Microbiol. 2009, 160 (4): 278-287. 10.1016/j.resmic.2009.03.004.PubMedView ArticleGoogle Scholar
- Henke JM, Bassler BL: Bacterial social engagements. Trends Cell Biol. 2004, 14 (11): 648-656. 10.1016/j.tcb.2004.09.012.PubMedView ArticleGoogle Scholar
- Mahoney TF, Silhavy TJ: The Cpx stress response confers resistance to some, but not all, bactericidal antibiotics. J Bacteriol. 2013, 195 (9): 1869-1874. 10.1128/JB.02197-12.PubMed CentralPubMedView ArticleGoogle Scholar
- Kint CI, Verstraeten N, Wens I, Liebens VR, Hofkens J, Versees W, Fauvart M, Michiels J: The Escherichia coli GTPase ObgE modulates hydroxyl radical levels in response to DNA replication fork arrest. FEBS J. 2012, 279 (19): 3692-3704. 10.1111/j.1742-4658.2012.08731.x.PubMedView ArticleGoogle Scholar
- Altuvia S, Weinstein-Fischer D, Zhang A, Postow L, Storz G: A small, stable RNA induced by oxidative stress: role as a pleiotropic regulator and antimutator. Cell. 1997, 90 (1): 43-53. 10.1016/S0092-8674(00)80312-8.PubMedView ArticleGoogle Scholar
- Bibi E: Early targeting events during membrane protein biogenesis in Escherichia coli. Biochim Biophys Acta. 2011, 1808 (3): 841-850. 10.1016/j.bbamem.2010.07.025.PubMedView ArticleGoogle Scholar
- Vogel J, Papenfort K: Small non-coding RNAs and the bacterial outer membrane. Curr Opin Microbiol. 2006, 9 (6): 605-611. 10.1016/j.mib.2006.10.006.PubMedView ArticleGoogle Scholar
- Faith JJ, Driscoll ME, Fusaro VA, Cosgrove EJ, Hayete B, Juhn FS, Schneider SJ, Gardner TS: Many microbe microarrays database: uniformly normalized affymetrix compendia with structured experimental metadata. Nucleic Acids Res. 2008, 36 (Database issue): D866-D870.PubMed CentralPubMedGoogle Scholar
- Cheadle C, Vawter MP, Freed WJ, Becker KG: Analysis of microarray data using Z score transformation. J Mol Diagn. 2003, 5 (2): 73-81. 10.1016/S1525-1578(10)60455-2.PubMed CentralPubMedView ArticleGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.