RiboMiner: a toolset for mining multi-dimensional features of the translatome with ribosome profiling data

Alternative translation initiation produces synaptic organizer proteoforms with distinct localization and functions

While many mRNAs contain more than one translation initiation site (TIS), the functions of most alternative TISs and their corresponding protein isoforms (proteoforms) remain undetermined. Here, we showed that alternative usage of CUG and AUG TISs in neuronal pentraxin receptor (NPR) mRNA produced two proteoforms, of which the ratio was regulated by RNA secondary structure and neuronal activity. Downstream AUG initiation truncated the N-terminal transmembrane domain and produced a secreted NPR proteoform sufficient in promoting synaptic clustering of AMPA-type glutamate receptors. Mutations that altered the ratio of NPR proteoforms reduced AMPA receptors in parvalbumin-positive interneurons and affected learning behaviors in mice. In addition to NPR, upstream AUU-initiated N-terminal extension of C1q-like synaptic organizers anchored these otherwise secreted factors to the membrane. Together, these results uncovered the plasticity of N-terminal signal sequences regulated by alternative TIS usage as a potentially widespread mechanism in diversifying protein localization and functions.

RPFdb v3.0: an enhanced repository for ribosome profiling data and related content

RPFdb (http://www.rpfdb.org or http://sysbio.gzzoc.com/rpfdb/) is a comprehensive repository dedicated to hosting ribosome profiling (Ribo-seq) data and related content. Herein, we present RPFdb v3.0, a significant update featuring expanded data content and improved functionality. Key enhancements include (i) increased data coverage, now encompassing 5018 Ribo-seq datasets and 2343 matched RNA-seq datasets from 496 studies across 34 species; (ii) implementation of translation efficiency, combining Ribo-seq and RNA-seq data to provide gene-specific translation efficiency; (iii) addition of pausing score, facilitating the identification of condition-specific triplet amino acid motifs with enhanced ribosome enrichment; (iv) refinement of open reading frame (ORF) annotation, leveraging RibORF v2.0 for more sensitive detection of actively translated ORFs; (v) introduction of a resource hub, curating advances in translatome sequencing techniques and data analytics tools to support a panoramic overview of the field; and (vi) redesigned web interface, providing intuitive navigation with dedicated pages for streamlined data retrieval, comparison and visualization. These enhancements make RPFdb a more powerful and user-friendly resource for researchers in the field of translatomics. The database is freely accessible and regularly updated to ensure its continued relevance to the scientific community.

SMYD5 is a ribosomal methyltransferase that catalyzes RPL40 lysine methylation to enhance translation output and promote hepatocellular carcinoma

While lysine methylation is well-known for regulating gene expression transcriptionally, its implications in translation have been largely uncharted. Trimethylation at lysine 22 (K22me3) on RPL40, a core ribosomal protein located in the GTPase activation center, was first reported 27 years ago. Yet, its methyltransferase and role in translation remain unexplored. Here, we report that SMYD5 has robust in vitro activity toward RPL40 K22 and primarily catalyzes RPL40 K22me3 in cells. The loss of SMYD5 and RPL40 K22me3 leads to reduced translation output and disturbed elongation as evidenced by increased ribosome collisions. SMYD5 and RPL40 K22me3 are upregulated in hepatocellular carcinoma (HCC) and negatively correlated with patient prognosis. Depleting SMYD5 renders HCC cells hypersensitive to mTOR inhibition in both 2D and 3D cultures. Additionally, the loss of SMYD5 markedly inhibits HCC development and growth in both genetically engineered mouse and patient-derived xenograft (PDX) models, with the inhibitory effect in the PDX model further enhanced by concurrent mTOR suppression. Our findings reveal a novel role of the SMYD5 and RPL40 K22me3 axis in translation elongation and highlight the therapeutic potential of targeting SMYD5 in HCC, particularly with concurrent mTOR inhibition. This work also conceptually broadens the understanding of lysine methylation, extending its significance from transcriptional regulation to translational control.

Widespread changes to the translational landscape in a maize microRNA biogenesis mutant

MicroRNAs are short, non-coding RNAs that repress gene expression in both plants and animals and have diverse functions related to growth, development, and stress responses. The ribonuclease, DICER-LIKE1 (DCL1) is required for two steps in plant miRNA biogenesis: cleavage of the primary miRNAs (pri-miRNAs) to release a hairpin structure, called the precursor miRNA (pre-miRNA) and cleavage of the pre-miRNA to generate the miRNA/miRNA* duplex. The mature miRNA guides the RNA-induced silencing complex to target RNAs with complementary sequences, resulting in translational repression and/or RNA cleavage of target mRNAs. However, the relative contribution of translational repression versus mRNA degradation by miRNAs remains unknown at the genome-level in crops, especially in maize. The maize fuzzy tassel (fzt) mutant contains a hypomorphic mutation in DCL1 resulting in broad developmental defects. While most miRNAs are reduced in fzt, the levels of miRNA-targeted mRNAs are not dramatically increased, suggesting that translational regulation by miRNAs may be common. To gain insight into the repression mechanism of plant miRNAs, we combined ribosome profiling and RNA-sequencing to globally survey miRNA activities in maize. Our data indicate that translational repression contributes significantly to regulation of most miRNA targets and that approximately one-third of miRNA targets are regulated primarily at the translational level. Surprisingly, ribosomes appear altered in fzt mutants suggesting that DCL1 may also have a role in ribosome biogenesis. Thus, DICER-LIKE1 shapes the translational landscape in plants through both miRNA-dependent and miRNA-independent mechanisms.

Next generation sequencing technologies to address aberrant mRNA translation in cancer

In this review, we explore the transformative impact of next generation sequencing technologies in the realm of translatomics (the study of how translational machinery acts on a genome-wide scale). Despite the expectation of a direct correlation between mRNA and protein content, the complex regulatory mechanisms that affect this relationship remark the limitations of standard RNA-seq approaches. Then, the review characterizes crucial techniques such as polysome profiling, ribo-seq, trap-seq, proximity-specific ribosome profiling, rnc-seq, tcp-seq, qti-seq and scRibo-seq. All these methods are summarized within the context of cancer research, shedding light on their applications in deciphering aberrant translation in cancer cells. In addition, we encompass databases and bioinformatic tools essential for researchers that want to address translatome analysis in the context of cancer biology.

MetageneCluster: a Python package for filtering conflicting signal trends in metagene plots

BACKGROUND

Metagene plots provide a visualization of biological signal trends over subsections of the genome and are used to perform high-level analysis of experimental data by aggregating genome-level data to create an average profile. The generation of metagene plots is useful for summarizing the results of many sequencing-based applications. Despite their prevalence and utility, the standard metagene plot is blind to conflicting signals within data. If multiple distinct trends occur, they can interact destructively, creating a plot that does not accurately represent any of the underlying trends.

RESULTS

We present MetageneCluster, a Python tool to generate a collection of representative metagene plots based on k-means clustering of genomic regions of interest. Clustering the data by similarity allows us to identify patterns within the features of interest. We are then able to summarize each pattern present in the data, rather than averaging across the entire feature space. We show that our method performs well when used to identify conflicting signals in real-world genome-level data.

CONCLUSIONS

Overall, MetageneCluster is a user-friendly tool for the creation of metagene plots that capture distinct patterns in underlying sequence data.

Mechanisms of readthrough mitigation reveal principles of GCN1-mediated translational quality control

Readthrough into the 3' untranslated region (3' UTR) of the mRNA results in the production of aberrant proteins. Metazoans efficiently clear readthrough proteins, but the underlying mechanisms remain unknown. Here, we show in Caenorhabditis elegans and mammalian cells that readthrough proteins are targeted by a coupled, two-level quality control pathway involving the BAG6 chaperone complex and the ribosome-collision-sensing protein GCN1. Readthrough proteins with hydrophobic C-terminal extensions (CTEs) are recognized by SGTA-BAG6 and ubiquitylated by RNF126 for proteasomal degradation. Additionally, cotranslational mRNA decay initiated by GCN1 and CCR4/NOT limits the accumulation of readthrough products. Unexpectedly, selective ribosome profiling uncovered a general role of GCN1 in regulating translation dynamics when ribosomes collide at nonoptimal codons, enriched in 3' UTRs, transmembrane proteins, and collagens. GCN1 dysfunction increasingly perturbs these protein classes during aging, resulting in mRNA and proteome imbalance. Our results define GCN1 as a key factor acting during translation in maintaining protein homeostasis.

