Dynamic changes in translational efficiency are deduced from codon usage of the transcriptome

tRNA expression and modification landscapes, and their dynamics during zebrafish embryo development

tRNA genes exist in multiple copies in the genome of all organisms across the three domains of life. Besides the sequence differences across tRNA copies, extensive post-transcriptional modification adds a further layer to tRNA diversification. Whilst the crucial role of tRNAs as adapter molecules in protein translation is well established, whether all tRNAs are actually expressed, and whether the differences across isodecoders play any regulatory role is only recently being uncovered. Here we built upon recent developments in the use of NGS-based methods for RNA modification detection and developed tRAM-seq, an experimental protocol and in silico analysis pipeline to investigate tRNA expression and modification. Using tRAM-seq, we analysed the full ensemble of nucleo-cytoplasmic and mitochondrial tRNAs during embryonic development of the model vertebrate zebrafish. We show that the repertoire of tRNAs changes during development, with an apparent major switch in tRNA isodecoder expression and modification profile taking place around the start of gastrulation. Taken together, our findings suggest the existence of a general reprogramming of the expressed tRNA pool, possibly gearing the translational machinery for distinct stages of the delicate and crucial process of embryo development.

MatK impacts differential chloroplast translation by limiting spliced tRNA-K(UUU) abundance

The protein levels of chloroplast photosynthetic genes and genes related to the chloroplast genetic apparatus vary to adapt to different conditions. However, the underlying mechanisms governing these variations remain unclear. The chloroplast intron Maturase K is encoded within the trnK intron and has been suggested to be required for splicing several group IIA introns, including the trnK intron. In this study, we used RNA immunoprecipitation followed by high-throughput sequencing (RIP-Seq) to identify MatK's preference for binding to group IIA intron domains I and VI within target transcripts. Importantly, these domains are crucial for splice site selection, and we discovered alternative 5'-splice sites in three MatK target introns. The resulting alternative trnK lariat structure showed increased accumulation during heat acclimation. The cognate codon of tRNA-K(UUU) is highly enriched in mRNAs encoding ribosomal proteins and a trnK-matK over-expressor exhibited elevated levels of the spliced tRNA-K(UUU). Ribosome profiling analysis of the overexpressor revealed a significant up-shift in the translation of ribosomal proteins compared to photosynthetic genes. Our findings suggest the existence of a novel regulatory mechanism linked to the abundance of tRNA-K(UUU), enabling the differential expression of functional chloroplast gene groups.

Codon usage bias in yeasts and its correlation with gene expression, growth temperature, and protein structure

Codon usage bias (CUB) has been described in viruses, prokaryotes, and eukaryotes and has been linked to several cellular and environmental factors, such as the organism's growth temperature, gene expression levels, and regulation of protein synthesis and folding. Most of the studies in this area have been conducted in bacteria and higher eukaryotes, in some cases with different results. In this study, a comparative analysis of CUB in yeasts isolated from cold and template environments was performed in order to evaluate the correlation of CUB with yeast optimal temperature of growth (OTG), gene expression levels, cellular function, and structure of encoded proteins. Among the main findings, highly expressed ORFs tend to have a more similar CUB within and between yeasts, and a direct correlation between codons ending in C and expression level was generally found. A low correspondence between CUB and OTG was observed, with an inverse correlation for some codons ending in C. The clustering of yeasts based on their CUB partially aligns with their OTG, being more consistent for yeasts with lower OTG. In most yeasts, the abundance of preferred codons was generally lower at the 5' end of ORFs, higher in segments encoding beta strand, lower in segments encoding extracellular and transmembrane regions, and higher in "translation" and "energy metabolism" pathways, especially in highly expressed ORFs. Based on our findings, it is suggested that the abundance and distribution of preferred and non-preferred codons along mRNAs contribute to proper protein folding and functionality by regulating protein synthesis rates, becoming a more important factor under conditions that require faster protein synthesis in yeasts.

tRNA epitranscriptome determines pathogenicity of the opportunistic pathogen Pseudomonas aeruginosa

The success of bacterial pathogens depends on the coordinated expression of virulence determinants. Regulatory circuits that drive pathogenesis are complex, multilayered, and incompletely understood. Here, we reveal that alterations in tRNA modifications define pathogenic phenotypes in the opportunistic pathogen Pseudomonas aeruginosa. We demonstrate that the enzymatic activity of GidA leads to the introduction of a carboxymethylaminomethyl modification in selected tRNAs. Modifications at the wobble uridine base (cmnm5U34) of the anticodon drives translation of transcripts containing rare codons. Specifically, in P. aeruginosa the presence of GidA-dependent tRNA modifications modulates expression of genes encoding virulence regulators, leading to a cellular proteomic shift toward pathogenic and well-adapted physiological states. Our approach of profiling the consequences of chemical tRNA modifications is general in concept. It provides a paradigm of how environmentally driven tRNA modifications govern gene expression programs and regulate phenotypic outcomes responsible for bacterial adaption to challenging habitats prevailing in the host niche.

Identification of hyperosmotic stress-responsive genes in Chinese hamster ovary cells via genome-wide virus-free CRISPR/Cas9 screening

Chinese hamster ovary (CHO) cells are the preferred mammalian host cells for therapeutic protein production that have been extensively engineered to possess the desired attributes for high-yield protein production. However, empirical approaches for identifying novel engineering targets are laborious and time-consuming. Here, we established a genome-wide CRISPR/Cas9 screening platform for CHO-K1 cells with 111,651 guide RNAs (gRNAs) targeting 21,585 genes using a virus-free recombinase-mediated cassette exchange-based gRNA integration method. Using this platform, we performed a positive selection screening under hyperosmotic stress conditions and identified 180 genes whose perturbations conferred resistance to hyperosmotic stress in CHO cells. Functional enrichment analysis identified hyperosmotic stress responsive gene clusters, such as tRNA wobble uridine modification and signaling pathways associated with cell cycle arrest. Furthermore, we validated 32 top-scoring candidates and observed a high rate of hit confirmation, demonstrating the potential of the screening platform. Knockout of the novel target genes, Zfr and Pnp, in monoclonal antibody (mAb)-producing recombinant CHO (rCHO) cells and bispecific antibody (bsAb)-producing rCHO cells enhanced their resistance to hyperosmotic stress, thereby improving mAb and bsAb production. Overall, the collective findings demonstrate the value of the screening platform as a powerful tool to investigate the functions of genes associated with hyperosmotic stress and to discover novel targets for rational cell engineering on a genome-wide scale in CHO cells.

