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

Long-read RNA sequencing reveals allele-specific N 6-methyladenosine modifications

Long-read sequencing technology enables highly accurate detection of allele-specific RNA expression, providing insights into the effects of genetic variation on splicing and RNA abundance. Furthermore, the ability to directly sequence RNA enables the detection of RNA modifications in tandem with ascertaining the allelic origin of each molecule. Here, we leverage these advantages to determine allele-biased patterns of N 6-methyladenosine (m6A) modifications in native mRNA. We used human and mouse cells with known genetic variants to assign the allelic origin of each mRNA molecule combined with a supervised machine learning model to detect read-level m6A modification ratios. Our analyses reveal the importance of sequences adjacent to the DRACH motif in determining m6A deposition, in addition to allelic differences that directly alter the motif. Moreover, we discover allele-specific m6A modification events with no genetic variants in close proximity to the differentially modified nucleotide, demonstrating the unique advantage of using long-reads and surpassing the capabilities of antibody-based short-read approaches. This technological advance will further our understanding of the role of genetics in determining mRNA modifications.

Landscape and regulation of mRNA translation in the early C. elegans embryo

Animal embryos rely on regulated translation of maternally deposited mRNAs to drive early development. Using low-input ribosome profiling combined with RNA sequencing on precisely staged embryos, we measured mRNA translation during the first four cell cycles of C. elegans development. We uncovered stage-specific patterns of developmentally coordinated translational regulation. We confirmed that mRNA localization correlates with translational eLiciency, though initial translational repression in germline precursors occurs before P-granule association. Our analysis suggests that the RNA-binding protein OMA-1 represses the translation of its target mRNAs in a stage-specific manner, while indirectly promoting the translational eLiciency of other transcripts. These findings illuminate how post-transcriptional mechanisms shape the embryonic proteome to direct cell diLerentiation, with implications for understanding similar regulation across species where maternal factors guide early development.

Differential impacts of ribosomal protein haploinsufficiency on mitochondrial function

The interplay between ribosomal protein (RP) composition and mitochondrial function is essential for energy homeostasis. Balanced RP production optimizes protein synthesis while minimizing energy costs, but its impact on mitochondrial functionality remains unclear. Here, we investigated haploinsufficiency for RP genes (rps-10, rpl-5, rpl-33, and rps-23) in Caenorhabditis elegans and corresponding reductions in human lymphoblast cells. Significant mitochondrial morphological differences, upregulation of glutathione transferases, and SKN-1-dependent oxidative stress resistance were observed across mutants. Loss of a Datasingle rps-10 copy reduced mitochondrial activity, energy levels, and oxygen consumption, mirrored by similar reductions in mitochondrial activity and energy levels in lymphoblast cells with 50% lower RPS10 transcripts. Both systems exhibited altered translation efficiency (TE) of mitochondrial electron transport chain components, suggesting a conserved mechanism to adjust mitochondrial protein synthesis under ribosomal stress. Finally, mitochondrial membrane and cytosolic RPs showed significant RNA and TE covariation in lymphoblastoid cells, highlighting the interplay between protein synthesis machinery and mitochondrial energy production.

The assembly factor Reh1 is released from the ribosome during its initial round of translation

Assembly of functional ribosomal subunits and successfully delivering them to the translating pool is a prerequisite for protein synthesis and cell growth. In S. cerevisiae, the ribosome assembly factor Reh1 binds to pre-60S subunits at a late stage during their cytoplasmic maturation. Previous work shows that the C-terminus of Reh1 inserts into the polypeptide exit tunnel of the pre-60S subunit. Here, we show that Reh1-bound nascent 60S subunits associate with 40S subunits to form actively translating ribosomes. Using selective ribosome profiling, we found that Reh1-bound ribosomes populate open reading frames near start codons. Reh1-bound ribosomes are also strongly enriched for initiator tRNA, indicating they are associated with early elongation. Using cryo-electron microscopy to image Reh1-bound 80S ribosomes, we found they contain A site peptidyl tRNA, P site tRNA and eIF5A, indicating that Reh1 does not dissociate from 60S until translation elongation. We propose that Reh1 is displaced by the elongating peptide chain, making it the last assembly factor released from the nascent 60S subunit during its initial round of translation.

Translation efficiency covariation across cell types is a conserved organizing principle of mammalian transcriptomes

