Codon usage and protein length-dependent feedback from translation elongation regulates translation initiation and elongation speed

Essential cellular functions require efficient production of many large proteins but synthesis of large proteins encounters many obstacles in cells. Translational control is mostly known to be regulated at the initiation step. Whether translation elongation process can feedback to regulate initiation efficiency is unclear. Codon usage bias, a universal feature of all genomes, plays an important role in determining gene expression levels. Here, we discovered that there is a conserved but codon usage-dependent genome-wide negative correlation between protein abundance and CDS length. The codon usage effects on protein expression and ribosome flux on mRNAs are influenced by CDS length; optimal codon usage preferentially promotes production of large proteins. Translation of mRNAs with long CDS and non-optimal codon usage preferentially induces phosphorylation of initiation factor eIF2α, which inhibits translation initiation efficiency. Deletion of the eIF2α kinase CPC-3 (GCN2 homolog) in Neurospora preferentially up-regulates large proteins encoded by non-optimal codons. Surprisingly, CPC-3 also inhibits translation elongation rate in a codon usage and CDS length-dependent manner, resulting in slow elongation rates for long CDS mRNAs. Together, these results revealed a codon usage and CDS length-dependent feedback mechanism from translation elongation to regulate both translation initiation and elongation kinetics.     

Translation rate modification by preferential codon usage: intragenic position effects

We present a model for calculating the protein production rate as a function of the translation rate. The model takes into account that the elongation rate along an mRNA molecule is non-uniform as a result of different tRNA availabilities for different codons. Initiation of ribosomes on an mRNA is normally the rate-limiting step in the translation process, and blocking of the initiation site can be avoided if the codons closest to this site allow fast translation by the ribosome. Hence, different selective forces may act on the choice of synonymous codons in the initiation region than elsewhere on a given mRNA. We show that the elongation rate along the whole mRNA influences the production rate of abundant proteins, whereas only the elongation rate in the initiation region is of importance for the production rate of rare proteins. We also present an analysis of the codon distribution along known mRNAs coding for abundant and rare proteins.     

Effect of codon adaptation on codon-level and gene-level translation efficiency in vivo

Background

There is a significant difference between synonymous codon usage in many organisms, and it is known that codons used more frequently generally showed efficient decoding rate. At the gene level, however, there are conflicting reports on the existence of a correlation between codon adaptation and translation efficiency, even in the same organism.

Results

To resolve this issue, we cultured Escherichia coli under conditions designed to maintain constant levels of mRNA and protein and subjected the cells to ribosome profiling (RP) and mRNA-seq analyses. We showed that the RP results correlated more closely with protein levels generated under similar culture conditions than with the mRNA abundance from the mRNA-seq. Our result indicated that RP/mRNA ratio could be used as a measure of translation efficiency at gene level. On the other hand, the RP data showed that codon-specific ribosome density at the decoding site negatively correlated with codon usage, consistent with the hypothesis that preferred codons display lower ribosome densities due to their faster decoding rate. However, highly codon-adapted genes showed higher ribosome densities at the gene level, indicating that the efficiency of translation initiation, rather than higher elongation efficiency of preferred codons, exerted a greater effect on ribosome density and thus translation efficiency.

Conclusions

These findings indicate that evolutionary pressure on highly expressed genes influenced both codon bias and translation initiation efficiency and therefore explains contradictory findings that codon usage bias correlates with translation efficiency of native genes, but not with the artificially created gene pool, which was not subjected to evolution pressure.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1115) contains supplementary material, which is available to authorized users.     

Trade-offs between tRNA abundance and mRNA secondary structure support smoothing of translation elongation rate

