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

RiboParser/RiboShiny: an integrated platform for comprehensive analysis and visualization of Ribo-seq data

Translation is a crucial step in gene expression. Over the past decade, the development and application of Ribosome profiling (Ribo-seq) have significantly advanced our understanding of translational regulation in vivo. However, the analysis and visualization of Ribo-seq data remain challenging. Despite the availability of various analytical pipelines, improvements in comprehensiveness, accuracy, and user-friendliness are still necessary. In this study, we develop RiboParser/RiboShiny, a robust framework for analyzing and visualizing Ribo-seq data. Building on published methods, we optimize ribosome structure-based and start/stop-based models to improve the accuracy and stability of P-site detection, even in species with a high proportion of leaderless transcripts. Leveraging these improvements, RiboParser offers comprehensive analyses, including quality control, gene-level analysis, codon-level analysis, and the analysis of Ribo-seq variants. Meanwhile, RiboShiny provides a user-friendly and adaptable platform for data visualization, facilitating deeper insights into the translational landscape. Furthermore, the integration of standardized genome annotation renders our platform universally applicable to various organisms with sequenced genomes. This framework has the potential to significantly improve the precision and efficiency of Ribo-seq data interpretation, thereby deepening our understanding of translational regulation.

Small open reading frame-encoded microproteins in cancer: identification, biological functions and clinical significance

The human genome harbors approximately twenty thousand protein-coding genes, and a significant portion of life science research focuses on elucidating their functions and the underlying mechanisms. Recent studies have revealed that small open reading frame (sORF), originating from non-coding RNAs or the 5' leader sequences of messenger RNAs, can be translated into small peptides called microproteins through cap-dependent or cap-independent mechanisms. These microproteins interact with diverse molecular partners to modulate gene expression at multiple regulatory levels, thereby playing critical roles in various biological processes. Notably, sORF-encoded microproteins exhibit aberrant expression patterns in cancer and are implicated in tumor initiation and progression, expanding our understanding of cancer biology. In this review, we introduce the translational mechanisms and identification methods of microproteins, summarize their dysregulation in cancer and their biological functions in regulating gene expression, and emphasize their roles in driving hallmark events of cancer. Furthermore, we discuss their clinical significance as diagnostic and prognostic biomarkers, as well as therapeutic targets.

Extensive location bias of the GPCR-dependent translatome via site-selective activation of mTOR

G protein-coupled receptors (GPCRs) modulate various physiological functions by rewiring cellular gene expression in response to extracellular signals. Control of gene expression by GPCRs has been studied almost exclusively at the transcriptional level, neglecting an extensive amount of regulation that takes place translationally. Hence, little is known about the nature and mechanisms of gene-specific posttranscriptional regulation downstream of receptor activation. Here, we apply an unbiased multiomics approach to delineate an extensive translational regulatory program initiated by the prototypical beta2-adrenergic receptor (β2-AR) and provide mechanistic insights into how these processes are orchestrated. Using ribosome profiling (Ribo-seq), we identify nearly 120 gene targets of adrenergic receptor activity for which expression is exclusively regulated at the level of translation. We next show that all translational changes are induced selectively by endosomal β2-ARs and report that this proceeds through activation of the mammalian target of rapamycin (mTOR) pathway. Specifically, within the set of translational GPCR targets, we find significant enrichment of genes with 5' terminal oligopyrimidine (TOP) motifs, a gene class classically known to be translationally regulated by mTOR. We then demonstrate that endosomal β2-ARs are required for mTOR activation and subsequent mTOR-dependent TOP mRNA translation. This site-selective crosstalk between the pathways is observed in multiple cell models with native β2-ARs, across a range of endogenous and synthetic adrenergic agonists, and for other GPCRs with intracellular activity. Together, this comprehensive analysis of drug-induced translational regulation establishes a critical role for location-biased GPCR signaling in fine-tuning the cellular protein landscape.

Decoding Codon Bias: The Role of tRNA Modifications in Tissue-Specific Translation

The tRNA epitranscriptome has been recognized as an important player in mRNA translation regulation. Our knowledge of the role of the tRNA epitranscriptome in fine-tuning translation via codon decoding at tissue or cell levels remains incomplete. We analyzed tRNA expression and modifications as well as codon optimality across seven mouse tissues. Our analysis revealed distinct enrichment patterns of tRNA modifications in different tissues. Queuosine (Q) tRNA modification was most enriched in the brain compared to other tissues, while mitochondrial tRNA modifications and tRNA expression were highest in the heart. Using this observation, we synthesized, and delivered in vivo, codon-mutated EGFP for Q-codons, where the C-ending Q-codons were replaced with U-ending codons. The protein levels of mutant EGFP were downregulated in liver, which is poor in Q, while in brain EGFP, levels did not change. These data show that understanding tRNA modification enrichments across tissues is not only essential for understanding codon decoding and bias but can also be utilized for optimizing gene and mRNA therapeutics to be more tissue-, cell-, or condition-specific.

The transcriptional and translational landscape of HCoV-OC43 infection

The coronavirus HCoV-OC43 circulates continuously in the human population and is a frequent cause of the common cold. Here, we generated a high-resolution atlas of the transcriptional and translational landscape of OC43 during a time course following infection of human lung fibroblasts. Using ribosome profiling, we quantified the relative expression of the canonical open reading frames (ORFs) and identified previously unannotated ORFs. These included several potential short upstream ORFs and a putative ORF nested inside the M gene. In parallel, we analyzed the cellular response to infection. Endoplasmic reticulum (ER) stress response genes were transcriptionally and translationally induced beginning 12 and 18 hours post infection, respectively. By contrast, conventional antiviral genes mostly remained quiescent. At the same time points, we observed accumulation and increased translation of noncoding transcripts normally targeted by nonsense mediated decay (NMD), suggesting NMD is suppressed during the course of infection. This work provides resources for deeper understanding of OC43 gene expression and the cellular responses during infection.

Phosphorylation of P-stalk proteins defines the ribosomal state for interaction with auxiliary protein factors

