m1A within cytoplasmic mRNAs at single nucleotide resolution: a reconciled transcriptome-wide map

Abundant mRNA m1A modification in dinoflagellates: a new layer of gene regulation

Dinoflagellates, a class of unicellular eukaryotic phytoplankton, exhibit minimal transcriptional regulation, representing a unique model for exploring gene expression. The biosynthesis, distribution, regulation, and function of mRNA N1-methyladenosine (m1A) remain controversial due to its limited presence in typical eukaryotic mRNA. This study provides a comprehensive map of m1A in dinoflagellate mRNA and shows that m1A, rather than N6-methyladenosine (m6A), is the most prevalent internal mRNA modification in various dinoflagellate species, with an asymmetric distribution along mature transcripts. In Amphidinium carterae, we identify 6549 m1A sites characterized by a non-tRNA T-loop-like sequence motif within the transcripts of 3196 genes, many of which are involved in regulating carbon and nitrogen metabolism. Enriched within 3'UTRs, dinoflagellate mRNA m1A levels negatively correlate with translation efficiency. Nitrogen depletion further decreases mRNA m1A levels. Our data suggest that distinctive patterns of m1A modification might influence the expression of metabolism-related genes through translational control.

MePMe-seq: antibody-free simultaneous m6A and m5C mapping in mRNA by metabolic propargyl labeling and sequencing

Internal modifications of mRNA have emerged as widespread and versatile regulatory mechanism to control gene expression at the post-transcriptional level. Most of these modifications are methyl groups, making S-adenosyl-L-methionine (SAM) a central metabolic hub. Here we show that metabolic labeling with a clickable metabolic precursor of SAM, propargyl-selenohomocysteine (PSH), enables detection and identification of various methylation sites. Propargylated A, C, and G nucleosides form at detectable amounts via intracellular generation of the corresponding SAM analogue. Integration into next generation sequencing enables mapping of N6-methyladenosine (m6A) and 5-methylcytidine (m5C) sites in mRNA with single nucleotide precision (MePMe-seq). Analysis of the termination profiles can be used to distinguish m6A from 2'-O-methyladenosine (Am) and N1-methyladenosine (m1A) sites. MePMe-seq overcomes the problems of antibodies for enrichment and sequence-motifs for evaluation, which was limiting previous methodologies. Metabolic labeling via clickable SAM facilitates the joint evaluation of methylation sites in RNA and potentially DNA and proteins.

nf-core/clipseq - a robust Nextflow pipeline for comprehensive CLIP data analysis

Crosslinking and immunoprecipitation (CLIP) technologies have become a central component of the molecular biologists' toolkit to study protein-RNA interactions and thus to uncover core principles of RNA biology. There has been a proliferation of CLIP-based experimental protocols, as well as computational tools, especially for peak-calling. Consequently, there is an urgent need for a well-documented bioinformatic pipeline that enshrines the principles of robustness, reproducibility, scalability, portability and flexibility while embracing the diversity of experimental and computational CLIP tools. To address this, we present nf-core/clipseq - a robust Nextflow pipeline for quality control and analysis of CLIP sequencing data. It is part of the international nf-core community effort to develop and curate a best-practice, gold-standard set of pipelines for data analysis. The standards enabled by Nextflow and nf-core, including workflow management, version control, continuous integration and containerisation ensure that these key needs are met. Furthermore, multiple tools are implemented ( e.g. for peak-calling), alongside visualisation of quality control metrics to empower the user to make their own informed decisions based on their data. nf-core/clipseq remains under active development, with plans to incorporate newly released tools to ensure that pipeline remains up-to-date and relevant for the community. Engagement with users and developers is encouraged through the nf-core GitHub repository and Slack channel to promote collaboration. It is available at https://nf-co.re/clipseq.

N1-Methyladenosine modification of mRNA regulates neuronal gene expression and oxygen glucose deprivation/reoxygenation induction

N¹-Methyladenosine (m1A) is an abundant modification of transcripts, plays important roles in regulating mRNA structure and translation efficiency, and is dynamically regulated under stress. However, the characteristics and functions of mRNA m1A modification in primary neurons and oxygen glucose deprivation/reoxygenation (OGD/R) induced remain unclear. We first constructed a mouse cortical neuron OGD/R model and then used methylated RNA immunoprecipitation (MeRIP) and sequencing technology to demonstrate that m1A modification is abundant in neuron mRNAs and dynamically regulated during OGD/R induction. Our study suggests that Trmt10c, Alkbh3, and Ythdf3 may be m1A-regulating enzymes in neurons during OGD/R induction. The level and pattern of m1A modification change significantly during OGD/R induction, and differential methylation is closely associated with the nervous system. Our findings show that m1A peaks in cortical neurons aggregate at both the 5' and 3' untranslated regions. m1A modification can regulate gene expression, and peaks in different regions have different effects on gene expression. By analysing m1A-seq and RNA-seq data, we show a positive correlation between differentially methylated m1A peaks and gene expression. The correlation was verified by using qRT-PCR and MeRIP-RT-PCR. Moreover, we selected human tissue samples from Parkinson's disease (PD) and Alzheimer's disease (AD) patients from the Gene Expression Comprehensive (GEO) database to analyse the selected differentially expressed genes (DEGs) and differential methylation modification regulatory enzymes, respectively, and found similar differential expression results. We highlight the potential relationship between m1A modification and neuronal apoptosis following OGD/R induction. Furthermore, by mapping mouse cortical neurons and OGD/R-induced modification characteristics, we reveal the important role of m1A modification in OGD/R and gene expression regulation, providing new ideas for research on neurological damage.

Biological roles of adenine methylation in RNA

N6-Methyladenosine (m6A) is one of the most abundant modifications of the epitranscriptome and is found in cellular RNAs across all kingdoms of life. Advances in detection and mapping methods have improved our understanding of the effects of m6A on mRNA fate and ribosomal RNA function, and have uncovered novel functional roles in virtually every species of RNA. In this Review, we explore the latest studies revealing roles for m6A-modified RNAs in chromatin architecture, transcriptional regulation and genome stability. We also summarize m6A functions in biological processes such as stem-cell renewal and differentiation, brain function, immunity and cancer progression.

Determining RNA Natural Modifications and Nucleoside Analog-Labeled Sites by a Chemical/Enzyme-Induced Base Mutation Principle

The natural chemical modifications of messenger RNA (mRNA) in living organisms have shown essential roles in both physiology and pathology. The mapping of mRNA modifications is critical for interpreting their biological functions. In another dimension, the synthesized nucleoside analogs can enable chemical labeling of cellular mRNA through a metabolic pathway, which facilitates the study of RNA dynamics in a pulse-chase manner. In this regard, the sequencing tools for mapping both natural modifications and nucleoside tags on mRNA at single base resolution are highly necessary. In this work, we review the progress of chemical sequencing technology for determining both a variety of naturally occurring base modifications mainly on mRNA and a few on transfer RNA and metabolically incorporated artificial base analogs on mRNA, and further discuss the problems and prospects in the field.

Developing an m5C regulator-mediated RNA methylation modification signature to predict prognosis and immunotherapy efficacy in rectal cancer

Background

Currently, a very small number of patients with colorectal cancer (CRC) respond to immune checkpoint inhibitor (ICI) treatment. Therefore, there is an urgent need to investigate effective biomarkers to determine the responsiveness to ICI treatment. Recently, aberrant 5-methylcytosine (m⁵C) RNA modification has emerged as a key player in the pathogenesis of cancer. Thus, we aimed to explore the predictive signature based on m⁵C regulator-related genes for characterizing the immune landscapes and predicting the prognosis and response to therapies.

Methods

The Cancer Genome Atlas (TCGA) cohort was used as the training set, while GEO data sets, real-time quantitative PCR (RT-qPCR) analysis from paired frozen tissues, and immunohistochemistry (IHC) data from tissue microarray (TMA) were used for validation. We constructed a novel signature based on three m⁵C regulator-related genes in patients with rectal adenocarcinoma (READ) using a least absolute shrinkage and selection operator (LASSO)-Cox regression and unsupervised consensus clustering analyses. Additionally, we correlated the three-gene signature risk model with the tumor immune microenvironment, immunotherapy efficiency, and potential applicable drugs.

Results

The m⁵C methylation-based signature was an independent prognostic factor, where low-risk patients showed a stronger immunoreactivity phenotype and a superior response to ICI therapy. Conversely, the high-risk patients had enriched pathways of cancer hallmarks and presented immune-suppressive state, which demonstrated that they are more insensitive to immunotherapy. Additionally, the signature markedly correlated with drug susceptibility.

Conclusions

We developed a reliable m⁵C regulator-based risk model to predict the prognosis, clarify the molecular and tumor microenvironment status, and identify patients who would benefit from immunotherapy or chemotherapy. Our study could provide vital guidance to improve prognostic stratification and optimize personalized therapeutic strategies for patients with rectal cancer.

Exploring the epitranscriptome by native RNA sequencing

Chemical RNA modifications, collectively referred to as the "epitranscriptome," are essential players in fine-tuning gene expression. Our ability to analyze RNA modifications has improved rapidly in recent years, largely due to the advent of high-throughput sequencing methodologies, which typically consist of coupling modification-specific reagents, such as antibodies or enzymes, to next-generation sequencing. Recently, it also became possible to map RNA modifications directly by sequencing native RNAs using nanopore technologies, which has been applied for the detection of a number of RNA modifications, such as N6-methyladenosine (m⁶A), pseudouridine (Ψ), and inosine (I). However, the signal modulations caused by most RNA modifications are yet to be determined. A global effort is needed to determine the signatures of the full range of RNA modifications to avoid the technical biases that have so far limited our understanding of the epitranscriptome.

m1A and m6A modifications function cooperatively to facilitate rapid mRNA degradation

N⁶-Methyladenosine (m⁶A), the most abundant internal mRNA modification, affects multiple steps in gene expression. Mechanistically, the binding of YTHDF2 to m⁶A on mRNAs elicits rapid mRNA degradation by recruiting several RNA degrading enzymes. Here, we show that N¹-methyladenosine (m¹A), another type of RNA modification, accelerates rapid m⁶A RNA degradation. We identify HRSP12 as an RNA-binding protein that recognizes m¹A. The binding of HRSP12 to m¹A promotes efficient interaction of YTHDF2 with m⁶A, consequently facilitating endoribonucleolytic cleavage via the RNase P/MRP complex. Transcriptome-wide analyses also reveal that mRNAs harboring both m¹A and m⁶A are downregulated in an HRSP12-dependent manner compared with mRNAs harboring m⁶A only. Accordingly, a subset of endogenous circular RNAs that harbor m⁶A and associate with YTHDF2 in an HRSP12-dependent manner is also subjected to m¹A-facilitated rapid degradation. Together, our observations provide compelling evidence for crosstalk between different RNA modifications.

RNA nucleotide methylation: 2021 update

Among RNA modifications, transfer of methylgroups from the typical cofactor S-adenosyl-l-methionine by methyltransferases (MTases) to RNA is by far the most common reaction. Since our last review about a decade ago, the field has witnessed the re-emergence of mRNA methylation as an important mechanism in gene regulation. Attention has then spread to many other RNA species; all being included into the newly coined concept of the "epitranscriptome." The focus moved from prokaryotes and single cell eukaryotes as model organisms to higher eukaryotes, in particular to mammals. The perception of the field has dramatically changed over the past decade. A previous lack of phenotypes in knockouts in single cell organisms has been replaced by the apparition of MTases in numerous disease models and clinical investigations. Major driving forces of the field include methylation mapping techniques, as well as the characterization of the various MTases, termed "writers." The latter term has spilled over from DNA modification in the neighboring epigenetics field, along with the designations "readers," applied to mediators of biological effects upon specific binding to a methylated RNA. Furthermore "eraser" enzymes effect the newly discovered oxidative removal of methylgroups. A sense of reversibility and dynamics has replaced the older perception of RNA modification as a concrete-cast, irreversible part of RNA maturation. A related concept concerns incompletely methylated residues, which, through permutation of each site, lead to inhomogeneous populations of numerous modivariants. This review recapitulates the major developments of the past decade outlined above, and attempts a prediction of upcoming trends. This article is categorized under: RNA Processing > RNA Editing and Modification.

RNA methylation in chloroplasts or mitochondria in plants

Recent advances in our understanding of epitranscriptomic RNA methylation have expanded the complexity of gene expression regulation beyond epigenetic regulation involving DNA methylation and histone modifications. The instalment, removal, and interpretation of methylation marks on RNAs are carried out by writers (methyltransferases), erasers (demethylases), and readers (RNA-binding proteins), respectively. Contrary to an emerging body of evidence demonstrating the importance of RNA methylation in the diverse fates of RNA molecules, including splicing, export, translation, and decay in the nucleus and cytoplasm, their roles in plant organelles remain largely unclear and are only now being discovered. In particular, extremely high levels of methylation marks in chloroplast and mitochondrial RNAs suggest that RNA methylation plays essential roles in organellar biogenesis and functions in plants that are crucial for plant development and responses to environmental stimuli. Thus, unveiling the cellular components involved in RNA methylation in cell organelles is essential to better understand plant biology.

Multiplexed profiling facilitates robust m6A quantification at site, gene and sample resolution

N⁶-methyladenosine (m6A) is the most prevalent modification of messenger RNA in mammals. To interrogate its functions and dynamics, there is a critical need to quantify m6A at three levels: site, gene and sample. Current approaches address these needs in a limited manner. Here we develop m6A-seq2, relying on multiplexed m6A-immunoprecipitation of barcoded and pooled samples. m6A-seq2 allows a big increase in throughput while reducing technical variability, requirements of input material and cost. m6A-seq2 is furthermore uniquely capable of providing sample-level relative quantitations of m6A, serving as an orthogonal alternative to mass spectrometry-based approaches. Finally, we develop a computational approach for gene-level quantitation of m6A. We demonstrate that using this metric, roughly 30% of the variability in RNA half life in mouse embryonic stem cells can be explained, establishing m6A as a main driver of RNA stability. m6A-seq2 thus provides an experimental and analytic framework for dissecting m6A-mediated regulation at three different levels.

Analysis of pseudouridines and other RNA modifications using HydraPsiSeq protocol

Detection of RNA modified nucleotides using deep sequencing can be performed by several approaches, including antibody-driven enrichment and natural or chemically induced RT signatures. However, only very few RNA modified nucleotides generate natural RT signatures and antibody-driven enrichment heavily depends on the quality of antibodies used and may be highly biased. Thus, the use of chemically-induced RT signatures is now considered as the most trusted experimental approach. In addition, the use of chemical reagents allows inclusion of simple "mock-treated" controls, to exclude spontaneous RT arrests, SNPs and other misincorporation-prone sites. Hydrazine is a well-known RNA-specific reagent, already extensively used in the past for RNA sequencing and structural probing. Hydrazine is highly reactive to U and shows low reaction rates with ψ residues, allowing their distinction by deep sequencing-based protocols. However, other modified RNA residues also show particular behavior upon hydrazine treatment. Here we present methodological developments allowing to use HydraPsiSeq for precise quantification of RNA pseudouridylation and also detection and quantification of some other RNA modifications, in addition to ψ residues.

Integrated analyses of m1A regulator-mediated modification patterns in tumor microenvironment-infiltrating immune cells in colon cancer

Emerging evidence has revealed the crucial role of transcriptional RNA methyladenosine modification in immune response. However, the potential role of RNA N1-methyladenosine (m¹A) modification of immune cells in the tumor microenvironment (TME) still remains unclear. In this study, we identified three distinct m¹A modification patterns based on the integrated analyses of nine m¹A regulators, which are significantly related to Relapse-free survival (RFS), Overall survival (OS), and TME infiltration cells in colon cancer patients. Furthermore, the m1AScore was generated by using principal components analysis (PCA) of expression of the 71 m¹A-related genes to further demonstrate the characteristics of m¹A patterns in colon cancer. In summary, a low m1AScore could be characterized by lower EMT, pan-F TBRS, and TNM stages, as well as less presence of lymphatic invasion, and, hence, good prognosis. At the same time, a low m1AScore could also be linked to CD8 + T effector proliferation, in addition to high microsatellite instability (MSI), neoantigen burden and PD-L1 expression, showing prolonged survival and better response after undergoing an anti-PD-L1 immunotherapy regimen in the public immunotherapy cohort. Our work reveals that m¹A modification patterns play a key role in the formation of TME complexity and diversity in the context of immune cell infiltration. Accordingly, this m1AScore system provides an efficient method by which to identify and characterize TME immune cell infiltration, thereby allowing for more personalized and effective antitumor immunotherapy strategies.

Mapping messenger RNA methylations at single base resolution

The messenger RNA (mRNA) methylations in mammalian cells have been found to contain N⁶-methyladenosine (m⁶A), N⁶-2'-O-dimethyladenosine (m⁶Am), 7-methylguanosine (m⁷G), 1-methyladenosine (m¹A), 5-methylcytosine (m⁵C), and 2'-O-methylation (2'-OMe). Their regulatory functions in control of mRNA fate and gene expression are being increasingly uncovered. To unambiguously understand the critical roles of mRNA methylations in physiological and pathological processes, mapping these methylations at single base resolution is highly required. Here, we will review the progresses made in methylation sequencing methodologies developed mainly in recent two years, with an emphasis on chemical labeling-assisted single base resolution methods, and discuss the problems and prospects as well.

The Regulation of RNA Modification Systems: The Next Frontier in Epitranscriptomics?

RNA modifications, long considered to be molecular curiosities embellishing just abundant and non-coding RNAs, have now moved into the focus of both academic and applied research. Dedicated research efforts (epitranscriptomics) aim at deciphering the underlying principles by determining RNA modification landscapes and investigating the molecular mechanisms that establish, interpret and modulate the information potential of RNA beyond the combination of four canonical nucleotides. This has resulted in mapping various epitranscriptomes at high resolution and in cataloguing the effects caused by aberrant RNA modification circuitry. While the scope of the obtained insights has been complex and exciting, most of current epitranscriptomics appears to be stuck in the process of producing data, with very few efforts to disentangle cause from consequence when studying a specific RNA modification system. This article discusses various knowledge gaps in this field with the aim to raise one specific question: how are the enzymes regulated that dynamically install and modify RNA modifications? Furthermore, various technologies will be highlighted whose development and use might allow identifying specific and context-dependent regulators of epitranscriptomic mechanisms. Given the complexity of individual epitranscriptomes, determining their regulatory principles will become crucially important, especially when aiming at modifying specific aspects of an epitranscriptome both for experimental and, potentially, therapeutic purposes.

Analysis of RNA Modifications by Second- and Third-Generation Deep Sequencing: 2020 Update

The precise mapping and quantification of the numerous RNA modifications that are present in tRNAs, rRNAs, ncRNAs/miRNAs, and mRNAs remain a major challenge and a top priority of the epitranscriptomics field. After the keystone discoveries of massive m⁶A methylation in mRNAs, dozens of deep sequencing-based methods and protocols were proposed for the analysis of various RNA modifications, allowing us to considerably extend the list of detectable modified residues. Many of the currently used methods rely on the particular reverse transcription signatures left by RNA modifications in cDNA; these signatures may be naturally present or induced by an appropriate enzymatic or chemical treatment. The newest approaches also include labeling at RNA abasic sites that result from the selective removal of RNA modification or the enhanced cleavage of the RNA ribose-phosphate chain (perhaps also protection from cleavage), followed by specific adapter ligation. Classical affinity/immunoprecipitation-based protocols use either antibodies against modified RNA bases or proteins/enzymes, recognizing RNA modifications. In this survey, we review the most recent achievements in this highly dynamic field, including promising attempts to map RNA modifications by the direct single-molecule sequencing of RNA by nanopores.

Deciphering RNA modifications at base resolution: from chemistry to biology

Nearly 200 distinct chemical modifications of RNAs have been discovered to date. Their analysis via direct methods has been possible in abundant RNA species, such as ribosomal, transfer or viral RNA, since several decades. However, their analysis in less abundant RNAs species, especially cellular messenger RNAs, was rendered possible only recently with the advent of high throughput sequencing techniques. Given the growing biomedical interest of the proteins that write, erase and read RNA modifications, ingenious new methods to enrich and identify RNA modifications at base resolution have been implemented, and more efforts are underway to render them more quantitative. Here, we review several crucial modification-specific (bio)chemical approaches and discuss their advantages and shortcomings for exploring the epitranscriptome.

Effects of mRNA Modifications on Translation: An Overview

The mRNA epitranscriptome imparts diversity to gene expression by installing chemical modifications. Advances in detection methods have identified chemical modifications in eukaryotic, bacterial, and viral messenger RNAs (mRNAs). The biological functions of modifications in mRNAs still remain to be understood. Chemical modifications are introduced in synthetic mRNAs meant for therapeutic applications to maximize expression from the synthetic mRNAs and to evade the host immune response. This overview provides a background of chemical modifications found in mRNAs, with an emphasis on pseudouridine and its known effects on the mRNA life cycle, its potential applications in synthetic mRNA, and the methods used to assess its effects on mRNA translation.

New Twists in Detecting mRNA Modification Dynamics

Modified nucleotides in mRNA are an essential addition to the standard genetic code of four nucleotides in animals, plants, and their viruses. The emerging field of epitranscriptomics examines nucleotide modifications in mRNA and their impact on gene expression. The low abundance of nucleotide modifications and technical limitations, however, have hampered systematic analysis of their occurrence and functions. Selective chemical and immunological identification of modified nucleotides has revealed global candidate topology maps for many modifications in mRNA, but further technical advances to increase confidence will be necessary. Single-molecule sequencing introduced by Oxford Nanopore now promises to overcome such limitations, and we summarize current progress with a particular focus on the bioinformatic challenges of this novel sequencing technology.

MicroRNA-214 Inhibits Chicken Myoblasts Proliferation, Promotes Their Differentiation, and Targets the TRMT61A Gene

The proliferation and differentiation of myoblasts is an important process of skeletal muscle development. In this process, microRNAs (miRNAs) play an important role in the proliferation and differentiation of chicken primary myoblasts (CPMs). Our previous study found that miR-214 and the tRNA methyltransferase 61A (TRMT61A) gene were differentially expressed in different stages of proliferation and differentiation. Therefore, this study aimed to explore the effect of miR-214 on the proliferation and differentiation of CPMs and the functional relationship between miR-214 and TRMT61A. In this study, we detected the effect of miR-214 on the proliferation of CPMs by qPCR, flow cytometry, CCK-8, and EdU after the overexpression and interference of miR-214. qPCR, Western blotting, and indirect immunofluorescence were used to detect the effect of miR-214 on the differentiation of the CPMs. The expression patterns of miR-214 and TRMT61A were observed at different time points of differentiation induced by the CPMs. The results show that miR-214 inhibited the proliferation of the CPMs and promoted the differentiation of the CPMs. The Dual-Luciferase Reporter assay and the expression pattern of miR-214 and TRMT61A suggested that they had a negative regulatory target relationship. This study revealed the function and regulatory mechanism of miR-214 in the proliferation and differentiation of CPMs.

The epitranscriptome beyond m6A

Following its transcription, RNA can be modified by >170 chemically distinct types of modifications - the epitranscriptome. In recent years, there have been substantial efforts to uncover and characterize the modifications present on mRNA, motivated by the potential of such modifications to regulate mRNA fate and by discoveries and advances in our understanding of N ⁶-methyladenosine (m⁶A). Here, we review our knowledge regarding the detection, distribution, abundance, biogenesis, functions and possible mechanisms of action of six of these modifications - pseudouridine (Ψ), 5-methylcytidine (m⁵C), N ¹-methyladenosine (m¹A), N ⁴-acetylcytidine (ac⁴C), ribose methylations (N_m) and N ⁷-methylguanosine (m⁷G). We discuss the technical and analytical aspects that have led to inconsistent conclusions and controversies regarding the abundance and distribution of some of these modifications. We further highlight shared commonalities and important ways in which these modifications differ with respect to m⁶A, based on which we speculate on their origin and their ability to acquire functions over evolutionary timescales.

HydraPsiSeq: a method for systematic and quantitative mapping of pseudouridines in RNA

Developing methods for accurate detection of RNA modifications remains a major challenge in epitranscriptomics. Next-generation sequencing-based mapping approaches have recently emerged but, often, they are not quantitative and lack specificity. Pseudouridine (ψ), produced by uridine isomerization, is one of the most abundant RNA modification. ψ mapping classically involves derivatization with soluble carbodiimide (CMCT), which is prone to variation making this approach only semi-quantitative. Here, we developed 'HydraPsiSeq', a novel quantitative ψ mapping technique relying on specific protection from hydrazine/aniline cleavage. HydraPsiSeq is quantitative because the obtained signal directly reflects pseudouridine level. Furthermore, normalization to natural unmodified RNA and/or to synthetic in vitro transcripts allows absolute measurements of modification levels. HydraPsiSeq requires minute amounts of RNA (as low as 10-50 ng), making it compatible with high-throughput profiling of diverse biological and clinical samples. Exploring the potential of HydraPsiSeq, we profiled human rRNAs, revealing strong variations in pseudouridylation levels at ∼20-25 positions out of total 104 sites. We also observed the dynamics of rRNA pseudouridylation throughout chondrogenic differentiation of human bone marrow stem cells. In conclusion, HydraPsiSeq is a robust approach for the systematic mapping and accurate quantification of pseudouridines in RNAs with applications in disease, aging, development, differentiation and/or stress response.

The epitranscriptome in stem cell biology and neural development

The blossoming field of epitranscriptomics has recently garnered attention across many fields by findings that chemical modifications on RNA have immense biological consequences. Methylation of nucleotides in RNA, including N6-methyladenosine (m⁶A), 2-O-dimethyladenosine (m⁶A_m), N1-methyladenosine (m¹A), 5-methylcytosine (m⁵C), and isomerization of uracil to pseudouridine (Ψ), have the potential to alter RNA processing events and contribute to developmental processes and different diseases. Though the abundance and roles of some RNA modifications remain contentious, the epitranscriptome is thought to be especially relevant in stem cell biology and neurobiology. In particular, m⁶A occurs at the highest levels in the brain and plays major roles in embryonic stem cell differentiation, brain development, and neurodevelopmental disorders. However, studies in these areas have reported conflicting results on epitranscriptomic regulation of stem cell pluripotency and mechanisms in neural development. In this review we provide an overview of the current understanding of several RNA modifications and disentangle the various findings on epitranscriptomic regulation of stem cell biology and neural development.

Naturally occurring modified ribonucleosides

The chemical identity of RNA molecules beyond the four standard ribonucleosides has fascinated scientists since pseudouridine was characterized as the "fifth" ribonucleotide in 1951. Since then, the ever-increasing number and complexity of modified ribonucleosides have been found in viruses and throughout all three domains of life. Such modifications can be as simple as methylations, hydroxylations, or thiolations, complex as ring closures, glycosylations, acylations, or aminoacylations, or unusual as the incorporation of selenium. While initially found in transfer and ribosomal RNAs, modifications also exist in messenger RNAs and noncoding RNAs. Modifications have profound cellular outcomes at various levels, such as altering RNA structure or being essential for cell survival or organism viability. The aberrant presence or absence of RNA modifications can lead to human disease, ranging from cancer to various metabolic and developmental illnesses such as Hoyeraal-Hreidarsson syndrome, Bowen-Conradi syndrome, or Williams-Beuren syndrome. In this review article, we summarize the characterization of all 143 currently known modified ribonucleosides by describing their taxonomic distributions, the enzymes that generate the modifications, and any implications in cellular processes, RNA structure, and disease. We also highlight areas of active research, such as specific RNAs that contain a particular type of modification as well as methodologies used to identify novel RNA modifications. This article is categorized under: RNA Processing > RNA Editing and Modification.

A Census and Categorization Method of Epitranscriptomic Marks

In the past few years, thorough investigation of chemical modifications operated in the cells on ribonucleic acid (RNA) molecules is gaining momentum. This new field of research has been dubbed “epitranscriptomics”, in analogy to best-known epigenomics, to stress the potential of ensembles of RNA modifications to constitute a post-transcriptional regulatory layer of gene expression orchestrated by writer, reader, and eraser RNA-binding proteins (RBPs). In fact, epitranscriptomics aims at identifying and characterizing all functionally relevant changes involving both non-substitutional chemical modifications and editing events made to the transcriptome. Indeed, several types of RNA modifications that impact gene expression have been reported so far in different species of cellular RNAs, including ribosomal RNAs, transfer RNAs, small nuclear RNAs, messenger RNAs, and long non-coding RNAs. Supporting functional relevance of this largely unknown regulatory mechanism, several human diseases have been associated directly to RNA modifications or to RBPs that may play as effectors of epitranscriptomic marks. However, an exhaustive epitranscriptome’s characterization, aimed to systematically classify all RNA modifications and clarify rules, actors, and outcomes of this promising regulatory code, is currently not available, mainly hampered by lack of suitable detecting technologies. This is an unfortunate limitation that, thanks to an unprecedented pace of technological advancements especially in the sequencing technology field, is likely to be overcome soon. Here, we review the current knowledge on epitranscriptomic marks and propose a categorization method based on the reference ribonucleotide and its rounds of modifications (“stages”) until reaching the given modified form. We believe that this classification scheme can be useful to coherently organize the expanding number of discovered RNA modifications.

Cytoplasmic m1A reader YTHDF3 inhibits trophoblast invasion by downregulation of m1A-methylated IGF1R

N1-methyladenosine (m1A) is one of the important post-transcriptional modifications in RNA and plays an important role in promoting translation or decay of m1A-methylated messenger RNA (mRNA), but the "reader" protein and the exact biological role of m1A remain to be determined. Here, we identified that nine potential m1A "reader" proteins including YTH domain family and heterogeneous nuclear ribonucleoprotein by mass spectrometry, and among them, YTH domain-containing protein 3 (YTHDF3), could bind directly to m1A-carrying RNA. YTHDF3 was then identified to negatively regulate invasion and migration of trophoblast. Mechanistically, we found that the m1A "reader" YTHDF3 bound to certain m1A-methylated transcripts, such as insulin-like growth factor 1 receptor (IGF1R), with the combination of iCLIP-seq (individual-nucleotide resolution ultraviolet crosslinking and immunoprecipitation high-throughput sequencing) and m1A-seq. Furthermore, YTHDF3 could promote IGF1R mRNA degradation and thus inhibit IGF1R protein expression along with its downstream matrix metallopeptidase 9 signaling pathway, consequently decreasing migration and invasion of trophoblast. Thus, we demonstrated that YTHDF3 as an m1A reader decreased invasion and migration of trophoblast by inhibiting IGF1R expression. Our study outlines a new m1A epigenetic way to regulate the trophoblast activity, which suggests a novel therapeutic target for trophoblast-associated pregnancy disorders.

A Benzophenone-Based Photocaging Strategy for the N7 Position of Guanosine

Selective modification of nucleobases with photolabile caging groups enables the study and control of processes and interactions of nucleic acids. Numerous positions on nucleobases have been targeted, but all involve formal substitution of a hydrogen atom with a photocaging group. Nature, however, also uses ring-nitrogen methylation, such as m7 G and m1 A, to change the electronic structure and properties of RNA and control biomolecular interactions essential for translation and turnover. We report that aryl ketones such as benzophenone and α-hydroxyalkyl ketone are photolabile caging groups if installed at the N7 position of guanosine or the N1 position of adenosine. Common photocaging groups derived from the ortho-nitrobenzyl moiety were not suitable. Both chemical and enzymatic methods for site-specific modification of N7G in nucleosides, dinucleotides, and RNA were developed, thereby opening the door to studying the molecular interactions of m7 G and m1 A with spatiotemporal control.

Holistic Optimization of Bioinformatic Analysis Pipeline for Detection and Quantification of 2'-O-Methylations in RNA by RiboMethSeq

A major trend in the epitranscriptomics field over the last 5 years has been the high-throughput analysis of RNA modifications by a combination of specific chemical treatment(s), followed by library preparation and deep sequencing. Multiple protocols have been described for several important RNA modifications, such as 5-methylcytosine (m⁵C), pseudouridine (ψ), 1-methyladenosine (m¹A), and 2′-O-methylation (Nm). One commonly used method is the alkaline cleavage-based RiboMethSeq protocol, where positions of reads' 5'-ends are used to distinguish nucleotides protected by ribose methylation. This method was successfully applied to detect and quantify Nm residues in various RNA species such as rRNA, tRNA, and snRNA. Such applications require adaptation of the initially published protocol(s), both at the wet bench and in the bioinformatics analysis. In this manuscript, we describe the optimization of RiboMethSeq bioinformatics at the level of initial read treatment, alignment to the reference sequence, counting the 5′- and 3′- ends, and calculation of the RiboMethSeq scores, allowing precise detection and quantification of the Nm-related signal. These improvements introduced in the original pipeline permit a more accurate detection of Nm candidates and a more precise quantification of Nm level variations. Applications of the improved RiboMethSeq treatment pipeline for different cellular RNA types are discussed.

Pan-cancer analysis of RNA methyltransferases identifies FTSJ3 as a potential regulator of breast cancer progression

RNA methylation, catalysed by a set of RNA methyltransferases (RNMTs), modulates RNA structures, properties, and biological functions. RNMTs are increasingly documented to be dysregulated in various human diseases, particularly developmental disorders and cancer. However, the genomic and transcriptomic alterations of RNMTs, as well as their functional roles in human cancer, are limited. In this study, we utilized an unbiased approach to examine copy number alterations and mutation rates of 58 RNMTs in more than 10,000 clinical samples across 32 human cancer types. We also investigated these alterations and RNMT expression level as they related to clinical features such as tumour subtype, grade, and survival in a large cohort of tumour samples, focusing on breast cancer. Loss-of-function analysis was performed to examine RNMT candidates with important roles in growth and viability of breast cancer cells. We identified a subset of RNMTs, notably TRMT12, NSUN2, TARBP1, and FTSJ3, that were amplified or mutated in a subset of human cancers. Several RNMTs were significantly associated with breast cancer aggressiveness and poor prognosis. Loss-of-function analysis indicated FTSJ3, a 2'-O-Me methyltransferase, as a candidate RNMT with functional roles in promoting cancer growth and survival. A subset of RNMTs, like FTSJ3, represents promising novel targets for anticancer drug discovery. Our findings provide a framework for further study of the functional consequences of RNMT alterations in human cancer and for developing therapies that target cancer-promoting RNMTs in the future.

YTHDF2 Recognition of N1-Methyladenosine (m1A)-Modified RNA Is Associated with Transcript Destabilization

Epitranscriptomic modifications play an important role in RNA function and can impact gene expression. Here, we apply a chemical proteomics approach to investigate readers of N¹-methyladenosine (m¹A), a poorly characterized modification on mammalian mRNA. We find that YTHDF proteins, known m⁶A readers, recognize m¹A-modified sequences in a methylation-specific manner. We characterize binding of recombinant YTHDF1/2 proteins to m¹A-modified oligonucleotides to demonstrate that these interactions can exhibit comparable affinity to m⁶A-recognition events and occur in diverse sequence contexts. Further, we demonstrate YTHDF2 interacts specifically with endogenously modified m¹A transcripts. Finally, we deplete cellular YTHDF2 to show that the abundance of m¹A-modified transcripts is increased in its absence. Similarly, increasing m¹A levels through depletion of ALKBH3, an m¹A eraser protein, destabilizes known m¹A-containing RNAs. Our results shed light on the function of m¹A on mRNA and provide a mechanistic framework to further evaluate the role of m¹A in biological processes.

Limited antibody specificity compromises epitranscriptomic analyses

A controversial discussion on the occurrence of the RNA modification m¹A in mRNA takes a new turn, as an antibody with a central role in modification mapping was shown to also bind mRNA cap structures.

Antibody cross-reactivity accounts for widespread appearance of m1A in 5'UTRs

N1-methyladenosine (m1A) was proposed to be a highly prevalent modification in mRNA 5'UTRs based on mapping studies using an m1A-binding antibody. We developed a bioinformatic approach to discover m1A and other modifications in mRNA throughout the transcriptome by analyzing preexisting ultra-deep RNA-Seq data for modification-induced misincorporations. Using this approach, we detected appreciable levels of m1A only in one mRNA: the mitochondrial MT-ND5 transcript. As an alternative approach, we also developed an antibody-based m1A-mapping approach to detect m1A at single-nucleotide resolution, and confirmed that the commonly used m1A antibody maps sites to the transcription-start site in mRNA 5'UTRs. However, further analysis revealed that these were false-positives caused by binding of the antibody to the m7G-cap. A different m1A antibody that lacks cap-binding cross-reactivity does not show enriched binding in 5'UTRs. These results demonstrate that high-stoichiometry m1A sites are exceedingly rare in mRNAs and that previous mappings of m1A to 5'UTRs were the result of antibody cross-reactivity to the 5' cap.

Comprehensive review and assessment of computational methods for predicting RNA post-transcriptional modification sites from RNA sequences

RNA post-transcriptional modifications play a crucial role in a myriad of biological processes and cellular functions. To date, more than 160 RNA modifications have been discovered; therefore, accurate identification of RNA-modification sites is fundamental for a better understanding of RNA-mediated biological functions and mechanisms. However, due to limitations in experimental methods, systematic identification of different types of RNA-modification sites remains a major challenge. Recently, more than 20 computational methods have been developed to identify RNA-modification sites in tandem with high-throughput experimental methods, with most of these capable of predicting only single types of RNA-modification sites. These methods show high diversity in their dataset size, data quality, core algorithms, features extracted and feature selection techniques and evaluation strategies. Therefore, there is an urgent need to revisit these methods and summarize their methodologies, in order to improve and further develop computational techniques to identify and characterize RNA-modification sites from the large amounts of sequence data. With this goal in mind, first, we provide a comprehensive survey on a large collection of 27 state-of-the-art approaches for predicting N1-methyladenosine and N6-methyladenosine sites. We cover a variety of important aspects that are crucial for the development of successful predictors, including the dataset quality, operating algorithms, sequence and genomic features, feature selection, model performance evaluation and software utility. In addition, we also provide our thoughts on potential strategies to improve the model performance. Second, we propose a computational approach called DeepPromise based on deep learning techniques for simultaneous prediction of N1-methyladenosine and N6-methyladenosine. To extract the sequence context surrounding the modification sites, three feature encodings, including enhanced nucleic acid composition, one-hot encoding, and RNA embedding, were used as the input to seven consecutive layers of convolutional neural networks (CNNs), respectively. Moreover, DeepPromise further combined the prediction score of the CNN-based models and achieved around 43% higher area under receiver-operating curve (AUROC) for m1A site prediction and 2–6% higher AUROC for m6A site prediction, respectively, when compared with several existing state-of-the-art approaches on the independent test. In-depth analyses of characteristic sequence motifs identified from the convolution-layer filters indicated that nucleotide presentation at proximal positions surrounding the modification sites contributed most to the classification, whereas those at distal positions also affected classification but to different extents. To maximize user convenience, a web server was developed as an implementation of DeepPromise and made publicly available at http://DeepPromise.erc.monash.edu/, with the server accepting both RNA sequences and genomic sequences to allow prediction of two types of putative RNA-modification sites.

Dynamic methylome of internal mRNA N7-methylguanosine and its regulatory role in translation

Over 150 types of RNA modifications are identified in RNA molecules. Transcriptome profiling is one of the key steps in decoding the epitranscriptomic panorama of these chemical modifications and their potential functions. N⁷-methylguanosine (m⁷G) is one of the most abundant modifications present in tRNA, rRNA and mRNA 5'cap, and has critical roles in regulating RNA processing, metabolism and function. Besides its presence at the cap position in mRNAs, m⁷G is also identified in internal mRNA regions. However, its transcriptome-wide distribution and dynamic regulation within internal mRNA regions remain unknown. Here, we have established m⁷G individual-nucleotide-resolution cross-linking and immunoprecipitation with sequencing (m⁷G miCLIP-seq) to specifically detect internal mRNA m⁷G modification. Using this approach, we revealed that m⁷G is enriched at the 5'UTR region and AG-rich contexts, a feature that is well-conserved across different human/mouse cell lines and mouse tissues. Strikingly, the internal m⁷G modification is dynamically regulated under both H₂O₂ and heat shock treatments, with remarkable accumulations in the CDS and 3'UTR regions, and functions in promoting mRNA translation efficiency. Consistently, a PCNA 3'UTR minigene reporter harboring the native m⁷G modification site displays both enriched m⁷G modification and increased mRNA translation upon H₂O₂ treatment compared to the m⁷G site-mutated minigene reporter (G to A). Taken together, our findings unravel the dynamic profiles of internal mRNA m⁷G methylome and highlight m⁷G as a novel epitranscriptomic marker with regulatory roles in translation.

m1A Regulated Genes Modulate PI3K/AKT/mTOR and ErbB Pathways in Gastrointestinal Cancer

BACKGROUND: Gene expression can be posttranscriptionally regulated by a complex network of proteins. N1-methyladenosine (m1A) is a newly validated RNA modification. However, little is known about both its influence and biogenesis in tumor development. METHODS: This study analyzed TCGA data of patients with five kinds of gastrointestinal (GI) cancers. Using data from cBioPortal, molecular features of the nine known m1A-related enzymes in GI cancers were investigated. Using a variety of bioinformatics approach, the impact of m1A regulators on its downstream signaling pathway was studied. To further confirm this regulation, the effect of m1A writer ALKBH3 knockdown was studied using RNA-seq data from published database. RESULTS: Dysregulation and multiple types of genetic alteration of putative m1A-related enzymes in tumor samples were observed. The ErbB and mTOR pathways with ErbB2, mTOR, and AKT1S1 hub genes were identified as being regulated by m1A-related enzymes. The expression of both ErbB2 and AKT1S1 was decreased after m1A writer ALKBH3 knockdown. Furthermore, Gene Ontology analysis revealed that m1A downstream genes were associated with cell proliferation, and the results showed that m1A genes are reliably linked to mTOR. CONCLUSION: This study demonstrated for the first time the dysregulation of m1A regulators in GI cancer and its signaling pathways and will contribute to the understanding of RNA modification in cancer.

The mammalian mitochondrial epitranscriptome

Correct expression of the mitochondrially-encoded genes is critical for the production of the components of the oxidative phosphorylation machinery. Post-transcriptional modifications of mitochondrial transcripts have been emerging as an important regulatory feature of mitochondrial gene expression. Here we review the current knowledge on how the mammalian mitochondrial epitranscriptome participates in regulating mitochondrial homeostasis. In particular, we focus on the latest breakthroughs made towards understanding the roles of the modified nucleotides in mitochondrially-encoded ribosomal and transfer RNAs, the enzymes responsible for introducing these modifications and on recent transcriptome-wide studies reporting modifications to mitochondrial messenger RNAs. This article is part of a Special Issue entitled: mRNA modifications in gene expression control edited by Dr. Matthias Soller and Dr. Rupert Fray.

Eukaryotic Translation Elongation is Modulated by Single Natural Nucleotide Derivatives in the Coding Sequences of mRNAs

RNA modifications are crucial factors for efficient protein synthesis. All classes of RNAs that are involved in translation are modified to different extents. Recently, mRNA modifications and their impact on gene regulation became a focus of interest because they can exert a variety of effects on the fate of mRNAs. mRNA modifications within coding sequences can either directly or indirectly interfere with protein synthesis. In order to investigate the roles of various natural occurring modified nucleotides, we site-specifically introduced them into the coding sequence of reporter mRNAs and subsequently translated them in HEK293T cells. The analysis of the respective protein products revealed a strong position-dependent impact of RNA modifications on translation efficiency and accuracy. Whereas a single 5-methylcytosine (m⁵C) or pseudouridine (Ψ) did not reduce product yields, N¹-methyladenosine (m¹A) generally impeded the translation of the respective modified mRNA. An inhibitory effect of 2′O-methlyated nucleotides (Nm) and N⁶-methyladenosine (m⁶A) was strongly dependent on their position within the codon. Finally, we could not attribute any miscoding potential to the set of mRNA modifications tested in HEK293T cells.

Methods for RNA Modification Mapping Using Deep Sequencing: Established and New Emerging Technologies

New analytics of post-transcriptional RNA modifications have paved the way for a tremendous upswing of the biological and biomedical research in this field. This especially applies to methods that included RNA-Seq techniques, and which typically result in what is termed global scale modification mapping. In this process, positions inside a cell’s transcriptome are receiving a status of potential modification sites (so called modification calling), typically based on a score of some kind that issues from the particular method applied. The resulting data are thought to represent information that goes beyond what is contained in typical transcriptome data, and hence the field has taken to use the term “epitranscriptome”. Due to the high rate of newly published mapping techniques, a significant number of chemically distinct RNA modifications have become amenable to mapping, albeit with variegated accuracy and precision, depending on the nature of the technique. This review gives a brief overview of known techniques, and how they were applied to modification calling.

Misincorporation signatures for detecting modifications in mRNA: Not as simple as it sounds

Post-transcriptional modification on mRNA has become a field of intense interest in recent years, and next-generation sequencing based technologies are constantly emerging to detect an increasing number of modifications at a transcriptome-wide level. Some of these approaches are based on identification of misincorporation events induced by reverse transcriptase at modified sites. Although conceptually trivial, sensitive and specific identification of such events is a challenge prone to a surprising number of artifacts, which can lead to substantially inflated estimates of the abundance of diverse modifications. Here we discuss the sources of some of these artifacts and delineate approaches to overcome them.

Epitranscriptomics: mapping methods and beyond

The field of epigenetics has changed the way we look at both DNA and proteins. Could the same be true for RNA? Joseph Martin discusses.

[Formula: see text]