REVIEW article

Front. Cell. Neurosci., 19 December 2022

Sec. Cellular Neuropathology

Volume 16 - 2022 | https://doi.org/10.3389/fncel.2022.1058083

Dynamic regulation and key roles of ribonucleic acid methylation

  • 1. Community Health Service Center, Geriatric Hospital Affiliated to Wuhan University of Science and Technology, Wuhan, China

  • 2. Brain Science and Advanced Technology Institute, School of Medicine, Wuhan University of Science and Technology, Wuhan, China

  • 3. School of Public Health, Wuhan University of Science and Technology, Wuhan, China

  • 4. Community Health Service Center, Wuchang Hospital, Wuhan, China

Abstract

Ribonucleic acid (RNA) methylation is the most abundant modification in biological systems, accounting for 60% of all RNA modifications, and affects multiple aspects of RNA (including mRNAs, tRNAs, rRNAs, microRNAs, and long non-coding RNAs). Dysregulation of RNA methylation causes many developmental diseases through various mechanisms mediated by N6-methyladenosine (m6A), 5-methylcytosine (m5C), N1-methyladenosine (m1A), 5-hydroxymethylcytosine (hm5C), and pseudouridine (Ψ). The emerging tools of RNA methylation can be used as diagnostic, preventive, and therapeutic markers. Here, we review the accumulated discoveries to date regarding the biological function and dynamic regulation of RNA methylation/modification, as well as the most popularly used techniques applied for profiling RNA epitranscriptome, to provide new ideas for growth and development.

1 Introduction

Ribonucleic acid (RNA) modification plays an important role in linking deoxyribonucleic acid (DNA) to proteins during the transmission of genetic information. More than 100 post-transcriptional RNA modifications have been identified in important biological processes in viruses, archaea, bacteria, and eukaryotes. Methylation accounts for 60% of the total RNA modifications. RNA methylation was discovered as early as the 1970s (), and it is widely distributed in messenger RNA (mRNA), transfer RNA (tRNA), ribosomal RNA (rRNA), small nuclear RNA (snRNA), small nucleolar RNA (snoRNA), and micro RNA (El Yacoubi et al., 2012). Currently, various distinct modifications in natural RNA have been characterized: N1-methyladenosine (m1A), N6-methyladenosine (m6A), N6, 2′-O-dimethyladenosine (m6Am), 5-methylcytosine (m5C), 5-hydroxymethylcytosine (hm5C), pseudouridine (Ψ), and Adenosine to inosine editing (A-to-I editing), etc. (Roundtree et al., 2017a; Walkley and Li, 2017; Zhao et al., 2020). Among these, m6A, m5C, m1A, and Ψ have been well studied (Figure 1). Compared with DNA methylation, RNA methylation is more complex and diverse, participates in the regulation of many more biological processes, and determines the diversification of RNA post-transcriptional modifications. Mutations in approximately half of the currently known RNA methylation enzymes have been linked to human diseases, including cancer, cardiovascular diseases, congenital genetic disabilities, metabolic diseases, neurological disorders, and mitochondrial-related defects (Jonkhout et al., 2017; ). Here, we review the discoveries to date regarding the dynamic regulation and key roles of RNA methylation.

FIGURE 1

2 Biological patterns of RNA methylation/modification

2.1 N6-methyladenosine (m6A)

N6-methyladenosine (m6A), the most abundant methylation modification in eukaryotic mRNA, is the most thoroughly studied type of RNA modification (Schwartz et al., 2013). m6A has been found in many eukaryotes, ranging from yeast, Arabidopsis, and Drosophila to mammals and even in viruses (Krug et al., 1976; Schwartz et al., 2013; Luo et al., 2014; Lin et al., 2017; Guo et al., 2018). In mammals, m6A is widely distributed in multiple tissues, with the highest level of m6A in the brain, kidney, and liver (Meyer et al., 2012). Although the existence of m6A on RNA was identified as early as the 1970s, little was known about its precise location, temporal dynamics, and regulation until 2011 (Desrosiers et al., 1974). With the identification of the first demethylase fat mass and obesity-associated protein (FTO) of RNA m6A and the development of antibody enrichment and high-throughput sequencing technology, accurate mapping of the transcriptome-level m6A distribution has been achieved (Jia et al., 2011). The sequence near the m6A methylation site on mRNA is highly conserved and mainly occurs on the adenine of RRm6ACH, where R = G/A (G > A) and H = U/A/C (U > A > C) (Harper et al., 1990). In addition, whole-transcriptome m6A sequencing and comprehensive bioinformatics analysis of human and mouse samples revealed that m6A modifications are species-specific (Liu et al., 2020a). m6A methylation is a reversible dynamic modification in many species and is jointly regulated by methyltransferases (writers), demethylases (erasers), and binding proteins (readers). The main components of the methyltransferase complex that have been identified include methyltransferase 3 (METTL3), methyltransferase 14 (METTL14), Wilms’ tumor 1-associating protein (WTAP), and KIAA1429 (vir-like m6A methyltransferase associated) (Oerum et al., 2021). The demethylases FTO and AlkB homolog 5 (ALKBH5) can reverse methylation (Figure 2; Liu et al., 2013; Zheng et al., 2013).

FIGURE 2

The biological function of m6A is mainly through the post-transcriptional regulation of RNA by m6A-binding proteins. Currently known binding proteins include the YT521-B homology (YTH) domain proteins (YTHDF1, YTHDF2, YTHDF3, YTHDC1, and YTHDC2) and nuclear heterogeneous ribonucleoprotein HNRNP family proteins (HNRNPA2B1, HNRNPC, and HNRNPG) (; Roundtree et al., 2017b). With the development of RNA modification detection technology, m6A modifications have been determined to functionally regulate the transcriptome of eukaryotes and processes such as mRNA stability, splicing, nucleation, localization, and translation (Figure 2). Furthermore, m6A is involved in a variety of biological processes such as stem cell differentiation, cell division, gametogenesis, and biological rhythms. Under the catalytic regulation of relevant enzymes, m6A participates in various diseases, including tumors, obesity, and infertility (Jiang et al., 2021).

Knocking out the m6A demethylase gene ALKBH5 accelerates the nuclear export of mRNA (Zheng et al., 2013), whereas RNAi silencing of the core of the m6A methyltransferase complex METTL3 inhibits the nuclear transport process (Fustin et al., 2013), thereby demonstrating that m6A modification can promote nuclear transfer. The nuclear export of mRNA is a bridge connecting mRNA processing in the nucleus and translation in the cytoplasm. Further, nuclear export is coupled with each step of pre-mRNA processing, and only properly processed mRNA is nucleated and translated into protein.

As an exocyclic ring amine involved in Watson Crick base pairing, the rotation direction of Watson Crick base pairing of m6A is opposite to the U force of the carbon-nitrogen bond, showing the methyl group in the reverse conformation, which destabilizes the RNA double strand to local unstructured transcripts (Roost et al., 2015). m6A tends to weaken these structures. This modification is enriched in the 3′-untranslated region (UTR) and alternative splicing exons and introns and is related to mRNA preprocessing, such as splicing regulation and polyadenylation (Meyer et al., 2012). Structural switches, including helices with buried binding sites, can be refolded to allow access to their respective protein partners (Liu et al., 2015). The formation of hybrid structures, such as miRNA-target interactions, is also affected by m6A ().

2.2 5-methylcytosine (m5C)

RNA m5C modifications were discovered more than 40 years ago. Similar to m6A, RNA m5C modifications are dynamically reversible. Methyltransferase uses S-adenosyl methionine (SAM) as a methyl donor to methylate cytosine (C) to form m5C (Figure 3). The distribution and function of m5C may be species- and tissue-specific. Bisulfite treatment combined with transcriptome sequencing has identified m5C modification sites in thousands of mRNAs in HeLa cells (Squires et al., 2012). m5C antibody immunoprecipitation combined with bisulfite sequencing allowed the identification of multiple m5C modifications in archaeal mRNA. The conserved sequence of AU (m5C) GANGU was consistent with the m5C conserved sequence in archaeal rRNA (Edelheit et al., 2013). The discovery of the m5C methyltransferase NSUN2 (NOP2/Sun RNA methyltransferase family member 2) and the binding protein ALYREF (Aly/REF export factor) also proved that RNA m5C modifications have dynamic reversibility (Yang et al., 2017b). RNA m5C methylation is widespread in cells and plays an important role in various physiological (Flores et al., 2017) and pathological processes, such as tumors (Huang et al., 2021; Nombela et al., 2021), neurological disorders (Van Haute et al., 2019), viral infections (Wnuk et al., 2020), and organism ontogeny (Flores et al., 2017). However, research on RNA m5C methylation is still in its infancy.

FIGURE 3

RNA m5C modification can be recognized by the transporter linker protein ALYREF, which promotes the transport of related mRNA out of the nucleus. Knockout of the m5C-modified methyltransferase gene NSUN2 blocks the shuttling of the ALYREF protein between the nucleus and cytoplasm (Yang et al., 2017b). In addition, the patterns of m5C distribution on mRNA in relation to cis-acting regulatory motifs and miRNA/RISC-binding sites suggest that this modification may be involved in the post-transcriptional regulation of mRNA metabolism (Squires et al., 2012). NSUN2-mediated methylation is required to process non-coding vault RNAs (vtRNAs) into small vault RNAs (svRNAs); however, consequences in the downstream coding transcripts have not emerged because of the defect of RNA m5C modification (Hussain et al., 2013). ALYREF can recognize m5C in mRNA via a methyl-specific RNA-binding motif and regulate the export of bound transcripts in an NSUN2-dependent manner (Yang et al., 2017b), whereas hm5C, derived from the Tet-dependent oxidation of m6C, preferentially marks mRNAs within coding regions and favors the translation of Drosophila transcripts (Fu et al., 2014; Delatte et al., 2016).

2.3 Other RNA methylation/modifications

In addition to m6A and m5C, other RNA modifications, including m1A, m6Am, hm5C, and pseudouridine (Ψ), have been a hot field of research in recent years.

2.3.1 N1-methyladenosine (m1A)

RNA m1A methylation was first discovered in non-coding RNAs such as rRNA and tRNA and is widely present in prokaryotic and eukaryotic mRNAs. Unlike m6A methylation, m1A methylation occurs at the N1 position of the adenosine base group. It carries a positive charge by blocking Watson-Crick base-pairing under physiological conditions. Thus, it can drastically alter protein-RNA interactions and RNA secondary structures through electrostatic effects. m1A maps uniquely to positions near the translation start and first splice sites in coding transcripts and correlates with the upregulation of translation (Dominissini et al., 2016). m1A can be removed by ALKBH3 and is responsive to various types of cellular stress (Li et al., 2015a). However, the methyltransferases and binding proteins of m1A remain unknown; hence, its specific function and mechanism of action need to be further explored.

In mRNA, m1A exists in the highly structured 5′UTRs, indicating that it may change the predicted secondary structure (Dominissini et al., 2016; Li et al., 2016c). In the loop structure, this positive charge may stabilize interactions with the RNA phosphate backbone. The methylation level of m1A in transcripts is related to increased translation, which may be due to the availability or direct recruitment of initiation factors and extension factors. The positive charge of this modification enables it to adapt to specific protein RNA and unique RNA-RNA interactions, and its biological effects are still unclear.

2.3.2 N6, 2′-O-dimethyladenosine (m6Am)

The N6, 2′-O-dimethyladenosine (m6Am) modification was found at the first nucleotide of certain mRNAs (Crain et al., 1978). m6Am is formed by the combination of 2′-O-methyltransferase (2′-O-MTase) and 2′-O-methyladenosine-N6-methyltransferase. Adenosine is methylated by 2′-O-methyltransferase (2′-O-MTase) to form Am, which can then be methylated by 2′-O-methyladenosine-N6-methyltransferase to form m6Am (Wei et al., 1975). Compared to m6A, the m6Am level in RNA is very low. Studies in H1-ESCs and GM12878 cells found that H1-ESC poly(A) + RNAs contained approximately three m6A nucleotides per 105 nucleotides compared to ∼100 m6A nucleotides per 105 nucleotides, revealing 33 times more m6A compared to m6Am (Molinie et al., 2016). m6Am and m6A share similar chemical characteristics; m6Am can be detected by m6A-seq, and it can be demethylated by FTO and is involved in the regulation of mRNA stability (Mauer et al., 2016). demonstrated that METTL4 mediates the N6-methylating process of Am30 on U2 small nuclear RNA (snRNA) under an AAG motif in vitro and in vivo (). However, the binding protein of m6Am and molecular mechanisms involved in regulating its biological functions require further study.

2.3.3 5-hydroxymethylcytosine (hm5C)

Similar to m5C in DNA, m5C in RNA can be oxidized by ten-eleven translocation (Tet)-family enzymes to hm5C (Figure 4; Fu et al., 2014). The Tet family of Fe(II)- and 2-oxoglutarate-dependent dioxygenases can induce the oxidation of m5C to yield hm5C (Kriaucionis and Heintz, 2009; Ito et al., 2011). In Drosophila melanogaster, which lacks DNA hydroxymethylation, hm5C is present in greater than 1,500 mRNAs, and hm5C modification is mainly distributed in the exons of mRNA and contains CU-enriched conserved motifs (Delatte et al., 2016). hm5C modification can promote the translation efficiency of mRNA, which is significantly higher than that of RNA without hm5C modifications. After TET knockout, the level of hm5C modification in RNA is reduced, resulting in the abnormal development of the Drosophila brain (Delatte et al., 2016).

FIGURE 4

2.3.4 Pseudouridine (Ψ)

In addition to methylated modifications, there are some unmethylated modifications of RNA, such as ψ (sometimes referred to as pseudouracil). It is one of the most abundant forms of post-transcriptional RNA modifications, widely present in cellular RNA, and highly conserved among species. ψ is formed by the sequence-specific isomerization of uracil (U) (Figure 5), which is abundant in tRNA and rRNA (Goodman et al., 1968). Ψ is known to affect the secondary structure of RNA, and its function in altering stop codon read-through may also be biologically relevant (Karijolich and Yu, 2010; Fernandez et al., 2013).

FIGURE 5

3 Dynamic regulation of RNA methylation/modifications

3.1 Methyltransferase-writers

3.1.1 Core m6A methyltransferase complex components: METTL3, METTL14, and WTAP

The RNA m6A methyltransferase holoenzyme complex consists of at least two multicomponent factors of the whole methyltransferase called MT-A (∼ 200 kDa) and MT-B (∼ 800 kDa). However, only the METTL3 (methyltransferase like 3, ∼ 70 kDa) protein has been identified, and the subunit alone had no enzymatic activity (Tuck, 1992). METTL3 is widely present in various human tissues, especially in the testes (Leach and Tuck, 2001). METTL3 has two key domains that are used to combine SAM and catalyze the formation of m6A (Figure 6). Knocking down METTL3 causes apoptosis in human HeLa and HepG2 cells, accompanied by a significant decrease in m6A levels (Dominissini et al., 2012). Apart from the METTL3 component, other components of the methyltransferase complex have not been comprehensively studied. In 2014, an evolutionary analysis of the METTL3 family revealed the METTL14 and METTL4 genes, which are highly homologous to different families of METTL3 (). METTL14 is another component in the m6A methyltransferase complex with enzymatic activity (Liu et al., 2013). In HeLa and HEK 293FT cells, knocking down METTL14 resulted in a decrease in total mRNA of m6A content, while METTL14 and METTL3 interacted with each other. In vitro size exclusion chromatography (gel filtration) experiments and two-dimensional gel electrophoresis analysis showed that METTL14 and METTL3 can form a stable complex at a 1:1 ratio. In vitro experiments on the activity of enzymes involved in m6A formation showed that although a single METTL14 had slightly higher enzymatic activity than METTL3, the heterodimer formed by METTL14 and METTL3 with a strong preference for substrates had the highest enzymatic activity for the RNA of the stem-loop structure, which has no obvious secondary structures. It preferentially methylates GGACU, which is consistent with the conserved sequence of m6A distribution reported previously (Liu et al., 2013).

FIGURE 6

The third component of the mammalian m6A methyltransferase complex is WTAP (Ping et al., 2014b). WTAP plays a key role in transcriptional and post-transcriptional regulation and is an oncogene whose expression is elevated in several types of human tumors. Through a yeast two-hybrid experiment, FIP37 (FKBP12 interacting protein 37, At3g54170, which encodes a human WTAP homologous gene) was found in Arabidopsis as a protein that interacts with methylthioadenosine (MTA) (Zhong et al., 2008). In HeLa cells, WTAP can co-precipitate with the METTL3-METTL14 heterodimer, and their combined effect is weaker than the interaction between METTL3 and METTL14 (Liu et al., 2013). WTAP lacks a catalytic region and has no methyltransferase activity, but the interaction between WTAP and METTL3-METTL14 can localize the methyltransferase complex in nuclear spots (Ping et al., 2014a), thereby regulating the binding of the methyltransferase complex to the target RNA, which ultimately affects the m6A level. Downregulating WTAP can suppress the localization of METTL3 in the nucleus and reduce m6A levels. The field has demonstrated that the consistency of the RNA substrates of the WTAP, METTL3, and METTL14 combination reached 36% using immunoprecipitation (PAR-CLIP) experiments (Liu et al., 2013). The conserved sequence of RNA binding is the same as the previously reported conserved sequence of m6A, and the three sites of RNA binding are mainly located in the intergenic and intron regions, indicating that m6A is generated from the precursor RNA. Silencing of the methyltransferase complex increases the expression of bound RNA, indicating that m6A is inversely related to gene expression (; Schwartz et al., 2014a). Because WTAP is a splicing factor, knocking down WTAP or METTL3 can result in the formation of different isoforms of RNA containing m6A (Dominissini et al., 2012; Ping et al., 2014a), indicating that m6A affects alternative splicing of RNA.

