Pseudouridine profiling reveals regulated mRNA pseudouridylation in yeast and human cells

Locus-specific detection of pseudouridine with CRISPR-Cas13a

We combined the CRISPR-Cas13a system with CMC chemical labeling, developing an approach that enables precise identification of pseudouridine (Ψ) sites at specific loci within ribosomal RNA (rRNA), messenger RNA (mRNA) and small nuclear RNAs (snRNA). This method, with good efficiency and simplicity, detects Ψ sites through fluorescence measurement, providing a straightforward and fast validation for targeted Ψ sites of interest.

PseUpred-ELPSO Is an Ensemble Learning Predictor with Particle Swarm Optimizer for Improving the Prediction of RNA Pseudouridine Sites

RNA pseudouridine modification exists in different RNA types of many species, and it has a significant role in regulating the expression of biological processes. To understand the functional mechanisms for RNA pseudouridine sites, the accurate identification of pseudouridine sites in RNA sequences is essential. Although several fast and inexpensive computational methods have been proposed, the challenge of improving recognition accuracy and generalization still exists. This study proposed a novel ensemble predictor called PseUpred-ELPSO for improved RNA pseudouridine site prediction. After analyzing the nucleotide composition preferences between RNA pseudouridine site sequences, two feature representations were determined and fed into the stacking ensemble framework. Then, using five tree-based machine learning classifiers as base classifiers, 30-dimensional RNA profiles are constructed to represent RNA sequences, and using the PSO algorithm, the weights of the RNA profiles were searched to further enhance the representation. A logistic regression classifier was used as a meta-classifier to complete the final predictions. Compared to the most advanced predictors, the performance of PseUpred-ELPSO is superior in both cross-validation and the independent test. Based on the PseUpred-ELPSO predictor, a free and easy-to-operate web server has been established, which will be a powerful tool for pseudouridine site identification.

Directing RNA-modifying machineries towards endogenous RNAs: opportunities and challenges

Over 170 chemical modifications can be naturally installed on RNA, all of which are catalyzed by dedicated machineries. These modifications can alter RNA sequence structure, stability, and translation as well as serving as quality control marks that record aspects of RNA processing. The diverse roles played by RNAs within cells has motivated endeavors to exogenously introduce RNA modifications at target sites for diverse purposes ranging from recording RNA:protein interactions to therapeutic applications. Here, we discuss these applications and the approaches that have been employed to engineer RNA-modifying machineries, and highlight persisting challenges and perspectives.

RNA modification-mediated mRNA translation regulation in liver cancer: mechanisms and clinical perspectives

Malignant liver cancer is characterized by rapid tumour progression and a high mortality rate, whereas the molecular mechanisms underlying liver cancer initiation and progression are still poorly understood. The dynamic and reversible RNA modifications have crucial functions in gene expression regulation by modulating RNA processing and mRNA translation. Emerging evidence has revealed that alterations in RNA modifications facilitate the selective translation of oncogenic transcripts and promote the diverse tumorigenic processes of liver cancer. In this Review, we first highlight the current progress on the functions and mechanisms underlying RNA modifications in the regulation of mRNA translation and then summarize the exciting discoveries on aberrant RNA modification-mediated mRNA translation in the regulation of tumour initiation, metastasis, metabolism, tumour microenvironment, and drug and radiotherapy resistance in liver cancer. Finally, we discuss the diagnostic and therapeutic potentials of targeting RNA modifications and mRNA translation for the clinical management of liver cancer.

RNA methylation-related inhibitors: Biological basis and therapeutic potential for cancer therapy

RNA methylation is widespread in nature. Abnormal expression of proteins associated with RNA methylation is strongly associated with a number of human diseases including cancer. Increasing evidence suggests that targeting RNA methylation holds promise for cancer treatment. This review specifically describes several common RNA modifications, such as the relatively well-studied N6-methyladenosine, as well as 5-methylcytosine and pseudouridine (Ψ). The regulatory factors involved in these modifications and their roles in RNA are also comprehensively discussed. We summarise the diverse regulatory functions of these modifications across different types of RNAs. Furthermore, we elucidate the structural characteristics of these modifications along with the development of specific inhibitors targeting them. Additionally, recent advancements in small molecule inhibitors targeting RNA modifications are presented to underscore their immense potential and clinical significance in enhancing therapeutic efficacy against cancer. KEY POINTS: In this paper, several important types of RNA modifications and their related regulatory factors are systematically summarised. Several regulatory factors related to RNA modification types were associated with cancer progression, and their relationships with cancer cell migration, invasion, drug resistance and immune environment were summarised. In this paper, the inhibitors targeting different regulators that have been proposed in recent studies are summarised in detail, which is of great significance for the development of RNA modification regulators and cancer treatment in the future.

Fuzzy kernel evidence Random Forest for identifying pseudouridine sites

Pseudouridine is an RNA modification that is widely distributed in both prokaryotes and eukaryotes, and plays a critical role in numerous biological activities. Despite its importance, the precise identification of pseudouridine sites through experimental approaches poses significant challenges, requiring substantial time and resources.Therefore, there is a growing need for computational techniques that can reliably and quickly identify pseudouridine sites from vast amounts of RNA sequencing data. In this study, we propose fuzzy kernel evidence Random Forest (FKeERF) to identify pseudouridine sites. This method is called PseU-FKeERF, which demonstrates high accuracy in identifying pseudouridine sites from RNA sequencing data. The PseU-FKeERF model selected four RNA feature coding schemes with relatively good performance for feature combination, and then input them into the newly proposed FKeERF method for category prediction. FKeERF not only uses fuzzy logic to expand the original feature space, but also combines kernel methods that are easy to interpret in general for category prediction. Both cross-validation tests and independent tests on benchmark datasets have shown that PseU-FKeERF has better predictive performance than several state-of-the-art methods. This new method not only improves the accuracy of pseudouridine site identification, but also provides a certain reference for disease control and related drug development in the future.

RNA modifications in cellular metabolism: implications for metabolism-targeted therapy and immunotherapy

Cellular metabolism is an intricate network satisfying bioenergetic and biosynthesis requirements of cells. Relevant studies have been constantly making inroads in our understanding of pathophysiology, and inspiring development of therapeutics. As a crucial component of epigenetics at post-transcription level, RNA modification significantly determines RNA fates, further affecting various biological processes and cellular phenotypes. To be noted, immunometabolism defines the metabolic alterations occur on immune cells in different stages and immunological contexts. In this review, we characterize the distribution features, modifying mechanisms and biological functions of 8 RNA modifications, including N6-methyladenosine (m6A), N6,2'-O-dimethyladenosine (m6Am), N1-methyladenosine (m1A), 5-methylcytosine (m5C), N4-acetylcytosine (ac4C), N7-methylguanosine (m7G), Pseudouridine (Ψ), adenosine-to-inosine (A-to-I) editing, which are relatively the most studied types. Then regulatory roles of these RNA modification on metabolism in diverse health and disease contexts are comprehensively described, categorized as glucose, lipid, amino acid, and mitochondrial metabolism. And we highlight the regulation of RNA modifications on immunometabolism, further influencing immune responses. Above all, we provide a thorough discussion about clinical implications of RNA modification in metabolism-targeted therapy and immunotherapy, progression of RNA modification-targeted agents, and its potential in RNA-targeted therapeutics. Eventually, we give legitimate perspectives for future researches in this field from methodological requirements, mechanistic insights, to therapeutic applications.

Pseudouridine synthase 1 regulates erythropoiesis via transfer RNAs pseudouridylation and cytoplasmic translation

