Conflicting selection pressures on synonymous codon use in yeast suggest selection on mRNA secondary structures

BMC Evolutionary Biology20088:224

DOI: 10.1186/1471-2148-8-224

Received: 02 October 2007

Accepted: 31 July 2008

Published: 31 July 2008



Eukaryotic mRNAs often contain secondary structures in their untranslated regions that are involved in expression regulation. Whether secondary structures in the protein coding regions are of functional importance remains unclear: laboratory studies suggest stable secondary structures within the protein coding sequence interfere with translation, while several bioinformatic studies indicate stable mRNA structures are more frequent than expected.


In contrast to several studies testing for unexpected structural stabilities, I directly compare the selective constraint of sites that differ in their structural importance. I.e. for each nucleotide, I identify whether it is paired with another nucleotide, or unpaired, in the predicted secondary structure. I assume paired sites are more important for the predicted secondary structure than unpaired sites. I look at protein coding yeast sequences and use optimal codons and synonymous substitutions to test for structural constraints. As expected under selection for secondary structures, paired sites experience higher constraint than unpaired sites, i.e. significantly lower numbers of conserved optimal codons and consistently lower numbers of synonymous substitutions. This is true for structures predicted by different algorithms.


The results of this study are consistent with purifying selection on mRNA secondary structures in yeast protein coding sequences and suggest their biological importance. One should be aware, however, that accuracy of structure prediction is unknown for mRNAs and interrelated selective forces may contribute as well. Note that if selection pressures alternative to translational selection affect synonymous (and optimal) codon use, this may lead to under- or over-estimates of selective strength on optimal codon use depending on strength and direction of translational selection.


Messenger RNA (mRNA) sequences encode the amino acid sequence of the protein but may also bear additional information. For example, certain synonymous codons may improve translation [13] and a variety of motifs may regulate expression at the level of translation, cellular localization, decay or splicing [49]. Many of these motifs are secondary structures, and eukaryotic mRNAs contain regulatory structures in their 5' and 3' UTRs [1015], or introns [16, 17]. However, it remains unclear whether secondary structures in the coding regions are of functional importance. Laboratory studies suggest that local secondary structures within coding regions can interfere with translation [18, 19], and one may therefore expect selection against structures that are too stable. Surprisingly, however, several bioinformatic studies find that RNA structures within the protein coding regions are more stable than expected by chance [2023] (but see [24] for opposing result). These studies used various algorithms to predict the secondary structures of mRNA sequences, and then compared the free energy values of these structures to the values for randomized sequences.

Here, I test for selection on mRNA secondary structure using another approach. Instead of testing for unexpected structural stabilities, I directly compare the selective constraint of sites that differ in their importance for the predicted secondary structure. I.e. I predict the secondary structure of coding yeast sequences using different algorithms, and for each nucleotide, I identify whether it is paired with another nucleotide, or unpaired. I assume paired sites are more important for the predicted secondary structure than unpaired sites. If there is selection for secondary structures, one might expect higher structural constraint at paired than at unpaired sites. Such constraint would affect synonymous codon use and substitution rates. In S. cerevisiae a relationship between codon use, tRNA abundance and expression level indicates that codon use is affected by selection for translationally optimal codons [1]. If there is selection for mRNA structure, structurally important sites may be under conflicting selection pressures: a codon might support the preferred mRNA structure that is translationally non-optimal. Under structural selection, one might expect lower numbers of optimal codons at paired than at unpaired sites. If mRNA structure is conserved across species, one might further expect lower numbers of synonymous substitutions at paired than at unpaired sites; possible compensatory substitutions however may make the latter test predictions less clear-cut. When a mutation occurs at a paired site and disrupts the pairing ability, a second compensatory mutation on the corresponding paired site may restore the pairing ability [25, 26]. Compensatory mutations may increase substitution numbers at paired sites. Innan and Stephan [27] show however, that unless selection against deleterious intermediates is very small, substitutions should occur only very slowly in paired regions [27].

Accurate structure prediction is obviously crucial for these tests. In several studies [2022], mRNA structures are predicted by thermodynamic properties using the minimum free energy (MFE) algorithm [28] only although taking the whole ensemble of possible structures and comparative information into account is known to increase predictive accuracy [2932]. I therefore predict the secondary structures by thermodynamic and comparative information (RNA- and ALIfold [33]), using the minimum free energy (MFE) algorithm and McCaskill's partition function of the thermodynamic equilibrium [34].

Results of this study are consistent with selection upon mRNA structures: numbers of conserved optimal codons and synonymous substitutions are reduced at structurally important sites.


Choice of study organism & data

