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.     

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.     

Deciphering the rules by which dynamics of mRNA secondary structure affect translation efficiency in Saccharomyces cerevisiae

Messenger RNA (mRNA) secondary structure decreases the elongation rate, as ribosomes must unwind every structure they encounter during translation. Therefore, the strength of mRNA secondary structure is assumed to be reduced in highly translated mRNAs. However, previous studies in vitro reported a positive correlation between mRNA folding strength and protein abundance. The counterintuitive finding suggests that mRNA secondary structure affects translation efficiency in an undetermined manner. Here, we analyzed the folding behavior of mRNA during translation and its effect on translation efficiency. We simulated translation process based on a novel computational model, taking into account the interactions among ribosomes, codon usage and mRNA secondary structures. We showed that mRNA secondary structure shortens ribosomal distance through the dynamics of folding strength. Notably, when adjacent ribosomes are close, mRNA secondary structures between them disappear, and codon usage determines the elongation rate. More importantly, our results showed that the combined effect of mRNA secondary structure and codon usage in highly translated mRNAs causes a short ribosomal distance in structural regions, which in turn eliminates the structures during translation, leading to a high elongation rate. Together, these findings reveal how the dynamics of mRNA secondary structure coupling with codon usage affect translation efficiency.     

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.     

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.     

A Major Controversy in Codon-Anticodon Adaptation Resolved by a New Codon Usage Index

Two alternative hypotheses attribute different benefits to codon-anticodon adaptation. The first assumes that protein production is rate limited by both initiation and elongation and that codon-anticodon adaptation would result in higher elongation efficiency and more efficient and accurate protein production, especially for highly expressed genes. The second claims that protein production is rate limited only by initiation efficiency but that improved codon adaptation and, consequently, increased elongation efficiency have the benefit of increasing ribosomal availability for global translation. To test these hypotheses, a recent study engineered a synthetic library of 154 genes, all encoding the same protein but differing in degrees of codon adaptation, to quantify the effect of differential codon adaptation on protein production in Escherichia coli. The surprising conclusion that “codon bias did not correlate with gene expression” and that “translation initiation, not elongation, is rate-limiting for gene expression” contradicts the conclusion reached by many other empirical studies. In this paper, I resolve the contradiction by reanalyzing the data from the 154 sequences. I demonstrate that translation elongation accounts for about 17% of total variation in protein production and that the previous conclusion is due to the use of a codon adaptation index (CAI) that does not account for the mutation bias in characterizing codon adaptation. The effect of translation elongation becomes undetectable only when translation initiation is unrealistically slow. A new index of translation elongation I_TE is formulated to facilitate studies on the efficiency and evolution of the translation machinery.     

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.     

Selective constraints in yeast genes with differential expressivity: codon pair usage and mRNA stability perspectives

Protein translation has been elucidated to be dictated by evolutionary constraints, namely, variations in tRNA availabilities and/or variations in codon-anticodon binding that is manifested in biased codon usage. Taking advantage of publicly available mRNA expression and protein abundance data for Saccharomyces cerevisiae, we have performed a comprehensive analysis of the diverse factors guiding translation leading to desired protein levels irrespective of the corresponding high or low mRNA levels. It has been elucidated in this study that different combinations of most abundant/non abundant tRNA isoacceptors are selected for in S. cerevisiae that helps in achieving the optimum speed and accuracy in the protein translation process. This is also accompanied by the strategic location of codon pairs in coherence to mRNA secondary structure folding stability for the above mentioned combinations of tRNA isoacceptors. We thus find that codon pair contextual effects; in addition to tRNA abundance and mRNA folding stability during translation elongation process play plausible roles in maintaining translation accuracy and speed that can achieve desired protein levels.     

Genetic depletion of the RNA helicase DDX3 leads to impaired elongation of translating ribosomes triggering co-translational quality control of newly synthesized polypeptides

DDX3 is a multifaceted RNA helicase of the DEAD-box family that plays central roles in all aspects of RNA metabolism including translation initiation. Here, we provide evidence that the Leishmania DDX3 ortholog functions in post-initiation steps of translation. We show that genetic depletion of DDX3 slows down ribosome movement resulting in elongation-stalled ribosomes, impaired translation elongation and decreased de novo protein synthesis. We also demonstrate that the essential ribosome recycling factor Rli1/ABCE1 and termination factors eRF3 and GTPBP1 are less recruited to ribosomes upon DDX3 loss, suggesting that arrested ribosomes may be inefficiently dissociated and recycled. Furthermore, we show that prolonged ribosome stalling triggers co-translational ubiquitination of nascent polypeptide chains and a higher recruitment of E3 ubiquitin ligases and proteasome components to ribosomes of DDX3 knockout cells, which further supports that ribosomes are not elongating optimally. Impaired elongation of translating ribosomes also results in the accumulation of cytoplasmic protein aggregates, which implies that defects in translation overwhelm the normal quality controls. The partial recovery of translation by overexpressing Hsp70 supports this possibility. Collectively, these results suggest an important novel contribution of DDX3 to optimal elongation of translating ribosomes by preventing prolonged translation stalls and stimulating recycling of arrested ribosomes.     

