Distinct stages of the translation elongation cycle revealed by sequencing ribosome-protected mRNA fragments

An Aggregation Method to Identify the RNA Meta-Stable Secondary Structure and its Functionally Interpretable Structure Ensemble

RNA can provide vital cellular functions through its secondary or tertiary structure. Due to the low-throughput nature of experimental approaches, studies on RNA structures mainly resort to computational methods. However, current existing tools fail to consider RNA structure ensembles and do not provide ways to decipher functional hypotheses for the new predictions. In this research, a novel method was proposed to identify the functionally interpretable structure ensemble of a given RNA sequence and provide the meta-stable structure, or the most frequently observed functional RNA cellular conformation, based on the ensemble. In the prediction of meta-stable structures, the proposed method outperformed existing tools on a yeast test set. The inferred functional aspects were then manually checked and demonstrated a micro-averaging F1 value of 0.92. Further, a biological example of the yeast ASH1-E1 element was discussed to articulate that these functional aspects can also suggest testable hypotheses. Then the proposed method was verified to be well applicable to other species through a human test set. Finally, the proposed method was demonstrated to show resistance to sequence length-dependent performance deterioration.

High-Resolution Ribosome Profiling for Determining Ribosome Functional States During Translation Elongation

Translation elongation is a highly choreographed process that involves substantial conformational changes of the ribosome to accommodate aminoacyl-tRNAs and traverse along the mRNA template. To capture distinct functional states of the ribosome, a high-resolution ribosome profiling-based approach has been developed. By deep-sequencing differently sized ribosome-protected mRNA fragments, this approach captures not only ribosome positions but also their functional states in vivo across the Saccharomyces cerevisiae transcriptome with codon resolution. This chapter presents a condensed and step-by-step protocol for preserving ribosomes in their functional states using a cocktail of antibiotics that traps distinct steps of elongating ribosomes and for constructing a cDNA library derived from the ribosome-protected mRNA fragments for deep sequencing.

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.

Investigating molecular mechanisms of 2A-stimulated ribosomal pausing and frameshifting in Theilovirus

The 2A protein of Theiler's murine encephalomyelitis virus (TMEV) acts as a switch to stimulate programmed -1 ribosomal frameshifting (PRF) during infection. Here, we present the X-ray crystal structure of TMEV 2A and define how it recognises the stimulatory RNA element. We demonstrate a critical role for bases upstream of the originally predicted stem-loop, providing evidence for a pseudoknot-like conformation and suggesting that the recognition of this pseudoknot by beta-shell proteins is a conserved feature in cardioviruses. Through examination of PRF in TMEV-infected cells by ribosome profiling, we identify a series of ribosomal pauses around the site of PRF induced by the 2A-pseudoknot complex. Careful normalisation of ribosomal profiling data with a 2A knockout virus facilitated the identification, through disome analysis, of ribosome stacking at the TMEV frameshifting signal. These experiments provide unparalleled detail of the molecular mechanisms underpinning Theilovirus protein-stimulated frameshifting.

The highs and lows of ionizing radiation and its effects on protein synthesis

Ionizing radiation (IR) is a constant feature of our environment and one that can dramatically affect organismal health and development. Although the impacts of high-doses of IR on mammalian cells and systems have been broadly explored, there are still challenges in accurately quantifying biological responses to IR, especially in the low-dose range to which most individuals are exposed in their lifetime. The resulting uncertainty has led to the entrenchment of conservative radioprotection policies around the world. Thus, uncovering long-sought molecular mechanisms and tissue responses that are targeted by IR could lead to more informed policymaking and propose new therapeutic avenues for a variety of pathologies. One often overlooked target of IR is mRNA translation, a highly regulated cellular process that consumes more than 40% of the cell's energy. In response to environmental stimuli, regulation of mRNA translation allows for precise and rapid changes to the cellular proteome, and unsurprisingly high-dose of IR was shown to trigger a severe reprogramming of global protein synthesis allowing the cell to conserve energy by preventing the synthesis of unneeded proteins. Nonetheless, under these conditions, certain mRNAs encoding specific proteins are translationally favoured to produce the factors essential to repair the cell or send it down the path of no return through programmed cell death. Understanding the mechanisms controlling protein synthesis in response to varying doses of IR could provide novel insights into how this stress-mediated cellular adaptation is regulated and potentially uncover novel targets for radiosensitization or radioprotection. Here, we review the current literature on the effects of IR at both high- and low-dose on the mRNA translation machinery.

Modifications of Ribosome Profiling that Provide New Data on the Translation Regulation

Ribosome profiling (riboseq) has opened the possibilities for the genome-wide studies of translation in all living organisms. This method is based on deep sequencing of mRNA fragments protected by the ribosomes from hydrolysis by ribonucleases, the so-called ribosomal footprints (RFPs). Ribosomal profiling together with RNA sequencing allows not only to identify with a reasonable accuracy translated reading frames in the transcriptome, but also to track changes in gene expression in response to various stimuli. Notably, ribosomal profiling in its classical version has certain limitations. The size of the selected mRNA fragments is 25-35 nts, while RFPs of other sizes are usually omitted from analysis. Also, ribosomal profiling "averages" the data from all ribosomes and does not allow to study specific ribosomal complexes associated with particular translation factors. However, recently developed modifications of ribosomal profiling provide answers to a number of questions. Thus, it has become possible to analyze not only elongating, but also scanning and reinitiating ribosomes, to study events associated with the collision of ribosomes during mRNA translation, to discover new ways of cotranslational assembly of multisubunit protein complexes during translation, and to selectively isolate ribosomal complexes associated with certain protein factors. New data obtained using these modified approaches provide a better understanding of the mechanisms of translation regulation and the functional roles of translational apparatus components.

Charcot-Marie-Tooth mutation in glycyl-tRNA synthetase stalls ribosomes in a pre-accommodation state and activates integrated stress response

Toxic gain-of-function mutations in aminoacyl-tRNA synthetases cause a degeneration of peripheral motor and sensory axons, known as Charcot-Marie-Tooth (CMT) disease. While these mutations do not disrupt overall aminoacylation activity, they interfere with translation via an unknown mechanism. Here, we dissect the mechanism of function of CMT mutant glycyl-tRNA synthetase (CMT-GARS), using high-resolution ribosome profiling and reporter assays. We find that CMT-GARS mutants deplete the pool of glycyl-tRNAGly available for translation and inhibit the first stage of elongation, the accommodation of glycyl-tRNA into the ribosomal A-site, which causes ribosomes to pause at glycine codons. Moreover, ribosome pausing activates a secondary repression mechanism at the level of translation initiation, by inducing the phosphorylation of the alpha subunit of eIF2 and the integrated stress response. Thus, CMT-GARS mutant triggers translational repression via two interconnected mechanisms, affecting both elongation and initiation of translation.

