The epitranscriptome toolbox.

In the last decade, the notion that mRNA modifications are involved in regulation of gene expression was demonstrated in thousands of studies. To date, new technologies and methods allow accurate identification, transcriptome-wide mapping, and functional characterization of a growing number of RNA modifications, providing important insights into the biology of these marks. Most of the methods and approaches were developed for studying m6A, the most prevalent internal mRNA modification. However, unique properties of other RNA modifications stimulated the development of additional approaches. In this technical primer, we will discuss the available tools and approaches for detecting and studying different RNA modifications.

Temporal Control of Mammalian Cortical Neurogenesis by m6A Methylation.

N6-methyladenosine (m6A), installed by the Mettl3/Mettl14 methyltransferase complex, is the most prevalent internal mRNA modification. Whether m6A regulates mammalian brain development is unknown. Here, we show that m6A depletion by Mettl14 knockout in embryonic mouse brains prolongs the cell cycle of radial glia cells and extends cortical neurogenesis into postnatal stages. m6A depletion by Mettl3 knockdown also leads to a prolonged cell cycle and maintenance of radial glia cells. m6A sequencing of embryonic mouse cortex reveals enrichment of mRNAs related to transcription factors, neurogenesis, the cell cycle, and neuronal differentiation, and m6A tagging promotes their decay. Further analysis uncovers previously unappreciated transcriptional prepatterning in cortical neural stem cells. m6A signaling also regulates human cortical neurogenesis in forebrain organoids. Comparison of m6A-mRNA landscapes between mouse and human cortical neurogenesis reveals enrichment of human-specific m6A tagging of transcripts related to brain-disorder risk genes. Our study identifies an epitranscriptomic mechanism in heightened transcriptional coordination during mammalian cortical neurogenesis.

2'-O-methylation at internal sites on mRNA promotes mRNA stability.

2'-O-methylation (Nm) is a prominent RNA modification well known in noncoding RNAs and more recently also found at many mRNA internal sites. However, their function and base-resolution stoichiometry remain underexplored. Here, we investigate the transcriptome-wide effect of internal site Nm on mRNA stability. Combining nanopore sequencing with our developed machine learning method, NanoNm, we identify thousands of Nm sites on mRNAs with a single-base resolution. We observe a positive effect of FBL-mediated Nm modification on mRNA stability and expression level. Elevated FBL expression in cancer cells is associated with increased expression levels for 2'-O-methylated mRNAs of cancer pathways, implying the role of FBL in post-transcriptional regulation. Lastly, we find that FBL-mediated 2'-O-methylation connects to widespread 3' UTR shortening, a mechanism that globally increases RNA stability. Collectively, we demonstrate that FBL-mediated Nm modifications at mRNA internal sites regulate gene expression by enhancing mRNA stability.

Base-resolution m5C profiling across the mammalian transcriptome by bisulfite-free enzyme-assisted chemical labeling approach.

5-methylcytosine (m5C) is a prevalent RNA modification crucial for gene expression regulation. However, accurate and sensitive m5C sites identification remains challenging due to severe RNA degradation and reduced sequence complexity during bisulfite sequencing (BS-seq). Here, we report m5C-TAC-seq, a bisulfite-free approach combining TET-assisted m5C-to-f5C oxidation with selective chemical labeling, therefore enabling direct base-resolution m5C detection through pre-enrichment and C-to-T transitions at m5C sites. With m5C-TAC-seq, we comprehensively profiled the m5C methylomes in human and mouse cells, identifying a substantially larger number of confident m5C sites. Through perturbing potential m5C methyltransferases, we deciphered the responsible enzymes for most m5C sites, including the characterization of NSUN5's involvement in mRNA m5C deposition. Additionally, we characterized m5C dynamics during mESC differentiation. Notably, the mild reaction conditions and preservation of nucleotide composition in m5C-TAC-seq allow m5C detection in chromatin-associated RNAs. The accurate and robust m5C-TAC-seq will advance research into m5C methylation functional investigation.

Single-nucleus multiomic mapping of m6A methylomes and transcriptomes in native populations of cells with sn-m6A-CT.

N6-methyladenosine (m6A) RNA modification plays important roles in the governance of gene expression and is temporally regulated in different cell states. In contrast to global m6A profiling in bulk sequencing, single-cell technologies for analyzing m6A heterogeneity are not extensively established. Here, we developed single-nucleus m6A-CUT&Tag (sn-m6A-CT) for simultaneous profiling of m6A methylomes and transcriptomes within a single nucleus using mouse embryonic stem cells (mESCs). m6A-CT is capable of enriching m6A-marked RNA molecules in situ, without isolating RNAs from cells. We adapted m6A-CT to the droplet-based single-cell omics platform and demonstrated high-throughput performance in analyzing nuclei isolated from thousands of cells from various cell types. We show that sn-m6A-CT profiling is sufficient to determine cell identity and allows the generation of cell-type-specific m6A methylome landscapes from heterogeneous populations. These indicate that sn-m6A-CT provides additional dimensions to multimodal datasets and insights into epitranscriptomic landscape in defining cell fate identity and states.