Pan-cancer analyses of synonymous mutations based on tissue-specific codon optimality

•

Developed tissue-specific codon optimality in 29 human tissues.

•

Applied these to analyze synonymous mutations in ∼10,000 tumor and normal samples.

•

Synonymous mutations frequently increase optimal codons in most cancer types.

•

Synonymous mutations frequently increase optimal codons cell cycle-related genes.

•

Frequency of optimal codon gain relates to proliferation, DDR deficiency, and survival.

Codon optimality has been demonstrated to be an important determinant of mRNA stability and expression levels in multiple model organisms and human cell lines. However, tissue-specific codon optimality has not been developed to investigate how codon optimality is usually perturbed by somatic synonymous mutations in human cancers. Here, we determined tissue-specific codon optimality in 29 human tissues based on mRNA expression data from the Genotype-Tissue Expression project. We found that optimal codons were associated with differentiation, whereas non-optimal codons were correlated with proliferation. Furthermore, codons biased toward differentiation displayed greater tissue specificity in codon optimality, and the tissue specificity of codon optimality was primarily present in amino acids with high degeneracy of the genetic code. By applying tissue-specific codon optimality to somatic synonymous mutations in 8532 tumor samples across 24 cancer types and to those in 416 normal cells across six human tissues, we found that synonymous mutations frequently increased optimal codons in tumor cells and cancer-related genes (e.g., genes involved in cell cycle). Furthermore, an elevated frequency of optimal codon gain was found to promote tumor cell proliferation in three cancer types characterized by DNA damage repair deficiency and could act as a prognostic biomarker for patients with triple-negative breast cancer. In summary, this study profiled tissue-specific codon optimality in human tissues, revealed alterations in codon optimality caused by synonymous mutations in human cancers, and highlighted the non-negligible role of optimal codon gain in tumorigenesis and therapeutics.