Causal signals between codon bias,mRNA structure,and the efficiency of translation and elongation

Ribosome profiling data report on the distribution of translating ribosomes, at steady‐state, with codon‐level resolution. We present a robust method to extract codon translation rates and protein synthesis rates from these data, and identify causal features associated with elongation and translation efficiency in physiological conditions in yeast. We show that neither elongation rate nor translational efficiency is improved by experimental manipulation of the abundance or body sequence of the rare AGG tRNA. Deletion of three of the four copies of the heavily used ACA tRNA shows a modest efficiency decrease that could be explained by other rate‐reducing signals at gene start. This suggests that correlation between codon bias and efficiency arises as selection for codons to utilize translation machinery efficiently in highly translated genes. We also show a correlation between efficiency and RNA structure calculated both computationally and from recent structure probing data, as well as the Kozak initiation motif, which may comprise a mechanism to regulate initiation.     

YaeJ is a novel ribosome-associated protein in Escherichia coli that can hydrolyze peptidyl-tRNA on stalled ribosomes

In bacteria, ribosomes often become stalled and are released by a trans-translation process mediated by transfer-messenger RNA (tmRNA). In the absence of tmRNA, however, there is evidence that stalled ribosomes are released from non-stop mRNAs. Here, we show a novel ribosome rescue system mediated by a small basic protein, YaeJ, from Escherichia coli, which is similar in sequence and structure to the catalytic domain 3 of polypeptide chain release factor (RF). In vitro translation experiments using the E. coli-based reconstituted cell-free protein synthesis system revealed that YaeJ can hydrolyze peptidyl–tRNA on ribosomes stalled by both non-stop mRNAs and mRNAs containing rare codon clusters that extend downstream from the P-site and prevent Ala-tmRNA•SmpB from entering the empty A-site. In addition, YaeJ had no effect on translation of a normal mRNA with a stop codon. These results suggested a novel tmRNA-independent rescue system for stalled ribosomes in E. coli. YaeJ was almost exclusively found in the 70S ribosome and polysome fractions after sucrose density gradient sedimentation, but was virtually undetectable in soluble fractions. The C-terminal basic residue-rich extension was also found to be required for ribosome binding. These findings suggest that YaeJ functions as a ribosome-attached rescue device for stalled ribosomes.     

Generation of pokeweed antiviral protein mutations in Saccharomyces cerevisiae: evidence that ribosome depurination is not sufficient for cytotoxicity

Pokeweed antiviral protein (PAP) is a ribosome-inactivating protein that depurinates the highly conserved α-sarcin/ricin loop in the large rRNA. Here, using site-directed mutagenesis and systematic deletion analysis from the 5′ and the 3′ ends of the PAP cDNA, we identified the amino acids important for ribosome depurination and cytotoxicity of PAP. Truncating the first 16 amino acids of PAP eliminated its cytotoxicity and the ability to depurinate ribosomes. Ribosome depurination gradually decreased upon the sequential deletion of C-terminal amino acids and was abolished when a stop codon was introduced at Glu-244. Cytotoxicity of the C-terminal deletion mutants was lost before their ability to depurinate ribosomes. Mutations in Tyr-123 at the active site affected cytotoxicity without altering the ribosome depurination ability. Total translation was not inhibited in yeast expressing the non-toxic Tyr-123 mutants, although ribosomes were depurinated. These mutants depurinated ribosomes only during their translation and could not depurinate ribosomes in trans in a translation-independent manner. A mutation in Leu-71 in the central domain affected cytotoxicity without altering the ability to depurinate ribosomes in trans and inhibit translation. These results demonstrate that the ability to depurinate ribosomes in trans in a catalytic manner is required for the inhibition of translation, but is not sufficient for cytotoxicity.     

Ribosome profiling reveals sequence-independent post-initiation pausing as a signature of translation

The journey of a newly synthesized polypeptide starts in the peptidyltransferase center of the ribosome, from where it traverses the exit tunnel. The interior of the ribosome exit tunnel is neither straight nor smooth. How the ribosome dynamics in vivo is influenced by the exit tunnel is poorly understood. Genome-wide ribosome profiling in mammalian cells reveals elevated ribosome density at the start codon and surprisingly the downstream 5th codon position as well. We found that the highly focused ribosomal pausing shortly after initiation is attributed to the geometry of the exit tunnel, as deletion of the loop region from ribosome protein L4 diminishes translational pausing at the 5th codon position. Unexpectedly, the ribosome variant undergoes translational abandonment shortly after initiation, suggesting that there exists an obligatory step between initiation and elongation commitment. We propose that the post-initiation pausing of ribosomes represents an inherent signature of the translation machinery to ensure productive translation.     

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.     

The h subunit of eIF3 promotes reinitiation competence during translation of mRNAs harboring upstream open reading frames

Upstream open reading frames (uORFs) are protein coding elements in the 5′ leader of messenger RNAs. uORFs generally inhibit translation of the main ORF because ribosomes that perform translation elongation suffer either permanent or conditional loss of reinitiation competence. After conditional loss, reinitiation competence may be regained by, at the minimum, reacquisition of a fresh methionyl-tRNA. The conserved h subunit of Arabidopsis eukaryotic initiation factor 3 (eIF3) mitigates the inhibitory effects of certain uORFs. Here, we define more precisely how this occurs, by combining gene expression data from mutated 5′ leaders of Arabidopsis AtbZip11 (At4g34590) and yeast GCN4 with a computational model of translation initiation in wild-type and eif3h mutant plants. Of the four phylogenetically conserved uORFs in AtbZip11, three are inhibitory to translation, while one is anti-inhibitory. The mutation in eIF3h has no major effect on uORF start codon recognition. Instead, eIF3h supports efficient reinitiation after uORF translation. Modeling suggested that the permanent loss of reinitiation competence during uORF translation occurs at a faster rate in the mutant than in the wild type. Thus, eIF3h ensures that a fraction of uORF-translating ribosomes retain their competence to resume scanning. Experiments using the yeast GCN4 leader provided no evidence that eIF3h fosters tRNA reaquisition. Together, these results attribute a specific molecular function in translation initiation to an individual eIF3 subunit in a multicellular eukaryote.     

Transient idling of posttermination ribosomes ready to reinitiate protein synthesis

The fate of ribosomes between termination and initiation during protein synthesis is very basic, yet poorly understood. Here we found that translational reinitiation of the alkaline phosphatase gene occurs in Escherichia coli from an internal methionine codon when the authentic translation is prematurely terminated at a nonsense codon that is within seven codons upstream of the reinitiation codon (which we refer to as "reinitiation window"). Changing the reading frame downstream of the stop codon did not abolish the reinitiation, while inactivating the upstream initiation codon abolished the reinitiation. Moreover, depletion of the ribosome recycling factor (RRF), which disassembles posttermination ribosomes in conjunction with elongation factor G, did not influence the observed reinitiation. These findings suggest that posttermination ribosomes can undergo a transient idling state ready to reinitiate protein synthesis even in the absence of the Shine-Dalgarno (SD) sequence within the reinitiation window by evading disengagement from the mRNA.     

Base 2661 in Escherichia coli 23S rRNA influences the binding of elongation factor Tu during protein synthesis in vivo.

The binding of the EF-Tu.GTP.aminoacyl-tRNA ternary complex (EF, elongation factor) to the ribosome is known to be strengthened by a 2661G-to-C mutation in 23S ribosomal RNA, whereas the binding to normal ribosomes is weakened if the factor is in an appropriate mutant form (Aa). In this report we describe the mutual effects by the 2661C alteration in 23S rRNA and EF-Tu(Aa) on bacterial viability and translation efficiency in strains with normal or mutationally altered ribosomes. The rrnB(2661C) allele on a multicopy plasmid was introduced by transformation into Escherichia coli K-12 strains, harbouring either the wild-type or the mutant gene (tufA) for EF-Tu as well as normal or mutant ribosomal protein S12 (rpsL). Together with wild-type EF-Tu, the 2661C mutant ribosomes decreased the translation elongation rate in a rpsL+ strain or a non-restrictive rpsL224 strain. This reduction was not seen in strains which harbored EF-Tu(Aa) instead of EF-Tu(As) (As, wild-type form). Nonsense codon suppression by tyrT(Su3) suppressor tRNA was reduced by 2661C in a rpsL224 strain in the presence of EF-Tu(As) but not in the presence of EF-Tu(Aa). The lethal effect obtained by the combination of 2661C and a restrictive ribosomal protein S12 mutation (rpsL282) disappeared if EF-Tu(As) was replaced by EF-Tu(Aa) in the strain. In such a viable strain, 2661C had no effect on either the translation elongation rate or nonsense codon suppression. Our data suggest that the G base at position 2661 in 23S rRNA is important for binding of EF-Tu during protein synthesis in vivo. The interaction between this base and EF-Tu is strongly influenced by the structure of ribosomal protein S12.