Advances in RNA cytosine-5 methylation: detection, regulatory mechanisms, biological functions and links to cancer
Biomarker Research volume 8, Article number: 43 (2020)
As an important posttranscriptional modification of RNA, 5-methylcytosine (m5C) has attracted increasing interest recently, with accumulating evidence suggesting the involvement of RNA m5C modification in multiple cellular processes as well as tumorigenesis. Cooperatively, advances in m5C detection techniques have enabled transcriptome mapping of RNA methylation at single-nucleotide resolution, thus stimulating m5C-based investigations. In this review, we summarize currently available approaches for detecting m5C distribution in RNA as well as the advantages and disadvantages of these techniques. Moreover, we elucidate the regulatory mechanisms of RNA m5C modification by introducing the molecular structure, catalytic substrates, cellular distributions and biological functions of RNA m5C regulators. The functional consequences of m5C modification on mRNAs, tRNAs, rRNAs and other RNA species, including viral RNAs and vault RNAs, are also discussed. Finally, we review the role of RNA m5C modification in cancer pathogenesis and progression, in hopes of providing new insights into cancer treatment.
Epigenetic modifications, mainly including DNA methylation and histone modifications, are chemical alterations in nucleic acids that do not change the DNA sequence but play key roles in heredity, growth, longevity, aging and disease [1, 2]. DNA 5-methylcytosine (m5C) has been found to be the most abundant DNA modification in mammalian cells and is characterized by the addition of a methyl group at the carbon-5 position of the cytosine base . In eukaryotes, the role of m5C in DNA and its oxidized derivatives, including 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC), has been extensively studied (reviewed in ). In recent years, great progress has been made in the research of RNA modifications, which were originally regarded as fine-tuning chemostructural features of non-protein-coding RNAs but are now considered to be dynamically regulated, reversible and widespread posttranscriptional regulators in diverse cellular processes (reviewed in ). Predominant mammalian RNA methylation modifications include N6-methyladenosine (m6A), N1-methyladenosine (m1A), pseudouridine (Ψ) and m5C, but previous studies on RNA methylation have focused on m6A [5, 6]; however, emerging evidence has gradually uncovered the role of RNA m5C in posttranscriptional regulation.
With the advances in RNA m5C detection techniques, including bisulfite sequencing, m5C RNA immunoprecipitation sequencing (m5C-RIP-seq), 5-AZAcytidine-mediated RNA immunoprecipitation sequencing (Aza-IP-seq) and methylation-individual nucleotide resolution crosslinking immunoprecipitation sequencing (miCLIP-seq), over 10,000 potential sites of m5C modification were detected within the whole human transcriptome , and the existence of m5C was found in diverse RNA species from various organisms, not only in ribosomal RNA (rRNA), messenger RNA (mRNA) and transfer RNA (tRNA) but also in viral RNA, vault RNA, small nuclear RNA, small nucleolar RNA, pseudogenes, natural antisense transcripts, enhancer RNA, long noncoding RNA (lncRNA) and microRNA (miRNA) [8,9,10,11,12,13,14,15]. The distribution of m5C differs among different organisms, as it has been common to find m5C methylation substrates in the tRNA and mRNA of eukaryotes and archaea, while no m5C methylation substrates have been detected in bacteria . It has been demonstrated that RNA m5C is regulated by “writers”, “erasers”, and “readers”, referring to RNA m5C methyltransferases (RCMTs), which are enzymes of the ten-eleven translocation (TET) family and RNA binding proteins . The alteration of m5C modification as well as its regulators has been proven to influence RNA stability, gene expression, and protein synthesis and thus has a great impact on various cellular and physical processes [15, 17,18,19,20,21,22]. Additionally, to date, RNA m5C modification and its regulators have been confirmed to play crucial roles in the pathogenesis of bladder cancer , hepatocellular carcinoma (HCC) , glioblastoma multiforme (GBM)  and leukemia , indicating the promising prospect of m5C modification in cancer treatment.
Herein, we introduce the principles of m5C mapping methods as well as their advantages and disadvantages. Next, we summarize the mechanisms and biological functions of m5C modification of RNAs. Finally, by discussing the role of m5C in tumorigenesis and cancer progression, we propose m5C-based approaches for cancer diagnosis, prognosis and clinical treatment. Although some statements in our study may briefly touch on other reviews [8, 16, 25], which discussed the detection methods and biological functions of RNA m5C methylation in depth, we focus on literature published in the past five years and are the first to thoroughly discuss the clinical prospect of RNA m5C in cancer treatment.
Approaches in detecting m5C distribution in RNA
Although m5C was discovered in RNA almost 60 years ago, little progress has been made regarding the distribution and biological roles of RNA m5C due to a lack of reliable and sensitive methods for RNA methylation detection . With the development of next-generation sequencing (NGS) techniques, RNA methylation can be mapped transcriptome-wide at single-nucleotide resolution. Herein, we summarize the principles of currently available approaches for the characterization of RNA m5C patterns, including bisulfite sequencing, m5C-RIP-seq, Aza-IP-seq and miCLIP-seq. Additionally, their advantages and limitations are discussed to provide a more comprehensive view of RNA m5C mapping (Table 1).
Sodium bisulfite has been applied as a robust and reliable method for DNA m5C analysis since the 1990s . In an acidic pH, sodium bisulfite reacts with methylated and unmethylated C, resulting in deamination to uracil sulfonate/5-methyluracil sulfonate, which can be further converted to uracil/thymine at basic pH. It was not originally considered an available technique for RNA methylation detection, as the stability of phosphodiester bonds in RNA might be disrupted due to harsh chemical and temperature treatment, and the methylation status of individual RNA sequences failed to be properly interrogated . In 2009, a versatile bisulfite deamination protocol for RNA m5C detection was developed based on PCR that allowed amplification of cDNA from low-abundance cellular RNAs . In the presence of sodium bisulfite, unmethylated cytosines (Cs) can be chemically deaminated to uracils (Us), which are ultimately replaced by thymines (Ts) during subsequent PCR amplification, while methylated Cs remain unchanged in this process; following this mechanism, bisulfite sequencing is designed to differentiate unmethylated and methylated Cs and thus investigate the prevalence of cytosine methylation in RNA . Coupled with NGS techniques, bisulfite conversion is now regarded as the gold standard for mapping modified cytosine in RNA at single nucleotide resolution and with high specificity. However, the broader application of bisulfite sequencing is hindered by the degradation of low-abundance RNAs during the conversion process . Additionally, it fails to distinguish m5C from other types of cytosine modifications, such as 5hmC . Further efforts should be made to optimize the validity and scope of bisulfite sequencing.
Originally applied to m6A detection in randomly fragmented RNA transcripts , m5C-RIP was later adapted to detect RNA m5C in bacterial, archaeal, yeast and plant transcriptomes . Conducted with a m5C-specific antibody combined with a deep sequencing technique, this transcriptome-wide protocol is available for mapping m5C sites in low-abundance RNAs; it has high specificity and no harsh reaction conditions, as presented in bisulfite sequencing. However, as the antibody and sequencing method cannot detect modifications with single-nucleotide resolution, as the sequence read by antibodies is usually 100–150 nt in length . Given that bisulfite sequencing cannot distinguish m5C from other types of cytosine modifications, such as 5hmC, m5C-RIP-seq can be used to confirm the methylation of target RNA recognized by bisulfite to improve the stability and accuracy of bisulfite sequencing.
Aza-IP-seq and miCLIP-seq
Aza-IP-seq and miCLIP-seq are designed to investigate specific catalytic sites of RCMTs by immunoprecipitation. All characterized mammalian RCMTs form a reversible covalent intermediate with their cytosine substrate—a covalent linkage that is trapped when formed with the cytosine analog 5-AZAcytidine (5-AZA-C) . Replaced by the cytidine analog 5-AZA-C, an RNA m5C site traps its target methyltransferase by forming a covalent enzyme-RNA adduct, which is then immunoprecipitated with enzyme- or tag-targeted antibodies. During reverse transcription and sequencing, the covalent bond between RNA and methyltransferase is broken, and the release of RNA triggers a hydrolytic opening of the 5-AZA-C ring, which is later read as a guanine (G). Not surprisingly, Aza-IP enables identification of RCMT-specific m5C sites at single nucleotide resolution. However, it is also limited for some reasons, including the toxicity of 5-AZAcytidine to cells, the possible absence of RNAs that are not stably replaced by 5-AZA-C, and random incorporation of cytidine analog into DNAs . Remarkably, a cytosine-to-guanine (C → G) transversion occurs specifically at target cytosines, allowing the simultaneous identification of the precise target cytosine within each RNA.
In contrast to Aza-IP, due to the mutation of the conserved cysteine at the catalytic domain (C271A) of the RNA methyltransferase, an irreversible covalent bond is formed between the enzyme and catalytic site during the miCLIP process . Following immunoprecipitation and deep sequencing, m5C sites are mapped at the + 1 site of the sequencing reads, which is distinct from the C-to-G transversion signatures in Aza-IP. Despite its high specificity, miCLIP requires a large amount of time and money.
RNA m5C regulatory mechanisms
RNA 5-cytosine methylation is a dynamic and reversible process modified by 3 major regulators, referred to as writers, erasers and readers. RCMTs, which include NOL1/NOP2/sun (NSUN) subgroups of methyltransferases and DNA methyltransferase homologue DNMT2, are recognized as writers for catalyzing cytosine-5 methylation . Additionally, RNA m5C binding proteins (readers), such as ALYREF and YBX1, exert biological effects by recognizing and binding to m5C sites . The eraser protein demethylases, including enzymes of the TET family, exhibit reversible effects by mediating written RNA degradation. Herein, we summarize the RNA m5C regulators in detail, emphasizing their molecular structures, catalytic substrates, cellular distributions and biological functions (Table 2).
The DNA methyltransferase homologue DNMT2 contains a Rossman-fold catalytic domain and an S-adenosyl methionine (SAM) binding site. DNMT2 and NSUN2 have complementary target specificities . DNMT2 mainly mediates the m5C modification of tRNA. Studies have confirmed that after knocking out DNMT2, there is no detectable change in mRNA levels, but the steady-state level of unmethylated tRNAs is significantly reduced, resulting in a significant decrease in the protein synthesis rate. These findings indicate that the complete loss of cytosine-5-tRNA methylation has a significant impact on the stability of certain tRNAs and establishes the role of tRNA methylation in protein translation .
Similar to DNMT2, the NSUN family of proteins was also found to contain a Rossman-fold catalytic domain and a SAM binding site. Using SAM as a methyl group donor, RCMTs covalently bind with the RNA cytosine pyrimidine ring through the conserved cysteine residue and catalyze 5-cytosine methylation. Structurally, the NSUN family of proteins contains two cysteine residues in motif IV and motif VI within the catalytic domain . The thiol group of the cysteine in motif VI, together with threonine (TC), forms a conserved dipeptide motif to catalyze the covalent binding between the enzyme and target cytosine, thus promoting the methylation of the carbon-5 atom by electrophilic SAM . Continually, the dipeptide motif formed by cysteine in motif IV and proline (PC) catalyzes covalent bond breakage as well as the release of methylated RNA and enzymes.
The NSUN family of proteins comprises 7 members, including NSUN1, NSUN2, NSUN3, NSUN4, NSUN5, NSUN6 and NSUN7, among which NSUN2 was the first-discovered and is the most well-elucidated member. Dynamic cellular locations of NSUN2 were demonstrated to be associated with cell cycle progression of human epidermal cells, and it was observed to be enriched in the nucleolus in G1, move to the position between nucleoli and nucleoplasm during S phase, and localize to cytoplasm in G2 and centrioles during M phase . Moreover, recent evidence has also shown the mapping of m5C in mitochondrial transfer RNAs (mt-tRNAs) at position C48–50 within mammalian mitochondria . Originally identified to methylate cytoplasmic tRNAs at the junction of the variable loop and the T stem spanning positions 34, 40 and 48–50, NSUN2 has gradually been found to have diverse RNA substrates, including rRNAs, tRNAs, mRNAs, mt-tRNAs and viral RNAs [9,10,11, 43,44,45]. Additionally, NSUN2 exerts multiple functions in the regulation of plant root development, mitochondrial oxidative phosphorylation, protein synthesis and cell cycle progression in response to oxidative stress stimuli , and it has also been reported to regulate HIV replication , Epstein-Barr virus degradation , epidermal differentiation  and tumorigenesis [23, 48, 49]. In addition, studies have shown that NSUN2 plays a role in spindle assembly during mitosis and chromosome segregation. A conserved residue of this gene undergoes a missense change in autosomal recessive mental retardation, causing NSUN2 to fail to localize to the cerebellar Purkinje cell nucleus. Rennai showed that the involvement of NSUN2 highlights the role of RNA methyltransferase in human neurocognitive development . NSUN1, NSUN4 and NSUN5 have been confirmed to play crucial roles in rRNA methylation, contributing to ribosome biogenesis and assembly. In the 1990s, NSUN1 was demonstrated to be involved in 60S ribosomal subunit biogenesis and the inhibition of tumor progression. Predominantly localized to nucleoli, Nop2, the yeast homologue of NSUN1, was reported to catalyze 25S rRNA methylation at position C2870 in domain V, thus affecting 60S large ribosome subunit biogenesis . Recent studies have indicated the role of NSUN1 as a multifunctional protein that is crucially associated with RNA modification, tumor aggressiveness , cell cycle progression , chromatin organization  and HIV-1 latency . NSUN4 is also identified as a rRNA-specific RCMT and is transported into mitochondria in an N-terminal 26-amino acid motif-dependent manner . By interacting with the mitochondrial regulatory factor MTERF4, NSUN4 is recruited to the large subunit of the mitochondrial ribosome, thus facilitating mitochondrial ribosome assembly by modifying 12S rRNA methylation at position C911 . Notably, NSUN4 is now considered an essential m5C regulator that is highly involved in HCC progression . Another rRNA-specific RCMT is NSUN5 (also called Rcm1 in yeast), which is localized to nucleoli, and the methylation substrate of NSUN5 is at C2278 within the conserved domain IV of 25S rRNA . The knockout of NSUN5 results in the extension of life span and increased resistance to multiple stresses in yeast, worms and files, indicating the important role of NSUN5 in cell senescence and stress response .
RCMTs exert broad regulation not only of rRNAs but also of tRNAs, which are predominantly modified by NSUN3, NSUN6 and DNMT2. NSUN3 is localized to mitochondria and is required for the methylation of mitochondrially encoded transfer RNA methionine (mt-tRNAMet) at the “wobble base C34” within the anticodon loop, and the methylated cytosine is further oxidized to 5fC for normal mitochondrial translation of the respiratory chain complex and oxidative phosphorylation. The mutation of NSUN3 is tightly linked to mitochondrial disease and skewed embryonic stem cell differentiation towards the meso- and endoderm lineages instead of neuroectoderm [57, 58]. Another tRNA methylation regulator is NSUN6, which partially resides in the Golgi apparatus and pericentriolar matrix in the cytoplasm and mediates specific methylation on tRNACys and tRNAThr at position C72, which affects tRNAs in a late step in their biogenesis. Both tRNA substrates share a 3′-CCA sequence in the structure, which can be recognized by NSUN6 to carry out cytosine methylation . Distinct from other RCMTs, DNMT2, which was renamed tRNA methyltransferase 1 (TRMDT1), was previously classified as a DNA methyltransferase, as it has the typical sequence and structural features of a DNA methyltransferase . However, further investigations indicated that DNMT2 acts more like an RCMT, since the enrichment of DNMT2 was observed in the cytoplasm and since knockdown of DNMT2 in mice failed to change the DNA methylation level . The reported tRNA substrates of DNMT2 are tRNAAsp-GTC, tRNAGly-GCC and tRNAVal-AAC, all of which contain 3′-GCG and 5′-CA sequences around the methylated C38 position within the anticodon loop . A recent study showed that DNMT2 depletion led to a decrease in the m5C methylation rate in mRNA, indicating that DNMT2 also catalyzes the methylation of other RNA types aside from tRNA . Similar to NSUN2, DNMT2 is also a multifunctional enzyme that broadly participates in tumorigenesis , protein synthesis , cell differentiation , malarial parasite pathogenicity  and HIV-1 RNA survival .
Enzymes of the TET family (TET1, UniProtKB Q8NFU7; TET2, UniProtKB Q6N021; TET3, UniProtKB O43151) are Fe (II)- and α-ketoglutarate (α-KG)-dependent dioxygenases and were initially found to catalyze the oxidation of DNA 5-methyl-2′-deoxycytidine (5mdC) to form 5-hydroxymethyl-2′-deoxycytidine (5hmdC), 5-formyl-2′-deoxycytidine (5fdC), and 5-carboxyl-2′-deoxycytidine (5cadC) . Members of the TET family have distinct cellular distributions, and it was reported that TET3 is located in the nucleus and cytoplasm, while TET1 and TET2 are primarily detected in the nucleus . Recently, the TET family of enzymes has also been found to exhibit activity on RNA 5-methylcytidine (5mrC), and oxidized analogs of 5mrC, including 5-hydroxymethylcytidine (5hmrC), 5-formylcytidine (5frC) and 5-carboxycytidine (5carC), were detected in both coding and noncoding RNAs . It was demonstrated that with the preference for dsDNA, demethylases were tolerant to diverse nucleic acid substrates including dsDNA, ssDNA, ssRNA and DNA-RNA hybrids, and questions remain open regarding how to enhance the specificity and selectivity for TET-mediated oxidation and regarding their structures and biological functions.
Aly/REF export factor (ALYREF), the THO subcomplex and the RNA helicase UAP56, are the 3 major components of the transcription-export complex, which plays an essential role in mRNA nuclear export . Mechanistically, the mRNA export adaptor ALYREF is recruited and binds to the 5′ end region of mature RNA in a CBP80-dependent manner as well as to the 3′ end region in a PABPN1-dependent manner, and by directly interacting with the 3′ processing factor CstF64, the overall binding of ALYREF and mRNA is enhanced [9, 19, 20]. Intriguingly, ALYREF protein shows a clear enrichment in mRNA and retroviral RNA m5C sites [9, 20]. In human HeLa cells and multiple mouse tissues, ALYREF directly binds to mRNA m5C sites and facilitates mRNA nuclear-cytoplasmic shuttling, and the RNA-ALYREF binding affinity as well as the export process are mediated by NSUN2 . In addition, an interaction between ALYREF and methylated retroviral RNA transcripts was observed, and ALYREF was confirmed to play an important role in promoting viral RNA export and replication . However, the molecular structure of the binding adduct and the regulatory mechanisms of ALYREF on RNA m5C remain largely unknown; thus, further investigations are urgently needed.
In addition to ALYREF, Y-box binding protein 1 (YBX1) is a newly discovered m5C reader protein that was confirmed to regulate mRNA stability in the cytoplasm [23, 69]. In human bladder cancer, YBX1 recognizes and binds with m5C-modified mRNAs through an indole ring of W65 within its cold-shock domain (CDS) . Then, by recruiting YBX1-interacting partner ELAV-like protein 1 (ELAVL1), YBX1 regulates mRNA metabolism through the stabilization of m5C-modified mRNAs. Similarly, it was also reported in zebrafish that YBX1 contains a crucial residue, Trp45, in its CDS, through which it can recognize and bind to m5C-modified mRNAs via π-π interactions, thus stabilizing mRNA by recruiting poly(A)-binding protein cytoplasmic 1a .
Functional consequences of m5C modification on different RNA species
Posttranscriptional methylation of 5-cytosine was previously mapped only in tRNAs and rRNAs; however, with the development of high-throughput sequencing technology, m5C is now regarded to broadly occur in mRNAs as well as other noncoding RNAs, including viral RNAs, vault RNAs, lncRNAs and microRNAs. Herein, we summarize the functional consequences of m5C modification on coding and noncoding RNAs, elucidating the underlying mechanisms and providing a novel view for further investigations.
The regulatory effects of m5C modifications on mRNA metabolism have recently become a prevalent issue because of the high abundance of m5C sites in mRNAs. The functional role of m5C regulation on mRNA can be classified into 3 aspects, which refer to the influence on mRNA translation, transport and stability.
The effect of m5C modification on mRNA translation is variable due to distinct the methylation profiles of mRNAs (Fig. 1). In Arabidopsis thaliana, m5C-RIP-seq analysis showed that m5C was predominantly enriched in the coding sequence (CDS), and the sequences located immediately after start codons and before stop codons seemed to have peak m5C modifications . Further investigations suggested that a high abundance of m5C in the CDS was negatively associated with mRNA translation efficiency, as the results showed that a high m5C level was correlated with the presence of ribosomal subunits (40S and 60S) and monosomes (80S), which strongly impaired mRNA translation efficiency . Similarly, it was also shown in the HeLa cell transcriptome that candidate m5C site contents in the CDS were negatively correlated with mRNA translation efficiency . However, another report based on TP53-deficient human colon carcinoma cells (HCT116 p53−/−) and HeLa cells exhibited opposite results . NSUN2-induced m5C methylation, together with methyltransferase like protein 3 (METTL3)/METTL14-induced m6A methylation, synergistically mediated p21 mRNA methylation in the 3′-UTR and facilitated p21 mRNA translation efficiency. Furthermore, the increase in p21 expression profoundly affected the outcome of oxidative stress-induced cellular senescence . The contradiction may be related to distinct m5C profiles, as a recent study showed that the negative correlation between m5C modification and mRNA translation efficiency was significantly associated with methylation site enrichment in the CDS, whereas enhanced methylation in the 3′-UTR was associated with a positive correlation between these features .
Cytosine methylation not only enables mRNA export from the nucleus to cytoplasm in human HeLa cells but also triggers mRNA transport to distant plant regions in A. thaliana (Fig. 2). Following bisulfite sequencing, m5C modification, which was found to be enriched in CG-rich regions as well as regions immediately after start codons, was confirmed to play a key role in mRNA nuclear-cytoplasmic shuttling . Mechanistically, m5C methylation of mRNA was mediated by NSUN2, which increased the binding affinity of mRNA and ALYREF and ultimately modulated nuclear-cytoplasmic export of the ALYREF-mRNA binding adduct . Moreover, m5C modification was also reported to be involved in transcript transport to regulate the development and growth of distant parts of plants. By generating mRNA methylation-deficient DNMT2-NSUN2 double mutants, researchers observed the diminished m5C modification and decreased mobility of graft-mobile transcripts, which ultimately inhibited root growth of A. thaliana .
Recent studies have also uncovered the significant role of m5C modification in affecting mRNA stability. It was reported in A. thaliana that m5C methylation modulated by TRM4B (also named NSUN2 in yeast and animals) plays an essential role in preventing the decay of root development-associated transcripts, thus maintaining mRNA stability and promoting plant growth . Additionally, m5C modification is also involved in the early embryogenesis of zebrafish . During the maternal-to-zygotic transition (MZT) of zebrafish, maternal mRNAs with m5C modification exhibit higher stability than unmodified mRNAs. This could be explained by the following mechanism: YBX1, together with its potential partner Pabpc 1a, specifically recognize and interact with m5C-modified mRNAs, thus maintaining maternal mRNA stability during MZT . In contrast, some studies indicate that no correlation or a negative correlation was observed between the m5C level and mRNA stability [15, 72], and further investigations are urgently needed to clarify the molecular mechanisms of methylation-associated mRNA stabilization.
tRNAs were the first identified substrates of RCMTs, with methylation sites mapped at positions C34, 38, 40, 48–50 and 72 within the anticodon loop. NSUN3 mediates cytosine methylation at position C34 of mitochondrial tRNAs, while NSUN6, NSUN2 and DNMT2 catalyze cytoplasmic tRNA methylation at the following positions: NSUN7-C72, DNMT2-C38, NSUN2-C34, 40, 48–49 [8, 22]. Notably, m5C modifications mediated by NSUN2 at position C48–50 in mt-tRNA have also been detected [42, 73]. tRNA m5C modifications mediated by NSUN2 and DNMT2 have been broadly demonstrated to maintain tRNA stability and regulate cell metabolism [15, 21, 22, 58, 72], as summarized in Fig. 3. In yeast, the proportion of TRM4/NSUN2-modified m5C sites in tRNALeu-CCA is increased to facilitate the translation of survival proteins in response to hydrogen peroxide, suggesting that reprogramming of tRNA modifications is required for the cellular stress response . In humans, mice and plants, it was confirmed that m5C methylation by TRM4/NSUN2 prevents tRNA degradation from oxidative stress, and reduced tRNA stability and hypersensitivity to oxidative stress were observed in TRM4/NSUN2 mutants [15, 74]. Similarly, DNMT2-dependent methylation of tRNA protects it from endonucleolytic cleavage and modulates substrate tRNAAsp-GTC and tRNAGly-GCC stability for accurate polypeptide synthesis during hematopoiesis . Additionally, NSUN2 and DNMT2 were reported to cooperatively stabilize the common substrates tRNAAsp-GTC and tRNAGly-GCC, as double-knockout cells exhibited a robust reduction in abundance in comparison with either DNMT2−/− or NSUN2−/− single-knockout cells .
rRNA methylation can affect the structure of ribonucleoproteins and is essential for ribosome biogenesis. It was confirmed that Nop2/NSUN1 and Rcm1/NSUN5 mediates the methylation of 5-cytosine in domains V and IV of 25S rRNA in yeast  and that NSUN4 catalyzes the formation of m5C at position 911 of mammalian 12S rRNA . The functional role of RCMTs in maintaining rRNA stability was also demonstrated to be similar to that in maintaining tRNA and mRNA stability. The loss of m5C2278, together with the loss of ribose methylation at G2288, resulted in a structural change in 25S RNA and the consequent loss of various ribosome proteins from the 60S large ribosomal subunit under low-salt conditions, thus impairing ribosome biogenesis . Moreover, the effect of Rcm1/NSUN5 alone on ribosome conformation was also elucidated. When cells were exposed to hydrogen peroxide, the loss of Rcm1/NSUN5 led to a more relaxed 25S rRNA fold in the sequence near C2278, indicating that Rcm1/NSUN5 plays an essential role in maintaining rRNA stability in the presence of oxidative stress . Metodiev et al.  proved that the m5C methyltransferase NSUN4 methylates cytosine 911 (m5C911) in small subunit (SSU) 12S rRNA and forms a complex with MTERF4 to ensure that only mature large subunits and SSUs are assembled into monomers. The above evidence proves that NSUN4 plays a key role in controlling the last step of ribosomal biogenesis.
Other RNA species
Cytosine-5 methylation sites are also detected in viral RNAs; however, the regulatory role of m5C modification in viral RNAs is quite controversial. It was reported that DNMT2-mediated m5C methylation promotes HIV-1 genome RNA stability and ensures its survival in the host cell . Moreover, NSUN2-dependent methylation of HIV-1 genomic RNA not only stimulates the translation of HIV-1 mRNA to produce proteins required for HIV-1 replication but also affects alternative HIV-1 RNA splicing . In contrast, the adverse effects of m5C modification on viral RNA replication and stability have also been reported recently. NSUN1 catalyzes the methylation of HIV-1 TAR RNA at the 5′-UTR by competing with HIV-1 Tat protein, thus suppressing HIV-1 replication and facilitating its latency . Additionally, by modulating m5C methylation of Epstein-Barr virus-encoded RNA 2 (EBER2), NSUN2 promotes RNase angiogenin-induced cleavage of EBER2 and leads to a decrease in its level . Therefore, further studies are required to determine the reason different results were observed when investigating the role of m5C in viral RNAs.
Vault RNAs were originally derived from large ribonucleoprotein vault particles and participate in various biological processes ; NSUN2-catalyzed m5C was found in VTRNA1.1 and VTRNA1.3 . NSUN2-mediated methylation of cytosine 69 stimulates the processing of VTRNA1.1 into small vault RNAs, thus promoting epidermal differentiation . Cytosine-5 methylation is also enriched in lncRNAs and was demonstrated to modulate lncRNA functions [9, 15, 77]. The m5C methylation sites were mapped within or near the region mediating interaction with chromatin-modifying complexes in lncRNA HOTAIR and XIST and affected XIST function by influencing the binding of XIST to chromatin-associated protein complex PRC2 . These results showed broad m5C profiles in diverse RNA species and highlighted their essential regulatory roles.
The role of RNA m5C modification in human cancers
Since accumulating evidence has suggested the emerging role of RNA m5C modification in cancer pathogenesis and progression, below, we discuss the essential role of m5C modification as well as the underlying mechanisms in human cancers and propose m5C-based clinical approaches for cancer treatment.
The RNA m5C modification landscape provides insights into dissecting the etiology of bladder cancer pathogenesis. It was reported that the overexpression of NSUN2 and YBX1 in cancer tissues promotes tumorigenesis by maintaining mRNA stability and is significantly correlated with poor prognosis in bladder cancer patients . More specifically, the mRNA of the oncogenic gene heparin binding growth factor (HDGF) is methylated by NSUN2, with high enrichment in its 3′-UTR, and is then recognized by the RNA binding protein YBX1 . By binding to m5C methylated sites and recruiting ELAVL1, YBX1 maintains HDGF mRNA stability and therefore exerts an oncogenic role in bladder cancer development through the activation of HDGF . In addition, the enrichment of m5C mRNAs in tumor-related pathways, such as the PI3K-AKT  and ERK-MAPK signaling pathways , has also been observed, indicating that m5C hypermethylation might initiate tumorigenesis through the activation of oncogenic signaling pathways .
Hepatocellular carcinoma (HCC)
As one of the most common malignances worldwide, HCC remains poorly understood in terms of its pathogenesis and development. Recently, the regulatory role of RNA m5C on HCC progression has been uncovered, and evidence has indicated that the elevated m5C regulators NSUN4 and ALYREF are negatively correlated with poor survival outcome in HCC patients and exhibit great value as significant diagnostic and prognostic biomarkers for HCC . Using 6 pairs of HCC tissues and adjacent nontumor tissues, the same study group demonstrated that the number of mRNA m5C peaks in HCC tissues was significantly higher than that in adjacent tissues, with a more widespread distribution . Together with the bioinformatics analysis results, the study group indicated that the alteration of m5C mRNA is correlated with the pathogenesis of HCC. However, the role of m5C regulators in cancer cell malignant phenotypes and the related molecular mechanisms remain to be further investigated.
Glioblastoma multiforme (GBM)
miRNAs are also common m5C substrates, and the methylation of mature miR-181a-5p results in the abolishment of its tumor suppressor effect and is negatively correlated with poor prognosis in GBM. Mechanistically, mediated by the complex comprising DNMT3a and AGO4 (a miRNA-induced silencing complex), miR-181a-5p methylation impedes the formation of the miRNA-181a-5p/mRNA BIM duplex, which was originally defined to induce apoptosis by interacting with antiapoptotic Bcl-2 or Bcl-xl . As a result, reduced apoptosis, as well as enhanced proliferation and invasion, was observed in cancer cells. In addition, as estimated by Kaplan-Meier analysis, the methylation level of miR-181a-5p is associated with a poor survival rate, suggesting the profound potential for attenuating m5C methylation and restoring normal miRNA function for cancer treatment .
RNA m5C and RCMTs have been reported to participate in drug response/resistance not only in solid tumors but also in leukemia by interacting with distinct partners to regulate the formation of different chromatin structures. NSUN1 specifically interacts with BRD4 and then directly binds to RNA-pol-II CTD-S2P to form a unique NSUN1/BRD4/RNA-pol-II CTD-S2P complex in 5-azacitidine (5-AZA)-resistant leukemia cells, thus mediating the formation of a 5-AZA-resistant chromatin structure and inducing 5-AZA resistance in leukemia . In contrast, NSUN3 and DNMT2 were confirmed to exhibit opposite effects on 5-AZA-sensitive leukemia cells . Mechanistically, the RNA binding protein hnRNPK directly interacts with RCMTs (NSUN3 and DNMT2), lineage-determining transcription factors (GATA1 and SPI1/PU.1), and CDK9/P-TEFb and recruits RNA-polymerase-II at nascent RNA to form a distinct complex that results in the formation of a 5-AZA-sensitive chromatin structure . Using clinical bone marrow specimens from 5-AZA-resistant and 5-AZA-sensitive leukemia cases, a clinical trial showed a significant increase in m5C mRNA in 5-AZA-resistant bone marrow compared with the 5-AZA-sensitive bone marrow. The expression levels of hnRNPK, NSUN1 and BRD4 were correlated with the course of leukemia and were involved in 5-AZA resistance and cancer development . These results predict the potential role of m5C application in the treatment of hematologic malignancies.
Searching for efficient biomarkers for the diagnosis and prognosis of cancer and developing suitable targets for cancer treatment have always been the mainstream approaches against cancer. In HCC, it was confirmed that cancer tissues have higher RNA m5C levels than adjacent normal tissues . Similarly, the number of mRNA m5C peaks was higher in drug-resistant cancer tissues than in drug-sensitive cancer tissues . These results suggest that detecting the level of RNA m5C modification may facilitate cancer diagnosis. In addition, there is evidence suggesting the feasibility of using m5C and m5C regulators as prognostic biomarkers, which are critical for treatment planning and patient life expectancy [20, 23, 51]. Furthermore, as the elevated expression of RNA m5C regulators and the enrichment of m5C in RNA were correlated with the growth, proliferation, invasion, and drug resistance of cancer cells [18, 23, 69], targeting RNA m5C is available for developing new treatment strategies. However, although RNA m5C exhibits promising prospects in the diagnosis, prognosis and treatment of cancer, studies on RNA m5C are limited, and the cancer types studied are too few to construct a thorough review.
In summary, existing evidence has revealed that RNA m5C methylation is essential for maintaining the normal physiological functions of cells and organisms, and its abnormal distribution and abundance are closely related to various diseases, including cancers. Therefore, it is proposed that RNA m5C methylation has great potential in the diagnosis, prognosis and treatment of cancer. Additionally, the mechanisms leading to the dysregulation of various RNA m5C regulators and the pathogenic mechanisms of RNA m5C methylation need to be further explored.
Availability of data and materials
m5C RNA immunoprecipitation sequencing
5-azacytidine–mediated RNA immunoprecipitation sequencing
methylation-individual nucleotide resolution crosslinking immunoprecipitation sequencing
respectively referring to RNA methyltransferases
enzymes of ten-eleven translocation
Y-box binding protein 1
Epstein-Barr virus-encoded RNA 2
heparin binding growth factor
Chen Y, Hong T, Wang S, Mo J, Tian T, Zhou X. Epigenetic modification of nucleic acids: from basic studies to medical applications. Chem Soc Rev. 2017;46:2844–72.
Sen P, Shah PP, Nativio R, Berger SL. Epigenetic mechanisms of longevity and aging. Cell. 2016;166:822–39.
Bestor TH. Cloning of a mammalian DNA methyltransferase. Gene. 1988;74:9–12.
Skvortsova K, Stirzaker C, Taberlay P. The DNA methylation landscape in cancer. Essays Biochem. 2019;63:797–811.
Jonkhout N, Tran J, Smith MA, Schonrock N, Mattick JS, Novoa EM. The RNA modification landscape in human disease. RNA. 2017;23:1754–69.
Roundtree IA, Evans ME, Pan T, He C. Dynamic RNA modifications in gene expression regulation. Cell. 2017;169:1187–200.
Bourgeois G, Ney M, Gaspar I, Aigueperse C, Schaefer M, Kellner S, et al. Eukaryotic rRNA modification by yeast 5-methylcytosine-methyltransferases and human proliferation-associated antigen p120. PLoS One. 2015;10:e0133321.
Garcia-Vilchez R, Sevilla A, Blanco S. Post-transcriptional regulation by cytosine-5 methylation of RNA. Biochim Biophys Acta Gene Regul Mech. 1862;2019:240–52.
Yang X, Yang Y, Sun BF, Chen YS, Xu JW, Lai WY, et al. 5-methylcytosine promotes mRNA export - NSUN2 as the methyltransferase and ALYREF as an m (5) C reader. Cell Res. 2017;27:606–25.
Van Haute L, Lee SY, McCann BJ, Powell CA, Bansal D, Vasiliauskaite L, et al. NSUN2 introduces 5-methylcytosines in mammalian mitochondrial tRNAs. Nucleic Acids Res. 2019;47:8720–33.
Winans S, Beemon K. m (5) C goes viral. Cell Host Microbe. 2019;26:154–5.
Sharma S, Yang J, Watzinger P, Kotter P, Entian KD. Yeast Nop2 and Rcm1 methylate C2870 and C2278 of the 25S rRNA, respectively. Nucleic Acids Res. 2013;41:9062–76.
Xue S, Xu H, Sun Z, Shen H, Chen S, Ouyang J, et al. Depletion of TRDMT1 affects 5-methylcytosine modification of mRNA and inhibits HEK293 cell proliferation and migration. Biochem Biophys Res Commun. 2019;520:60–6.
Cheray M, Etcheverry A, Jacques C, Pacaud R, Bougras-Cartron G, Aubry M, et al. Cytosine methylation of mature microRNAs inhibits their functions and is associated with poor prognosis in glioblastoma multiforme. Mol Cancer. 2020;19:36.
David R, Burgess A, Parker B, Li J, Pulsford K, Sibbritt T, et al. Transcriptome-wide mapping of RNA 5-methylcytosine in Arabidopsis mRNAs and noncoding RNAs. Plant Cell. 2017;29:445–60.
Trixl L, Lusser A. The dynamic RNA modification 5-methylcytosine and its emerging role as an epitranscriptomic mark. Wiley Interdiscip Rev RNA. 2019;10:e1510.
King MY, Redman KL. RNA methyltransferases utilize two cysteine residues in the formation of 5-methylcytosine. Biochemistry. 2002;41:11218–25.
Cheng JX, Chen L, Li Y, Cloe A, Yue M, Wei J, et al. RNA cytosine methylation and methyltransferases mediate chromatin organization and 5-azacytidine response and resistance in leukaemia. Nat Commun. 2018;9:1163.
Shi M, Zhang H, Wu X, He Z, Wang L, Yin S, et al. ALYREF mainly binds to the 5′ and the 3′ regions of the mRNA in vivo. Nucleic Acids Res. 2017;45:9640–53.
Eckwahl M, Xu R, Michalkiewicz J, Zhang W, Patel P, Cai Z, et al. 5-methylcytosine RNA modifications promote retrovirus replication in an ALYREF reader protein-dependent manner. J Virol. 2020;94. https://doi.org/10.1128/JVI.00544-20.
Huang T, Chen W, Liu J, Gu N, Zhang R. Genome-wide identification of mRNA 5-methylcytosine in mammals. Nat Struct Mol Biol. 2019;26:380–8.
Muller M, Samel-Pommerencke A, Legrand C, Tuorto F, Lyko F, Ehrenhofer-Murray AE. Division of labour: tRNA methylation by the NSun2 tRNA methyltransferases Trm4a and Trm4b in fission yeast. RNA Biol. 2019;16:249–56.
Chen X, Li A, Sun BF, Yang Y, Han YN, Yuan X, et al. 5-methylcytosine promotes pathogenesis of bladder cancer through stabilizing mRNAs. Nat Cell Biol. 2019;21:978–90.
He Y, Yu X, Li J, Zhang Q, Zheng Q, Guo W. Role of m (5) C-related regulatory genes in the diagnosis and prognosis of hepatocellular carcinoma. Am J Transl Res. 2020;12:912–22.
Popis MC, Blanco S, Frye M. Posttranscriptional methylation of transfer and ribosomal RNA in stress response pathways, cell differentiation, and cancer. Curr Opin Oncol. 2016;28:65–71.
Amos H, Korn M. 5-methyl cytosine in the RNA of Escherichia coli. Biochim Biophys Acta. 1958;29:444–5.
Frommer M, McDonald LE, Millar DS, Collis CM, Watt F, Grigg GW, et al. A genomic sequencing protocol that yields a positive display of 5-methylcytosine residues in individual DNA strands. Proc Natl Acad Sci U S A. 1992;89:1827–31.
Chakraburtty K. Effect of sodium bisulfite modification on the arginine acceptance of E. coli tRNA Arg. Nucleic Acids Res. 1975;2:1793–804.
Schaefer M, Pollex T, Hanna K, Lyko F. RNA cytosine methylation analysis by bisulfite sequencing. Nucleic Acids Res. 2009;37:e12.
Cui X, Liang Z, Shen L, Zhang Q, Bao S, Geng Y, et al. 5-methylcytosine RNA methylation in Arabidopsis thaliana. Mol Plant. 2017;10:1387–99.
Huber SM, van Delft P, Mendil L, Bachman M, Smollett K, Werner F, et al. Formation and abundance of 5-hydroxymethylcytosine in RNA. Chembiochem. 2015;16:752–5.
Dominissini D, Moshitch-Moshkovitz S, Schwartz S, Salmon-Divon M, Ungar L, Osenberg S, et al. Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq. Nature. 2012;485:201–6.
Edelheit S, Schwartz S, Mumbach MR, Wurtzel O, Sorek R. Transcriptome-wide mapping of 5-methylcytidine RNA modifications in bacteria, archaea, and yeast reveals m5C within archaeal mRNAs. PLoS Genet. 2013;9:e1003602.
Khoddami V, Cairns BR. Transcriptome-wide target profiling of RNA cytosine methyltransferases using the mechanism-based enrichment procedure Aza-IP. Nat Protoc. 2014;9:337–61.
Khoddami V, Cairns BR. Identification of direct targets and modified bases of RNA cytosine methyltransferases. Nat Biotechnol. 2013;31:458–64.
Hussain S, Sajini AA, Blanco S, Dietmann S, Lombard P, Sugimoto Y, et al. NSun2-mediated cytosine-5 methylation of vault noncoding RNA determines its processing into regulatory small RNAs. Cell Rep. 2013;4:255–61.
Bohnsack KE, Hobartner C, Bohnsack MT. Eukaryotic 5-methylcytosine (m (5) C) RNA methyltransferases: Mechanisms, cellular functions, and links to disease. Genes (Basel). 2019;10:102.
Jeltsch A, Ehrenhofer-Murray A, Jurkowski TP, Lyko F, Reuter G, Ankri S, et al. Mechanism and biological role of Dnmt2 in nucleic acid methylation. RNA Biol. 2017;14:1108–23.
Tuorto F, Liebers R, Musch T, Schaefer M, Hofmann S, Kellner S, et al. RNA cytosine methylation by Dnmt2 and NSun2 promotes tRNA stability and protein synthesis. Nat Struct Mol Biol. 2012;19:900–5.
Liu Y, Santi DV. m5C RNA and m5C DNA methyl transferases use different cysteine residues as catalysts. Proc Natl Acad Sci U S A. 2000;97:8263–5.
Frye M, Watt FM. The RNA methyltransferase Misu (NSun2) mediates Myc-induced proliferation and is upregulated in tumors. Curr Biol. 2006;16:971–81.
Shinoda S, Kitagawa S, Nakagawa S, Wei FY, Tomizawa K, Araki K, et al. Mammalian NSUN2 introduces 5-methylcytidines into mitochondrial tRNAs. Nucleic Acids Res. 2019;47:8734–45.
Squires JE, Preiss T. Function and detection of 5-methylcytosine in eukaryotic RNA. Epigenomics. 2010;2:709–15.
Courtney DG, Tsai K, Bogerd HP, Kennedy EM, Law BA, Emery A, et al. Epitranscriptomic addition of m (5) C to HIV-1 transcripts regulates viral gene expression. Cell Host Microbe. 2019;26:217–27 e6.
Henry BA, Kanarek JP, Kotter A, Helm M, Lee N. 5-methylcytosine modification of an Epstein-Barr virus noncoding RNA decreases its stability. RNA. 2020;26:1038–48.
Gkatza NA, Castro C, Harvey RF, Heiss M, Popis MC, Blanco S, et al. Cytosine-5 RNA methylation links protein synthesis to cell metabolism. PLoS Biol. 2019;17:e3000297.
Sajini AA, Choudhury NR, Wagner RE, Bornelov S, Selmi T, Spanos C, et al. Loss of 5-methylcytosine alters the biogenesis of vault-derived small RNAs to coordinate epidermal differentiation. Nat Commun. 2019;10:2550.
Sato G, Saijo Y, Uchiyama B, Kumano N, Sugawara S, Fujimura S, et al. Prognostic value of nucleolar protein p120 in patients with resected lung adenocarcinoma. J Clin Oncol. 1999;17:2721–7.
Jhiang SM, Yaneva M, Busch H. Expression of human proliferation-associated nucleolar antigen p120. Cell Growth Differ. 1990;1:319–24.
Khan MA, Rafiq MA, Noor A, Hussain S, Flores JV, Rupp V, et al. Mutation in NSUN2, which encodes an RNA methyltransferase, causes autosomal-recessive intellectual disability. Am J Hum Genet. 2012;90:856–63.
Kosi N, Alic I, Kolacevic M, Vrsaljko N, Jovanov Milosevic N, Sobol M, et al. Nop2 is expressed during proliferation of neural stem cells and in adult mouse and human brain. Brain Res. 2015;1597:65–76.
Hong J, Lee JH, Chung IK. Telomerase activates transcription of cyclin D1 gene through an interaction with NOL1. J Cell Sci. 2016;129:1566–79.
Kong W, Biswas A, Zhou D, Fiches G, Fujinaga K, Santoso N, et al. Nucleolar protein NOP2/NSUN1 suppresses HIV-1 transcription and promotes viral latency by competing with tat for TAR binding and methylation. PLoS Pathog. 2020;16:e1008430.
Camara Y, Asin-Cayuela J, Park CB, Metodiev MD, Shi Y, Ruzzenente B, et al. MTERF4 regulates translation by targeting the methyltransferase NSUN4 to the mammalian mitochondrial ribosome. Cell Metab. 2011;13:527–39.
Gigova A, Duggimpudi S, Pollex T, Schaefer M, Kos M. A cluster of methylations in the domain IV of 25S rRNA is required for ribosome stability. RNA. 2014;20:1632–44.
Schosserer M, Minois N, Angerer TB, Amring M, Dellago H, Harreither E, et al. Methylation of ribosomal RNA by NSUN5 is a conserved mechanism modulating organismal lifespan. Nat Commun. 2015;6:6158.
Van Haute L, Dietmann S, Kremer L, Hussain S, Pearce SF, Powell CA, et al. Deficient methylation and formylation of mt-tRNA (met) wobble cytosine in a patient carrying mutations in NSUN3. Nat Commun. 2016;7:12039.
Trixl L, Amort T, Wille A, Zinni M, Ebner S, Hechenberger C, et al. RNA cytosine methyltransferase Nsun3 regulates embryonic stem cell differentiation by promoting mitochondrial activity. Cell Mol Life Sci. 2018;75:1483–97.
Haag S, Warda AS, Kretschmer J, Gunnigmann MA, Hobartner C, Bohnsack MT. NSUN6 is a human RNA methyltransferase that catalyzes formation of m5C72 in specific tRNAs. RNA. 2015;21:1532–43.
Dong A, Yoder JA, Zhang X, Zhou L, Bestor TH, Cheng X. Structure of human DNMT2, an enigmatic DNA methyltransferase homolog that displays denaturant-resistant binding to DNA. Nucleic Acids Res. 2001;29:439–48.
Goll MG, Kirpekar F, Maggert KA, Yoder JA, Hsieh CL, Zhang X, et al. Methylation of tRNAAsp by the DNA methyltransferase homolog Dnmt2. Science. 2006;311:395–8.
Shanmugam R, Fierer J, Kaiser S, Helm M, Jurkowski TP, Jeltsch A. Cytosine methylation of tRNA-asp by DNMT2 has a role in translation of proteins containing poly-asp sequences. Cell Discov. 2015;1:15010.
Tuorto F, Herbst F, Alerasool N, Bender S, Popp O, Federico G, et al. The tRNA methyltransferase Dnmt2 is required for accurate polypeptide synthesis during haematopoiesis. EMBO J. 2015;34:2350–62.
Govindaraju G, Jabeena CA, Sethumadhavan DV, Rajaram N, Rajavelu A. DNA methyltransferase homologue TRDMT1 in plasmodium falciparum specifically methylates endogenous aspartic acid tRNA. Biochim Biophys Acta Gene Regul Mech. 1860;2017:1047–57.
Dev RR, Ganji R, Singh SP, Mahalingam S, Banerjee S, Khosla S. Cytosine methylation by DNMT2 facilitates stability and survival of HIV-1 RNA in the host cell during infection. Biochem J. 2017;474:2009–26.
Traube FR, Carell T. The chemistries and consequences of DNA and RNA methylation and demethylation. RNA Biol. 2017;14:1099–107.
Huang Y, Wang G, Liang Z, Yang Y, Cui L, Liu CY. Loss of nuclear localization of TET2 in colorectal cancer. Clin Epigenetics. 2016;8:9.
Huang W, Lan MD, Qi CB, Zheng SJ, Wei SZ, Yuan BF, et al. Formation and determination of the oxidation products of 5-methylcytosine in RNA. Chem Sci. 2016;7:5495–502.
Yang Y, Wang L, Han X, Yang WL, Zhang M, Ma HL, et al. RNA 5-methylcytosine facilitates the maternal-to-zygotic transition by preventing maternal mRNA decay. Mol Cell. 2019;75:1188–202 e11.
Li Q, Li X, Tang H, Jiang B, Dou Y, Gorospe M, et al. NSUN2-mediated m5C methylation and METTL3/METTL14-mediated m6A methylation cooperatively enhance p21 translation. J Cell Biochem. 2017;118:2587–98.
Schumann U, Zhang HN, Sibbritt T, Pan A, Horvath A, Gross S, et al. Multiple links between 5-methylcytosine content of mRNA and translation. BMC Biol. 2020;18:40.
Yang L, Perrera V, Saplaoura E, Apelt F, Bahin M, Kramdi A, et al. M (5) C methylation guides systemic transport of messenger RNA over graft junctions in plants. Curr Biol. 2019;29:2465–76 e5.
Chan CT, Pang YL, Deng W, Babu IR, Dyavaiah M, Begley TJ, et al. Reprogramming of tRNA modifications controls the oxidative stress response by codon-biased translation of proteins. Nat Commun. 2012;3:937.
Blanco S, Dietmann S, Flores JV, Hussain S, Kutter C, Humphreys P, et al. Aberrant methylation of tRNAs links cellular stress to neuro-developmental disorders. EMBO J. 2014;33:2020–39.
Metodiev MD, Spahr H, Loguercio Polosa P, Meharg C, Becker C, Altmueller J, et al. NSUN4 is a dual function mitochondrial protein required for both methylation of 12S rRNA and coordination of mitoribosomal assembly. PLoS Genet. 2014;10:e1004110.
Stadler PF, Chen JJ, Hackermuller J, Hoffmann S, Horn F, Khaitovich P, et al. Evolution of vault RNAs. Mol Biol Evol. 2009;26:1975–91.
Amort T, Rieder D, Wille A, Khokhlova-Cubberley D, Riml C, Trixl L, et al. Distinct 5-methylcytosine profiles in poly(a) RNA from mouse embryonic stem cells and brain. Genome Biol. 2017;18:1.
Amort T, Souliere MF, Wille A, Jia XY, Fiegl H, Worle H, et al. Long non-coding RNAs as targets for cytosine methylation. RNA Biol. 2013;10:1003–8.
Tewari D, Patni P, Bishayee A, Sah AN, Bishayee A. Natural products targeting the PI3K-Akt-mTOR signaling pathway in cancer: A novel therapeutic strategy. Semin Cancer Biol. 2019;S1044-579X(19):30405. https://doi.org/10.1016/j.semcancer.2019.12.008.
Guo YJ, Pan WW, Liu SB, Shen ZF, Xu Y, Hu LL. ERK/MAPK signalling pathway and tumorigenesis. Exp Ther Med. 2020;19:1997–2007.
Zhang Q, Zheng Q, Yu X, He Y, Guo W. Overview of distinct 5-methylcytosine profiles of messenger RNA in human hepatocellular carcinoma and paired adjacent non-tumor tissues. J Transl Med. 2020;18:245.
We thanks for the authors of references for providing the data of the final manuscript.
This work is supported by the National Natural Science Foundation of China (81790631); and the National Key Research and Development Program of China (2018YFC2000500).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Xue, C., Zhao, Y. & Li, L. Advances in RNA cytosine-5 methylation: detection, regulatory mechanisms, biological functions and links to cancer. Biomark Res 8, 43 (2020). https://doi.org/10.1186/s40364-020-00225-0
- RNA modification
- Detection techniques
- Biological functions