Characterization of shared patterns of RNA expression between genes across conditions has led to the discovery of regulatory networks and novel biological functions. However, it is unclear if such coordination extends to translation, a critical step in gene expression. Here, we uniformly analyzed 3,819 ribosome profiling datasets from 117 human and 94 mouse tissues and cell lines. We introduce the concept of Translation Efficiency Covariation (TEC), identifying coordinated translation patterns across cell types. We nominate potential mechanisms driving shared patterns of translation regulation. TEC is conserved across human and mouse cells and helps uncover gene functions. Moreover, our observations indicate that proteins that physically interact are highly enriched for positive covariation at both translational and transcriptional levels. Our findings establish translational covariation as a conserved organizing principle of mammalian transcriptomes.

Predicting the translation efficiency of messenger RNA in mammalian cells

The degree to which translational control is specified by mRNA sequence is poorly understood in mammalian cells. Here, we constructed and leveraged a compendium of 3,819 ribosomal profiling datasets, distilling them into a transcriptome-wide atlas of translation efficiency (TE) measurements encompassing >140 human and mouse cell types. We subsequently developed RiboNN, a multitask deep convolutional neural network, and classic machine learning models to predict TEs in hundreds of cell types from sequence-encoded mRNA features, achieving state-of-the-art performance (r=0.79 in human and r=0.78 in mouse for mean TE across cell types). While the majority of earlier models solely considered 5' UTR sequence1, RiboNN integrates contributions from the full-length mRNA sequence, learning that the 5' UTR, CDS, and 3' UTR respectively possess ~67%, 31%, and 2% per-nucleotide information density in the specification of mammalian TEs. Interpretation of RiboNN revealed that the spatial positioning of low-level di- and tri-nucleotide features (i.e., including codons) largely explain model performance, capturing mechanistic principles such as how ribosomal processivity and tRNA abundance control translational output. RiboNN is predictive of the translational behavior of base-modified therapeutic RNA, and can explain evolutionary selection pressures in human 5' UTRs. Finally, it detects a common language governing mRNA regulatory control and highlights the interconnectedness of mRNA translation, stability, and localization in mammalian organisms.

Long-read RNA sequencing reveals allele-specific N6-methyladenosine modifications

Long-read sequencing technology enables highly accurate detection of allele-specific RNA expression, providing insights into the effects of genetic variation on splicing and RNA abundance. Furthermore, the ability to directly sequence RNA promises the detection of RNA modifications in tandem with ascertaining the allelic origin of each molecule. Here, we leverage these advantages to determine allele-biased patterns of N6-methyladenosine (m6A) modifications in native mRNA. We utilized human and mouse cells with known genetic variants to assign allelic origin of each mRNA molecule combined with a supervised machine learning model to detect read-level m6A modification ratios. Our analyses revealed the importance of sequences adjacent to the DRACH-motif in determining m6A deposition, in addition to allelic differences that directly alter the motif. Moreover, we discovered allele-specific m6A modification (ASM) events with no genetic variants in close proximity to the differentially modified nucleotide, demonstrating the unique advantage of using long reads and surpassing the capabilities of antibody-based short-read approaches. This technological advancement promises to advance our understanding of the role of genetics in determining mRNA modifications.

RiboGraph: an interactive visualization system for ribosome profiling data at read length resolution

MOTIVATION

Ribosome profiling is a widely-used technique for measuring ribosome occupancy at nucleotide resolution. However, the need to analyze this data at nucleotide resolution introduces unique challenges in data visualization and analyses.

RESULTS

In this study, we introduce RiboGraph, a dedicated visualization tool designed to work with .ribo files, a specialized and efficient format for ribosome occupancy data. Unlike existing solutions that rely on large alignment files and time-consuming preprocessing steps, RiboGraph operates on a purpose designed compact file type. This efficiency allows for interactive, real-time visualization at ribosome-protected fragment length resolution. By providing an integrated toolset, RiboGraph empowers researchers to conduct comprehensive visual analysis of ribosome occupancy data.

AVAILABILITY AND IMPLEMENTATION

Source code, step-by-step installation instructions and links to documentation are available on GitHub: https://github.com/ribosomeprofiling/ribograph. On the same page, we provide test files and a step-by-step tutorial highlighting the key features of RiboGraph.

Cytosolic Ribosomal Protein Haploinsufficiency affects Mitochondrial Morphology and Respiration