Pseudouridylation plays a regulatory role in various physiological and pathological processes. A prime example is the mitochondrial myopathy, lactic acidosis, and sideroblastic anemia syndrome (MLASA), characterized by defective pseudouridylation resulting from genetic mutations in pseudouridine synthase 1 (PUS1). However, the roles and mechanisms of pseudouridylation in normal erythropoiesis and MLASA-related anemia remain elusive. We established a mouse model carrying a point mutation (R110W) in the enzymatic domain of PUS1, mimicking the common mutation in human MLASA. Pus1-mutant mice exhibited anemia at 4 weeks old. Impaired mitochondrial oxidative phosphorylation was also observed in mutant erythroblasts. Mechanistically, mutant erythroblasts showed defective pseudouridylation of targeted tRNAs, altered tRNA profiles, decreased translation efficiency of ribosomal protein genes, and reduced globin synthesis, culminating in ineffective erythropoiesis. Our study thus provided direct evidence that pseudouridylation participates in erythropoiesis in vivo. We demonstrated the critical role of pseudouridylation in regulating tRNA homeostasis, cytoplasmic translation, and erythropoiesis.

Emerging Roles and Mechanisms of RNA Modifications in Neurodegenerative Diseases and Glioma

Despite a long history of research, neurodegenerative diseases and malignant brain tumor gliomas are both considered incurable, facing challenges in the development of treatments. Recent evidence suggests that RNA modifications, previously considered as static components of intracellular RNAs, are in fact dynamically regulated across various RNA species in cells and play a critical role in major biological processes in the nervous system. Innovations in next-generation sequencing have enabled the accurate detection of modifications on bases and sugars within various RNA molecules. These RNA modifications influence the stability and transportation of RNA, and crucially affect its translation. This review delves into existing knowledge on RNA modifications to offer a comprehensive inventory of these modifications across different RNA species. The detailed regulatory functions and roles of RNA modifications within the nervous system are discussed with a focus on neurodegenerative diseases and gliomas. This article presents a comprehensive overview of the fundamental mechanisms and emerging roles of RNA modifications in these diseases, which can facilitate the creation of innovative diagnostics and therapeutics for these conditions.

RNA-binding proteins in pain

RNAs are meticulously controlled by proteins. Through direct and indirect associations, every facet in the brief life of an mRNA is subject to regulation. RNA-binding proteins (RBPs) permeate biology. Here, we focus on their roles in pain. Chronic pain is among the largest challenges facing medicine and requires new strategies. Mounting pharmacologic and genetic evidence obtained in pre-clinical models suggests fundamental roles for a broad array of RBPs. We describe their diverse roles that span RNA modification, splicing, stability, translation, and decay. Finally, we highlight opportunities to expand our understanding of regulatory interactions that contribute to pain signaling. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications Translation > Regulation RNA in Disease and Development > RNA in Disease.

Interpretable Multi-Scale Deep Learning for RNA Methylation Analysis across Multiple Species

RNA modification plays a crucial role in cellular regulation. However, traditional high-throughput sequencing methods for elucidating their functional mechanisms are time-consuming and labor-intensive, despite extensive research. Moreover, existing methods often limit their focus to specific species, neglecting the simultaneous exploration of RNA modifications across diverse species. Therefore, a versatile computational approach is necessary for interpretable analysis of RNA modifications across species. A multi-scale biological language-based deep learning model is proposed for interpretable, sequential-level prediction of diverse RNA modifications. Benchmark comparisons across species demonstrate the model's superiority in predicting various RNA methylation types over current state-of-the-art methods. The cross-species validation and attention weight visualization also highlight the model's capability to capture sequential and functional semantics from genomic backgrounds. Our analysis of RNA modifications helps us find the potential existence of "biological grammars" in each modification type, which could be effective for mapping methylation-related sequential patterns and understanding the underlying biological mechanisms of RNA modifications.

How snoRNAs can contribute to cancer at multiple levels

snoRNAs are a class of non-coding RNAs known to guide site specifically RNA modifications such as 2'-O-methylation and pseudouridylation. Recent results regarding snoRNA alterations in cancer has been made available and suggest their potential evaluation as diagnostic and prognostic biomarkers. A large part of these data, however, was not consistently confirmed and failed to provide mechanistic insights on the contribution of altered snoRNA expression to the neoplastic process. Here, we aim to critically review the available literature on snoRNA in cancer focusing on the studies elucidating the functional consequences of their deregulation. Beyond the canonical guide function in RNA processing and modification we also considered additional roles in which snoRNA, in various forms and through different modalities, are involved and that have been recently reported.

General Principles and Limitations for Detection of RNA Modifications by Sequencing

Among the many analytical methods applied to RNA modifications, a particularly pronounced surge has occurred in the past decade in the field of modification mapping. The occurrence of modifications such as m6A in mRNA, albeit known since the 1980s, became amenable to transcriptome-wide analyses through the advent of next-generation sequencing techniques in a rather sudden manner. The term "mapping" here refers to detection of RNA modifications in a sequence context, which has a dramatic impact on the interpretation of biological functions. As a consequence, an impressive number of mapping techniques were published, most in the perspective of what now has become known as "epitranscriptomics". While more and more different modifications were reported to occur in mRNA, conflicting reports and controversial results pointed to a number of technical and theoretical problems rooted in analytics, statistics, and reagents. Rather than finding the proverbial needle in a haystack, the tasks were to determine how many needles of what color in what size of a haystack one was looking at.As the authors of this Account, we think it important to outline the limitations of different mapping methods since many life scientists freshly entering the field confuse the accuracy and precision of modification mapping with that of normal sequencing, which already features numerous caveats by itself. Indeed, we propose here to qualify a specific mapping method by the size of the transcriptome that can be meaningfully analyzed with it.We here focus on high throughput sequencing by Illumina technology, referred to as RNA-Seq. We noted with interest the development of methods for modification detection by other high throughput sequencing platforms that act directly on RNA, e.g., PacBio SMRT and nanopore sequencing, but those are not considered here.In contrast to approaches relying on direct RNA sequencing, current Illumina RNA-Seq protocols require prior conversion of RNA into DNA. This conversion relies on reverse transcription (RT) to create cDNA; thereafter, the cDNA undergoes a sequencing-by-synthesis type of analysis. Thus, a particular behavior of RNA modified nucleotides during the RT-step is a prerequisite for their detection (and quantification) by deep sequencing, and RT properties have great influence on the detection efficiency and reliability. Moreover, the RT-step requires annealing of a synthetic primer, a prerequisite with a crucial impact on library preparation. Thus, all RNA-Seq protocols must feature steps for the introduction of primers, primer landing sites, or adapters on both the RNA 3'- and 5'-ends.

BID-seq for transcriptome-wide quantitative sequencing of mRNA pseudouridine at base resolution

Pseudouridine (Ψ) is an abundant RNA modification that is present in and affects the functions of diverse non-coding RNA species, including rRNA, tRNA and small nuclear RNA. Ψ also exists in mammalian mRNA and probably exhibits functional roles; however, functional investigations of mRNA Ψ modifications in mammals have been hampered by the lack of a quantitative method that detects Ψ at base precision. We have recently developed bisulfite-induced deletion sequencing (BID-seq), which provides the community with a quantitative method to map RNA Ψ distribution transcriptome-wide at single-base resolution. Here, we describe an optimized BID-seq protocol for mapping Ψ distribution across cellular mRNAs, which includes fast steps in both library preparation and data analysis. This protocol generates highly reproducible results by inducing high deletion ratios at Ψ modification within diverse sequence contexts, and meanwhile displayed almost zero background deletions at unmodified uridines. When used for transcriptome-wide Ψ profiling in mouse embryonic stem cells, the current protocol uncovered 8,407 Ψ sites from as little as 10 ng of polyA+ RNA input. This optimized BID-seq workflow takes 5 days to complete and includes four main sections: RNA preparation, library construction, next-generation sequencing (NGS) and data analysis. Library construction can be completed by researchers who have basic knowledge and skills in molecular biology and genetics. In addition to the experimental protocol, we provide BID-pipe ( https://github.com/y9c/pseudoU-BIDseq ), a user-friendly data analysis pipeline for Ψ site detection and modification stoichiometry quantification, requiring only basic bioinformatic and computational skills to uncover Ψ signatures from BID-seq data.

RNA modifications in physiology and disease: towards clinical applications

The ability of chemical modifications of single nucleotides to alter the electrostatic charge, hydrophobic surface and base pairing of RNA molecules is exploited for the clinical use of stable artificial RNAs such as mRNA vaccines and synthetic small RNA molecules - to increase or decrease the expression of therapeutic proteins. Furthermore, naturally occurring biochemical modifications of nucleotides regulate RNA metabolism and function to modulate crucial cellular processes. Studies showing the mechanisms by which RNA modifications regulate basic cell functions in higher organisms have led to greater understanding of how aberrant RNA modification profiles can cause disease in humans. Together, these basic science discoveries have unravelled the molecular and cellular functions of RNA modifications, have provided new prospects for therapeutic manipulation and have led to a range of innovative clinical approaches.

In silico λ-dynamics predicts protein binding specificities to modified RNAs

RNA modifications shape gene expression through a smorgasbord of chemical changes to canonical RNA bases. Although numbering in the hundreds, only a few RNA modifications are well characterized, in part due to the absence of methods to identify modification sites. Antibodies remain a common tool to identify modified RNA and infer modification sites through straightforward applications. However, specificity issues can result in off-target binding and confound conclusions. This work utilizes in silico λ-dynamics to efficiently estimate binding free energy differences of modification-targeting antibodies between a variety of naturally occurring RNA modifications. Crystal structures of inosine and N6-methyladenosine (m6A) targeting antibodies bound to their modified ribonucleosides were determined and served as structural starting points. λ-Dynamics was utilized to predict RNA modifications that permit or inhibit binding to these antibodies. In vitro RNA-antibody binding assays supported the accuracy of these in silico results. High agreement between experimental and computed binding propensities demonstrated that λ-dynamics can serve as a predictive screen for antibody specificity against libraries of RNA modifications. More importantly, this strategy is an innovative way to elucidate how hundreds of known RNA modifications interact with biological molecules without the limitations imposed by in vitro or in vivo methodologies.

2.7 Å cryo-EM structure of human telomerase H/ACA ribonucleoprotein

Telomerase is a ribonucleoprotein (RNP) enzyme that extends telomeric repeats at eukaryotic chromosome ends to counterbalance telomere loss caused by incomplete genome replication. Human telomerase is comprised of two distinct functional lobes tethered by telomerase RNA (hTR): a catalytic core, responsible for DNA extension; and a Hinge and ACA (H/ACA) box RNP, responsible for telomerase biogenesis. H/ACA RNPs also have a general role in pseudouridylation of spliceosomal and ribosomal RNAs, which is critical for the biogenesis of the spliceosome and ribosome. Much of our structural understanding of eukaryotic H/ACA RNPs comes from structures of the human telomerase H/ACA RNP. Here we report a 2.7 Å cryo-electron microscopy structure of the telomerase H/ACA RNP. The significant improvement in resolution over previous 3.3 Å to 8.2 Å structures allows us to uncover new molecular interactions within the H/ACA RNP. Many disease mutations are mapped to these interaction sites. The structure also reveals unprecedented insights into a region critical for pseudouridylation in canonical H/ACA RNPs. Together, our work advances understanding of telomerase-related disease mutations and the mechanism of pseudouridylation by eukaryotic H/ACA RNPs.

Role of RNA Modifications, Especially m6A, in Aflatoxin Biosynthesis of Aspergillus flavus

RNA modifications play key roles in eukaryotes, but the functions in Aspergillus flavus are still unknown. Temperature has been reported previously to be a critical environmental factor that regulates the aflatoxin production of A. flavus, but much remains to be learned about the molecular networks. Here, we demonstrated that 12 kinds of RNA modifications in A. flavus were significantly changed under 29 °C compared to 37 °C incubation; among them, m6A was further verified by a colorimetric method. Then, the transcriptome-wide m6A methylome and m6A-altered genes were comprehensively illuminated through methylated RNA immunoprecipitation sequencing and RNA sequencing, from which 22 differentially methylated and expressed transcripts under 29 °C were screened out. It is especially notable that AFCA_009549, an aflatoxin biosynthetic pathway gene (aflQ), and the m6A methylation of its 332nd adenine in the mRNA significantly affect aflatoxin biosynthesis in A. flavus both on media and crop kernels. The content of sterigmatocystin in both ΔaflQ and aflQA332C strains was significantly higher than that in the WT strain. Together, these findings reveal that RNA modifications are associated with secondary metabolite biosynthesis of A. flavus.

Pseudouridylation-mediated gene expression modulation

RNA-guided pseudouridylation, a widespread post-transcriptional RNA modification, has recently gained recognition for its role in cellular processes such as pre-mRNA splicing and the modulation of premature termination codon (PTC) readthrough. This review provides insights into its mechanisms, functions, and potential therapeutic applications. It examines the mechanisms governing RNA-guided pseudouridylation, emphasizing the roles of guide RNAs and pseudouridine synthases in catalyzing uridine-to-pseudouridine conversion. A key focus is the impact of RNA-guided pseudouridylation of U2 small nuclear RNA on pre-mRNA splicing, encompassing its influence on branch site recognition and spliceosome assembly. Additionally, the review discusses the emerging role of RNA-guided pseudouridylation in regulating PTC readthrough, impacting translation termination and genetic disorders. Finally, it explores the therapeutic potential of pseudouridine modifications, offering insights into potential treatments for genetic diseases and cancer and the development of mRNA vaccine.

Detecting m6A RNA modification from nanopore sequencing using a semi-supervised learning framework

Direct nanopore-based RNA sequencing can be used to detect post-transcriptional base modifications, such as m6A methylation, based on the electric current signals produced by the distinct chemical structures of modified bases. A key challenge is the scarcity of adequate training data with known methylation modifications. We present Xron, a hybrid encoder-decoder framework that delivers a direct methylation-distinguishing basecaller by training on synthetic RNA data and immunoprecipitation-based experimental data in two steps. First, we generate data with more diverse modification combinations through in silico cross-linking. Second, we use this dataset to train an end-to-end neural network basecaller followed by fine-tuning on immunoprecipitation-based experimental data with label-smoothing. The trained neural network basecaller outperforms existing methylation detection methods on both read-level and site-level prediction scores. Xron is a standalone, end-to-end m6A-distinguishing basecaller capable of detecting methylated bases directly from raw sequencing signals, enabling de novo methylome assembly.

RMBase v3.0: decode the landscape, mechanisms and functions of RNA modifications

Although over 170 chemical modifications have been identified, their prevalence, mechanism and function remain largely unknown. To enable integrated analysis of diverse RNA modification profiles, we have developed RMBase v3.0 (http://bioinformaticsscience.cn/rmbase/), a comprehensive platform consisting of eight modules. These modules facilitate the exploration of transcriptome-wide landscape, biogenesis, interactome and functions of RNA modifications. By mining thousands of epitranscriptome datasets with novel pipelines, the 'RNA Modifications' module reveals the map of 73 RNA modifications of 62 species. the 'Genes' module allows to retrieve RNA modification profiles and clusters by gene and transcript. The 'Mechanisms' module explores 23 382 enzyme-catalyzed or snoRNA-guided modified sites to elucidate their biogenesis mechanisms. The 'Co-localization' module systematically formulates potential correlations between 14 histone modifications and 6 RNA modifications in various cell-lines. The 'RMP' module investigates the differential expression profiles of 146 RNA-modifying proteins (RMPs) in 18 types of cancers. The 'Interactome' integrates the interactional relationships between 73 RNA modifications with RBP binding events, miRNA targets and SNPs. The 'Motif' illuminates the enriched motifs for 11 types of RNA modifications identified from epitranscriptome datasets. The 'Tools' introduces a novel web-based 'modGeneTool' for annotating modifications. Overall, RMBase v3.0 provides various resources and tools for studying RNA modifications.

Base-Resolution Sequencing Methods for Whole-Transcriptome Quantification of mRNA Modifications

RNA molecules are not merely a combination of four bases of A, C, G, and U. Chemical modifications occur in almost all RNA species and play diverse roles in gene expression regulation. The abundant cellular RNAs, such as ribosomal RNA (rRNA) and transfer RNA (tRNA), are known to have the highest density of RNA modifications, which exert critical functions in rRNA and tRNA biogenesis, stability, and subsequent translation. In recent years, modifications on low-abundance RNA species in mammalian cells, such as messenger RNA (mRNA), regulatory noncoding RNA (ncRNA), and chromatin-associated RNA (caRNA), have been shown to contain multiple different chemical modifications with functional significance. As the most abundant mRNA modification in mammals, N6-methyladenosine (m6A) affects nearly every stage of mRNA processing and metabolism, with the antibody-based m6A-MeRIP-seq (methylated RNA immunoprecipitation sequencing) followed by high-throughput sequencing widely employed in mapping m6A distribution transcriptome-wide in diverse biological systems. In addition to m6A, other chemical modifications such as pseudouridine (Ψ), 2'-O-methylation (Nm), 5-methylcytidine (m5C), internal N7-methylguanosine (m7G), N1-methyladenosine (m1A), N4-acetylcytidine (ac4C), etc. also exist in polyA-tailed RNA in mammalian cells, requiring effective mapping approaches for whole-transcriptome profiling of these non-m6A mRNA modifications. Like m6A, the antibody-based enrichment followed by sequencing has been the primary method to study distributions of these modifications. Methods to more quantitatively map these modifications would dramatically improve our understanding of distributions and modification density of these chemical marks on RNA, thereby bettering informing functional implications. In this Account, aimed at both single-base resolution and modification fraction quantification, we summarize our recent advances in developing a series of chemistry- or biochemistry-based methods to quantitatively map RNA modifications, including m6A, Ψ, m5C, m1A, 2'-O-methylation (Nm), and internal m7G, in mammalian mRNA at base resolution. These new methods, including m6A-SAC-seq, eTAM-seq, BID-seq, UBS-seq, DAMM-seq, m1A-quant-seq, Nm-Mut-seq, and m7G-quant-seq, promise to conduct base-resolution mapping of most major mRNA modifications with low RNA input and uncover dynamic changes in modification stoichiometry during biological and physiological processes, facilitating future investigations on these RNA modifications in regulating cellular gene expression and as potential biomarkers for clinical diagnosis and prognosis. These quantitative sequencing methods allow the mapping of most mRNA modifications with limited input sample requirements. The same modifications on diverse RNA species, such as caRNA, ncRNA, nuclear nascent RNA, mitochondrial RNA, cell-free RNA (cfRNA), etc., could be sequenced using the same methods.

Why U matters: detection and functions of pseudouridine modifications in mRNAs

The uridine modifications pseudouridine (Ψ), dihydrouridine, and 5-methyluridine are present in eukaryotic mRNAs. Many uridine-modifying enzymes are associated with human disease, underscoring the importance of uncovering the functions of uridine modifications in mRNAs. These modified uridines have chemical properties distinct from those of canonical uridines, which impact RNA structure and RNA-protein interactions. Ψ, the most abundant of these uridine modifications, is present across (pre-)mRNAs. Recent work has shown that many Ψs are present at intermediate to high stoichiometries that are likely conducive to function and at locations that are poised to influence pre-/mRNA processing. Technological innovations and mechanistic investigations are unveiling the functions of uridine modifications in pre-mRNA splicing, translation, and mRNA stability, which are discussed in this review.

Quantitative analysis of tRNA abundance and modifications by nanopore RNA sequencing

Transfer RNAs (tRNAs) play a central role in protein translation. Studying them has been difficult in part because a simple method to simultaneously quantify their abundance and chemical modifications is lacking. Here we introduce Nano-tRNAseq, a nanopore-based approach to sequence native tRNA populations that provides quantitative estimates of both tRNA abundances and modification dynamics in a single experiment. We show that default nanopore sequencing settings discard the vast majority of tRNA reads, leading to poor sequencing yields and biased representations of tRNA abundances based on their transcript length. Re-processing of raw nanopore current intensity signals leads to a 12-fold increase in the number of recovered tRNA reads and enables recapitulation of accurate tRNA abundances. We then apply Nano-tRNAseq to Saccharomyces cerevisiae tRNA populations, revealing crosstalks and interdependencies between different tRNA modification types within the same molecule and changes in tRNA populations in response to oxidative stress.

Single-cell analysis of the epitranscriptome: RNA modifications under the microscope

The identification of mechanisms capable of modifying genetic information by the addition of covalent RNA modifications distinguishes a level of complexity in gene expression which challenges key long-standing concepts of RNA biology. One of the current challenges of molecular biology is to properly understand the molecular functions of these RNA modifications, with more than 170 different ones having been identified so far. However, it has not been possible to map specific RNA modifications at a single-cell resolution until very recently. This review will highlight the technological advances in single-cell methodologies aimed at assessing and testing the biological function of certain RNA modifications, focusing on m6A. These advances have allowed for the development of novel strategies that enable the study of the 'epitranscriptome'. Nevertheless, despite all these improvements, many challenges and difficulties still need fixing for these techniques to work efficiently.

Emerging Role of Environmental Epitranscriptomics and RNA Modifications in Parkinson's Disease

Environmental risk factors and gene-environment interactions play a critical role in Parkinson's disease (PD). However, the relatively large contribution of environmental risk factors in the overwhelming majority of PD cases has been widely neglected in the field. A "PD prevention agenda" proposed in this journal laid out a set of research priorities focused on preventing PD through modification of environmental risk factors. This agenda includes a call for preclinical studies to employ new high-throughput methods for analyzing transcriptomics and epigenomics to provide a deeper understanding of the effects of exposures linked to PD. Here, we focus on epitranscriptomics as a novel area of research with the potential to add to our understanding of the interplay between genes and environmental exposures in PD. Both epigenetics and epitranscriptomics have been recognized as potential mediators of the complex relationship between genes, environment, and disease. Multiple studies have identified epigenetic alterations, such as DNA methylation, associated with PD and PD-related exposures in human studies and preclinical models. In addition, recent technological advancements have made it possible to study epitranscriptomic RNA modifications, such as RNA N6-methyladenosine (m6A), and a handful of recent studies have begun to explore epitranscriptomics in PD-relevant exposure models. Continued exploration of epitranscriptomic mechanisms in environmentally relevant PD models offers the opportunity to identify biomarkers, pre-degenerative changes that precede symptom onset, and potential mitigation strategies for disease prevention and treatment.

Analyzing RNA posttranscriptional modifications to decipher the epitranscriptomic code

The discovery of RNA silencing has revealed that non-protein-coding sequences (ncRNAs) can cover essential roles in regulatory networks and their malfunction may result in severe consequences on human health. These findings have prompted a general reassessment of the significance of RNA as a key player in cellular processes. This reassessment, however, will not be complete without a greater understanding of the distribution and function of the over 170 variants of the canonical ribonucleotides, which contribute to the breathtaking structural diversity of natural RNA. This review surveys the analytical approaches employed for the identification, characterization, and detection of RNA posttranscriptional modifications (rPTMs). The merits of analyzing individual units after exhaustive hydrolysis of the initial biopolymer are outlined together with those of identifying their position in the sequence of parent strands. Approaches based on next generation sequencing and mass spectrometry technologies are covered in depth to provide a comprehensive view of their respective merits. Deciphering the epitranscriptomic code will require not only mapping the location of rPTMs in the various classes of RNAs, but also assessing the variations of expression levels under different experimental conditions. The fact that no individual platform is currently capable of meeting all such demands implies that it will be essential to capitalize on complementary approaches to obtain the desired information. For this reason, the review strived to cover the broadest possible range of techniques to provide readers with the fundamental elements necessary to make informed choices and design the most effective possible strategy to accomplish the task at hand.

Genome-wide identification, characterization, and expression analysis unveil the roles of pseudouridine synthase (PUS) family proteins in rice development and stress response

Pseudouridylation, the conversion of uridine (U) to pseudouridine (Ѱ), is one of the most prevalent and evolutionary conserved RNA modifications, which is catalyzed by pseudouridine synthase (PUS) enzymes. Ѱs play a crucial epitranscriptomic role by regulating attributes of cellular RNAs across diverse organisms. However, the precise biological functions of PUSs in plants remain largely elusive. In this study, we identified and characterized 21 members in the rice PUS family which were categorized into six distinct subfamilies, with RluA and TruA emerging as the most extensive. A comprehensive analysis of domain structures, motifs, and homology modeling revealed that OsPUSs possess all canonical features of true PUS proteins, essential for substrate recognition and catalysis. The exploration of OsPUS promoters revealed presence of cis-acting regulatory elements associated with hormone and abiotic stress responses. Expression analysis of OsPUS genes showed differential expression at developmental stages and under stress conditions. Notably, OsTruB3 displayed high expression in salt, heat, and drought stresses. Several OsRluA members showed induction in heat stress, while a significant decline in expression was observed for various OsTruA members in drought and salinity. Furthermore, miRNAs predicted to target OsPUSs were themselves responsive to variable stresses, adding an additional layer of regulatory complexity of OsPUSs. Study of protein-protein interaction networks provided substantial support for the potential regulatory role of OsPUSs in numerous cellular and stress response pathways. Conclusively, our study provides functional insights into the OsPUS family, contributing to a better understanding of their crucial roles in shaping the development and stress adaptation in rice.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01396-4.

tRNA modifications: insights into their role in human cancers

Transfer RNA (tRNA) plays a central role in translation by functioning as a biological link between messenger RNA (mRNA) and proteins. One prominent feature of the tRNA molecule is its heavily modified status, which greatly affects its biogenesis and function. Modifications within the anticodon loop are crucial for translation efficiency and accuracy, whereas other modifications in the body region affect tRNA structure and stability. Recent research has revealed that these diverse modifications are critical regulators of gene expression. They are involved in many important physiological and pathological processes, including cancers. In this review we focus on six different tRNA modifications to delineate their functions and mechanisms in tumorigenesis and tumor progression, providing insights into their clinical potential as biomarkers and therapeutic targets.

Noncoding snoRNA host genes are a distinct subclass of long noncoding RNAs

Mammalian genomes are pervasively transcribed into different noncoding (nc)RNA classes, each one with its own hallmarks and exceptions. Some of them are nested into each other, such as host genes for small nucleolar RNAs (snoRNAs), which were long believed to simply act as molecular containers strictly facilitating snoRNA biogenesis. However, recent findings show that noncoding snoRNA host genes (ncSNHGs) display features different from those of 'regular' long ncRNAs (lncRNAs) and, more importantly, they can exert independent and unrelated functions to those of the encoded snoRNAs. Here, we review and summarize past and recent evidence that ncSNHGs form a defined subclass among the plethora of lncRNAs, and discuss future research that can further elucidate their biological relevance.

RNA pseudouridine modification in plants

Pseudouridine is one of the well-known chemical modifications in various RNA species. Current advances to detect pseudouridine show that the pseudouridine landscape is dynamic and affects multiple cellular processes. Although our understanding of this post-transcriptional modification mainly depends on yeast and human models, the recent findings provide strong evidence for the critical role of pseudouridine in plants. Here, we review the current knowledge of pseudouridine in plant RNAs, including its synthesis, degradation, regulatory mechanisms, and functions. Moreover, we propose future areas of research on pseudouridine modification in plants.

A single pseudouridine on rRNA regulates ribosome structure and function in the mammalian parasite Trypanosoma brucei

Trypanosomes are protozoan parasites that cycle between insect and mammalian hosts and are the causative agent of sleeping sickness. Here, we describe the changes of pseudouridine (Ψ) modification on rRNA in the two life stages of the parasite using four different genome-wide approaches. CRISPR-Cas9 knock-outs of all four snoRNAs guiding Ψ on helix 69 (H69) of the large rRNA subunit were lethal. A single knock-out of a snoRNA guiding Ψ530 on H69 altered the composition of the 80S monosome. These changes specifically affected the translation of only a subset of proteins. This study correlates a single site Ψ modification with changes in ribosomal protein stoichiometry, supported by a high-resolution cryo-EM structure. We propose that alteration in rRNA modifications could generate ribosomes preferentially translating state-beneficial proteins.

MRT-ModSeq - Rapid Detection of RNA Modifications with MarathonRT

Chemical modifications are essential regulatory elements that modulate the behavior and function of cellular RNAs. Despite recent advances in sequencing-based RNA modification mapping, methods combining accuracy and speed are still lacking. Here, we introduce MRT-ModSeq for rapid, simultaneous detection of multiple RNA modifications using MarathonRT. MRT-ModSeq employs distinct divalent cofactors to generate 2-D mutational profiles that are highly dependent on nucleotide identity and modification type. As a proof of concept, we use the MRT fingerprints of well-studied rRNAs to implement a general workflow for detecting RNA modifications. MRT-ModSeq rapidly detects positions of diverse modifications across a RNA transcript, enabling assignment of m1acp3Y, m1A, m3U, m7G and 2'-OMe locations through mutation-rate filtering and machine learning. m1A sites in sparsely modified targets, such as MALAT1 and PRUNE1 could also be detected. MRT-ModSeq can be trained on natural and synthetic transcripts to expedite detection of diverse RNA modification subtypes across targets of interest.

Do RNA modifications contribute to modulation of immune responses in allergic diseases?

RNA modifications have emerged as a fundamental mechanism of post-transcriptional gene regulation, playing vital roles in cellular physiology and the development of various diseases. While the investigation of RNA modifications has seen significant advancements, the exploration of their implication in allergic diseases has been comparatively overlooked. Allergic reactions, including hay fever, asthma, eczema and food allergies, result from hypersensitive immune responses, affecting a considerable population worldwide. Despite the high prevalence, the molecular mechanisms underlying these responses remain partially understood. The potential role of RNA modifications in modulating the hypersensitive immune responses has yet to be fully elucidated. This mini-review seeks to shed light on potential connections between RNA modifications and allergy, highlighting recent findings and potential future research directions. By expanding our understanding of the complex interplay between RNA modifications and allergic responses, we hope to unlock new avenues for allergy diagnosis, prognosis, and therapeutic intervention.

The structural basis of mRNA recognition and binding by yeast pseudouridine synthase PUS1

The chemical modification of RNA bases represents a ubiquitous activity that spans all domains of life. Pseudouridylation is the most common RNA modification and is observed within tRNA, rRNA, ncRNA and mRNAs. Pseudouridine synthase or 'PUS' enzymes include those that rely on guide RNA molecules and others that function as 'stand-alone' enzymes. Among the latter, several have been shown to modify mRNA transcripts. Although recent studies have defined the structural requirements for RNA to act as a PUS target, the mechanisms by which PUS1 recognizes these target sequences in mRNA are not well understood. Here we describe the crystal structure of yeast PUS1 bound to an RNA target that we identified as being a hot spot for PUS1-interaction within a model mRNA at 2.4 Å resolution. The enzyme recognizes and binds both strands in a helical RNA duplex, and thus guides the RNA containing the target uridine to the active site for subsequent modification of the transcript. The study also allows us to show the divergence of related PUS1 enzymes and their corresponding RNA target specificities, and to speculate on the basis by which PUS1 binds and modifies mRNA or tRNA substrates.

Quantifying m6A and Ψ Modifications in the Transcriptome via Chemical-Assisted Approaches

Since the discovery of the first chemically modified RNA nucleotide in 1951, more than 170 types of chemical modifications have been characterized in RNA so far. Since the discovery of the reversible and dynamic nature of N6-methyladenosine (m6A) in mRNA modification, researchers have identified about ten modifications in eukaryotic mRNA; together with modifications on the noncoding RNAs, the term "epitranscriptome" has been coined to describe the ensemble of various chemical RNA modifications. The past decade has witnessed the discovery of many novel molecular functions of mRNA modifications, demonstrating their crucial roles in gene expression regulation. As the most abundant modifications in mRNA, the study of m6A and Ψ has been facilitated by innovative high-throughput sequencing technologies, which can be based on antibodies, enzymes, or novel chemistry. Among them, chemical-assisted methods utilize selective chemistry that can discriminate modified versus unmodified nucleotides, enabling the transcriptome-wide mapping of m6A and Ψ modifications and functional studies.Our group has developed several sequencing technologies to investigate these epitranscriptomic marks including m6A, Ψ, m1A, and m6Am. Among them, we have recently developed two methods for absolute quantification of m6A and Ψ in the transcriptome based on chemical reactivity to distinguish and measure the two modifications. In GLORI, we used glyoxal and nitrite to mediate efficient deamination of regular adenosine, while m6A remained unaffected, thereby enabling efficient and unbiased detection of single-base resolution and absolute quantification of m6A modification. In CeU-seq and PRAISE, we used different chemistry to achieve selective labeling and detection of transcriptome-wide Ψ. CeU-seq is based on an azido-derivatized carbodiimide compound, while PRAISE utilizes the unique activity of bisulfite to Ψ. PRAISE results in the formation of ring-opening Ψ-bisulfite adduct and quantitatively detects Ψ as 1-2 nt deletion signatures during sequencing. The resulting base-resolution and stoichiometric information expanded our understanding to the profiles of RNA modifications in the transcriptome. In particular, the quantitative information on RNA methylome is critical for characterizing the dynamic and reversible nature of RNA modifications, for instance, under environmental stress or during development. Additionally, base-resolution and stoichiometric information can greatly facilitate the analysis and characterization of functional modification sites that are important for gene expression regulation, especially when one modification type may have multiple or even opposing functions within a specific transcript. Together, the quantitative profiling methods provide the modification stoichiometry information, which is critical to study the regulatory roles of RNA modifications.In this Account, we will focus on the quantitative sequencing technologies of m6A and Ψ developed in our group, review recent advances in chemical-assisted reactions for m6A and Ψ detection, and discuss the challenges and future opportunities of transcriptome-wide mapping technologies for RNA modifications.

Pseudouridine guides germline small RNA transport and epigenetic inheritance

Epigenetic modifications that arise during plant and animal development, such as DNA and histone modification, are mostly reset during gamete formation, but some are inherited from the germline including those marking imprinted genes1. Small RNAs guide these epigenetic modifications, and some are also inherited by the next generation2,3. In C. elegans, these inherited small RNAs have poly (UG) tails4, but how inherited small RNAs are distinguished in other animals and plants is unknown. Pseudouridine (Ψ) is the most abundant RNA modification but has not been explored in small RNAs. Here, we develop novel assays to detect Ψ in short RNA sequences, demonstrating its presence in mouse and Arabidopsis microRNAs and their precursors. We also detect substantial enrichment in germline small RNAs, namely epigenetically activated siRNAs (easiRNAs) in Arabidopsis pollen, and piwi-interacting piRNAs in mouse testis. In pollen, pseudouridylated easiRNAs are localized to sperm cells, and we found that PAUSED/HEN5 (PSD), the plant homolog of Exportin-t, interacts genetically with Ψ and is required for transport of easiRNAs into sperm cells from the vegetative nucleus. We further show that Exportin-t is required for the triploid block: chromosome dosage-dependent seed lethality that is epigenetically inherited from pollen. Thus, Ψ has a conserved role in marking inherited small RNAs in the germline.

Integrative analysis of multi-omics data reveals a pseudouridine-related lncRNA signature for prediction of glioma prognosis and chemoradiotherapy sensitivity

BACKGROUND

Glioblastoma is the most common type of glioma with a high incidence and poor prognosis, and effective medical treatment remains challenging. Pseudouridine (Ψ) is the first post-transcriptional modification discovered and one of the most abundant modifications to RNA. However, the prognostic value of Ψ-related lncRNAs (ΨrLs) for glioma patients has never been systematically evaluated. This study aims to construct a risk model based on ΨrLs signature and to validate the predictive efficiency of the model.

METHOD

Transcriptomic data, genomic data, and relevant clinical data of glioma patients were extracted from the Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA). ΨrLs with significant correlation with Ψ-related genes were identified, and univariate Cox regression, least absolute shrinkage and selection operator (LASSO) regression, and multivariate Cox regression were used to further select biomarkers and construct a ΨrLs signature risk model. Then, the expression of lncRNAs of ΨrLs signature in multiple glioma cell lines was detected by qPCR. Further, ROC analysis, stratification analysis, correlation analysis, survival analysis, nomogram, enrichment analysis, immune infiltration analysis, chemoradiotherapy sensitivity analysis, somatic mutation, and recurrent copy number variation (CNV) analysis were used to validate the predictive efficiency of ΨrLs signature in TCGA and CGGA datasets.

RESULTS

A four-lncRNA ΨrLs signature (DNAJC27-AS1, GDNF-AS1, ZBTB20-AS4, and DNMBP-AS1) risk model was constructed. By ROC analysis, stratified analysis, correlation analysis, survival analysis, and nomogram, the signature showed satisfactory predictive efficiency. Functional enrichment analysis revealed the differences in immune-related biological processes between high- and low-risk groups. Immune infiltration analysis showed that the high-risk group had lower tumor purity and higher stromal, immune and ESTIMATE scores. Mitoxantrone was identified as effective drug for low-risk group of glioma patients. Key genes in glioma development, including IDH1, EGFR, PTEN, etc., were differentially mutated between risk groups. The main recurrent CNVs in low-risk groups were 19q13.42 deletion and 7q34 amplification; 10q23.31 deletion and 12q14.1 in the high-risk group.

CONCLUSIONS

Our study identified a four-lncRNA ΨrLs signature that effectively predicts the prognosis of glioma patients and may serve as a diagnostic tool. Risk scores of glioma patients generated by the signature is associated with immune-related biological processes and chemoradiotherapy sensitivity. These findings may inform the development of more targeted and effective therapies for glioma patients.

Advantages and challenges associated with bisulfite-assisted nanopore direct RNA sequencing for modifications

Nanopore direct RNA sequencing is a technology that allows sequencing for epitranscriptomic modifications with the possibility of a quantitative assessment. In the present work, pseudouridine (Ψ) was sequenced with the nanopore before and after the pH 7 bisulfite reaction that yields stable ribose adducts at C1' of Ψ. The adducted sites produced greater base call errors in the form of deletion signatures compared to Ψ. Sequencing studies on E. coli rRNA and tmRNA before and after the pH 7 bisulfite reaction demonstrated that using chemically-assisted nanopore sequencing has distinct advantages for minimization of false positives and false negatives in the data. The rRNA from E. coli has 19 known U/C sequence variations that give similar base call signatures as Ψ, and therefore, are false positives when inspecting base call data; however, these sites are refractory to reacting with bisulfite as is easily observed in nanopore data. The E. coli tmRNA has a low occupancy Ψ in a pyrimidine-rich sequence context that is called a U representing a false negative; partial occupancy by Ψ is revealed after the bisulfite reaction. In a final study, 5-methylcytidine (m5C) in RNA can readily be observed after the pH 5 bisulfite reaction in which the parent C deaminates to U and the modified site does not react. This locates m5C when using bisulfite-assisted nanopore direct RNA sequencing, which is otherwise challenging to observe. The advantages and challenges of the overall approach are discussed.

The potential role of RNA modification in skin diseases, as well as the recent advances in its detection methods and therapeutic agents

RNA modification is considered as an epigenetic modification that plays an indispensable role in biological processes such as gene expression and genome editing without altering nucleotide sequence, but the molecular mechanism of RNA modification has not been discussed systematically in the development of skin diseases. This article mainly presents the whole picture of theoretical achievements on the potential role of RNA modification in dermatology. Furthermore, this article summarizes the latest advances in clinical practice related with RNA modification, including its detection methods and drug development. Based on this comprehensive review, we aim to illustrate the current blind spots and future directions of RNA modification, which may provide new insights for researchers in this field.

Regulatory role of RNA modifications in the treatment of pancreatic ductal adenocarcinoma (PDAC)

Pancreatic ductal adenocarcinoma (PDAC) is an extremely life-threatening malignancy with a relatively unfavorable prognosis. The early occurrence of metastasis and local recurrence subsequent to surgery contribute to the poor survival rates of PDAC patients, thereby limiting the effectiveness of surgical intervention. Additionally, the desmoplastic and immune-suppressive tumor microenvironment of PDAC diminishes its responsiveness to conventional treatment modalities such as chemotherapy, radiotherapy, and immunotherapy. Therefore, it is imperative to identify novel therapeutic targets for PDAC treatment. Chemical modifications are prevalent in various types of RNA and exert significant influence on their structure and functions. RNA modifications, exemplified by m6A, m5C, m1A, and Ψ, have been identified as general regulators of cellular functions. The abundance of specific modifications, such as m6A, has been correlated with cell proliferation, invasion, migration, and patient prognosis in PDAC. Pre-clinical data has indicated that manipulating RNA modification regulators could enhance the efficacy of chemotherapy, radiotherapy, and immunotherapy. Therefore, targeting RNA modifications in conjunction with current adjuvant or neoadjuvant therapy holds promise. The objective of this review is to provide a comprehensive overview of RNA modifications in PDAC treatment, encompassing their behaviors, mechanisms, and potential treatment targets. Therefore, it aims to stimulate the development of novel therapeutic approaches and future clinical trials.

RNA modification in cardiovascular disease: implications for therapeutic interventions

Cardiovascular disease (CVD) is the leading cause of death in the world, with a high incidence and a youth-oriented tendency. RNA modification is ubiquitous and indispensable in cell, maintaining cell homeostasis and function by dynamically regulating gene expression. Accumulating evidence has revealed the role of aberrant gene expression in CVD caused by dysregulated RNA modification. In this review, we focus on nine common RNA modifications: N6-methyladenosine (m6A), N1-methyladenosine (m1A), 5-methylcytosine (m5C), N7-methylguanosine (m7G), N4-acetylcytosine (ac4C), pseudouridine (Ψ), uridylation, adenosine-to-inosine (A-to-I) RNA editing, and modifications of U34 on tRNA wobble. We summarize the key regulators of RNA modification and their effects on gene expression, such as RNA splicing, maturation, transport, stability, and translation. Then, based on the classification of CVD, the mechanisms by which the disease occurs and progresses through RNA modifications are discussed. Potential therapeutic strategies, such as gene therapy, are reviewed based on these mechanisms. Herein, some of the CVD (such as stroke and peripheral vascular disease) are not included due to the limited availability of literature. Finally, the prospective applications and challenges of RNA modification in CVD are discussed for the purpose of facilitating clinical translation. Moreover, we look forward to more studies exploring the mechanisms and roles of RNA modification in CVD in the future, as there are substantial uncultivated areas to be explored.

m5U-GEPred: prediction of RNA 5-methyluridine sites based on sequence-derived and graph embedding features

5-Methyluridine (m5U) is one of the most common post-transcriptional RNA modifications, which is involved in a variety of important biological processes and disease development. The precise identification of the m5U sites allows for a better understanding of the biological processes of RNA and contributes to the discovery of new RNA functional and therapeutic targets. Here, we present m5U-GEPred, a prediction framework, to combine sequence characteristics and graph embedding-based information for m5U identification. The graph embedding approach was introduced to extract the global information of training data that complemented the local information represented by conventional sequence features, thereby enhancing the prediction performance of m5U identification. m5U-GEPred outperformed the state-of-the-art m5U predictors built on two independent species, with an average AUROC of 0.984 and 0.985 tested on human and yeast transcriptomes, respectively. To further validate the performance of our newly proposed framework, the experimentally validated m5U sites identified from Oxford Nanopore Technology (ONT) were collected as independent testing data, and in this project, m5U-GEPred achieved reasonable prediction performance with ACC of 91.84%. We hope that m5U-GEPred should make a useful computational alternative for m5U identification.

Epitranscriptomics as a New Layer of Regulation of Gene Expression in Skeletal Muscle: Known Functions and Future Perspectives

Epitranscriptomics refers to post-transcriptional regulation of gene expression via RNA modifications and editing that affect RNA functions. Many kinds of modifications of mRNA have been described, among which are N6-methyladenosine (m6A), N1-methyladenosine (m1A), 7-methylguanosine (m7G), pseudouridine (Ψ), and 5-methylcytidine (m5C). They alter mRNA structure and consequently stability, localization and translation efficiency. Perturbation of the epitranscriptome is associated with human diseases, thus opening the opportunity for potential manipulations as a therapeutic approach. In this review, we aim to provide an overview of the functional roles of epitranscriptomic marks in the skeletal muscle system, in particular in embryonic myogenesis, muscle cell differentiation and muscle homeostasis processes. Further, we explored high-throughput epitranscriptome sequencing data to identify RNA chemical modifications in muscle-specific genes and we discuss the possible functional role and the potential therapeutic applications.

Bisulfite and Nanopore Sequencing for Pseudouridine in RNA

Nucleophilic addition of bisulfite to pyrimidine bases has been known for a half century, and the reaction has been in use for at least a quarter of a century for identifying 5-methylcytidine in DNA. This account focuses on the chemistry of bisulfite with pseudouridine, an isomer of the RNA nucleoside uridine in which the uracil base is connected to C1' of ribose via C5 instead of N1. Pseudouridine, Ψ, is the most common nucleotide modification found in cellular RNA overall, in part due to its abundance in rRNAs and tRNAs. It has a stabilizing influence on RNA structure because N1 is now available for additional hydrogen bonding and because the heterocycle is slightly better at π stacking. The isomerization of U to Ψ in RNA strands is catalyzed by 13 different enzymes in humans and 11 in E. coli; some of these enzymes are implicated in disease states which is testament to the biological importance of pseudouridine in cells. Recently, pseudouridine came into the limelight as the key modification that, after N1 methylation, enables mRNA vaccines to be delivered efficiently into human tissue with minimal generation of a deleterious immunogenic response. Here we describe the bisulfite reaction with pseudouridine which gives rise to a chemical sequencing method to map the modified base in the epitranscriptome. Unlike the reaction with cytidine, the addition of bisulfite to Ψ leads irreversibly to form an adduct that is bypassed during cDNA synthesis by reverse transcriptases yielding a characteristic deletion signature. Although there were hints to the structure of the bisulfite adduct(s) 30 to 50 years ago, it took modern spectroscopic and computational methods to solve the mystery. Raman spectroscopy along with extensive NMR, ECD, and computational work led to the assignment of the major product as the (R) diastereomer of an oxygen adduct at C1' of a ring-opened pseudouridine. Mechanistically, this arose from a succession of conjugate addition, E2 elimination, and a [2,3] sigmatropic rearrangement, all of which are stereodefined reactions. A minor reaction with excess bisulfite led to the (S) isomer of a S-adducted SO3- group. Understanding structure and mechanism aided the design of a Ψ-specific sequencing reaction and guided attempts to improve the utility and specificity of the method. Separately, we have been investigating the use of nanopore direct RNA sequencing, a single-molecule method that directly analyzes RNA strands isolated from cells after end-ligation of adaptor sequences. By combining the electrical current and base-calling data from the nanopore with dwell-time analysis from the helicase employed to deliver RNA to the nanopore, we were able to map Ψ sites in nearly all sequence contexts. This analysis was employed to find Ψ residues in the SARS-CoV-2 vRNA, to analyze the sequence context effects of mRNA vaccine synthesis via in vitro transcription, and to evaluate the impact of stress on chemical modifications in the E. coli ribosome. Most recently, we found that bisulfite treatment of RNA leading to Ψ adducts could modulate the nanopore signal to help in mapping modifications of low occupancy.

RNA modification in mRNA cancer vaccines

RNA modification is manifested as chemically altered nucleotides, widely exists in diverse natural RNAs, and is closely related to RNA structure and function. Currently, mRNA-based vaccines have received great attention and rapid development as novel and mighty fighters against various diseases including cancer. The achievement of RNA vaccines in clinical application is largely attributed to some methodological innovations including the incorporation of modified nucleotides into the synthetic RNA. The selection of optimal RNA modifications aimed at reducing the instability and immunogenicity of RNA molecules is a very critical task to improve the efficacy and safety of mRNA vaccines. This review summarizes the functions of RNA modifications and their application in mRNA vaccines, highlights recent advances of mRNA vaccines in cancer immunotherapy, and provides perspectives for future development of mRNA vaccines in the context of personalized tumor therapy.

Regulation and functions of non-m6A mRNA modifications

Nucleobase modifications are prevalent in eukaryotic mRNA and their discovery has resulted in the emergence of epitranscriptomics as a research field. The most abundant internal (non-cap) mRNA modification is N6-methyladenosine (m6A), the study of which has revolutionized our understanding of post-transcriptional gene regulation. In addition, numerous other mRNA modifications are gaining great attention because of their major roles in RNA metabolism, immunity, development and disease. In this Review, we focus on the regulation and function of non-m6A modifications in eukaryotic mRNA, including pseudouridine (Ψ), N6,2'-O-dimethyladenosine (m6Am), N1-methyladenosine (m1A), inosine, 5-methylcytidine (m5C), N4-acetylcytidine (ac4C), 2'-O-methylated nucleotide (Nm) and internal N7-methylguanosine (m7G). We highlight their regulation, distribution, stoichiometry and known roles in mRNA metabolism, such as mRNA stability, translation, splicing and export. We also discuss their biological consequences in physiological and pathological processes. In addition, we cover research techniques to further study the non-m6A mRNA modifications and discuss their potential future applications.

Quantitative profiling of pseudouridylation landscape in the human transcriptome

Pseudouridine (Ψ) is an abundant post-transcriptional RNA modification in ncRNA and mRNA. However, stoichiometric measurement of individual Ψ sites in human transcriptome remains unaddressed. Here we develop 'PRAISE', via selective chemical labeling of Ψ by bisulfite to induce nucleotide deletion signature during reverse transcription, to realize quantitative assessment of the Ψ landscape in the human transcriptome. Unlike traditional bisulfite treatment, our approach is based on quaternary base mapping and revealed an ~10% median modification level for 2,209 confident Ψ sites in HEK293T cells. By perturbing pseudouridine synthases, we obtained differential mRNA targets of PUS1, PUS7, TRUB1 and DKC1, with TRUB1 targets showing the highest modification stoichiometry. In addition, we quantified known and new Ψ sites in mitochondrial mRNA catalyzed by PUS1. Collectively, we provide a sensitive and convenient method to measure transcriptome-wide Ψ; we envision this quantitative approach would facilitate emerging efforts to elucidate the function and mechanism of mRNA pseudouridylation.

The emerging role of epitranscriptome in shaping stress responses in plants

KEY MESSAGE

RNA modifications and editing changes constitute 'epitranscriptome' and are crucial in regulating the development and stress response in plants. Exploration of the epitranscriptome and associated machinery would facilitate the engineering of stress tolerance in crops. RNA editing and modifications post-transcriptionally decorate almost all classes of cellular RNAs, including tRNAs, rRNAs, snRNAs, lncRNAs and mRNAs, with more than 170 known modifications, among which m6A, Ψ, m5C, 8-OHG and C-to-U editing are the most abundant. Together, these modifications constitute the "epitranscriptome", and contribute to changes in several RNA attributes, thus providing an additional structural and functional diversification to the "cellular messages" and adding another layer of gene regulation in organisms, including plants. Numerous evidences suggest that RNA modifications have a widespread impact on plant development as well as in regulating the response of plants to abiotic and biotic stresses. High-throughput sequencing studies demonstrate that the landscapes of m6A, m5C, Am, Cm, C-to-U, U-to-G, and A-to-I editing are remarkably dynamic during stress conditions in plants. GO analysis of transcripts enriched in Ψ, m6A and m5C modifications have identified bonafide components of stress regulatory pathways. Furthermore, significant alterations in the expression pattern of genes encoding writers, readers, and erasers of certain modifications have been documented when plants are grown in challenging environments. Notably, manipulating the expression levels of a few components of RNA editing machinery markedly influenced the stress tolerance in plants. We provide updated information on the current understanding on the contribution of RNA modifications in shaping the stress responses in plants. Unraveling of the epitranscriptome has opened new avenues for designing crops with enhanced productivity and stress resilience in view of global climate change.