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

Translation Rates and Protein Folding

The mRNA coding sequence defines not only the amino acid sequence of the protein, but also the speed at which the ribosomes move along the mRNA while making the protein. The non-uniform local kinetics - denoted as translational rhythm - is similar among mRNAs coding for related protein folds. Deviations from this conserved rhythm can result in protein misfolding. In this review we summarize the experimental evidence demonstrating how local translation rates affect cotranslational protein folding, with the focus on the synonymous codons and patches of charged residues in the nascent peptide as best-studied examples. Alterations in nascent protein conformations due to disturbed translational rhythm can persist off the ribosome, as demonstrated by the effects of synonymous codon variants of several disease-related proteins. Charged amino acid patches in nascent chains also modulate translation and cotranslational protein folding, and can abrogate translation when placed at the N-terminus of the nascent peptide. During cotranslational folding, incomplete nascent chains navigate through a unique conformational landscape in which earlier intermediate states become inaccessible as the nascent peptide grows. Precisely tuned local translation rates, as well as interactions with the ribosome, guide the folding pathway towards the native structure, whereas deviations from the natural translation rhythm may favor pathways leading to trapped misfolded states. Deciphering the 'folding code' of the mRNA will contribute to understanding the diseases caused by protein misfolding and to rational protein design.

The Giardia lamblia ribosome structure reveals divergence in several biological pathways and the mode of emetine function

Giardia lamblia is a deeply branching protist and a human pathogen. Its unusual biology presents the opportunity to explore conserved and fundamental molecular mechanisms. We determined the structure of the G. lamblia 80S ribosome bound to tRNA, mRNA, and the antibiotic emetine by cryo-electron microscopy, to an overall resolution of 2.49 Å. The structure reveals rapidly evolving protein and nucleotide regions, differences in the peptide exit tunnel, and likely altered ribosome quality control pathways. Examination of translation initiation factor binding sites suggests these interactions are conserved despite a divergent initiation mechanism. Highlighting the potential of G. lamblia to resolve conserved biological principles; our structure reveals the interactions of the translation inhibitor emetine with the ribosome and mRNA, thus providing insight into the mechanism of action for this widely used antibiotic. Our work defines key questions in G. lamblia and motivates future experiments to explore the diversity of eukaryotic gene regulation.

Ribosome profiling: a powerful tool in oncological research

Neoplastic cells need to adapt their gene expression pattern to survive in an ever-changing or unfavorable tumor microenvironment. Protein synthesis (or mRNA translation), an essential part of gene expression, is dysregulated in cancer. The emergence of distinct translatomic technologies has revolutionized oncological studies to elucidate translational regulatory mechanisms. Ribosome profiling can provide adequate information on diverse aspects of translation by aiding in quantitatively analyzing the intensity of translating ribosome-protected fragments. Here, we review the primary currently used translatomics techniques and highlight their advantages and disadvantages as tools for translatomics studies. Subsequently, we clarified the areas in which ribosome profiling could be applied to better understand translational control. Finally, we summarized the latest advances in cancer studies using ribosome profiling to highlight the extensive application of this powerful and promising translatomic tool.

Streamlined and sensitive mono- and di-ribosome profiling in yeast and human cells

Ribosome profiling has unveiled diverse regulation and perturbations of translation through a transcriptome-wide survey of ribosome occupancy, read out by sequencing of ribosome-protected messenger RNA fragments. Generation of ribosome footprints and their conversion into sequencing libraries is technically demanding and sensitive to biases that distort the representation of physiological ribosome occupancy. We address these challenges by producing ribosome footprints with P1 nuclease rather than RNase I and replacing RNA ligation with ordered two-template relay, a single-tube protocol for sequencing library preparation that incorporates adaptors by reverse transcription. Our streamlined approach reduced sequence bias and enhanced enrichment of ribosome footprints relative to ribosomal RNA. Furthermore, P1 nuclease preserved distinct juxtaposed ribosome complexes informative about yeast and human ribosome fates during translation initiation, stalling and termination. Our optimized methods for mRNA footprint generation and capture provide a richer translatome profile with low input and fewer technical challenges.

choros: correction of sequence-based biases for accurate quantification of ribosome profiling data

Ribosome profiling quantifies translation genome-wide by sequencing ribosome-protected fragments, or footprints. Its single-codon resolution allows identification of translation regulation, such as ribosome stalls or pauses, on individual genes. However, enzyme preferences during library preparation lead to pervasive sequence artifacts that obscure translation dynamics. Widespread over- and under-representation of ribosome footprints can dominate local footprint densities and skew estimates of elongation rates by up to five fold. To address these biases and uncover true patterns of translation, we present choros, a computational method that models ribosome footprint distributions to provide bias-corrected footprint counts. choros uses negative binomial regression to accurately estimate two sets of parameters: (i) biological contributions from codon-specific translation elongation rates; and (ii) technical contributions from nuclease digestion and ligation efficiencies. We use these parameter estimates to generate bias correction factors that eliminate sequence artifacts. Applying choros to multiple ribosome profiling datasets, we are able to accurately quantify and attenuate ligation biases to provide more faithful measurements of ribosome distribution. We show that a pattern interpreted as pervasive ribosome pausing near the beginning of coding regions is likely to arise from technical biases. Incorporating choros into standard analysis pipelines will improve biological discovery from measurements of translation.