Ribosomal action is facilitated by the orchestrated work of trans-acting factors and ribosomal elements, which are subject to regulatory events, often involving phosphorylation. One such element is the ribosomal P-stalk, which plays a dual function: it activates translational GTPases, which support basic ribosomal functions, and interacts with the Gcn2 kinase, linking the ribosomes to the ISR pathway. We show that P-stalk proteins, which form a pentamer, exist in the cell exclusively in a phosphorylated state at five C-terminal domains (CTDs), ensuring optimal translation (speed and accuracy) and may play a role in the timely regulation of the Gcn2-dependent stress response. Phosphorylation of the CTD induces a structural transition from a collapsed to a coil-like structure, and the CTD gains conformational freedom, allowing specific but transient binding to various protein partners, optimizing the ribosome action. The report reveals a unique feature of the P-stalk proteins, indicating that, unlike most ribosomal proteins, which are regulated by phosphorylation in an on/off manner, the P-stalk proteins exist in a constantly phosphorylated state, which optimizes their interaction with auxiliary factors.

Transposable Elements Contribute to the Regulation of Long Noncoding RNAs in Drosophila melanogaster

Background: Transposable elements (TEs) and noncoding sequences are major components of the genome, yet their functional contributions to long noncoding RNAs (lncRNAs) are not well understood. Although many lncRNAs originating from TEs (TE-lncRNAs) have been identified across various organisms, their characteristics and regulatory roles, particularly in insects, remain largely unexplored. This study integrated multi-omics data to investigate TE-lncRNAs in D. melanogaster, focusing on the influence of transposons across different omics levels. Results: We identified 16,118 transposons overlapping with lncRNA sequences that constitute 2119 TE-lncRNAs (40.4% of all lncRNAs) using 256 public RNA-seq samples and 15 lncRNA-seq samples of Drosophila S2 cells treated with heavy metals. Of these, 67.2% of TE-lncRNAs contain more than one TE. The LTR/Gypsy family was the most common transposon insertion. Transposons preferred to insert into promoters, transcription starting sites, and intronic regions, especially in chromosome ends. Compared with lncRNAs, TE-lncRNAs showed longer lengths, a lower conservation, and lower levels but a higher specificity of expression. Multi-omics data analysis revealed positive correlations between transposon insertions and chromatin openness at the pre-transcriptional level. Notably, a total of 516 TE-lncRNAs provided transcriptional factor binding sites through transposon insertions. The regulatory network of a key transcription factor was rewired by transposons, potentially recruiting other transcription factors to exert regulatory functions under heavy metal stress. Additionally, 99 TE-lncRNAs were associated with m6A methylation modification sites, and 115 TE-lncRNAs potentially provided candidate small open reading frames through transposon insertions. Conclusions: Our data analysis demonstrated that TEs contribute to the regulation of lncRNAs. TEs not only promote the transcriptional regulation of lncRNAs, but also facilitate their post-transcriptional and epigenetic regulation.

Exploring the potential of large language model-based chatbots in challenges of ribosome profiling data analysis: a review

Ribosome profiling (Ribo-seq) provides transcriptome-wide insights into protein synthesis dynamics, yet its analysis poses challenges, particularly for nonbioinformatics researchers. Large language model-based chatbots offer promising solutions by leveraging natural language processing. This review explores their convergence, highlighting opportunities for synergy. We discuss challenges in Ribo-seq analysis and how chatbots mitigate them, facilitating scientific discovery. Through case studies, we illustrate chatbots' potential contributions, including data analysis and result interpretation. Despite the absence of applied examples, existing software underscores the value of chatbots and the large language model. We anticipate their pivotal role in future Ribo-seq analysis, overcoming limitations. Challenges such as model bias and data privacy require attention, but emerging trends offer promise. The integration of large language models and Ribo-seq analysis holds immense potential for advancing translational regulation and gene expression understanding.

Disruption of tRNA biogenesis enhances proteostatic resilience, improves later-life health, and promotes longevity

tRNAs are evolutionarily ancient molecular decoders essential for protein translation. In eukaryotes, tRNAs and other short, noncoding RNAs are transcribed by RNA polymerase (Pol) III, an enzyme that promotes ageing in yeast, worms, and flies. Here, we show that a partial reduction in Pol III activity specifically disrupts tRNA levels. This effect is conserved across worms, flies, and mice, where computational models indicate that it impacts mRNA decoding. In all 3 species, reduced Pol III activity increases proteostatic resilience. In worms, it activates the unfolded protein response (UPR) and direct disruption of tRNA metabolism is sufficient to recapitulate this. In flies, decreasing Pol III's transcriptional initiation on tRNA genes by a loss-of-function in the TFIIIC transcription factor robustly extends lifespan, improves proteostatic resilience and recapitulates the broad-spectrum benefits to late-life health seen following partial Pol III inhibition. We provide evidence that a partial reduction in Pol III activity impacts translation, quantitatively or qualitatively, in both worms and flies, indicating a potential mode of action. Our work demonstrates a conserved and previously unappreciated role of tRNAs in animal ageing.

Metabolic rewiring during bone development underlies tRNA m7G-associated primordial dwarfism

Translation of mRNA to protein is tightly regulated by transfer RNAs (tRNAs), which are subject to various chemical modifications that maintain structure, stability, and function. Deficiency of tRNA N7-methylguanosine (m7G) modification in patients causes a type of primordial dwarfism, but the underlying mechanism remains unknown. Here we report that the loss of m7G rewires cellular metabolism, leading to the pathogenesis of primordial dwarfism. Conditional deletion of the catalytic enzyme Mettl1 or missense mutation of the scaffold protein Wdr4 severely impaired endochondral bone formation and bone mass accrual. Mechanistically, Mettl1 knockout decreased abundance of m7G-modified tRNAs and inhibited translation of mRNAs relating to cytoskeleton and Rho GTPase signaling. Meanwhile, Mettl1 knockout enhanced cellular energy metabolism despite incompetent proliferation and osteogenic commitment. Further exploration revealed that impairment of Rho GTPase signaling upregulated the level of branched-chain amino acid transaminase 1 (BCAT1) that rewired cell metabolism and restricted intracellular α-ketoglutarate (αKG). Supplementation of αKG ameliorated the skeletal defect of Mettl1-deficient mice. In addition to the selective translation of metabolism-related mRNAs, we further revealed that Mettl1 knockout globally regulated translation via integrated stress response (ISR) and mammalian target of rapamycin complex 1 (mTORC1) signaling. Restoring translation by targeting either ISR or mTORC1 aggravated bone defects of Mettl1-deficient mice. Overall, our study unveils a critical role of m7G tRNA modification in bone development by regulation of cellular metabolism and indicates suspension of translation initiation as a quality control mechanism in response to tRNA dysregulation.

