Skip to main content
  • Research article
  • Open access
  • Published:

Codon usage and modular interactions between messenger RNA coding regions and small RNAs in Escherichia coli

Abstract

Background

Small RNAs (sRNAs) are key regulators of gene expression in bacteria. In addition to modulating translation initiation, sRNAs can interact with mRNA coding regions to regulate mRNA stability and translation efficiency, enhancing or impeding progression of the ribosome along the mRNA. Since most amino acids are decoded by more than one codon (synonymous) we asked as to whether there is a codon bias in the interaction of sRNAs with coding regions of mRNAs. Therefore, we explored whether there are differences in codon usage or tRNA availability according to whether an mRNA is regulated by sRNAs or not. We also explored these parameters in the coding interaction regions in mRNAs. We focused our analysis on sRNAs that regulate multiple mRNAs.

Results

We found differences in codon adaptation index and tRNA adaptation index between sRNA-regulated and non-sRNA-regulated mRNAs. Interestingly, the sRNA-mRNA interacting regions tended to be enriched in unpreferred codons decoded by scarce tRNAs. We also found that sRNAs with multiple targets often contained modular segments capable of recognizing conserved motifs among these mRNAs.

Conclusions

Our results show that sRNAs in E. coli tend to recognize mRNA coding regions in which the ribosome is predicted to advance at low speeds. Identified motifs in interacting regions are conserved among mRNAs that are recognized by the same sRNA.

Background

Small non-coding RNAs (sRNAs) are crucial regulators of gene expression in bacteria. These 50- to 500-nucleotide-long RNAs are encoded in the genomes of all known bacteria and play major roles in many cellular processes, particularly in response to stress conditions [1, 2]. The major role of sRNAs is regulation of translation initiation, achieved by interfering with interactions between ribosomes and their binding sites at the 5’untranslated region (UTR) of the target mRNA and/or the very first codons of the open reading frame. sRNAs can also modulate mRNA stability and processing. These functions are exerted through imperfect base pairing with the target mRNA, initiated by the interaction of a few seed nucleotides followed by other base pairings. The RNA chaperone Hfq plays a major role in facilitating and stabilizing many of these interactions between these two RNA partners [3]. This Sm-like, homohexameric protein interacts with RNAs primarily through recognition of U-rich sequences and RNA motifs in the sRNA and target mRNA, respectively. Hfq-assisted recruitment of RNAse E to the base-pairing region results in degradation of the mRNA. Conversely, other interactions between an sRNA and its target block the action of RNAses, thereby stabilizing the mRNA.

The role of sRNAs in regulatory networks, particularly in modulation of transcription, has gained increasing recognition in recent years. How transcription and translation of the genetic information are coordinated with the participation of sRNAs in response to cellular cues is a topic of intense research [4]. It has been reported that a single mRNA can be the target of multiple sRNAs with different regulatory modes. This discovery shed critical light on the conformation of regulatory networks, particularly given that some mRNA targets may serve as regulatory molecules and transcription factors. The finding that one mRNA might be the target of several sRNAs also suggests that cross-talk among various cues may influence the fate of an mRNA and therefore translation of the encoded protein. Understanding the global effects of each sRNA might contribute to elucidating the general mechanisms of cellular response to changing environmental conditions.

Predictive and experimental approaches have been used to identify the mRNA targets of sRNAs. A number of computational tools to detect the complementary regions of sRNAs and mRNAs have been developed, including tools that make use of energy as well as phylogenetic criteria to discriminate between true and false predictions [5,6,7,8]. Experimental global approaches have also been developed to identify Hfq-mediated sRNA-mRNA pairings [9]. Most of the experimentally-reported target sites for sRNAs are located in the 5’UTR of the mRNAs, which is consistent with the well-established role of sRNAs in regulating translation initiation. However, an increasing number of interactions between sRNAs and mRNA coding regions have been reported recently as well [9, 10], including experimentally-confirmed examples. Most reported cases involve complexes consisting of sRNA-mRNA duplexes and Hfq, indicating that this chaperone plays a major role in modulating interactions involving mRNA coding regions.

As interactions between sRNAs and target mRNAs require base complementarity, interactions with mRNA coding regions may be influenced not only by sRNA-mRNA sequence complementarity but also by the synonymous codons encoding the amino acids of the protein.

Therefore, codon usage bias may add another level of restriction to this interaction. Most amino acids are encoded by more than one codon (often 2, 4, or 6 codons); however, these codons are not evenly distributed throughout the genome, the chromosome, or even the gene. This non-random distribution represents the so-called codon usage bias [11]. The basis of this bias at different structural levels is not yet understood. At the gene level, codon bias may be implicated in translation efficiency and therefore protein folding [12]. Therefore, both the amino acids encoded and the codon bias in the interaction region of the mRNA may influence pairings with sRNAs.

To assess this hypothesis, we developed a global prediction of mRNA-sRNAs pairings occurring in mRNA coding regions in E. coli. Codon usage and predicted availability of tRNA were determined. Special emphasis was placed on sRNAs that interact with multiple mRNAs, as conservation of interaction regions is expected, constraining the locations and characteristics of interactions potentially involved in network conformation.

Methods

sRNA and mRNA target sequences from E. coli K12 MG1655