Oxygen level regulates N-terminal translation elongation of selected proteins through deoxyhypusine hydroxylation

Hypusine is a post-translational modification on eukaryotic translation initiation factor 5A (eIF5A). The last step of hypusine biosynthesis, deoxyhypusine hydroxylation, is an oxygen-dependent reaction. Here we show that deletion of the deoxyhypusine hydroxylase Lia1 compromises yeast respiration through translation downregulation of selected proteins in the respiration pathway. The translation suppression, because of the lack of deoxyhypusine hydroxylation, mainly affects translation of the N termini of the proteins, independent of the presence of proline residues but likely dependent on the interaction between the N-terminal nascent peptide and the ribosomal peptide exit tunnel. Proteomics and biochemical studies reveal that Lia1 deletion decreases N-terminal translation of proteins involved in mitochondrial respiration, oxidative stress response, and protein folding. Our work uncovers functions of the hypusine modification by considering the substrate requirement of the post-translational modification, highlights the unique challenges of translating the N termini of proteins, and reveals an oxygen-sensing mechanism in eukaryotic cells.

Peptidyl transferase center decompaction and structural constraints during early protein elongation on the ribosome

Peptide bond formation on the ribosome requires that aminoacyl-tRNAs and peptidyl-tRNAs are properly positioned on the A site and the P site of the peptidyl transferase center (PTC) so that nucleophilic attack can occur. Here we analyse some constraints associated with the induced-fit mechanism of the PTC, that promotes this positioning through a compaction around the aminoacyl ester orchestrated by U2506. The physical basis of PTC decompaction, that allows the elongated peptidyl-tRNA to free itself from that state and move to the P site of the PTC, is still unclear. From thermodynamics considerations and an analysis of published ribosome structures, the present work highlights the rational of this mechanism, in which the free-energy released by the new peptide bond is used to kick U2506 away from the reaction center. Furthermore, we show the evidence that decompaction is impaired when the nascent peptide is not yet anchored inside the exit tunnel, which may contribute to explain why the first rounds of elongation are inefficient, an issue that has attracted much interest for about two decades. Results in this field are examined in the light of the present analysis and a physico-chemical correlation in the genetic code, which suggest that elementary constraints associated with the size of the side-chain of the amino acids penalize early elongation events.

Translation complex stabilization on messenger RNA and footprint profiling to study the RNA responses and dynamics of protein biosynthesis in the cells

During protein biosynthesis, ribosomes bind to messenger (m)RNA, locate its protein-coding information, and translate the nucleotide triplets sequentially as codons into the corresponding sequence of amino acids, forming proteins. Non-coding mRNA features, such as 5' and 3' untranslated regions (UTRs), start sites or stop codons of different efficiency, stretches of slower or faster code and nascent polypeptide interactions can alter the translation rates transcript-wise. Most of the homeostatic and signal response pathways of the cells converge on individual mRNA control, as well as alter the global translation output. Among the multitude of approaches to study translational control, one of the most powerful is to infer the locations of translational complexes on mRNA based on the mRNA fragments protected by these complexes from endonucleolytic hydrolysis, or footprints. Translation complex profiling by high-throughput sequencing of the footprints allows to quantify the transcript-wise, as well as global, alterations of translation, and uncover the underlying control mechanisms by attributing footprint locations and sizes to different configurations of the translational complexes. The accuracy of all footprint profiling approaches critically depends on the fidelity of footprint generation and many methods have emerged to preserve certain or multiple configurations of the translational complexes, often in challenging biological material. In this review, a systematic summary of approaches to stabilize translational complexes on mRNA for footprinting is presented and major findings are discussed. Future directions of translation footprint profiling are outlined, focusing on the fidelity and accuracy of inference of the native in vivo translation complex distribution on mRNA.

SMN-primed ribosomes modulate the translation of transcripts related to spinal muscular atrophy