The interplay between ribosomal protein composition and mitochondrial function is essential for sustaining energy homeostasis. Precise stoichiometric production of ribosomal proteins is crucial to maximize protein synthesis efficiency while reducing the energy costs to the cell. However, the impact of this balance on mitochondrial ATP generation, morphology and function remains unclear. Particularly, the loss of a single copy ribosomal protein gene is observed in Mendelian disorders like Diamond Blackfan Anemia and is common in somatic tumors, yet the implications of this imbalance on mitochondrial function and energy dynamics are still unclear. In this study, we investigated the impact of haploinsufficiency for four ribosomal protein genes implicated in ribosomopathy disorders (rps-10, rpl-5, rpl-33, rps-23) in Caenorhabditis elegans and corresponding reductions in human lymphoblast cells. Our findings uncover significant, albeit variably penetrant, mitochondrial morphological differences across these mutants, alongside an upregulation of glutathione transferases, and SKN-1 dependent increase in oxidative stress resistance, indicative of increased ROS production. Specifically, loss of a single copy of rps-10 in C. elegans led to decreased mitochondrial activity, characterized by lower energy levels and reduced oxygen consumption. A similar reduction in mitochondrial activity and energy levels was observed in human leukemia cells with a 50% reduction in RPS10 transcript levels. Importantly, we also observed alterations in the translation efficiency of nuclear and mitochondrial electron transport chain components in response to reductions in ribosomal protein genes' expression in both C. elegans and human cells. This suggests a conserved mechanism whereby the synthesis of components vital for mitochondrial function are adjusted in the face of compromised ribosomal machinery. Finally, mitochondrial membrane and cytosolic ribosomal components exhibited significant covariation at the RNA and translation efficiency level in lymphoblastoid cells across a diverse group of individuals, emphasizing the interplay between the protein synthesis machinery and mitochondrial energy production. By uncovering the impact of ribosomal protein haploinsufficiency on the translation efficiency of electron transport chain components, mitochondrial physiology, and the adaptive stress responses, we provide evidence for an evolutionarily conserved strategy to safeguard cellular functionality under genetic stress.

Integrated multiplexed assays of variant effect reveal determinants of catechol-O-methyltransferase gene expression

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase or decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

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.

Translatome and Transcriptome Analyses Reveal the Mechanism that Underlies the Enhancement of Salt Stress by the Small Peptide Ospep5 in Plants

Salt stress significantly impedes plant growth and the crop yield. This study utilized de novo transcriptome assembly and ribosome profiling to explore mRNA translation's role in rice salt tolerance. We identified unrecognized translated open reading frames (ORFs), including 42 upstream transcripts and 86 unannotated transcripts. A noteworthy discovery was the role of a small ORF, Ospep5, in conferring salt tolerance. Overexpression of Ospep5 in plants increased salt tolerance, while its absence led to heightened sensitivity. This hypothesis was corroborated by the findings that exogenous application of the synthetic small peptide Ospep5 bolstered salt tolerance in both rice and Arabidopsis. We found that the mechanism underpinning the Ospep5-mediated salt tolerance involves the maintenance of intracellular Na+/K+ homeostasis, facilitated by upregulation of high-affinity potassium transporters (HKT) and Na+/H+ exchangers (SOS1). Furthermore, a comprehensive multiomics approach, particularly ribosome profiling, is instrumental in uncovering unannotated ORFs and elucidating their functions in plant stress responses.

Cellular targets and lysine selectivity of the HERC5 ISG15 ligase

ISG15 is a type I interferon-induced ubiquitin-like modifier that functions in innate immune responses. The major human ISG15 ligase is hHERC5, a ribosome-associated HECT E3 that broadly ISGylates proteins cotranslationally. Here, we characterized the hHERC5-dependent ISGylome and identified over 2,000 modified lysines in over 1,100 proteins in IFN-β-stimulated cells. In parallel, we compared the substrate selectivity hHERC5 to the major mouse ISG15 ligase, mHERC6, and analysis of sequences surrounding ISGylation sites revealed that hHERC5 and mHERC6 have distinct preferences for amino acid sequence context. Several features of the datasets were consistent with ISGylation of ribosome-tethered nascent chains, and mHERC6, like hHERC5, cotranslationally modified nascent polypeptides. The ISGylome datasets presented here represent the largest numbers of protein targets and modification sites attributable to a single Ub/Ubl ligase and the lysine selectivities of the hHERC5 and mHERC6 enzymes may have implications for the activities of HECT domain ligases, generally.