The tRNA Gm18 methyltransferase TARBP1 promotes hepatocellular carcinoma progression via metabolic reprogramming of glutamine

Cancer cells rely on metabolic reprogramming to sustain the prodigious energetic requirements for rapid growth and proliferation. Glutamine metabolism is frequently dysregulated in cancers and is being exploited as a potential therapeutic target. Using CRISPR/Cas9 interference (CRISPRi) screening, we identified TARBP1 (TAR (HIV-1) RNA Binding Protein 1) as a critical regulator involved in glutamine reliance of cancer cell. Consistent with this discovery, TARBP1 amplification and overexpression are frequently observed in various cancers. Knockout of TARBP1 significantly suppresses cell proliferation, colony formation and xenograft tumor growth. Mechanistically, TARBP1 selectively methylates and stabilizes a small subset of tRNAs, which promotes efficient protein synthesis of glutamine transporter-ASCT2 (also known as SLC1A5) and glutamine import to fuel the growth of cancer cell. Moreover, we found that the gene expression of TARBP1 and ASCT2 are upregulated in combination in clinical cohorts and their upregulation is associated with unfavorable prognosis of HCC (hepatocellular carcinoma). Taken together, this study reveals the unexpected role of TARBP1 in coordinating the tRNA availability and glutamine uptake during HCC progression and provides a potential target for tumor therapy.

m3C32 tRNA modification controls serine codon-biased mRNA translation, cell cycle, and DNA-damage response

The epitranscriptome includes a diversity of RNA modifications that influence gene expression. N3-methylcytidine (m3C) mainly occurs in the anticodon loop (position C32) of certain tRNAs yet its role is poorly understood. Here, using HAC-Seq, we report comprehensive METTL2A/2B-, METTL6-, and METTL2A/2B/6-dependent m3C profiles in human cells. METTL2A/2B modifies tRNA-arginine and tRNA-threonine members, whereas METTL6 modifies the tRNA-serine family. However, decreased m3C32 on tRNA-Ser-GCT isodecoders is only observed with combined METTL2A/2B/6 deletion. Ribo-Seq reveals altered translation of genes related to cell cycle and DNA repair pathways in METTL2A/2B/6-deficient cells, and these mRNAs are enriched in AGU codons that require tRNA-Ser-GCT for translation. These results, supported by reporter assays, help explain the observed altered cell cycle, slowed proliferation, and increased cisplatin sensitivity phenotypes of METTL2A/2B/6-deficient cells. Thus, we define METTL2A/2B/6-dependent methylomes and uncover a particular requirement of m3C32 tRNA modification for serine codon-biased mRNA translation of cell cycle, and DNA repair genes.

Extensive location bias of the GPCR-dependent translatome via site-selective activation of mTOR

G protein-coupled receptors (GPCRs) modulate various physiological functions by re-wiring cellular gene expression in response to extracellular signals. Control of gene expression by GPCRs has been studied almost exclusively at the transcriptional level, neglecting an extensive amount of regulation that takes place translationally. Hence, little is known about the nature and mechanisms of gene-specific post-transcriptional regulation downstream of receptor activation. Here, we apply an unbiased multiomics approach to delineate an extensive translational regulatory program initiated by the prototypical beta2-adrenergic receptor (β2-AR) and provide mechanistic insights into how these processes are orchestrated. Using ribosome profiling (Ribo-seq), we identify nearly 120 novel gene targets of adrenergic receptor activity which expression is exclusively regulated at the level of translation. We next show that all translational changes are induced selectively by endosomal β2-ARs. We further report that this proceeds through activation of the mammalian target of rapamycin (mTOR) pathway. Specifically, within the set of translational GPCR targets we discover significant enrichment of genes with 5' terminal oligopyrimidine (TOP) motifs, a gene class classically known to be translationally regulated by mTOR. We then demonstrate that endosomal β2-ARs are required for mTOR activation and subsequent mTOR-dependent TOP mRNA translation. Together, this comprehensive analysis of drug-induced translational regulation establishes a critical role for location-biased GPCR signaling in fine-tuning the cellular protein landscape.

Comparative transcriptome analysis reveals molecular damage associated with cryopreservation in Crassostrea angulata D-larvae rather than to cryoprotectant exposure

BACKGROUND

The Portuguese oyster Crassostrea angulata, a bivalve of significant economic and ecological importance, has faced a decline in both production and natural populations due to pathologies, climate change, and anthropogenic factors. To safeguard its genetic diversity and improve reproductive management, cryopreservation emerges as a valuable strategy. However, the cryopreservation methodologies lead to some damage in structures and functions of the cells and tissues that can affect post-thaw quality. Transcriptomics may help to understand the molecular consequences related to cryopreservation steps and therefore to identify different freezability biomarkers. This study investigates the molecular damage induced by cryopreservation in C. angulata D-larvae, focusing on two critical steps: exposure to cryoprotectant solution and the freezing/thawing process.

RESULTS

Expression analysis revealed 3 differentially expressed genes between larvae exposed to cryoprotectant solution and fresh larvae and 611 differentially expressed genes in cryopreserved larvae against fresh larvae. The most significantly enriched gene ontology terms were "carbohydrate metabolic process", "integral component of membrane" and "chitin binding" for biological processes, cellular components and molecular functions, respectively. Kyoto Encyclopedia of Genes and Genomes enrichment analysis identified the "neuroactive ligand receptor interaction", "endocytosis" and "spliceosome" as the most enriched pathways. RNA sequencing results were validate by quantitative RT-PCR, once both techniques presented the same gene expression tendency and a group of 11 genes were considered important molecular biomarkers to be used in further studies for the evaluation of cryodamage.

CONCLUSIONS

The current work provided valuable insights into the molecular repercussions of cryopreservation on D-larvae of Crassostrea angulata, revealing that the freezing process had a more pronounced impact on larval quality compared to any potential cryoprotectant-induced toxicity. Additionally, was identify 11 genes serving as biomarkers of freezability for D-larvae quality assessment. This research contributes to the development of more effective cryopreservation protocols and detection methods for cryodamage in this species.

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.

Metabolic reprogramming driven by METTL1-mediated tRNA m7G modification promotes acquired anlotinib resistance in oral squamous cell carcinoma

Tyrosine kinase inhibitors (TKIs) are frequently utilized in the management of malignant tumors. Studies have indicated that anlotinib has a significant inhibitory effect on oral squamous cell carcinoma (OSCC). However, the mechanisms underlying the development of resistance with long-term anlotinib treatment remain obscure. Our research found that METTL1 expression was heightened in anlotinib-resistant OSCC cells. We observed that METTL1 played a role in fostering resistance to anlotinib in both transgenic mouse models and in vitro. Mechanistically, the elevated METTL1 levels in anlotinib-resistant OSCC cells contributed to enhanced global mRNA translation and stimulated oxidative phosphorylation (OXPHOS) through m7G tRNA modification. Bioenergetic profiling demonstrated that METTL1 drived a metabolic shift from glycolysis to OXPHOS in anlotinib-resistant OSCC cells. Additionally, inhibition of OXPHOS biochemically negated METTL1's impact on anlotinib resistance. Overall, this study underscores the pivotal role of METTL1-mediated m7G tRNA modification in anlotinib resistance and lays the groundwork for novel therapeutic interventions to counteract resistance in OSCC.

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.

Translational response to mitochondrial stresses is orchestrated by tRNA modifications

Mitochondrial stress and dysfunction play important roles in many pathologies. However, how cells respond to mitochondrial stress is not fully understood. Here, we examined the translational response to electron transport chain (ETC) inhibition and arsenite induced mitochondrial stresses. Our analysis revealed that during mitochondrial stress, tRNA modifications (namely f5C, hm5C, queuosine and its derivatives, and mcm5U) dynamically change to fine tune codon decoding, usage, and optimality. These changes in codon optimality drive the translation of many pathways and gene sets, such as the ATF4 pathway and selenoproteins, involved in the cellular response to mitochondrial stress. We further examined several of these modifications using targeted approaches. ALKBH1 knockout (KO) abrogated f5C and hm5C levels and led to mitochondrial dysfunction, reduced proliferation, and impacted mRNA translation rates. Our analysis revealed that tRNA queuosine (tRNA-Q) is a master regulator of the mitochondrial stress response. KO of QTRT1 or QTRT2, the enzymes responsible for tRNA-Q synthesis, led to mitochondrial dysfunction, translational dysregulation, and metabolic alterations in mitochondria-related pathways, without altering cellular proliferation. In addition, our analysis revealed that tRNA-Q loss led to a domino effect on various tRNA modifications. Some of these changes could be explained by metabolic profiling. Our analysis also revealed that utilizing serum deprivation or alteration with Queuine supplementation to study tRNA-Q or stress response can introduce various confounding factors by altering many other tRNA modifications. In summary, our data show that tRNA modifications are master regulators of the mitochondrial stress response by driving changes in codon decoding.

Transcriptional programming of translation by BCL6 controls skeletal muscle proteostasis

Skeletal muscle is dynamically controlled by the balance of protein synthesis and degradation. Here we discover an unexpected function for the transcriptional repressor B cell lymphoma 6 (BCL6) in muscle proteostasis and strength in mice. Skeletal muscle-specific Bcl6 ablation in utero or in adult mice results in over 30% decreased muscle mass and force production due to reduced protein synthesis and increased autophagy, while it promotes a shift to a slower myosin heavy chain fibre profile. Ribosome profiling reveals reduced overall translation efficiency in Bcl6-ablated muscles. Mechanistically, tandem chromatin immunoprecipitation, transcriptomic and translational analyses identify direct BCL6 repression of eukaryotic translation initiation factor 4E-binding protein 1 (Eif4ebp1) and activation of insulin-like growth factor 1 (Igf1) and androgen receptor (Ar). Together, these results uncover a bifunctional role for BCL6 in the transcriptional and translational control of muscle proteostasis.

Ribosome Profiling of Plants

Translation is a key step in control of gene expression, yet most analyses of global responses to a stimulus focus on transcription and the transcriptome. For RNA viruses in particular, which have no DNA-templated transcriptional control, control of viral and host translation is crucial. Here, we describe the method of ribosome profiling (ribo-seq) in plants, applied to virus infection. Ribo-seq is a deep sequencing technique that reveals the translatome by presenting a snapshot of the positions and relative amounts of translating ribosomes on all mRNAs in the cell. In contrast to RNA-seq, a crude cell extract is first digested with ribonuclease to degrade all mRNA not protected by a translating 80S ribosome. The resulting ribosome-protected fragments (RPFs) are deep sequenced. The number of reads mapping to a specific mRNA compared to the standard RNA-seq reads reveals the translational efficiency of that mRNA. Moreover, the precise positions of ribosome pause sites, previously unknown translatable open reading frames, and noncanonical translation events can be characterized quantitatively using ribo-seq. As this technique requires meticulous technique, here we present detailed step-by-step instructions for cell lysate preparation by flash freezing of samples, nuclease digestion of cell lysate, monosome collection by sucrose cushion ultracentrifugation, size-selective RNA extraction and rRNA depletion, library preparation for sequencing and finally quality control of sequenced data. These experimental methods apply to many plant systems, with minor nuclease digestion modifications depending on the plant tissue and species. This protocol should be valuable for studies of plant virus gene expression, and the global translational response to virus infection, or any other biotic or abiotic stress, by the host plant.

Ribosome biogenesis controls cranial suture MSC fate via the complement pathway in mouse and human iPSC models

Disruption of global ribosome biogenesis selectively affects craniofacial tissues with unclear mechanisms. Craniosynostosis is a congenital craniofacial disorder characterized by premature fusion of cranial suture(s) with loss of suture mesenchymal stem cells (MSCs). Here we focused on ribosomopathy disease gene Snord118, which encodes a small nucleolar RNA (snoRNA), to genetically disturb ribosome biogenesis in suture MSCs using mouse and human induced pluripotent stem cell (iPSC) models. Snord118 depletion exhibited p53 activation, increased cell death, reduced proliferation, and premature osteogenic differentiation of MSCs, leading to suture growth and craniosynostosis defects. Mechanistically, Snord118 deficiency causes translational dysregulation of ribosomal proteins and downregulation of complement pathway genes. Further complement pathway disruption by knockout of complement C3a receptor 1 (C3ar1) exacerbated MSC and suture defects in mutant mice, whereas activating the complement pathway rescued MSC cell fate and suture growth defects. Thus, ribosome biogenesis controls MSC fate via the complement pathway to prevent craniosynostosis.

A Zpr1 co-chaperone mediates folding of eukaryotic translation elongation factor 1A via a GTPase cycle

General protein folding is mediated by chaperones that utilize ATP hydrolysis to regulate client binding and release. Zinc-finger protein 1 (Zpr1) is an essential ATP-independent chaperone dedicated to the biogenesis of eukaryotic translation elongation factor 1A (eEF1A), a highly abundant GTP-binding protein. How Zpr1-mediated folding is regulated to ensure rapid Zpr1 recycling remains an unanswered question. Here, we use yeast genetics and microscopy analysis, biochemical reconstitution, and structural modeling to reveal that folding of eEF1A by Zpr1 requires GTP hydrolysis. Furthermore, we identify the highly conserved altered inheritance of mitochondria 29 (Aim29) protein as a Zpr1 co-chaperone that recognizes eEF1A in the GTP-bound, pre-hydrolysis conformation. This interaction dampens Zpr1⋅eEF1A GTPase activity and facilitates client exit from the folding cycle. Our work reveals that a bespoke ATP-independent chaperone system has mechanistic similarity to ATPase chaperones but unexpectedly relies on client GTP hydrolysis to regulate the chaperone-client interaction.

Small Open Reading Frame-Encoded Micro-Peptides: An Emerging Protein World

Small open reading frames (sORFs) are often overlooked features in genomes. In the past, they were labeled as noncoding or "transcriptional noise". However, accumulating evidence from recent years suggests that sORFs may be transcribed and translated to produce sORF-encoded polypeptides (SEPs) with less than 100 amino acids. The vigorous development of computational algorithms, ribosome profiling, and peptidome has facilitated the prediction and identification of many new SEPs. These SEPs were revealed to be involved in a wide range of basic biological processes, such as gene expression regulation, embryonic development, cellular metabolism, inflammation, and even carcinogenesis. To effectively understand the potential biological functions of SEPs, we discuss the history and development of the newly emerging research on sORFs and SEPs. In particular, we review a range of recently discovered bioinformatics tools for identifying, predicting, and validating SEPs as well as a variety of biochemical experiments for characterizing SEP functions. Lastly, this review underlines the challenges and future directions in identifying and validating sORFs and their encoded micropeptides, providing a significant reference for upcoming research on sORF-encoded peptides.

Variational inference for detecting differential translation in ribosome profiling studies

Translational efficiency change is an important mechanism for regulating protein synthesis. Experiments with paired ribosome profiling (Ribo-seq) and mRNA-sequencing (RNA-seq) allow the study of translational efficiency by simultaneously quantifying the abundances of total transcripts and those that are being actively translated. Existing methods for Ribo-seq data analysis either ignore the pairing structure in the experimental design or treat the paired samples as fixed effects instead of random effects. To address these issues, we propose a hierarchical Bayesian generalized linear mixed effects model which incorporates a random effect for the paired samples according to the experimental design. We provide an analytical software tool, "riboVI," that uses a novel variational Bayesian algorithm to fit our model in an efficient way. Simulation studies demonstrate that "riboVI" outperforms existing methods in terms of both ranking differentially translated genes and controlling false discovery rate. We also analyzed data from a real ribosome profiling experiment, which provided new biological insight into virus-host interactions by revealing changes in hormone signaling and regulation of signal transduction not detected by other Ribo-seq data analysis tools.

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.

Comparative oxidation proteomics analyses suggest redox regulation of cytosolic translation in rice leaves upon Magnaporthe oryzae infection

Pathogen attack can increase plant levels of reactive oxygen species (ROS), which act as signaling molecules to activate plant defense mechanisms. Elucidating these processes is crucial for understanding redox signaling pathways in plant defense responses. Using an iodo-tandem mass tag (TMT)-based quantitative proteomics approach, we mapped 3362 oxidized cysteine sites in 2275 proteins in rice leaves. Oxidized proteins were involved in gene expression, peptide biosynthetic processes, stress responses, ROS metabolic processes, and translation pathways. Magnaporthe oryzae infection led to increased oxidative modification levels of 512 cysteine sites in 438 proteins, including many transcriptional regulators and ribosomal proteins. Ribosome profiling (Ribo-seq) analysis revealed that the oxidative modification of ribosomal proteins promoted the translational efficiency of many mRNAs involved in defense response pathways, thereby affecting rice immunity. Our results suggest that increased oxidative modification of ribosomal proteins in rice leaves promotes cytosolic translation, thus revealing a novel function of post-translational modifications. Furthermore, the oxidation-sensitive proteins identified here provide a valuable resource for research on protein redox regulation and can guide future mechanistic studies.

Effective extraction of polyribosomes exposes gene expression strategies in primary astrocytes

Regulation of mRNA translation in astrocytes gains a growing interest. However, until now, successful ribosome profiling of primary astrocytes has not been reported. Here, we optimized the standard 'polysome profiling' method and generated an effective protocol for polyribosome extraction, which enabled genome-wide assessment of mRNA translation dynamics along the process of astrocyte activation. Transcriptome (RNAseq) and translatome (Riboseq) data generated at 0, 24 and 48 h after cytokines treatment, revealed dynamic genome-wide changes in the expression level of ∼12 000 genes. The data clarify whether a change in protein synthesis rate results from a change in mRNA level or translation efficiency per se. It exhibit different expression strategies, based on changes in mRNA abundance and/or translation efficiency, which are specifically assigned to gene subsets depending on their function. Moreover, the study raises an important take-home message related to the possible presence of 'difficult to extract' polyribosome sub-groups, in all cell types, thus illuminating the impact of ribosomes extraction methodology on experiments addressing translation regulation.

Velcrin-induced selective cleavage of tRNALeu(TAA) by SLFN12 causes cancer cell death

Velcrin compounds kill cancer cells expressing high levels of phosphodiesterase 3A (PDE3A) and Schlafen family member 12 (SLFN12) by inducing complex formation between these two proteins, but the mechanism of cancer cell killing by the PDE3A-SLFN12 complex is not fully understood. Here, we report that the physiological substrate of SLFN12 RNase is tRNALeu(TAA). SLFN12 selectively digests tRNALeu(TAA), and velcrin treatment promotes the cleavage of tRNALeu(TAA) by inducing PDE3A-SLFN12 complex formation in vitro. We found that distinct sequences in the variable loop and acceptor stem of tRNALeu(TAA) are required for substrate digestion. Velcrin treatment of sensitive cells results in downregulation of tRNALeu(TAA), ribosome pausing at Leu-TTA codons and global inhibition of protein synthesis. Velcrin-induced cleavage of tRNALeu(TAA) by SLFN12 and the concomitant global inhibition of protein synthesis thus define a new mechanism of apoptosis initiation.

A Dynamic rRNA Ribomethylome Drives Stemness in Acute Myeloid Leukemia

The development and regulation of malignant self-renewal remain unresolved issues. Here, we provide biochemical, genetic, and functional evidence that dynamics in ribosomal RNA (rRNA) 2'-O-methylation regulate leukemia stem cell (LSC) activity in vivo. A comprehensive analysis of the rRNA 2'-O-methylation landscape of 94 patients with acute myeloid leukemia (AML) revealed dynamic 2'-O-methylation specifically at exterior sites of ribosomes. The rRNA 2'-O-methylation pattern is closely associated with AML development stage and LSC gene expression signature. Forced expression of the 2'-O-methyltransferase fibrillarin (FBL) induced an AML stem cell phenotype and enabled engraftment of non-LSC leukemia cells in NSG mice. Enhanced 2'-O-methylation redirected the ribosome translation program toward amino acid transporter mRNAs enriched in optimal codons and subsequently increased intracellular amino acid levels. Methylation at the single site 18S-guanosine 1447 was instrumental for LSC activity. Collectively, our work demonstrates that dynamic 2'-O-methylation at specific sites on rRNAs shifts translational preferences and controls AML LSC self-renewal.

SIGNIFICANCE

We establish the complete rRNA 2'-O-methylation landscape in human AML. Plasticity of rRNA 2'-O-methylation shifts protein translation toward an LSC phenotype. This dynamic process constitutes a novel concept of how cancers reprogram cell fate and function. This article is highlighted in the In This Issue feature, p. 247.

A novel active transposon creates allelic variation through altered translation rate to influence protein abundance

Protein translation is tightly and precisely controlled by multiple mechanisms including upstream open reading frames (uORFs), but the origins of uORFs and their role in maize are largely unexplored. In this study, an active transposition event was identified during the propagation of maize inbred line B73. The transposon, which was named BTA for 'B73 active transposable element hAT', creates a novel dosage-dependent hypomorphic allele of the hexose transporter gene ZmSWEET4c through insertion within the coding sequence in the first exon, and results in reduced kernel size. The BTA insertion does not affect transcript abundance but reduces protein abundance of ZmSWEET4c, probably through the introduction of a uORF. Furthermore, the introduction of BTA sequence in the exon of other genes can regulate translation efficiency without affecting their mRNA levels. A transposon capture assay revealed 79 novel insertions for BTA and BTA-like elements. These insertion sites have typical euchromatin features, including low levels of DNA methylation and high levels of H3K27ac. A putative autonomous element that mobilizes BTA and BTA-like elements was identified. Together, our results suggest a transposon-based origin of uORFs and document a new role for transposable elements to influence protein abundance and phenotypic diversity by affecting the translation rate.

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.

Discovering microproteins: making the most of ribosome profiling data

Building a reference set of protein-coding open reading frames (ORFs) has revolutionized biological process discovery and understanding. Traditionally, gene models have been confirmed using cDNA sequencing and encoded translated regions inferred using sequence-based detection of start and stop combinations longer than 100 amino-acids to prevent false positives. This has led to small ORFs (smORFs) and their encoded proteins left un-annotated. Ribo-seq allows deciphering translated regions from untranslated irrespective of the length. In this review, we describe the power of Ribo-seq data in detection of smORFs while discussing the major challenge posed by data-quality, -depth and -sparseness in identifying the start and end of smORF translation. In particular, we outline smORF cataloguing efforts in humans and the large differences that have arisen due to variation in data, methods and assumptions. Although current versions of smORF reference sets can already be used as a powerful tool for hypothesis generation, we recommend that future editions should consider these data limitations and adopt unified processing for the community to establish a canonical catalogue of translated smORFs.

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.

Short open reading frame genes in innate immunity: from discovery to characterization

Next-generation sequencing (NGS) technologies have greatly expanded the size of the known transcriptome. Many newly discovered transcripts are classified as long noncoding RNAs (lncRNAs) which are assumed to affect phenotype through sequence and structure and not via translated protein products despite the vast majority of them harboring short open reading frames (sORFs). Recent advances have demonstrated that the noncoding designation is incorrect in many cases and that sORF-encoded peptides (SEPs) translated from these transcripts are important contributors to diverse biological processes. Interest in SEPs is at an early stage and there is evidence for the existence of thousands of SEPs that are yet unstudied. We hope to pique interest in investigating this unexplored proteome by providing a discussion of SEP characterization generally and describing specific discoveries in innate immunity.

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.

Hidden in plain sight: challenges in proteomics detection of small ORF-encoded polypeptides

Genomic studies of bacteria have long pointed toward widespread prevalence of small open reading frames (sORFs) encoding for short proteins, <100 amino acids in length. Despite the mounting genomic evidence of their robust expression, relatively little progress has been made in their mass spectrometry-based detection and various blanket statements have been used to explain this observed discrepancy. In this study, we provide a large-scale riboproteogenomics investigation of the challenging nature of proteomic detection of such small proteins as informed by conditional translation data. A panel of physiochemical properties alongside recently developed mass spectrometry detectability metrics was interrogated to provide a comprehensive evidence-based assessment of sORF-encoded polypeptide (SEP) detectability. Moreover, a large-scale proteomics and translatomics compendium of proteins produced by Salmonella Typhimurium (S. Typhimurium), a model human pathogen, across a panel of growth conditions is presented and used in support of our in silico SEP detectability analysis. This integrative approach is used to provide a data-driven census of small proteins expressed by S. Typhimurium across growth phases and infection-relevant conditions. Taken together, our study pinpoints current limitations in proteomics-based detection of novel small proteins currently missing from bacterial genome annotations.

The Identity of the Constriction Region of the Ribosomal Exit Tunnel Is Important to Maintain Gene Expression in Escherichia coli

Mutational changes in bacterial ribosomes often affect gene expression and consequently cellular fitness. Understanding how mutant ribosomes disrupt global gene expression is critical to determining key genetic factors that affect bacterial survival. Here, we describe gene expression and phenotypic changes presented in Escherichia coli cells carrying an uL22(K90D) mutant ribosomal protein, which displayed alterations during growth. Ribosome profiling analyses revealed reduced expression of operons involved in catabolism, indole production, and lysine-dependent acid resistance. In general, translation initiation of proximal genes in several of these affected operons was substantially reduced. These reductions in expression were accompanied by increases in the expression of acid-induced membrane proteins and chaperones, the glutamate-decarboxylase regulon, and the autoinducer-2 metabolic regulon. In agreement with these changes, uL22(K90D) mutant cells had higher glutamate decarboxylase activity, survived better in extremely acidic conditions, and generated more biofilm in static cultures compared to their parental strain. Our work demonstrates that a single mutation in a non-conserved residue of a ribosomal protein affects a substantial number of genes to alter pH resistance and the formation of biofilms. IMPORTANCE All newly synthesized proteins must pass through a channel in the ribosome named the exit tunnel before emerging into the cytoplasm, membrane, and other compartments. The structural characteristics of the tunnel could govern protein folding and gene expression in a species-specific manner but how the identity of tunnel elements influences gene expression is less well-understood. Our global transcriptomics and translatome profiling demonstrate that a single substitution in a non-conserved amino acid of the E. coli tunnel protein uL22 has a profound impact on catabolism, cellular signaling, and acid resistance systems. Consequently, cells bearing the uL22 mutant ribosomes had an increased ability to survive acidic conditions and form biofilms. This work reveals a previously unrecognized link between tunnel identity and bacterial stress adaptation involving pH response and biofilm formation.

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.

N7-methylguanosine tRNA modification promotes esophageal squamous cell carcinoma tumorigenesis via the RPTOR/ULK1/autophagy axis

Mis-regulated RNA modifications promote the processing and translation of oncogenic mRNAs to facilitate cancer progression, while the molecular mechanisms remain unclear. Here we reveal that tRNA m⁷G methyltransferase complex proteins METTL1 and WDR4 are significantly up-regulated in esophageal squamous cell carcinoma (ESCC) tissues and associated with poor ESCC prognosis. In addition, METTL1 and WDR4 promote ESCC progression via the tRNA m⁷G methyltransferase activity in vitro and in vivo. Mechanistically, METTL1 or WDR4 knockdown leads to decreased expression of m⁷G-modified tRNAs and reduces the translation of a subset of oncogenic transcripts enriched in RPTOR/ULK1/autophagy pathway. Furthermore, ESCC models using Mettl1 conditional knockout and knockin mice uncover the essential function of METTL1 in promoting ESCC tumorigenesis in vivo. Our study demonstrates the important oncogenic function of mis-regulated tRNA m⁷G modification in ESCC, and suggest that targeting METTL1 and its downstream signaling axis could be a promising therapeutic target for ESCC treatment.

Deregulation of METTL1-mediated N7- methylguanosine tRNA modification can promote oncogenesis. Here, the authors report that this modification regulates the translation of proteins in both the mTOR and negative regulators of autophagy pathways, resulting in the progression of esophageal squamous cell carcinoma.

Aberrant translation regulated by METTL1/WDR4-mediated tRNA N7-methylguanosine modification drives head and neck squamous cell carcinoma progression

Background

Cancer cells selectively promote the translation of oncogenic transcripts to stimulate cancer progression. Although growing evidence has revealed that tRNA modifications and related genes participate in this process, their roles in head and neck squamous cell carcinoma (HNSCC) remain largely uncharacterized. Here, we sought to investigate the function and mechanisms of the transfer RNA (tRNA) N7‐methylguanosine (m⁷G) modification in regulating the occurrence and development of HNSCC.

Methods

Cell lost‐of‐function and gain‐of‐function assays, xenograft models, conditional knockout and knockin mouse models were used to study the physiological functions of tRNA m⁷G modification in HNSCC tumorigenesis. tRNA modification and expression profiling, mRNA translation profiling and rescue assays were performed to uncover the underlying molecular mechanisms. Single‐cell RNA sequencing (scRNA‐seq) was conducted to explore the tumor microenvironment changes.

Results

The tRNA m⁷G methyltransferase complex components Methyltransferase‐like 1 (METTL1)/WD repeat domain 4 (WDR4) were upregulated in HNSCC and associated with a poor prognosis. Functionally, METTL1/WDR4 promoted HNSCC progression and metastasis in cell‐based and transgenic mouse models. Mechanistically, ablation of METTL1 reduced the m⁷G levels of 16 tRNAs, inhibiting the translation of a subset of oncogenic transcripts, including genes related to the phosphatidylinositol‐3‐kinase/protein kinase B/mammalian target of rapamycin (PI3K/AKT/mTOR) signaling pathway. In addition, chemical modulators of the PI3K/Akt/mTOR signaling pathway reversed the effects of Mettl1 in mouse HNSCC. Furthermore, scRNA‐seq results revealed that Mettl1 knockout in mouse tumor cells altered the immune landscape and cell‐cell interaction between the tumor and stromal compartment.

Conclusions

The tRNA m⁷G methyltransferase METTL1 was found to promote the development and malignancy of HNSCC through regulating global mRNA translation, including the PI3K/AKT/mTOR signaling pathway, and found to alter immune landscape. METTL1 could be a promising treatment target for HNSCC patients.

RiboChat: a chat-style web interface for analysis and annotation of ribosome profiling data

The increasing volume of ribosome profiling (Ribo-seq) data, computational complexity of its data processing and operational handicap of related analytical procedures present a daunting set of informatics challenges. These impose a substantial barrier to researchers particularly with no or limited bioinformatics expertise in analyzing and decoding translation information from Ribo-seq data, thus driving the need for a new research paradigm for data computation and information extraction. In this knowledge base, we herein present a novel interactive web platform, RiboChat (https://db.cngb.org/ribobench/chat.html), for direct analyzing and annotating Ribo-seq data in the form of a chat conversation. It consists of a user-friendly web interface and a backend cloud-computing service. When typing a data analysis question into the chat window, the object-text detection module will be run to recognize relevant keywords from the input text. Based on the features identified in the input, individual analytics modules are then scored to find the perfect-matching candidate. The corresponding analytics module will be further executed after checking the completion status of the uploading of datasets and configured parameters. Overall, RiboChat represents an important step forward in the emerging direction of next-generation data analytics and will enable the broad research community to conveniently decipher translation information embedded within Ribo-seq data.

The METTL5-TRMT112 N6-methyladenosine methyltransferase complex regulates mRNA translation via 18S rRNA methylation

Ribosomal RNAs (rRNAs) have long been known to carry chemical modifications, including 2'O-methylation, pseudouridylation, N⁶-methyladenosine (m⁶A), and N^6,6-dimethyladenosine. While the functions of many of these modifications are unclear, some are highly conserved and occur in regions of the ribosome critical for mRNA decoding. Both 28S rRNA and 18S rRNA carry single m⁶A sites, and while the methyltransferase ZCCHC4 has been identified as the enzyme responsible for the 28S rRNA m⁶A modification, the methyltransferase responsible for the 18S rRNA m⁶A modification has remained unclear. Here, we show that the METTL5-TRMT112 methyltransferase complex installs the m⁶A modification at position 1832 of human 18S rRNA. Our work supports findings that TRMT112 is required for METTL5 stability and reveals that human METTL5 mutations associated with microcephaly and intellectual disability disrupt this interaction. We show that loss of METTL5 in human cancer cell lines and in mice regulates gene expression at the translational level; additionally, Mettl5 knockout mice display reduced body size and evidence of metabolic defects. While recent work has focused heavily on m⁶A modifications in mRNA and their roles in mRNA processing and translation, we demonstrate here that deorphanizing putative methyltransferase enzymes can reveal previously unappreciated regulatory roles for m⁶A in noncoding RNAs.

Accurate Prediction of Protein Sequences for Proteogenomics Data Integration

This book chapter discusses proteogenomics data integration and provides an overview into the different omics layer involved in defining the proteome of a living organism. Various aspects of genome variability affecting either the sequence or abundance level of proteins are discussed in this book chapter, such as the effect of single-nucleotide variants or larger genomic structural variants on the proteome. Next, various sequencing technologies are introduced and discussed from a proteogenomics data integration perspective such as those providing short- and long-read sequencing and listing their respective advantages and shortcomings for accurate protein variant prediction using genomic/transcriptomics sequencing data. Finally, the various bioinformatics tools used to process and analyze DNA/RNA sequencing data are discussed with the ultimate goal of obtaining accurately predicted sample-specific protein sequences that can be used as a drop-in replacement in existing approaches for peptide and protein identification using popular database search engines such as MSFragger, SearchGUI/PeptideShaker.

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.

Analysis of codon-specific translation by ribosome profiling

Chemical modifications of RNA molecules can affect translation in multiple ways. Therefore, it is critical to understand how their absence changes cellular translation dynamics and in particular codon-specific translation. In this chapter, we discuss the application of ribosome profiling to analyze changes in codon-specific translation and differential translation in Saccharomyces cerevisiae and human cells.

METTL1-mediated m7G modification of Arg-TCT tRNA drives oncogenic transformation

The emerging "epitranscriptomics" field is providing insights into the biological and pathological roles of different RNA modifications. The RNA methyltransferase METTL1 catalyzes N7-methylguanosine (m7G) modification of tRNAs. Here we find METTL1 is frequently amplified and overexpressed in cancers and is associated with poor patient survival. METTL1 depletion causes decreased abundance of m7G-modified tRNAs and altered cell cycle and inhibits oncogenicity. Conversely, METTL1 overexpression induces oncogenic cell transformation and cancer. Mechanistically, we find increased abundance of m7G-modified tRNAs, in particular Arg-TCT-4-1, and increased translation of mRNAs, including cell cycle regulators that are enriched in the corresponding AGA codon. Accordingly, Arg-TCT expression is elevated in many tumor types and is associated with patient survival, and strikingly, overexpression of this individual tRNA induces oncogenic transformation. Thus, METTL1-mediated tRNA modification drives oncogenic transformation through a remodeling of the mRNA "translatome" to increase expression of growth-promoting proteins and represents a promising anti-cancer target.

N7-Methylguanosine tRNA modification enhances oncogenic mRNA translation and promotes intrahepatic cholangiocarcinoma progression

Cancer cells selectively promote translation of specific oncogenic transcripts to facilitate cancer survival and progression, but the underlying mechanisms are poorly understood. Here, we find that N⁷-methylguanosine (m⁷G) tRNA modification and its methyltransferase complex components, METTL1 and WDR4, are significantly upregulated in intrahepatic cholangiocarcinoma (ICC) and associated with poor prognosis. We further reveal the critical role of METTL1/WDR4 in promoting ICC cell survival and progression using loss- and gain-of-function assays in vitro and in vivo. Mechanistically, m⁷G tRNA modification selectively regulates the translation of oncogenic transcripts, including cell-cycle and epidermal growth factor receptor (EGFR) pathway genes, in m⁷G-tRNA-decoded codon-frequency-dependent mechanisms. Moreover, using overexpression and knockout mouse models, we demonstrate the crucial oncogenic function of Mettl1-mediated m⁷G tRNA modification in promoting ICC tumorigenesis and progression in vivo. Our study uncovers the important physiological function and mechanism of METTL1-mediated m⁷G tRNA modification in the regulation of oncogenic mRNA translation and cancer progression.

Large-scale translatome profiling annotates the functional genome and reveals the key role of genic 3' untranslated regions in translatomic variation in plants

The translatome, a profile of the translational status of genetic information within cells, provides a new perspective on gene expression. Although many plant genomes have been sequenced, comprehensive translatomic annotations are not available for plants due to a lack of efficient translatome profiling techniques. Here, we developed a new technique termed 3' ribosome-profiling sequencing (3'Ribo-seq) for reliable, robust translatomic profiling. 3'Ribo-seq combines polysome profiling and 3' selection with a barcoding and pooling strategy. Systematic translatome profiling of different tissues of Arabidopsis, rice, and maize using conventional ribosome profiling (Ribo-seq) and 3'Ribo-seq revealed many novel translational genomic loci, thereby complementing functional genome annotation in plants. Using the low-cost, efficient 3'Ribo-seq technique and genome-wide association mapping of translatome expression (eGWAS), we performed a population-level dissection of the translatomes of 159 diverse maize inbred lines and identified 1,777 translational expression quantitative trait loci (eQTLs). Notably, local eQTLs are significantly enriched in the 3' untranslated regions of genes. Detailed eQTL analysis suggested that sequence variation around the polyadenylation (polyA) signal motif plays a key role in translatomic variation. Our study provides a comprehensive translatome annotation of plant functional genomes and introduces 3'Ribo-seq, which paves the way for deep translatomic analysis at the population level.

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.