Sex Differences in Primary HIV Infection: Revisiting the Role of TLR7-Driven Type 1 IFN Production by Plasmacytoid Dendritic Cells in Women

Plasmacytoid dendritic cells (pDCs) produce type I interferon (IFN-I) during HIV-1 infection in response to TLR7 stimulation. However, IFN-I-signaling has been shown to play opposite effects in HIV-1 and SIV infection. TLR7-driven type I interferon production in pDCs is higher in women than in men due to the cell-intrinsic actions of estrogen and X-chromosome complement. Indeed, TLR7 is encoded on the X-chromosome, and the TLR7 gene escapes the X-chromosome inactivation in immune cells of women which express significantly higher levels of TLR7 protein than male cells. Following HIV infection, women have a lower viremia during acute infection and exhibit stronger antiviral responses than men, which has been attributed to the increased capacity of female pDCs to produce IFN-α upon TLR7-stimulation. However, a deleterious functional impact of an excessive TLR7 response on acute viremia in women has been recently revealed by the analysis of the frequent rs179008 c.32A>T SNP of TLR7. This SNP was identified as a sex-specific protein abundance quantitative trait locus (pQTL) causing a difference in the TLR7 protein dosage and effector function in females only. T allele expression was associated with a lower TLR7 protein synthesis, blunted production of IFN-α by pDCs upon TLR7 stimulation, and an unexpectedly lower viral load during primary HIV-1 infection in women. In the present review, the author will revisit the role of TLR7-driven pDC innate function in the context of HIV-1 infection to discuss at what stage of primary HIV-1 infection the TLR7 rs179008 T allele is likely to be protective in women.

tRNA overexpression rescues peripheral neuropathy caused by mutations in tRNA synthetase

[Figure: see text].

Small-scale sequencing enables quality assessment of Ribo-Seq data: an example from Arabidopsis cell culture

BACKGROUND

Translation is a tightly regulated process, controlling the rate of protein synthesis in cells. Ribosome sequencing (Ribo-Seq) is a recently developed tool for studying actively translated mRNA and can thus directly address translational regulation. Ribo-Seq libraries need to be sequenced to a great depth due to high contamination by rRNA and other contaminating nucleic acid fragments. Deep sequencing is expensive, and it generates large volumes of data, making data analysis complicated and time consuming.

METHODS AND RESULTS

Here we developed a platform for Ribo-Seq library construction and data analysis to enable rapid quality assessment of Ribo-Seq libraries with the help of a small-scale sequencer. Our data show that several qualitative features of a Ribo-Seq library, such as read length distribution, P-site distribution, reading frame and triplet periodicity, can be effectively evaluated using only the data generated by a benchtop sequencer with a very limited number of reads.

CONCLUSION

Our pipeline enables rapid evaluation of Ribo-Seq libraries, opening up possibilities for optimization of Ribo-Seq library construction from difficult samples, and leading to better decision making prior to more costly deep sequencing.

Humans and other commonly used model organisms are resistant to cycloheximide-mediated biases in ribosome profiling experiments

Ribosome profiling measures genome-wide translation dynamics at sub-codon resolution. Cycloheximide (CHX), a widely used translation inhibitor to arrest ribosomes in these experiments, has been shown to induce biases in yeast, questioning its use. However, whether such biases are present in datasets of other organisms including humans is unknown. Here we compare different CHX-treatment conditions in human cells and yeast in parallel experiments using an optimized protocol. We find that human ribosomes are not susceptible to conformational restrictions by CHX, nor does it distort gene-level measurements of ribosome occupancy, measured decoding speed or the translational ramp. Furthermore, CHX-induced codon-specific biases on ribosome occupancy are not detectable in human cells or other model organisms. This shows that reported biases of CHX are species-specific and that CHX does not affect the outcome of ribosome profiling experiments in most settings. Our findings provide a solid framework to conduct and analyze ribosome profiling experiments.

Ribosome profiling has become the gold standard to analyze mRNA translation dynamics, and the translation inhibitor cycloheximide (CHX) is often used in its application. Here the authors systematically demonstrate that CHX does not bias the outcome of ribosome profiling experiments in most organisms.

Combinatorial analysis of translation dynamics reveals eIF2 dependence of translation initiation at near-cognate codons

Although ribosome-profiling and translation initiation sequencing (TI-seq) analyses have identified many noncanonical initiation codons, the precise detection of translation initiation sites (TISs) remains a challenge, mainly because of experimental artifacts of such analyses. Here, we describe a new method, TISCA (TIS detection by translation Complex Analysis), for the accurate identification of TISs. TISCA proved to be more reliable for TIS detection compared with existing tools, and it identified a substantial number of near-cognate codons in Kozak-like sequence contexts. Analysis of proteomics data revealed the presence of methionine at the NH₂-terminus of most proteins derived from near-cognate initiation codons. Although eukaryotic initiation factor 2 (eIF2), eIF2A and eIF2D have previously been shown to contribute to translation initiation at near-cognate codons, we found that most noncanonical initiation events are most probably dependent on eIF2, consistent with the initial amino acid being methionine. Comprehensive identification of TISs by TISCA should facilitate characterization of the mechanism of noncanonical initiation.

Persistence of Ambigrammatic Narnaviruses Requires Translation of the Reverse Open Reading Frame

Narnaviruses are RNA viruses detected in diverse fungi, plants, protists, arthropods, and nematodes. Though initially described as simple single-gene nonsegmented viruses encoding RNA-dependent RNA polymerase (RdRp), a subset of narnaviruses referred to as "ambigrammatic" harbor a unique genomic configuration consisting of overlapping open reading frames (ORFs) encoded on opposite strands. Phylogenetic analysis supports selection to maintain this unusual genome organization, but functional investigations are lacking. Here, we establish the mosquito-infecting Culex narnavirus 1 (CxNV1) as a model to investigate the functional role of overlapping ORFs in narnavirus replication. In CxNV1, a reverse ORF without homology to known proteins covers nearly the entire 3.2-kb segment encoding the RdRp. Additionally, two opposing and nearly completely overlapping novel ORFs are found on the second putative CxNV1 segment, the 0.8-kb "Robin" RNA. We developed a system to launch CxNV1 in a naive mosquito cell line and then showed that functional RdRp is required for persistence of both segments, and an intact reverse ORF is required on the RdRp segment for persistence. Mass spectrometry of persistently CxNV1-infected cells provided evidence for translation of this reverse ORF. Finally, ribosome profiling yielded a striking pattern of footprints for all four CxNV1 RNA strands that was distinct from actively translating ribosomes on host mRNA or coinfecting RNA viruses. Taken together, these data raise the possibility that the process of translation itself is important for persistence of ambigrammatic narnaviruses, potentially by protecting viral RNA with ribosomes, thus suggesting a heretofore undescribed viral tactic for replication and transmission. IMPORTANCE Fundamental to our understanding of RNA viruses is a description of which strand(s) of RNA are transmitted as the viral genome relative to which encode the viral proteins. Ambigrammatic narnaviruses break the mold. These viruses, found broadly in fungi, plants, and insects, have the unique feature of two overlapping genes encoded on opposite strands, comprising nearly the full length of the viral genome. Such extensive overlap is not seen in other RNA viruses and comes at the cost of reduced evolutionary flexibility in the sequence. The present study is motivated by investigating the benefits which balance that cost. We show for the first time a functional requirement for the ambigrammatic genome configuration in Culex narnavirus 1, which suggests a model for how translation of both strands might benefit this virus. Our work highlights a new blueprint for viral persistence, distinct from strategies defined by canonical definitions of the coding strand.