Taken together, the discovery of METTL3, METTL14, and WTAP, which are important active center components in the m6A methyltransferase complex (Yao et al., 2020), has provided many new insights into m6A, laying an important foundation for revealing the biological function of m6A in RNA.

3.1.2 Other m6A methyltransferase complex components

3.1.2.1 KIAA1429

KIAA1429 is an important methyltransferase that participates in m6A modifications (Schwartz et al., 2014b). Based on the protein immunoprecipitation-mass spectrometry analysis of the core component of m6A methyltransferase, KIAA1429 may be a new subunit component of the m6A methyltransferase complex to catalyze the formation of m6A in mRNA (Schwartz et al., 2014). In Drosophila, the homologous gene of KIAA1429 interacts with the homologous gene of WTAP and regulates the pre-mRNA selective splicing of the important gene sex lethal (Sxl), which controls sex determination (Ortega et al., 2003). In human A549 cells, deletion of KIAA1429 can cause a sharp drop in the m6A peak, indicating that KIAA1429 plays an important role in the methyltransferase complex (Schwartz et al., 2014). Another study found that KIAA1429 can recruit METTL3/METTL14/WTAP, the catalytic core of the m6A methyltransferase complex, to achieve site-specific regulation of m6A levels in mRNA (Yue et al., 2018).

3.1.2.2 RNA binding motif protein 15(B) (RBM15/RBM15B)

RBM15/15B, a protein-coding gene, is a part of the WTAP-METTL3 m6A methyltransferase complex that interacts with WTAP to recruit the complex to target mRNAs. RBM15 and RBM15B can participate in catalyzing the formation of m6A on mRNA and non-coding RNA X-inactive specific transcripts (XIST). Further, immunoprecipitation experiments showed that RBM15 and RBM15B could bind and recruit the WTAP-METTL3 complex to specific target sites. RBM15/RBM15B consists of an RNA recognition domain, as per iCLIP (individual nucleotide resolution crosslinking and immunoprecipitation) sequencing and single miCLIP (m6A individual nucleotide resolution crosslinking and immunoprecipitation). Sequencing analysis of the base resolution revealed that the RBM15/RBM15B-binding site was significantly enriched near the m6A methylation modification site (Patil et al., 2016). In addition, the homologous protein Nito of RBM15 was also identified as a component of the methylase complex in Drosophila (Lence et al., 2016).

3.1.2.3 METTL16

METTL16, a homologous protein of METTL3, can mediate the formation of m6A in U6 snRNA. A 2017 study showed that METTL16 could also mediate the methylation of mRNA and regulate the intracellular levels of SAM. METTL16 contains methylation marks within the conserved UACm6AGAGAA sequence in the 3′UTR hairpins of the mRNA MAT2A and a special type of RNA structure (Figure 7; Pendleton et al., 2017).

FIGURE 7

3.1.2.4 Zinc finger CCCH-type containing 13 (Zc3h13)

Zc3h13, a critical RNA m6A regulator, plays an important role in modulating RNA m6A methylation in the nucleus (Wen et al., 2018). In 2018, three studies almost simultaneously described Zc3h13 as a new component of the m6A methyltransferase complex that regulates m6A modification. Wen et al. (2018) found that mouse Zc3h13 can stabilize the nuclear localization of the Zc3h13-WTAP-Virilizer-Hakai complex and regulate the self-renewal of mouse embryonic stem cells (mESCs) by promoting m6A methylation (Figure 8). Knuckles et al. (2018) further clarified that mouse Zc3h13 and its Drosophila homologous gene CG7358 (named Flacc) mediate the interaction between RBM15/Nito and WTAP/Fl (2)d, thereby promoting m6A modification of Mrna. They also found that Flacc regulates sex determination and dosage compensation in Drosophila by regulating alternative splicing of the Sxl gene. Guo et al. (2018) identified the CG7358 gene of Drosophila (named Xio) and its involvement in the sex determination pathway of Drosophila through regulating the alternative splicing of Sxl via m6A modification.

FIGURE 8

In addition to the above-mentioned members, there are other methyltransferase components involved in the selective recognition of methylation sites and division of labor to ensure fine post-transcriptional regulation.

3.1.3 m5C methyltransferase: NSUN2 and DNMT2

The identified RNA m5C-specific methyltransferases include the NSUN (NOL1/NOP2/SUN) family proteins and DNMT2 (DNA methyltransferase 2), which belong to the superfamily of Rosman-folded enzymes and are characterized by the conserved cysteine residue SAM. The RNA m5C methyltransferase that has been studied the most recently is the NSUN family, which acts as a methyl donor to catalyze the transfer of methyl groups to cytosine residues of different RNA substrates (). Among the nine members of this family, many have catalytic and release sites for methyltransferases, particularly NSUN2 (also known as Trm4). NSUN2, encoded by the NSUN2 gene on chromosome 5p15.31-33, is a nucleolar RNA methyltransferase that catalyzes the m5C methylation of various RNAs such as mRNA, tRNA, and ncRNA (Figure 9). Ultrahigh-performance liquid chromatography-triple quadrupole mass spectrometry with multiple-reaction monitoring (UHPLC-QQQ-MRM-MS/MS) studies have found that NSUN2 is an mRNA-specific m5C methyltransferase, and its catalytic activity depends on the C271 (cysteine 271) and C321 (cysteine 321) sites, where C321 catalyzes the methylation of cytosine by binding to the cytosine pyrimidine ring to form a covalent bond, and C271 mediates the release of RNA. High-throughput sequencing analysis also showed that NSUN2 knockdown significantly reduced mRNA m5C modification (Yang et al., 2017b).

FIGURE 9

In addition to NSUN2, DNMT2 is an RNA m5C methyltransferase that plays a role in multiple species. It was initially considered a DNA methyltransferase, but the latest research shows that human DNMT2 cannot catalyze DNA methylation but can catalyze miRNA methylation and the changes in the C38 methyl group of tRNAAsp (aspartatetRNA) (Figure 9; Goll et al., 2006).

3.2 Demethylase-erasers

3.2.1 m6A demethylase: FTO and ALKBH5

During brain development, a dynamic change in m6A levels in RNA has been observed (Meyer et al., 2012), suggesting the existence of RNA m6A demethylases. Two novel mRNA m6A demethylases (FTO and ALKBH5) have recently been identified, confirming the dynamic regulation of m6A (Jia et al., 2011; Zheng et al., 2013).

The FTO gene is a major regulator of metabolism and energy utilization (, ; Fischer et al., 2009). It is a member of the Fe(II)- and oxoglutarate-dependent AlkB oxygenase family and was originally shown to catalyze the oxidative demethylation of methylated thymidine and uracil (Gerken et al., 2007; Jia et al., 2008). Overexpressing FTO in mice causes excessive energy intake and fattening; knocking out FTO makes mice thinner and grow slowly and causes teratogenic and lethal phenomena. In addition, FTO is closely related to diseases such as type 2 diabetes, cancer, and dementia (; Fu et al., 2013; Zhao et al., 2014). FTO can remove 3mT and m3U from single-stranded DNA and RNA (Jia et al., 2008), but its activity is very weak, and these nucleic acid modifications are rarely found in the body. It is generally believed that these are not the actual substrates for FTO. In 2011, m6A was reported to be an enzyme substrate for FTO (Jia et al., 2011). By simulating the demethylation reaction conditions in the physiological environment in vitro, the wild-type and mutant proteins of the obesity gene FTO were incubated with methylated substrates, and various methylation forms were analyzed using mass spectrometry and high-performance liquid phase technology. FTO had a demethylation function for m6A on single-stranded RNA, and the enzyme kinetic data indicated that the activity of FTO on m6A was much higher than that on m3U. In vivo experiments showed that the m6A content of mRNA in cells with knockdown of the FTO gene increased by approximately 23%; overexpression of FTO reduced the m6A content of mRNA by approximately 18% (Jia et al., 2011). Both in vivo and in vitro data indicate that the actual demethylase substrate of FTO is m6A on nuclear RNA. MeRIP-seq analysis of FTO-knockout brain tissue showed that the m6A level was almost unaffected, but the m6A level of certain specific mRNAs was significantly increased, indicating that FTO acts only on specific mRNA substrates (Hess et al., 2013). Epidemiological fluorescence experiments indicated that FTO is partially located in nuclear spots, suggesting that nuclear RNA is the main substrate of FTO (Jia et al., 2011).

Soon after FTO was discovered, a study reported a second m6A demethylase, ALKBH5, in mammals (Zheng et al., 2013). Integrating mass spectrometry, cell biology, genomics, bioinformatics, and model biology revealed the catalytic m6A demethylation activity of ALKBH5 (RNA demethylase). As the second identified m6A demethylase, ALKBH5 has m6A demethylase activity comparable to that of FTO and preferentially acts on the m6A site of conserved motifs. Unlike the oxidative demethylation of FTO, ALKBH5 can directly catalyze the removal of m6A without generating intermediate products (Zheng et al., 2013). In addition to FTO and ALKBH5, more m6A demethylases are yet to be discovered.

3.3 Binding proteins-readers

3.3.1 m6A readers

The discovery of m6A methyltransferase and demethylase in RNA means that m6A is a reversible chemical modification that can dynamically regulate physiological development processes. It is currently known that m6A-binding proteins in mammals are mainly of families containing YTH (YT521-B homology) domains, including YTHDF1, YTHDF2, YTHDF3, YTHDC1, YTHDC2, and heterogeneous nuclear ribonucleoproteins (HNRNP). Readers mainly work by directly or indirectly binding RNA. Direct RNA binding is mediated by a reader-specific recognition modification domain such as the YTH family proteins, including YTHDF1, YTHDF2, YTHDF3, YTHDC1, and YTHDC2, which can directly bind to methylation modifications of the target genes in the nucleus and cytoplasm (Figure 2). It can affect mRNA modifications, including mRNA splicing, nuclear export, translation, and mRNA degradation (Nachtergaele and He, 2017; Zhao et al., 2017). Indirect RNA binding relies on methylation-induced RNA unwinding, thereby exposing hidden protein-binding motifs and making proteins more likely to bind to RNA. HNRNP, known as heterogeneous nuclear ribonucleoproteins and indirect RNA-binding proteins, is a complex of proteins and RNA found in the nucleus of cells during gene transcription and subsequent post-transcriptional modification of newly synthesized RNA. It is involved in various cellular functions, including mRNA stability (Zhao et al., 2009), mRNA transport (Ma et al., 2018), miRNA maturation (Guil and Cáceres, 2007), and telomere biogenesis (Labranche et al., 1998).

3.3.1.1 YT521-B homology domain family members

YTH is a newly discovered domain that can bind to short, denatured single-stranded RNA (Zhang et al., 2010). The human genome contains at least five YTH domain-containing proteins: YTHDF1–3 and YTHDC1–2. Some YTH proteins contain a low-complexity region (LCR) that can form a phase separation of the protein (Ries et al., 2019), causing the cytoplasmic YTHDF-m6A-RNA polymer to become processing bodies (P-bodies). P-bodies are cytoplasmic ribonucleoprotein (RNP) granules, which are aggregates of multiple functional proteins and RNA in the cytoplasm (Wang et al., 2014). The binding capacity of the YTH domain to RNA containing conserved m6A sequences is approximately 1 μmol/L (Wang et al., 2014; Xu et al., 2014, 2015; Zhu et al., 2014).

In 2014, the m6A-binding protein YTHDF2 (YTH domain family, member 2) was first reported to mediate the degradation of m6A-modified mRNA (Wang et al., 2014). Photoactivatable ribonucleoside–enhanced crosslinking and PAR-CLIP sequencing results indicated that YTHDF2 mainly binds to mRNA and some long non-coding RNA, and the binding site is mainly in the 3′UTR rich in GAC sequences, which is very similar to the distribution characteristics of m6A. In general, YTHDF2 colocalizes with the deadenylation and decapsidase complex in the cell, which can direct the RNA substrate to the cytoplasmic subcellular organelle P-body that mediates RNA degradation (Wang et al., 2014). Subsequent studies further showed that YTHDF2 accelerated the degradation of m6A-modified transcripts by recruiting CCR4-NOT adenylate complexes (Du et al., 2016), primarily comprising subunits with deadenylase activity (CNOT6, CNOT6 L, CNOT7, or CNOT8) and regulatory NOT modules (CNOT1, CNOT2, CNOT3, CNOT9, CNOT10, and CNOT11).

Similar to YTHDF2, the second identified m6A-binding protein is YTHDF1 (YTH domain family, member 1), which is also a cytoplasmic protein that can promote the translation efficiency of the m6A-mRNA it binds. Ribosomal map analysis showed that YTHDF1 promotes the binding of substrate RNA to ribosomes, which is conducive to RNA translation (Wang et al., 2015). In addition, YTHDF1 promotes translation initiation, which, in turn, promotes protein synthesis. Downregulation of YTHDF1 will result in a reduction in substrate RNA translation efficiency, and the magnitude of the reduction is directly related to how much YTHDF1 is downregulated. The discovery of the YTHDF1 function indicates that m6A, as a dynamic RNA marker, can effectively regulate protein production.

YTHDC1, an m6A-binding protein in the nucleus, recognizes and binds m6A-containing RNAs and can mediate m6A-regulated mRNA splicing (Xiao et al., 2016). YTHDC1 is located at the YT body of the nucleus. In vitro electrophoretic mobility shift assay (EMSA) experiments showed significant binding of YTHDC1 to RNA containing m6A modifications, while in vivo photoactivatable ribonucleoside enhanced crosslinking and PAR-CLIP sequencing showed that the binding motif of YTHDC1 is consistent with the conserved motif RRACH of m6A. Most of the binding sites were located near the stop codon, which is also consistent with the distribution characteristics of m6A. YTHDC1 can interact with the splicing factor SR protein to regulate the alternative splicing of m6A-containing exons (Xiao et al., 2016). YTHDC1 also promotes XIST-mediated X-chromosome silencing by identifying the m6A modification site on the non-coding gene XIST (Patil et al., 2016). In HeLa cells, YTHDC1 interacts with SRSF3 and the RNA export factor 1 (NXF1) to promote the export of m6A-modified mRNA (Roundtree et al., 2017b).

YTHDF3 interacts with YTHDF1 and YTHDF2 to enhance the binding ability of YTHDF1 or YTHDF2 to RNA containing an m6A-modified substrate, thereby promoting RNA translation or degradation. YTHDF3 may also interact with other proteins to exert cell-specific regulatory functions (Li et al., 2017; Ni et al., 2019).

YTHDC2 enhances the translation efficiency of target mRNA by binding to the m6A-conserved motif (Wojtas et al., 2017). YTHDC2 also plays an important role in spermatogenesis by interacting with the meiosis-specific protein MEIOC to affect substrate stability (; Soh et al., 2017). A study using infertile male and female mice with YTHDC2 knockout presented a defective phenotype during meiosis phase I (Tanabe et al., 2016). In addition to containing the YTH domain, YTHDC2 also contains an RNA-binding domain, a helicase domain, and two ankyrin repeats, which may play a role in the recruitment of RNA secondary structures, RNA binding proteins, and other interacting proteins (; Tanabe et al., 2016). This also suggests that YTHDC2 may play multiple roles in various biological processes.

3.3.1.2 HNRNP family members

The HNRNP protein family is mainly distributed in the nucleus and is an RNA-binding protein that participates in processes such as precursor RNA splicing, transportation, and translation (Han et al., 2010). HNRNPA2B1, a member of the HNRNP family, can bind m6A-modified RNA transcripts both in vivo and in vitro. However, it remains controversial whether HNRNPA2B1 is an m6A-binding protein. found that HNRNPA2B1 can directly bind to m6A modification sites and regulate alternative splicing and pri miRNA processing, while Wu et al. (2018) found that HNRNPA2B1 does not directly bind m6A but acts as a converter based on structural studies. In addition, two other HNRNP family proteins, HNRNPC and HNRNPG, can regulate the processing of RNA transcripts containing m6A modifications. Unlike HNRNPA2B1, HNRNPC and HNRNPG do not directly bind to the m6A site but mediate the alternative splicing process of transcripts containing m6A by recognizing and binding m6A-dependent structural switches (Liu et al., 2015, 2017).