I focus on Saccharomyces cerevisiae, as this model eukaryote is well studied, with genome sequences available for it and several related species. Importantly yeast allows using optimal codon numbers to investigate alternative selective constraints while controlling for effects of base composition. This is because (i) translational selection has been investigated extensively and supported in yeast [13]: certain translationally "optimal" codons increase in frequency with expression level and correspond to the most abundant tRNAs in the cell or to the tRNA with which they form the strongest binding. (ii) Crucially, translationally optimal codons in yeast are not biased towards GC-ending codons, as in many other Eukaryotic organisms. In yeast 12 optimal codons end with G or C (-GC), 12 with A or T (-AT). To control for base composition is important as RNA secondary structure predictions are – at least partly- based on thermodynamic properties and will therefore be affected by GC content: GC nucleotides form the most stable binding with three hydrogen bonds and will consequently more likely be paired in the structure. From the yeast alignments provided by Kellis et al. [35] comparing Saccharomyces cerevisiae with S. paradoxus, S. mikatae and S. bayanus, I use 492 genes that have start and stop codons but no premature stop codons or frame-shifting indels in all four species.

Secondary structure

I predict the secondary structure of the coding sequences using the below methods and identify for each nucleotide whether it is paired with another nucleotide, or unpaired. I assume paired sites are more important for the predicted secondary structure than unpaired sites. Note however, that unpaired sites may well be important for maintaining the mRNA's tertiary structure.

Secondary structure prediction methods

The thermodynamic stability of a secondary structure is measured as the amount of free energy released or used by forming base pairs. Positive free energy requires work to form a structure, negative free energy releases stored work. Free energy parameters are estimated from chemical melting experiments. The widely used Minimum Free Energy (MFE) algorithm [28] computes the one single structure with the most negative energy value, that thermodynamically is hence the most likely to be formed. The MFE algorithm seems fairly accurate for short RNA sequences, for which ~73% of paired sites are accurately predicted. mRNAs however are likely to be present in a population of structures [36, 37]. Often 5–10% of structures share very similar free energy values [38], and the predicted MFE structure might just be one out of many thermodynamically similar structures. Taking all possible secondary structures of the thermodynamic equilibrium into account, McCaskill's algorithm [34] computes the most probable structure and calculates the probability that each site is paired. When taking base pairings with high probabilities, the accuracy of the prediction increases [29]. Another benefit of McCaskill's algorithm is that it is less affected by small but reasonable variations in the underlying energy parameters – while the MFE prediction is very sensitive [39, 40]. I used the RNAfold (Vienna RNA Secondary Structure [33, 41]) package to predict structures of the four yeasts separately using the MFE and McCaskill's algorithms. When using McCaskill's algorithm, I consider sites to be paired that pair with high probability (>2/3) across the structure ensemble; all other sites are considered as unpaired. With increasing sequence lengths predictive accuracy decreases presumably because of the enormous increase in the number of potential base pairings that can be made as sequence length increases [42]. I therefore look at both the complete set of genes, and at the subset of genes shorter than 800 bp.

To predict the secondary structure, one can also assume structural conservation, and compute the one consensus structure that allows the largest amount of structural conservation across homologous sequences. Especially supportive of structural conservation are sites that vary at the sequence level but retain potential of Watson-Crick pairings in the structure (co-variations). Structures predicted with the aid of comparative data appear to be more accurate than those based on thermodynamic properties alone [3032]. I use the ALIfold package [33, 43] that integrates comparative information in the prediction made with either MFE or McCaskill's algorithm and predict the consensus structures of the four yeasts together using the ALIfold default settings for co-variation weight (Φ1 = 1, and Φ2 = 1).

Optimal codon use

Codon identification is based on the S. cerevisiae sequence. Optimal codons are defined as in Kliman et al. (2003) [44]. The relative frequency of optimal codons (Fop[45]) is the ratio of optimal codons to synonymous codons. I compute the relative frequency of optimal codons for each amino acid and gene separately. For amino acids with both one AT- as well as one GC-ending optimal codon (thr, val, ile, ser), I compute the relative optimal codon frequencies of the two optimal codons per amino acid separately. Throughout the paper, the terms "optimal" and "suboptimal" will refer to translational selection.


If there is selection for secondary structures, one may expect higher constraint at structurally important (paired) than at structurally less important (unpaired) sites.

(1) Under translational selection one may expect lower numbers of translationally optimal codons at paired compared to unpaired sites. Note that the analysis is restricted to those codons that are conserved across the four yeast species and are likely to experience stronger selection pressures. Restricting the analysis to conserved sites is crucial for the ALIfold measure, as it incorporates substitutions in its prediction: ALIfold may tend to pair conserved sites, and under translational selection conserved sites tend to have higher optimal codon use than non-optimal sites. This could generate an artificial positive correlation between optimal codon numbers and structure when considering all codons. As GC-ending optimal codons are more likely to be paired, I look at GC- and AT-ending optimal codons separately. I do this for the four yeast species separately (using RNAfold) as well as for their consensus structure (using ALIfold), using MFE as well as McCaskill's algorithm for both methods.