RiboGraph: An interactive visualization system for ribosome profiling data at read length resolution

Ribosome profiling is a widely-used technique for measuring ribosome occupancy at nucleotide resolution. However, the need to analyze this data at nucleotide resolution introduces unique challenges in data visualization and analyses. In this study, we introduce RiboGraph, a dedicated visualization tool designed to work with .ribo files, a specialized and efficient format for ribosome occupancy data. Unlike existing solutions that rely on large alignment files and time-consuming preprocessing steps, RiboGraph operates on a purpose designed compact file type and eliminates the need for data preprocessing. This efficiency allows for interactive, real-time visualization at ribosome-protected fragment length resolution. By providing an integrated toolset, RiboGraph empowers researchers to conduct comprehensive visual analysis of ribosome occupancy data.

Availability and Implementation

Source code, step-by-step installation instructions and links to documentation are available on GitHub: https://github.com/ribosomeprofiling/ribograph. On the same page, we provide test files and a step-by-step tutorial highlighting the key features of RiboGraph.

Integrated multiplexed assays of variant effect reveal cis-regulatory determinants of catechol-O-methyltransferase gene expression

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase and decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

The Ribosome Assembly Factor Reh1 is Released from the Polypeptide Exit Tunnel in the Pioneer Round of Translation

Assembly of functional ribosomal subunits and successfully delivering them to the translating pool is a prerequisite for protein synthesis and cell growth. In S. cerevisiae, the ribosome assembly factor Reh1 binds to pre-60S subunits at a late stage during their cytoplasmic maturation. Previous work shows that the C-terminus of Reh1 inserts into the polypeptide exit tunnel (PET) of the pre-60S subunit. Unlike canonical assembly factors, which associate exclusively with pre-60S subunits, we observed that Reh1 sediments with polysomes in addition to free 60S subunits. We therefore investigated the intriguing possibility that Reh1 remains associated with 60S subunits after the release of the anti-association factor Tif6 and after subunit joining. Here, we show that Reh1-bound nascent 60S subunits associate with 40S subunits to form actively translating ribosomes. Using selective ribosome profiling, we found that Reh1-bound ribosomes populate open reading frames near start codons. Reh1-bound ribosomes are also strongly enriched for initiator tRNA, indicating they are associated with early elongation events. Using single particle cryo-electron microscopy to image cycloheximide-arrested Reh1-bound 80S ribosomes, we found that Reh1-bound 80S contain A site peptidyl tRNA, P site tRNA and eIF5A indicating that Reh1 does not dissociate from 60S until early stages of translation elongation. We propose that Reh1 is displaced by the elongating peptide chain. These results identify Reh1 as the last assembly factor released from the nascent 60S subunit during its pioneer round of translation.

Inhibition of ribosome biogenesis in the epidermis is sufficient to trigger organism-wide growth quiescence independently of nutritional status in C. elegans

Interorgan communication is crucial for multicellular organismal growth, development, and homeostasis. Cell nonautonomous inhibitory cues, which limit tissue-specific growth alterations, are not well characterized due to cell ablation approach limitations. In this study, we employed the auxin-inducible degradation system in C. elegans to temporally and spatially modulate ribosome biogenesis, through depletion of essential factors (RPOA-2, GRWD-1, or TSR-2). Our findings reveal that embryo-wide inhibition of ribosome biogenesis induces a reversible early larval growth quiescence, distinguished by a unique gene expression signature that is different from starvation or dauer stages. When ribosome biogenesis is inhibited in volumetrically similar tissues, including body wall muscle, epidermis, pharynx, intestine, or germ line, it results in proportionally stunted growth across the organism to different degrees. We show that specifically inhibiting ribosome biogenesis in the epidermis is sufficient to trigger an organism-wide growth quiescence. Epidermis-specific ribosome depletion leads to larval growth quiescence at the L3 stage, reduces organism-wide protein synthesis, and induced cell nonautonomous gene expression alterations. Further molecular analysis reveals overexpression of secreted proteins, suggesting an organism-wide regulatory mechanism. We find that UNC-31, a dense-core vesicle (DCV) pathway component, plays a significant role in epidermal ribosome biogenesis-mediated growth quiescence. Our tissue-specific knockdown experiments reveal that the organism-wide growth quiescence induced by epidermal-specific ribosome biogenesis inhibition is suppressed by reducing unc-31 expression in the epidermis, but not in neurons or body wall muscles. Similarly, IDA-1, a membrane-associated protein of the DCV, is overexpressed, and its knockdown in epidermis suppresses the organism-wide growth quiescence in response to epidermal ribosome biogenesis inhibition. Finally, we observe an overall increase in DCV puncta labeled by IDA-1 when epidermal ribosome biogenesis is inhibited, and these puncta are present in or near epidermal cells. In conclusion, these findings suggest a novel mechanism of nutrition-independent multicellular growth coordination initiated from the epidermis tissue upon ribosome biogenesis inhibition.