Additionally, studies based on RNA pull-down experiments have detected other potential m6A binding proteins, including ELAV-like RNA binding protein 1 (ELAVL1, also called HuR), FMRP translational regulator 1 (FMR1), leucine-rich pentatricopeptide repeat containing (LRPPRC), and insulin-like growth factor 2 mRNA-binding proteins (IGF2BPs).

3.3.2 m5C-binding proteins: ALYREF

To identify the binding protein of m5C, two oligonucleotide substrates, with and without m5C, were designed. Through oligonucleotide enrichment combined with protein profiling, the mRNA was identified as a nuclear functional complex component. ALYREF binds to m5C-modified RNA oligonucleotides. EMSA and RIP combined with HPLC and sequencing technology demonstrated that ALYREF’s K171 (lysine 171) mutation can significantly reduce its ability to bind m5C, thereby reducing its ability to bind RNA. Further, fluorescence in situ hybridization technology (FISH) showed that the efficiency of mRNA nucleation decreased with the knockdown of NSUN2 and ALYREF. However, backfilling experiments demonstrated that only the wild-type NSUN2 and ALYREF could reverse the reduction in mRNA nucleation caused by the knockdown of NSUN2 and ALYREF, suggesting that m5C plays an important regulatory role in the process of mRNA nucleation (Shi et al., 2017; Yang et al., 2017a; Fan et al., 2019).

4 Biological processing of RNA methylation

4.1 m6A and neural stem cells

m6A is distributed in distinct developmental stages in the brain and controls the self-renewal of neural stem cells (NSCs) and their differentiation into neurons (Yoon et al., 2017; Du et al., 2019; Zhu et al., 2020). In METTL3/METTL14 gene knockout mice, m6A can promote the proliferation of NSCs and prolong the cell cycle of radial glial cells (Figure 6; Yoon et al., 2017; Wang et al., 2018b). Such regulation may promote the development of stem cell and gene-targeted therapies for Alzheimer’s, Parkinson’s, and cognitive-related neurological diseases.

4.2 m6A and synaptic functions

m6A modifications also contribute to the local regulation of synaptic functions (Flamand and Meyer, 2019). Experimental studies have shown that FTO protein in neurons can shift between the nucleus, cell body, and dendrites (including synapses) and can lead to changes in local RNA methylation kinetics (Meyer et al., 2012; ). Methylated transcripts are highly biased toward neuronal genes and functions such as synaptic functions. The m6A transcriptome is spatially regulated in different brain regions. At the single neuron level, m6A-modified RNA and its interaction groups diffuse to specific structures such as axons, dendrites, pre-synaptic nerve endings, and dendritic spines. This spatial distribution supports the m6A’s functional library to control synapses, touch transmission, and plasticity. Transcriptomic analysis of m6A in a mammalian brain revealed a specific bias of m6A toward neuronal genes rather than glial cell genes. The functional classification of m6A-seq analysis in the entire mouse brain, midbrain, cortex, and cerebellum has been identified. Most m6A-target genes are involved in nervous system development, synaptic transmission, and post-synaptic function. Other studies have reported the localization of m6A in axons and its role in axon growth. Transcripts encoding the axon elongation factor growth-associated protein-43 (GAP-43), a specific marker for axonal regeneration, were identified as targets for m6A modifications. The local translation of these transcripts is negatively regulated by m6A and can be regulated by FTO in axons (Geuens et al., 2016; Yu et al., 2018), showing that RNA methylation also plays an important role in the development and function of synapses.

4.3 m6A regulates nervous system development

Among the various tissues and organs, the brain has the most abundant expression of m6A and its recognition proteins. m6A-related regulatory proteins perform important functions in the cerebral cortex, such as synaptic function, axon regeneration, neural stem cell self-renewal, and cerebellar development (Yoon et al., 2017). The m6A modification plays an irreplaceable regulatory role in the development and function of the nervous system. Drosophila with a knockout of a yeast m6A methyltransferase encoded by the Ime4 gene can be born and survive until adulthood; however, the life span of the fruit fly is shortened and shows obvious abnormal behavior, suggesting that the deletion of m6A modification affects the function of the fruit fly nervous system (Lence et al., 2016). A specific knockout of METTL14 in the central nervous system of mice can affect the development of the mouse cerebral cortex (Yoon et al., 2017). Moreover, the deletion of the YTHDF2 gene in mice leads to an increase in the overall level of m6A, preventing the differentiation of neural stem cells and the formation of neuronal axis dendrites from entering the RNA degradation pathway, resulting in the asymmetric division of neural stem cells in the cerebral cortex (Li et al., 2018b). When neural precursor cells are missing in mice, the differentiation of neurons is affected, resulting in the slower development of the forebrain cortex.

Except for the cerebral cortex, cerebellar RNA m6A patterns and levels are particularly prominent. The dynamic process of m6A methylation and demethylation occurs throughout the entire developmental process of the cerebellum after birth, and the lack of ALKBH5 in a low-pressure and low-oxygen environment causes the m6A level of genes involved in the regulation of cerebellar development to be disordered, speeding up the process of RNA export and leading to a marked lag in cerebellar development (Ma et al., 2018). The specific knockout of METTL3 in the CNS causes severe motor dysfunction in mice during lactation and leads to death. Anatomical and pathological examination revealed that the deletion of m6A methylation caused by METTL3 knockout severely affected the development of the cerebral cortex and cerebellum, resulting in thinning of the cerebral cortex and dysplasia of the cerebellum (Wang et al., 2018a). The absence of m6A also causes disordered gene expression regulation during the differentiation and maturation of granule neurons in the cerebellum, resulting in severe apoptosis of newborn granule neurons, revealing that METTL3-mediated m6A modification plays an important role in the development of the mammalian central nervous system (Wang et al., 2018a).

4.4 m5C is involved in brain development

m5C-modified methyltransferase NOP2/Sun RNA methyltransferase 2 (NSUN2), a protein-coding gene, also regulates neuronal development. In a developing mouse brain, the absence of NSUN2 did not affect RG cells but delayed their differentiation into superior neurons. NSUN2 is expressed in early neuroectodermal cells and can differentiate into various region-specific neuronal and glial cell types (Li et al., 2005). Mutations or deletions in the m5C methylase genes NSUN2 and NSUN3 can cause defects in the nervous system. NSUN2 mutations are associated with autosomal recessive mental retardation (Khan et al., 2012). In the mouse brain, NSUN2 is localized in the nucleolus of Purkinje cells in the cerebellum. The mutated NSUN2 cannot function normally because it cannot converge in the nucleolus. Further research revealed that the methylation level of tRNA in the brains of NSUN2-deficient mice decreased, the cerebral cortex thickened, the number of intermediate progenitors increased, and the number of upper-layer neurons decreased. In this process, the loss of NSUN2 may result in a reduction in the size of the hindbrain due to the failure to produce sufficient numbers of differentiated neurons, while the migration ability of human neuroepithelial stem cells with NSUN2 loss is significantly reduced, which also suggests that NSUN2-dependent tRNA methylation is essential for brain development. The differentiation and migration of neural progenitor cells thus play an important role (Flores et al., 2017). In addition, the inactivation of NSUN3 in mouse embryonic stem cells can further lead to impaired differentiation of the neuroectodermal lineage (Trixl et al., 2018). The summary of RNA modifications is presented in Table 1.

TABLE 1

Modification typesRNAs presentedModifying enzymes and modification readersBiological function
N6-methyladenosine (m6A)mRNAs,
rRNAs,
ncRNAs,
lncRNAs.
Write: METTL3
METTL14
WTAP
KIAA1429
RBM15
RBM15B
METTL16
Reader: YTHDF1
YTHDF2
YTHDF3
YTHDC1
YTHDC2
eIF3
HNRNPA2B1
HNRNPC
Eraser: ALKBH5
FTO
Regulating the transcriptome of eukaryotes and processes such as mRNA stability, splicing, nucleation, localization, and translation.
5-methylcytosine (m5C)mRNAs,
tRNAs,
rRNAs,
miRNAs,
lncRNAs,
circRNAs.
Writer: NSUN1
NSUN2 NSUN3
NSUN4 NSUN5
NSUN6
NSUN7
DNMT2
Reader: ALYREF
YBX1
Eraser: still under debate
1. Stabilizing the secondary structure.
2. Influencing the anticodon stem-loop conformation.
3. Affecting translational fidelity.
N1-methyladenosine (m1A)tRNAs,
rRNAs,
mRNAs.
Writes: TRMT61A TRMT6 TRMT10C
TRMT61B
Reader: remaining controversial
Eraser: ALKBH3 ALKBH1
1. Facilitating translation.
2. Affect tRNA structure folding.
3. Stabilizing the tertiary structures of the tRNA molecules.
Pseudouridine (Ψ)rRNAs,
tRNAs,
snRNAs.
Writes:
PUS1
PUSL1
PUS3
PUS7
PUS7L
PUS10 RPUSD1 RPUSD2 RPUSD3
RPUSD4 TRUB1 TRUB2
DKC1
Reader: none identified
Eraser: none identified
Stabilizing secondary structures to alter translation efficiency, RNA localization, and RNA stability.
N6, 2′-O-dimethyladenosine (m6Am)mRNAs,
rRNAs,
snoRNAs
Writes: METTL4
Reader: none identified
Eraser: FTO
1. Inhibiting adenosine deamination, increase mRNA stability.
2. Controlling RNA stability and translation.

The summary of RNA modifications.

m6A, N6-methyladenosine; m5C, 5-methylcytidine; m1A, N1-methyladenosine; Ψ, pseudouridine; mRNA, messenger RNA; rRNAs, ribosomal RNA; ncRNA, non-coding RNA; lncRNA, long non-coding RNA; tRNAs, transfer RNAs; miRNA, micro RNA; circRNA, circular RNA; snRNAs, small nuclear RNAs; snoRNAs, small nucleolar RNAs.

5 Methods to detect RNA modifications

5.1 Wet-lab approaches applied for profiling RNA epitranscriptome

Targeted approaches for detecting RNA modifications have been well established. These include roughly four types (Vandivier and Gregory, 2017): (1) direct sequencing through nucleotide labeling and chromatography, (2) mass spectrometry, (3) detecting the stalling and termination of reverse transcriptase (RT), and (4) direct measurement of changes in base pairing. We have summarized these targeted technologies for detecting specific modifications (Li et al., 2016b; Helm and Motorin, 2017; Lv et al., 2021). Recently, the progress in combining existing biochemical techniques and high-throughput sequencing has been rapid in covalent RNA modification studies. These technologies include: (1) high-throughput sequencing + antibody pull-down and (2) high-throughput sequencing + chemical conversion and chemical adduct coupling methods (Vandivier and Gregory, 2017). (3) Single-molecule real time (SMRT) technology using nanowells and nanopore sequencing (Furlan et al., 2021). Recently developed in silico methods that detect RT errors in high-throughput RNA sequencing data and high-throughput single molecule sequencing data may read transcriptome and extratranscriptome information simultaneously (Helm and Motorin, 2017; Vandivier and Gregory, 2017). These methods are summarized in Table 2.

TABLE 2

MethodsDetectionProsStrengthsDrawbacks
Direct sequencingVia nucleotide labeling and chromatographyVarious modifications● Powerful and able to produce high quality data● Labor intensive
● Restricting its use to highly abundant RNAs
Direct sequencing with SCARLETSite-specific cleavage and radioactive-labeling followed by ligation-assisted extraction and thin-layer chromatographym6A
Ψ
● Purifing less abundant species of RNA
● Including sequence information
● Quantification possible
● No specialized equipment
● Single-site query
● No high throughput
● Labor intensive
Direct sequencing with mass spectrometryAnalyze fragmented and whole RNAVarious modifications● Unbiased manner
● Highly accurate quantification
● Expertise is reasonably widespread
● Labor-intensive
● No sequence information
● Requires specialized equipment
● Methodological and computational challenges
Reverse transcriptase-based methodsDetecting the stalling and terminationVarious modifications● Targeting transcripts in a heterogeneous pool of RNA
● Ideal approach for studying less abundant RNAs
● Straightforward protocol
● Precise single-nucleotide mapping
● Adaptable to different types of modification
● Semi-quantitative
High-resolution meltingDNA polymorphisms
DNA methylation
covalent RNA modifications
Various modifications● Performing with any existing set of PCR probes
● Covering a putative modification site
● Putative modifications
● Relying on a shift in melting temperature.
Global methods for detecting RNA modifications
Antibody-based enrichment coupled to high-throughput sequencingMethyl and hydroxymethyl RIP-seq, which rely upon antibodies recognizing modified ribonucleotide epitopesm6A
m1A
Ψ
m5C
hm5C
● Unbiased surveys● Modification sites cannot be defined with singlenucleotide resolution
High-throughput sequencing with chemical-based methodsSpecifically target or exclude modified ribonucleotides with high-throughput sequencingm6A
m5C
Ψ
● Determining the location of modification sites
● Single-nucleotide resolution techniques
● Potential false negatives
● Apparent mismatches from the expected sequence
High-throughput single-molecule sequencingDirectly measuring changes in base pairing like SMRT and Nanopore sequencingm6A, m5C, hm5C
Inosine
Ψ
● Much longer sequencing-length
● Allow direct readout of modification sites
● Providing unbiased views of both the transcriptome and epitranscriptome
● Allowing direct quantitation of modification abundance
● Prone to noise and sequencing error
● Statistics problems
● Unmatured base-calling
In silico methodsHigh-throughput analysis of modified ribonucleotidesVarious modifications● Identifying modifications transcriptome-wide with single nucleotide resolution
● Retrospectively and can be readily applied to existing data and in meta analyses.
● Surveying multiple
● modification subtypes simultaneously
● Artifacts not be properly controlled Multiplesteps
● Limited to diploid and haploid organisms

Methods for detecting RNA modifications.

m6A, N6-methyladenosine; m5C, 5-methylcytidine; m1A, N1-methyladenosine; Ψ, pseudouridine; hm5C, 5-hydroxymethylcytosine; SCARLET, site-specific cleavage and radioactive-labeling followed by ligation-assisted extraction and thin-layer chromatography; SMRT, single molecule real-time; snRNAs; small nuclear RNAs; snoRNAs, small nucleolarRNAs.

These high-throughput mapping methods have enabled the characterization of the epitranscriptome in diverse cellular environments. However, further technical improvements are needed to improve the resolution and sensitivity of these methods. Following the recent discovery of more writers, readers, and erasers of the epitranscriptome, as well as the detailed analysis of the known epitranscriptome, new research directions will emerge, which may lead to new treatment strategies. In the future, improved high-throughput technology should help improve spatial resolution, provide chemometric data, and possibly detect new mRNA modifications in the transcriptome.

5.2 Applied computational approaches for profiling RNA epitranscriptome

Although wet-lab approaches can obtain relatively accurate mRNA modification information, they are time-consuming, costly, and difficult to conduct. Given that RNA sequence numbers show explosive growth in the post-genomic era, wet-lab approaches are obviously not suitable for systematic and in-depth analysis of the relevant mechanisms and functions of RNA methylation modification. Therefore, many researchers have developed predictive and computational tools for identifying epigenetic modifications (; Guo et al., 2021), which have evolved rapidly in recent years (Rehman et al., 2021b,2022b; Liao et al., 2022). These tools are mainly based on machine learning (ML) or deep learning (DL) algorithms (). The information for most detection tools and databases has been described in recent reviews published by , El Allali et al. (2021) and Wang et al. (2022b). We have integrated the information from those studies with the most recent advances; this information is listed in Table 3. Furthermore, we categorized emerging tools into two types that identify single modification and multiple RNA methylation sites, respectively, based on the year that the tool developed (Table 3).

TABLE 3