Race with virus evolution: The development and application of mRNA vaccines against SARS-CoV-2

Since the COVID-19 pandemic was declared, vaccines against SARS-CoV-2 have been urgently developed around the world. On the basis of the mRNA vaccine technology developed previously, COVID-19 mRNA vaccines were promptly tested in animals, advanced to clinical trials, and then authorized for emergency use in humans. The administration of COVID-19 mRNA vaccines has successfully reduced the hospitalization and mortality caused by the viral infection, although the virus continuously evolves with its transmission. Therefore, the development of mRNA vaccine technology, including RNA modification and delivery systems, is well recognized for its contribution to moderating the harms caused by the COVID-19 pandemic. The scientists who developed these technologies, Katalin Karikó, Drew Weissman, and Pieter Cullis, were awarded the 2022 Tang Prize in Biopharmaceutical Science. In this review, we summarize the principles, safety and efficacy of as well as the immune response to COVID-19 mRNA vaccines. Since mRNA vaccine approaches could be practical for the prevention of infectious diseases, we also briefly describe mRNA vaccines against other human viral pathogens in clinical trials.

The Ramp Atlas: facilitating tissue and cell-specific ramp sequence analyses through an intuitive web interface

Ramp sequences occur when the average translational efficiency of codons near the 5′ end of highly expressed genes is significantly lower than the rest of the gene sequence, which counterintuitively increases translational efficiency by decreasing downstream ribosomal collisions. Here, we show that the relative codon adaptiveness within different tissues changes the existence of a ramp sequence without altering the underlying genetic code. We present the first comprehensive analysis of tissue and cell type-specific ramp sequences and report 3108 genes with ramp sequences that change between tissues and cell types, which corresponds with increased gene expression within those tissues and cells. The Ramp Atlas (https://ramps.byu.edu/) allows researchers to query precomputed ramp sequences in 18 388 genes across 62 tissues and 66 cell types and calculate tissue-specific ramp sequences from user-uploaded FASTA files through an intuitive web interface. We used The Ramp Atlas to identify seven SARS-CoV-2 genes and seven human SARS-CoV-2 entry factor genes with tissue-specific ramp sequences that may help explain viral proliferation within those tissues. We anticipate that The Ramp Atlas will facilitate personalized and creative tissue-specific ramp sequence analyses for both human and viral genes that will increase our ability to utilize this often-overlooked regulatory region.

Upregulation of tRNA-Ser-AGA-2-1 Promotes Malignant Behavior in Normal Bronchial Cells

Serine tRNAs (tRNA^Ser) are frequently overexpressed in tumors and associated with poor prognosis and increased risk of recurrence in breast cancer. Impairment of tRNA biogenesis and abundance also impacts proteome homeostasis, and activates protein quality control systems. Herein, we aimed at testing whether increasing tRNA^Ser abundance could foster tumor establishment through activation of the UPR. In order to do so, firstly we confirmed that the expression of tRNA-Ser-AGA-2-1 [hereafter tRNA^Ser(AGA)] was upregulated by 1.79-fold in Stage I NSCLC tumors when compared to normal adjacent tissue. To study the impact of tRNA^Ser(AGA) in early stage tumorigenesis, we induced its upregulation in a non-tumoral bronchial cell line, BEAS-2B. Upregulation of this tRNA increased cellular proliferation and protein synthesis rate, driven by eIF2α dephosphorylation and ATF4 activation downstream of PERK signaling. Futhermore, tRNA^Ser(AGA) enhanced transformation potential in vitro, and promoted the establishment of slow growing tumors with aggressive features in nude mice. Our work highlights the importance of studying tRNA deregulation on early stage tumorigenesis, as they may be potential malignancy and aggressiveness biomarkers.

Cold Regulation of Genes Encoding Ion Transport Systems in the Oligotrophic Bacterium Caulobacter crescentus

In this study, we characterize the response of the free-living oligotrophic alphaproteobacterium Caulobacter crescentus to low temperatures by global transcriptomic analysis. Our results showed that 656 genes were upregulated and 619 were downregulated at least 2-fold after a temperature downshift. The identified differentially expressed genes (DEG) belong to several functional categories, notably inorganic ion transport and metabolism, and a subset of these genes had their expression confirmed by reverse transcription quantitative real-time PCR (RT-qPCR). Several genes belonging to the ferric uptake regulator (Fur) regulon were downregulated, indicating that iron homeostasis is relevant for adaptation to cold. Several upregulated genes encode proteins that interact with nucleic acids, particularly RNA: cspA, cspB, and the DEAD box RNA helicases rhlE, dbpA, and rhlB. Moreover, 31 small regulatory RNAs (sRNAs), including the cell cycle-regulated noncoding RNA (ncRNA) CcnA, were upregulated, indicating that posttranscriptional regulation is important for the cold stress response. Interestingly, several genes related to transport were upregulated under cold stress, including three AcrB-like cation/multidrug efflux pumps, the nitrate/nitrite transport system, and the potassium transport genes kdpFABC. Further characterization showed that kdpA is upregulated in a potassium-limited medium and at a low temperature in a SigT-independent way. kdpA mRNA is less stable in rho and rhlE mutant strains, but while the expression is positively regulated by RhlE, it is negatively regulated by Rho. A kdpA-deleted strain was generated, and its viability in response to osmotic, acidic, or cold stresses was determined. The implications of such variation in the gene expression for cold adaptation are discussed. IMPORTANCE Low-temperature stress is an important factor for nucleic acid stability and must be circumvented in order to maintain the basic cell processes, such as transcription and translation. The oligotrophic lifestyle presents further challenges to ensure the proper nutrient uptake and osmotic balance in an environment of slow nutrient flow. Here, we show that in Caulobacter crescentus, the expression of the genes involved in cation transport and homeostasis is altered in response to cold, which could lead to a decrease in iron uptake and an increase in nitrogen and high-affinity potassium transport by the Kdp system. This previously uncharacterized regulation of the Kdp transporter has revealed a new mechanism for adaptation to low temperatures that may be relevant for oligotrophic bacteria.

The genomic loci of specific human tRNA genes exhibit ageing-related DNA hypermethylation

The epigenome has been shown to deteriorate with age, potentially impacting on ageing-related disease. tRNA, while arising from only ˜46 kb (<0.002% genome), is the second most abundant cellular transcript. tRNAs also control metabolic processes known to affect ageing, through core translational and additional regulatory roles. Here, we interrogate the DNA methylation state of the genomic loci of human tRNA. We identify a genomic enrichment for age-related DNA hypermethylation at tRNA loci. Analysis in 4,350 MeDIP-seq peripheral-blood DNA methylomes (16-82 years), identifies 44 and 21 hypermethylating specific tRNAs at study-and genome-wide significance, respectively, contrasting with none hypomethylating. Validation and replication (450k array and independent targeted Bisuphite-sequencing) supported the hypermethylation of this functional unit. Tissue-specificity is a significant driver, although the strongest consistent signals, also independent of major cell-type change, occur in tRNA-iMet-CAT-1-4 and tRNA-Ser-AGA-2-6. This study presents a comprehensive evaluation of the genomic DNA methylation state of human tRNA genes and reveals a discreet hypermethylation with advancing age.

Adaptation of codon and amino acid use for translational functions in highly expressed cricket genes

BACKGROUND

For multicellular organisms, much remains unknown about the dynamics of synonymous codon and amino acid use in highly expressed genes, including whether their use varies with expression in different tissue types and sexes. Moreover, specific codons and amino acids may have translational functions in highly transcribed genes, that largely depend on their relationships to tRNA gene copies in the genome. However, these relationships and putative functions are poorly understood, particularly in multicellular systems.

RESULTS

Here, we studied codon and amino acid use in highly expressed genes from reproductive and nervous system tissues (male and female gonad, somatic reproductive system, brain and ventral nerve cord, and male accessory glands) in the cricket Gryllus bimaculatus. We report an optimal codon, defined as the codon preferentially used in highly expressed genes, for each of the 18 amino acids with synonymous codons in this organism. The optimal codons were mostly shared among tissue types and both sexes. However, the frequency of optimal codons was highest in gonadal genes. Concordant with translational selection, a majority of the optimal codons had abundant matching tRNA gene copies in the genome, but sometimes obligately required wobble tRNAs. We suggest the latter may comprise a mechanism for slowing translation of abundant transcripts, particularly for cell-cycle genes. Non-optimal codons, defined as those least commonly used in highly transcribed genes, intriguingly often had abundant tRNAs, and had elevated use in a subset of genes with specialized functions (gametic and apoptosis genes), suggesting their use promotes the translational upregulation of particular mRNAs. In terms of amino acids, we found evidence suggesting that amino acid frequency, tRNA gene copy number, and amino acid biosynthetic costs (size/complexity) had all interdependently evolved in this insect model, potentially for translational optimization.

CONCLUSIONS

Collectively, the results suggest a model whereby codon use in highly expressed genes, including optimal, wobble, and non-optimal codons, and their tRNA abundances, as well as amino acid use, have been influenced by adaptation for various functional roles in translation within this cricket. The effects of expression in different tissue types and the two sexes are discussed.

Whole-Genome Transformation Promotes tRNA Anticodon Suppressor Mutations under Stress

tRNAs are encoded by a large gene family, usually with several isogenic tRNAs interacting with the same codon. Mutations in the anticodon region of other tRNAs can overcome specific tRNA deficiencies. Phylogenetic analysis suggests that such mutations have occurred in evolution, but the driving force is unclear. We show that in yeast suppressor mutations in other tRNAs are able to overcome deficiency of the essential TRT2-encoded tRNAThrCGU at high temperature (40°C). Surprisingly, these tRNA suppressor mutations were obtained after whole-genome transformation with DNA from thermotolerant Kluyveromyces marxianus or Ogataea polymorpha strains but from which the mutations did apparently not originate. We suggest that transient presence of donor DNA in the host facilitates proliferation at high temperature and thus increases the chances for occurrence of spontaneous mutations suppressing defective growth at high temperature. Whole-genome sequence analysis of three transformants revealed only four to five nonsynonymous mutations of which one causing TRT2 anticodon stem stabilization and two anticodon mutations in non-threonyl-tRNAs, tRNALysCUU and tRNAeMetCAU, were causative. Both anticodon mutations suppressed lethality of TRT2 deletion and apparently caused the respective tRNAs to become novel substrates for threonyl-tRNA synthetase. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) data could not detect any significant mistranslation, and reverse transcription-quantitative PCR results contradicted induction of the unfolded protein response. We suggest that stress conditions have been a driving force in evolution for the selection of anticodon-switching mutations in tRNAs as revealed by phylogenetic analysis.IMPORTANCE In this work, we have identified for the first time the causative elements in a eukaryotic organism introduced by applying whole-genome transformation and responsible for the selectable trait of interest, i.e., high temperature tolerance. Surprisingly, the whole-genome transformants contained just a few single nucleotide polymorphisms (SNPs), which were unrelated to the sequence of the donor DNA. In each of three independent transformants, we have identified a SNP in a tRNA, either stabilizing the essential tRNAThrCGU at high temperature or switching the anticodon of tRNALysCUU or tRNAeMetCAU into CGU, which is apparently enough for in vivo recognition by threonyl-tRNA synthetase. LC-MS/MS analysis indeed indicated absence of significant mistranslation. Phylogenetic analysis showed that similar mutations have occurred throughout evolution and we suggest that stress conditions may have been a driving force for their selection. The low number of SNPs introduced by whole-genome transformation may favor its application for improvement of industrial yeast strains.