Ultrasensitive Ribo-seq reveals translational landscapes during mammalian oocyte-to-embryo transition and pre-implantation development

In mammals, translational control plays critical roles during oocyte-to-embryo transition (OET) when transcription ceases. However, the underlying regulatory mechanisms remain challenging to study. Here, using low-input Ribo-seq (Ribo-lite), we investigated translational landscapes during OET using 30-150 mouse oocytes or embryos per stage. Ribo-lite can also accommodate single oocytes. Combining PAIso-seq to interrogate poly(A) tail lengths, we found a global switch of translatome that closely parallels changes of poly(A) tails upon meiotic resumption. Translation activation correlates with polyadenylation and is supported by polyadenylation signal proximal cytoplasmic polyadenylation elements (papCPEs) in 3' untranslated regions. By contrast, translation repression parallels global de-adenylation. The latter includes transcripts containing no CPEs or non-papCPEs, which encode many transcription regulators that are preferentially re-activated before zygotic genome activation. CCR4-NOT, the major de-adenylation complex, and its key adaptor protein BTG4 regulate translation downregulation often independent of RNA decay. BTG4 is not essential for global de-adenylation but is required for selective gene de-adenylation and production of very short-tailed transcripts. In sum, our data reveal intimate interplays among translation, RNA stability and poly(A) tail length regulation underlying mammalian OET.

Cellular responses to halofuginone reveal a vulnerability of the GCN2 branch of the integrated stress response