Genome-wide Survey of Ribosome Collision

Ribosome movement is not always smooth and is rather often impeded. For ribosome pauses, fundamental issues remain to be addressed, including where ribosomes pause on mRNAs, what kind of RNA/amino acid sequence causes this pause, and the physiological significance of this attenuation of protein synthesis. Here, we survey the positions of ribosome collisions caused by ribosome pauses in humans and zebrafish using modified ribosome profiling. Collided ribosomes, i.e., disomes, emerge at various sites: Pro-Pro/Gly/Asp motifs; Arg-X-Lys motifs; stop codons; and 3' untranslated regions. The electrostatic interaction between the charged nascent chain and the ribosome exit tunnel determines the eIF5A-mediated disome rescue at the Pro-Pro sites. In particular, XBP1u, a precursor of endoplasmic reticulum (ER)-stress-responsive transcription factor, shows striking queues of collided ribosomes and thus acts as a degradation substrate by ribosome-associated quality control. Our results provide insight into the causes and consequences of ribosome pause by dissecting collided ribosomes.

Deep conservation of ribosome stall sites across RNA processing genes

The rate of translation can vary depending on the mRNA template. During the elongation phase the ribosome can transiently pause or permanently stall. A pause can provide the nascent protein with the time to fold or be transported, while stalling can serve as quality control and trigger degradation of aberrant mRNA and peptide. Ribosome profiling has allowed for the genome-wide detection of such pauses and stalls, but due to library-specific biases, these predictions are often unreliable. Here, we take advantage of the deep conservation of protein synthesis machinery, hypothesizing that similar conservation could exist for functionally important locations of ribosome slowdown, here collectively called stall sites. We analyze multiple ribosome profiling datasets from phylogenetically diverse eukaryotes: yeast, fruit fly, zebrafish, mouse and human to identify conserved stall sites. We find thousands of stall sites across multiple species, with the enrichment of proline, glycine and negatively charged amino acids around conserved stalling. Many of the sites are found in RNA processing genes, suggesting that stalling might have a conserved role in RNA metabolism. In summary, our results provide a rich resource for the study of conserved stalling and indicate possible roles of stalling in gene regulation.

RiboDoc: A Docker-based package for ribosome profiling analysis