Targets of RNA modificationTools for identifying RNA methylation sites by years
N6-Methyladenosine (m6A)● iRNA-Methyl (), m6Apred () (2015)
● pRNAm-PC (Liu et al., 2016), RNA-MethylPred (Jia et al., 2016), AthMethPre (Xiang et al., 2016b),
RNAMethPre (Xiang et al., 2016a), SRAMP (Zhou et al., 2016), TargetM6A (Li et al., 2016a), M6A-HPCS (Zhang et al., 2016), M6ATH () (2016)
● MethyRNA (), RAM-ESVM (), RAM-NPPS (Xing et al., 2017) (2017)
● iMethyl-STTNC (), M6APred-EL (Wei et al., 2018), RFAthM6A (Wang and Yan, 2018), BERMP (Huang et al., 2018), HMpre (Zhao et al., 2018), Zhang et al.’s (2018a) method, M6AMRFS (Qiang et al., 2018), DeepM6ASeq (Zhang and Hamada, 2018), iRNA (m6A)-PseDNC (), M6Apred-EL (Wei et al., 2018), m6ASNP (Jiang et al., 2018) (2018)
● Gene2vec (Zou et al., 2019), iN6-Methyl (5-step) (Nazari et al., 2019), WHISTLE (), DeepM6APred (Wei et al., 2019), FunDMDeep-m6A (Zhang et al., 2019), iRNA-Freq (Zhuang et al., 2019) (2019)
● LITHOPHONE (Liu et al., 2020e), WITMSG (Liu et al., 2020g), m6A-pred (Khan et al., 2020a), iRNA-m6A (Dao et al., 2020), im6A-TS-CNN (Liu et al., 2020c), iMethyl-deep (Mahmoudi et al., 2020), Pm6A-CNN (), M6A-word2vec (Tahir et al., 2020), m6A Reader (Zhen et al., 2020) (2020)
● TS-m6A-DL (), M6A-GSMS (Zhang et al., 2021c), m6AGE (Wang et al., 2021), m6Aboost (Kortel et al., 2021),
DNN-m6A (Zhang et al., 2021b), EDLm6APred (Zhang et al., 2021a), m6A-NeuralTool (Rehman et al., 2021a) (2021)
● DL-m6A (Rehman et al., 2022a), m6A-TSFinder (Song et al., 2022a) (2022)
5-methylcytosine (m5C)Feng et al.’s (2016) method, m5C-PseDNC (Feng et al., 2016) (2016)
● iRNAm5C-PseDNC (Qiu et al., 2017) (2017)
● pM5CS-Comp-mRMR (Sabooh et al., 2018), M5C-HPCR (Zhang et al., 2018b), PEA-m5C (Song et al., 2018), RNAm5Cfinder (Li et al., 2018a), M5C–HPCR (Zhang et al., 2018b) (2018)
● RNAm5CPred (Fang et al., 2019) (2019)
● iRNA-m5C_SVM (Dou et al., 2020b), m5CPred-SVM (), iRNAm5C_NB (Dou et al., 2020a) (2020)
● Staem5 () (2021)
Li et al.’s (2022b) method (2022)
N1-methyladenosine (m1A)● RAMPred () (2016)
● ISGm1A (Liu et al., 2020f) (2020)
● m1ARegpred (Yao et al., 2022), m1A-Pred (Suleman and Khan, 2022) (2022)
N6, 2′-O-dimethyladenosine (m6Am)● m6AmPred (Jiang et al., 2022), DLm6Am (Luo et al., 2022) (2022)
5-hydroxymethylcytosine (hm5C)● iRNA5hmC (Liu et al., 2020i), iRNA5hmC-PS () (2020)
● iRhm5CNN (), iR5hmcSC (Zhang and Shi, 2021) (2021)
● R5hmCFDV (Shi et al., 2022) (2022)
Pseudouridine (Ψ)● tRNAmod (Panwar and Raghava, 2014) (2014)
● PPUS (Li et al., 2015b) (2015)
● iRNA-PseU () (2016)
● PseUI (He et al., 2018) (2018)
● iPseU-CNN (Tahir et al., 2019), iPseU-NCP (Nguyen-Vo et al., 2019) (2019)
● RF-PseU (Lv et al., 2020b), iPseU–Layer (Mu et al., 2020), PIANO (Song et al., 2020b), EnsemPseU (), PSI-MOUSE (Song et al., 2020a), MixedCNN-PseUI (), MU-PseUDeep (Khan et al., 2020b) (2020)
● PA-PseU (Wang and Zhang, 2021), XG–PseU (Liu et al., 2020d), Aziz et al.’s model (), Porpoise (Li et al., 2021a), PseUdeep (Zhuang et al., 2021) (2021)
Multi-modification type prediction tools.● HAMR (Ryvkin et al., 2013) (2013)
● iRNA-PseColl (Feng et al., 2017) (2017)
● iRNA-3typeA () (2018)
● DeepMRMP (Sun et al., 2019) (2019)
● iMRM (Liu and Chen, 2020), DeepPromise () (2020)
● iRNA-Mod-CNN (Tahir et al., 2021), MultiRM (Song et al., 2021b) (2021)
● EMDLP (Wang et al., 2022a), ZayyuNet () (2022)

Summary of tools for identifying RNA methylation sites.

5.2.1 Prediction tools for identifying m6A sites

Dao et al. (2020) established the iRNA-m6A tool based on the Support Vector Machine (SVM) with fivefold cross-validation test. It can accurate identify m6A sites with the data from high-throughput sequencing techniques in multiple human, mouse, and rat tissues. Liu et al. (2020c) further developed the im6A-TS-CNN tool and improved the results of iRNA-m6A using a convolutional neural network (CNN). Most recently, deep neural network (DNN)-based m6A site computation models, such as TS-m6A-DL, iMethyl-Deep, DNN-M6A, EDLm6APred, and DL-M6A, have been developed. These tools can achieve the identification of m6A methylation sites across species, and are specific to tissue and even RNA types (Mahmoudi et al., 2020; ; Zhang et al., 2021a,b; Rehman et al., 2022a). Kortel et al. (2021) developed m6Aboost by combining ML with miCLIP (m6A individual-nucleotide resolution UV crosslinking and immunoprecipitation) to significantly improve the detection of m6A sites. miCLIP is a kind of antibody-based approach for m6A site mapping with single-nucleotide resolution. Zhang et al. (2021c) developed a novel predictor named M6A-GSMS based on the GBDT (Gradient Boosting Decision Tree) and stacking learning to identify m6A sites (Zhang et al., 2021c). At the same time, Wang et al. (2021) proposed the m6AGE predictor that combines sequence-derived features and graph embeddings for m6A site prediction. Recently, Rehman et al. (2021a) has made use of artificial intelligence to produce an effective model, the m6A-NeuralTool, which can be utilized for speedy and efficient identification of N6-methyladenosine sites (Rehman et al., 2021a). m6A TSHub is a comprehensive online platform established by Song et al. (2022a), which provides a web server constructed by multi-instance deep neural networks with gated attention for high-accuracy prediction of m6A methylation sites, named as m6A-TSFinder (Song et al., 2022a).

5.2.2 Prediction tools for identifying m5C sites

Lv et al. (2020a) established a typical predictor named as iRNA-m5C based on the best features and random forest algorithm, to identify m5C sites in Homo sapiens, Mus musculus, Saccharomyces cerevisiae, and Arabidopsis thaliana (Lv et al., 2020a). Li et al. (2022b) developed a model to identify m5C based on a deep fusion approach with an improved residual network via 10-fold cross-validation and independent set testing, which shows a considerable improvement compared to previous tools (Li et al., 2022b).

5.2.3 Prediction tools for identifying m1A sites

Yao et al. (2022) built the framework m1ARegpred (m1A regulators substrate prediction), m1ARegpred was achieved based on ML and the combination of sequence-derived and genome-derived features. Suleman and Khan (2022) developed an extreme gradient boost predictor named as m1A-Pred for the prediction of modified m1A sites.

5.2.4 Prediction tools for identifying m6Am sites

Jiang et al. (2022) presented the m6AmPred, the first web server, for in silico identification of m6Am sites from the primary sequences of RNA. m6AmPred was built upon the XgbDart (eXtreme Gradient Boosting with Dart algorithm) and EIIP-PseEIIP encoding scheme. Luo et al. (2022) proposed an ensemble DL framework, named as DLm6Am, to identify m6Am sites. DLm6Am consists of three similar base classifiers, each of which contains a multi-head attention module, an embedding module with two parallel DL sub-modules, a CNN and a BiLSTM (Bidirectional long short-term memory), and a prediction module. Compared with the existing state-of-the-art m6Am prediction tool, m6AmPred and MultiRM show superior performance (Luo et al., 2022).

5.2.5 Prediction tools for identifying hm5C sites

Zhang and Shi (2021) designed a novel and powerful model called iR5hmcSC for identifying hm5C. iR5hmcSC can achieve high-throughput identification of hm5C (Zhang and Shi, 2021). Shi et al. (2022) designed a model called R5hmCFDV. R5hmCFDV showed higher accuracy than iR5hmcSC does in the 10-fold cross-validation.

5.2.6 Prediction tools for identifying Ψ sites

Zhuang et al. (2021) built PseUdeep, an RNA Pseudouridine Site Identification framework with DL Algorithm. PseUdeep outperformed the best traditional ML model available, which was evaluated through 10-fold cross-validation and two independent testing data sets (Zhuang et al., 2021).

5.2.7 Prediction tools for multi-modification sites

Wang et al. (2022a) combined convolutional CNN and BiLSTM, and developed an ensemble multiscale DL predictor, EMDLP. It can identify RNA methylation sites by NLP (natural language processing) and DL way, and also take better advantage of the local and global information for site prediction. proposed a unified DL model named as ZayyuNet, which can efficiently receive large inputs and achieve better performance based on its SpinalNet architecture that was inspired by the human somatosensory system.

5.3 Benchmark datasets related to RNA methylation

It is important to understand the methylation-related databases involved in developing these tools. We summarized the database information involved in methylation modification (Table 4) and sorted it into either single site or multiple site modification based on the year that the database was established in. Databases related to methylation modification are under development from comprehensive databases to more detailed and specialized ones. The following are the emerged databases of methylation modifications associated with some specific diseases and fields: (1) After m6A methylation-related genes based on The Cancer Genome Atlas (TCGA) were used to predict the prognosis of hepatocellular carcinoma (Group et al., 2020; Liu et al., 2020b; Li et al., 2021b), cancer-related methylation modification databases have begun to emerge, including, Lnc2Cancer 3.0 and OncoDB (Gao et al., 2021; Tang et al., 2022); (2) Osteoarthritis-omics and molecular biomarkers (OAOB), which are a group of database containing differential molecular biomarkers related to osteoarthritis (Li et al., 2022a); (3) Other than disease, Pan et al. (2022) established an integrative multi-omic database (iMOMdb) of Asian pregnant women providing the first blood-based multi-omic analysis of pregnant women in Asia. This database contains high-resolution genotypes, DNA methylation, and transcriptome profiles, and fills the knowledge gap of complex traits in populations of Asian ancestry; (4) Gao et al. (2022) developed AgingBank, an experimentally supported multiomics database of information related to aging in multiple species; (5) ProMetheusDB, a database generated by analyzing and sorting cell culture experiments data using ML tools from the protein perspective (Massignani et al., 2022); (6) compendium of protein lysine modifications (CPLM 4.0), a post-translational modification (PTMs) database (Zhang et al., 2022); (7) tRNA-related databases containing high-throughput tRNA sequencing data (Sajek et al., 2020); (8) RNAWRE and RM2 Target are two databases focusing on information of writers, readers, and erasers (Nie et al., 2020; ); and (9) SyStemCell, a multiple-levels experimental database for stem cell research (Yu et al., 2012). With the emergence of these specialized and multi-angle RNA methylation related databases, the traditional databases are also constantly updated and developed. Song et al. (2022a) established a comprehensive online platform, m6A-TSHub, to reveal context-specific m6A methylation and gene mutations that may regulate m6A epigenetic markers. Ma et al. (2022) proposed the M5C Atlas, a database for the comprehensive collection and annotation of RNA5 methylcytosine. Zhou et al. (2022) built ASMdb, a DNA modification database containing RNA sequencing data. Nevertheless, at present, specific databases for methylation modification such as m1A, m6Am, hm5C, and Ψ are still relatively lacking, which requires further improvement and development by researchers.

TABLE 4

Targets of RNA methylation/modificationDatabase by years
N6-Methyladenosine (m6A)● MeT-DB V2.0 (Liu et al., 2018), m6Avar (Zheng et al., 2018) (2018)
● CVm6A (Han et al., 2019) (2019)
● REPIC (Liu et al., 2020h) (2020)
● ConsRM (Song et al., 2021a), M6A2Target (Deng et al., 2021) (2021)
● m6A-TSHub (Song et al., 2022a) (2022)
5-methylcytosine (m5C)● SyStemCell (Yu et al., 2012) (2012)
● m5C-Atlas (Ma et al., 2022) (2022)
N1-methyladenosine (m1A)/
N6, 2′-O-dimethyladenosine (m6Am)/
5-hydroxymethylcytosine (hm5C)/
Pseudouridine (Ψ)/
Database containing multiple types of RNA methylation/modification● MODOMICS (Dunin-Horkawicz et al., 2006)
● TCGA (Wang et al., 2016) (2006)
● REACTOME (Croft et al., 2011), RNAMDB () (2011)
● Gene-Expression Omnibus (GEO) (), DARNED (Kiran et al., 2013) (2013)
● RCAS (Uyar et al., 2017) (2017)
● RMBase V2 (Xuan et al., 2018), REDIdb 3.0 (Lo Giudice et al., 2018) (2018)
● RNAmod (Liu and Gregory, 2019) (2019)
● RNAWRE (Nie et al., 2020), T-psi-C (Sajek et al., 2020) (2020)
● m6A-Atlas (Tang et al., 2021), RMVar (Luo et al., 2021), Lnc2Cancer 3.0 (Gao et al., 2021) (2021)
● RMDisease V2.0 (Song et al., 2022b), AgingBank (Gao et al., 2022), CPLM 4.0 (Zhang et al., 2022), OncoDB (Tang et al., 2022), ASMdb (Zhou et al., 2022), iMOMdb (Pan et al., 2022), OAOB (Li et al., 2022a), ProMetheusDB (Massignani et al., 2022), RM2Target () (2022)

RNA modification related database.

6 Conclusion and outlook

To provide biochemical researchers with the latest progress in RNA methylation, we reviewed the findings to date on the dynamic regulation and key roles of RNA methylation. We found that combining biochemical technology with high-throughput sequencing has made rapid progress in understanding the form and function of RNA modification, especially in mRNA and lncRNA. We mainly discussed m6A, m5C, m1A, pseudouridine, and 2’OMe. Other types of modifications that can be detected but are not clearly defined were not included. We also discussed the advantages and disadvantages of the detecting method based on high-throughput sequencing. It is worth mentioning that the databases related to RNA methylation modification, as well as the prediction and identification tools for RNA methylation sites developed with the help of these databases, have played an increasingly important role in RNA methylation modification research. We found that establishing a powerful method to investigate transcriptome modification is key to understanding post-transcriptional regulation. We emphasized that RNA methylation modification plays a key role in neural stem cells, synaptic functions, nervous system development, and brain development. Future research should focus on the role and mechanism of RNA methylation in neurodevelopmental disorders, which will greatly contribute to the prevention and treatment of developmental diseases.

Statements

Author contributions

JZ and HL were responsible for gathering information and drawing the illustrations. YZ contributed to information interpretation, editing, and critical revision of the manuscript. WT, Y-QC, JD, S-YB, and Z-XW were responsible for revision of the manuscript. All authors read and approved the final manuscript.

Funding