Pyruvate Kinase M (PKM) binds ribosomes in a poly-ADP ribosylation dependent manner to induce translational stalling

In light of the numerous studies identifying post-transcriptional regulators on the surface of the endoplasmic reticulum (ER), we asked whether there are factors that regulate compartment specific mRNA translation in human cells. Using a proteomic survey of spatially regulated polysome interacting proteins, we identified the glycolytic enzyme Pyruvate Kinase M (PKM) as a cytosolic (i.e. ER-excluded) polysome interactor and investigated how it influences mRNA translation. We discovered that the PKM-polysome interaction is directly regulated by ADP levels-providing a link between carbohydrate metabolism and mRNA translation. By performing enhanced crosslinking immunoprecipitation-sequencing (eCLIP-seq), we found that PKM crosslinks to mRNA sequences that are immediately downstream of regions that encode lysine- and glutamate-enriched tracts. Using ribosome footprint protection sequencing, we found that PKM binding to ribosomes causes translational stalling near lysine and glutamate encoding sequences. Lastly, we observed that PKM recruitment to polysomes is dependent on poly-ADP ribosylation activity (PARylation)-and may depend on co-translational PARylation of lysine and glutamate residues of nascent polypeptide chains. Overall, our study uncovers a novel role for PKM in post-transcriptional gene regulation, linking cellular metabolism and mRNA translation.

Dynamic interplay between RPL3- and RPL3L-containing ribosomes modulates mitochondrial activity in the mammalian heart

The existence of naturally occurring ribosome heterogeneity is now a well-acknowledged phenomenon. However, whether this heterogeneity leads to functionally diverse 'specialized ribosomes' is still a controversial topic. Here, we explore the biological function of RPL3L (uL3L), a ribosomal protein (RP) paralogue of RPL3 (uL3) that is exclusively expressed in skeletal muscle and heart tissues, by generating a viable homozygous Rpl3l knockout mouse strain. We identify a rescue mechanism in which, upon RPL3L depletion, RPL3 becomes up-regulated, yielding RPL3-containing ribosomes instead of RPL3L-containing ribosomes that are typically found in cardiomyocytes. Using both ribosome profiling (Ribo-seq) and a novel orthogonal approach consisting of ribosome pulldown coupled to nanopore sequencing (Nano-TRAP), we find that RPL3L modulates neither translational efficiency nor ribosome affinity towards a specific subset of transcripts. In contrast, we show that depletion of RPL3L leads to increased ribosome-mitochondria interactions in cardiomyocytes, which is accompanied by a significant increase in ATP levels, potentially as a result of fine-tuning of mitochondrial activity. Our results demonstrate that the existence of tissue-specific RP paralogues does not necessarily lead to enhanced translation of specific transcripts or modulation of translational output. Instead, we reveal a complex cellular scenario in which RPL3L modulates the expression of RPL3, which in turn affects ribosomal subcellular localization and, ultimately, mitochondrial activity.

Single-cell quantification of ribosome occupancy in early mouse development

Translation regulation is critical for early mammalian embryonic development¹. However, previous studies had been restricted to bulk measurements², precluding precise determination of translation regulation including allele-specific analyses. Here, to address this challenge, we developed a novel microfluidic isotachophoresis (ITP) approach, named RIBOsome profiling via ITP (Ribo-ITP), and characterized translation in single oocytes and embryos during early mouse development. We identified differential translation efficiency as a key mechanism regulating genes involved in centrosome organization and N⁶-methyladenosine modification of RNAs. Our high-coverage measurements enabled, to our knowledge, the first analysis of allele-specific ribosome engagement in early development. These led to the discovery of stage-specific differential engagement of zygotic RNAs with ribosomes and reduced translation efficiency of transcripts exhibiting allele-biased expression. By integrating our measurements with proteomics data, we discovered that ribosome occupancy in germinal vesicle-stage oocytes is the predominant determinant of protein abundance in the zygote. The Ribo-ITP approach will enable numerous applications by providing high-coverage and high-resolution ribosome occupancy measurements from ultra-low input samples including single cells.