(2) If mRNA structures are conserved across species one may further expect lower numbers of substitutions at paired compared to unpaired sites. As ALIfold incorporates comparative information, this test is only meaningful for structures predicted by RNAfold. Codons experiencing non-synonymous substitutions are excluded from this analysis as one may expect possible selection on mRNA structure will mainly affect synonymous substitutions, while non-synonymous substitutions will be more constrained for other reasons. To check the structural similarity and potential conservation of predicted structures across species, I first compute the relative number of base pairings per gene that are consistently, i.e. unambiguously, predicted to be paired or unpaired across species. To estimate structural constraint at synonymous sites, I count for each synonymous optimal and non-optimal codon how often the respective third codon position is paired and unpaired in the S. cerevisiae structure (RNAfold) and how often the codon is conserved or experiences a synonymous substitution compared to S. paravensis. Note that translational selection and structural selection may be counter-balancing with respect to synonymous substitution numbers. I.e. unpaired sites with high numbers of optimal codons may experience reduced synonymous substitution numbers due to translational selection while paired sites with high numbers of non-optimal codons may experience reduced synonymous substitution numbers due to structural selection. To disentangle structural selection from translational selection, I look at optimal and non-optimal codons separately as. I further look at GC- and AT-ending codons separately as mutational processes and gene conversion events may be compositionally biased [46].


Each of our analyses generates a set of 2 × 2 contingency tables per gene and per amino acid or codon. These are divided according to whether the site is paired or unpaired in the predicted secondary structure, and whether (1) the codon is optimal or non-optimal, and whether (2) the codon is conserved or synonymous polymorphic across the four species. To combine these independent 2 × 2 tables, I use the Mantel-Haenszel Z statistic according to Sokal and Rohlf [47]. I compute joint probabilities for all tables or certain subsets. To disentangle an effect of GC content on synonymous codon use at paired sites, I combine amino acids with AT-ending ending and amino acids with GC-ending optimal codons. I exclude contingency tables when expected values were zero, tested for homogeneity and computed the joint odds ratio (WMH) and its significance, including the continuity correction. I orient the odds ratio such that selection in favour of mRNA secondary structure is indicated by WMH <1: i.e. lower numbers of optimal codons, and lower numbers of synonymous substitutions at paired sites.