Translational adaptation of human viruses to the tissues they infect

Viruses need to hijack the translational machinery of the host cell for a productive infection to happen. However, given the dynamic landscape of tRNA pools among tissues, it is unclear whether different viruses infecting different tissues have adapted their codon usage toward their tropism. Here, we collect the coding sequences of 502 human-infecting viruses and determine that tropism explains changes in codon usage. Using the tRNA abundances across 23 human tissues from The Cancer Genome Atlas (TCGA), we build an in silico model of translational efficiency that validates the correspondence of the viral codon usage with the translational machinery of their tropism. For instance, we detect that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is specifically adapted to the upper respiratory tract and alveoli. Furthermore, this correspondence is specifically defined in early viral proteins. The observed tissue-specific translational efficiency could be useful for the development of antiviral therapies and vaccines.

The effects of codon bias and optimality on mRNA and protein regulation

The central dogma of molecular biology entails that genetic information is transferred from nucleic acid to proteins. Notwithstanding retro-transcribing genetic elements, DNA is transcribed to RNA which in turn is translated into proteins. Recent advancements have shown that each stage is regulated to control protein abundances for a variety of essential physiological processes. In this regard, mRNA regulation is essential in fine-tuning or calibrating protein abundances. In this review, we would like to discuss one of several mRNA-intrinsic features of mRNA regulation that has been gaining traction of recent-codon bias and optimality. Specifically, we address the effects of codon bias with regard to codon optimality in several biological processes centred on translation, such as mRNA stability and protein folding among others. Finally, we examine how different organisms or cell types, through this system, are able to coordinate physiological pathways to respond to a variety of stress or growth conditions.