The identities and sequences of the sRNAs described for E. coli K12 MG1655 were downloaded from the Bacterial Small Regulatory RNA Database (BSRD database) [13]. The identities of the target mRNAs reported for these sRNAs in E. coli K12 MG1655 were also obtained from the BSRD database. The sequences of the coding regions, 5’UTR (untranslated region), and 3’UTR were downloaded from EcoGene 3.0 (https://www.re3data.org/repository/r3d100010546) [14]. The UTRs for each gene were defined as the 50 nucleotides upstream or downstream from the start or stop codon, respectively. Targets of sRNAs not assigned in the BSRD database were predicted using RNApredator [6].

Identification and prediction of sRNA-mRNA interacting regions

The sequences of empirically-determined sRNA-mRNA interaction regions were obtained from sRNATarBase 3.0 [15]. Uncharacterized interaction regions were predicted using IntaRNA software [16]. For all interactions, the ΔG of binding was predicted using IntaRNA software.

Determination of CAI and tAI

The codon adaptation index (CAI) was calculated for each coding sequence using JEmboss [17]. To calculate the CAI for the region within the mRNA coding sequence that interacts with the sRNA, the region was identified using IntaRNA, and the minimal set of extra nucleotides required to conserve the original open reading frame was extracted. These sequences were then analyzed using JEmboss. To calculate adaption to the tRNA pool, we calculated the tRNA adaptation index (tAI) using a Python script [18]. To calculate tAI for the regions within the mRNA coding sequence that interact with sRNA, we extracted the sequences as described above and analyzed the sequences using a tAI calculation script.

Analyses of RIL-seq data

Information about the region of interaction and genes involved in the Hfq complex was obtained from the Additional file 1 from Melamed et al. report [9]. Coordinates and sense from interaction regions reported in this file were used to extract the sequences from E. coli MG1655 genome with a script written in Python language. Then to analyze only the region that encodes for complete codons, the sequences were analyzed using BLASTX against the protein sequences encoded in E. coli MG1655 genome to obtain the coordinates of coding sequences. This information was used to re-edit the sequences taking into account that all sequences must have a size that is a multiple of 3. CAI and tAI values were calculated according to the methods described above.

Interaction density maps for sRNAs with multiple targets

For each sRNA that recognized more than four mRNA targets, we developed a density map of interactions with mRNA targets along the sRNA sequence. In brief, each sRNA nucleotide was assigned a value of 1 when it interacted with an mRNA and a value of 0 when it did not. The sums of all positive interactions for each sRNA position were used to develop interaction profiles, in which a higher total score indicated that the nucleotide interacted with a greater number of mRNAs, and a total score of zero meant that the nucleotide did not interact with any mRNA.

Identification of conserved mRNA regions that interact with the same sRNA

Conserved regions among the different mRNAs that recognized a common sRNA were identified using the interaction density maps. In brief, the density maps indicated which mRNAs primarily interacted with a single sRNA subzone or with two or more subzones. The mRNA sequences that interacted with each subzone were determined using IntaRna software and extracted manually. The groups of sequences that recognized a single sRNA zone or multiple subzones were analyzed using MEME Package to identify conserved motifs (http://meme-suite.org/tools/meme) [19]. We used a p-value of < 0.05 as cutoff value for a statistically-significant motif.

Network identification

The sRNA-mRNA networks present in E. coli K12 MG1655 were identified and visualized using Cytoscape software [20]. sRNAs and mRNAs were represented as nodes. Edges connecting sRNA and mRNA nodes represented interactions between sRNAs and mRNAs. For networks involving an sRNA with multiple target mRNAs, edges connecting the sRNA zone and subzone to other subzones or the target mRNA also represented interactions.

Results

Identification of interacting regions in mRNAs regulated by sRNAs

Current information about sRNA and their targets are deposited in databases such as BSRD (http://www.bac-srna.org/BSRD/index.jsp) [13]. According to BSRD, a total of 108 sRNAs are reported in E. coli K-12 MG1655, and 96 of which are known to interact with 304 mRNAs, resulting in 412 sRNA-mRNA interactions. Our analysis of this database identified three types of sRNA: a) sRNAs with no identified targets, b) sRNAs with identified targets, but no characterized sRNA-mRNA interactions, and c) sRNAs with identified target(s) and characterized sRNA-mRNA interaction(s). Approximately 26.0% of the interactions (113) involved complexes in which the mRNA-sRNA interaction region had been localized, while 74.0% (191) involved interactions in which only the existence of sRNA-mRNA pair was known. These interactions were classified as “characterized” or “uncharacterized,” respectively. A total of 64 sRNAs had uncharacterized and 32 had characterized interactions. As only 26% of the interactions were characterized at the level of the sRNA-mRNA, we chose to characterize the remaining interactions. We predicted the binding regions for all uncharacterized interactions using IntaRNA software [16] (Additional file 2). Interactions were then classified according to the location of the binding region in the target mRNA. Empirically-documented interactions were most commonly found in the 5’UTR or adjacent to the 5’ end of the coding DNA sequence (CDS). Although some predicted interactions also occurred in these regions, most were found to occur in the middle of the coding regions (Fig. 1). This observation is consistent with the RNA-seq characterization of Hfq-dependent sRNA-mRNA interactions described by Melamed et al. [9, 21], in which approximately 60% of all described interactions occur in the middle of the coding regions.

Fig. 1
figure 1

Status of sRNA-mRNA interaction data in E. coli K12 MG1655. a Schematic representation of the arbitrary subdivision of mRNAs used in this work. 5’-UTR, coding DNA sequence (CDS-5’, middle (CDS-M), CDS-3’), 3’-UTR). b Upper panel: sRNAs archived in BSRB were classified according to the status of knowledge regarding their interactions with target mRNAs: characterized interaction (Char), uncharacterized interaction (Non-char), or unknown (sRNA) (Fig. 1B). Lower panel: interactions found in BSRD were also classified as empirical or predicted according to the method used for their determination. c Comparison of the distributions of predicted (this work) and characterized (BSDR) interactions of the various segments of the mRNAs. d Distributions of empirical and predicted binding energies of sRNA-mRNA interactions. Nt indicates for nucleotides

To verify that the predicted interactions represented plausible models, we compared the predicted and empirically-determined ΔG of interaction. The results showed that the binding energy of both types of interactions were within the same range (Fig. 1d).

Modular interactions of sRNAs with multiple mRNA targets

To understand how sRNAs interact with multiple targets, we identified sRNAs with more than four targets and mapped the nucleotides that interact with the targets for each case (Fig. 2). This analysis revealed that a total of 16 sRNAs interacted with at least four mRNAs (Table 1). The density of interactions (measured as the number of interactions per nucleotide) along the sRNA sequence allowed for identification of zones and subzones that interact with defined groups of mRNAs. This analysis also showed that some mRNAs interacted with more than one subzone (Fig. 2c). Most sRNAs contained at least two zones, one or more interacting with mRNAs and another that had a high probability of forming a stem-loop structure (Fig. 2b). Only RyhB was found to interact with mRNAs along its entire sequence. RygD, RprA, SgrS, DsrA, OxyS, and GcvB showed more than one zone with a high probability of forming a stem-loop secondary structure that would not interact with any mRNA (Table 1, Additional file 1). Analyses focused on interactions within the coding regions produced similar profiles (Additional files 3 and 4). To determine whether the different zones and subzones of each sRNA have a preferred interaction region, we localized the interactions within the target mRNA sequences. Our results showed that each sRNA was capable of interacting with different regions on the target mRNAs (5’UTR, 5’UTR/ -CDS, CDS-5, CDS-M, 3’CDS CDS-3, and 3’UTR); however, sRNA zones and subzones located at the ends of the sequence (5’ or 3’) tended to recognize particular regions of the mRNA. The sRNA OmrA provides an example of the first pattern of interaction. About 89% of sRNA-mRNA interactions (25/28) involving subzone OmrA-A1 occurred in the CDS-M, while 60% (6/10) of interactions involving subzone OmrA-A2 occurred in the 5’UTR of the target mRNA. The sRNAs RyhB, OmrB, GcvB, and CyaR showed a similar pattern. OxyS exemplified the opposite pattern; 67% of interactions involving OxyS-A occurred in the 5’UTR, while 70% of interactions involving OxyS-C2 occurred in the CDS of the target mRNA. Spot42, FnrS, DsrA, MicF, and RygD shared this pattern (Additional file 5). These results suggest that sRNAs with multiple targets have a modular structure, in which each subzone primarily recognizes the same type of region in their target mRNAs.

Fig. 2
figure 2

Modular structure of sRNAs with multiple targets. a The numbers of interactions for Spot42 (top panel) and OmrA (bottom panel) were plotted as a function of the position of each nucleotide in the primary sequence of the sRNA, allowing for construction of a density map that was then used to identify the zones (top line) and subzones (brackets) that interact with mRNAs. b Interaction zones and subzones were located in the predicted secondary structures of the sRNAs. The probability of forming a base pair is indicated using a color scale ranging from 0 (green) to 1 (red). c A network of the interactions between the mRNA and the sRNA interaction zones was then developed (Panel c). Blue circles connected by blue lines represent the sRNA structure. Boxes represent the mRNA, and edges represent interactions with the corresponding sRNA zone. Repression (red lines), empirical (solid line), and predicted (dotted line) interactions are indicated. The CAI values for each target are indicated using a color scale ranging from green (low CAI values) to red (high CAI values)

Table 1 Distribution of interactions of sRNAs that recognize multiple mRNA targets

sRNAs that interact with multiple targets recognize regions that are conserved among sRNA-regulated mRNAs

The results presented in the previous section suggest that sRNAs with multiple targets have a modular structure with zones and subzones that are specialized for interaction with particular mRNA regions, including the coding sequences. To determine whether the subzones identified recognized conserved regions of target mRNAs, we analyzed the sequences of the mRNAs that interacted with each subzone (Fig. 3). We found conserved sequences (statistically significant with p value < 0.05) in target mRNAs that were recognized by each of the following sRNAs: OmrA, RyhB, GcvB, OxyS, and Spot42. (Table 2). We evaluated whether these conserved sequence motifs were complementary to the cognate sRNA subzones. Our analysis shows that these motifs were complementary to the sRNA subzones, indicating that interaction with sRNA is likely the major role of these regions. This result supports the idea that sRNAs with multiple targets interact with target mRNAs through modules conserved in sequences, which allow for sRNA-mRNA interactions based on complementarity. These motifs were located mainly in coding regions (5’UTR (17%), 5’UTR-CDS (21.3%), CDS-5’ (21.3%), CDS-M (36.2%), CDS.3’ (1%), and 3’UTR (3.2%)), suggesting that constraints induced by the coding sequences influence sRNA-mRNA interactions through conserved motifs.

Fig. 3
figure 3

Strategy for identifying conserved regions in mRNAs recognized by a common sRNA. As an example, the interactions between OmrA and mRNAs containing one conserved motif (gntP and malK) or two conserved motifs (xylH) are shown. a To identify conserved mRNA regions that recognize a common sRNA, we first constructed an interaction density map. The mRNA sequence segments recognized by zones or subzones of a particular sRNA were then determined based on these density maps. These short sequences were aligned and analyzed using MEME to identify conserved motifs. Significantly conserved motifs are represented by a sequence logo. b Complimentary interactions between sRNAs and mRNAs were predicted, and the conserved motifs were mapped in the mRNA. A network approach was used to graphically represent the presence of conserved motifs in the mRNAs that interacted with the sRNA zone and subzones (nodes and edges in blue). mRNAs with conserved motifs complementary to the A1 subzone, the A2 subzone, or both are shown as red, cyan, and green nodes, respectively

Table 2 Conserved mRNA motifs that recognize modular sRNA sequences

Our analysis found 5 sRNAs (OmrA, RhyB, GcvB, Spot42, and OxyS) that interacted with significantly conserved mRNA motifs. Melamed et al. [9] previously described 16 sRNAs with target motifs, but their list did not include OmrA, Spot42, or OxyS. Despite the different approaches, we identified the same target motifs reported by Melamed et al. for GcvB and RyhB. Additionally, our study included an analysis of interactions between specific mRNA regions and sRNA subzones. In fact, splitting the analysis into sRNA subzones increased the number of motifs identified (Fig. 3, Table 2).

This analysis also showed that mRNAs that interact with two or more sRNA subzones generally do so via a single motif. In fact, the mRNAs encoding for XylH, GalK, and MetI were the only identified examples of mRNAs containing two motifs, each recognizing a specific sRNA subzone.

For example, the OmrA-A1 subzone interacted with 39 mRNAs, and 21/39 contained the 6 bases motif CUCUGG. On the other hand, the OmrA2 subzone interacted with 20 mRNAs, all of which contained the 6 bases consensus motif GCGUAC. Of the 10 mRNAs that interacted with these two subzones, xylH was the only mRNA with interactions involving both consensus motifs.

RyhB was the sRNA with the second-largest number of interactions. We identified a motif in approximately 50% of the mRNAs regulated by this sRNA. This consensus sequence was located in the RyhB-A2 subzone, where the MEME package also localized a motif in 19/33 mRNAs that interacted with this region. RyhB also shows a conserved motif in 5/8 mRNAs that interacted with the RyhB-A3 subzone. metI was the only mRNA that interacted with RyhB through the RyhB-A2 and RyhB-A 3 motifs.

For Spot42, we identified a motif in 100% of the mRNAs that interacted with the Spot42-A1 subzone. A consensus motif was identified in 71% of the mRNAs that interacted with Spot42-A3. GalK was the only gene that showed a conserved motif for interactions with the Spot42-A2 and Spot42-A3 subzones (Table 2).

Our analysis also showed that there are two classes of sRNA subzones, a “primary” subzone that contains the “seed” sequence for interaction with the conserved mRNA motif and a “secondary” subzone that stabilizes the interaction (data not shown).

Codon usage of sRNA-regulated mRNAs and adaptation to the predicted tRNA pool

Once the target regions in the predicted sRNA-regulated mRNAs were identified, we asked whether these mRNAs might show a different codon composition pattern than non-sRNA-regulated mRNAs. To perform this analysis, we compared the codon adaptation index (CAI) of both groups of mRNAs. No significant differences between the groups were found for CAI values calculated using the complete coding sequence (data not shown). However, splitting the coding sequence into 3 regions (CDS-5’, CDS-M, and CDS-3’, shown in Fig. 1a), a small, statistically significative difference between the groups (p = 0.0217) was identified in the C-terminal coding region, where the sRNA-regulated mRNAs showed slightly higher CAI values (Fig. 4a). A similar result was observed for adaptation to the tRNA pool (tAI) (Fig. 4b). Our results showed that sRNA-regulated mRNAs did not show any special codon preference pattern as compared to mRNAs not regulated by sRNAs. The absence of a significant difference could be the consequence of the reduced number of interactions analyzed. Currently, the data presented by Melamed et al. [9] represent the most extensive empirical information reported for sRNA-mRNA interactions in E. coli. Based on this information concerning the sRNA-mRNA interactions mediated by the RNA chaperone Hfq, we found that genes that are part of Hfq complexes tend to have slightly higher values of CAI (0.7186 ± 0.0015, n = 1155) compared to genes that do not form complex (0.7034 ± 0.001054, n = 3164) with a p value < 0.0001 using Mann-Whitney test. Similar results were obtained for tAI values of genes forming complexes with Hfq (0.2602 ± 0.001) compared to those that do not form complex (0.2536 ± 0.0001).

Fig. 4
figure 4

Codon preferences in sRNA-mRNA interactions. The mRNAs regulated by sRNAs were characterized according to codon usage (Codon Adaptation Index) (a) and availability of tRNAs to decode the mRNA (tRNA Adaptation Index) (b). sRNA-mRNA interaction regions (IR) located within the coding sequence were also characterized according to codon usage (c) and tRNA availability (d). Statistical significance was assessed using the Mann-Whitney test, *0.05 > p > 0.01, **0.01 > p > 0.001, ***p < 0.001

As the free energy of the sRNA-mRNA interaction depends on the length of the sequence and GC content of the interacting segments, it is possible that certain codon biases could be involved in the determination of interacting regions (IR). The CAI values for sRNA-mRNA interacting regions were compared to CAI values for the entire CDS (or other segments, as shown in Fig. 1a) in sRNA-regulated mRNAs. We found that the interacting regions tended to use the most-abundant codons (higher CAI values) (Fig. 4c). Although small differences were observed between the interacting regions and the entire coding region, statistically-significant differences were observed only for the first 150 (CDS-5’) or the last 30 nucleotides (CDS-3’) of the coding region. This result could reflect the previously-described translational ramp, which prevents the ribosome stalling [22]. When adaptation to tRNA availability was predicted, (Fig. 4d), the results showed that the interacting regions had lower tAI values vs. the coding sequence as a whole and vs. other segments (according to Fig. 1a). These results suggest that sRNAs tend to interact with mRNA coding regions with limited availability of tRNAs to decode them. When these analyses were performed using the information reported by Melamed et al. [9], we found similar low tAI values for the interacting regions (0.2062 ± 0.001), however CAI values were lower (0.6160 ± 0.001) (Additional files 6 and 7). Since these values were obtained using a much larger sample (n = 3580), they strongly support the information provided by the BSRD platform. Lower CAI values are more consistent with low tAI values found for interacting regions using either BRSD or Melamed et al. data, reinforcing the idea that these regions have a codon bias deviated to unpreferred codons.

Altogether, these findings suggest that ribosomes may tend to pause at the sRNA-mRNA interacting regions. Whether the ribosome protects the mRNA against degradation by RNAse E by blocking the sRNA-mRNA interaction remains unknown, but this could be a mechanism for avoiding degradation at low concentrations of sRNA [23]. Alternatively, the predicted ribosomal pauses at these sites might increase the chances of interactions between the sRNA and mRNA, inhibiting elongation of translation either by hitting an on/off switch or by enhancing mRNA degradation [24]. Whether pauses in translation are required for sRNAs to exert their regulatory actions within the mRNA coding region might be approached performing ribosome profiling experiments under the same culture conditions as Melamed et al. [9] analysis of sRNA-mRNA interactions were carried out.

Discussion

sRNAs are small, noncoding RNAs that modify translation by preventing or enhancing the interaction between the ribosome and its binding site, or by modifying the translation efficiency of the target mRNA [25,26,27]. sRNA-mediated translational repression results in rapid degradation of the target mRNA, through Hfq-dependent recruitment of RNAse E. While the translational control mechanisms of sRNAs that target the 5’UTR of mRNAs are well understood, less is known about sRNAs that target the coding regions of mRNAs. Experimental characterization of Hfq-bound mRNAs and sRNAs has shown that coding regions are important targets of sRNAs [9, 28]. sRNAs that target the coding regions of mRNAs can modify mRNA stability without altering translation efficiency, through a RNAse E-dependent mechanism [29], as well as modulating translation efficiency without promoting mRNA degradation [30, 31].

Translation efficiency also can be influenced by the frequency and disposition of synonymous codons across the coding region, determining the rate and rhythm of translation [22, 32,33,34]. As both types of regulation are exerted in the coding sequence, our hypothesis was that these two mechanisms for altering elongation of translation might be related and complementary. To address this hypothesis, we used the information available in BSRD [13], which collects the targets of sRNA predicted by the IntaRNA [16] and RNAplex [35] computational algorithms, as well as the experimentally-validated sRNAs targets in sRNATarBase [15]. To our knowledge, this is the most complete database for E. coli. However, to carry out a broad analysis of sRNA-mRNA interactions in E. coli, it was necessary to complete the uncharacterized interactions, which represented 75% of the sRNA-mRNA pairs in BSRD. We used a bioinformatics approach to identify the sequences within the mRNAs that were targeted by sRNAs. We also used for a comparative analysis, the empirical data reported by Melamed et al. [9] on the Hfq-mediated interacting regions between sRNA and mRNAs. Our results indicated that most predicted interactions occurred near the middle of the mRNA coding regions, similar to findings obtained using empirical methods in other studies (9), reinforcing our initial hypothesis. However, comparisons between the codon preferences (as measured by CAI) of sRNA-regulated and non-sRNA-regulated mRNAs showed no significant differences between the two groups, even when this analysis was constrained to the mRNAs with interaction regions located in the middle of their coding sequences. When we split the coding sequence into separate subzones, we found that codon preference and adaptation to the tRNA pool were different in each zone, increasing from the coding DNA sequence near the 5’ end towards the central portion of the coding sequence. These results are consistent with the translational ramp previously described by Tuller et al. [22]. Our results showed that interaction regions tended to be rich in codons with limited tRNA available to decode them, which could suggest that sRNAs bind to the regions where the ribosomes tend to pause. These results could be biased by the interactions of sRNAs in the five codon window [36], which are also located in the translational ramp, rich in scarce codons with limited tRNA availability. sRNA interactions at the first five codons (CDS-5’) allow for repression of translation in the absence of any identifiable interaction by complementarity at the Shine-Dalgarno sequence or start codon, the classical mechanism of sRNA-mediated translation inhibition [37]. However, most interactions were located in the central region of the coding sequence, where CAI and tAI values are higher than in the CDS-5’ region. Further experiments are required to establish whether sRNA-mRNA interactions in coding regions require a slowing of ribosome translation. Whether the coding sequence constrains the sRNA-mRNA interaction is unknown. If this constraint exists, it would have a higher impact when sRNAs interact with multiple mRNAs.

Previous studies on sRNAs that recognize multiple targets have shown that interactions occur in different regions of the sRNA and that structural flexibility of the sRNA is required for interactions with different targets [38, 39]. To analyze multiple interactions, we built density maps of the interactions. These maps allowed for identification of sRNA regions that preferentially interact with groups of mRNAs. The profiles showed that most sRNAs have a modular structure, with some regions that interact preferentially with mRNAs. Although these characteristics have been previously described for Spot42 and FnrS [40, 41], our analyses showed that this feature is a general property of sRNAs that interact with multiple mRNAs. Our results also demonstrate the presence of conserved motifs in target mRNAs. Although previous studies have also identified conserved motifs in target mRNAs, our analysis shows that most of these motifs are present in coding regions, suggesting that the coding sequence imposes a constraining factor.

Conclusions

Based on our analysis of the information available in BSRD, complemented with our own predictions of sRNA-mRNA interactions occurring within the mRNA coding sequences and the empirical data reported by Melamed et al. [9], we have shown that sRNAs preferentially recognize regions where tRNA availability is limited. As the binding of sRNAs to mRNA coding regions implies competition with the translating ribosome, this preference may have a dual effect. A reduced translation rate in these regions might induce pauses in the advancement or progression of the translating ribosome, preventing the sRNA-mRNA interaction and therefore RNAse E recruitment and mRNA degradation. Alternatively, such pauses might facilitate the sRNA-mRNA interaction, enhancing recruitment of RNase E and negatively regulating translation. These two alternative mechanisms may be strongly influenced by the availability of the aa-tRNAs that translate the interaction regions.

Abbreviations

BSRD:

Small Regulatory RNA Database

CAI:

Codon Adaptation Index

CDS:

Coding sequences

IR:

Interaction regions

mRNA:

Messenger RNA

sRNA:

Small RNA

tAI:

tRNA Adaptation Index

tRNA:

Transfer RNA

UTR:

Untranslated region

References

  1. Storz G, Vogel J, Wassarman KM. Regulation by Small RNAs in Bacteria: Expanding Frontiers. Mol Cell. 2011;43:880–91. https://doi.org/10.1016/j.molcel.2011.08.022.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. Wagner EGH, Romby P. Small RNAs in Bacteria and Archaea. Advances in genetics. 2015:133–208. https://doi.org/10.1016/bs.adgen.2015.05.001.

  3. Updegrove TB, Zhang A, Storz G. Hfq: the flexible RNA matchmaker. Curr Opin Microbiol. 2016;30:133–8. https://doi.org/10.1016/j.mib.2016.02.003.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  4. Nitzan M, Rehani R, Margalit H. Integration of Bacterial Small RNAs in Regulatory Networks. Annu Rev. Biophys. 2017;46:131–48. https://doi.org/10.1146/annurev-biophys-070816-034058.

    Article  PubMed  CAS  Google Scholar 

  5. Rehmsmeier M, Steffen P, Hochsmann M, Giegerich R. Fast and effective prediction of microRNA/target duplexes. RNA. 2004;10:1507–17. https://doi.org/10.1261/rna.5248604.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Eggenhofer F, Tafer H, Stadler PF, Hofacker IL. RNApredator: fast accessibility-based prediction of sRNA targets. Nucleic Acids Res. 2011;39 Web Server issue:W149–W154. doi:https://doi.org/10.1093/nar/gkr467.

  7. Busch A, Richter AS, Backofen R. IntaRNA: efficient prediction of bacterial sRNA targets incorporating target site accessibility and seed regions. Bioinformatics. 2008;24:2849–56. https://doi.org/10.1093/bioinformatics/btn544.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Wright PR, Richter AS, Papenfort K, Mann M, Vogel J, Hess WR, et al. Comparative genomics boosts target prediction for bacterial small RNAs. Proc Natl Acad Sci. 2013;110:E3487–96. https://doi.org/10.1073/pnas.1303248110.

    Article  PubMed  Google Scholar 

  9. Melamed S, Peer A, Faigenbaum-Romm R, Gatt YE, Reiss N, Bar A, et al. Global Mapping of Small RNA-Target Interactions in Bacteria. Mol Cell. 2016;63:884–97. https://doi.org/10.1016/j.molcel.2016.07.026.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Wroblewska Z, Olejniczak M. Hfq assists small RNAs in binding to the coding sequence of ompD mRNA and in rearranging its structure. RNA. 2016;22:979–94. https://doi.org/10.1261/rna.055251.115.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. Quax TEF, Claassens NJ, Söll D, van der Oost J. Codon Bias as a Means to Fine-Tune Gene Expression. Mol Cell. 2015;59:149–61. https://doi.org/10.1016/j.molcel.2015.05.035.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Rodnina MV. The ribosome in action: Tuning of translational efficiency and protein folding. Protein Sci. 2016;25:1390–406. https://doi.org/10.1002/pro.2950.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Li L, Huang D, Cheung MK, Nong W, Huang Q, Kwan HS. BSRD: a repository for bacterial small regulatory RNA. Nucleic Acids Res. 2013;41(Database issue):D233–8. https://doi.org/10.1093/nar/gks1264.

    Article  PubMed  CAS  Google Scholar 

  14. Zhou J, Rudd KE. EcoGene 3.0. Nucleic Acids Res. 2013;41 Database issue: D613–24. doi:https://doi.org/10.1093/nar/gks1235.

  15. Wang J, Liu T, Zhao B, Lu Q, Wang Z, Cao Y, et al. sRNATarBase 3.0: an updated database for sRNA-target interactions in bacteria. Nucleic Acids Res. 2016;44:D248–53. https://doi.org/10.1093/nar/gkv1127.

    Article  PubMed  CAS  Google Scholar 

  16. Mann M, Wright PR, Backofen R. IntaRNA 2.0: enhanced and customizable prediction of RNA–RNA interactions. Nucleic Acids Res. 2017;45:W435–9. https://doi.org/10.1093/nar/gkx279.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Carver T, Bleasby A. The design of Jemboss: a graphical user interface to EMBOSS. Bioinformatics. 2003;19:1837–43.

    Article  PubMed  CAS  Google Scholar 

  18. Tello M, Saavedra JM, Spencer E. Analysis of the use of codon pairs in the HE gene of the ISA virus shows a correlation between bias in HPR codon-pair use and mortality rates caused by the virus. Virol J. 2013; https://doi.org/10.1186/1743-422X-10-180.

  19. Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, et al. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37(Web Server):W202–8. https://doi.org/10.1093/nar/gkp335.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. Cline MS, Smoot M, Cerami E, Kuchinsky A, Landys N, Workman C, et al. Integration of biological networks and gene expression data using Cytoscape. Nat Protoc. 2007;2:2366–82. https://doi.org/10.1038/nprot.2007.324.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Melamed S, Faigenbaum-Romm R, Peer A, Reiss N, Shechter O, Bar A, et al. Mapping the small RNA interactome in bacteria using RIL-seq. Nat Protoc. 2017;13:1–33. https://doi.org/10.1038/nprot.2017.115.

    Article  PubMed  Google Scholar 

  22. Tuller T, Carmi A, Vestsigian K, Navon S, Dorfan Y, Zaborske J, et al. An evolutionarily conserved mechanism for controlling the efficiency of protein translation. Cell. 2010;141:344–54. https://doi.org/10.1016/j.cell.2010.03.031.

    Article  PubMed  CAS  Google Scholar 

  23. Bossi L, Figueroa-Bossi N. Competing endogenous RNAs: a target-centric view of small RNA regulation in bacteria. Nat Rev. Microbiol. 2016;14:775–84. https://doi.org/10.1038/nrmicro.2016.129.

    Article  PubMed  CAS  Google Scholar 

  24. Lavi-Itzkovitz A, Peterman N, Jost D, Levine E. Quantitative effect of target translation on small RNA efficacy reveals a novel mode of interaction. Nucleic Acids Res. 2014;42:12200–11. https://doi.org/10.1093/nar/gku889.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Wang J, Rennie W, Liu C, Carmack CS, Prévost K, Caron M-P, et al. Identification of bacterial sRNA regulatory targets using ribosome profiling. Nucleic Acids Res. 2015;43:10308–20. https://doi.org/10.1093/nar/gkv1158.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Duval M, Simonetti A, Caldelari I, Marzi S. Multiple ways to regulate translation initiation in bacteria: Mechanisms, regulatory circuits, dynamics. Biochimie. 2015;114:18–29. https://doi.org/10.1016/j.biochi.2015.03.007.

    Article  PubMed  CAS  Google Scholar 

  27. Morita T, Mochizuki Y, Aiba H. Translational repression is sufficient for gene silencing by bacterial small noncoding RNAs in the absence of mRNA destruction. Proc Natl Acad Sci U S A. 2006;103:4858–63. https://doi.org/10.1073/pnas.0509638103.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Holmqvist E, Wright PR, Li L, Bischler T, Barquist L, Reinhardt R, et al. Global RNA recognition patterns of post-transcriptional regulators Hfq and CsrA revealed by UV crosslinking in vivo. EMBO J. 2016;35:991–1011. https://doi.org/10.15252/embj.201593360.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Desnoyers G, Bouchard M-P, Massé E. New insights into small RNA-dependent translational regulation in prokaryotes. Trends Genet. 2013;29:92–8. https://doi.org/10.1016/j.tig.2012.10.004.

    Article  PubMed  CAS  Google Scholar 

  30. Rice JB, Vanderpool CK. The small RNA SgrS controls sugar–phosphate accumulation by regulating multiple PTS genes. Nucleic Acids Res. 2011;39:3806–19. https://doi.org/10.1093/nar/gkq1219.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Bobrovskyy M, Vanderpool CK. Diverse mechanisms of post-transcriptional repression by the small RNA regulator of glucose-phosphate stress. Mol Microbiol. 2016;99:254–73. https://doi.org/10.1111/mmi.13230.

    Article  PubMed  CAS  Google Scholar 

  32. Tuller T, Waldman YY, Kupiec M, Ruppin E. Translation efficiency is determined by both codon bias and folding energy. Proc Natl Acad Sci U S A. 2010;107:3645–50. https://doi.org/10.1073/pnas.0909910107.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Dana A, Tuller T. The effect of tRNA levels on decoding times of mRNA codons. Nucleic Acids Res. 2014;42:9171–81. https://doi.org/10.1093/nar/gku646.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Angov E. Codon usage: nature’s roadmap to expression and folding of proteins. Biotechnol J. 2011;6:650–9. https://doi.org/10.1002/biot.201000332.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. Tafer H, Hofacker IL. RNAplex: a fast tool for RNA–RNA interaction search. Bioinformatics. 2008;24:2657–63. https://doi.org/10.1093/bioinformatics/btn193.

    Article  PubMed  CAS  Google Scholar 

  36. Wroblewska Z, Olejniczak M. Contributions of the Hfq protein to translation regulation by small noncoding RNAs binding to the mRNA coding sequence. Acta Biochim Pol. 2016;63:701–7. https://doi.org/10.18388/abp.2016_1362.

    Article  PubMed  CAS  Google Scholar 

  37. Bouvier M, Sharma CM, Mika F, Nierhaus KH, Vogel J. Small RNA binding to 5’ mRNA coding region inhibits translational initiation. Mol Cell. 2008;32:827–37. https://doi.org/10.1016/j.molcel.2008.10.027.

    Article  PubMed  CAS  Google Scholar 

  38. Wu P, Liu X, Yang L, Sun Y, Gong Q, Wu J, et al. The important conformational plasticity of DsrA sRNA for adapting multiple target regulation. Nucleic Acids Res. 2017;45:9625–39. https://doi.org/10.1093/nar/gkx570.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Bækkedal C, Haugen P. The Spot 42 RNA: A regulatory small RNA with roles in the central metabolism. RNA Biol. 2015;12:1071–7. https://doi.org/10.1080/15476286.2015.1086867.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Beisel CL, Storz G. The Base-Pairing RNA Spot 42 Participates in a Multioutput Feedforward Loop to Help Enact Catabolite Repression in Escherichia coli. Mol Cell. 2011;41:286–97. https://doi.org/10.1016/j.molcel.2010.12.027.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Durand S, Storz G. Reprogramming of anaerobic metabolism by the FnrS small RNA. Mol Microbiol. 2010;75:1215–31. https://doi.org/10.1111/j.1365-2958.2010.07044.x.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

We thank to FONDECYT, Chile, grant 1150834 to OO and grant USA1799 to MT, Universidad de Chile and Universidad de Santiago de Chile for providing the logistic support for this work.

Availability of data and material

The datasets used and analysed during the current study are available from the corresponding author on reasonable request.

Funding

This work was supported by FONDECYT, Chile, grant 1150834 to OO and grant USA1799 to MT.

Author information

Authors and Affiliations

Authors

Contributions

MT contributed with conceptual and technical expertise in bioinformatics, drafting and critically revising the manuscript. FA contributed with the performance of most computational applications and data bases constructions. OO made substantial contributions to conception, design of fundamental questions, analysis and interpretation of data and drafting and review or revision of manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Omar Orellana.

Ethics declarations

Ethics approval and consent to participate

NA

Consent for publication

NA

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figures in Powerpoint format (ppt) showing density maps of the interaction of each sRNA with more than four targets. (PPTX 153 kb)

Additional file 2:

Table in Excel format containing the information about sRNA, target mRNA, sRNA subzone, site of sRNA-mRNA interaction, CAI values, tAI values, and ΔG of sRNA-mRNA interaction using the data from BRSD. (XLSX 111 kb)

Additional file 3:

Table in Word format (doc) showing sRNA zones and number of interaction in coding regions by zone and subzones for each sRNA with more than four interactions. (DOCX 14 kb)

Additional file 4:

Figures in Powerpoint format (ppt) showing density maps of the interaction for each sRNA with multiple interactions in coding region. (PPTX 107 kb)

Additional file 5:

Table in Excel format showing the number of interactions for each sRNA zones and for each site of interaction in target mRNAs. (XLSX 45 kb)

Additional file 6:

Figure in Powerpoint format (ppt) showing the average CAI (A) and tAI (B) values of sRNA-intercating regions on mRNA coding sequences obtained from BSRDatabase (IR), from the Hfq-mediated sRNA-interacting regions reported by Melamed et al. (IE_E), and CDS-5’, (CDS-M), CDS-3’, and complete (CDS) coding sequences (as in Fig. 1). Comparison of CAI and tAI values were performed using non-parametric tests against IR (*) and IR_E (*). ***/*** p < 0.001. (PPTX 168 kb)

Additional file 7:

Table in Excel format containing the information about target mRNAs, CAI and tAI values of mRNA interacting regions using the data from Melamed et al. (2016) [21]. (XLSX 538 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tello, M., Avalos, F. & Orellana, O. Codon usage and modular interactions between messenger RNA coding regions and small RNAs in Escherichia coli. BMC Genomics 19, 657 (2018). https://doi.org/10.1186/s12864-018-5038-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-018-5038-6

Keywords