Translation of protein from mRNA is a complex multi-step process that occurs at a non-uniform rate. Variability in ribosome speed along an mRNA enables refinement of the proteome and plays a critical role in protein biogenesis. Detailed single protein studies have found both tRNA abundance and mRNA secondary structure as key modulators of translation elongation rate, but recent genome-wide ribosome profiling experiments have not observed significant influence of either on translation efficiency. Here we provide evidence that this results from an inherent trade-off between these factors. We find codons pairing to high-abundance tRNAs are preferentially used in regions of high secondary structure content, while codons read by significantly less abundant tRNAs are located in lowly structured regions. By considering long stretches of high and low mRNA secondary structure in Saccharomyces cerevisiae and Escherichia coli and comparing them to randomized-gene models and experimental expression data, we were able to distinguish clear selective pressures and increased protein expression for specific codon choices. The trade-off between secondary structure and tRNA-concentration based codon choice allows for compensation of their independent effects on translation, helping to smooth overall translational speed and reducing the chance of potentially detrimental points of excessively slow or fast ribosome movement.     

Efficient translation initiation dictates codon usage at gene start

The genetic code is degenerate; thus, protein evolution does not uniquely determine the coding sequence. One of the puzzles in evolutionary genetics is therefore to uncover evolutionary driving forces that result in specific codon choice. In many bacteria, the first 5–10 codons of protein‐coding genes are often codons that are less frequently used in the rest of the genome, an effect that has been argued to arise from selection for slowed early elongation to reduce ribosome traffic jams. However, genome analysis across many species has demonstrated that the region shows reduced mRNA folding consistent with pressure for efficient translation initiation. This raises the possibility that unusual codon usage is a side effect of selection for reduced mRNA structure. Here we discriminate between these two competing hypotheses, and show that in bacteria selection favours codons that reduce mRNA folding around the translation start, regardless of whether these codons are frequent or rare. Experiments confirm that primarily mRNA structure, and not codon usage, at the beginning of genes determines the translation rate.     

Ribosome profiling of mouse embryonic stem cells reveals the complexity and dynamics of mammalian proteomes

The ability to sequence genomes has far outstripped approaches for deciphering the information they encode. Here we present a suite of techniques, based on ribosome profiling (the deep sequencing of ribosome-protected mRNA fragments), to provide genome-wide maps of protein synthesis as well as a pulse-chase strategy for determining rates of translation elongation. We exploit the propensity of harringtonine to cause ribosomes to accumulate at sites of translation initiation together with a machine learning algorithm to define protein products systematically. Analysis of translation in mouse embryonic stem cells reveals thousands of strong pause sites and unannotated translation products. These include amino-terminal extensions and truncations and upstream open reading frames with regulatory potential, initiated at both AUG and non-AUG codons, whose translation changes after differentiation. We also define a class of short, polycistronic ribosome-associated coding RNAs (sprcRNAs) that encode small proteins. Our studies reveal an unanticipated complexity to mammalian proteomes.     

Isolation of Ribosome Bound Nascent Polypeptides in vitro to Identify Translational Pause Sites Along mRNA

The rate of translational elongation is non-uniform. mRNA secondary structure, codon usage and mRNA associated proteins may alter ribosome movement on the message^{for review see 1}. However, it''s now widely accepted that synonymous codon usage is the primary cause of non-uniform translational elongation rates¹. Synonymous codons are not used with identical frequency. A bias exists in the use of synonymous codons with some codons used more frequently than others². Codon bias is organism as well as tissue specific^2,3. Moreover, frequency of codon usage is directly proportional to the concentrations of cognate tRNAs⁴. Thus, a frequently used codon will have higher multitude of corresponding tRNAs, which further implies that a frequent codon will be translated faster than an infrequent one. Thus, regions on mRNA enriched in rare codons (potential pause sites) will as a rule slow down ribosome movement on the message and cause accumulation of nascent peptides of the respective sizes^5-8. These pause sites can have functional impact on the protein expression, mRNA stability and protein folding^{for review see 9}. Indeed, it was shown that alleviation of such pause sites can alter ribosome movement on mRNA and subsequently may affect the efficiency of co-translational (in vivo) protein folding^1,7,10,11. To understand the process of protein folding in vivo, in the cell, that is ultimately coupled to the process of protein synthesis it is essential to gain comprehensive insights into the impact of codon usage/tRNA content on the movement of ribosomes along mRNA during translational elongation.Here we describe a simple technique that can be used to locate major translation pause sites for a given mRNA translated in various cell-free systems^6-8. This procedure is based on isolation of nascent polypeptides accumulating on ribosomes during in vitro translation of a target mRNA. The rationale is that at low-frequency codons, the increase in the residence time of the ribosomes results in increased amounts of nascent peptides of the corresponding sizes. In vitro transcribed mRNA is used for in vitro translational reactions in the presence of radioactively labeled amino acids to allow the detection of the nascent chains. In order to isolate ribosome bound nascent polypeptide complexes the translation reaction is layered on top of 30% glycerol solution followed by centrifugation. Nascent polypeptides in polysomal pellet are further treated with ribonuclease A and resolved by SDS PAGE. This technique can be potentially used for any protein and allows analysis of ribosome movement along mRNA and the detection of the major pause sites. Additionally, this protocol can be adapted to study factors and conditions that can alter ribosome movement and thus potentially can also alter the function/conformation of the protein.     