Ribo-uORF: a comprehensive data resource of upstream open reading frames (uORFs) based on ribosome profiling

Upstream open reading frames (uORFs) are typically defined as translation sites located within the 5' untranslated region upstream of the main protein coding sequence (CDS) of messenger RNAs (mRNAs). Although uORFs are prevalent in eukaryotic mRNAs and modulate the translation of downstream CDSs, a comprehensive resource for uORFs is currently lacking. We developed Ribo-uORF (http://rnainformatics.org.cn/RiboUORF) to serve as a comprehensive functional resource for uORF analysis based on ribosome profiling (Ribo-seq) data. Ribo-uORF currently supports six species: human, mouse, rat, zebrafish, fruit fly, and worm. Ribo-uORF includes 501 554 actively translated uORFs and 107 914 upstream translation initiation sites (uTIS), which were identified from 1495 Ribo-seq and 77 quantitative translation initiation sequencing (QTI-seq) datasets, respectively. We also developed mRNAbrowse to visualize items such as uORFs, cis-regulatory elements, genetic variations, eQTLs, GWAS-based associations, RNA modifications, and RNA editing. Ribo-uORF provides a very intuitive web interface for conveniently browsing, searching, and visualizing uORF data. Finally, uORFscan and UTR5var were developed in Ribo-uORF to precisely identify uORFs and analyze the influence of genetic mutations on uORFs using user-uploaded datasets. Ribo-uORF should greatly facilitate studies of uORFs and their roles in mRNA translation and posttranscriptional control of gene expression.

Developing high-affinity decoy receptors to treat multiple myeloma and diffuse large B cell lymphoma

Disease relapse and treatment-induced immunotoxicity pose significant clinical challenges for patients with hematological cancers. Here, we reveal distinctive requirements for neutralizing TNF receptor ligands APRIL and BAFF and their receptor activity in MM and DLBCL, impacting protein translation and production in MM cells and modulating the translation efficiency of the ATM interactor (ATMIN/ACSIZ). Therapeutically, we investigated the use of BCMA decoy receptor (sBCMA-Fc) as an inhibitor of APRIL and BAFF. While wild-type sBCMA-Fc effectively blocked APRIL signaling in MM, it lacked activity in DLBCL due to its weak BAFF binding. To expand the therapeutic utility of sBCMA-Fc, we engineered an affinity-enhanced mutant sBCMA-Fc fusion molecule (sBCMA-Fc V3) 4- and 500-fold stronger in binding to APRIL and BAFF, respectively. The mutant sBCMA-Fc V3 clone significantly enhanced antitumor activity against both MM and DLBCL. Importantly, we also demonstrated an adequate toxicity profile and on-target mechanism of action in nonhuman primate studies.

Genes with 5' terminal oligopyrimidine tracts preferentially escape global suppression of translation by the SARS-CoV-2 Nsp1 protein

Viruses rely on the host translation machinery to synthesize their own proteins. Consequently, they have evolved varied mechanisms to co-opt host translation for their survival. SARS-CoV-2 relies on a nonstructural protein, Nsp1, for shutting down host translation. However, it is currently unknown how viral proteins and host factors critical for viral replication can escape a global shutdown of host translation. Here, using a novel FACS-based assay called MeTAFlow, we report a dose-dependent reduction in both nascent protein synthesis and mRNA abundance in cells expressing Nsp1. We perform RNA-seq and matched ribosome profiling experiments to identify gene-specific changes both at the mRNA expression and translation levels. We discover that a functionally coherent subset of human genes is preferentially translated in the context of Nsp1 expression. These genes include the translation machinery components, RNA binding proteins, and others important for viral pathogenicity. Importantly, we uncovered a remarkable enrichment of 5' terminal oligo-pyrimidine (TOP) tracts among preferentially translated genes. Using reporter assays, we validated that 5' UTRs from TOP transcripts can drive preferential expression in the presence of Nsp1. Finally, we found that LARP1, a key effector protein in the mTOR pathway, may contribute to preferential translation of TOP transcripts in response to Nsp1 expression. Collectively, our study suggests fine-tuning of host gene expression and translation by Nsp1 despite its global repressive effect on host protein synthesis.