Halofuginone (HF) is a phase 2 clinical compound that inhibits the glutamyl-prolyl-tRNA synthetase (EPRS) thereby inducing the integrated stress response (ISR). Here, we report that halofuginone indeed triggers the predicted canonical ISR adaptations, consisting of attenuation of protein synthesis and gene expression reprogramming. However, the former is surprisingly atypical and occurs to a similar magnitude in wild-type cells, cells lacking GCN2 and those incapable of phosphorylating eIF2α. Proline supplementation rescues the observed HF-induced changes indicating that they result from inhibition of EPRS. The failure of the GCN2-to-eIF2α pathway to elicit a measurable protective attenuation of translation initiation allows translation elongation defects to prevail upon HF treatment. Exploiting this vulnerability of the ISR, we show that cancer cells with increased proline dependency are more sensitive to halofuginone. This work reveals that the consequences of EPRS inhibition are more complex than anticipated and provides novel insights into ISR signaling, as well as a molecular framework to guide the targeted development of halofuginone as a therapeutic.

riboviz 2: a flexible and robust ribosome profiling data analysis and visualization workflow

MOTIVATION

Ribosome profiling, or Ribo-seq, is the state-of-the-art method for quantifying protein synthesis in living cells. Computational analysis of Ribo-seq data remains challenging due to the complexity of the procedure, as well as variations introduced for specific organisms or specialized analyses.

RESULTS

We present riboviz 2, an updated riboviz package, for the comprehensive transcript-centric analysis and visualization of Ribo-seq data. riboviz 2 includes an analysis workflow built on the Nextflow workflow management system for end-to-end processing of Ribo-seq data. riboviz 2 has been extensively tested on diverse species and library preparation strategies, including multiplexed samples. riboviz 2 is flexible and uses open, documented file formats, allowing users to integrate new analyses with the pipeline.

AVAILABILITY AND IMPLEMENTATION

riboviz 2 is freely available at github.com/riboviz/riboviz.

HRIBO: high-throughput analysis of bacterial ribosome profiling data

MOTIVATION

Ribosome profiling (Ribo-seq) is a powerful approach based on deep sequencing of cDNA libraries generated from ribosome-protected RNA fragments to explore the translatome of a cell, and is especially useful for the detection of small proteins (50-100 amino acids) that are recalcitrant to many standard biochemical and in silico approaches. While pipelines are available to analyze Ribo-seq data, none are designed explicitly for the automatic processing and analysis of data from bacteria, nor are they focused on the discovery of unannotated open reading frames (ORFs).

RESULTS

We present HRIBO (High-throughput annotation by Ribo-seq), a workflow to enable reproducible and high-throughput analysis of bacterial Ribo-seq data. The workflow performs all required pre-processing and quality control steps. Importantly, HRIBO outputs annotation-independent ORF predictions based on two complementary bacteria-focused tools, and integrates them with additional feature information and expression values. This facilitates the rapid and high-confidence discovery of novel ORFs and their prioritization for functional characterization.

AVAILABILITY AND IMPLEMENTATION

HRIBO is a free and open source project available under the GPL-3 license at: https://github.com/RickGelhausen/HRIBO.

ORFLine: a bioinformatic pipeline to prioritise small open reading frames identifies candidate secreted small proteins from lymphocytes

MOTIVATION

The annotation of small open reading frames (smORFs) of less than 100 codons (<300 nucleotides) is challenging due to the large number of such sequences in the genome.

RESULTS

In this study, we developed a computational pipeline, which we have named ORFLine, that stringently identifies smORFs and classifies them according to their position within transcripts. We identified a total of 5744 unique smORFs in datasets from mouse B and T lymphocytes and systematically characterized them using ORFLine. We further searched smORFs for the presence of a signal peptide, which predicted known secreted chemokines as well as novel micropeptides. Four novel micropeptides show evidence of secretion and are therefore candidate mediators of immunoregulatory functions.

AVAILABILITY

Freely available on the web at https://github.com/boboppie/ORFLine.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.