Manipulation of the human tRNA pool reveals distinct tRNA sets that act in cellular proliferation or cell cycle arrest

Different subsets of the tRNA pool in human cells are expressed in different cellular conditions. The 'proliferation-tRNAs' are induced upon normal and cancerous cell division, while the 'differentiation-tRNAs' are active in non-dividing, differentiated cells. Here we examine the essentiality of the various tRNAs upon cellular growth and arrest. We established a CRISPR-based editing procedure with sgRNAs that each target a tRNA family. We measured tRNA essentiality for cellular growth and found that most proliferation-tRNAs are essential compared to differentiation- tRNAs in rapidly growing cell lines. Yet in more slowly dividing lines, the differentiation-tRNAs were more essential. In addition, we measured the essentiality of each tRNA family upon response to cell cycle arresting signals. Here we detected a more complex behavior with both proliferation-tRNAs and differentiation tRNAs showing various levels of essentiality. These results provide the so-far most comprehensive functional characterization of human tRNAs with intricate roles in various cellular states.

Translational quality control and reprogramming during stress adaptation

Organisms encounter stress throughout their lives, and therefore require the ability to respond rapidly to environmental changes. Although transcriptional responses are crucial for controlling changes in gene expression, regulation at the translational level often allows for a faster response at the protein levels which permits immediate adaptation. The fidelity and robustness of protein synthesis are actively regulated under stress. For example, mistranslation can be beneficial to cells upon environmental changes and also alters cellular stress responses. Additionally, stress modulates both global and selective translational regulation through mechanisms including the change of aminoacyl-tRNA activity, tRNA pool reprogramming and ribosome heterogeneity. In this review, we draw on studies from both the prokaryotic and eukaryotic systems to discuss current findings of cellular adaptation at the level of translation, specifically translational fidelity and activity changes in response to a wide array of environmental stressors including oxidative stress, nutrient depletion, temperature variation, antibiotics and host colonization.

Does proteostasis get lost in translation? Implications for protein aggregation across the lifespan

Protein aggregation is a phenomenon of major relevance in neurodegenerative and neuromuscular disorders, cataracts, diabetes and many other diseases. Research has unveiled that proteins also aggregate in multiple tissues during healthy aging yet, the biological and biomedical relevance of this apparently asymptomatic phenomenon remains to be understood. It is known that proteome homeostasis (proteostasis) is maintained by a balanced protein synthesis rate, high protein synthesis accuracy, efficient protein folding and continual tagging of damaged proteins for degradation, suggesting that protein aggregation during healthy aging may be associated with alterations in both protein synthesis and the proteostasis network (PN) pathways. In particular, dysregulation of protein synthesis and alterations in translation fidelity are hypothesized to lead to the production of misfolded proteins which could explain the occurrence of age-related protein aggregation. Nevertheless, some data on this topic is controversial and the biological mechanisms that lead to widespread protein aggregation remain to be elucidated. We review the recent literature about the age-related decline of proteostasis, highlighting the need to build an integrated view of protein synthesis rate, fidelity and quality control pathways in order to better understand the proteome alterations that occur during aging and in age-related diseases.

Regulators of Viral Frameshifting: More Than RNA Influences Translation Events

Programmed ribosomal frameshifting (PRF) is a conserved translational recoding mechanism found in all branches of life and viruses. In bacteria, archaea, and eukaryotes PRF is used to downregulate protein production by inducing a premature termination of translation, which triggers messenger RNA (mRNA) decay. In viruses, PRF is used to drive the production of a new protein while downregulating the production of another protein, thus maintaining a stoichiometry optimal for productive infection. Traditionally, PRF motifs have been defined by the characteristics of two cis elements: a slippery heptanucleotide sequence followed by an RNA pseudoknot or stem-loop within the mRNA. Recently, additional cis and new trans elements have been identified that regulate PRF in both host and viral translation. These additional factors suggest PRF is an evolutionarily conserved process whose function and regulation we are just beginning to understand.

Evidence of multifaceted functions of codon usage in translation within the model beetle Tribolium castaneum

Synonymous codon use is non-random. Codons most used in highly transcribed genes, often called optimal codons, typically have high gene counts of matching tRNA genes (tRNA abundance) and promote accurate and/or efficient translation. Non-optimal codons, those least used in highly expressed genes, may also affect translation. In multicellular organisms, codon optimality may vary among tissues. At present, however, tissue specificity of codon use remains poorly understood. Here, we studied codon usage of genes highly transcribed in germ line (testis and ovary) and somatic tissues (gonadectomized males and females) of the beetle Tribolium castaneum. The results demonstrate that: (i) the majority of optimal codons were organism-wide, the same in all tissues, and had numerous matching tRNA gene copies (Opt-codon_↑tRNAs), consistent with translational selection; (ii) some optimal codons varied among tissues, suggesting tissue-specific tRNA populations; (iii) wobble tRNA were required for translation of certain optimal codons (Opt-codon_wobble), possibly allowing precise translation and/or protein folding; and (iv) remarkably, some non-optimal codons had abundant tRNA genes (Nonopt-codon_↑tRNAs), and genes using those codons were tightly linked to ribosomal and stress-response functions. Thus, Nonopt-codon_↑tRNAs codons may regulate translation of specific genes. Together, the evidence suggests that codon use and tRNA genes regulate multiple translational processes in T. castaneum.

Widespread non-modular overlapping codes in the coding regions

Messenger RNAs (mRNAs) consist of a coding region (open reading frame (ORF)) and two untranslated regions (UTRs), 5'UTR and 3'UTR. Ribosomes travel along the coding region, translating nucleotide triplets (called codons) to a chain of amino acids. The coding region was long believed to mainly encode the amino acid content of proteins, whereas regulatory signals reside in the UTRs and in other genomic regions. However, in recent years we have learned that the ORF is expansively populated with various regulatory signals, or codes, which are related to all gene expression steps and additional intracellular aspects. In this paper, we review the current knowledge related to overlapping codes inside the coding regions, such as the influence of synonymous codon usage on translation speed (and, in turn, the effect of translation speed on protein folding), ribosomal frameshifting, mRNA stability, methylation, splicing, transcription and more. All these codes come together and overlap in the ORF sequence, ensuring production of the right protein at the right time.

Translational efficiency across healthy and tumor tissues is proliferation-related

