In this work, we have analyzed several translation parameters and their relevance to the decay of ΔermC mRNA, a constitutively transcribed, 254-nt mRNA. The factors that were considered were the RBS, start codon, extent of translation of the ΔermC ORF, and location of the RBS or start codon relative to the 5′ end. Although it would be difficult to extrapolate from a study of ΔermC mRNA alone to B. subtilis mRNAs in general, we believe that intensive study of a representative mRNA species can be informative as to the primary mechanism of decay initiation.
The 8.2-min half-life of Δ
ermC mRNA is at the high end expected for bacterial mRNAs. In a recent global analysis of the half-lives of mRNAs from
E.
coli grown in minimal medium, it was found that 80% of individual mRNA half-lives were between 2 and 8 min, with a mean half-life of 5.7 min (
5). The Δ
ermC mRNA RBS consists of a 7-nt sequence complementary to the 3′-end sequence of 16S rRNA (Fig. ) and is surrounded by AU-rich sequences on both sides (Fig. , upper portion). The absence of G and C residues around the RBS is thought to be optimal for ribosome binding (
21). In addition, the small size of Δ
ermC mRNA may offer fewer internal decay initiation sites. Although we believe that the initiation of decay in
B.
subtilis is controlled primarily by access to the 5′ end (see below), alternative mechanisms operating through direct internal entry are also likely to occur. Such alternative mechanisms would explain the decreased half-life of Δ
ermC-
lacZ mRNA (6.3 min, as opposed to 8.2 min for Δ
ermC mRNA). Nevertheless, the fact that Δ
ermC-
lacZ mRNA is about 12 times the size of Δ
ermC mRNA but showed only a 23% decrease in half-life demonstrates the overriding importance of the 5′ end. Furthermore, the fact that there was no difference in the half-lives of translated and untranslated Δ
ermC mRNA and Δ
ermC-
lacZ mRNA suggests that any minor degradation pathway is not affected by ribosome flow.
Our study goes well beyond two earlier reports about the stability of
B.
subtilis mRNAs with regard to translational signals. In one report (
25), the 20-min half-life of
gsiB mRNA, encoding a stress response protein, was shown to be dependent on a strong RBS. Changing 3 nt of the
gsiB RBS resulted in a fourfold decrease in the half-life. The 5′ end of
gsiB is at −25; thus, the results obtained with
gsiB mRNA are very similar to our results, in which the RBS mutation resulted in a fourfold decrease in the half-life. Other translational effects on stability were not determined for
gsiB mRNA. Figure shows that the
gsiB RBS is predicted to have a higher free energy for binding 16S rRNA than the Δ
ermC RBS. Whether the difference in the free energy of RBS interactions between Δ
ermC and
gsiB (−8.6 and −9.5 Kcal/mol, respectively) accounts for the large difference in the mRNA half-lives (8 and 20 min) is a matter of speculation. This notion could be tested by substituting the
gsiB RBS sequence with that of Δ
ermC and vice versa.
The other report (
22), which appeared as our studies of Δ
ermC mRNA were in progress, looked at the stability of
aprE mRNA, which has a half-life of ≥25 min.
aprE mRNA also showed a fourfold decrease in stability when one or two residues of its RBS were mutated. Another result similar to ours was that the insertion of a stop codon after
aprE codon 4 did not affect stability, demonstrating that translation of the coding sequence was not needed for wild-type stability. On the other hand, mutating the
aprE start codon resulted in only a small decrease in half-life, from ≥25 min to 18 min, a finding which is very different from our finding that mutating the Δ
ermC mRNA start codon had as great a negative effect on the half-life as mutating the RBS. This difference can be explained by the fact that the
aprE start codon is GUG, which we have shown gives less stability than the AUG start codon when present in Δ
ermC mRNA (Fig. ). Thus, the
aprE start codon may not contribute as much to the relative stability of
aprE mRNA as the AUG start codon does to Δ
ermC mRNA. More importantly, the
aprE mRNA has a 5′-terminal secondary structure, which is very important for stability; weakening of this structure resulted in a fivefold decrease in stability (
22). It is likely that the stability of
aprE mRNA is determined primarily by this structure, since the strength of the RBS interaction, the GUG start codon, and the suboptimal spacing between the RBS and the start codon (Fig. ) all suggest weak ribosome binding. Thus, the negative effect of interfering with ternary complex formation on
aprE mRNA may be mitigated by the presence of a 5′-terminal secondary structure. Analysis of the 5′ untranslated region of Δ
ermC mRNA with the GCG FOLD program (Wisconsin Package version 10.0; Genetics Computer Group, Madison, Wis.) predicted no secondary structure for this sequence.
The recent report of Yusupova et al. on the path of mRNA through the ribosome (
43) is germane to thinking about mRNA protection achieved by ribosome binding. The ribosome was shown to protect a region of mRNA from −15 to +16 (the A residue of the AUG start codon is at +1), similar to much earlier estimates of −16 to +16, which were based on the protection of ribosome-bound mRNA from endonuclease digestion (
37). According to these data, we can estimate that the binding of a ribosome in a ternary complex on Δ
ermC mRNA would leave about 16 or 17 nt unprotected at the 5′ end. Based on earlier studies of RNase E-dependent decay in
E.
coli, one might expect that a 5′-terminal, single-stranded stretch of 16 or 17 nt would render the mRNA susceptible to rapid decay. In
E.
coli studies, it was demonstrated that, although 5′-terminal secondary structures could function as 5′ stabilizers, single-stranded extensions of 4 nt or more upstream of the structure decreased RNA stability (
6,
16). To explain why the predicted 5′-terminal, single-stranded sequence of Δ
ermC mRNA does not lead to instability, one could propose that, although the 5′-terminal 16 or 17 nt may not be in the mRNA tunnel of the ribosome, they could be in contact with the surface of the ribosome, thus limiting accessibility to a decay-initiating RNase. A model for our results, based on this view of ribosome binding and mRNA stability, is presented in Fig. . According to this model, the mRNA half-life is determined primarily by access of the 5′ end to a 5′-end-dependent endoribonuclease. A flow of ribosomes that bind near the 5′ end protects the mRNA from RNase binding. Formation of the ternary complex alone is sufficient for the ribosome to compete with the RNase; the protective effect of a bound ternary complex would be evident even in the complete absence of downstream translation. Earlier conclusions regarding the ability of erythromycin-induced ribosome stalling to stabilize downstream RNA sequences (
4,
11) can thus be viewed as a special case of the protective effect of ribosome binding.
An alternative hypothesis for explaining why an unbound 5′-terminal stretch might not render the mRNA susceptible to rapid decay has to do with the nature of the putative RNase E-like endoribonuclease that is thought to exist in
B.
subtilis (
10). From comparisons of 5′-end-dependent decay of mRNAs in
E.
coli and
B.
subtilis, it has been proposed that the RNase E-like RNase of
B.
subtilis may have a different mechanism for translocating from the 5′ end to an internal cleavage site (
9,
24). While
E.
coli RNase E can loop around a 5′-proximal bound ribosome or an RNA secondary structure to cleave the mRNA at a downstream site, the putative, functionally equivalent
B.
subtilis RNase may track along the mRNA to its cleavage site. Thus, the initiation of decay is hindered by a 5′-proximal ribosome to a greater degree in
B.
subtilis. Testing of such models awaits the identification of a decay-initiating enzyme in
B.
subtilis.