Ribosome profiling (RiboSeq) has emerged as a powerful technique for studying the genome-wide regulation of translation in various cells. Several steps in the biological protocol have been improved, but the bioinformatics part of RiboSeq suffers from a lack of standardization, preventing the straightforward and complete reproduction of published results. Too many published studies provide insufficient detail about the bioinformatics pipeline used. The broad range of questions that can be asked with RiboSeq makes it difficult to use a single bioinformatics tool. Indeed, many scripts have been published for addressing diverse questions. Here (https://github.com/equipeGST/RiboDoc), we propose a unique tool (for use with multiple operating systems, OS) to standardize the general steps that must be performed systematically in RiboSeq analysis, together with the statistical analysis and quality control of the sample. The data generated can then be exploited with more specific tools. We hope that this tool will help to standardize bioinformatics analyses pipelines in the field of translation.

Synonymous mutations that regulate translation speed might play a non-negligible role in liver cancer development

Background

Synonymous mutations do not change the protein sequences. Automatically, they have been regarded as neutral events and are ignored in the mutation-based cancer studies. However, synonymous mutations will change the codon optimality, resulting in altered translational velocity.

Methods

We fully utilized the transcriptome and translatome of liver cancer and normal tissue from ten patients. We profiled the mutation spectrum and examined the effect of synonymous mutations on translational velocity.

Results

Synonymous mutations that increase the codon optimality significantly enhanced the translational velocity, and were enriched in oncogenes. Meanwhile, synonymous mutations decreasing codon optimality slowed down translation, and were enriched in tumor suppressor genes. These synonymous mutations significantly contributed to the translational changes in tumor samples compared to normal samples.

Conclusions

Synonymous mutations might play a role in liver cancer development by altering codon optimality and translational velocity. Synonymous mutations should no longer be ignored in the genome-wide studies.

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.

Tools for Assessing Translation in Cryptococcus neoformans

Cryptococcus neoformans is a ubiquitous environmental fungus capable of establishing an infection in a human host. Rapid changes in environments and exposure to the host immune system results in a significant amount of cellular stress, which is effectively combated at the level of translatome reprogramming. Repression of translation following stress allows for the specific reallocation of limited resources. Understanding the mechanisms involved in regulating translation in C. neoformans during host infection is critical in the development of new antifungal drugs. In this review, we discuss the main tools available for assessing changes in translation state and translational output during cellular stress.

Influence of nascent polypeptide positive charges on translation dynamics

Protein segments with a high concentration of positively charged amino acid residues are often used in reporter constructs designed to activate ribosomal mRNA/protein decay pathways, such as those involving nonstop mRNA decay (NSD), no-go mRNA decay (NGD) and the ribosome quality control (RQC) complex. It has been proposed that the electrostatic interaction of the positively charged nascent peptide with the negatively charged ribosomal exit tunnel leads to translation arrest. When stalled long enough, the translation process is terminated with the degradation of the transcript and an incomplete protein. Although early experiments made a strong argument for this mechanism, other features associated with positively charged reporters, such as codon bias and mRNA and protein structure, have emerged as potent inducers of ribosome stalling. We carefully reviewed the published data on the protein and mRNA expression of artificial constructs with diverse compositions as assessed in different organisms. We concluded that, although polybasic sequences generally lead to lower translation efficiency, it appears that an aggravating factor, such as a nonoptimal codon composition, is necessary to cause translation termination events.

Positive epistasis between disease-causing missense mutations and silent polymorphism with effect on mRNA translation velocity

Epistasis refers to the dependence of a mutation on other mutation(s) and the genetic context in general. In the context of human disorders, epistasis complicates the spectrum of disease symptoms and has been proposed as a major contributor to variations in disease outcome. The nonadditive relationship between mutations and the lack of complete understanding of the underlying physiological effects limit our ability to predict phenotypic outcome. Here, we report positive epistasis between intragenic mutations in the cystic fibrosis transmembrane conductance regulator (CFTR)-the gene responsible for cystic fibrosis (CF) pathology. We identified a synonymous single-nucleotide polymorphism (sSNP) that is invariant for the CFTR amino acid sequence but inverts translation speed at the affected codon. This sSNP in cis exhibits positive epistatic effects on some CF disease-causing missense mutations. Individually, both mutations alter CFTR structure and function, yet when combined, they lead to enhanced protein expression and activity. The most robust effect was observed when the sSNP was present in combination with missense mutations that, along with the primary amino acid change, also alter the speed of translation at the affected codon. Functional studies revealed that synergistic alteration in ribosomal velocity is the underlying mechanism; alteration of translation speed likely increases the time window for establishing crucial domain-domain interactions that are otherwise perturbed by each individual mutation.

Translational landscape and protein biogenesis demands of the early secretory pathway in Komagataella phaffii

Background

Eukaryotes use distinct networks of biogenesis factors to synthesize, fold, monitor, traffic, and secrete proteins. During heterologous expression, saturation of any of these networks may bottleneck titer and yield. To understand the flux through various routes into the early secretory pathway, we quantified the global and membrane-associated translatomes of Komagataella phaffii.

Results

By coupling Ribo-seq with long-read mRNA sequencing, we generated a new annotation of protein-encoding genes. By using Ribo-seq with subcellular fractionation, we quantified demands on co- and posttranslational translocation pathways. During exponential growth in rich media, protein components of the cell-wall represent the greatest number of nascent chains entering the ER. Transcripts encoding the transmembrane protein PMA1 sequester more ribosomes at the ER membrane than any others. Comparison to Saccharomyces cerevisiae reveals conservation in the resources allocated by gene ontology, but variation in the diversity of gene products entering the secretory pathway.

Conclusion

A subset of host proteins, particularly cell-wall components, impose the greatest biosynthetic demands in the early secretory pathway. These proteins are potential targets in strain engineering aimed at alleviating bottlenecks during heterologous protein production.

A stress-induced tyrosine-tRNA depletion response mediates codon-based translational repression and growth suppression

Eukaryotic transfer RNAs can become selectively fragmented upon various stresses, generating tRNA-derived small RNA fragments. Such fragmentation has been reported to impact a small fraction of the tRNA pool and thus presumed to not directly impact translation. We report that oxidative stress can rapidly generate tyrosine-tRNAGUA fragments in human cells-causing significant depletion of the precursor tRNA. Tyrosine-tRNAGUA depletion impaired translation of growth and metabolic genes enriched in cognate tyrosine codons. Depletion of tyrosine tRNAGUA or its translationally regulated targets USP3 and SCD repressed proliferation-revealing a dedicated tRNA-regulated growth-suppressive pathway for oxidative stress response. Tyrosine fragments are generated in a DIS3L2 exoribonuclease-dependent manner and inhibit hnRNPA1-mediated transcript destabilization. Moreover, tyrosine fragmentation is conserved in C. elegans. Thus, tRNA fragmentation can coordinately generate trans-acting small RNAs and functionally deplete a tRNA. Our findings reveal the existence of an underlying adaptive codon-based regulatory response inherent to the genetic code.

Insights Into Translatomics in the Nervous System

Most neurological disorders are caused by abnormal gene translation. Generally, dysregulation of elements involved in the translational process disrupts homeostasis in neurons and neuroglia. Better understanding of how the gene translation process occurs requires detailed analysis of transcriptomic and proteomic profile data. However, a lack of strictly direct correlations between mRNA and protein levels limits translational investigation by combining transcriptomic and proteomic profiling. The much better correlation between proteins and translated mRNAs than total mRNAs in abundance and insufficiently sensitive proteomics approach promote the requirement of advances in translatomics technology. Translatomics which capture and sequence the mRNAs associated with ribosomes has been effective in identifying translational changes by genetics or projections, ribosome stalling, local translation, and transcript isoforms in the nervous system. Here, we place emphasis on the main three translatomics methods currently used to profile mRNAs attached to ribosome-nascent chain complex (RNC-mRNA). Their prominent applications in neurological diseases including glioma, neuropathic pain, depression, fragile X syndrome (FXS), neurodegenerative disorders are outlined. The content reviewed here expands our understanding on the contributions of aberrant translation to neurological disease development.

Disome-seq reveals widespread ribosome collisions that promote cotranslational protein folding

BACKGROUND

The folding of proteins is challenging in the highly crowded and sticky environment of a cell. Regulation of translation elongation may play a crucial role in ensuring the correct folding of proteins. Much of our knowledge regarding translation elongation comes from the sequencing of mRNA fragments protected by single ribosomes by ribo-seq. However, larger protected mRNA fragments have been observed, suggesting the existence of an alternative and previously hidden layer of regulation.

RESULTS

In this study, we performed disome-seq to sequence mRNA fragments protected by two stacked ribosomes, a product of translational pauses during which the 5'-elongating ribosome collides with the 3'-paused one. We detected widespread ribosome collisions that are related to slow ribosome release when stop codons are at the A-site, slow peptide bond formation from proline, glycine, asparagine, and cysteine when they are at the P-site, and slow leaving of polylysine from the exit tunnel of ribosomes. The structure of disomes obtained by cryo-electron microscopy suggests a different conformation from the substrate of the ribosome-associated protein quality control pathway. Collisions occurred more frequently in the gap regions between α-helices, where a translational pause can prevent the folding interference from the downstream peptides. Paused or collided ribosomes are associated with specific chaperones, which can aid in the cotranslational folding of the nascent peptides.

CONCLUSIONS

Therefore, cells use regulated ribosome collisions to ensure protein homeostasis.

Performing Ribosome Profiling to Assess Translation in Vegetative and Meiotic Yeast Cells

Ribosome profiling, first developed in 2009, is the gold standard for quantifying and qualifying changes to translation genome-wide (Ingolia et al., Science, 2009). Though first designed and optimized in vegetative budding yeast, it has since been modified and specialized for use in diverse cellular states in yeast, as well as in bacteria, plants, human cells, and many other organisms (Ingolia et al. Science, 2009, reviewed in (Ingolia et al., Cold Spring Harb Perspect Biol, 2019; Brar and Weissman, Nat Rev Mol Cell Biol, 2015)). Here we report the current ribosome profiling protocol used in our lab to study genome-wide changes to translation in budding yeast undergoing the developmental process of meiosis (Brar et al., Science, 2012; Cheng et al., Cell, 2018). We describe this protocol in detail, including the following steps: collection and flash freezing samples, cell lysis and extract preparation, sucrose gradient centrifugation and monosome collection, RNA extraction, library preparation, and library quality control. Almost every step presented here should be directly applicable to performing ribosome profiling in other eukaryotic cell types or cell states.

Nuclease-mediated depletion biases in ribosome footprint profiling libraries

Ribosome footprint profiling is a high-throughput sequencing-based technique that provides detailed and global views of translation in living cells. An essential part of this technology is removal of unwanted, normally very abundant, ribosomal RNA sequences that dominate libraries and increase sequencing costs. The most effective commercial solution (Ribo-Zero) has been discontinued as a standalone product and a number of new, experimentally distinct commercial applications have emerged on the market. Here we evaluated several commercially available alternatives designed for RNA-seq of human samples and find them generally unsuitable for ribosome footprint profiling. We instead recommend the use of custom-designed biotinylated oligos, which were widely used in early ribosome profiling studies. Importantly, we warn that depletion solutions based on targeted nuclease cleavage significantly perturb the high-resolution information that can be derived from the data, and thus do not recommend their use for any applications that require precise determination of the ends of RNA fragments.

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.

Control of translation by eukaryotic mRNA transcript leaders-Insights from high-throughput assays and computational modeling

Eukaryotic gene expression is tightly regulated during translation of mRNA to protein. Mis-regulation of translation can lead to aberrant proteins which accumulate in cancers and cause neurodegenerative diseases. Foundational studies on model genes established fundamental roles for mRNA 5' transcript leader (TL) sequences in controlling ribosome recruitment, scanning, and initiation. TL cis-regulatory elements and their corresponding trans-acting factors control cap-dependent initiation under unstressed conditions. Under stress, cap-dependent initiation is suppressed, and specific mRNA structures and sequences promote translation of stress-responsive transcripts to remodel the proteome. In this review, we summarize current knowledge of TL functions in translation initiation. We focus on insights from high-throughput analyses of ribosome occupancy, mRNA structure, RNA Binding Protein occupancy, and massively parallel reporter assays. These data-driven approaches, coupled with computational analyses and modeling, have paved the way for a comprehensive understanding of TL functions. Finally, we will discuss areas of future research on the roles of mRNA sequences and structures in translation. This article is categorized under: Translation > Translation Mechanisms RNA Evolution and Genomics > Computational Analyses of RNA RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

Selective Translation Complex Profiling Reveals Staged Initiation and Co-translational Assembly of Initiation Factor Complexes

Translational control targeting the initiation phase is central to the regulation of gene expression. Understanding all of its aspects requires substantial technological advancements. Here we modified yeast translation complex profile sequencing (TCP-seq), related to ribosome profiling, and adapted it for mammalian cells. Human TCP-seq, capable of capturing footprints of 40S subunits (40Ss) in addition to 80S ribosomes (80Ss), revealed that mammalian and yeast 40Ss distribute similarly across 5'TRs, indicating considerable evolutionary conservation. We further developed yeast and human selective TCP-seq (Sel-TCP-seq), enabling selection of 40Ss and 80Ss associated with immuno-targeted factors. Sel-TCP-seq demonstrated that eIF2 and eIF3 travel along 5' UTRs with scanning 40Ss to successively dissociate upon AUG recognition; notably, a proportion of eIF3 lingers on during the initial elongation cycles. Highlighting Sel-TCP-seq versatility, we also identified four initiating 48S conformational intermediates, provided novel insights into ATF4 and GCN4 mRNA translational control, and demonstrated co-translational assembly of initiation factor complexes.

Ribosomes: An Exciting Avenue in Stem Cell Research

Stem cell research has focused on genomic studies. However, recent evidence has indicated the involvement of epigenetic regulation in determining the fate of stem cells. Ribosomes play a crucial role in epigenetic regulation, and thus, we focused on the role of ribosomes in stem cells. Majority of living organisms possess ribosomes that are involved in the translation of mRNA into proteins and promote cellular proliferation and differentiation. Ribosomes are stable molecular machines that play a role with changes in the levels of RNA during translation. Recent research suggests that specific ribosomes actively regulate gene expression in multiple cell types, such as stem cells. Stem cells have the potential for self-renewal and differentiation into multiple lineages and, thus, require high efficiency of translation. Ribosomes induce cellular transdifferentiation and reprogramming, and disrupted ribosome synthesis affects translation efficiency, thereby hindering stem cell function leading to cell death and differentiation. Stem cell function is regulated by ribosome-mediated control of stem cell-specific gene expression. In this review, we have presented a detailed discourse on the characteristics of ribosomes in stem cells. Understanding ribosome biology in stem cells will provide insights into the regulation of stem cell function and cellular reprogramming.

Recommendations for bacterial ribosome profiling experiments based on bioinformatic evaluation of published data

Ribosome profiling (RIBO-Seq) has improved our understanding of bacterial translation, including finding many unannotated genes. However, protocols for RIBO-Seq and corresponding data analysis are not yet standardized. Here, we analyzed 48 RIBO-Seq samples from nine studies of Escherichia coli K12 grown in lysogeny broth medium and particularly focused on the size-selection step. We show that for conventional expression analysis, a size range between 22 and 30 nucleotides is sufficient to obtain protein-coding fragments, which has the advantage of removing many unwanted rRNA and tRNA reads. More specific analyses may require longer reads and a corresponding improvement in rRNA/tRNA depletion. There is no consensus about the appropriate sequencing depth for RIBO-Seq experiments in prokaryotes, and studies vary significantly in total read number. Our analysis suggests that 20 million reads that are not mapping to rRNA/tRNA are required for global detection of translated annotated genes. We also highlight the influence of drug-induced ribosome stalling, which causes bias at translation start sites. The resulting accumulation of reads at the start site may be especially useful for detecting weakly expressed genes. As different methods suit different questions, it may not be possible to produce a "one-size-fits-all" ribosome profiling data set. Therefore, experiments should be carefully designed in light of the scientific questions of interest. We propose some basic characteristics that should be reported with any new RIBO-Seq data sets. Careful attention to the factors discussed should improve prokaryotic gene detection and the comparability of ribosome profiling data sets.

A molecular-level perspective on the frequency, distribution, and consequences of messenger RNA modifications

Cells use chemical modifications to alter the sterics, charge, and conformations of large biomolecules, modulating their biogenesis, function, and stability. Until recently post-transcriptional RNA modifications were thought to be largely limited to nonprotein coding RNA species. However, this dogma has rapidly transformed with the discovery of a host of modifications in protein coding messenger RNAs (mRNAs). Recent advancements in genome-wide sequencing technologies have enabled the identification of mRNA modifications as a potential new frontier in gene regulation-leading to the development of the epitranscriptome field. As a result, there has been a flurry of multiple groundbreaking discoveries, including new modifications, nucleoside modifying enzymes ("writers" and "erasers"), and RNA binding proteins that recognize chemical modifications ("readers"). These discoveries opened the door to understanding how post-transcriptional mRNA modifications can modulate the mRNA lifecycle, and established a link between the epitranscriptome and human health and disease. Despite a rapidly growing recognition of their importance, fundamental questions regarding the identity, prevalence, and functional consequences of mRNA modifications remain to be answered. Here, we highlight quantitative studies that characterize mRNA modification abundance, frequency, and interactions with cellular machinery. As the field progresses, we see a need for the further integration of quantitative and reductionist approaches to complement transcriptome wide studies in order to establish a molecular-level framework for understanding the consequences of mRNA chemical modifications on biological processes. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Processing > RNA Editing and Modification.

Disome and Trisome Profiling Reveal Genome-wide Targets of Ribosome Quality Control

The ribosome-associated protein quality control (RQC) system that resolves stalled translation events is activated when ribosomes collide and form disome, trisome, or higher-order complexes. However, it is unclear whether this system distinguishes collision complexes formed on defective mRNAs from those with functional roles on endogenous transcripts. Here, we performed disome and trisome footprint profiling in yeast and found collisions were enriched on diverse sequence motifs known to slow translation. When 60S recycling was inhibited, disomes accumulated at stop codons and could move into the 3' UTR to reinitiate translation. The ubiquitin ligase and RQC factor Hel2/ZNF598 generally recognized collisions but did not induce degradation of endogenous transcripts. However, loss of Hel2 triggered the integrated stress response, via phosphorylation of eIF2α, thus linking these pathways. Our results suggest that Hel2 has a role in sensing ribosome collisions on endogenous mRNAs, and such events may be important for cellular homeostasis.

Slowing ribosome velocity restores folding and function of mutant CFTR

Cystic fibrosis (CF) is caused by mutations in the CF transmembrane conductance regulator (CFTR), with approximately 90% of patients harboring at least one copy of the disease-associated variant F508del. We utilized a yeast phenomic system to identify genetic modifiers of F508del-CFTR biogenesis, from which ribosomal protein L12 (RPL12/uL11) emerged as a molecular target. In the present study, we investigated mechanism(s) by which suppression of RPL12 rescues F508del protein synthesis and activity. Using ribosome profiling, we found that rates of translation initiation and elongation were markedly slowed by RPL12 silencing. However, proteolytic stability and patch-clamp assays revealed RPL12 depletion significantly increased F508del-CFTR steady-state expression, interdomain assembly, and baseline open-channel probability. We next evaluated whether Rpl12-corrected F508del-CFTR could be further enhanced with concomitant pharmacologic repair (e.g., using clinically approved modulators lumacaftor and tezacaftor) and demonstrated additivity of these treatments. Rpl12 knockdown also partially restored maturation of specific CFTR variants in addition to F508del, and WT Cftr biogenesis was enhanced in the pancreas, colon, and ileum of Rpl12 haplosufficient mice. Modulation of ribosome velocity therefore represents a robust method for understanding both CF pathogenesis and therapeutic response.

Detection and Degradation of Stalled Nascent Chains via Ribosome-Associated Quality Control

Stalled protein synthesis produces defective nascent chains that can harm cells. In response, cells degrade these nascent chains via a process called ribosome-associated quality control (RQC). Here, we review the irregularities in the translation process that cause ribosomes to stall as well as how cells use RQC to detect stalled ribosomes, ubiquitylate their tethered nascent chains, and deliver the ubiquitylated nascent chains to the proteasome. We additionally summarize how cells respond to RQC failure.

TLR7 dosage polymorphism shapes interferogenesis and HIV-1 acute viremia in women

Type I IFN (IFN-I) production by plasmacytoid DCs (pDCs) occurs during acute HIV-1 infection in response to TLR7 stimulation, but the role of pDC-derived IFN-I in controlling or promoting HIV-1 infection is ambiguous. We report here a sex-biased interferogenic phenotype for a frequent single-nucleotide polymorphism of human TLR7, rs179008, displaying an impact on key parameters of acute HIV-1 infection. We show allele rs179008 T to determine lower TLR7 protein abundance in cells from women, specifically - likely by diminishing TLR7 mRNA translation efficiency through codon usage. The hypomorphic TLR7 phenotype is mirrored by decreased TLR7-driven IFN-I production by female pDCs. Among women from the French ANRS PRIMO cohort of acute HIV-1 patients, carriage of allele rs179008 T associated with lower viremia, cell-associated HIV-1 DNA, and CXCL10 (IP-10) plasma concentrations. RNA viral load was decreased by 0.85 log10 (95% CI, -1.51 to -0.18) among T/T homozygotes, who also exhibited a lower frequency of acute symptoms. TLR7 emerges as an important control locus for acute HIV-1 viremia, and the clinical phenotype for allele rs179008 T, carried by 30%-50% of European women, supports a beneficial effect of toning down TLR7-driven IFN-I production by pDCs during acute HIV-1 infection.

Ribosome profiling in archaea reveals leaderless translation, novel translational initiation sites, and ribosome pausing at single codon resolution

High-throughput methods, such as ribosome profiling, have revealed the complexity of translation regulation in Bacteria and Eukarya with large-scale effects on cellular functions. In contrast, the translational landscape in Archaea remains mostly unexplored. Here, we developed ribosome profiling in a model archaeon, Haloferax volcanii, elucidating, for the first time, the translational landscape of a representative of the third domain of life. We determined the ribosome footprint of H. volcanii to be comparable in size to that of the Eukarya. We linked footprint lengths to initiating and elongating states of the ribosome on leadered transcripts, operons, and on leaderless transcripts, the latter representing 70% of H. volcanii transcriptome. We manipulated ribosome activity with translation inhibitors to reveal ribosome pausing at specific codons. Lastly, we found that the drug harringtonine arrested ribosomes at initiation sites in this archaeon. This drug treatment allowed us to confirm known translation initiation sites and also reveal putative novel initiation sites in intergenic regions and within genes. Ribosome profiling revealed an uncharacterized complexity of translation in this archaeon with bacteria-like, eukarya-like, and potentially novel translation mechanisms. These mechanisms are likely to be functionally essential and to contribute to an expanded proteome with regulatory roles in gene expression.

Deficiency of mitoribosomal S10 protein affects translation and splicing in Arabidopsis mitochondria

The ribosome is not only a protein-making machine, but also a regulatory element in protein synthesis. This view is supported by our earlier data showing that Arabidopsis mitoribosomes altered due to the silencing of the nuclear RPS10 gene encoding mitochondrial ribosomal protein S10 differentially translate mitochondrial transcripts compared with the wild-type. Here, we used ribosome profiling to determine the contribution of transcriptional and translational control in the regulation of protein synthesis in rps10 mitochondria compared with the wild-type ones. Oxidative phosphorylation system proteins are preferentially synthesized in wild-type mitochondria but this feature is lost in the mutant. The rps10 mitoribosomes show slightly reduced translation efficiency of most respiration-related proteins and at the same time markedly more efficiently synthesize ribosomal proteins and MatR and TatC proteins. The mitoribosomes deficient in S10 protein protect shorter transcript fragments which exhibit a weaker 3-nt periodicity compared with the wild-type. The decrease in the triplet periodicity is particularly drastic for genes containing introns. Notably, splicing is considerably less effective in the mutant, indicating an unexpected link between the deficiency of S10 and mitochondrial splicing. Thus, a shortage of the mitoribosomal S10 protein has wide-ranging consequences on mitochondrial gene expression.

Ribosome assembly defects subvert initiation Factor3 mediated scrutiny of bona fide start signal

In bacteria, the assembly factors tightly orchestrate the maturation of ribosomes whose competency for protein synthesis is validated by translation machinery at various stages of translation cycle. However, what transpires to the quality control measures when the ribosomes are produced with assembly defects remains enigmatic. In Escherichia coli, we show that 30S ribosomes that harbour assembly defects due to the lack of assembly factors such as RbfA and KsgA display suboptimal initiation codon recognition and bypass the critical codon–anticodon proofreading steps during translation initiation. These premature ribosomes on entering the translation cycle compromise the fidelity of decoding that gives rise to errors during initiation and elongation. We show that the assembly defects compromise the binding of initiation factor 3 (IF3), which in turn appears to license the rapid transition of 30S (pre) initiation complex to 70S initiation complex by tempering the validation of codon–anticodon interaction during translation initiation. This suggests that the premature ribosomes harbouring the assembly defects subvert the IF3 mediated proofreading of cognate initiation codon to enter the translation cycle.

Computational methods for ribosome profiling data analysis

Since the introduction of the ribosome profiling technique in 2009 its popularity has greatly increased. It is widely used for the comprehensive assessment of gene expression and for studying the mechanisms of regulation at the translational level. As the number of ribosome profiling datasets being produced continues to grow, so too does the need for reliable software that can provide answers to the biological questions it can address. This review describes the computational methods and tools that have been developed to analyze ribosome profiling data at the different stages of the process. It starts with initial routine processing of raw data and follows with more specific tasks such as the identification of translated open reading frames, differential gene expression analysis, or evaluation of local or global codon decoding rates. The review pinpoints challenges associated with each step and explains the ways in which they are currently addressed. In addition it provides a comprehensive, albeit incomplete, list of publicly available software applicable to each step, which may be a beneficial starting point to those unexposed to ribosome profiling analysis. The outline of current challenges in ribosome profiling data analysis may inspire computational biologists to search for novel, potentially superior, solutions that will improve and expand the bioinformatician's toolbox for ribosome profiling data analysis. This article is characterized under: Translation > Ribosome Structure/Function RNA Evolution and Genomics > Computational Analyses of RNA Translation > Translation Mechanisms Translation > Translation Regulation.

RiboFlow, RiboR and RiboPy: an ecosystem for analyzing ribosome profiling data at read length resolution

SUMMARY

Ribosome occupancy measurements enable protein abundance estimation and infer mechanisms of translation. Recent studies have revealed that sequence read lengths in ribosome profiling data are highly variable and carry critical information. Consequently, data analyses require the computation and storage of multiple metrics for a wide range of ribosome footprint lengths. We developed a software ecosystem including a new efficient binary file format named 'ribo'. Ribo files store all essential data grouped by ribosome footprint lengths. Users can assemble ribo files using our RiboFlow pipeline that processes raw ribosomal profiling sequencing data. RiboFlow is highly portable and customizable across a large number of computational environments with built-in capabilities for parallelization. We also developed interfaces for writing and reading ribo files in the R (RiboR) and Python (RiboPy) environments. Using RiboR and RiboPy, users can efficiently access ribosome profiling quality control metrics, generate essential plots and carry out analyses. Altogether, these components create a software ecosystem for researchers to study translation through ribosome profiling.

AVAILABILITY AND IMPLEMENTATION

For a quickstart, please see https://ribosomeprofiling.github.io. Source code, installation instructions and links to documentation are available on GitHub: https://github.com/ribosomeprofiling.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Robust landscapes of ribosome dwell times and aminoacyl-tRNAs in response to nutrient stress in liver

Translation depends on messenger RNA (mRNA)-specific initiation, elongation, and termination rates. While translation elongation is well studied in bacteria and yeast, less is known in higher eukaryotes. Here we combined ribosome and transfer RNA (tRNA) profiling to investigate the relations between translation elongation rates, (aminoacyl-) tRNA levels, and codon usage in mammals. We modeled codon-specific ribosome dwell times from ribosome profiling, considering codon pair interactions between ribosome sites. In mouse liver, the model revealed site- and codon-specific dwell times that differed from those in yeast, as well as pairs of adjacent codons in the P and A site that markedly slow down or speed up elongation. While translation efficiencies vary across diurnal time and feeding regimen, codon dwell times were highly stable and conserved in human. Measured tRNA levels correlated with codon usage and several tRNAs showed reduced aminoacylation, which was conserved in fasted mice. Finally, we uncovered that the longest codon dwell times could be explained by aminoacylation levels or high codon usage relative to tRNA abundance.

Pervasive functional translation of noncanonical human open reading frames

Ribosome profiling has revealed pervasive but largely uncharacterized translation outside of canonical coding sequences (CDSs). In this work, we exploit a systematic CRISPR-based screening strategy to identify hundreds of noncanonical CDSs that are essential for cellular growth and whose disruption elicits specific, robust transcriptomic and phenotypic changes in human cells. Functional characterization of the encoded microproteins reveals distinct cellular localizations, specific protein binding partners, and hundreds of microproteins that are presented by the human leukocyte antigen system. We find multiple microproteins encoded in upstream open reading frames, which form stable complexes with the main, canonical protein encoded on the same messenger RNA, thereby revealing the use of functional bicistronic operons in mammals. Together, our results point to a family of functional human microproteins that play critical and diverse cellular roles.

RiboVIEW: a computational framework for visualization, quality control and statistical analysis of ribosome profiling data

Recently, newly developed ribosome profiling methods based on high-throughput sequencing of ribosome-protected mRNA footprints allow to study genome-wide translational changes in detail. However, computational analysis of the sequencing data still represents a bottleneck for many laboratories. Further, specific pipelines for quality control and statistical analysis of ribosome profiling data, providing high levels of both accuracy and confidence, are currently lacking. In this study, we describe automated bioinformatic and statistical diagnoses to perform robust quality control of ribosome profiling data (RiboQC), to efficiently visualize ribosome positions and to estimate ribosome speed (RiboMine) in an unbiased way. We present an R pipeline to setup and undertake the analyses that offers the user an HTML page to scan own data regarding the following aspects: periodicity, ligation and digestion of footprints; reproducibility and batch effects of replicates; drug-related artifacts; unbiased codon enrichment including variability between mRNAs, for A, P and E sites; mining of some causal or confounding factors. We expect our pipeline to allow an optimal use of the wealth of information provided by ribosome profiling experiments.

Ribosome profiling of HEK293T cells overexpressing codon optimized coagulation factor IX

Ribosome profiling provides the opportunity to evaluate translation kinetics at codon level resolution. Here, we describe ribosome profiling data, generated from two HEK293T cell lines. The ribosome profiling data are composed of Ribo-seq (mRNA sequencing data from ribosome protected fragments) and RNA-seq data (total RNA sequencing). The two HEK293T cell lines each express a version of the F9 gene, both of which are translated into identical proteins in terms of their amino acid sequences. However, these F9 genes vary drastically in their codon usage and predicted mRNA structure. We also provide the pipeline that we used to analyze the data. Further analyzing this dataset holds great potential as it can be used i) to unveil insights into the composition and regulation of the transcriptome, ii) for comparison with other ribosome profiling datasets, iii) to measure the rate of protein synthesis across the proteome and identify differences in elongation rates, iv) to discover previously unidentified translation of peptides, v) to explore the effects of codon usage or codon context in translational kinetics and vi) to investigate cotranslational folding. Importantly, a unique feature of this dataset, compared to other available ribosome profiling data, is the presence of the F9 gene in two very distinct coding sequences.

Elucidation of Codon Usage Signatures across the Domains of Life

Because of the degeneracy of the genetic code, multiple codons are translated into the same amino acid. Despite being “synonymous,” these codons are not equally used. Selective pressures are thought to drive the choice among synonymous codons within a genome, while GC content, which is typically attributed to mutational drift, is the major determinant of variation across species. Here, we find that in addition to GC content, interspecies codon usage signatures can also be detected. More specifically, we show that a single amino acid, arginine, is the major contributor to codon usage bias differences across domains of life. We then exploit this finding and show that domain-specific codon bias signatures can be used to classify a given sequence into its corresponding domain of life with high accuracy. We then wondered whether the inclusion of codon usage codon autocorrelation patterns, which reflects the nonrandom distribution of codon occurrences throughout a transcript, might improve the classification performance of our algorithm. However, we find that autocorrelation patterns are not domain-specific, and surprisingly, are unrelated to tRNA reusage, in contrast to previous reports. Instead, our results suggest that codon autocorrelation patterns are a by-product of codon optimality throughout a sequence, where highly expressed genes display autocorrelated “optimal” codons, whereas lowly expressed genes display autocorrelated “nonoptimal” codons.

A role for the Saccharomyces cerevisiae ABCF protein New1 in translation termination/recycling

Translation is controlled by numerous accessory proteins and translation factors. In the yeast Saccharomyces cerevisiae, translation elongation requires an essential elongation factor, the ABCF ATPase eEF3. A closely related protein, New1, is encoded by a non-essential gene with cold sensitivity and ribosome assembly defect knock-out phenotypes. Since the exact molecular function of New1 is unknown, it is unclear if the ribosome assembly defect is direct, i.e. New1 is a bona fide assembly factor, or indirect, for instance due to a defect in protein synthesis. To investigate this, we employed yeast genetics, cryo-electron microscopy (cryo-EM) and ribosome profiling (Ribo-Seq) to interrogate the molecular function of New1. Overexpression of New1 rescues the inviability of a yeast strain lacking the otherwise strictly essential translation factor eEF3. The structure of the ATPase-deficient (EQ₂) New1 mutant locked on the 80S ribosome reveals that New1 binds analogously to the ribosome as eEF3. Finally, Ribo-Seq analysis revealed that loss of New1 leads to ribosome queuing upstream of 3′-terminal lysine and arginine codons, including those genes encoding proteins of the cytoplasmic translational machinery. Our results suggest that New1 is a translation factor that fine-tunes the efficiency of translation termination or ribosome recycling.

Effects of codon optimization on coagulation factor IX translation and structure: Implications for protein and gene therapies

Synonymous codons occur with different frequencies in different organisms, a phenomenon termed codon usage bias. Codon optimization, a common term for a variety of approaches used widely by the biopharmaceutical industry, involves synonymous substitutions to increase protein expression. It had long been presumed that synonymous variants, which, by definition, do not alter the primary amino acid sequence, have no effect on protein structure and function. However, a critical mass of reports suggests that synonymous codon variations may impact protein conformation. To investigate the impact of synonymous codons usage on protein expression and function, we designed an optimized coagulation factor IX (FIX) variant and used multiple methods to compare its properties to the wild-type FIX upon expression in HEK293T cells. We found that the two variants differ in their conformation, even when controlling for the difference in expression levels. Using ribosome profiling, we identified robust changes in the translational kinetics of the two variants and were able to identify a region in the gene that may have a role in altering the conformation of the protein. Our data have direct implications for codon optimization strategies, for production of recombinant proteins and gene therapies.

Computational design and interpretation of single-RNA translation experiments

Advances in fluorescence microscopy have introduced new assays to quantify live-cell translation dynamics at single-RNA resolution. We introduce a detailed, yet efficient sequence-based stochastic model that generates realistic synthetic data for several such assays, including Fluorescence Correlation Spectroscopy (FCS), ribosome Run-Off Assays (ROA) after Harringtonine application, and Fluorescence Recovery After Photobleaching (FRAP). We simulate these experiments under multiple imaging conditions and for thousands of human genes, and we evaluate through simulations which experiments are most likely to provide accurate estimates of elongation kinetics. Finding that FCS analyses are optimal for both short and long length genes, we integrate our model with experimental FCS data to capture the nascent protein statistics and temporal dynamics for three human genes: KDM5B, β-actin, and H2B. Finally, we introduce a new open-source software package, RNA Sequence to NAscent Protein Simulator (rSNAPsim), to easily simulate the single-molecule translation dynamics of any gene sequence for any of these assays and for different assumptions regarding synonymous codon usage, tRNA level modifications, or ribosome pauses. rSNAPsim is implemented in Python and is available at: https://github.com/MunskyGroup/rSNAPsim.git.