The contribution of ribosome heterogeneity and ribosome-associated proteins to the molecular control of proteomes in health and disease remains unclear. Here, we demonstrate that survival motor neuron (SMN) protein-the loss of which causes the neuromuscular disease spinal muscular atrophy (SMA)-binds to ribosomes and that this interaction is tissue-dependent. SMN-primed ribosomes are preferentially positioned within the first five codons of a set of mRNAs that are enriched for translational enhancer sequences in the 5' untranslated region (UTR) and rare codons at the beginning of their coding sequence. These SMN-specific mRNAs are associated with neurogenesis, lipid metabolism, ubiquitination, chromatin regulation and translation. Loss of SMN induces ribosome depletion, especially at the beginning of the coding sequence of SMN-specific mRNAs, leading to impairment of proteins that are involved in motor neuron function and stability, including acetylcholinesterase. Thus, SMN plays a crucial role in the regulation of ribosome fluxes along mRNAs encoding proteins that are relevant to SMA pathogenesis.

Nutrient Control of mRNA Translation

The emergence of genome-wide analyses to interrogate cellular DNA, RNA, and protein content has revolutionized the study of control networks that mediate cellular homeostasis. mRNA translation represents the last step of genetic flow and primarily defines the proteome. Translational regulation is thus critical for gene expression, in particular under nutrient excess or deficiency. Until recently, it was unclear how the global effects of translational control are orchestrated by nutrient signaling pathways. An emerging concept of translational reprogramming addresses how to maintain the expression of specific proteins during nutrient stress by translation of selective mRNAs. In this review, we describe recent advances in our understanding of translational control principles; nutrient-sensing mechanisms; and their dysregulation in human diseases such as diabetes, cancer, and aging. The mechanistic understanding of translational regulation in response to different nutrient conditions may help identify potential dietary and therapeutic targets to improve human health.

Transcriptome-wide sites of collided ribosomes reveal principles of translational pausing

Translation initiation is the major regulatory step defining the rate of protein production from an mRNA. Meanwhile, the impact of nonuniform ribosomal elongation rates is largely unknown. Using a modified ribosome profiling protocol based on footprints from two closely packed ribosomes (disomes), we have mapped ribosomal collisions transcriptome-wide in mouse liver. We uncover that the stacking of an elongating onto a paused ribosome occurs frequently and scales with translation rate, trapping ∼10% of translating ribosomes in the disome state. A distinct class of pause sites is indicative of deterministic pausing signals. Pause site association with specific amino acids, peptide motifs, and nascent polypeptide structure is suggestive of programmed pausing as a widespread mechanism associated with protein folding. Evolutionary conservation at disome sites indicates functional relevance of translational pausing. Collectively, our disome profiling approach allows unique insights into gene regulation occurring at the step of translation elongation.

Hallmarks and Determinants of Oncogenic Translation Revealed by Ribosome Profiling in Models of Breast Cancer

Gene expression is extensively and dynamically modulated at the level of translation. How cancer cells prioritize the translation of certain mRNAs over others from a pool of competing mRNAs remains an open question. Here, we analyze translation in cell line models of breast cancer and normal mammary tissue by ribosome profiling. We identify key recurrent themes of oncogenic translation: higher ribosome occupancy, greater variance of translational efficiencies, and preferential translation of transcriptional regulators and signaling proteins in malignant cells as compared with their nonmalignant counterpart. We survey for candidate RNA interacting proteins that could associate with the 5′untranslated regions of the transcripts preferentially translated in breast tumour cells. We identify SRSF1, a prototypic splicing factor, to have a pervasive direct and indirect impact on translation. In a representative estrogen receptor–positive and estrogen receptor–negative cell line, we find that protein synthesis relies heavily on SRSF1. SRSF1 is predominantly intranuclear. Under certain conditions, SRSF1 translocates from the nucleus to the cytoplasm where it associates with MYC and CDK1 mRNAs and upregulates their internal ribosome entry site–mediated translation. Our results point to a synergy between splicing and translation and unveil how certain RNA-binding proteins modulate the translational landscape in breast cancer.

A short translational ramp determines the efficiency of protein synthesis

Translation initiation is a major rate-limiting step for protein synthesis. However, recent studies strongly suggest that the efficiency of protein synthesis is additionally regulated by multiple factors that impact the elongation phase. To assess the influence of early elongation on protein synthesis, we employed a library of more than 250,000 reporters combined with in vitro and in vivo protein expression assays. Here we report that the identity of the amino acids encoded by codons 3 to 5 impact protein yield. This effect is independent of tRNA abundance, translation initiation efficiency, or overall mRNA structure. Single-molecule measurements of translation kinetics revealed pausing of the ribosome and aborted protein synthesis on codons 4 and 5 of distinct amino acid and nucleotide compositions. Finally, introduction of preferred sequence motifs only at specific codon positions improves protein synthesis efficiency for recombinant proteins. Collectively, our data underscore the critical role of early elongation events in translational control of gene expression.

Several factors contribute to the efficiency of protein expression. Here the authors show that the identity of amino acids encoded by codons at position 3–5 significantly impact translation efficiency and protein expression levels.

From reporters to endogenous genes: the impact of the first five codons on translation efficiency in Escherichia coli

Translation initiation is a critical step in the regulation of protein synthesis, and it is subjected to different control mechanisms, such as 5' UTR secondary structure and initiation codon context, that can influence the rates at which initiation and consequentially translation occur. For some genes, translation elongation also affects the rate of protein synthesis. With a GFP library containing nearly all possible combinations of nucleotides from the 3rd to the 5th codon positions in the protein coding region of the mRNA, it was previously demonstrated that some nucleotide combinations increased GFP expression up to four orders of magnitude. While it is clear that the codon region from positions 3 to 5 can influence protein expression levels of artificial constructs, its impact on endogenous proteins is still unknown. Through bioinformatics analysis, we identified the nucleotide combinations of the GFP library in Escherichia coli genes and examined the correlation between the expected levels of translation according to the GFP data with the experimental measures of protein expression. We observed that E. coli genes were enriched with the nucleotide compositions that enhanced protein expression in the GFP library, but surprisingly, it seemed to affect the translation efficiency only marginally. Nevertheless, our data indicate that different enterobacteria present similar nucleotide composition enrichment as E. coli, suggesting an evolutionary pressure towards the conservation of short translational enhancer sequences.

Decomposition of time-dependent fluorescence signals reveals codon-specific kinetics of protein synthesis

During protein synthesis, the nascent peptide chain traverses the peptide exit tunnel of the ribosome. We monitor the co-translational movement of the nascent peptide using a fluorescent probe attached to the N-terminus of the nascent chain. Due to fluorophore quenching, the time-dependent fluorescence signal emitted by an individual peptide is determined by co-translational events, such as secondary structure formation and peptide-tunnel interactions. To obtain information on these individual events, the measured ensemble fluorescence signal has to be decomposed into position-dependent intensities. Here, we describe mRNA translation as a Markov process with specific fluorescence intensities assigned to the different states of the process. Combining the computed stochastic time evolution of the translation process with a sequence of observed ensemble fluorescence time courses, we compute the unknown position-specific intensities and obtain detailed information on the kinetics of the translation process. In particular, we find that translation of poly(U) mRNAs dramatically slows down at the fourth UUU codon. The method presented here detects subtle differences in the position-specific fluorescence intensities and thus provides a novel approach to study translation kinetics in ensemble experiments.

Ribosome Profiling: Global Views of Translation

The translation of messenger RNA (mRNA) into protein and the folding of the resulting protein into an active form are prerequisites for virtually every cellular process and represent the single largest investment of energy by cells. Ribosome profiling-based approaches have revolutionized our ability to monitor every step of protein synthesis in vivo, allowing one to measure the rate of protein synthesis across the proteome, annotate the protein coding capacity of genomes, monitor localized protein synthesis, and explore cotranslational folding and targeting. The rich and quantitative nature of ribosome profiling data provides an unprecedented opportunity to explore and model complex cellular processes. New analytical techniques and improved experimental protocols will provide a deeper understanding of the factors controlling translation speed and its impact on protein function and cell physiology as well as the role of ribosomal RNA and mRNA modifications in regulating translation.

Identifying A- and P-site locations on ribosome-protected mRNA fragments using Integer Programming

Identifying the A- and P-site locations on ribosome-protected mRNA fragments from Ribo-Seq experiments is a fundamental step in the quantitative analysis of transcriptome-wide translation properties at the codon level. Many analyses of Ribo-Seq data have utilized heuristic approaches applied to a narrow range of fragment sizes to identify the A-site. In this study, we use Integer Programming to identify the A-site by maximizing an objective function that reflects the fact that the ribosome's A-site on ribosome-protected fragments must reside between the second and stop codons of an mRNA. This identifies the A-site location as a function of the fragment's size and its 5' end reading frame in Ribo-Seq data generated from S. cerevisiae and mouse embryonic stem cells. The correctness of the identified A-site locations is demonstrated by showing that this method, as compared to others, yields the largest ribosome density at established stalling sites. By providing greater accuracy and utilization of a wider range of fragment sizes, our approach increases the signal-to-noise ratio of underlying biological signals associated with translation elongation at the codon length scale.

SENP3-mediated host defense response contains HBV replication and restores protein synthesis

Certain organs are capable of containing the replication of various types of viruses. In the liver, infection of Hepatitis B virus (HBV), the etiological factor of Hepatitis B and hepatocellular carcinoma (HCC), often remains asymptomatic and leads to a chronic carrier state. Here we investigated how hepatocytes contain HBV replication and promote their own survival by orchestrating a translational defense mechanism via the stress-sensitive SUMO-2/3-specific peptidase SENP3. We found that SENP3 expression level decreased in HBV-infected hepatocytes in various models including HepG2-NTCP cell lines and a humanized mouse model. Downregulation of SENP3 reduced HBV replication and boosted host protein translation. We also discovered that IQGAP2, a Ras GTPase-activating-like protein, is a key substrate for SENP3-mediated de-SUMOylation. Downregulation of SENP3 in HBV infected cells facilitated IQGAP2 SUMOylation and degradation, which leads to suppression of HBV gene expression and restoration of global translation of host genes via modulation of AKT phosphorylation. Thus, The SENP3-IQGAP2 de-SUMOylation axis is a host defense mechanism of hepatocytes that restores host protein translation and suppresses HBV gene expression.

A Coding Sequence-Embedded Principle Governs Translational Reading Frame Fidelity

Upon initiation at a start codon, the ribosome must maintain the correct reading frame for hundreds of codons in order to produce functional proteins. While some sequence elements are able to trigger programmed ribosomal frameshifting (PRF), very little is known about how the ribosome normally prevents spontaneous frameshift errors that can have dire consequences if uncorrected. Using high resolution ribosome profiling data sets, we discovered that the translating ribosome uses the 3' end of 18S rRNA to scan the AUG-like codons after the decoding process. The postdecoding mRNA:rRNA interaction not only contributes to predominant translational pausing, but also provides a retrospective mechanism to safeguard the ribosome in the correct reading frame. Partially eliminating the AUG-like "sticky" codons in the reporter message leads to increased +1 frameshift errors. Remarkably, mutating the highly conserved CAU triplet of 18S rRNA globally changes the codon "stickiness". Further supporting the role of "sticky" sequences in reading frame maintenance, the codon composition of open reading frames is highly optimized across eukaryotic genomes. These results suggest an important layer of information embedded within the protein-coding sequences that instructs the ribosome to ensure reading frame fidelity during translation.

riboWaltz: Optimization of ribosome P-site positioning in ribosome profiling data

Ribosome profiling is a powerful technique used to study translation at the genome-wide level, generating unique information concerning ribosome positions along RNAs. Optimal localization of ribosomes requires the proper identification of the ribosome P-site in each ribosome protected fragment, a crucial step to determine the trinucleotide periodicity of translating ribosomes, and draw correct conclusions concerning where ribosomes are located. To determine the P-site within ribosome footprints at nucleotide resolution, the precise estimation of its offset with respect to the protected fragment is necessary. Here we present riboWaltz, an R package for calculation of optimal P-site offsets, diagnostic analysis and visual inspection of ribosome profiling data. Compared to existing tools, riboWaltz shows improved accuracies for P-site estimation and neat ribosome positioning in multiple case studies. riboWaltz was implemented in R and is available as an R package at https://github.com/LabTranslationalArchitectomics/RiboWaltz.

Effects of single amino acid deficiency on mRNA translation are markedly different for methionine versus leucine

Although amino acids are known regulators of translation, the unique contributions of specific amino acids are not well understood. We compared effects of culturing HEK293T cells in medium lacking either leucine, methionine, histidine, or arginine on eIF2 and 4EBP1 phosphorylation and measures of mRNA translation. Methionine starvation caused the most drastic decrease in translation as assessed by polysome formation, ribosome profiling, and a measure of protein synthesis (puromycin-labeled polypeptides) but had no significant effect on eIF2 phosphorylation, 4EBP1 hyperphosphorylation or 4EBP1 binding to eIF4E. Leucine starvation suppressed polysome formation and was the only tested condition that caused a significant decrease in 4EBP1 phosphorylation or increase in 4EBP1 binding to eIF4E, but effects of leucine starvation were not replicated by overexpressing nonphosphorylatable 4EBP1. This suggests the binding of 4EBP1 to eIF4E may not by itself explain the suppression of mRNA translation under conditions of leucine starvation. Ribosome profiling suggested that leucine deprivation may primarily inhibit ribosome loading, whereas methionine deprivation may primarily impair start site recognition. These data underscore our lack of a full understanding of how mRNA translation is regulated and point to a unique regulatory role of methionine status on translation initiation that is not dependent upon eIF2 phosphorylation.

Along the Central Dogma-Controlling Gene Expression with Small Molecules

The central dogma of molecular biology, that DNA is transcribed into RNA and RNA translated into protein, was coined in the early days of modern biology. Back in the 1950s and 1960s, bacterial genetics first opened the way toward understanding life as the genetically encoded interaction of macromolecules. As molecular biology progressed and our knowledge of gene control deepened, it became increasingly clear that expression relied on many more levels of regulation. In the process of dissecting mechanisms of gene expression, specific small-molecule inhibitors played an important role and became valuable tools of investigation. Small molecules offer significant advantages over genetic tools, as they allow inhibiting a process at any desired time point, whereas mutating or altering the gene of an important regulator would likely result in a dead organism. With the advent of modern sequencing technology, it has become possible to monitor global cellular effects of small-molecule treatment and thereby overcome the limitations of classical biochemistry, which usually looks at a biological system in isolation. This review focuses on several molecules, especially natural products, that have played an important role in dissecting gene expression and have opened up new fields of investigation as well as clinical venues for disease treatment.

Ribosome signatures aid bacterial translation initiation site identification

BACKGROUND

While methods for annotation of genes are increasingly reliable, the exact identification of translation initiation sites remains a challenging problem. Since the N-termini of proteins often contain regulatory and targeting information, developing a robust method for start site identification is crucial. Ribosome profiling reads show distinct patterns of read length distributions around translation initiation sites. These patterns are typically lost in standard ribosome profiling analysis pipelines, when reads from footprints are adjusted to determine the specific codon being translated.

RESULTS

Utilising these signatures in combination with nucleotide sequence information, we build a model capable of predicting translation initiation sites and demonstrate its high accuracy using N-terminal proteomics. Applying this to prokaryotic translatomes, we re-annotate translation initiation sites and provide evidence of N-terminal truncations and extensions of previously annotated coding sequences. These re-annotations are supported by the presence of structural and sequence-based features next to N-terminal peptide evidence. Finally, our model identifies 61 novel genes previously undiscovered in the Salmonella enterica genome.

CONCLUSIONS

Signatures within ribosome profiling read length distributions can be used in combination with nucleotide sequence information to provide accurate genome-wide identification of translation initiation sites.

Codon usage bias in 5' terminal coding sequences reveals distinct enrichment of gene functions

Codon bias at the 5' terminal of coding sequence (CDS) is known to be distinct from the rest of the CDS. A number of events occur in this short region to regulate early translation elongation and co-translational translocation. In the genes encoding secretory proteins, there is a special signal sequence which has a higher occurrence of rare codons. In this study, we analyzed codon bias of secretory genes in several eukaryotes. The results showed that secretory genes in the species except mammals had a higher occurrence of rare codons in the 5' terminal of CDS, and the bias was greater than the same region of non-secretory genes. GO analysis revealed that secretory genes containing rare codon clusters in different regions were responsible for various roles in gene functions. Moreover, codon bias in the region encoding the hydrophobic region of protein is similar in secretory and non-secretory genes, indicating that codon bias in secretory genes was partly influenced by amino acid bias. Rare codon clusters are found more frequently in specific regions, and continuous rare codons are not favoured probably because they will increase the probability of ribosome collision and drop-off. Based on ribosome profiling data, there is no significant difference in the average translation efficiencies between rare and optimal codons. Higher ribosomal density in the 5' terminal may result from ribosome pausing which could be involved in different translation events. These findings collectively provided rich information on codon bias in secretory genes, which may shed light on the co-effect of codon bias, mRNA structure and tRNA abundance in translational regulations.

Efficient in situ detection of mRNAs using the Chlorella virus DNA ligase for padlock probe ligation

Padlock probes are single-stranded DNA molecules that are circularized upon hybridization to their target sequence by a DNA ligase. In the following, the circulated padlock probes are amplified and detected with fluorescently labeled probes complementary to the amplification product. The hallmark of padlock probe assays is a high detection specificity gained by the ligation reaction. Concomitantly, the ligation reaction is the largest drawback for a quantitative in situ detection of mRNAs due to the low affinities of common DNA or RNA ligases to RNA-DNA duplex strands. Therefore, current protocols require that mRNAs be reverse transcribed to DNA before detection with padlock probes. Recently, it was found that the DNA ligase from Paramecium bursaria Chlorella virus 1 (PBCV-1) is able to efficiently ligate RNA-splinted DNA. Hence, we designed a padlock probe assay for direct in situ detection of mRNAs using the PBCV-1 DNA ligase. Experimental single-cell data were used to optimize and characterize the efficiency of mRNA detection with padlock probes. Our results demonstrate that the PBCV-1 DNA ligase overcomes the efficiency limitation of current protocols for direct in situ mRNA detection, making the PBCV-1 DNA ligase an attractive tool to simplify in situ ligation sequencing applications.

Predictive biophysical modeling and understanding of the dynamics of mRNA translation and its evolution

mRNA translation is the fundamental process of decoding the information encoded in mRNA molecules by the ribosome for the synthesis of proteins. The centrality of this process in various biomedical disciplines such as cell biology, evolution and biotechnology, encouraged the development of dozens of mathematical and computational models of translation in recent years. These models aimed at capturing various biophysical aspects of the process. The objective of this review is to survey these models, focusing on those based and/or validated on real large-scale genomic data. We consider aspects such as the complexity of the models, the biophysical aspects they regard and the predictions they may provide. Furthermore, we survey the central systems biology discoveries reported on their basis. This review demonstrates the fundamental advantages of employing computational biophysical translation models in general, and discusses the relative advantages of the different approaches and the challenges in the field.