Protein Synthesis in E. coli: Dependence of Codon-Specific Elongation on tRNA Concentration and Codon Usage

To synthesize a protein, a ribosome moves along a messenger RNA (mRNA), reads it codon by codon, and takes up the corresponding ternary complexes which consist of aminoacylated transfer RNAs (aa-tRNAs), elongation factor Tu (EF-Tu), and GTP. During this process of translation elongation, the ribosome proceeds with a codon-specific rate. Here, we present a general theoretical framework to calculate codon-specific elongation rates and error frequencies based on tRNA concentrations and codon usages. Our theory takes three important aspects of in-vivo translation elongation into account. First, non-cognate, near-cognate and cognate ternary complexes compete for the binding sites on the ribosomes. Second, the corresponding binding rates are determined by the concentrations of free ternary complexes, which must be distinguished from the total tRNA concentrations as measured in vivo. Third, for each tRNA species, the difference between total tRNA and ternary complex concentration depends on the codon usages of the corresponding cognate and near-cognate codons. Furthermore, we apply our theory to two alternative pathways for tRNA release from the ribosomal E site and show how the mechanism of tRNA release influences the concentrations of free ternary complexes and thus the codon-specific elongation rates. Using a recently introduced method to determine kinetic rates of in-vivo translation from in-vitro data, we compute elongation rates for all codons in Escherichia coli. We show that for some tRNA species only a few tRNA molecules are part of ternary complexes and, thus, available for the translating ribosomes. In addition, we find that codon-specific elongation rates strongly depend on the overall codon usage in the cell, which could be altered experimentally by overexpression of individual genes.     

Mutation and selection on the anticodon of tRNA genes in vertebrate mitochondrial genomes

The H-strand of vertebrate mitochondrial DNA is left single-stranded for hours during the slow DNA replication. This facilitates C-->U mutations on the H-strand (and consequently G-->A mutations on the L-strand) via spontaneous deamination which occurs much more frequently on single-stranded than on double-stranded DNA. For the 12 coding sequences (CDS) collinear with the L-strand, NNY synonymous codon families (where N stands for any of the four nucleotides and Y stands for either C or U) end mostly with C, and NNR and NNN codon families (where R stands for either A or G) end mostly with A. For the lone ND6 gene on the other strand, the codon bias is the opposite, with NNY codon families ending mostly with U and NNR and NNN codon families ending mostly with G. These patterns are consistent with the strand-specific mutation bias. The codon usage biased towards C-ending and A-ending in the 12 CDS sequences affects the codon-anticodon adaptation. The wobble site of the anticodon is always G for NNY codon families dominated by C-ending codons and U for NNR and NNN codon families dominated by A-ending codons. The only, but consistent, exception is the anticodon of tRNA-Met which consistently has a 5'-CAU-3' anticodon base-pairing with the AUG codon (the translation initiation codon) instead of the more frequent AUA. The observed CAU anticodon (matching AUG) would increase the rate of translation initiation but would reduce the rate of peptide elongation because most methionine codons are AUA, whereas the unobserved UAU anticodon (matching AUA) would increase the elongation rate at the cost of translation initiation rate. The consistent CAU anticodon in tRNA-Met suggests the importance of maximizing the rate of translation initiation.     