SRSF2 plays an unexpected role as reader of m5C on mRNA, linking epitranscriptomics to cancer.

A common mRNA modification is 5-methylcytosine (m5C), whose role in gene-transcript processing and cancer remains unclear. Here, we identify serine/arginine-rich splicing factor 2 (SRSF2) as a reader of m5C and impaired SRSF2 m5C binding as a potential contributor to leukemogenesis. Structurally, we identify residues involved in m5C recognition and the impact of the prevalent leukemia-associated mutation SRSF2P95H. We show that SRSF2 binding and m5C colocalize within transcripts. Furthermore, knocking down the m5C writer NSUN2 decreases mRNA m5C, reduces SRSF2 binding, and alters RNA splicing. We also show that the SRSF2P95H mutation impairs the ability of the protein to read m5C-marked mRNA, notably reducing its binding to key leukemia-related transcripts in leukemic cells. In leukemia patients, low NSUN2 expression leads to mRNA m5C hypomethylation and, combined with SRSF2P95H, predicts poor outcomes. Altogether, we highlight an unrecognized mechanistic link between epitranscriptomics and a key oncogenesis driver.

RNA m6A methylation across the transcriptome.

Since the early days of foundational studies of nucleic acids, many chemical moieties have been discovered to decorate RNA and DNA in diverse organisms. In mammalian cells, one of these chemical modifications, N6-methyl adenosine (m6A), is unique in a way that it is highly abundant not only on RNA polymerase II (RNAPII) transcribed, protein-coding transcripts but also on non-coding RNAs, such as ribosomal RNAs and snRNAs, mediated by distinct, evolutionarily conserved enzymes. Here, we review RNA m6A modification in the light of the recent appreciation of nuclear roles for m6A in regulating chromatin states and gene expression, as well as the recent discoveries of the evolutionarily conserved methyltransferases, which catalyze methylation of adenosine on diverse sets of RNAs. Considering that the substrates of these enzymes are involved in many important biological processes, this modification warrants further research to understand the molecular mechanisms and functions of m6A in health and disease.

Transcription-wide mapping of dihydrouridine reveals that mRNA dihydrouridylation is required for meiotic chromosome segregation.

The epitranscriptome has emerged as a new fundamental layer of control of gene expression. Nevertheless, the determination of the transcriptome-wide occupancy and function of RNA modifications remains challenging. Here we have developed Rho-seq, an integrated pipeline detecting a range of modifications through differential modification-dependent rhodamine labeling. Using Rho-seq, we confirm that the reduction of uridine to dihydrouridine (D) by the Dus reductase enzymes targets tRNAs in E. coli and fission yeast. We find that the D modification is also present on fission yeast mRNAs, particularly those encoding cytoskeleton-related proteins, which is supported by large-scale proteome analyses and ribosome profiling. We show that the α-tubulin encoding mRNA nda2 undergoes Dus3-dependent dihydrouridylation, which affects its translation. The absence of the modification on nda2 mRNA strongly impacts meiotic chromosome segregation, resulting in low gamete viability. Applying Rho-seq to human cells revealed that tubulin mRNA dihydrouridylation is evolutionarily conserved.

The Mettl3 epitranscriptomic writer amplifies p53 stress responses.

The p53 transcription factor drives anti-proliferative gene expression programs in response to diverse stressors, including DNA damage and oncogenic signaling. Here, we seek to uncover new mechanisms through which p53 regulates gene expression using tandem affinity purification/mass spectrometry to identify p53-interacting proteins. This approach identified METTL3, an m6A RNA-methyltransferase complex (MTC) constituent, as a p53 interactor. We find that METTL3 promotes p53 protein stabilization and target gene expression in response to DNA damage and oncogenic signals, by both catalytic activity-dependent and independent mechanisms. METTL3 also enhances p53 tumor suppressor activity in in vivo mouse cancer models and human cancer cells. Notably, METTL3 only promotes tumor suppression in the context of intact p53. Analysis of human cancer genome data further supports the notion that the MTC reinforces p53 function in human cancer. Together, these studies reveal a fundamental role for METTL3 in amplifying p53 signaling in response to cellular stress.

Enlightening epigenetics: optochemical tools illuminate the path.