Trips-Viz: an environment for the analysis of public and user-generated ribosome profiling data

Trips-Viz (https://trips.ucc.ie/) is an interactive platform for the analysis and visualization of ribosome profiling (Ribo-Seq) and shotgun RNA sequencing (RNA-seq) data. This includes publicly available and user generated data, hence Trips-Viz can be classified as a database and as a server. As a database it provides access to many processed Ribo-Seq and RNA-seq data aligned to reference transcriptomes which has been expanded considerably since its inception. Here, we focus on the server functionality of Trips-viz which also has been greatly improved. Trips-viz now enables visualisation of proteomics data from a large number of processed mass spectrometry datasets. It can be used to support translation inferred from Ribo-Seq data. Users are now able to upload a custom reference transcriptome as well as data types other than Ribo-Seq/RNA-Seq. Incorporating custom data has been streamlined with RiboGalaxy (https://ribogalaxy.ucc.ie/) integration. The other new functionality is the rapid detection of translated open reading frames (ORFs) through a simple easy to use interface. The analysis of differential expression has been also improved via integration of DESeq2 and Anota2seq in addition to a number of other improvements of existing Trips-viz features.

ORFik: a comprehensive R toolkit for the analysis of translation

BACKGROUND

With the rapid growth in the use of high-throughput methods for characterizing translation and the continued expansion of multi-omics, there is a need for back-end functions and streamlined tools for processing, analyzing, and characterizing data produced by these assays.

RESULTS

Here, we introduce ORFik, a user-friendly R/Bioconductor API and toolbox for studying translation and its regulation. It extends GenomicRanges from the genome to the transcriptome and implements a framework that integrates data from several sources. ORFik streamlines the steps to process, analyze, and visualize the different steps of translation with a particular focus on initiation and elongation. It accepts high-throughput sequencing data from ribosome profiling to quantify ribosome elongation or RCP-seq/TCP-seq to also quantify ribosome scanning. In addition, ORFik can use CAGE data to accurately determine 5'UTRs and RNA-seq for determining translation relative to RNA abundance. ORFik supports and calculates over 30 different translation-related features and metrics from the literature and can annotate translated regions such as proteins or upstream open reading frames (uORFs). As a use-case, we demonstrate using ORFik to rapidly annotate the dynamics of 5' UTRs across different tissues, detect their uORFs, and characterize their scanning and translation in the downstream protein-coding regions.

CONCLUSION

In summary, ORFik introduces hundreds of tested, documented and optimized methods. ORFik is designed to be easily customizable, enabling users to create complete workflows from raw data to publication-ready figures for several types of sequencing data. Finally, by improving speed and scope of many core Bioconductor functions, ORFik offers enhancement benefiting the entire Bioconductor environment.

AVAILABILITY

http://bioconductor.org/packages/ORFik .

Genes with 5' terminal oligopyrimidine tracts preferentially escape global suppression of translation by the SARS-CoV-2 Nsp1 protein

Viruses rely on the host translation machinery to synthesize their own proteins. Consequently, they have evolved varied mechanisms to co-opt host translation for their survival. SARS-CoV-2 relies on a non-structural protein, Nsp1, for shutting down host translation. However, it is currently unknown how viral proteins and host factors critical for viral replication can escape a global shutdown of host translation. Here, using a novel FACS-based assay called MeTAFlow, we report a dose-dependent reduction in both nascent protein synthesis and mRNA abundance in cells expressing Nsp1. We perform RNA-Seq and matched ribosome profiling experiments to identify gene-specific changes both at the mRNA expression and translation level. We discover a functionally-coherent subset of human genes are preferentially translated in the context of Nsp1 expression. These genes include the translation machinery components, RNA binding proteins, and others important for viral pathogenicity. Importantly, we uncovered a remarkable enrichment of 5' terminal oligo-pyrimidine (TOP) tracts among preferentially translated genes. Using reporter assays, we validated that 5' UTRs from TOP transcripts can drive preferential expression in the presence of NSP1. Finally, we found that LARP1, a key effector protein in the mTOR pathway may contribute to preferential translation of TOP transcripts in response to Nsp1 expression. Collectively, our study suggests fine tuning of host gene expression and translation by Nsp1 despite its global repressive effect on host protein synthesis.

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.