Understanding Biases in Ribosome Profiling Experiments Reveals Signatures of Translation Dynamics in Yeast

Ribosome profiling produces snapshots of the locations of actively translating ribosomes on messenger RNAs. These snapshots can be used to make inferences about translation dynamics. Recent ribosome profiling studies in yeast, however, have reached contradictory conclusions regarding the average translation rate of each codon. Some experiments have used cycloheximide (CHX) to stabilize ribosomes before measuring their positions, and these studies all counterintuitively report a weak negative correlation between the translation rate of a codon and the abundance of its cognate tRNA. In contrast, some experiments performed without CHX report strong positive correlations. To explain this contradiction, we identify unexpected patterns in ribosome density downstream of each type of codon in experiments that use CHX. These patterns are evidence that elongation continues to occur in the presence of CHX but with dramatically altered codon-specific elongation rates. The measured positions of ribosomes in these experiments therefore do not reflect the amounts of time ribosomes spend at each position in vivo. These results suggest that conclusions from experiments in yeast using CHX may need reexamination. In particular, we show that in all such experiments, codons decoded by less abundant tRNAs were in fact being translated more slowly before the addition of CHX disrupted these dynamics.     

Codon replacement in the PGK1 gene of Saccharomyces cerevisiae: experimental approach to study the role of biased codon usage in gene expression.

The coding sequences of genes in the yeast Saccharomyces cerevisiae show a preference for 25 of the 61 possible coding triplets. The degree of this biased codon usage in each gene is positively correlated to its expression level. Highly expressed genes use these 25 major codons almost exclusively. As an experimental approach to studying biased codon usage and its possible role in modulating gene expression, systematic codon replacements were carried out in the highly expressed PGK1 gene. The expression of phosphoglycerate kinase (PGK) was studied both on a high-copy-number plasmid and as a single copy gene integrated into the chromosome. Replacing an increasing number (up to 39% of all codons) of major codons with synonymous minor ones at the 5' end of the coding sequence caused a dramatic decline of the expression level. The PGK protein levels dropped 10-fold. The steady-state mRNA levels also declined, but to a lesser extent (threefold). Our data indicate that this reduction in mRNA levels was due to destabilization caused by impaired translation elongation at the minor codons. By preventing translation of the PGK mRNAs by the introduction of a stop codon 3' and adjacent to the start codon, the steady-state mRNA levels decreased dramatically. We conclude that efficient mRNA translation is required for maintaining mRNA stability in S. cerevisiae. These findings have important implications for the study of the expression of heterologous genes in yeast cells.     

同义密码子使用模式对多肽链共翻译折叠的影响研究进展

同义密码子使用模式作为核苷酸与氨基酸的纽带,其多样性介导了核糖体扫描速率,同时扩充了基因的遗传信息存储量。随着新型技术的应用,发现特异性密码子和密码子结合力可调节核糖体扫描速率并影响蛋白质构象。同义密码子使用模式通过多种方式在不同环节影响着核糖体扫描速率,同时还影响着自身mRNA的稳定性。本文简述了密码子使用模式如何在核糖体扫描翻译mRNA的过程中实现对多肽链翻译延伸的调控,为今后生物工程学领域如何优化蛋白高效表达提供可参考的思路与理念。     

The functional half-life of an mRNA depends on the ribosome spacing in an early coding region

Bacterial mRNAs are translated by closely spaced ribosomes and degraded from the 5′-end, with half-lives of around 2 min at 37 °C in most cases. Ribosome-free or “naked” mRNA is known to be readily degraded, but the initial event that inactivates the mRNA functionally has not been fully described. Here, we characterize a determinant of the functional stability of an mRNA, which is located in the early coding region. Using literature values for the mRNA half-lives of variant lacZ mRNAs in Escherichia coli, we modeled how the ribosome spacing is affected by the translation rate of the individual codons. When comparing the ribosome spacing at various segments of the mRNA to its functional half-life, we found a clear correlation between the functional mRNA half-life and the ribosome spacing in the mRNA region approximately between codon 20 and codon 45. From this finding, we predicted that inserts of slowly translated codons before codon 20 or after codon 45 should shorten or prolong, respectively, the functional mRNA half-life by altering the ribosome density in the important region. These predictions were tested on eight new lacZ variants, and their experimentally determined mRNA half-lives all supported the model. We thus suggest that translation-rate-mediated differences in the spacing between ribosomes in this early coding region is a parameter that determines the mRNAs functional half-life. We present a model that is in accordance with many earlier observations and that allows a prediction of the functional half-life of a given mRNA sequence.     

Silent substitutions predictably alter translation elongation rates and protein folding efficiencies

Genetic code redundancy allows most amino acids to be encoded by multiple codons that are non-randomly distributed along coding sequences. An accepted theory explaining the biological significance of such non-uniform codon selection is that codons are translated at different speeds. Thus, varying codon placement along a message may confer variable rates of polypeptide emergence from the ribosome, which may influence the capacity to fold toward the native state. Previous studies report conflicting results regarding whether certain codons correlate with particular structural or folding properties of the encoded protein. This is partly due to different criteria traditionally utilized for predicting translation speeds of codons, including their usage frequencies and the concentration of tRNA species capable of decoding them, which do not always correlate. Here, we developed a metric to predict organism-specific relative translation rates of codons based on the availability of tRNA decoding mechanisms: Watson-Crick, non-Watson-Crick or both types of interactions. We determine translation rates of messages by pulse-chase analyses in living Escherichia coli cells and show that sequence engineering based on these concepts predictably modulates translation rates in a manner that is superior to codon usage frequency, which occur during the elongation phase, and significantly impacts folding of the encoded polypeptide. Finally, we demonstrate that sequence harmonization based on expression host tRNA pools, designed to mimic ribosome movement of the original organism, can significantly increase the folding of the encoded polypeptide. These results illuminate how genetic code degeneracy may function to specify properties beyond amino acid encoding, including folding.     

A model for protein translation: polysome self-organization leads to maximum protein synthesis rates

The genetic information in DNA is transcribed to mRNA and then translated to proteins, which form the building blocks of life. Translation, or protein synthesis, is hence a central cellular process. We have developed a gene-sequence-specific mechanistic model for the translation machinery, which accounts for all the elementary steps of the translation mechanism. We performed a sensitivity analysis to determine the effects of kinetic parameters and concentrations of the translational components on protein synthesis rate. Utilizing our mathematical framework and sensitivity analysis, we investigated the translational kinetic properties of a single mRNA species in Escherichia coli. We propose that translation rate at a given polysome size depends on the complex interplay between ribosomal occupancy of elongation phase intermediate states and ribosome distributions with respect to codon position along the length of the mRNA, and this interplay leads to polysome self-organization that drives translation rate to maximum levels.     

Effects of codon distributions and tRNA competition on protein translation

Translation is a central cellular process and the complexity of its mechanism necessitates mathematical frameworks to better understand system properties and make quantitative predictions. We have developed a gene sequence-specific mechanistic model for translation which accounts for all the elementary steps of translation elongation. Included in our model is the nonspecific binding of tRNAs to the ribosomal A site, and we find that the competitive, nonspecific binding of the tRNAs is the rate-limiting step in the elongation cycle for every codon. By introducing our model in terms of the Michaelis-Menten kinetic framework, we determine that these results are due to the tRNAs that do not recognize the ribosomal A site codon acting as competitive inhibitors to the tRNAs that do recognize the ribosomal A site codon. We present the results of a sensitivity analysis to determine the contribution of elongation cycle kinetic parameters of each codon on the overall translation rate, and observe that the translation rates of mRNAs are controlled by segments of rate-limiting codons that are sequence-specific. Along these lines, we find that the relative position of codons along the mRNA determines the optimal protein synthesis rate.