Optochemical tools have become potent instruments for understanding biological processes at the molecular level, and the past decade has witnessed their use in epigenetics and epitranscriptomics (also known as RNA epigenetics) for deciphering gene expression regulation. By using photoresponsive molecules such as photoswitches and photocages, researchers can achieve precise control over when and where specific events occur. Therefore, these are invaluable for studying both histone and nucleotide modifications and exploring disease-related mechanisms. We systematically report and assess current examples in the field, and identify open challenges and future directions. These outstanding proof-of-concept investigations will inspire other chemical biologists to participate in these emerging fields given the potential of photochromic molecules in research and biomedicine.

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.

PCIF1 Catalyzes m6Am mRNA Methylation to Regulate Gene Expression.

mRNA modifications play important roles in regulating gene expression. One of the most abundant mRNA modifications is N6,2-O-dimethyladenosine (m6Am). Here, we demonstrate that m6Am is an evolutionarily conserved mRNA modification mediated by the Phosphorylated CTD Interacting Factor 1 (PCIF1), which catalyzes m6A methylation on 2-O-methylated adenine located at the 5' ends of mRNAs. Furthermore, PCIF1 catalyzes only 5' m6Am methylation of capped mRNAs but not internal m6A methylation in vitro and in vivo. To study the biological role of m6Am, we developed a robust methodology (m6Am-Exo-Seq) to map its transcriptome-wide distribution, which revealed no global crosstalk between m6Am and m6A under assayed conditions, suggesting that m6Am is functionally distinct from m6A. Importantly, we find that m6Am does not alter mRNA transcription or stability but negatively impacts cap-dependent translation of methylated mRNAs. Together, we identify the only human mRNA m6Am methyltransferase and demonstrate a mechanism of gene expression regulation through PCIF1-mediated m6Am mRNA methylation.

Transcriptome-wide Mapping of Internal N7-Methylguanosine Methylome in Mammalian mRNA.

N7-methylguanosine (m7G) is a positively charged, essential modification at the 5' cap of eukaryotic mRNA, regulating mRNA export, translation, and splicing. m7G also occurs internally within tRNA and rRNA, but its existence and distribution within eukaryotic mRNA remain to be investigated. Here, we show the presence of internal m7G sites within mammalian mRNA. We then performed transcriptome-wide profiling of internal m7G methylome using m7G-MeRIP sequencing (MeRIP-seq). To map this modification at base resolution, we developed a chemical-assisted sequencing approach that selectively converts internal m7G sites into abasic sites, inducing misincorporation at these sites during reverse transcription. This base-resolution m7G-seq enabled transcriptome-wide mapping of m7G in human tRNA and mRNA, revealing distribution features of the internal m7G methylome in human cells. We also identified METTL1 as a methyltransferase that installs a subset of m7G within mRNA and showed that internal m7G methylation could affect mRNA translation.

tRNA epitranscriptome determines pathogenicity of the opportunistic pathogen Pseudomonas aeruginosa.

The success of bacterial pathogens depends on the coordinated expression of virulence determinants. Regulatory circuits that drive pathogenesis are complex, multilayered, and incompletely understood. Here, we reveal that alterations in tRNA modifications define pathogenic phenotypes in the opportunistic pathogen Pseudomonas aeruginosa. We demonstrate that the enzymatic activity of GidA leads to the introduction of a carboxymethylaminomethyl modification in selected tRNAs. Modifications at the wobble uridine base (cmnm5U34) of the anticodon drives translation of transcripts containing rare codons. Specifically, in P. aeruginosa the presence of GidA-dependent tRNA modifications modulates expression of genes encoding virulence regulators, leading to a cellular proteomic shift toward pathogenic and well-adapted physiological states. Our approach of profiling the consequences of chemical tRNA modifications is general in concept. It provides a paradigm of how environmentally driven tRNA modifications govern gene expression programs and regulate phenotypic outcomes responsible for bacterial adaption to challenging habitats prevailing in the host niche.

2'-O-methylation (Nm) in RNA: progress, challenges, and future directions.

RNA 2'-O-methylation (Nm) is highly abundant in noncoding RNAs including ribosomal RNA (rRNA), transfer RNA (tRNA), and small nuclear RNA (snRNA), and occurs in the 5' cap of virtually all messenger RNAs (mRNAs) in higher eukaryotes. More recently, Nm has also been reported to occur at internal sites in mRNA. High-throughput methods have been developed for the transcriptome-wide detection of Nm. However, these methods have mostly been applied to abundant RNAs such as rRNA, and the validity of the internal mRNA Nm sites detected with these approaches remains controversial. Nonetheless, Nm in both coding and noncoding RNAs has been demonstrated to impact cellular processes, including translation and splicing. In addition, Nm modifications at the 5' cap and possibly at internal sites in mRNA serve to prevent the binding of nucleic acid sensors, thus preventing the activation of the innate immune response by self-mRNAs. Finally, Nm has been implicated in a variety of diseases including cancer, cardiovascular diseases, and neurologic syndromes. In this review, we discuss current challenges in determining the distribution, regulation, function, and disease relevance of Nm, as well as potential future directions for the field.