1) Conserved optimal codon numbers are significantly lower at paired compared to unpaired sites, irrespective of the method (RNA- and ALIfold) and algorithm (MFE and McCaskill's) used to predict secondary structure. Crucially, the tendency remains whether I consider amino acids with AT- or GC-ending optimal codons (Tables 1, 2).
Table 1

Comparison of conserved optimal codon numbers at paired and unpaired sites.






GC leu & lys

Genes shorter than 800 bp

RNAfold S. cerevisiae


0.646 ***

0.624 **

0.503 ***

1.353 ***



0.542 ***

0.567 ***

0.422 ***

0.910 ***

RNAfold S. paravensis


0.667 ***

0.608 ***

0.542 NS

1.137 ***



0.560 ***

0.571 ***

0.407 **

1.119 ***

RNAfold S. mikitae


0.653 ***

0.590 ***

0.544 *

1.172 ***



0.572 ***

0.623 ***

0.411 ***

1.131 ***

RNAfold S. bayanus


0.640 ***

0.597 **

0.502 ***

1.181 ***



0.537 **

0.557 ***

0.401 ***

1.003 ***



0.638 ***

0.577 ***

0.465 NS

1.499 ***



0.468 **

0.444 ***

0.326 ***

0.997 ***

All genes

RNAfold S. cerevisiae


0.920 ***

0.863 ***

0.751 ***

1.584 ***



0.841 ***

0.866 ***

0.676 ***

1.436 ***

RNAfold S. paravensis


0.878 ***

0.826 ***

0.733 ***

1.497 ***



0.819 ***

0.822 ***

0.659 ***

1.460 ***

RNAfold S. mikitae


0.912 ***

0.814 ***

0.790 NS

1.160 ***



0.839 ***

0.869 ***

0.740 NS

1.404 ***

RNAfold S. bayanus


0.904 ***

0.855 ***

0.742 ***

1.590 ***



0.833 ***

0.840 ***

0.675 ***

1.455 ***



0.899 ***

0.759 **

0.739 NS

1.937 ***



0.770 ***

0.724 ***

0.645 ***

1.529 ***

I combine contingency tables for all amino acids and genes (ALL) and subsets of amino acids with GC- and AT- ending optimal codons (leu and lys are treated separately, as these two GC-ending amino acids behave very opposing, see below Table 2). Mantel Haenzsel estimators and significances are presented, WMH <1 = lower optimal codon use at paired than at unpaired sites.

* < 0.05, ** < 0.01, *** < 0.005, NS = not significant. Structure prediction is based on ALIfold and RNAfold using MFE and McCaskill's (Mc) algorithm.

Table 2

Comparison of conserved optimal codon numbers at paired and unpaired sites.


RNAfold (S. cerevisiae)







Amino acids with GC-ending optimal codons


1.380 ***

1.273 ***

1.811 ***

1.456 ***


1.772 ***

1.654 ***

2.019 ***

1.654 ***


1.075 ***

1.120 ***

0.786 NS

0.864 **


0.809 NS

0.720 NS

0.741 NS

0.631 NS


0.707 NS

0.600 NS

0.568 NS

0.657 *


0.813 NS

0.838 NS

0.656 ***

0.746 **


0.891 *

0.831 NS

0.725 **

0.637 NS

Amino acids with one GC- and one AT-ending optimal codon


1.318 ***

1.254 ***

1.555 ***

1.049 ***


1.272 ***

0.978 ***

1.561 ***

1.009 ***


0.585 NS

0.852 ***

0.734 ***

0.560 ***


0.705 *

0.794 NS

0.806 *

0.686 NS


0.838 **

0.921 ***

0.811 ***

0.750 ***


0.940 ***

0.876 ***

0.932 ***

0.655 ***


0.614 NS

0.668 NS

0.583 NS

0.591 NS


0.794 NS

0.816 NS

0.816 ***

0.985 ***

Amino acids with AT-ending optimal codons


0.968 **

0.904 ***

1.073 ***

0.883 ***


0.588 ***

0.545 ***

0.518 ***

0.460 ***


0.894 NS

0.779 NS

0.901 ***

0.752 NS


0.442 ***

0.353 ***

0.293 ***

0.278 ***


0.946 ***

0.350 ***

0.386 ***

0.315 ***


0.851 NS

0.734 NS

0.729 **

0.689 *


0.500 NS

0.382 NS

0.755 ***

0.403 NS

Separately for each amino acid, I combine contingency tables of the different genes. Mantel Haenzsel estimators and significances are presented, with WMH <1 = lower optimal codon use at paired than at unpaired sites. Structure prediction is based on ALIfold and RNAfold using MFE and McCaskill's (Mc) algorithm.

* < 0.05, ** < 0.01, *** < 0.005, NS = not significant

The tendency is true for most amino acids separately; even GC-ending optimal codons that are more likely to be paired for thermodynamic reasons tend to be less frequent at paired sites (Table 2). Notable exceptions however are leu, lys, ile, and for RNAfold additionally phe (Table 2). One explanation for these exceptions may be that selection strength for translationally optimal codons is stronger in these amino acids, for example translational errors may be more likely or more costly. Considering prediction accuracy may decrease with gene length, I first restricted the data to genes shorter than 800 bp; including all genes however does not change the result.

(2) I first check the similarity and potential conservation of structures predicted by RNAfold. The major parts of mRNAs do not seem conserved in structure across species or prediction accuracy is low: 75% of sites are ambiguous, i.e. predicted to be paired in one or more species, but predicted to be unpaired in the remaining species (Table 3). When looking pairwise on average 41% of sites are ambiguous; number of ambiguous sites is only slightly lower for short genes. The ambiguity of predicted structural status will introduce considerable noise and may cause non-significant results. Despite high ambiguity in structure prediction the numbers of synonymous substitutions are consistently lower (WMH<1) at paired sites (Table 4). The tendencies remain when restricting the data to genes shorter than 800 bp. Other species comparisons lead to similar results (data not presented). Results become significant for GC-ending (optimal and non-optimal) codons (when structure is predicted using McCaskill's algorithm). G and C nucleotides do not only form stronger bonds and are more likely to be paired and structurally important, they are also more likely to be unambiguously predicted paired than A and T nucleotides (means GC: 0.263, AT: 0.153, t = 29.8409, df = 866.05, ***). This could reduce the level of noise and cause the significance of results for GC-ending codons.
Table 3

Similarity of predicted structures for species pairs.

Species comparison

Prediction Method




Genes shorter than 800 bp

Across all 4 yeasts


27% ± 0.6

17% ± 0.3

10% ± 0.5



27% ± 0.4

9% ± 0.2

18% ± 0.3

S. cerevisiae – S. paravensis


63% ± 0.6

36% ± 0.3

27% ± 0.3



64% ± 0.5

22% ± 0.4

42% ± 0.5

S. cerevisiae – S. mikitae


59% ± 0.3

34% ± 0.2

25% ± 0.2



61% ± 0.4

20% ± 0.4

41% ± 0.6

S. cerevisiae – S. bayanus


42% ± 0.4

23% ± 0.5

19% ± 0.1



59% ± 0.3

20% ± 0.0

40% ± 0.5

All genes

Across all 4 yeasts


25% ± 0.3

8% ± 0.1

16% ± 0.2



25% ± 0.4

16% ± 0.2

9% ± 0.1

S. cerevisiae – S. paravensis


61% ± 0.4

36% ± 0.2

25% ± 0.2



63% ± 0.3

22% ± 0.3

41% ± 0.4

S. cerevisiae – S. mikitae


58% ± 0.2

35% ± 0.2

23% ± 0.1



60% ± 0.3

20% ± 0.3

40% ± 0.4

S. cerevisiae – S. bayanus


57% ± 0.2

34% ± 0.1

23% ± 0.1



58% ± 0.2

20% ± 0.2

39% ± 0.4

The average percentages of sites (± variances) unambiguously predicted to be paired (P/all) and/or unpaired ((P+U)/all, U/all) for the respective species comparison using RNAfold MFE and McCaskill's (Mc) algorithm are presented.

Table 4

Comparison of synonymous substitution numbers at paired and unpaired sites.






































































Looking at S. cerevisiae and S. paravensis, I compare numbers of each codon in S. cerevisiae being either synonymous non-conserved or conserved at paired or unpaired sites. Structure prediction is based on RNAfold upon the S. cerevisiae sequence using MFE and McCaskill's (Mc) algorithm. Mantel Haenzsel estimators and significances are presented. WMH<1 = lower numbers of synonymous substitutions at paired sites. (1) All genes, (2) Genes that are shorter than 800 bp.

* < 0.05, ** < 0.01, *** < 0.005, NS = not significant


I tested for evidence of selective constraint acting on mRNA secondary structures in protein coding yeast genes. Predicted secondary structures differ greatly according to the prediction method used and between species. Nevertheless, there are significantly fewer conserved optimal codons and consistently fewer synonymous substitutions at paired sites for all predicted secondary structures. The results of this study are consistent with purifying selection on mRNA secondary structures.

Similar tendencies of codon use have been reported for Drosophila and humans: mRNA stability seems high when optimal codon use is low in Drosophila [48] and paired sites contain an excess of rare codons in humans [49]. Note that in this study, the comparison of optimal codon use is restricted to conserved sites. Besides the methodological need for ALIfold (see Material and Methods), the restriction to conserved sites restricts the analysis to sites potentially under considerable strong selection. For RNAfold structures, results become non-significant when not restricting the data to these conserved sites (data not presented). Strong conflicting selection pressures seem to act on certain sites while the remaining sites seem less constrained for structure. Selection on local and not global structures may explain these results and contribute to the low structural similarity across species. Selection on local mRNA structures in coding regions of eukaryotic genes has been suggested before [49]. Beside the low structural similarity also compensatory substitutions may contribute to the non-significant results when comparing substitution numbers at paired and unpaired sites.

Previous bioinformatic studies that focussed on whether or not the thermodynamic stability of mRNA structures of various organisms is selected for or against [2023, 49, 55] lead to partly inconsistent results and controversies about the accurate randomization procedure. In these studies, the observed MFE is compared to the expected MFE, which is estimated by taking the mean MFE of randomized versions of the same sequence, and a significant deviation is taken as evidence for selection for or against thermodynamic stability of the structure. The randomization of sequences can be performed in a number of different ways holding various properties of the sequence constant, while randomizing others. The properties are of biological importance; variables that are affected by forces other than selection for mRNA structure – for example the amino acid sequence – should be fixed. Which variables should remain free to vary however may not always be obvious, while the results are very sensitive to them. Di-nucleotide content for example might be selected for its effect on stability and should be allowed to vary for randomized sequences argue Chamary and Hurst [22]. However, di-nucleotides might well be affected by mutation bias, or selected for some other reason [21], in which case, di-nucleotide content should be kept fixed. The control of di-nucleotides in fact renders significant results non-significant [2023, 55].

Note that in contrast to comparing observed and expected MFE values, the comparison of constraint at paired and unpaired sites does not indicate that selection acts for or against the thermodynamic stability of the structure, but that the very predicted structure is under selection. With respect to selection for or against stability of structures, ALIfold results indicate that the thermodynamically most stable global structure is not conserved across the four yeasts: ALIfold consensus energy value is much higher i.e. less stable compared to the average energy value of the single sequences [see also Washietl et al. [50] for approach]. This is conform with results of Babak et al. [24] which support selection against stability of structures in coding regions. It is reasonable to expect selection on mRNA structures may act against too stable structures because too stable and un-flexible mRNA structures may interfere for instance with translation [18] and some mRNAs flexibility may allow their specific and dynamic complexes with other factors. mRNAs lead a complex life [51] and besides thermodynamic stability, selection on mRNA structure may also exist to maintain specific local or global mRNA structures that allow binding and interaction with other factors and thus effect biological functioning. Not only structural targets may be of effect, also accessibility of sequence targets may depend on global or local mRNA structures.

While results of this study are consistent with selection upon mRNA structures in coding regions and support laboratory studies that report synonymous substitutions are functionally important with respect to mRNA structure and translation in humans [5254], two considerations should be made. First, we do not know whether thermodynamic mRNA structure predictions predict the mRNA structures that are formed in the cell. mRNAs are generally associated with other factors [51], and effects of mRNA-associated microRNAs and proteins on the structure are hard to predict. Also, kinetics of mRNA folding and pseudo-knots are not considered here. Even with the comparative method, mRNA structures may remain at best approximations of the real mRNA structures in the cell. Secondly, the predicted and also the real structure will be affected by certain DNA patterns – however whether or not the respective DNA patterns are selected for mRNA structure or another reason may be hard to judge. There are several DNA patterns one may consider. (i) Di-nucleotide content of naturally occurring sequences leads to higher than expected thermodynamic stability [e.g. [21, 23, 55]]. Di-nucleotide content may be selected for its effect on mRNA structure but it may also be affected by mutation bias, or selected for some other reason, for example for nucleosome positioning [5658] or transcription pause sites [59]. (ii) Frequency of polypurine tracts is increased in exons and may affect thermodynamic structure. Again, polypurine tracts may be selected with respect to mRNA structure but also for other reasons such as enhancing splicing [60]. (iii) Translational protein folding into alpha-helix and beta sheets may affect synonymous codon use [61] and periodic DNA patterns may affect mRNA structure. If thermodynamic predictions correspond to any other force such as selection on nucleosome positioning and transcription or co-translational pause sites, the observed patterns may be a consequence of that and inference of selection acting directly upon on mRNA structure may be incorrect.

Alternative selection upon mRNA structures (or any other selective target) may counterbalance translational selection and explain why the bias towards translationally optimal codon is never complete and even in highly expressed genes non-optimal codons are used. Alternative selection may also contribute to the discrepancy between expected and observed codon bias [58], and may lead to systematic underestimates of selection strength for optimal codons. As selection for mRNA structures may be acting stronger on GC-ending codons, in organisms in which potential translationally optimal codons are biased towards GC, such as Drosophila, mammals, C. elegans, estimates of selective strength for optimal codons may also be overestimated. It will be worth considering effects of alternative selection and disentangling the different targets of selection.


I tested for evidence of selective constraint acting on mRNA secondary structures in yeast. Predicted structures differ greatly according to the prediction method used and between species. Nevertheless, there are significantly fewer conserved optimal codons and consistently fewer synonymous substitutions at paired sites for all predicted secondary structures. These results are consistent with purifying selection on mRNA secondary structures in protein coding yeast sequences and suggest their biological importance. One should consider however that accuracy of structure prediction is unknown for mRNAs and interrelated selective forces may contribute. Selective pressures alternative to translational selection seem affect synonymous and optimal codon use in yeast. Depending on strength and direction of translational selection in an organism, such alternative selective forces may lead to under- or over-estimates of selective strength on optimal codon use.



I first would like to thank Adam Eyre-Walker very much; he has provided plenty valuable advice, many helpful comments and suggestions at various stages of the manuscript; I am very grateful for his continuous presence for discussions and encouragement. I thank John Welch very much for an early helpful discussion and encouragement, as for many helpful comments later on the manuscript. I thank Wolfgang Stephan for friendly welcome me in his RNA discussion group, including Andreas Brunnert, Mihaela Martis and Sebastian Strempel who are all acknowledged for helpful discussions on RNA prediction methods and Perl codes. I am very grateful for support from Joachim Hermisson and his respect for my independent PhD projects. Finally, I thank both anonymous reviewers for heir helpful comments. NS was partly supported by an Emmy-Noether grant to J. Hermisson.

Authors’ Affiliations

Ludwig-Maximilan Universität, Biocenter
Centre for the Study of Evolution, School of Life Sciences, University of Sussex


  1. Ikemura T: Review Codon Usage and tRNA Content in Unicellular and Multicellular Organisms. Mol Biol Evol 1985, 2: 13–34.PubMed
  2. Kanaya S, Yamada Y, Kinouchi M, Kudo Y, Ikemura T: Codon Usage and tRNA genes in eukarotes: correlation of codon usage diversity with translational efficiency and with CG–dinucleotide usage as assessed by multivariate analysis. J Mol Evol 2001, 53 (4–5) : 290–298.View ArticlePubMed
  3. Percudani R, Pavesi A, Ottonello S: Transfer RNA gene redundancy and translational selection in Saccharomyces cerevisiae . J Mol Biol 1997, 268: 322–330.View ArticlePubMed
  4. Baker KE, Parker R: Nonsense–mediated mRNA decay: terminating erroneous gene expression. Cur Opinion Cell Biol 2004, 16: 293–299.View Article
  5. Buratti E, Baralle FE: Influence of RNA secondary structure on the Pre–mRNA splicing process. Mol Cell Biol 2004, 24: 10505–10514.View ArticlePubMed
  6. Gonsalvez GB, Urbinati CR, Long RM: RNA localization in yeast: moving towards a mechanism. Biol Cell 2005, 97: 75–86.View ArticlePubMed
  7. Kozak M: Review Regulation of translation via mRNA structure in prokaryotes and eukaryotes. Gene 2005, 361: 13–37.View ArticlePubMed
  8. StJohnston D: Moving messages: the intracellular localization of mRNAs. Nature Reviews Molecular Cell Biology 2005, 6: 363–375.View Article
  9. Tekotte H, Davis I: Intracellular mRNA localization: motors move messages. Trends in Genetics 2002, 18: 636–642.View ArticlePubMed
  10. Babendure J, Babendure JL, Ding J–H, Tsien RY: Control of mammalian translation by mRNA structure near caps. RNA 2006, 12: 1–11.View Article
  11. Pelletier J, Sonnenberg N: The involvment of mRNA secondary structure in protein synthesis. Biochem Cell Biol 1987, 65: 576–581.View ArticlePubMed
  12. Batey RT: Structures of regulatory elements in mRNAs. Cur Opinion Struct Biol 2006, 16: 299–306.View Article
  13. Ringner M, Krogh M: Folding Free Energies of 5'–UTRs Impact Post–Transcriptional Regulation on a Genomic Scale in Yeast. PLOS Comp Biol 2005, 1, 7: e72.
  14. Mignone F, Gissi C, LIuni S, Pesole G: Untranslated regions of mRNAs. Genome Biol 2002, 3: Review S0004.View Article
  15. Parsch J, Stephan W, Tanda S: A Highly Conserved Sequence in the 3'–Untranslated Region of the Drosophila Adh Gene Plays a Functional Role in Adh Expression. Genetics 1999, 151: 667–674.PubMed
  16. Kirby DA, Muse SV, Stephan W: Maintenance of pre–mRNA secondary structure by epistatic selection. Proc Natl Acad Sci USA 1995, 92: 9047–9051.View ArticlePubMed
  17. Chen Y, Stephan W: Compensatory evolution of a precursor messenger RNA secondary structure in the Drosophila melanogaster Adh gene. Proc Natl Acad Sci USA 2003, 100: 11499–11504.View ArticlePubMed
  18. Doma MK, Parker R: Endonucleotic cleavage of eukaryotic mRNAs with stalls in translation elongation. Nature 2006, 440: 561–564.View ArticlePubMed
  19. Liebhaber SA, Cash F, Eshleman SS: Translation inhibition by an mRNA coding region secondary structure is determined by its proximity to the AUG initiation codon. J Mol Biol 1992, 226: 609–621.View ArticlePubMed
  20. Seffens W, Digby D: mRNAs have greater negative folding free energies than shuffled or codon choice randomized sequences. Nucleic Acids Res 1999, 27: 1578–1584.View ArticlePubMed
  21. Katz L, Burge CB: Widespread selection for Local RNA secondary Structure in Coding Regions of Bacterial Genes. Genome Research 2003, 13: 2042–2051.View ArticlePubMed
  22. Chamary JV, Hurst LD: Evidence for selection on synonymous mutations affecting stability of mRNA secondary structure in mammals. Genome Biol 2005, 6: R75.View ArticlePubMed
  23. Shabalina SA, Ogurtsov AY, Spiridonov NA: A periodic pattern of mRNA secondary structure created by the genetic code. Nucleic Acids Res 2006, 34: 2428–2437.View ArticlePubMed
  24. Babak T, Blencowe BJ, Hughes TR: Considerations in the identification of functional RNA structural elements in genomic alignments. BMC Bioinformatics 2007, 8: 33.View ArticlePubMed
  25. Kimura M: The role of compensatory neutral mutations in molecular evolution. J Genet 1985, 64: 7–19.View Article
  26. Chen Y, Carlini DB, Baines JF, Parsch J, Braverman J, Tanda S, Stephan W: RNA secondary structure and compensatory evolution. Genes Genet Syst 1999, 74: 271–286.View ArticlePubMed
  27. Innan H, Stephan W: Selection Intensity Against Deleterious Mutations in RNA Secondary Structures and Rate of Compensatory Nucleotide Substitutions. Genetics 2001, 159: 389–399.PubMed
  28. Zuker M, Stiegler P: Optimal computer folding of large RNA sequences using thermodynamic and auxiliary information. Nucleic Acids Res 1981, 9: 133–148.View ArticlePubMed
  29. Mathews DH: Using an RNA secondary structure partition function to determine confidence in base pairs predicted by free energy minimization. RNA 2004, 10 (8) : 1178–1190.View ArticlePubMed
  30. Gutell RR, JC Lee, JJ Cannone: The accuracy of ribosomal RNA comparative structure models. Cur Opinion Struct Biol 2002, 12: 301–310.View Article
  31. Gardner PP, Griegerich R: A comprehensive comparison of comparative RNA structure prediction methods. BMC Bioinformatics 2004, 5: 140–158.View ArticlePubMed
  32. Mathews DH, Turner DH: Prediction of RNA secondary structure by free energy minimization. Cur Opinion Struct Biol 2006, 260: 270–278.View Article
  33. Hofacker IL: Vienna RNA secondary structure server. Nucleic Acids Res 2003, 31: 3429–3431.View ArticlePubMed
  34. McCaskill JS: The Equilibrium Partition Function and Base Pair Binding Probabilities for RNA Secondary Structure. Biopolymers 1990, 29: 1105–1119.View ArticlePubMed
  35. Kellis M, Patterson N, Endrizzi M, Lander ES: Sequencing and comparison of yeast species to identify genes and regulatory elements. Nature 2003, 423: 241–254.View ArticlePubMed
  36. Christofferson RE, McSwiegen JA, Konings D: Application of computational technologies to ribozyme biotechnology products. J Mol Struct THEOCHEM 1994, 311: 273–284.View Article
  37. Betts L, Spremulli LL: Analysis of the role of the Shine–Dalgarno sequence and mRNA secondary structure on the efficiency of translational initiation in the Euglena gracilis chloroplast atpH mRNA. J Biol Chem 1994, 269: 26456–26463.PubMed
  38. Eddy SR: How do RNA folding algorithms work? Nature Biotechnology 2004, 22: 1457–1458.View ArticlePubMed
  39. Zuker M: Calculating nucleic acid secondary structure. Cur Opinion Struct Biol 2000, 10: 303–310.View Article
  40. Layton DM, Bundschuh R: A statistical analysis of RNA folding algorithms through thermodynamic parameter perturbation. Nucleic Acids Res 2005, 33: 519–524.View ArticlePubMed
  41. Hofacker IL, Fontana W, Stadler PF, Bonhoeffer LS, Tacker M, Schuster P: Fast folding and comparison of RNA secondary structures. Monatshefte Chemie 1994, 125: 167–188.View Article
  42. Mathews DH, Sabina J, Zuker M, Turner DH: Expanded sequence dependence of thermodynamic parameters provides robust prediction of RNA secondary structure. J Mol Biol 1999, 288: 911–940.View ArticlePubMed
  43. Hofacker IL, Fakete M, Stadler PF: Secondary Structure Prediction for Aligned RNA Sequences. J Mol Biol 2002, 319: 1059–66.View ArticlePubMed
  44. Kliman RM, Naheelah I, Santiago M: Selection conflicts, gene expression, and codon usage trends in yeast. J Mol Evol 2003, 57: 98–109.View ArticlePubMed
  45. Ikemura T: Correlation between the abundance of Escherichia coli transfer RNAs and the occurrence of the respective codons in its protein genes: a proposal for a synonymous codon choice that is optimal for the E. coli translational system. J Mol Biol 1981, 151: 389–409.View ArticlePubMed
  46. Birdsell JA: Integrating Genomics, Bioinformatics, and Classical Genetics to Study the Effects of Recombination on Genome Evolution. Mol Biol Evol 2002, 19: 1181–1197.PubMed
  47. Sokal RR, Rohlf JF: Biometry. 3 Edition Freeman and Company. New York 1995.
  48. Carlini DB, Chen Y, Stephan W: The relationship between third–codon position nucleotide content, codon bias, mRNA secondary structure and gene expression in the Drosophilid Alcohol Dehydrogenase genes Adh and Adhr . Genetics 2001, 159: 623–633.PubMed
  49. Meyer IM, Miklos I: Statistical evidence for conserved local secondary structure in the coding regions of eukaryotic mRNAs and pre–mRNAs. Nucleic Acids Res 2005, 33: 6338–6348.View ArticlePubMed
  50. Washietl S, Hofacker IL, Stadler PF: Fast and reliable prediction of non–coding RNAs. Proc Natl Acad Sci 2005, 107: 2454–2459.View Article
  51. Moore MJ: From Birth to Death: The complex lifes of eukaryotic mRNAs. Science 2005., 309:
  52. Shen LX, Basilion JP, Stanton VP Jr: Single–nucleotide polymorphisms can cause different structural folds of mRNA. Proc Natl Acad Sci 1999, 96: 7871–7876.View ArticlePubMed
  53. Duan J, Wainwright MS, Comeron JM, Saitou N, Sanders AR, Gelernter J, Gejman PV: Synonymous mutation in the human dopamine receptor D2 (DrD2) affect mRNA stability and synthesis of the receptor. Hum Mol Genetics 2003, 12: 205–216.View Article
  54. Nackley AG, Shabalina SA, Tichileva IE, Satterfield K, Korchynskyi O, Makarov SS, Maixner W, Diatchenko L: Human Catechol–O–Methyltransferase Haplotypes modulate protein expression by altering secondary structure. Science 2006, 314: 1930–1933.View ArticlePubMed
  55. Workman C, Krogh A: No evidence that mRNAs have lower free energies than random sequences with the same di–nucleotide distribution. Nucleic Acids Res 1999, 27: 4816–4822.View ArticlePubMed
  56. Segal E, Fondufe–Mittendorf Y, Chen L, Thastrom AC, Field Y, Moore IK, Wang J–PZ, Widom J: A genomic code for nucleosome positioning. Nature 2006, 442: 772–778.View ArticlePubMed
  57. Ioshikhes IP, Albert I, Zanton SJ, Pugh F: Nucleosome positions predicted through comparative genomics. Nat Genet 2006, 38: 1104–5.View Article
  58. Ioshikhes I, Bolshoy A, Derenshteyn K, Borodovsky M, Trifonov EN: Nucleosome DNA Sequence Pattern Revealed by Multiple Alignment of Experimentally Mapped Sequences. J Mol Biol 1996, 262: 129–139.View ArticlePubMed
  59. Herbert KM, La Porta A, Wong BJ, Mooney RA, Neuman KC, Landick R, Block SM: Sequence–Resolved Detection of Pausing by Single RNA Polymerase Molecules. Cell 2006, 125: 1083–1094.View ArticlePubMed
  60. Paz A, Mester D, Nevo E, Karol A: Looking for Organization Patterns of Highly Extressed Genes: Purine–Pyrimidine Composition of Precursor mRNAs. J Mol Evol 2007, 64: 248–260.View ArticlePubMed
  61. Kimichi–Sarfaty C, Oh MM, Kim I–W, Sauna ZE, Calcangno AM, Ambudkar SV: A "Silent" Polymorphism in the MDR1 Gene changes substrate specificity. Science 2007, 315: 525–528.View Article
  62. Bulmer M: The selection–mutation–drift theory of synonymous codon usage. Genetics 1991, 129: 897–907.PubMed


© Stoletzki. 2008

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 cited.