Different tissues express genes with particular codon usage and anticodon tRNA repertoires. However, the codon-anticodon co-adaptation in humans is not completely understood, nor is its effect on tissue-specific protein levels. Here, we first validated the accuracy of small RNA-seq for tRNA quantification across five human cell lines. We then analyzed the tRNA abundance of more than 8,000 tumor samples from TCGA, together with their paired mRNA-seq and proteomics data, to determine the Supply-to-Demand Adaptation. We thereby elucidate that the dynamic adaptation of the tRNA pool is largely related to the proliferative state across tissues. The distribution of such tRNA pools over the whole cellular translatome affects the subsequent translational efficiency, which functionally determines a condition-specific expression program both in healthy and tumor states. Furthermore, the aberrant translational efficiency of some codons in cancer, exemplified by ProCCA and GlyGGT, is associated with poor patient survival. The regulation of these tRNA profiles is partly explained by the tRNA gene copy numbers and their promoter DNA methylation.

TissueCoCoPUTs: Novel Human Tissue-Specific Codon and Codon-Pair Usage Tables Based on Differential Tissue Gene Expression

Protein expression in multicellular organisms varies widely across tissues. Codon usage in the transcriptome of each tissue is derived from genomic codon usage and the relative expression level of each gene. We created a comprehensive computational resource that houses tissue-specific codon, codon-pair, and dinucleotide usage data for 51 Homo sapiens tissues (TissueCoCoPUTs: https://hive.biochemistry.gwu.edu/review/tissue_codon), using transcriptome data from the Broad Institute Genotype-Tissue Expression (GTEx) portal. Distances between tissue-specific codon and codon-pair frequencies were used to generate a dendrogram based on the unique patterns of codon and codon-pair usage in each tissue that are clearly distinct from the genomic distribution. This novel resource may be useful in unraveling the relationship between codon usage and tRNA abundance, which could be critical in determining translation kinetics and efficiency across tissues. Areas of investigation such as biotherapeutic development, tissue-specific genetic engineering, and genetic disease prediction will greatly benefit from this resource.

Genome-wide Translation Profiling by Ribosome-Bound tRNA Capture

In the ribosome complex, tRNA is a critical element of mRNA translation. A rich repertoire of cell regulation is hypothesized to occur during the recruitment of specific tRNAs in polypeptide formation. However, this basic question in nascent chain biology remains unaddressed due to the lack of technologies to report the complete tRNA complement inside ribosomes during active translation. Here, we characterize a technique for profiling ribosome-embedded tRNA and their modifications. With this method, we generated a comprehensive survey of the quantity and quality of intra-ribosomal tRNAs. In cells under environmental stress, we show that methionine tRNA inside ribosomes is a robust biomarker for the impairment of translation initiation or elongation steps. Concurrent tRNA/mRNA ribosome profiling revealed a stress-dependent incorporation of damaged and uncharged tRNAs into ribosomes causing translation arrest. Thus, tRNA ribosome profiling can provide insights on translation control mechanisms in diverse biological contexts.

Pathological changes are associated with shifts in the employment of synonymous codons at the transcriptome level

BACKGROUND

The usage of different synonymous codons reflects the genome organization and has been connected to parameters such as mRNA abundance and protein folding. It is also been established that mutations targeting specific synonymous codons can trigger disease.

RESULTS

We performed a systematic meta-analysis of transcriptome results from 75 datasets representing 40 pathologies. We found that a subset of codons was preferentially employed in abundant transcripts, while other codons were preferentially found in low-abundance transcripts. By comparing control and pathological transcriptomes, we observed a shift in the employment of synonymous codons for every analyzed disease. For example, cancerous tissue employed preferentially A- or U-ending codons, shifting from G- or C-ending codons, which were preferred by control tissues. This analysis was able to discriminate patients and controls with high specificity and sensitivity.

CONCLUSIONS

Here we show that the employment of specific synonymous codons, quantified at the whole transcriptome level, changes profoundly in many diseases. We propose that the changes in codon employment offer a novel perspective for disease studies, and could be used to design new diagnostic tools.

Translational and transcriptional responses in human primary hepatocytes under hypoxia

The liver is the primary source of a large number of plasma proteins and plays a critical role in multiple biological processes. Inadequate oxygen supply characterizing various clinical settings such as liver transplantation exposes the liver to hypoxic conditions. Studies assessing hypoxia-induced global translational changes in liver are lacking. Here, we employed a recently developed ribosome-profiling technique to assess global translational responses of human primary hepatocytes exposed to acute hypoxic stress (1% O₂) for the short term. In parallel, transcriptome profiling was performed to assess mRNA expression changes. We found that translational responses appeared earlier and were predominant over transcriptional responses. A significant decrease in translational efficiency of several ribosome genes indicated translational inhibition of new ribosome protein synthesis in hypoxia. Pathway enrichment analysis highlighted altered translational regulation of MAPK signaling, drug metabolism, oxidative phosphorylation, and nonalcoholic fatty liver disease pathways. Gene Ontology enrichment analysis revealed terms related to translation, metabolism, angiogenesis, apoptosis, and response to stress. Transcriptional induction of genes encoding heat shock proteins was observed within 30 min of hypoxia. Induction of genes encoding stress response mediators, metabolism regulators, and proangiogenic proteins was observed at 240 min. Despite the liver being the primary source of coagulation proteins and the implicated role of hypoxia in thrombosis, limited differences were observed in genes encoding coagulation-associated proteins. Overall, our study demonstrates the predominance of translational regulation over transcription and highlights differentially regulated pathways or biological processes in short-term hypoxic stress responses of human primary hepatocytes. NEW & NOTEWORTHY The novelty of this study lies in applying parallel ribosome- and transcriptome-profiling analyses to human primary hepatocytes in hypoxia. To our knowledge, this is the first study to assess global translational responses using ribosome profiling in hypoxic hepatocytes. Our results demonstrate the predominance of translational responses over transcriptional responses in early hepatic hypoxic stress responses. Furthermore, our study reveals multiple pathways and specific genes showing altered regulation in hypoxic hepatocytes.

Cells alter their tRNA abundance to selectively regulate protein synthesis during stress conditions

Decoding the information in mRNA during protein synthesis relies on tRNA adaptors, the abundance of which can affect the decoding rate and translation efficiency. To determine whether cells alter tRNA abundance to selectively regulate protein expression, we quantified changes in the abundance of individual tRNAs at different time points in response to diverse stress conditions in Saccharomyces cerevisiae. We found that the tRNA pool was dynamic and rearranged in a manner that facilitated selective translation of stress-related transcripts. Through genomic analysis of multiple data sets, stochastic simulations, and experiments with designed sequences of proteins with identical amino acids but altered codon usage, we showed that changes in tRNA abundance affected protein expression independently of factors such as mRNA abundance. We suggest that cells alter their tRNA abundance to selectively affect the translation rates of specific transcripts to increase the amounts of required proteins under diverse stress conditions.

Coping with stress by regulating tRNAs

Stress conditions curtail the energetically costly process of messenger RNA translation. In this issue of Science Signaling, Torrent et al. report key evidence for a direct link between codon usage and translation regulation in response to stress.

Repertoires of tRNAs: The Couplers of Genomics and Proteomics

The pool of transfer RNA (tRNA) molecules in cells allows the ribosome to decode genetic information. This repertoire of molecular decoders is positioned in the crossroad of the genome, the transcriptome, and the proteome. Omics and systems biology now allow scientists to explore the entire repertoire of tRNAs of many organisms, revealing basic exciting biology. The tRNA gene set of hundreds of species is now characterized, in addition to the tRNA genes of organelles and viruses. Genes encoding tRNAs for certain anticodon types appear in dozens of copies in a genome, while others are universally absent from any genome. Transcriptome measurement of tRNAs is challenging, but in recent years new technologies have allowed researchers to determine the dynamic expression patterns of tRNAs. These advances reveal that availability of ready-to-translate tRNA molecules is highly controlled by several transcriptional and posttranscriptional regulatory processes. This regulation shapes the proteome according to the cellular state. The tRNA pool profoundly impacts many aspects of cellular and organismal life, including protein expression level, translation accuracy, adequacy of folding, and even mRNA stability. As a result, the shape of the tRNA pool affects organismal health and may participate in causing conditions such as cancer and neurological conditions.

Analysis of codon usage patterns in Hirudinaria manillensis reveals a preference for GC-ending codons caused by dominant selection constraints

Background

Hirudinaria manillensis is an ephemeral, blood-sucking ectoparasite, possessing anticoagulant capacities with potential medical applications. Analysis of codon usage patterns would contribute to our understanding of the evolutionary mechanisms and genetic architecture of H. manillensis, which in turn would provide insight into the characteristics of other leeches. We analysed codon usage and related indices using 18,000 coding sequences (CDSs) retrieved from H. manillensis RNA-Seq data.

Results

We identified four highly preferred codons in H. manillensis that have G/C-endings. Points generated in an effective number of codons (ENC) plot distributed below the standard curve and the slope of a neutrality plot was less than 1. Highly expressed CDSs had lower ENC content and higher GC content than weakly expressed CDSs. Principal component analysis conducted on relative synonymous codon usage (RSCU) values divided CDSs according to GC content and divided codons according to ending bases. Moreover, by determining codon usage, we found that the majority of blood-diet related genes have undergone less adaptive evolution in H. manillensis, except for those with homologous sequences in the host species.

Conclusions

Codon usage in H. manillensis had an overall preference toward C-endings and indicated that codon usage patterns are mediated by differential expression, GC content, and biological function. Although mutation pressure effects were also notable, the majority of genetic evolution in H. manillensis was driven by natural selection.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4937-x) contains supplementary material, which is available to authorized users.

Functionality of tRNAs encoded in a mobile genetic element from an acidophilic bacterium

The genome of the acidophilic, bioleaching bacterium Acidithiobacillus ferrooxidans, strain ATCC 23270, contains 95 predicted tRNA genes. Thirty-six of these genes (all 20 species) are clustered within an actively excising integrative-conjugative element (ICEAfe1). We speculated that these tRNA genes might have a role in adapting the bacterial tRNA pool to the codon usage of ICEAfe1 genes. To answer this question, we performed theoretical calculations of the global tRNA adaptation index to the entire A. ferrooxidans genome with and without the ICEAfe1 encoded tRNA genes. Based on these calculations, we observed that tRNAs encoded in ICEAfe1 negatively contribute to adapt the tRNA pool to the codon use in A. ferrooxidans. Although some of the tRNAs encoded in ICEAfe1 are functional in aminoacylation or protein synthesis, we found that they are expressed at low levels. These findings, along with the identification of a tRNA-like RNA encoded in the same cluster, led us to speculate that tRNA genes encoded in the mobile genetic element ICEAfe1 might have acquired mutations that would result in either inactivation or the acquisition of new functions.

Beyond the Triplet Code: Context Cues Transform Translation

The elucidation of the genetic code remains among the most influential discoveries in biology. While innumerable studies have validated the general universality of the code and its value in predicting and analyzing protein coding sequences, established and emerging work has also suggested that full genome decryption may benefit from a greater consideration of a codon's neighborhood within an mRNA than has been broadly applied. This Review examines the evidence for context cues in translation, with a focus on several recent studies that reveal broad roles for mRNA context in programming translation start sites, the rate of translation elongation, and stop codon identity.

tRNA-mediated codon-biased translation in mycobacterial hypoxic persistence

Microbial pathogens adapt to the stress of infection by regulating transcription, translation and protein modification. We report that changes in gene expression in hypoxia-induced non-replicating persistence in mycobacteria—which models tuberculous granulomas—are partly determined by a mechanism of tRNA reprogramming and codon-biased translation. Mycobacterium bovis BCG responded to each stage of hypoxia and aerobic resuscitation by uniquely reprogramming 40 modified ribonucleosides in tRNA, which correlate with selective translation of mRNAs from families of codon-biased persistence genes. For example, early hypoxia increases wobble cmo⁵U in tRNA^Thr(UGU), which parallels translation of transcripts enriched in its cognate codon, ACG, including the DosR master regulator of hypoxic bacteriostasis. Codon re-engineering of dosR exaggerates hypoxia-induced changes in codon-biased DosR translation, with altered dosR expression revealing unanticipated effects on bacterial survival during hypoxia. These results reveal a coordinated system of tRNA modifications and translation of codon-biased transcripts that enhance expression of stress response proteins in mycobacteria.

Mycobacteria can adapt to the stress of human infection by entering a dormant state. Here the authors show that hypoxia-induced dormancy in M. bovis BCG involves the reprogramming of tRNA wobble modifications and copy numbers, coupled with biased use of synonymous codons in survival genes.

Dynamic translation regulation in Caulobacter cell cycle control

Progression of the Caulobacter cell cycle requires temporal and spatial control of gene expression, culminating in an asymmetric cell division yielding distinct daughter cells. To explore the contribution of translational control, RNA-seq and ribosome profiling were used to assay global transcription and translation levels of individual genes at six times over the cell cycle. Translational efficiency (TE) was used as a metric for the relative rate of protein production from each mRNA. TE profiles with similar cell cycle patterns were found across multiple clusters of genes, including those in operons or in subsets of operons. Collections of genes associated with central cell cycle functional modules (e.g., biosynthesis of stalk, flagellum, or chemotaxis machinery) have consistent but different TE temporal patterns, independent of their operon organization. Differential translation of operon-encoded genes facilitates precise cell cycle-timing for the dynamic assembly of multiprotein complexes, such as the flagellum and the stalk and the correct positioning of regulatory proteins to specific cell poles. The cell cycle-regulatory pathways that produce specific temporal TE patterns are separate from-but highly coordinated with-the transcriptional cell cycle circuitry, suggesting that the scheduling of translational regulation is organized by the same cyclical regulatory circuit that directs the transcriptional control of the Caulobacter cell cycle.

Tissue- and Time-Specific Expression of Otherwise Identical tRNA Genes

Codon usage bias affects protein translation because tRNAs that recognize synonymous codons differ in their abundance. Although the current dogma states that tRNA expression is exclusively regulated by intrinsic control elements (A- and B-box sequences), we revealed, using a reporter that monitors the levels of individual tRNA genes in Caenorhabditis elegans, that eight tryptophan tRNA genes, 100% identical in sequence, are expressed in different tissues and change their expression dynamically. Furthermore, the expression levels of the sup-7 tRNA gene at day 6 were found to predict the animal's lifespan. We discovered that the expression of tRNAs that reside within introns of protein-coding genes is affected by the host gene's promoter. Pairing between specific Pol II genes and the tRNAs that are contained in their introns is most likely adaptive, since a genome-wide analysis revealed that the presence of specific intronic tRNAs within specific orthologous genes is conserved across Caenorhabditis species.

Overcoming Challenges in Engineering the Genetic Code

Withstanding 3.5 billion years of genetic drift, the canonical genetic code remains such a fundamental foundation for the complexity of life that it is highly conserved across all three phylogenetic domains. Genome engineering technologies are now making it possible to rationally change the genetic code, offering resistance to viruses, genetic isolation from horizontal gene transfer, and prevention of environmental escape by genetically modified organisms. We discuss the biochemical, genetic, and technological challenges that must be overcome in order to engineer the genetic code.

A relay race on the evolutionary adaptation spectrum

Adaptation is the process in which organisms improve their fitness by changing their phenotype using genetic or non-genetic mechanisms. The adaptation toolbox consists of varied molecular and genetic means that we posit span an almost continuous "adaptation spectrum." Different adaptations are characterized by the time needed for organisms to attain them and by their duration. We suggest that organisms often adapt by progressing the adaptation spectrum, starting with rapidly attained physiological and epigenetic adaptations and culminating with slower long-lasting genetic ones. A tantalizing possibility is that earlier adaptations facilitate realization of later ones.

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.

Why should cancer biologists care about tRNAs? tRNA synthesis, mRNA translation and the control of growth

Transfer RNAs (tRNAs) are essential for mRNA translation. They are transcribed in the nucleus by RNA polymerase III and undergo many modifications before contributing to cytoplasmic protein synthesis. In this review I highlight our understanding of how tRNA biology may be linked to the regulation of mRNA translation, growth and tumorigenesis. First, I review how oncogenes and tumour suppressor signalling pathways, such as the PI3 kinase/TORC1, Ras/ERK, Myc, p53 and Rb pathways, regulate Pol III and tRNA synthesis. In several cases, this regulation contributes to cell, tissue and body growth, and has implications for our understanding of tumorigenesis. Second, I highlight some recent work, particularly in model organisms such as yeast and Drosophila, that shows how alterations in tRNA synthesis may be not only necessary, but also sufficient to drive changes in mRNA translation and growth. These effects may arise due to both absolute increases in total tRNA levels, but also changes in the relative levels of tRNAs in the overall pool. Finally, I review some recent studies that have revealed how tRNA modifications (amino acid acylation, base modifications, subcellular shuttling, and cleavage) can be regulated by growth and stress cues to selectively influence mRNA translation. Together these studies emphasize the importance of the regulation of tRNA synthesis and modification as critical control points in protein synthesis and growth. This article is part of a Special Issue entitled: Translation and Cancer.

Metabolic de-isotoping for improved LC-MS characterization of modified RNAs

Mapping, sequencing, and quantifying individual noncoding ribonucleic acids (ncRNAs), including post-transcriptionally modified nucleosides, by mass spectrometry is a challenge that often requires rigorous sample preparation prior to analysis. Previously, we have described a simplified method for the comparative analysis of RNA digests (CARD) that is applicable to relatively complex mixtures of ncRNAs. In the CARD approach for transfer RNA (tRNA) analysis, two complete sets of digestion products from total tRNA are compared using the enzymatic incorporation of (16)O/(18)O isotopic labels. This approach allows one to rapidly screen total tRNAs from gene deletion mutants or comparatively sequence total tRNA from two related bacterial organisms. However, data analysis can be challenging because of convoluted mass spectra arising from the natural (13)C and (15) N isotopes present in the ribonuclease-digested tRNA samples. Here, we demonstrate that culturing in (12)C-enriched/(13)C-depleted media significantly reduces the isotope patterns that must be interpreted during the CARD experiment. Improvements in data quality yield a 35 % improvement in detection of tRNA digestion products that can be uniquely assigned to particular tRNAs. These mass spectral improvements lead to a significant reduction in data processing attributable to the ease of spectral identification of labeled digestion products and will enable improvements in the relative quantification of modified RNAs by the (16)O/(18)O differential labeling approach.

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.

A comprehensive tRNA deletion library unravels the genetic architecture of the tRNA pool

Deciphering the architecture of the tRNA pool is a prime challenge in translation research, as tRNAs govern the efficiency and accuracy of the process. Towards this challenge, we created a systematic tRNA deletion library in Saccharomyces cerevisiae, aimed at dissecting the specific contribution of each tRNA gene to the tRNA pool and to the cell's fitness. By harnessing this resource, we observed that the majority of tRNA deletions show no appreciable phenotype in rich medium, yet under more challenging conditions, additional phenotypes were observed. Robustness to tRNA gene deletion was often facilitated through extensive backup compensation within and between tRNA families. Interestingly, we found that within tRNA families, genes carrying identical anti-codons can contribute differently to the cellular fitness, suggesting the importance of the genomic surrounding to tRNA expression. Characterization of the transcriptome response to deletions of tRNA genes exposed two disparate patterns: in single-copy families, deletions elicited a stress response; in deletions of genes from multi-copy families, expression of the translation machinery increased. Our results uncover the complex architecture of the tRNA pool and pave the way towards complete understanding of their role in cell physiology.

tRNA genes rapidly change in evolution to meet novel translational demands

Changes in expression patterns may occur when organisms are presented with new environmental challenges, for example following migration or genetic changes. To elucidate the mechanisms by which the translational machinery adapts to such changes, we perturbed the tRNA pool of Saccharomyces cerevisiae by tRNA gene deletion. We then evolved the deletion strain and observed that the genetic adaptation was recurrently based on a strategic mutation that changed the anticodon of other tRNA genes to match that of the deleted one. Strikingly, a systematic search in hundreds of genomes revealed that anticodon mutations occur throughout the tree of life. We further show that the evolution of the tRNA pool also depends on the need to properly couple translation to protein folding. Together, our observations shed light on the evolution of the tRNA pool, demonstrating that mutation in the anticodons of tRNA genes is a common adaptive mechanism when meeting new translational demands.

DOI:http://dx.doi.org/10.7554/eLife.01339.001

Genes contain the blueprints for the proteins that are essential for countless biological functions and processes, and the path that leads from a particular gene to the corresponding protein is long and complex. The genetic information stored in the DNA must first be transcribed to produce a messenger RNA molecule, which then has to be translated to produce a string of amino acids that fold to form a protein. The translation step is performed by a molecular machine called the ribosome, with transfer RNA molecules bringing the amino acids that are needed to make the protein.

The information in messenger RNA is stored as a series of letters, with groups of three letters called codons representing the different amino acids. Since there are four letters—A, C, G and U—it is possible to form 64 different codons. And since there are only 20 amino acids, two or more different codons can specify the same amino acid (for example, AGU and AGC both specify serine), and two or more different transfer RNA molecules can take this amino acid to the ribosome. Moreover, some codons are found more often than others in the messenger RNA molecules, so the genes that encode the related transfer RNA molecules are more common than the genes for other transfer RNA molecules.

Environmental pressures mean that organisms must adapt to survive, with some genes and proteins increasing in importance, and others becoming less important. Clearly the relative numbers of the different transfer RNA molecules will also need to change to reflect these evolutionary changes, but the details of how this happens were not understood.

Now Yona et al. have explored this issue by studying yeast cells that lack a gene for one of the less common transfer RNA molecules (corresponding to the codon AGG, which specifies the amino acid arginine). At first this mutation resulted in slower growth of the yeast cells, but after being allowed to evolve over 200 generations, the rate of growth matched that of a normal strain with all transfer RNA genes. Yona et al. found that the gene for a more common transfer RNA molecule, corresponding to the codon AGA, which also specifies arginine, had mutated to AGG. As a result, the mutated yeast was eventually able to produce proteins as quickly as wild type yeast. Moreover, further experiments showed that the levels of some transfer RNAs are kept deliberately low in order to slow down the production of proteins so as to ensure that the proteins assume their correct structure.

But does the way these cells evolved in the lab resemble what happened in nature? To address this question Yona et al. examined a database of transfer RNA sequences from more than 500 species, and found evidence for the same codon-based switching mechanism in many species across the tree of life.

DOI:http://dx.doi.org/10.7554/eLife.01339.002

Changes in translational efficiency is a dominant regulatory mechanism in the environmental response of bacteria

To understand how cell physiological state affects mRNA translation, we used Shewanella oneidensis MR-1 grown under steady state conditions at either 20% or 8.5% O2. Using a combination of quantitative proteomics and RNA-Seq, we generated high-confidence data on >1000 mRNA and protein pairs. By using a steady state model, we found that differences in protein-mRNA ratios were primarily due to differences in the translational efficiency of specific genes. When oxygen levels were lowered, 28% of the proteins showed at least a 2-fold change in expression. Transcription levels were sp. significantly altered for 26% of the protein changes; translational efficiency was significantly altered for 46% and a combination of both was responsible for the remaining 28%. Changes in translational efficiency were significantly correlated with the codon usage pattern of the genes and measurable tRNA pools changed in response to altered O2 levels. Our results suggest that changes in the translational efficiency of proteins, in part due to altered tRNA pools, is a major determinant of regulated alterations in protein expression levels in bacteria.

Genome-wide analysis of selective constraints on high stability regions of mRNA reveals multiple compensatory mutations in Escherichia coli

Message RNA (mRNA) carries a large number of local secondary structures, with structural stability to participate in the regulations of gene expression. A worthy question is how the local structural stability is maintained under the constraint that multiple selective pressures are imposed on mRNA local regions. Here, we performed the first genome-wide study of natural selection operating on high structural stability regions (HSRs) of mRNAs in Escherichia coli. We found that HSR tends to adjust the folded conformation to reduce the harm of mutations, showing a high level of mutational robustness. Moreover, guanine preference in HSR was observed, supporting the hypothesis that the selective constraint for high structural stability may partly account for the high percentage of G content in Escherichia coli genome. Notably, we found a substantially reduced synonymous substitution rate in HSRs compared with that in their adjacent regions. Surprisingly and interestingly, the non-key sites in HSRs, which have slight effect on structural stability, have synonymous substitution rate equivalent to background regions. To explain this result, we identified compensatory mutations in HSRs based on structural stability, and found that a considerable number of synonymous mutations occur to restore the structural stability decreased heavily by the mutations on key sites. Overall, these results suggest a significant role of local structural stability as a selective force operating on mRNA, which furthers our understanding of the constraints imposed on protein-coding RNAs.

Translational sensitivity of the Escherichia coli genome to fluctuating tRNA availability

The synthesis of protein from messenger RNA during translation is a highly dynamic process that plays a key role in controlling the efficiency and fidelity of genome-wide protein expression. The availability of aminoacylated transfer RNA (tRNA) is a major factor influencing the speed of ribosomal movement, which depending on codon choices, varies considerably along a transcript. Furthermore, it has been shown experimentally that tRNA availability can vary significantly under different growth and stress conditions, offering the cell a way to adapt translational dynamics across the genome. Existing models of translation have neglected fluctuations of tRNA pools, instead assuming fixed tRNA availabilities over time. This has lead to an incomplete understanding of this process. Here, we show for the entire Escherichia coli genome how and to what extent translational speed profiles, which capture local aspects of translational elongation, respond to measured shifts in tRNA availability. We find that translational profiles across the genome are affected to differing degrees, with genes that are essential or related to fundamental processes such as translation, being more robust than those linked to regulation. Furthermore, we reveal how fluctuating tRNA availability influences profiles of specific sequences known to play a significant role in translational control of gene expression.