Understanding the redundant functions of the m6A-binding YTHDF proteins.

N 6-methyladenosine (m6A) is the most prevalent modified nucleotide in mRNA, and it has important functions in mRNA regulation. However, our understanding of the specific functions of m6A along with its cytosolic readers, the YTHDF proteins, has changed substantially in recent years. The original view was that different m6A sites within an mRNA could have different functions depending on which YTHDF paralog was bound to it, with bound YTHDF1 inducing translation, while bound YTHDF2 induced mRNA degradation. As a result, each YTHDF was proposed to have unique physiologic roles that arise from their unique binding properties and regulatory effects on mRNA. More recent data have called much of this into question, showing that all m6A sites bind all YTHDF proteins with equal ability, with a single primary function of all three YTHDF proteins to mediate mRNA degradation. Here, we describe the diverse technical concerns that led to the original model being questioned and the newer data that overturned this model and led to the new understanding of m6A and YTHDF function. We also discuss how any remaining questions about the functions of the YTHDF proteins can be readily resolved.

Conserved reduction of m6A RNA modifications during aging and neurodegeneration is linked to changes in synaptic transcripts.

N6-methyladenosine (m6A) regulates mRNA metabolism. While it has been implicated in the development of the mammalian brain and in cognition, the role of m6A in synaptic plasticity, especially during cognitive decline, is not fully understood. In this study, we employed methylated RNA immunoprecipitation sequencing to obtain the m6A epitranscriptome of the hippocampal subregions CA1, CA3, and the dentate gyrus and the anterior cingulate cortex (ACC) in young and aged mice. We observed a decrease in m6A levels in aged animals. Comparative analysis of cingulate cortex (CC) brain tissue from cognitively intact human subjects and Alzheimer's disease (AD) patients showed decreased m6A RNA methylation in AD patients. m6A changes common to brains of aged mice and AD patients were found in transcripts linked to synaptic function including calcium/calmodulin-dependent protein kinase 2 (CAMKII) and AMPA-selective glutamate receptor 1 (Glua1). We used proximity ligation assays to show that reduced m6A levels result in decreased synaptic protein synthesis as exemplified by CAMKII and GLUA1. Moreover, reduced m6A levels impaired synaptic function. Our results suggest that m6A RNA methylation controls synaptic protein synthesis and may play a role in cognitive decline associated with aging and AD.

ADAR1 and its implications in cancer development and treatment.

The family of adenosine deaminases acting on RNA (ADARs) regulates global gene expression output by catalyzing adenosine-to-inosine (A-to-I) editing of double-stranded RNA (dsRNA) and through interacting with RNA and other proteins. ADARs play important roles in development and disease, including an increasing connection to cancer progression. ADAR1 has demonstrated a largely pro-oncogenic role in a growing list of cancer types, and its function in cancer has been attributed to diverse mechanisms. Here, we review existing literature on ADAR1 biology and function, its roles in human disease including cancer, and summarize known cancer-associated phenotypes and mechanisms. Lastly, we discuss implications and outstanding questions in the field, including strategies for targeting ADAR1 in cancer.

Intron retention: importance, challenges, and opportunities.

Recent landmark discoveries have underpinned the physiological importance of intron retention (IR) across multiple domains of life and revealed an unexpected breath of functions in a large variety of biological processes. Despite significant progress in the field, some challenges remain. Once solved, opportunities will arise for discovering more functions of IR.

Concepts and methods for transcriptome-wide prediction of chemical messenger RNA modifications with machine learning.

The expanding field of epitranscriptomics might rival the epigenome in the diversity of biological processes impacted. In recent years, the development of new high-throughput experimental and computational techniques has been a key driving force in discovering the properties of RNA modifications. Machine learning applications, such as for classification, clustering or de novo identification, have been critical in these advances. Nonetheless, various challenges remain before the full potential of machine learning for epitranscriptomics can be leveraged. In this review, we provide a comprehensive survey of machine learning methods to detect RNA modifications using diverse input data sources. We describe strategies to train and test machine learning methods and to encode and interpret features that are relevant for epitranscriptomics. Finally, we identify some of the current challenges and open questions about RNA modification analysis, including the ambiguity in predicting RNA modifications in transcript isoforms or in single nucleotides, or the lack of complete ground truth sets to test RNA modifications. We believe this review will inspire and benefit the rapidly developing field of epitranscriptomics in addressing the current limitations through the effective use of machine learning.