Redefining the Translational Status of 80S Monosomes

Fully assembled ribosomes exist in two populations: polysomes and monosomes. While the former has been studied extensively, to what extent translation occurs on monosomes and its importance for overall translational output remain controversial. Here, we used ribosome profiling to examine the translational status of 80S monosomes in Saccharomyces cerevisiae. We found that the vast majority of 80S monosomes are elongating, not initiating. Further, most mRNAs exhibit some degree of monosome occupancy, with monosomes predominating on nonsense-mediated decay (NMD) targets, upstream open reading frames (uORFs), canonical ORFs shorter than ∼ 590 nt, and ORFs for which the total time required to complete elongation is substantially shorter than that required for initiation. Importantly, mRNAs encoding low-abundance regulatory proteins tend to be enriched in the monosome fraction. Our data highlight the importance of monosomes for the translation of highly regulated mRNAs.

Genome-wide assessment of differential translations with ribosome profiling data

The closely regulated process of mRNA translation is crucial for precise control of protein abundance and quality. Ribosome profiling, a combination of ribosome foot-printing and RNA deep sequencing, has been used in a large variety of studies to quantify genome-wide mRNA translation. Here, we developed Xtail, an analysis pipeline tailored for ribosome profiling data that comprehensively and accurately identifies differentially translated genes in pairwise comparisons. Applied on simulated and real datasets, Xtail exhibits high sensitivity with minimal false-positive rates, outperforming existing methods in the accuracy of quantifying differential translations. With published ribosome profiling datasets, Xtail does not only reveal differentially translated genes that make biological sense, but also uncovers new events of differential translation in human cancer cells on mTOR signalling perturbation and in human primary macrophages on interferon gamma (IFN-γ) treatment. This demonstrates the value of Xtail in providing novel insights into the molecular mechanisms that involve translational dysregulations.

Ribosome hibernation factor promotes Staphylococcal survival and differentially represses translation

In opportunistic Gram-positive Staphylococcus aureus, a small protein called hibernation-promoting factor (HPF_Sa) is sufficient to dimerize 2.5-MDa 70S ribosomes into a translationally inactive 100S complex. Although the 100S dimer is observed in only the stationary phase in Gram-negative gammaproteobacteria, it is ubiquitous throughout all growth phases in S. aureus. The biological significance of the 100S ribosome is poorly understood. Here, we reveal an important role of HPF_Sa in preserving ribosome integrity and poising cells for translational restart, a process that has significant clinical implications for relapsed staphylococcal infections. We found that the hpf null strain is severely impaired in long-term viability concomitant with a dramatic loss of intact ribosomes. Genome-wide ribosome profiling shows that eliminating HPF_Sa drastically increased ribosome occupancy at the 5′ end of specific mRNAs under nutrient-limited conditions, suggesting that HPF_Sa may suppress translation initiation. The protective function of HPF_Sa on ribosomes resides at the N-terminal conserved basic residues and the extended C-terminal segment, which are critical for dimerization and ribosome binding, respectively. These data provide significant insight into the functional consequences of 100S ribosome loss for protein synthesis and stress adaptation.

High-Resolution Analysis of Coronavirus Gene Expression by RNA Sequencing and Ribosome Profiling

Members of the family Coronaviridae have the largest genomes of all RNA viruses, typically in the region of 30 kilobases. Several coronaviruses, such as Severe acute respiratory syndrome-related coronavirus (SARS-CoV) and Middle East respiratory syndrome-related coronavirus (MERS-CoV), are of medical importance, with high mortality rates and, in the case of SARS-CoV, significant pandemic potential. Other coronaviruses, such as Porcine epidemic diarrhea virus and Avian coronavirus, are important livestock pathogens. Ribosome profiling is a technique which exploits the capacity of the translating ribosome to protect around 30 nucleotides of mRNA from ribonuclease digestion. Ribosome-protected mRNA fragments are purified, subjected to deep sequencing and mapped back to the transcriptome to give a global "snap-shot" of translation. Parallel RNA sequencing allows normalization by transcript abundance. Here we apply ribosome profiling to cells infected with Murine coronavirus, mouse hepatitis virus, strain A59 (MHV-A59), a model coronavirus in the same genus as SARS-CoV and MERS-CoV. The data obtained allowed us to study the kinetics of virus transcription and translation with exquisite precision. We studied the timecourse of positive and negative-sense genomic and subgenomic viral RNA production and the relative translation efficiencies of the different virus ORFs. Virus mRNAs were not found to be translated more efficiently than host mRNAs; rather, virus translation dominates host translation at later time points due to high levels of virus transcripts. Triplet phasing of the profiling data allowed precise determination of translated reading frames and revealed several translated short open reading frames upstream of, or embedded within, known virus protein-coding regions. Ribosome pause sites were identified in the virus replicase polyprotein pp1a ORF and investigated experimentally. Contrary to expectations, ribosomes were not found to pause at the ribosomal frameshift site. To our knowledge this is the first application of ribosome profiling to an RNA virus.

Ribosome profiling reveals the what, when, where and how of protein synthesis

Ribosome profiling, which involves the deep sequencing of ribosome-protected mRNA fragments, is a powerful tool for globally monitoring translation in vivo. The method has facilitated discovery of the regulation of gene expression underlying diverse and complex biological processes, of important aspects of the mechanism of protein synthesis, and even of new proteins, by providing a systematic approach for experimental annotation of coding regions. Here, we introduce the methodology of ribosome profiling and discuss examples in which this approach has been a key factor in guiding biological discovery, including its prominent role in identifying thousands of novel translated short open reading frames and alternative translation products.

Dual transcriptional-translational cascade permits cellular level tuneable expression control

The ability to induce gene expression in a small molecule dependent manner has led to many applications in target discovery, functional elucidation and bio-production. To date these applications have relied on a limited set of protein-based control mechanisms operating at the level of transcription initiation. The discovery, design and reengineering of riboswitches offer an alternative means by which to control gene expression. Here we report the development and characterization of a novel tunable recombinant expression system, termed RiboTite, which operates at both the transcriptional and translational level. Using standard inducible promoters and orthogonal riboswitches, a multi-layered modular genetic control circuit was developed to control the expression of both bacteriophage T7 RNA polymerase and recombinant gene(s) of interest. The system was benchmarked against a number of commonly used E. coli expression systems, and shows tight basal control, precise analogue tunability of gene expression at the cellular level, dose-dependent regulation of protein production rates over extended growth periods and enhanced cell viability. This novel system expands the number of E. coli expression systems for use in recombinant protein production and represents a major performance enhancement over and above the most widely used expression systems.

Protein synthesis as an integral quality control mechanism during ageing

Ageing is manifested as functional and structural deterioration that affects cell and tissue physiology. mRNA translation is a central cellular process, supplying cells with newly synthesized proteins. Accumulating evidence suggests that alterations in protein synthesis are not merely a corollary but rather a critical factor for the progression of ageing. Here, we survey protein synthesis regulatory mechanisms and focus on the pre-translational regulation of the process exerted by non-coding RNA species, RNA binding proteins and alterations of intrinsic RNA properties. In addition, we discuss the tight relationship between mRNA translation and two central pathways that modulate ageing, namely the insulin/IGF-1 and TOR signalling cascades. A thorough understanding of the complex interplay between protein synthesis regulation and ageing will provide critical insights into the pathogenesis of age-related disorders, associated with impaired proteostasis and protein quality control.

Widespread Co-translational RNA Decay Reveals Ribosome Dynamics

It is generally assumed that mRNAs undergoing translation are protected from decay. Here, we show that mRNAs are, in fact, co-translationally degraded. This is a widespread and conserved process affecting most genes, where 5'-3' transcript degradation follows the last translating ribosome, producing an in vivo ribosomal footprint. By sequencing the ends of 5' phosphorylated mRNA degradation intermediates, we obtain a genome-wide drug-free measurement of ribosome dynamics. We identify general translation termination pauses in both normal and stress conditions. In addition, we describe novel codon-specific ribosomal pausing sites in response to oxidative stress that are dependent on the RNase Rny1. Our approach is simple and straightforward and does not require the use of translational inhibitors or in vitro RNA footprinting that can alter ribosome protection patterns.

Transcriptome-wide measurement of ribosomal occupancy by ribosome profiling

Gene expression profiling provides a tool to analyze the internal states of cells or organisms, and their responses to perturbations. While global measurements of mRNA levels have thus been widely used for many years, it is only through the recent development of the ribosome profiling technique that an analogous examination of global mRNA translation programs has become possible. Ribosome profiling reveals which RNAs are being translated to what extent and where the translated open reading frames are located. In addition, different modes of translation regulation can be distinguished and characterized. Here, we present an optimized, step-by-step protocol for ribosome profiling. Although established in Caenorhabditis elegans, our protocol and optimization approaches should be equally usable for other model organisms or cell culture with little adaptation. Next to providing a protocol, we compare two different methods for isolation of single ribosomes and two different library preparations, and describe strategies to optimize the RNase digest and to reduce ribosomal RNA contamination in the libraries. Moreover, we discuss bioinformatic strategies to evaluate the quality of the data and explain how the data can be analyzed for different applications. In sum, this article seeks to facilitate the understanding, execution, and optimization of ribosome profiling experiments.