This study was financially supported by the Ministry of Science and Technology of China (2020YFC2006000) and the National Natural Science Foundation of China (81870901 and 82071272). The funders had no role in the design, analysis, or writing of this manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1

    AbbasZ.TayaraH.ChongK. T. (2022). ZayyuNet - A unified deep learning model for the identification of epigenetic modifications using raw genomic sequences.IEEE ACM Trans. Comput. Biol. Bioinform.1925332544. 10.1109/Tcbb.2021.3083789

  • 2

    AbbasZ.TayaraH.ZouQ.ChongK. T. (2021). Ts-m6A-Dl: Tissue-specific identification of N6-methyladenosine sites using a universal deep learning model.Comput. Struct. Biotechnol. J.1946194625. 10.1016/j.csbj.2021.08.014

  • 3

    AbbyE.TourpinS.RibeiroJ.DanielK.MessiaenS.MoisonD.et al (2016). Implementation of meiosis prophase I programme requires a conserved retinoid-independent stabilizer of meiotic transcripts.Nat. Commun.7:10324. 10.1038/ncomms10324

  • 4

    AhmedS.HossainZ.UddinM.TaherzadehG.SharmaA.ShatabdaS.et al (2020). Accurate prediction of Rna 5-hydroxymethylcytosine modification by utilizing novel position-specific gapped k-mer descriptors.Comput. Struct. Biotechnol. J.1835283538. 10.1016/j.csbj.2020.10.032

  • 5

    AkbarS.HayatM. (2018). iMethyl-Sttnc: Identification of N(6)-methyladenosine sites by extending the idea of Saac into Chou’s Pseaac to formulate Rna sequences.J. Theory Biol.455205211. 10.1016/j.jtbi.2018.07.018

  • 6

    AlamW.AliS. D.TayaraH.ChongK. T. (2020). A Cnn-Based Rna N6-methyladenosine site predictor for multiple species using heterogeneous features representation.IEEE Access.8138203138209. 10.1109/access.2020.3002995

  • 7

    AlarcónC. R.GoodarziH.LeeH.LiuX.TavazoieS.TavazoieS. F. (2015). Hnrnpa2B1 Is a Mediator of m(6)A-Dependent nuclear rna processing events.Cell16212991308. 10.1016/j.cell.2015.08.011

  • 8

    AliS. D.KimJ. H.TayaraH.ChongK. T. (2021). Prediction of Rna 5-Hydroxymethylcytosine modifications using deep learning.IEEE Access.984918496. 10.1109/access.2021.3049146

  • 9

    AzizA. Z. B.HasanM. A. M. (2020). “A mixed convolution neural network for identifying rna pseudouridine sites,” in Ieee Region 10 Symposium (Tensymp), (Dhaka: Ieee), 799802.

  • 10

    AzizA. Z. B.HasanM. A. M.ShinJ. (2021). Identification of Rna pseudouridine sites using deep learning approaches.PLoS One16:e0247511. 10.1371/journal.pone.0247511

  • 11

    BaoX.ZhangY.LiH.TengY.MaL.ChenZ.et al (2022). Rm2Target: A comprehensive database for targets of writers, erasers and readers of Rna modifications.Nucleic Acids Res.10.1093/nar/gkac945

  • 12

    BarbieriI.KouzaridesT. (2020). Role of Rna modifications in cancer.Nat. Rev. Cancer20303322. 10.1038/s41568-020-0253-2

  • 13

    BarrettT.WilhiteS. E.LedouxP.EvangelistaC.KimI. F.TomashevskyM.et al (2013). Ncbi Geo: Archive for functional genomics data sets–update.Nucleic Acids Res.41D991D995. 10.1093/nar/gks1193

  • 14

    BatistaP. J.MolinieB.WangJ.QuK.ZhangJ.LiL.et al (2014). m(6)A Rna modification controls cell fate transition in mammalian embryonic stem cells.Cell Stem Cell15707719. 10.1016/j.stem.2014.09.019

  • 15

    BiY.JinD.JiaC. (2020). EnsemPseU: Identifying pseudouridine sites with an ensemble approach.IEEE Access.87937679382. 10.1109/access.2020.2989469

  • 16

    BoisselS.ReishO.ProulxK.Kawagoe-TakakiH.SedgwickB.YeoG. S.et al (2009). Loss-of-function mutation in the dioxygenase-encoding Fto gene causes severe growth retardation and multiple malformations.Am. J. Hum. Genet.85106111. 10.1016/j.ajhg.2009.06.002

  • 17

    BrummerA.YangY.ChanT. W.XiaoX. (2017). Structure-mediated modulation of mrna abundance by A-to-I editing.Nat. Commun.8:1255. 10.1038/s41467-017-01459-7

  • 18

    BujnickiJ. M.FederM.AyresC. L.RedmanK. L. (2004). Sequence-structure-function studies of trna:m5C methyltransferase Trm4p and its relationship to Dna:m5C and Rna:m5U methyltransferases.Nucleic Acids Res.3224532463. 10.1093/nar/gkh564

  • 19

    BujnickiJ.FederM.RadlinskaM.BlumenthalR. (2002). Structure prediction and phylogenetic analysis of a functionally diverse family of proteins homologous to the Mt-A70 subunit of the human mrna: m6A methyltran sferase.J. Mol. Evol.55431444. 10.1007/s00239-002-2339-8

  • 20

    CantaraW. A.CrainP. F.RozenskiJ.MccloskeyJ. A.HarrisK. A.ZhangX.et al (2011). The Rna Modification Database. Rnamdb: 2011 update.Nucleic Acids Res.39D195D201. 10.1093/nar/gkq1028

  • 21

    ChaiD.JiaC.ZhengJ.ZouQ.LiF. (2021). Staem5: A novel computational approachfor accurate prediction of m5C site.Mol. Ther. Nucleic Acids2610271034. 10.1016/j.omtn.2021.10.012

  • 22

    ChangM.LvH.ZhangW.MaC.HeX.ZhaoS.et al (2017). Region-specific Rna m(6)A methylation represents a new layer of control in the gene regulatory network in the mouse brain.Open Biol.7:170166. 10.1098/rsob.170166

  • 23

    ChenH.GuL.OrellanaE. A.WangY.GuoJ.LiuQ.et al (2020a). Mettl4 is an snrna m(6)Am methyltransferase that regulates Rna splicing.Cell Res.30544547. 10.1038/s41422-019-0270-4

  • 24

    ChenK.WeiZ.ZhangQ.WuX.RongR.LuZ.et al (2019a). Whistle: A high-accuracy map of the human N6-methyladenosine (m6A) epitranscriptome predicted using a machine learning approach.Nucleic Acids Res.47:e41. 10.1093/nar/gkz074

  • 25

    ChenW.DingH.ZhouX.LinH.ChouK. C. (2018a). irna(m6A)-Psednc: Identifying N(6)-methyladenosine sites using pseudo dinucleotide composition.Anal. Biochem.561-5625965. 10.1016/j.ab.2018.09.002

  • 26

    ChenW.FengP.DingH.LinH. (2016a). Identifying N (6)-methyladenosine sites in the Arabidopsis thaliana transcriptome.Mol. Genet. Genomics29122252229. 10.1007/s00438-016-1243-7

  • 27

    ChenW.FengP.DingH.LinH.ChouK. C. (2015a). irna-Methyl: Identifying N(6)-methyladenosine sites using pseudo nucleotide composition.Anal. Biochem.4902633. 10.1016/j.ab.2015.08.021

  • 28

    ChenW.FengP.TangH.DingH.LinH. (2016b). Rampred: Identifying the N(1)-methyladenosine sites in eukaryotic transcriptomes.Sci. Rep.6:31080. 10.1038/srep31080

  • 29

    ChenW.FengP.YangH.DingH.LinH.ChouK. C. (2018b). irna-3typeA: Identifying Three types of modification at Rna’s adenosine sites.Mol. Ther. Nucleic Acids11468474. 10.1016/j.omtn.2018.03.012

  • 30

    ChenW.TangH.LinH. (2017a). Methyrna: A web server for identification of N(6)-methyladenosine sites.J. Biomol. Struct. Dyn.35683687. 10.1080/07391102.2016.1157761

  • 31

    ChenW.TangH.YeJ.LinH.ChouK. C. (2016c). irna-PseU: Identifying Rna pseudouridine sites.Mol. Ther. Nucleic Acids5:e332. 10.1038/mtna.2016.37

  • 32

    ChenW.TranH.LiangZ.LinH.ZhangL. (2015b). Identification and analysis of the N(6)-methyladenosine in the Saccharomyces cerevisiae transcriptome.Sci. Rep.5:13859. 10.1038/srep13859

  • 33

    ChenW.XingP.ZouQ. (2017b). Detecting N6-methyladenosine sites from Rna transcriptomes using ensemble support vector machines.Sci. Rep.7:40242. 10.1038/srep40242

  • 34

    ChenX.SunY. Z.LiuH.ZhangL.LiJ. Q.MengJ. (2019b). Rna methylation and diseases: Experimental results, databases, Web servers and computational models.Brief Bioinform.20896917. 10.1093/bib/bbx142

  • 35

    ChenX.XiongY.LiuY.ChenY.BiS.ZhuX. (2020b). m5cpred-Svm: A novel method for predicting m5C sites of Rna.BMC Bioinform.21:489. 10.1186/s12859-020-03828-4

  • 36

    ChenZ.ZhaoP.LiF.WangY.SmithA. I.WebbG. I.et al (2020c). Comprehensive review and assessment of computational methods for predicting Rna post-transcriptional modification sites from Rna sequences.Brief Bioinform.2116761696. 10.1093/bib/bbz112

  • 37

    ChurchC.LeeS.BaggE. A.MctaggartJ. S.DeaconR.GerkenT.et al (2009). A mouse model for the metabolic effects of the human fat mass and obesity associated Fto gene.PLoS Genet.5:e1000599. 10.1371/journal.pgen.1000599

  • 38

    ChurchC.MoirL.McmurrayF.GirardC.BanksG. T.TeboulL.et al (2010). Overexpression of Fto leads to increased food intake and results in obesity.Nat. Genet.4210861092. 10.1038/ng.713

  • 39

    CrainP.ChoiY.BuschH.McCloskeyJ. (1978). Characterization of N6, O2′-dimethyladenosine from nuclear Rna of Novikoff hepatoma.Nucleic Acids Res.5771776. 10.1093/nar/5.3.771

  • 40

    CroftD.O’kellyG.WuG.HawR.GillespieM.MatthewsL.et al (2011). Reactome: A database of reactions, pathways and biological processes.Nucleic Acids Res.39D691D697. 10.1093/nar/gkq1018

  • 41

    DaoF. Y.LvH.YangY. H.ZulfiqarH.GaoH.LinH. (2020). Computational identification of N6-methyladenosine sites in multiple tissues of mammals.Comput. Struct. Biotechnol. J.1810841091. 10.1016/j.csbj.2020.04.015

  • 42

    DelatteB.WangF.NgocL.CollignonE.BonvinE.DeplusR.et al (2016). Transcriptome-wide distribution and function of Rna hydroxymethylcytosine.Science351282285. 10.1126/science.aac5253

  • 43

    DengS.ZhangH.ZhuK.LiX.YeY.LiR.et al (2021). M6A2Target: A comprehensive database for targets of m6A writers, erasers and readers.Brief Bioinform.22:bbaa055. 10.1093/bib/bbaa055

  • 44

    DesrosiersR.FridericiK.RottmanF. (1974). Identification of methylated nucleosides in messengerrna from Novikoff Hepatoma cells.Proc. Natl. Acad. Sci. U.S.A.7139713975. 10.1073/pnas.71.10.3971

  • 45

    DominissiniD.Moshitch-MoshkovitzS.SchwartzS.Salmon-DivonM.UngarL.OsenbergS.et al (2012). Topology of the human and mouse m6A Rna methylomes revealed by m6A-seq.Nature485201206. 10.1038/nature11112

  • 46

    DominissiniD.NachtergaeleS.Moshitch-MoshkovitzS.PeerE.KolN.Ben-HaimM. S.et al (2016). The dynamic N(1)-methyladenosine methylome in eukaryotic messenger Rna.Nature530441446. 10.1038/nature16998

  • 47

    DouL.LiX.DingH.XuL.XiangH. (2020a). irna-m5C_Nb: A novel predictor to identify rna 5-methylcytosine sites based on the naive bayes classifier.IEEE Access.88490684917. 10.1109/access.2020.2991477

  • 48

    DouL.LiX.DingH.XuL.XiangH. (2020b). Prediction of m5C modifications in rna sequences by combining multiple sequence features.Mol. Ther. Nucleic Acids21332342. 10.1016/j.omtn.2020.06.004

  • 49

    DuH.ZhaoY.HeJ.ZhangY.XiH.LiuM.et al (2016). Ythdf2 destabilizes m(6)A-containing Rna through direct recruitment of the Ccr4-Not deadenylase complex.Nat. Commun.7:12626. 10.1038/ncomms12626

  • 50

    DuK. Z.ZhangL. B.LeeT.SunT. (2019). m(6)A Rna methylation controls neural development and is involved in human diseases.Mol. Neurobiol.5615961606. 10.1007/s12035-018-1138-1

  • 51

    Dunin-HorkawiczS.CzerwoniecA.GajdaM. J.FederM.GrosjeanH.BujnickiJ. M. (2006). Modomics: A database of Rna modification pathways.Nucleic Acids Res.34D145D149. 10.1093/nar/gkj084

  • 52

    EdelheitS.SchwartzS.MumbachM. R.WurtzelO.SorekR. (2013). Transcriptome-wide mapping of 5-methylcytidine Rna modifications in bacteria, archaea, and yeast reveals m5C within archaeal mrnas.PLoS Genet.9:e1003602. 10.1371/journal.pgen.1003602

  • 53

    El AllaliA.ElhamraouiZ.DaoudR. (2021). Machine learning applications in Rna modification sites prediction.Comput. Struct. Biotechnol. J.1955105524. 10.1016/j.csbj.2021.09.025

  • 54

    El YacoubiB.BaillyM.De Crecy-LagardV. (2012). Biosynthesis and function of posttranscriptional modifications of transfer Rnas.Annu. Rev. Genet.466995. 10.1146/annurev-genet-110711-155641

  • 55

    FanJ.WangK.DuX.WangJ.ChenS.WangY.et al (2019). Alyref links 3′-end processing to nuclear export of non-polyadenylated mrnas.EMBO J.38:e99910. 10.15252/embj.201899910

  • 56

    FangT.ZhangZ.SunR.ZhuL.HeJ.HuangB.et al (2019). Rnam5cpred: Prediction of Rna 5-methylcytosine sites based on three different kinds of nucleotide composition.Mol. Ther. Nucleic Acids18739747. 10.1016/j.omtn.2019.10.008

  • 57

    FengP.DingH.ChenW.LinH. (2016). Identifying Rna 5-methylcytosine sites via pseudo nucleotide compositions.Mol. Biosyst.1233073311. 10.1039/c6mb00471g

  • 58

    FengP.DingH.YangH.ChenW.LinH.ChouK. C. (2017). irna-PseColl: Identifying the occurrence sites of different rna modifications by incorporating collective effects of nucleotides into pseknc.Mol. Ther. Nucleic Acids7155163. 10.1016/j.omtn.2017.03.006

  • 59

    FernandezI. S.NgC. L.KelleyA. C.WuG.YuY. T.RamakrishnanV. (2013). Unusual base pairing during the decoding of a stop codon by the ribosome.Nature500107110. 10.1038/nature12302

  • 60

    FischerJ.KochL.EmmerlingC.VierkottenJ.PetersT.BruningJ. C.et al (2009). Inactivation of the Fto gene protects from obesity.Nature458894898. 10.1038/nature07848

  • 61

    FlamandM. N.MeyerK. D. (2019). The epitranscriptome and synaptic plasticity.Curr. Opin. Neurobiol.594148. 10.1016/j.conb.2019.04.007

  • 62

    FloresJ. V.Cordero-EspinozaL.Oeztuerk-WinderF.Andersson-RolfA.SelmiT.BlancoS.et al (2017). Cytosine-5 Rna methylation regulates neural stem cell differentiation and motility.Stem Cell Rep.8112124. 10.1016/j.stemcr.2016.11.014

  • 63

    FuL.GuerreroC. R.ZhongN.AmatoN. J.LiuY.LiuS.et al (2014). Tet-mediated formation of 5-hydroxymethylcytosine in Rna.J. Am. Chem. Soc.1361158211585. 10.1021/ja505305z

  • 64

    FuY.JiaG.PangX.WangR. N.WangX.LiC. J.et al (2013). Fto-mediated formation of N6-hydroxymethyladenosine and N6-formyladenosine in mammalian Rna.Nat. Commun.4:1798. 10.1038/ncomms2822

  • 65

    FurlanM.Delgado-TejedorA.MulroneyL.PelizzolaM.NovoaE. M.LeonardiT. (2021). Computational methods for Rna modification detection from nanopore direct Rna sequencing data.RNA Biol.183140. 10.1080/15476286.2021.1978215

  • 66

    FustinJ. M.DoiM.YamaguchiY.HidaH.NishimuraS.YoshidaM.et al (2013). Rna-methylation-dependent Rna processing controls the speed of the circadian clock.Cell155793806. 10.1016/j.cell.2013.10.026

  • 67

    GaoY.ShangS.GuoS.LiX.ZhouH.LiuH.et al (2021). Lnc2Cancer 3.0: An updated resource for experimentally supported lncrna/circrna cancer associations and web tools based on Rna-seq and scrna-seq data.Nucleic Acids Res.49D1251D1258. 10.1093/nar/gkaa1006

  • 68

    GaoY.ShangS.GuoS.WangX.ZhouH.SunY.et al (2022). AgingBank: A manually curated knowledgebase and high-throughput analysis platform that provides experimentally supported multi-omics data relevant to aging in multiple species.Brief Bioinform.23:bbac438. 10.1093/bib/bbac438

  • 69

    GerkenT.GirardC. A.TungY. C.WebbyC. J.SaudekV.HewitsonK. S.et al (2007). The obesity-associated Fto gene encodes a 2-oxoglutarate-dependent nucleic acid demethylase.Science31814691472. 10.1126/science.1151710

  • 70

    GeuensT.BouhyD.TimmermanV. (2016). The hnrnp family: Insights into their role in health and disease.Hum. Genet.135851867. 10.1007/s00439-016-1683-5

  • 71

    GollM. G.KirpekarF.MaggertK. A.YoderJ. A.HsiehC. L.ZhangX.et al (2006). Methylation of trnaasp by the Dna methyltransferase homolog Dnmt2.Science311395398. 10.1126/science.1120976

  • 72

    GoodmanH. M.AbelsonJ.LandyA.BrennerS.SmithJ. D. (1968). Amber suppression: A nucleotide change in the anticodon of a tyrosine transfer Rna.Nature21710191024. 10.1038/2171019a0

  • 73

    GroupP. T. C.CalabreseC.DavidsonN. R.DemirciogluD.FonsecaN. A.HeY.et al (2020). Genomic basis for Rna alterations in cancer.Nature578129136. 10.1038/s41586-020-1970-0

  • 74

    GuilS.CáceresJ. F. (2007). The multifunctional Rna-binding protein hnrnp A1 is required for processing of miR-18a.Nat. Struct. Mol. Biol.14591596. 10.1038/nsmb1250

  • 75

    GuoG.PanK.FangS.YeL.TongX.WangZ.et al (2021). Advances in mrna 5-methylcytosine modifications: Detection, effectors, biological functions, and clinical relevance.Mol. Ther. Nucleic Acids26575593. 10.1016/j.omtn.2021.08.020

  • 76

    GuoJ.TangH. W.LiJ.PerrimonN.YanD. (2018). Xio is a component of the Drosophila sex determination pathway and Rna N(6)-methyladenosine methyltransferase complex.Proc. Natl. Acad. Sci. U.S.A.11536743679. 10.1073/pnas.1720945115

  • 77

    HanS. P.TangY. H.SmithR. (2010). Functional diversity of the hnrnps: Past, present and perspectives.Biochem. J.430379392. 10.1042/Bj20100396

  • 78

    HanY.FengJ.XiaL.DongX.ZhangX.ZhangS.et al (2019). Cvm6A: A visualization and exploration database for m(6)As in cell lines.Cells8:168. 10.3390/cells8020168

  • 79

    HarperJ. E.MiceliS. M.RobertsR. J.ManleyJ. L. (1990). Sequence specificity of the human mrna N6-adenosine methylase in vitro.Nucleic Acids Res.1857355741. 10.1093/nar/18.19.5735

  • 80

    HeJ.FangT.ZhangZ.HuangB.ZhuX.XiongY. (2018). Pseui: Pseudouridine sites identification based on Rna sequence information.BMC Bioinform.19:306. 10.1186/s12859-018-2321-0

  • 81

    HelmM.MotorinY. (2017). Detecting Rna modifications in the epitranscriptome: Predict and validate.Nat. Rev. Genet.18275291. 10.1038/nrg.2016.169

  • 82

    HessM. E.HessS.MeyerK. D.VerhagenL. A.KochL.BronnekeH. S.et al (2013). The fat mass and obesity associated gene (Fto) regulates activity of the dopaminergic midbrain circuitry.Nat. Neurosci.1610421048. 10.1038/nn.3449

  • 83

    HuangY.HeN.ChenY.ChenZ.LiL. (2018). Bermp: A cross-species classifier for predicting m(6)A sites by integrating a deep learning algorithm and a random forest approach.Int. J. Biol. Sci.1416691677. 10.7150/ijbs.27819

  • 84

    HuangZ.PanJ.WangH.DuX.XuY.WangZ.et al (2021). Prognostic significance and tumor immune microenvironment heterogenicity of m5C Rna methylation regulators in triple-negative breast cancer.Front. Cell Dev. Biol.9:657547. 10.3389/fcell.2021.657547

  • 85

    HussainS.SajiniA. A.BlancoS.DietmannS.LombardP.SugimotoY.et al (2013). Nsun2-mediated cytosine-5 methylation of vault noncoding Rna determines its processing into regulatory small Rnas.Cell Rep.4255261. 10.1016/j.celrep.2013.06.029

  • 86

    ItoS.ShenL.DaiQ.WuS. C.CollinsL. B.SwenbergJ. A.et al (2011). Tet proteins can convert 5-methylcytosine to 5-Formylcytosine and 5-Carboxylcytosine.Science33313001303. 10.1126/science.1210597

  • 87

    JiaC. Z.ZhangJ. J.GuW. Z. (2016). Rna-MethylPred: A high-accuracy predictor to identify N6-methyladenosine in Rna.Anal. Biochem.5107275. 10.1016/j.ab.2016.06.012

  • 88

    JiaG.FuY.ZhaoX.DaiQ.ZhengG.YangY.et al (2011). N6-methyladenosine in nuclear Rna is a major substrate of the obesity-associated Fto.Nat. Chem. Biol.7885887. 10.1038/nchembio.687

  • 89

    JiaG.YangC. G.YangS.JianX.YiC.ZhouZ.et al (2008). Oxidative demethylation of 3-methylthymine and 3-methyluracil in single-stranded Dna and Rna by mouse and human Fto.FEBS Lett.58233133319. 10.1016/j.febslet.2008.08.019

  • 90

    JiangJ.SongB.ChenK.LuZ.RongR.ZhongY.et al (2022). m6AmPred: Identifying Rna N6, 2′-O-dimethyladenosine (m(6)Am) sites based on sequence-derived information.Methods203328334. 10.1016/j.ymeth.2021.01.007

  • 91

    JiangS.XieY.HeZ.ZhangY.ZhaoY.ChenL.et al (2018). m6asnp: A tool for annotating genetic variants by m6A function.Gigascience7:giy035. 10.1093/gigascience/giy035

  • 92

    JiangX.LiuB.NieZ.DuanL.XiongQ.JinZ.et al (2021). The role of m6A modification in the biological functions and diseases.Signal Transduct Target Ther.6:74. 10.1038/s41392-020-00450-x

  • 93

    JonkhoutN.TranJ.SmithM. A.SchonrockN.MattickJ. S.NovoaE. M. (2017). The Rna modification landscape in human disease.RNA2317541769. 10.1261/rna.063503.117

  • 94

    KarijolichJ.YuY. T. (2010). Spliceosomal snrna modifications and their function.RNA Biol.7192204. 10.4161/rna.7.2.11207

  • 95

    KhanA.RehmanH. U.HabibU.IjazU. (2020a). Detecting N6-methyladenosine sites from Rna transcriptomes using random forest.J. Comput. Sci.47:40242. 10.1016/j.jocs.2020.101238

  • 96

    KhanM. A.RafiqM. A.NoorA.HussainS.FloresJ. V.RuppV.et al (2012). Mutation in Nsun2, which encodes an Rna methyltransferase, causes autosomal-recessive intellectual disability.Am. J. Hum. Genet.90856863. 10.1016/j.ajhg.2012.03.023

  • 97

    KhanS. M.HeF.WangD.ChenY.XuD. (2020b). Mu-Pseudeep: A deep learning method for prediction of pseudouridine sites.Comput. Struct. Biotechnol. J.1818771883. 10.1016/j.csbj.2020.07.010

  • 98

    KiranA. M.O’mahonyJ. J.SanjeevK.BaranovP. V. (2013). Darned in 2013: Inclusion of model organisms and linking with Wikipedia.Nucleic Acids Res.41D258D261. 10.1093/nar/gks961

  • 99

    KnucklesP.LenceT.HaussmannI. U.JacobD.KreimN.CarlS. H.et al (2018). Zc3h13/Flacc is required for adenosine methylation by bridging the mrna-binding factor Rbm15/Spenito to the m(6)A machinery component Wtap/Fl(2)d.Genes Dev.32415429. 10.1101/gad.309146.117

  • 100

    KortelN.RuckleC.ZhouY.BuschA.Hoch-KraftP.SutandyF. X. R.et al (2021). Deep and accurate detection of m6A Rna modifications using miclip2 and m6Aboost machine learning.Nucleic Acids Res.49:e92. 10.1093/nar/gkab485

  • 101

    KriaucionisS.HeintzN. (2009). The Nuclear Dna Base 5-Hydroxymethylcytosine is present in purkinje neurons and the brain.Science324929930. 10.1126/science.1169786

  • 102

    KrugR.MorganM.ShatkinA. (1976). Influenza viral mrna con- tains internal N6 -methyladenosine and 5′-terminal 7-methyl- guanosine in cap structures.J. Virol.204553. 10.1128/Jvi.20.1.45-53.1976

  • 103

    LabrancheH.DupuisS.Ben-DavidY.BaniM. R.WellingerR. J.ChabotB. (1998). Telomere elongation by hnrnp A1 and a derivative that interacts with telomeric repeats and telomerase.Nat. Genet.19199202. 10.1038/575

  • 104

    LeachR.TuckM. (2001). Expression of the mrna (N6-adeno sine)-methyltransferase S-adenosyl-L-methionine binding subunit mrna in cultured cells.Int. J. Biochem. Cell Biol.33984999. 10.1016/s1357-2725(01)00071-1

  • 105

    LenceT.AkhtarJ.BayerM.SchmidK.SpindlerL.HoC. H.et al (2016). m6A modulates neuronal functions and sex determination in Drosophila.Nature540242247. 10.1038/nature20568

  • 106

    LiA.ChenY. S.PingX. L.YangX.XiaoW.YangY.et al (2017). Cytoplasmic m(6)A reader Ythdf3 promotes mrna translation.Cell Res.27444447. 10.1038/cr.2017.10

  • 107

    LiF.GuoX.JinP.ChenJ.XiangD.SongJ.et al (2021a). Porpoise: A new approach for accurate prediction of Rna pseudouridine sites.Brief Bioinform.22:bbab245. 10.1093/bib/bbab245

  • 108

    LiG. Q.LiuZ.ShenH. B.YuD. J. (2016a). TargetM6A: Identifying N(6)-methyladenosine sites from Rna sequences via position-specific nucleotide propensities and a support vector machine.IEEE Trans. Nanobiosci.15674682. 10.1109/Tnb.2016.2599115

  • 109

    LiJ.HuangY.YangX.ZhouY.ZhouY. (2018a). Rnam5Cfinder: A web-server for predicting Rna 5-methylcytosine (m5C) sites based on random forest.Sci. Rep.8:17299. 10.1038/s41598-018-35502-4

  • 110

    LiJ.YangX.ChuQ.XieL.DingY.XuX.et al (2022a). Multi-omics molecular biomarkers and database of osteoarthritis.Database (Oxford)2022:baac052. 10.1093/database/baac052

  • 111

    LiM.ZhaoX.WangW.ShiH.PanQ.LuZ.et al (2018b). Ythdf2-mediated m(6)A mrna clearance modulates neural development in mice.Genome Biol.19:69. 10.1186/s13059-018-1436-y

  • 112

    LiX. J.DuZ. W.ZarnowskaE. D.PankratzM.HansenL. O.PearceR. A.et al (2005). Specification of motoneurons from human embryonic stem cells.Nat. Biotechnol.23215221. 10.1038/nbt1063

  • 113

    LiX.XiongX.YiC. (2016b). Epitranscriptome sequencing technologies: Decoding Rna modifications.Nat. Methods142331. 10.1038/nmeth.4110

  • 114

    LiX.XiongX.WangK.WangL.ShuX.MaS.et al (2016c). Transcriptome-wide mapping reveals reversible and dynamic N1-methyladenosine methylome.Nat. Chem. Biol.12311316. 10.1038/nchembio.2040

  • 115

    LiX.ZhangS.ShiH. (2022b). An improved residual network using deep fusion for identifying Rna 5-methylcytosine sites.Bioinformatics3842714277. 10.1093/bioinformatics/btac532

  • 116

    LiX.ZhuP.MaS.SongJ.BaiJ.SunF.et al (2015a). Chemical pulldown reveals dynamic pseudouridylation of the mammalian transcriptome.Nat. Chem. Biol.11592597. 10.1038/nchembio.1836

  • 117

    LiY.QiD.ZhuB.YeX. (2021b). Analysis of m6A Rna methylation-related genes in liver hepatocellular carcinoma and their correlation with survival.Int. J. Mol. Sci.22:1474. 10.3390/ijms22031474

  • 118

    LiY.-H.ZhangG.CuiQ. (2015b). Ppus: A web server to predict Pus-specific pseudouridine sites: Table 1.Bioinformatics3133623364. 10.1093/bioinformatics/btv366

  • 119

    LiaoJ.WangQ.WuF.HuangZ. (2022). In silico methods for identification of potential active sites of therapeutic targets.Molecules27:7103. 10.3390/molecules27207103

  • 120

    LinZ.HsuP. J.XingX.FangJ.LuZ.ZouQ.et al (2017). Mettl3-/Mettl14-mediated mrna N(6)-methyladenosine modulates murine spermatogenesis.Cell Res.2712161230. 10.1038/cr.2017.117

  • 121

    LiuH.WangH.WeiZ.ZhangS.HuaG.ZhangS. W.et al (2018). MeT-Db V2.0: Elucidating context-specific functions of N6-methyl-adenosine methyltranscriptome.Nucleic Acids Res.46D281D287. 10.1093/nar/gkx1080

  • 122

    LiuJ.LiK.CaiJ.ZhangM.ZhangX.XiongX.et al (2020a). Landscape and regulation of m(6)A and m(6)Am methylome across human and mouse tissues.Mol. Cell77426440e6. 10.1016/j.molcel.2019.09.032

  • 123

    LiuJ.SunG.PanS.QinM.OuyangR.LiZ.et al (2020b). The cancer genome atlas (Tcga) based m(6)A methylation-related genes predict prognosis in hepatocellular carcinoma.Bioengineered11759768. 10.1080/21655979.2020.1787764

  • 124

    LiuJ.YueY.HanD.WangX.FuY.ZhangL.et al (2013). A Mettl3–Mettl14 complex mediates mammalian nuclear Rna N6-adenosine methylation.Nat. Chem. Biol.109395. 10.1038/nchembio.1432

  • 125

    LiuK.ChenW. (2020). imrm: A platform for simultaneously identifying multiple kinds of Rna modifications.Bioinformatics3633363342. 10.1093/bioinformatics/btaa155

  • 126

    LiuK.CaoL.DuP.ChenW. (2020c). im6A-Ts-Cnn: Identifying the N(6)-Methyladenine site in multiple tissues by using the convolutional neural network.Mol. Ther. Nucleic Acids2110441049. 10.1016/j.omtn.2020.07.034

  • 127

    LiuK.ChenW.LinH. (2020d). Xg-PseU: An eXtreme Gradient Boosting based method for identifying pseudouridine sites.Mol. Genet. Genomics2951321. 10.1007/s00438-019-01600-9

  • 128

    LiuL.LeiX.FangZ.TangY.MengJ.WeiZ. (2020e). Lithophone: Improving lncrna methylation site prediction using an ensemble predictor.Front. Genet.11:545. 10.3389/fgene.2020.00545

  • 129

    LiuL.LeiX.MengJ.WeiZ. (2020f). Isgm1A: Integration of sequence features and genomic features to improve the prediction of human m1A Rna methylation sites.IEEE Access.88197181977. 10.1109/access.2020.2991070.

  • 130

    LiuL.LeiX.MengJ.WeiZ. (2020g). Witmsg: Large-scale prediction of human intronic m(6)A Rna methylation sites from sequence and genomic features.Curr. Genomics216776. 10.2174/1389202921666200211104140

  • 131

    LiuN.DaiQ.ZhengG.HeC.ParisienM.PanT. (2015). N(6)-methyladenosine-dependent Rna structural switches regulate Rna-protein interactions.Nature518560564. 10.1038/nature14234

  • 132

    LiuN.ZhouK. I.ParisienM.DaiQ.DiatchenkoL.PanT. (2017). N6-methyladenosine alters Rna structure to regulate binding of a low-complexity protein.Nucleic Acids Res.4560516063. 10.1093/nar/gkx141

  • 133

    LiuQ.GregoryR. I. (2019). Rnamod: An integrated system for the annotation of mrna modifications.Nucleic Acids Res.47W548W555. 10.1093/nar/gkz479

  • 134

    LiuS.ZhuA.HeC.ChenM. (2020h). Repic: A database for exploring the N(6)-methyladenosine methylome.Genome Biol.21:100. 10.1186/s13059-020-02012-4

  • 135

    LiuY.ChenD.SuR.ChenW.WeiL. (2020i). irna5hmC: The first predictor to identify Rna 5-hydroxymethylcytosine modifications using machine learning.Front. Bioeng. Biotechnol.8:227. 10.3389/fbioe.2020.00227

  • 136

    LiuZ.XiaoX.YuD. J.JiaJ.QiuW. R.ChouK. C. (2016). prnam-Pc: Predicting N(6)-methyladenosine sites in Rna sequences via physical-chemical properties.Anal. Biochem.4976067. 10.1016/j.ab.2015.12.017

  • 137

    Lo GiudiceC.PesoleG.PicardiE. (2018). Redidb 3.0: A comprehensive collection of Rna editing events in plant organellar genomes.Front. Plant Sci.9:482. 10.3389/fpls.2018.00482

  • 138

    LuoG. Z.MacqueenA.ZhengG.DuanH.DoreL. C.LuZ.et al (2014). Unique features of the m6A methylome in Arabidopsis thaliana.Nat. Commun.5:5630. 10.1038/ncomms6630

  • 139

    LuoX.LiH.LiangJ.ZhaoQ.XieY.RenJ.et al (2021). Rmvar: An updated database of functional variants involved in Rna modifications.Nucleic Acids Res.49D1405D1412. 10.1093/nar/gkaa811

  • 140

    LuoZ.SuW.LouL.QiuW.XiaoX.XuZ. (2022). Dlm6Am: A deep-learning-based tool for identifying N6,2′-O-Dimethyladenosine sites in rna sequences.Int. J. Mol. Sci.23:11026. 10.3390/ijms231911026

  • 141

    LvH.DaoF. Y.ZhangD.YangH.LinH. (2021). Advances in mapping the epigenetic modifications of 5-methylcytosine (5mC), N6-methyladenine (6mA), and N4-methylcytosine (4mC).Biotechnol. Bioeng.11842044216. 10.1002/bit.27911

  • 142

    LvH.ZhangZ. M.LiS. H.TanJ. X.ChenW.LinH. (2020a). Evaluation of different computational methods on 5-methylcytosine sites identification.Brief Bioinform.21982995. 10.1093/bib/bbz048

  • 143

    LvZ.ZhangJ.DingH.ZouQ. (2020b). Rf-PseU: A random forest predictor for rna pseudouridine sites.Front. Bioeng. Biotechnol.8:134. 10.3389/fbioe.2020.00134

  • 144

    MaC.ChangM.LvH.ZhangZ. W.ZhangW.HeX.et al (2018). Rna m(6)A methylation participates in regulation of postnatal development of the mouse cerebellum.Genome Biol.19:68. 10.1186/s13059-018-1435-z

  • 145

    MaJ.SongB.WeiZ.HuangD.ZhangY.SuJ.et al (2022). m5C-Atlas: A comprehensive database for decoding and annotating the 5-methylcytosine (m5C) epitranscriptome.Nucleic Acids Res.50D196D203. 10.1093/nar/gkab1075

  • 146

    MahmoudiO.WahabA.ChongK. T. (2020). iMethyl-Deep: N6 methyladenosine identification of yeast genome with automatic feature extraction technique by using deep learning algorithm.Genes (Basel)11:529. 10.3390/genes11050529

  • 147

    MassignaniE.GiambrunoR.ManiaciM.NicosiaL.YadavA.CuomoA.et al (2022). ProMetheusdb: An in-depth analysis of the high-quality human Methyl-proteome.Mol. Cell Proteomics21:100243. 10.1016/j.mcpro.2022.100243

  • 148

    MauerJ.LuoX.BlanjoieA.JiaoX.GrozhikA. V.PatilD. P.et al (2016). Reversible methylation of m6Am in the 5′ cap controls mrna stability.Nature541371375. 10.1038/nature21022

  • 149

    MeyerK. D.SaletoreY.ZumboP.ElementoO.MasonC. E.JaffreyS. R. (2012). Comprehensive analysis of mrna methylation reveals enrichment in 3′. Utrs and near stop codons.Cell14916351646. 10.1016/j.cell.2012.05.003

  • 150

    MolinieB.WangJ.LimK. S.HillebrandR.LuZ. X.Van WittenbergheN.et al (2016). m(6)A-Laic-seq reveals the census and complexity of the m(6)A epitranscriptome.Nat. Methods13692698. 10.1038/nmeth.3898

  • 151

    MuY.ZhangR.WangL.LiuX. (2020). iPseU-Layer: Identifying Rna pseudouridine sites using layered ensemble model.Interdiscip. Sci.12193203. 10.1007/s12539-020-00362-y

  • 152

    NachtergaeleS.HeC. (2017). The emerging biology of Rna post-transcriptional modifications.RNA Biol.14156163. 10.1080/15476286.2016.1267096

  • 153

    NazariI.TahirM.TayaraH.ChongK. T. (2019). iN6-Methyl (5-step): Identifying Rna N6-methyladenosine sites using deep learning mode via Chou’s 5-step rules and Chou’s general Pseknc.Chemometr. Intell. Lab. Syst.193:103811. 10.1016/j.chemolab.2019.103811

  • 154

    Nguyen-VoT.-H.NguyenQ. H.DoT. T. T.NguyenT.-N.RahardjaS.NguyenB. P. (2019). iPseU-Ncp: Identifying Rna pseudouridine sites using random forest and Ncp-encoded features.BMC Genomics20:971. 10.1186/s12864-019-6357-y

  • 155

    NiW.YaoS.ZhouY.LiuY.HuangP.ZhouA.et al (2019). Long noncoding Rna Gas5 inhibits progression of colorectal cancer by interacting with and triggering Yap phosphorylation and degradation and is negatively regulated by the m(6)A reader Ythdf3.Mol. Cancer18:143. 10.1186/s12943-019-1079-y

  • 156

    NieF.FengP.SongX.WuM.TangQ.ChenW. (2020). Rnawre: A resource of writers, readers and erasers of Rna modifications.Database (Oxford)2020:baaa049. 10.1093/database/baaa049

  • 157

    NombelaP.Miguel-LopezB.BlancoS. (2021). The role of m(6)A, m(5)C and Psi Rna modifications in cancer: Novel therapeutic opportunities.Mol. Cancer20:18. 10.1186/s12943-020-01263-w

  • 158

    OerumS.MeynierV.CatalaM.TisneC. (2021). A comprehensive review of m6A/m6Am Rna methyltransferase structures.Nucleic Acids Res.4972397255. 10.1093/nar/gkab378

  • 159

    OrtegaA.NiksicM.BachiA.WilmM.SanchezL.HastieN.et al (2003). Biochemical function of female-lethal (2)D/Wilms’ tumor suppressor-1-associated proteins in alternative pre-mrna splicing.J. Biol. Chem.27830403047. 10.1074/jbc.M210737200

  • 160

    PanH.TanP. F.LimI. Y.HuanJ.TehA. L.ChenL.et al (2022). Integrative multi-omics database (imomdb) of Asian pregnant women.Hum. Mol. Genet.3130513067. 10.1093/hmg/ddac079

  • 161

    PanwarB.RaghavaG. P. S. (2014). Prediction of uridine modifications in trna sequences.BMC Bioinform.15:326. 10.1186/1471-2105-15-326

  • 162

    PatilD. P.ChenC. K.PickeringB. F.ChowA.JacksonC.GuttmanM.et al (2016). m(6)A Rna methylation promotes Xist-mediated transcriptional repression.Nature537369373. 10.1038/nature19342

  • 163

    PendletonK. E.ChenB.LiuK.HunterO. V.XieY.TuB. P.et al (2017). The U6 snrna m(6)A Methyltransferase Mettl16 regulates sam synthetase intron retention.Cell169824835e14. 10.1016/j.cell.2017.05.003

  • 164

    PingX. L.SunB. F.WangL.XiaoW.YangX.WangW. J.et al (2014b). Mammalian Wtap is a regulatory subunit of the Rna N6-methyladenosine methyltransferase.Cell Res.24177189.

  • 165

    PingX.-L.SunB.-F.WangL.XiaoW.YangX.WangW.-J.et al (2014a). Mammalian Wtap is a regulatory subunit of the Rna N6-methyladenosine methyltransferase.Cell Res.24177189. 10.1038/cr.2014.3

  • 166

    QiangX.ChenH.YeX.SuR.WeiL. (2018). M6amrfs: Robust prediction of N6-methyladenosine sites with sequence-based features in multiple species.Front. Genet.9:495. 10.3389/fgene.2018.00495

  • 167

    QiuW. R.JiangS. Y.XuZ. C.XiaoX.ChouK. C. (2017). irnam5C-Psednc: Identifying Rna 5-methylcytosine sites by incorporating physical-chemical properties into pseudo dinucleotide composition.Oncotarget84117841188. 10.18632/oncotarget.17104

  • 168

    RehmanM. U.HongK. J.TayaraH.ChongK. T. (2021a). m6A-NeuralTool: Convolution neural tool for Rna N6-methyladenosine site identification in different species.IEEE Access.91777917786. 10.1109/access.2021.3054361

  • 169

    RehmanM. U.TayaraH.ChongK. T. (2021b). Dcnn-4mC: Densely connected neural network based N4-methylcytosine site prediction in multiple species.Comput. Struct. Biotechnol. J.1960096019. 10.1016/j.csbj.2021.10.034

  • 170

    RehmanM. U.TayaraH.ChongK. T. (2022a). “Dl-M6A: Identification of N6-methyladenosine sites in mammals using deep learning based on different encoding schemes,” in IEEE/ACM Transactions on computational biology and bioinformatics, (Piscataway: IEEE). 10.1109/Tcbb.2022.3192572

  • 171

    RehmanM. U.TayaraH.ZouQ.ChongK. T. (2022b). i6mA-Caps: A CapsuleNet-based framework for identifying Dna N6-methyladenine sites.Bioinformatics3838853891. 10.1093/bioinformatics/btac434

  • 172

    RiesR. J.ZaccaraS.KleinP.Olarerin-GeorgeA.NamkoongS.PickeringB. F.et al (2019). m(6)A enhances the phase separation potential of mrna.Nature571424428. 10.1038/s41586-019-1374-1

  • 173

    RoostC.LynchS. R.BatistaP. J.QuK.ChangH. Y.KoolE. T. (2015). Structure and thermodynamics of N6-methyladenosine in Rna: A spring-loaded base modification.J. Am. Chem. Soc.13721072115. 10.1021/ja513080v

  • 174

    RoundtreeI. A.EvansM. E.PanT.HeC. (2017a). Dynamic Rna modifications in gene expression regulation.Cell16911871200. 10.1016/j.cell.2017.05.045

  • 175

    RoundtreeI. A.LuoG. Z.ZhangZ.WangX.ZhouT.CuiY.et al (2017b). Ythdc1 mediates nuclear export of N(6)-methyladenosine methylated mrnas.Elife6:e31311. 10.7554/eLife.31311

  • 176

    RyvkinP.LeungY. Y.SilvermanI. M.ChildressM.ValladaresO.DragomirI.et al (2013). Hamr: High-throughput annotation of modified ribonucleotides.RNA1916841692. 10.1261/rna.036806.112

  • 177

    SaboohM. F.IqbalN.KhanM.KhanM.MaqboolH. F. (2018). Identifying 5-methylcytosine sites in Rna sequence using composite encoding feature into Chou’s Pseknc.J. Theory Biol.45219. 10.1016/j.jtbi.2018.04.037

  • 178

    SajekM. P.WozniakT.SprinzlM.JaruzelskaJ.BarciszewskiJ. (2020). T-psi-C: User friendly database of trna sequences and structures.Nucleic Acids Res.48D256D260. 10.1093/nar/gkz922

  • 179

    SchwartzS.AgarwalaS. D.MumbachM. R.JovanovicM.MertinsP.ShishkinA.et al (2013). High-resolution mapping reveals a conserved, widespread, dynamic mrna methylation program in yeast meiosis.Cell15514091421. 10.1016/j.cell.2013.10.047

  • 180

    SchwartzS.BernsteinD. A.MumbachM. R.JovanovicM.HerbstR. H.Leon-RicardoB. X.et al (2014a). Transcriptome-wide mapping reveals widespread dynamic-regulated pseudouridylation of ncrna and mrna.Cell159148162. 10.1016/j.cell.2014.08.028

  • 181

    SchwartzS.MumbachM. R.JovanovicM.WangT.MaciagK.BushkinG. G.et al (2014b). Perturbation of m6A writers reveals two distinct classes of mrna methylation at Internal and 5′ Sites.Cell Rep.8284296. 10.1016/j.celrep.2014.05.048

  • 182

    SchwartzS.MumbachM.JovanovicM.WangT.MaciagK.BushkinG.et al (2014). Perturbation of m6A writers reveals two distinct classes of mrna methylation at internal and 5′ sites.Cell Rep.8284296.

  • 183

    ShiH.ZhangS.LiX. (2022). R5hmcfdv: Computational identification of Rna 5-hydroxymethylcytosine based on deep feature fusion and deep voting.Brief Bioinform.23:bbac341. 10.1093/bib/bbac341

  • 184

    ShiM.ZhangH.WuX.HeZ.WangL.YinS.et al (2017). Alyref mainly binds to the 5′ and the 3′ regions of the mrna in vivo.Nucleic Acids Res.4596409653. 10.1093/nar/gkx597

  • 185

    SohY. Q. S.MikedisM. M.KojimaM.GodfreyA. K.De RooijD. G.PageD. C. (2017). Meioc maintains an extended meiotic prophase I in mice.PLoS Genet.13:e1006704. 10.1371/journal.pgen.1006704

  • 186

    SongB.ChenK.TangY.MaJ.MengJ.WeiZ. (2020a). Psi-Mouse: Predicting mouse pseudouridine sites from sequence and genome-derived features.Evol. Bioinform.16:1176934320925752. 10.1177/1176934320925752

  • 187

    SongB.ChenK.TangY.WeiZ.SuJ.De MagalhaesJ. P.et al (2021a). Consrm: Collection and large-scale prediction of the evolutionarily conserved Rna methylation sites, with implications for the functional epitranscriptome.Brief Bioinform.22:bbab088. 10.1093/bib/bbab088

  • 188

    SongB.HuangD.ZhangY.WeiZ.SuJ.Pedro de MagalhãesJ.et al (2022a). m6A-TSHub: Unveiling the context-specific m6A methylation and m6A-affecting mutations in 23 human tissues. Genom. Proteom. Bioinform.10.1016/j.gpb.2022.09.001

  • 189

    SongB.WangX.LiangZ.MaJ.HuangD.WangY.et al (2022b). Rmdisease V2.0: An updated database of genetic variants that affect Rna modifications with disease and trait implication.Nucleic Acids Res.10.1093/nar/gkac750

  • 190

    SongB.TangY.WeiZ.LiuG.SuJ.MengJ.et al (2020b). Piano: A web server for Pseudouridine-Site (Psi) identification and functional annotation.Front. Genet.11:88. 10.3389/fgene.2020.00088

  • 191

    SongJ.ZhaiJ.BianE.SongY.YuJ.MaC. (2018). Transcriptome-wide annotation of m(5)C Rna modifications using machine learning.Front. Plant Sci.9:519. 10.3389/fpls.2018.00519

  • 192

    SongZ.HuangD.SongB.ChenK.SongY.LiuG.et al (2021b). Attention-based multi-label neural networks for integrated prediction and interpretation of twelve widely occurring Rna modifications.Nat. Commun.12:4011. 10.1038/s41467-021-24313-3

  • 193

    SquiresJ. E.PatelH. R.NouschM.SibbrittT.HumphreysD. T.ParkerB. J.et al (2012). Widespread occurrence of 5-methylcytosine in human coding and non-coding Rna.Nucleic Acids Res.4050235033. 10.1093/nar/gks144

  • 194

    SulemanM. T.KhanY. D. (2022). m1A-pred: Prediction of Modified 1-methyladenosine Sites in Rna sequences through artificial intelligence.Comb. Chem. High. Throughput Screen2524732484. 10.2174/1386207325666220617152743

  • 195

    SunP. P.ChenY. B.LiuB.GaoY. X.HanY.HeF.et al (2019). Deepmrmp: A new predictor for multiple types of Rna modification sites using deep learning.Math. Biosci. Eng.1662316241. 10.3934/mbe.2019310

  • 196

    TahirM.HayatM.ChongK. T. (2020). Prediction of N6-methyladenosine sites using convolution neural network model based on distributed feature representations.Neural Netw.129385391. 10.1016/j.neunet.2020.05.027

  • 197

    TahirM.HayatM.ChongK. T. (2021). A convolution neural network-based computational model to identify the occurrence sites of various Rna modifications by fusing varied features.Chemometr. Intell. Lab. Syst.211:104233. 10.1016/j.chemolab.2021.104233

  • 198

    TahirM.TayaraH.ChongK. T. (2019). iPseU-Cnn: Identifying Rna pseudouridine sites using convolutional neural networks.Mol. Ther. Nucleic Acids16463470. 10.1016/j.omtn.2019.03.010

  • 199

    TanabeA.TanikawaK.TsunetomiM.TakaiK.IkedaH.KonnoJ.et al (2016). Rna helicase Ythdc2 promotes cancer metastasis via the enhancement of the efficiency by which Hif-1alpha mrna is translated.Cancer Lett.3763442. 10.1016/j.canlet.2016.02.022

  • 200

    TangG.ChoM.WangX. (2022). Oncodb: An interactive online database for analysis of gene expression and viral infection in cancer.Nucleic Acids Res.50D1334D1339. 10.1093/nar/gkab970

  • 201

    TangY.ChenK.SongB.MaJ.WuX.XuQ.et al (2021). m6A-Atlas: A comprehensive knowledgebase for unraveling the N6-methyladenosine (m6A) epitranscriptome.Nucleic Acids Res.49D134D143. 10.1093/nar/gkaa692

  • 202

    TrixlL.AmortT.WilleA.ZinniM.EbnerS.HechenbergerC.et al (2018). Rna cytosine methyltransferase Nsun3 regulates embryonic stem cell differentiation by promoting mitochondrial activity.Cell Mol. Life Sci.7514831497. 10.1007/s00018-017-2700-0

  • 203

    TuckM. T. (1992). Partial purification of a 6-methyladenine mrna methyltransferase which modifies internal adenine residues.Biochem. J.288 (Pt 1)233240. 10.1042/bj2880233

  • 204

    UyarB.YusufD.WurmusR.RajewskyN.OhlerU.AkalinA. (2017). Rcas: An Rna centric annotation system for transcriptome-wide regions of interest.Nucleic Acids Res.45:e91. 10.1093/nar/gkx120

  • 205

    Van HauteL.LeeS. Y.MccannB. J.PowellC. A.BansalD.VasiliauskaiteL.et al (2019). Nsun2 introduces 5-methylcytosines in mammalian mitochondrial trnas.Nucleic Acids Res.4787208733. 10.1093/nar/gkz559

  • 206

    VandivierL. E.GregoryB. D. (2017). Reading the Epitranscriptome: New techniques and perspectives.Enzymes41269298. 10.1016/bs.enz.2017.03.004

  • 207

    WalkleyC. R.LiJ. B. (2017). Rewriting the transcriptome: Adenosine-to-inosine Rna editing by Adars.Genome Biol.18:205. 10.1186/s13059-017-1347-3

  • 208

    WangC. X.CuiG. S.LiuX.XuK.WangM.ZhangX. X.et al (2018a). Mettl3-mediated m6A modification is required for cerebellar development.PLoS Biol.16:e2004880. 10.1371/journal.pbio.2004880

  • 209

    WangH.LiuH.HuangT.LiG.ZhangL.SunY. (2022a). Emdlp: Ensemble multiscale deep learning model for Rna methylation site prediction.BMC Bioinform.23:221. 10.1186/s12859-022-04756-1

  • 210

    WangH.WangS.ZhangY.BiS.ZhuX. (2022b). A brief review of machine learning methods for Rna methylation sites prediction.Methods203399421. 10.1016/j.ymeth.2022.03.001

  • 211

    WangJ.ZhangS. (2021). Pa-PseU: An incremental passive-aggressive based method for identifying Rna pseudouridine sites via Chou’s 5-steps rule.Chemometr. Intell. Lab. Syst.210:104250. 10.1016/j.chemolab.2021.104250

  • 212

    WangX.YanR. (2018). RfathM6A: A new tool for predicting m(6)A sites in Arabidopsis thaliana.Plant Mol. Biol.96327337. 10.1007/s11103-018-0698-9

  • 213

    WangX.LuZ.GomezA.HonG. C.YueY.HanD.et al (2014). N6-methyladenosine-dependent regulation of messenger Rna stability.Nature505117120. 10.1038/nature12730

  • 214

    WangX.ZhaoB. S.RoundtreeI. A.LuZ.HanD.MaH.et al (2015). N6-methyladenosine modulates messenger Rna translation efficiency.Cell16113881399. 10.1016/j.cell.2015.05.014

  • 215

    WangY.GuoR.HuangL.YangS.HuX.HeK. (2021). m6age: A Predictor for N6-methyladenosine sites identification utilizing sequence characteristics and graph embedding-based geometrical information.Front. Genet.12:670852. 10.3389/fgene.2021.670852

  • 216

    WangY.LiY.YueM.WangJ.KumarS.Wechsler-ReyaR. J.et al (2018b). N(6)-methyladenosine Rna modification regulates embryonic neural stem cell self-renewal through histone modifications.Nat. Neurosci.21195206. 10.1038/s41593-017-0057-1

  • 217

    WangZ.JensenM. A.ZenklusenJ. C. (2016). A practical guide to the cancer genome atlas (Tcga).Methods Mol. Biol.1418111141. 10.1007/978-1-4939-3578-9_6

  • 218

    WeiC.GershowitzA.MossB. (1975). N6, O2′-dimethyladenosine a novel methylated ribonucleoside next to the 5′ terminal of animal cell and virus mrnas.Nature257251253. 10.1038/257251a0

  • 219

    WeiL.ChenH.SuR. (2018). M6apred-El: A sequence-based predictor for identifying n6-methyladenosine sites using ensemble learning.Mol. Ther. Nucleic Acids12635644. 10.1016/j.omtn.2018.07.004

  • 220

    WeiL.SuR.WangB.LiX.ZouQ.GaoX. (2019). Integration of deep feature representations and handcrafted features to improve the prediction of N6-methyladenosine sites.Neurocomputing32439. 10.1016/j.neucom.2018.04.082

  • 221

    WenJ.LvR.MaH.ShenH.HeC.WangJ.et al (2018). Zc3h13 Regulates nuclear Rna m(6)A methylation and mouse embryonic stem cell self-renewal.Mol. Cell6910281038e6. 10.1016/j.molcel.2018.02.015

  • 222

    WnukM.SlipekP.DziedzicM.LewinskaA. (2020). The Roles of Host 5-methylcytosine Rna methyltransferases during viral infections.Int. J. Mol. Sci.21:8176. 10.3390/ijms21218176

  • 223

    WojtasM. N.PandeyR. R.MendelM.HomolkaD.SachidanandamR.PillaiR. S. (2017). Regulation of m(6)A Transcripts by the 3′–>5′ Rna Helicase Ythdc2 Is Essential for a successful meiotic program in the mammalian germline.Mol. Cell68374387e12. 10.1016/j.molcel.2017.09.021

  • 224

    WuB.SuS.PatilD. P.LiuH.GanJ.JaffreyS. R.et al (2018). Molecular basis for the specific and multivariant recognitions of Rna substrates by human hnrnp A2/B1.Nat. Commun.9:420. 10.1038/s41467-017-02770-z

  • 225

    XiangS.LiuK.YanZ.ZhangY.SunZ. (2016a). RnamethPre: A web server for the prediction and query of mrna m6A Sites.PLoS One11:e0162707. 10.1371/journal.pone.0162707

  • 226

    XiangS.YanZ.LiuK.ZhangY.SunZ. (2016b). AthMethPre: A web server for the prediction and query of mrna m(6)A sites in Arabidopsis thaliana.Mol. Biosyst.1233333337. 10.1039/c6mb00536e

  • 227

    XiaoW.AdhikariS.DahalU.ChenY. S.HaoY. J.SunB. F.et al (2016). Nuclear m(6)A reader ythdc1 regulates mrna splicing.Mol. Cell61507519. 10.1016/j.molcel.2016.01.012

  • 228

    XingP.SuR.GuoF.WeiL. (2017). Identifying N(6)-methyladenosine sites using multi-interval nucleotide pair position specificity and support vector machine.Sci. Rep.7:46757. 10.1038/srep46757

  • 229

    XuC.LiuK.AhmedH.LoppnauP.SchapiraM.MinJ. (2015). Structural Basis for the Discriminative Recognition of N6-Methyladenosine Rna by the Human Yt521-B homology domain family of proteins.J. Biol. Chem.2902490224913. 10.1074/jbc.M115.680389

  • 230

    XuC.WangX.LiuK.RoundtreeI. A.TempelW.LiY.et al (2014). Structural basis for selective binding of m6A Rna by the Ythdc1 Yth domain.Nat. Chem. Biol.10927929. 10.1038/nchembio.1654

  • 231

    XuanJ. J.SunW. J.LinP. H.ZhouK. R.LiuS.ZhengL. L.et al (2018). Rmbase v2.0: Deciphering the map of Rna modifications from epitranscriptome sequencing data.Nucleic Acids Res.46D327D334. 10.1093/nar/gkx934

  • 232

    YangX.YangY.SunB. F.ChenY. S.XuJ. W.LaiW. Y.et al (2017a). 5-methylcytosine promotes mrna export - Nsun2 as the methyltransferase and Alyref as an m(5)C reader.Cell Res.27606625. 10.1038/cr.2017.55

  • 233

    YangY.HsuP. J.ChenY.-S.YangY.-G. (2017b). 5-methylcytosine promotes mrna export — Nsun2 as the methyltransferase and Alyref as an m5C reader.Cell Res.28616624. 10.1038/s41422-018-0040-8

  • 234

    YaoJ. H.LinM. X.LiaoW. J.FanW. J.XuX. X.ShiH.et al (2022). M1aregpred: Epitranscriptome target prediction of N1-methyladenosine (m1A) regulators based on sequencing features and genomic features.Front. Biosci.27:269. 10.31083/j.fbl2709269

  • 235

    YaoW.HanX.GeM.ChenC.XiaoX.LiH.et al (2020). N(6)-methyladenosine (m(6)A) methylation in ischemia-reperfusion injury.Cell Death Dis.11:478. 10.1038/s41419-020-2686-7

  • 236

    YoonK. J.RingelingF. R.VissersC.JacobF.PokrassM.Jimenez-CyrusD.et al (2017). Temporal control of mammalian cortical neurogenesis by m(6)A METHYLATION.Cell171877889e17. 10.1016/j.cell.2017.09.003

  • 237

    YuJ.ChenM.HuangH.ZhuJ.SongH.ZhuJ.et al (2018). Dynamic m6A modification regulates local translation of mrna in axons.Nucleic Acids Res.4614121423. 10.1093/nar/gkx1182

  • 238

    YuJ.XingX.ZengL.SunJ.LiW.SunH.et al (2012). SyStemCell: A database populated with multiple levels of experimental data from stem cell differentiation research.PLoS One7:e35230. 10.1371/journal.pone.0035230

  • 239

    YueY.LiuJ.CuiX.CaoJ.LuoG.ZhangZ.et al (2018). Virma mediates preferential m(6)A mrna methylation in 3′utr and near stop codon and associates with alternative polyadenylation.Cell Discov.4:10. 10.1038/s41421-018-0019-0

  • 240

    ZhangJ.FengP.LinH.ChenW. (2018a). Identifying Rna N(6)-methyladenosine sites in Escherichia coli Genome.Front. Microbiol.9:955. 10.3389/fmicb.2018.00955

  • 241

    ZhangL.LiG.LiX.WangH.ChenS.LiuH. (2021a). Edlm(6)Apred: Ensemble deep learning approach for mrna m(6)A site prediction.BMC Bioinform.22:288. 10.1186/s12859-021-04206-4

  • 242

    ZhangL.QinX.LiuM.XuZ.LiuG. (2021b). Dnn-m6A: A cross-species method for identifying Rna N6-methyladenosine sites based on deep neural network with multi-information fusion.Genes (Basel)12:354. 10.3390/genes12030354

  • 243

    ZhangS.WangJ.LiX.LiangY. (2021c). M6A-Gsms: Computational identification of N(6)-methyladenosine sites with Gbdt and stacking learning in multiple species.J. Biomol. Struct. Dyn.112. 10.1080/07391102.2021.1970628

  • 244

    ZhangM.SunJ.-W.LiuZ.RenM.-W.ShenH.-B.YuD.-J. (2016). Improving N6-methyladenosine site prediction with heuristic selection of nucleotide physical–chemical properties.Anal. Biochem.508104113. 10.1016/j.ab.2016.06.001

  • 245

    ZhangM.XuY.LiL.LiuZ.YangX.YuD. J. (2018b). Accurate Rna 5-methylcytosine site prediction based on heuristic physical-chemical properties reduction and classifier ensemble.Anal. Biochem.5504148. 10.1016/j.ab.2018.03.027

  • 246

    ZhangS. Y.ZhangS. W.FanX. N.ZhangT.MengJ.HuangY. (2019). Fundmdeep-m6A: Identification and prioritization of functional differential m6A methylation genes.Bioinformatics35i90i98. 10.1093/bioinformatics/btz316

  • 247

    ZhangS.ShiH. (2021). iR5hmcsc: Identifying Rna 5-hydroxymethylcytosine with multiple features based on stacking learning.Comput. Biol. Chem.95:107583. 10.1016/j.compbiolchem.2021.107583

  • 248

    ZhangW.TanX.LinS.GouY.HanC.ZhangC.et al (2022). Cplm 4.0: An updated database with rich annotations for protein lysine modifications.Nucleic Acids Res.50D451D459. 10.1093/nar/gkab849

  • 249

    ZhangY.HamadaM. (2018). DeepM6aseq: Prediction and characterization of m6A-containing sequences using deep learning.BMC Bioinform.19:524. 10.1186/s12859-018-2516-4

  • 250

    ZhangZ.ThelerD.KaminskaK. H.HillerM.De La GrangeP.PudimatR.et al (2010). The Yth domain is a novel Rna binding domain.J. Biol. Chem.2851470114710. 10.1074/jbc.M110.104711

  • 251

    ZhaoB. S.RoundtreeI. A.HeC. (2017). Post-transcriptional gene regulation by mrna modifications.Nat. Rev. Mol. Cell Biol.183142. 10.1038/nrm.2016.132

  • 252

    ZhaoL.-Y.SongJ.LiuY.SongC.-X.YiC. (2020). Mapping the epigenetic modifications of Dna and Rna.Protein Cell11792808. 10.1007/s13238-020-00733-7

  • 253

    ZhaoT. T.GraberT. E.JordanL. E.CloutierM.LewisS. M.GouletI.et al (2009). hnrnp A1 regulates Uv-induced Nf-kappaB signalling through destabilization of ciap1 mrna.Cell Death Differ.16244252. 10.1038/cdd.2008.146

  • 254

    ZhaoX.YangY.SunB. F.ShiY.YangX.XiaoW.et al (2014). Fto-dependent demethylation of N6-methyladenosine regulates mrna splicing and is required for adipogenesis.Cell Res.2414031419. 10.1038/cr.2014.151

  • 255

    ZhaoZ.PengH.LanC.ZhengY.FangL.LiJ. (2018). Imbalance learning for the prediction of N(6)-Methylation sites in mrnas.BMC Genomics19:574. 10.1186/s12864-018-4928-y

  • 256

    ZhenD.WuY.ZhangY.ChenK.SongB.XuH.et al (2020). m(6)A Reader: Epitranscriptome target prediction and functional characterization of N (6)-Methyladenosine (m(6)A) Readers.Front. Cell Dev. Biol.8:741. 10.3389/fcell.2020.00741

  • 257

    ZhengG.DahlJ. A.NiuY.FedorcsakP.HuangC. M.LiC. J.et al (2013). Alkbh5 is a mammalian Rna demethylase that impacts Rna metabolism and mouse fertility.Mol. Cell491829. 10.1016/j.molcel.2012.10.015

  • 258

    ZhengY.NieP.PengD.HeZ.LiuM.XieY.et al (2018). m6avar: A database of functional variants involved in m6A modification.Nucleic Acids Res.46D139D145. 10.1093/nar/gkx895

  • 259

    ZhongS.LiH.BodiZ.ButtonJ.VespaL.HerzogM.et al (2008). Mta is an Arabidopsis messenger Rna adenosine methylase and interacts with a homolog of a sex-specific splicing factor.Plant Cell2012781288. 10.1105/tpc.108.058883

  • 260

    ZhouQ.GuanP.ZhuZ.ChengS.ZhouC.WangH.et al (2022). Asmdb: A comprehensive database for allele-specific Dna methylation in diverse organisms.Nucleic Acids Res.50D60D71. 10.1093/nar/gkab937

  • 261

    ZhouY.ZengP.LiY. H.ZhangZ.CuiQ. (2016). Sramp: Prediction of mammalian N6-methyladenosine (m6A) sites based on sequence-derived features.Nucleic Acids Res.44:e91. 10.1093/nar/gkw104

  • 262

    ZhuT.RoundtreeI. A.WangP.WangX.WangL.SunC.et al (2014). Crystal structure of the Yth domain of Ythdf2 reveals mechanism for recognition of N6-methyladenosine.Cell Res.2414931496. 10.1038/cr.2014.152

  • 263

    ZhuZ. M.HuoF. C.PeiD. S. (2020). Function and evolution of Rna N6-methyladenosine modification.Int. J. Biol. Sci.1619291940. 10.7150/ijbs.45231

  • 264

    ZhuangJ.LiuD.LinM.QiuW.LiuJ.ChenS. (2021). PseUdeep: Rna Pseudouridine site identification with deep learning algorithm.Front. Genet.12:773882. 10.3389/fgene.2021.773882

  • 265

    ZhuangY. Y.LiuH. J.SongX.JuY.PengH. (2019). A linear regression predictor for identifying N(6)-methyladenosine sites using frequent gapped K-mer pattern.Mol. Ther. Nucleic Acids18673680. 10.1016/j.omtn.2019.10.001

  • 266

    ZouQ.XingP.WeiL.LiuB. (2019). Gene2vec: Gene subsequence embedding for prediction of mammalian N(6)-methyladenosine sites from mrna.RNA25205218. 10.1261/rna.069112.118

Summary

Keywords

RNA methylation, N6-methyladenosine, 5-methylcytosine, 5-hydroxymethylcytosine, N1-methyladenosine, pseudouridine

Citation

Zou J, Liu H, Tan W, Chen Y, Dong J, Bai S, Wu Z and Zeng Y (2022) Dynamic regulation and key roles of ribonucleic acid methylation. Front. Cell. Neurosci. 16:1058083. doi: 10.3389/fncel.2022.1058083

Received

30 September 2022

Accepted

28 November 2022

Published

19 December 2022

Volume

16 - 2022

Edited by

Mario Gomes-Pereira, Institut National de la Santé et de la Recherche Médicale (INSERM), France

Reviewed by

Bowen Song, University of Liverpool, United Kingdom; Kil To Chong, Jeonbuk National University, South Korea

Updates

Copyright

*Correspondence: Yan Zeng,

†These authors share first authorship

This article was submitted to Cellular Neuropathology, a section of the journal Frontiers in Cellular Neuroscience

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics