Comparative survey of the relative impact of mRNA features on local ribosome profiling read density

Riboseq-flow: A streamlined, reliable pipeline for ribosome profiling data analysis and quality control

Ribosome profiling is a powerful technique to study translation at a transcriptome-wide level. However, ensuring good data quality is paramount for accurate interpretation, as is ensuring that the analyses are reproducible. We introduce a new Nextflow DSL2 pipeline, riboseq-flow, designed for processing and comprehensive quality control of ribosome profiling experiments. Riboseq-flow is user-friendly, versatile and upholds high standards in reproducibility, scalability, portability, version control and continuous integration. It enables users to efficiently analyse multiple samples in parallel and helps them evaluate the quality and utility of their data based on the detailed metrics and visualisations that are automatically generated. Riboseq-flow is available at https://github.com/iraiosub/riboseq-flow.

Riboformer: a deep learning framework for predicting context-dependent translation dynamics

Translation elongation is essential for maintaining cellular proteostasis, and alterations in the translational landscape are associated with a range of diseases. Ribosome profiling allows detailed measurements of translation at the genome scale. However, it remains unclear how to disentangle biological variations from technical artifacts in these data and identify sequence determinants of translation dysregulation. Here we present Riboformer, a deep learning-based framework for modeling context-dependent changes in translation dynamics. Riboformer leverages the transformer architecture to accurately predict ribosome densities at codon resolution. When trained on an unbiased dataset, Riboformer corrects experimental artifacts in previously unseen datasets, which reveals subtle differences in synonymous codon translation and uncovers a bottleneck in translation elongation. Further, we show that Riboformer can be combined with in silico mutagenesis to identify sequence motifs that contribute to ribosome stalling across various biological contexts, including aging and viral infection. Our tool offers a context-aware and interpretable approach for standardizing ribosome profiling datasets and elucidating the regulatory basis of translation kinetics.

Comprehensive quantitative modeling of translation efficiency in a genome-reduced bacterium

Translation efficiency has been mainly studied by ribosome profiling, which only provides an incomplete picture of translation kinetics. Here, we integrated the absolute quantifications of tRNAs, mRNAs, RNA half-lives, proteins, and protein half-lives with ribosome densities and derived the initiation and elongation rates for 475 genes (67% of all genes), 73 with high precision, in the bacterium Mycoplasma pneumoniae (Mpn). We found that, although the initiation rate varied over 160-fold among genes, most of the known factors had little impact on translation efficiency. Local codon elongation rates could not be fully explained by the adaptation to tRNA abundances, which varied over 100-fold among tRNA isoacceptors. We provide a comprehensive quantitative view of translation efficiency, which suggests the existence of unidentified mechanisms of translational regulation in Mpn.

Neoantigen identification: Technological advances and challenges

Neoantigens have emerged as promising targets for cutting-edge immunotherapies, such as cancer vaccines and adoptive cell therapy. These neoantigens are unique to tumors and arise exclusively from somatic mutations or non-genomic aberrations in tumor proteins. They encompass a wide range of alterations, including genomic mutations, post-transcriptomic variants, and viral oncoproteins. With the advancements in technology, the identification of immunogenic neoantigens has seen rapid progress, raising new opportunities for enhancing their clinical significance. Prediction of neoantigens necessitates the acquisition of high-quality samples and sequencing data, followed by mutation calling. Subsequently, the pipeline involves integrating various tools that can predict the expression, processing, binding, and recognition potential of neoantigens. However, the continuous improvement of computational tools is constrained by the availability of datasets which contain validated immunogenic neoantigens. This review article aims to provide a comprehensive summary of the current knowledge as well as limitations in neoantigen prediction and validation. Additionally, it delves into the origin and biological role of neoantigens, offering a deeper understanding of their significance in the field of cancer immunotherapy. This article thus seeks to contribute to the ongoing efforts to harness neoantigens as powerful weapons in the fight against cancer.

Riboformer: A Deep Learning Framework for Predicting Context-Dependent Translation Dynamics

Translation elongation is essential for maintaining cellular proteostasis, and alterations in the translational landscape are associated with a range of diseases. Ribosome profiling allows detailed measurement of translation at genome scale. However, it remains unclear how to disentangle biological variations from technical artifacts and identify sequence determinant of translation dysregulation. Here we present Riboformer, a deep learning-based framework for modeling context-dependent changes in translation dynamics. Riboformer leverages the transformer architecture to accurately predict ribosome densities at codon resolution. It corrects experimental artifacts in previously unseen datasets, reveals subtle differences in synonymous codon translation and uncovers a bottleneck in protein synthesis. Further, we show that Riboformer can be combined with in silico mutagenesis analysis to identify sequence motifs that contribute to ribosome stalling across various biological contexts, including aging and viral infection. Our tool offers a context-aware and interpretable approach for standardizing ribosome profiling datasets and elucidating the regulatory basis of translation kinetics.

RPL3L-containing ribosomes determine translation elongation dynamics required for cardiac function

Although several ribosomal protein paralogs are expressed in a tissue-specific manner, how these proteins affect translation and why they are required only in certain tissues have remained unclear. Here we show that RPL3L, a paralog of RPL3 specifically expressed in heart and skeletal muscle, influences translation elongation dynamics. Deficiency of RPL3L-containing ribosomes in RPL3L knockout male mice resulted in impaired cardiac contractility. Ribosome occupancy at mRNA codons was found to be altered in the RPL3L-deficient heart, and the changes were negatively correlated with those observed in myoblasts overexpressing RPL3L. RPL3L-containing ribosomes were less prone to collisions compared with RPL3-containing canonical ribosomes. Although the loss of RPL3L-containing ribosomes altered translation elongation dynamics for the entire transcriptome, its effects were most pronounced for transcripts related to cardiac muscle contraction and dilated cardiomyopathy, with the abundance of the encoded proteins being correspondingly decreased. Our results provide further insight into the mechanisms and physiological relevance of tissue-specific translational regulation.

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.

Evaluating data integrity in ribosome footprinting datasets through modelled polysome profiles

The assessment of transcriptome-wide ribosome binding to mRNAs is useful for studying the dynamic regulation of protein synthesis. Two methods frequently applied in eukaryotic cells that operate at different levels of resolution are polysome profiling, which reveals the distribution of ribosome loads across the transcriptome, and ribosome footprinting (also termed ribosome profiling or Ribo-Seq), which when combined with appropriate data on mRNA expression can reveal ribosome densities on individual transcripts. In this study we develop methods for relating the information content of these two methods to one another, by reconstructing theoretical polysome profiles from ribosome footprinting data. Our results validate both approaches as experimental tools. Although we show that both methods can yield highly consistent data, some published ribosome footprinting datasets give rise to reconstructed polysome profiles with non-physiological features. We trace these aberrant features to inconsistencies in RNA and Ribo-Seq data when compared to datasets yielding physiological polysome profiles, thereby demonstrating that modelled polysomes are useful for assessing global dataset properties such as its quality in a simple, visual approach. Aside from using polysome profile reconstructions on published datasets, we propose that this also provides a useful tool for validating new ribosome footprinting datasets in early stages of analyses.

Flow-Seq Method: Features and Application in Bacterial Translation Studies

The Flow-seq method is based on using reporter construct libraries, where a certain element regulating the gene expression of fluorescent reporter proteins is represented in many thousands of variants. Reporter construct libraries are introduced into cells, sorted according to their fluorescence level, and then subjected to next-generation sequencing. Therefore, it turns out to be possible to identify patterns that determine the expression efficiency, based on tens and hundreds of thousands of reporter constructs in one experiment. This method has become common in evaluating the efficiency of protein synthesis simultaneously by multiple mRNA variants. However, its potential is not confined to this area. In the presented review, a comparative analysis of the Flow-seq method and other alternative approaches used for translation efficiency evaluation of mRNA was carried out; the features of its application and the results obtained by Flow-seq were also considered.

Tracing Translational Footprint by Ribo-Seq: Principle, Workflow, and Applications to Understand the Mechanism of Human Diseases

RNA-seq has been widely used as a high-throughput method to characterize transcript dynamic changes in a broad context, such as development and diseases. However, whether RNA-seq-estimated transcriptional dynamics can be translated into protein level changes is largely unknown. Ribo-seq (Ribosome profiling) is an emerging technology that allows for the investigation of the translational footprint via profiling ribosome-bounded mRNA fragments. Ribo-seq coupled with RNA-seq will allow us to understand the transcriptional and translational control of the fundamental biological process and human diseases. This review focuses on discussing the principle, workflow, and applications of Ribo-seq to study human diseases.

Distinct elongation stalls during translation are linked with distinct pathways for mRNA degradation

Key protein adapters couple translation to mRNA decay on specific classes of problematic mRNAs in eukaryotes. Slow decoding on non-optimal codons leads to codon-optimality-mediated decay (COMD) and prolonged arrest at stall sites leads to no-go decay (NGD). The identities of the decay factors underlying these processes and the mechanisms by which they respond to translational distress remain open areas of investigation. We use carefully designed reporter mRNAs to perform genetic screens and functional assays in Saccharomyces cerevisiae. We characterize the roles of Hel2, Syh1, and Smy2 in coordinating translational repression and mRNA decay on NGD reporter mRNAs, finding that Syh1 and, to a lesser extent its paralog Smy2, act in a distinct pathway from Hel2. This Syh1/Smy2-mediated pathway acts as a redundant, compensatory pathway to elicit NGD when Hel2-dependent NGD is impaired. Importantly, we observe that these NGD factors are not involved in the degradation of mRNAs enriched in non-optimal codons. Further, we establish that a key factor previously implicated in COMD, Not5, contributes modestly to the degradation of an NGD-targeted mRNA. Finally, we use ribosome profiling to reveal distinct ribosomal states associated with each reporter mRNA that readily rationalize the contributions of NGD and COMD factors to degradation of these reporters. Taken together, these results provide new insight into the role of Syh1 and Smy2 in NGD and into the ribosomal states that correlate with the activation of distinct pathways targeting mRNAs for degradation in yeast.

Stress-induced perturbations in intracellular amino acids reprogram mRNA translation in osmoadaptation independently of the ISR

The integrated stress response (ISR) plays a pivotal role in adaptation of translation machinery to cellular stress. Here, we demonstrate an ISR-independent osmoadaptation mechanism involving reprogramming of translation via coordinated but independent actions of mTOR and plasma membrane amino acid transporter SNAT2. This biphasic response entails reduced global protein synthesis and mTOR signaling followed by translation of SNAT2. Induction of SNAT2 leads to accumulation of amino acids and reactivation of mTOR and global protein synthesis, paralleled by partial reversal of the early-phase, stress-induced translatome. We propose SNAT2 functions as a molecular switch between inhibition of protein synthesis and establishment of an osmoadaptive translation program involving the formation of cytoplasmic condensates of SNAT2-regulated RNA-binding proteins DDX3X and FUS. In summary, we define key roles of SNAT2 in osmotolerance.

The landscape of translational stall sites in bacteria revealed by monosome and disome profiling

Ribosome pauses are associated with various cotranslational events and determine the fate of mRNAs and proteins. Thus, the identification of precise pause sites across the transcriptome is desirable; however, the landscape of ribosome pauses in bacteria remains ambiguous. Here, we harness monosome and disome (or collided ribosome) profiling strategies to survey ribosome pause sites in Escherichia coli Compared to eukaryotes, ribosome collisions in bacteria showed remarkable differences: a low frequency of disomes at stop codons, collisions occurring immediately after 70S assembly on start codons, and shorter queues of ribosomes trailing upstream. The pause sites corresponded with the biochemical validation by integrated nascent chain profiling (iNP) to detect polypeptidyl-tRNA, an elongation intermediate. Moreover, the subset of those sites showed puromycin resistance, presenting slow peptidyl transfer. Among the identified sites, the ribosome pause at Asn586 of ycbZ was validated by biochemical reporter assay, tRNA sequencing (tRNA-seq), and cryo-electron microscopy (cryo-EM) experiments. Our results provide a useful resource for ribosome stalling sites in bacteria.

Low-input RNase footprinting for simultaneous quantification of cytosolic and mitochondrial translation

We describe a low-input RNase footprinting approach for the rapid quantification of ribosome-protected fragments with as few as 1000 cultured cells. The assay uses a simplified procedure to selectively capture ribosome footprints based on optimized RNase digestion. It simultaneously maps cytosolic and mitochondrial translation with single-nucleotide resolution. We applied it to reveal selective functions of the elongation factor TUFM in mitochondrial translation, as well as synchronized repression of cytosolic translation after TUFM perturbation. We show the assay is applicable to small amounts of primary tissue samples with low protein synthesis rates, including snap-frozen tissues and immune cells from an individual's blood draw. We showed its feasibility to characterize the personalized immuno-translatome. Our analyses revealed that thousands of genes show lower translation efficiency in monocytes compared with lymphocytes, and identified thousands of translated noncanonical open reading frames (ORFs). Altogether, our RNase footprinting approach opens an avenue to assay transcriptome-wide translation using low-input samples from a wide range of physiological conditions.

Polysome Fractionation for Transcriptome-Wide Studies of mRNA Translation

Protein synthesis and degradation determine the relationship between mRNA and corresponding protein amounts. This relationship can change in a dynamic and selective fashion when translational efficiencies of transcript subsets are altered downstream of, for example, translation factors and/or RNA binding proteins. Notably, even transcription factors such as estrogen receptor alpha (ERα) can modulate mRNA translation in a transcript-selective manner. Yet, despite ample evidence suggesting a key role for mRNA translation in shaping the proteome in health and disease, it remains largely unexplored. Here, we present a guide for the extraction of mRNA engaged in translation using polysome fractionation with linear and optimized sucrose gradients. The isolated polysome-associated RNA is then quantified, in parallel with total mRNA from the same conditions, using methods such as RNA sequencing; and the resulting data set is analyzed to derive transcriptome-wide insights into how mRNA translation is modulated. The methods we describe are applicable to cultured cells, small numbers of FACS-isolated primary cells, and small tissue samples from biobanks or animal studies. Accordingly, this approach can be applied to study in detail how ERα and other factors control gene expression by selectively modulating mRNA translation both in vitro and in vivo.

A critical period of translational control during brain development at codon resolution

Translation modulates the timing and amplification of gene expression after transcription. Brain development requires uniquely complex gene expression patterns, but large-scale measurements of translation directly in the prenatal brain are lacking. We measure the reactants, synthesis and products of mRNA translation spanning mouse neocortex neurogenesis, and discover a transient window of dynamic regulation at mid-gestation. Timed translation upregulation of chromatin-binding proteins like Satb2, which is essential for neuronal subtype differentiation, restricts protein expression in neuronal lineages despite broad transcriptional priming in progenitors. In contrast, translation downregulation of ribosomal proteins sharply decreases ribosome biogenesis, coinciding with a major shift in protein synthesis dynamics at mid-gestation. Changing activity of eIF4EBP1, a direct inhibitor of ribosome biogenesis, is concurrent with ribosome downregulation and affects neurogenesis of the Satb2 lineage. Thus, the molecular logic of brain development includes the refinement of transcriptional programs by translation. Modeling of the developmental neocortex translatome is provided as an open-source searchable resource at https://shiny.mdc-berlin.de/cortexomics .

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.

Ribo-DT: An automated pipeline for inferring codon dwell times from ribosome profiling data

Protein synthesis is an energy consuming process characterised as a pivotal and highly regulated step in gene expression. The net protein output is dictated by a combination of translation initiation, elongation and termination rates that have remained difficult to measure. Recently, the development of ribosome profiling has enabled the inference of translation parameters through modelling, as this method informs on the ribosome position along the mRNA. Here, we present an automated, reproducible and portable computational pipeline to infer relative single-codon and codon-pair dwell times as well as gene flux from raw ribosome profiling sequencing data. As a case study, we applied our workflow to a publicly available yeast ribosome profiling dataset consisting of 57 independent gene knockouts related to RNA and tRNA modifications. We uncovered the effects of those modifications on translation elongation and codon selection during decoding. In particular, knocking out mod5 and trm7 increases codon-specific dwell times which indicates their potential tRNA targets, and highlights effects of nucleotide modifications on ribosome decoding rate.

Exploring Evidence of Non-coding RNA Translation With Trips-Viz and GWIPS-Viz Browsers

Detection of translation in so-called non-coding RNA provides an opportunity for identification of novel bioactive peptides and microproteins. The main methods used for these purposes are ribosome profiling and mass spectrometry. A number of publicly available datasets already exist for a substantial number of different cell types grown under various conditions, and public data mining is an attractive strategy for identification of translation in non-coding RNAs. Since the analysis of publicly available data requires intensive data processing, several data resources have been created recently for exploring processed publicly available data, such as OpenProt, GWIPS-viz, and Trips-Viz. In this work we provide a detailed demonstration of how to use the latter two tools for exploring experimental evidence for translation of RNAs hitherto classified as non-coding. For this purpose, we use a set of transcripts with substantially different patterns of ribosome footprint distributions. We discuss how certain features of these patterns can be used as evidence for or against genuine translation. During our analysis we concluded that the MTLN mRNA, previously misannotated as lncRNA LINC00116, likely encodes only a short proteoform expressed from shorter RNA transcript variants.

Estimation of peptide elongation times from ribosome profiling spectra

Ribosome profiling spectra bear rich information on translation control and dynamics. Yet, due to technical biases in library generation, extracting quantitative measures of discrete translation events has remained elusive. Using maximum likelihood statistics and data set from Escherichia coli we develop a robust method for neutralizing technical biases (e.g. base specific RNase preferences in ribosome-protected mRNA fragments (RPF) generation), which allows for correct estimation of translation times at single codon resolution. Furthermore, we validated the method with available datasets from E. coli treated with antibiotic to inhibit isoleucyl-tRNA synthetase, and two datasets from Saccharomyces cerevisiae treated with two RNases with distinct cleavage signatures. We demonstrate that our approach accounts for RNase cleavage preferences and provides bias-corrected translation times estimates. Our approach provides a solution to the long-standing problem of extracting reliable information about peptide elongation times from highly noisy and technically biased ribosome profiling spectra.

Quantitative Modeling of Protein Synthesis Using Ribosome Profiling Data

Quantitative prediction on protein synthesis requires accurate translation initiation and codon translation rates. Ribosome profiling data, which provide steady-state distribution of relative ribosome occupancies along a transcript, can be used to extract these rate parameters. Various methods have been developed in the past few years to measure translation-initiation and codon translation rates from ribosome profiling data. In the review, we provide a detailed analysis of the key methods employed to extract the translation rate parameters from ribosome profiling data. We further discuss how these approaches were used to decipher the role of various structural and sequence-based features of mRNA molecules in the regulation of gene expression. The utilization of these accurate rate parameters in computational modeling of protein synthesis may provide new insights into the kinetic control of the process of gene expression.

CONCUR: quick and robust calculation of codon usage from ribosome profiling data

SUMMARY

CONCUR is a standalone tool for codon usage analysis in ribosome profiling experiments. CONCUR uses the aligned reads in BAM format to estimate codon counts at the ribosome E-, P- and A-sites and at flanking positions.

AVAILABILITY AND IMPLEMENTATION

CONCUR is written in Perl and is freely available at https://github.com/susbo/concur.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Full-length ribosome density prediction by a multi-input and multi-output model

Translation elongation is regulated by a series of complicated mechanisms in both prokaryotes and eukaryotes. Although recent advance in ribosome profiling techniques has enabled one to capture the genome-wide ribosome footprints along transcripts at codon resolution, the regulatory codes of elongation dynamics are still not fully understood. Most of the existing computational approaches for modeling translation elongation from ribosome profiling data mainly focus on local contextual patterns, while ignoring the continuity of the elongation process and relations between ribosome densities of remote codons. Modeling the translation elongation process in full-length coding sequence (CDS) level has not been studied to the best of our knowledge. In this paper, we developed a deep learning based approach with a multi-input and multi-output framework, named RiboMIMO, for modeling the ribosome density distributions of full-length mRNA CDS regions. Through considering the underlying correlations in translation efficiency among neighboring and remote codons and extracting hidden features from the input full-length coding sequence, RiboMIMO can greatly outperform the state-of-the-art baseline approaches and accurately predict the ribosome density distributions along the whole mRNA CDS regions. In addition, RiboMIMO explores the contributions of individual input codons to the predictions of output ribosome densities, which thus can help reveal important biological factors influencing the translation elongation process. The analyses, based on our interpretable metric named codon impact score, not only identified several patterns consistent with the previously-published literatures, but also for the first time (to the best of our knowledge) revealed that the codons located at a long distance from the ribosomal A site may also have an association on the translation elongation rate. This finding of long-range impact on translation elongation velocity may shed new light on the regulatory mechanisms of protein synthesis. Overall, these results indicated that RiboMIMO can provide a useful tool for studying the regulation of translation elongation in the range of full-length CDS.

Translation elongation is a process in which amino acids are linked into proteins by ribosomes in cells. Translation elongation rates along the mRNAs are not constant, and are regulated by a series of mechanisms, such as codon rarity and mRNA stability. In this study, we modeled the translation elongation process at a full-length coding sequence level and developed a deep learning based approach to predict the translation elongation rates from mRNA sequences, through extracting the regulatory codes of elongation rates from the contextual sequences. The analyses, based on our interpretable metric named codon impact score, for the first time (to the best of our knowledge), revealed that in addition to the neighboring codons of the ribosomal A sites, the remote codons may also have an important impact on the translation elongation rates. This new finding may stimulate additional experiments and shed light on the regulatory mechanisms of protein synthesis.

Inferring translational heterogeneity from Saccharomyces cerevisiae ribosome profiling

Translation of mRNAs into proteins by the ribosome is the most important step of protein biosynthesis. Accordingly, translation is tightly controlled and heavily regulated to maintain cellular homeostasis. Ribosome profiling (Ribo-seq) has revolutionized the study of translation by revealing many of its underlying mechanisms. However, equally many aspects of translation remain mysterious, in part also due to persisting challenges in the interpretation of data obtained from Ribo-seq experiments. Here, we show that some of the variability observed in Ribo-seq data has biological origins and reflects programmed heterogeneity of translation. Through a comparative analysis of Ribo-seq data from Saccharomyces cerevisiae, we systematically identify short 3-codon sequences that are differentially translated (DT) across mRNAs, that is, identical sequences that are translated sometimes fast and sometimes slowly beyond what can be attributed to variability between experiments. Remarkably, the thus identified DT sequences link to mechanisms known to regulate translation elongation and are enriched in genes important for protein and organelle biosynthesis. Our results thus highlight examples of translational heterogeneity that are encoded in the genomic sequences and tuned to optimizing cellular homeostasis. More generally, our work highlights the power of Ribo-seq to understand the complexities of translation regulation.

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.

Riboexp: an interpretable reinforcement learning framework for ribosome density modeling

Translation elongation is a crucial phase during protein biosynthesis. In this study, we develop a novel deep reinforcement learning-based framework, named Riboexp, to model the determinants of the uneven distribution of ribosomes on mRNA transcripts during translation elongation. In particular, our model employs a policy network to perform a context-dependent feature selection in the setting of ribosome density prediction. Our extensive tests demonstrated that Riboexp can significantly outperform the state-of-the-art methods in predicting ribosome density by up to 5.9% in terms of per-gene Pearson correlation coefficient on the datasets from three species. In addition, Riboexp can indicate more informative sequence features for the prediction task than other commonly used attribution methods in deep learning. In-depth analyses also revealed the meaningful biological insights generated by the Riboexp framework. Moreover, the application of Riboexp in codon optimization resulted in an increase of protein production by around 31% over the previous state-of-the-art method that models ribosome density. These results have established Riboexp as a powerful and useful computational tool in the studies of translation dynamics and protein synthesis. Availability: The data and code of this study are available on GitHub: https://github.com/Liuxg16/Riboexp. Contact:  zengjy321@tsinghua.edu.cn; songsen@tsinghua.edu.cn.

Codon Resolution Analysis of Ribosome Profiling Data

Translation is a central biological process in living cells. Ribosome profiling approach enables assessing translation on a global, cell-wide level. Extracting versatile information from the ribosome profiling data usually requires specialized expertise for handling the sequencing data that is not available to the broad community of experimentalists. Here, we provide an easy-to-use and modifiable workflow that uses a small set of commands and enables full data analysis in a standardized way, including precise positioning of the ribosome-protected fragments, for determining codon-specific translation features. The workflow is complemented with simple step-by-step explanations and is accessible to scientists with no computational background.

RiboDiPA: a novel tool for differential pattern analysis in Ribo-seq data

Ribosome profiling, also known as Ribo-seq, has become a popular approach to investigate regulatory mechanisms of translation in a wide variety of biological contexts. Ribo-seq not only provides a measurement of translation efficiency based on the relative abundance of ribosomes bound to transcripts, but also has the capacity to reveal dynamic and local regulation at different stages of translation based on positional information of footprints across individual transcripts. While many computational tools exist for the analysis of Ribo-seq data, no method is currently available for rigorous testing of the pattern differences in ribosome footprints. In this work, we develop a novel approach together with an R package, RiboDiPA, for Differential Pattern Analysis of Ribo-seq data. RiboDiPA allows for quick identification of genes with statistically significant differences in ribosome occupancy patterns for model organisms ranging from yeast to mammals. We show that differential pattern analysis reveals information that is distinct and complimentary to existing methods that focus on translational efficiency analysis. Using both simulated Ribo-seq footprint data and three benchmark data sets, we illustrate that RiboDiPA can uncover meaningful pattern differences across multiple biological conditions on a global scale, and pinpoint characteristic ribosome occupancy patterns at single codon resolution.

What determines eukaryotic translation elongation: recent molecular and quantitative analyses of protein synthesis

Protein synthesis from mRNA is an energy-intensive and tightly controlled cellular process. Translation elongation is a well-coordinated, multifactorial step in translation that undergoes dynamic regulation owing to cellular state and environmental determinants. Recent studies involving genome-wide approaches have uncovered some crucial aspects of translation elongation including the mRNA itself and the nascent polypeptide chain. Additionally, these studies have fuelled quantitative and mathematical modelling of translation elongation. In this review, we provide a comprehensive overview of the key determinants of translation elongation. We discuss consequences of ribosome stalling or collision, and how the cells regulate translation in case of such events. Next, we review theoretical approaches and widely used mathematical models that have become an essential ingredient to interpret complex molecular datasets and study translation dynamics quantitatively. Finally, we review recent advances in live-cell reporter and related analysis techniques, to monitor the translation dynamics of single cells and single-mRNA molecules in real time.

RiboToolkit: an integrated platform for analysis and annotation of ribosome profiling data to decode mRNA translation at codon resolution

Ribosome profiling (Ribo-seq) is a powerful technology for globally monitoring RNA translation; ranging from codon occupancy profiling, identification of actively translated open reading frames (ORFs), to the quantification of translational efficiency under various physiological or experimental conditions. However, analyzing and decoding translation information from Ribo-seq data is not trivial. Although there are many existing tools to analyze Ribo-seq data, most of these tools are designed for specific or limited functionalities and an easy-to-use integrated tool to analyze Ribo-seq data is lacking. Fortunately, the small size (26–34 nt) of ribosome protected fragments (RPFs) in Ribo-seq and the relatively small amount of sequencing data greatly facilitates the development of such a web platform, which is easy to manipulate for users with or without bioinformatic expertise. Thus, we developed RiboToolkit (http://rnabioinfor.tch.harvard.edu/RiboToolkit), a convenient, freely available, web-based service to centralize Ribo-seq data analyses, including data cleaning and quality evaluation, expression analysis based on RPFs, codon occupancy, translation efficiency analysis, differential translation analysis, functional annotation, translation metagene analysis, and identification of actively translated ORFs. Besides, easy-to-use web interfaces were developed to facilitate data analysis and intuitively visualize results. Thus, RiboToolkit will greatly facilitate the study of mRNA translation based on ribosome profiling.

Genetic circuit characterization by inferring RNA polymerase movement and ribosome usage

To perform their computational function, genetic circuits change states through a symphony of genetic parts that turn regulator expression on and off. Debugging is frustrated by an inability to characterize parts in the context of the circuit and identify the origins of failures. Here, we take snapshots of a large genetic circuit in different states: RNA-seq is used to visualize circuit function as a changing pattern of RNA polymerase (RNAP) flux along the DNA. Together with ribosome profiling, all 54 genetic parts (promoters, ribozymes, RBSs, terminators) are parameterized and used to inform a mathematical model that can predict circuit performance, dynamics, and robustness. The circuit behaves as designed; however, it is riddled with genetic errors, including cryptic sense/antisense promoters and translation, attenuation, incorrect start codons, and a failed gate. While not impacting the expected Boolean logic, they reduce the prediction accuracy and could lead to failures when the parts are used in other designs. Finally, the cellular power (RNAP and ribosome usage) required to maintain a circuit state is calculated. This work demonstrates the use of a small number of measurements to fully parameterize a regulatory circuit and quantify its impact on host.

Accurate detection of short and long active ORFs using Ribo-seq data

MOTIVATION

Ribo-seq, a technique for deep-sequencing ribosome-protected mRNA fragments, has enabled transcriptome-wide monitoring of translation in vivo. It has opened avenues for re-evaluating the coding potential of open reading frames (ORFs), including many short ORFs that were previously presumed to be non-translating. However, the detection of translating ORFs, specifically short ORFs, from Ribo-seq data, remains challenging due to its high heterogeneity and noise.

RESULTS

We present ribotricer, a method for detecting actively translating ORFs by directly leveraging the three-nucleotide periodicity of Ribo-seq data. Ribotricer demonstrates higher accuracy and robustness compared with other methods at detecting actively translating ORFs including short ORFs on multiple published datasets across species inclusive of Arabidopsis, Caenorhabditis elegans, Drosophila, human, mouse, rat, yeast and zebrafish.

AVAILABILITY AND IMPLEMENTATION

Ribotricer is available at https://github.com/smithlabcode/ribotricer. All analysis scripts and results are available at https://github.com/smithlabcode/ribotricer-results.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Translation initiation downstream from annotated start codons in human mRNAs coevolves with the Kozak context

Eukaryotic translation initiation involves preinitiation ribosomal complex 5′-to-3′ directional probing of mRNA for codons suitable for starting protein synthesis. The recognition of codons as starts depends on the codon identity and on its immediate nucleotide context known as Kozak context. When the context is weak (i.e., nonoptimal), leaky scanning takes place during which a fraction of ribosomes continues the mRNA probing. We explored the relationship between the context of AUG codons annotated as starts of protein-coding sequences and the next AUG codon occurrence. We found that AUG codons downstream from weak starts occur in the same frame more frequently than downstream from strong starts. We suggest that evolutionary selection on in-frame AUGs downstream from weak start codons is driven by the advantage of the reduction of wasteful out-of-frame product synthesis and also by the advantage of producing multiple proteoforms from certain mRNAs. We confirmed translation initiation downstream from weak start codons using ribosome profiling data. We also tested translation of alternative start codons in 10 specific human genes using reporter constructs. In all tested cases, initiation at downstream start codons was more productive than at the annotated ones. In most cases, optimization of Kozak context did not completely abolish downstream initiation, and in the specific example of CMPK1 mRNA, the optimized start remained unproductive. Collectively, our work reveals previously uncharacterized forces shaping the evolution of protein-coding genes and points to the plurality of translation initiation and the existence of sequence features influencing start codon selection, other than Kozak context.

Analysis of computational codon usage models and their association with translationally slow codons

Improved computational modeling of protein translation rates, including better prediction of where translational slowdowns along an mRNA sequence may occur, is critical for understanding co-translational folding. Because codons within a synonymous codon group are translated at different rates, many computational translation models rely on analyzing synonymous codons. Some models rely on genome-wide codon usage bias (CUB), believing that globally rare and common codons are the most informative of slow and fast translation, respectively. Others use the CUB observed only in highly expressed genes, which should be under selective pressure to be translated efficiently (and whose CUB may therefore be more indicative of translation rates). No prior work has analyzed these models for their ability to predict translational slowdowns. Here, we evaluate five models for their association with slowly translated positions as denoted by two independent ribosome footprint (RFP) count experiments from S. cerevisiae, because RFP data is often considered as a “ground truth” for translation rates across mRNA sequences. We show that all five considered models strongly associate with the RFP data and therefore have potential for estimating translational slowdowns. However, we also show that there is a weak correlation between RFP counts for the same genes originating from independent experiments, even when their experimental conditions are similar. This raises concerns about the efficacy of using current RFP experimental data for estimating translation rates and highlights a potential advantage of using computational models to understand translation rates instead.

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.

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.

Ribosome pausing, a dangerous necessity for co-translational events

In recent years translation elongation has emerged as an important contributor to the regulation of gene expression. There are multiple quality control checkpoints along the way of producing mature proteins and targeting them to the right cellular compartment, or associating them correctly with their partners. Ribosomes pause to allow co-translational protein folding, protein targeting or protein interactions, and the pausing is dictated by a combination of the mRNA sequence and structure, the tRNA availability and the nascent peptide. However, ribosome pausing can also lead to ribosome collisions and co-translational degradation of both mRNA and nascent chain. Understanding how the translating ribosome tunes the different maturation steps that nascent proteins must undergo, what the timing of these maturation events is, and how degradation can be avoided when pausing is needed, is now possible by the emergence of methods to follow ribosome dynamics in vivo. This review summarizes some of the recent studies that have advanced our knowledge about co-translational events using the power of ribosome profiling, and some of the questions that have emerged from these studies.

Generally applicable transcriptome-wide analysis of translation using anota2seq

mRNA translation plays an evolutionarily conserved role in homeostasis and when dysregulated contributes to various disorders including metabolic and neurological diseases and cancer. Notwithstanding that optimal and universally applicable methods are critical for understanding the complex role of translational control under physiological and pathological conditions, approaches to analyze translatomes are largely underdeveloped. To address this, we developed the anota2seq algorithm which outperforms current methods for statistical identification of changes in translation. Notably, in contrast to available analytical methods, anota2seq also allows specific identification of an underappreciated mode of gene expression regulation whereby translation acts as a buffering mechanism which maintains protein levels despite fluctuations in corresponding mRNA abundance (‘translational buffering’). Thus, the universal anota2seq algorithm allows efficient and hitherto unprecedented interrogation of translatomes which is anticipated to advance knowledge regarding the role of translation in homeostasis and disease.

XPRESSyourself: Enhancing, standardizing, and automating ribosome profiling computational analyses yields improved insight into data

Ribosome profiling, an application of nucleic acid sequencing for monitoring ribosome activity, has revolutionized our understanding of protein translation dynamics. This technique has been available for a decade, yet the current state and standardization of publicly available computational tools for these data is bleak. We introduce XPRESSyourself, an analytical toolkit that eliminates barriers and bottlenecks associated with this specialized data type by filling gaps in the computational toolset for both experts and non-experts of ribosome profiling. XPRESSyourself automates and standardizes analysis procedures, decreasing time-to-discovery and increasing reproducibility. This toolkit acts as a reference implementation of current best practices in ribosome profiling analysis. We demonstrate this toolkit’s performance on publicly available ribosome profiling data by rapidly identifying hypothetical mechanisms related to neurodegenerative phenotypes and neuroprotective mechanisms of the small-molecule ISRIB during acute cellular stress. XPRESSyourself brings robust, rapid analysis of ribosome-profiling data to a broad and ever-expanding audience and will lead to more reproducible and accessible measurements of translation regulation. XPRESSyourself software is perpetually open-source under the GPL-3.0 license and is hosted at https://github.com/XPRESSyourself, where users can access additional documentation and report software issues.

TDP-43 enhances translation of specific mRNAs linked to neurodegenerative disease

The RNA-binding protein TDP-43 is heavily implicated in neurodegenerative disease. Numerous patient mutations in TARDBP, the gene encoding TDP-43, combined with data from animal and cell-based models, imply that altered RNA regulation by TDP-43 causes Amyotrophic Lateral Sclerosis and Frontotemporal Dementia. However, underlying mechanisms remain unresolved. Increased cytoplasmic TDP-43 levels in diseased neurons suggest a possible role in this cellular compartment. Here, we examined the impact on translation of overexpressing human TDP-43 and the TDP-43A315T patient mutant protein in motor neuron-like cells and primary cultures of cortical neurons. In motor-neuron like cells, TDP-43 associates with ribosomes without significantly affecting global translation. However, ribosome profiling and additional assays revealed enhanced translation and direct binding of Camta1, Mig12, and Dennd4a mRNAs. Overexpressing either wild-type TDP-43 or TDP-43A315T stimulated translation of Camta1 and Mig12 mRNAs via their 5'UTRs and increased CAMTA1 and MIG12 protein levels. In contrast, translational enhancement of Dennd4a mRNA required a specific 3'UTR region and was specifically observed with the TDP-43A315T patient mutant allele. Our data reveal that TDP-43 can function as an mRNA-specific translational enhancer. Moreover, since CAMTA1 and DENND4A are linked to neurodegeneration, they suggest that this function could contribute to disease.

GWIPS-viz: 2018 update

The GWIPS-viz browser (http://gwips.ucc.ie/) is an on-line genome browser which is tailored for exploring ribosome profiling (Ribo-seq) data. Since its publication in 2014, GWIPS-viz provides Ribo-seq data for an additional 14 genomes bringing the current total to 23. The integration of new Ribo-seq data has been automated thereby increasing the number of available tracks to 1792, a 10-fold increase in the last three years. The increase is particularly substantial for data derived from human sources. Following user requests, we added the functionality to download these tracks in bigWig format. We also incorporated new types of data (e.g. TCP-seq) as well as auxiliary tracks from other sources that help with the interpretation of Ribo-seq data. Improvements in the visualization of the data have been carried out particularly for bacterial genomes where the Ribo-seq data are now shown in a strand specific manner. For higher eukaryotic datasets, we provide characteristics of individual datasets using the RUST program which includes the triplet periodicity, sequencing biases and relative inferred A-site dwell times. This information can be used for assessing the quality of Ribo-seq datasets. To improve the power of the signal, we aggregate Ribo-seq data from several studies into Global aggregate tracks for each genome.

PROTEOFORMER 2.0: Further Developments in the Ribosome Profiling-assisted Proteogenomic Hunt for New Proteoforms

PROTEOFORMER is a pipeline that enables the automated processing of data derived from ribosome profiling (RIBO-seq, i.e. the sequencing of ribosome-protected mRNA fragments). As such, genome-wide ribosome occupancies lead to the delineation of data-specific translation product candidates and these can improve the mass spectrometry-based identification. Since its first publication, different upgrades, new features and extensions have been added to the PROTEOFORMER pipeline. Some of the most important upgrades include P-site offset calculation during mapping, comprehensive data pre-exploration, the introduction of two alternative proteoform calling strategies and extended pipeline output features. These novelties are illustrated by analyzing ribosome profiling data of human HCT116 and Jurkat data. The different proteoform calling strategies are used alongside one another and in the end combined together with reference sequences from UniProt. Matching mass spectrometry data are searched against this extended search space with MaxQuant. Overall, besides annotated proteoforms, this pipeline leads to the identification and validation of different categories of new proteoforms, including translation products of up- and downstream open reading frames, 5' and 3' extended and truncated proteoforms, single amino acid variants, splice variants and translation products of so-called noncoding regions. Further, proof-of-concept is reported for the improvement of spectrum matching by including Prosit, a deep neural network strategy that adds extra fragmentation spectrum intensity features to the analysis. In the light of ribosome profiling-driven proteogenomics, it is shown that this allows validating the spectrum matches of newly identified proteoforms with elevated stringency. These updates and novel conclusions provide new insights and lessons for the ribosome profiling-based proteogenomic research field. More practical information on the pipeline, raw code, the user manual (README) and explanations on the different modes of availability can be found at the GitHub repository of PROTEOFORMER: https://github.com/Biobix/proteoformer.

Computational resources for ribosome profiling: from database to Web server and software

Ribosome profiling is emerging as a powerful technique that enables genome-wide investigation of in vivo translation at sub-codon resolution. The increasing application of ribosome profiling in recent years has achieved remarkable progress toward understanding the composition, regulation and mechanism of translation. This benefits from not only the awesome power of ribosome profiling but also an extensive range of computational resources available for ribosome profiling. At present, however, a comprehensive review on these resources is still lacking. Here, we survey the recent computational advances guided by ribosome profiling, with a focus on databases, Web servers and software tools for storing, visualizing and analyzing ribosome profiling data. This review is intended to provide experimental and computational biologists with a reference to make appropriate choices among existing resources for the question at hand.

RiboStreamR: a web application for quality control, analysis, and visualization of Ribo-seq data

BACKGROUND

Ribo-seq is a popular technique for studying translation and its regulation. A Ribo-seq experiment produces a snap-shot of the location and abundance of actively translating ribosomes within a cell's transcriptome. In practice, Ribo-seq data analysis can be sensitive to quality issues such as read length variation, low read periodicities, and contaminations with ribosomal and transfer RNA. Various software tools for data preprocessing, quality assessment, analysis, and visualization of Ribo-seq data have been developed. However, many of these tools require considerable practical knowledge of software applications, and often multiple different tools have to be used in combination with each other.

RESULTS

We present riboStreamR, a comprehensive Ribo-seq quality control (QC) platform in the form of an R Shiny web application. RiboStreamR provides visualization and analysis tools for various Ribo-seq QC metrics, including read length distribution, read periodicity, and translational efficiency. Our platform is focused on providing a user-friendly experience, and includes various options for graphical customization, report generation, and anomaly detection within Ribo-seq datasets.

CONCLUSIONS

RiboStreamR takes advantage of the vast resources provided by the R and Bioconductor environments, and utilizes the Shiny R package to ensure a high level of usability. Our goal is to develop a tool which facilitates in-depth quality assessment of Ribo-seq data by providing reference datasets and automatically highlighting quality issues and anomalies within datasets.

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.

Optimization of ribosome profiling using low-input brain tissue from fragile X syndrome model mice

Dysregulated protein synthesis is a major underlying cause of many neurodevelopmental diseases including fragile X syndrome. In order to capture subtle but biologically significant differences in translation in these disorders, a robust technique is required. One powerful tool to study translational control is ribosome profiling, which is based on deep sequencing of mRNA fragments protected from ribonuclease (RNase) digestion by ribosomes. However, this approach has been mainly applied to rapidly dividing cells where translation is active and large amounts of starting material are readily available. The application of ribosome profiling to low-input brain tissue where translation is modest and gene expression changes between genotypes are expected to be small has not been carefully evaluated. Using hippocampal tissue from wide type and fragile X mental retardation 1 (Fmr1) knockout mice, we show that variable RNase digestion can lead to significant sample batch effects. We also establish GC content and ribosome footprint length as quality control metrics for RNase digestion. We performed RNase titration experiments for low-input samples to identify optimal conditions for this critical step that is often improperly conducted. Our data reveal that optimal RNase digestion is essential to ensure high quality and reproducibility of ribosome profiling for low-input brain tissue.

Quantification and discovery of sequence determinants of protein-per-mRNA amount in 29 human tissues

Despite their importance in determining protein abundance, a comprehensive catalogue of sequence features controlling protein-to-mRNA (PTR) ratios and a quantification of their effects are still lacking. Here, we quantified PTR ratios for 11,575 proteins across 29 human tissues using matched transcriptomes and proteomes. We estimated by regression the contribution of known sequence determinants of protein synthesis and degradation in addition to 45 mRNA and 3 protein sequence motifs that we found by association testing. While PTR ratios span more than 2 orders of magnitude, our integrative model predicts PTR ratios at a median precision of 3.2-fold. A reporter assay provided functional support for two novel UTR motifs, and an immobilized mRNA affinity competition-binding assay identified motif-specific bound proteins for one motif. Moreover, our integrative model led to a new metric of codon optimality that captures the effects of codon frequency on protein synthesis and degradation. Altogether, this study shows that a large fraction of PTR ratio variation in human tissues can be predicted from sequence, and it identifies many new candidate post-transcriptional regulatory elements.

A systematically-revised ribosome profiling method for bacteria reveals pauses at single-codon resolution

In eukaryotes, ribosome profiling provides insight into the mechanism of protein synthesis at the codon level. In bacteria, however, the method has been more problematic and no consensus has emerged for how to best prepare profiling samples. Here, we identify the sources of these problems and describe new solutions for arresting translation and harvesting cells in order to overcome them. These improvements remove confounding artifacts and improve the resolution to allow analyses of ribosome behavior at the codon level. With a clearer view of the translational landscape in vivo, we observe that filtering cultures leads to translational pauses at serine and glycine codons through the reduction of tRNA aminoacylation levels. This observation illustrates how bacterial ribosome profiling studies can yield insight into the mechanism of protein synthesis at the codon level and how these mechanisms are regulated in response to changes in the physiology of the cell.

RiboProP: a probabilistic ribosome positioning algorithm for ribosome profiling

Motivation

Ribosome profiling has been widely used to study translation in a genome-wide fashion. It requires deep sequencing of ribosome protected mRNA fragments followed by mapping of fragments to the reference genome. For applications such as identification of ribosome pausing sites, it is not enough to map a fragment to a given gene, but the exact position of the ribosome represented by the fragment must be identified for each mRNA fragment. The assignment of the correct ribosome position is complicated by the broad length distribution of the ribosome protected fragments caused by the known sequence bias of micrococcal nuclease (MNase), the most widely used nuclease for digesting mRNAs in bacteria. Available mapping algorithms suffer from either MNase bias or low accuracy in characterizing the ribosome pausing kinetics.

Results

In this paper, we introduce a new computational method for mapping the ribosome protected fragments to ribosome locations. We first develop a mathematical model of the interplay between MNase digestion and ribosome protection of the mRNAs. We then use the model to reconstruct the ribosome occupancy profile on a per gene level. We demonstrate that our method has the capability of mitigating the sequence bias introduced by MNase and accurately locating ribosome pausing sites at codon resolution. We believe that our method can be broadly applied to ribosome profiling studies on bacteria where codon resolution is necessary.

Availability and implementation

Source code implementing our approach can be downloaded under GPL3 license at http://bioserv.mps.ohio-state.edu/RiboProP.

Supplementary information

Supplementary data are available at Bioinformatics online.