- 1Department of Neurosurgery, Zhongnan Hospital, Wuhan University, Wuhan, China
- 2Brain Research Center, Zhongnan Hospital, Wuhan University, Wuhan, China
Glioma is the most common malignant tumor in the central nervous system. The clinical treatment strategy is mainly surgery combined with concurrent temozolomide chemotherapy, but patients can develop drug resistance during treatment, which severely limits its therapeutic efficacy. Epigenetic regulation at the RNA level is plastic and adaptable, and it can induce a variety of tumor responses to drugs. The regulators of RNA modification include methyltransferases, demethylases, and methylation binding proteins; these are also considered to play an important role in the development, prognosis, and therapeutic response of gliomas, which provides a basis for finding new targets of epigenetic drugs and resetting the sensitivity of tumor cells to temozolomide. This review discusses the relationship between the development of adaptive drug resistance and RNA modification in glioma and summarizes the progress of several major RNA modification strategies in this field, especially RNA m6A modification, m5C modification, and adenosine-to-inosine editing.
1 Introduction
Gliomas are the most common primary central nervous system tumors and similar histological features to normal glial cells (i.e, astrocytes, oligodendrocytes, and ependymal cells), accounting for approximately 80% of malignant tumors of the central nervous system (Bauchet and Ostrom, 2019). According to the 2021 World Health Organization classification, glioblastoma (GBM) is defined as IDH-wildtype and assigned a grade 4 classification with high malignant potential, a short survival period, and poor reactions to treatment. GBM mostly occurs in the supratentorial area of the brain. The average age of patients with GBM in developed countries is 64 years old, and they show rapid progression and poor prognosis. The median overall survival time is 16–18 months, and the 5-year survival rate after diagnosis is 5.8% (Alexander and Cloughesy, 2017; Wang et al., 2021a; Ostrom et al., 2021).
In the clinical setting, the standard treatment for patients with primary GBM is surgical resection with a maximum safety margin, postoperative radiotherapy combined with concurrent temozolomide (TMZ) chemotherapy, and then 6-week of TMZ maintenance therapy (Tan et al., 2020). Surgical resection is the most important and effective treatment method for patients with gliomas of any grade. Over the past decade, in order to achieve maximum tumor resection, various methods have been used in neurosurgery to determine functional sites, such as intraoperative ultrasound imaging, fluorescent labeling, and nerve electrophysiology functional (Sanai and Berger, 2018). Despite continuous improvement in surgery techniques, gliomas are difficult to completely remove due to their invasive growth and unclear boundaries. A variety of adjuvant therapies are used in clinical treatment. A variety of new treatment options have been reported, but they are still underdeveloped, such as electric field therapy (Stupp et al., 2017), immunotherapy (Lim et al., 2018; Huang et al., 2020; Zhang et al., 2021a), targeted therapy (Yang et al., 2022a), and nanocarrier combination therapy (Zhao et al., 2020a). Chemotherapy remains an important component of postoperative adjuvant therapy, and TMZ is the most widely used drug. In addition to GBM, in the clinical treatment of patients with recurrent or grade II–III glioma, doctors often use TMZ cautiously according to the patient’s clinical manifestations, tumor pathology, and molecular subtypes, and patient drug sensitivity must be considered.
TMZ is a synthetic small-molecule alkylating agent (Newlands et al., 1997), which can be converted into active substances at physiological pH and precisely deliver methylated fragments to DNA. Clinical trials (Hart et al., 2008; Stupp et al., 2009; Wick et al., 2012) have demonstrated a significant increase in survival and time to tumor progression in patients with GBM after intervention with TMZ. However, the objective response rate of TMZ treatment is only 4–8%, and the final outcome of most patients is still tumor recurrence and death. An important reason for the limited therapeutic effect is that patients with GBM develop resistance to TMZ during long-term standardized chemotherapy.
The cytotoxicity and mutagenicity of TMZ mainly cause DNA damage to O6-MeG (methylation of the O6 position of guanine). Regulation of O6-MeG DNA methyltransferase (MGMT) (Hegi et al., 2005; Cancer Genome Atlas Research Network, 2008; Stupp et al., 2009; Jiang et al., 2014) and an absence of DNA mismatch repair (MMR) (Higuchi et al., 2020) are now widely recognized as the most important factors contributing to TMZ resistance in gliomas. A subset of gliomas with TMZ-induced hypermutations has also been identified (Choi et al., 2018; Daniel et al., 2019). However, these discoveries have not contributed a solution to the clinical problem of TMZ resistance.
There is a growing evidence suggesting that gliomas may evolve adaptive drug resistance through TMZ-induced epigenetic changes, of which Wu et al. (Wu et al., 2021) have conducted a detailed summary. Epigenetics aims at heritable changes in gene expression—without alterations in DNA sequence—using DNA methylation, histone modification, RNA modification, chromatin remodeling, chromatin three-dimensional structural composition, and post-transcriptional regulation of non-coding RNAs, which are all mechanisms that coordinate the regulation of chromatin structure and gene expression. Previous study on epigenetic regulation in cases of glioma drug resistance have mostly focused on genes and proteins (Belter et al., 2020). However, epigenetic drugs targeting glioma DNA-modifying proteins have not achieved satisfactory results. With the development of sequencing technology, the role of RNA modification in TMZ resistance has gradually attracted attention. Finding suitable epigenetic drugs that target the post-transcriptional level may be a novel way to break this deadlock.
There are more than 170 types of biochemical modifications of RNA in living organisms (Dong and Cui, 2020a). These RNA modifications determine the fate of mRNA by affecting RNA metabolic processes such as alternative splicing, nuclear export, maintenance of stability, induction of phase separation, and regulating translation in the mRNA life cycle (Zaccara et al., 2019). RNA modifications regulate gene expression at the post-transcriptional level, participate in maintaining a unique and sophisticated gene expression profile in cells, rapidly and reversibly affect cell phenotype, and play important roles in the healthy physiology and the pathological processes of disease (Lee et al., 2020). In recent years, RNA modification has been a focus of cancer biology; the disruption of epitranscriptomic homeostasis in the physiological state of cells may be a factor initiating carcinogenesis. Moreover, the epistatic regulation of RNA is plastic and adaptive, which can induce alterations in the responsiveness of many cancers to drugs, which is one of the important reasons for the formation of acquired drug resistance in gliomas.
The regulatory roles of RNA modifications in the occurrence, development, proliferation, differentiation, invasion, and metastasis of gliomas have been mostly summarized (Dong and Cui, 2020a; Galardi et al., 2020; Dixit et al., 2021). However, the question of how RNA modifications are involved in glioma drug resistance has not been answered, and related studies are not common and remain controversial. In this review, we discuss the relationship between adaptive drug resistance and RNA modifications in gliomas and summarize the roles of several major RNA methylation modifications in the development of TMZ resistance. This work may help unearth potential drug targets and develop new epitranscriptomic drugs, which may provide clues or directions for solving the issue of glioma chemotherapy resistance and prolonging the survival of glioma patients.
2 The formation of acquired drug resistance in glioma
Tumor chemotherapy resistance is generally divided into primary and acquired drug resistance. Acquired resistance refers to the responses of tumors in the long-term that initially respond to treatment. Tumor heterogeneity and clonal evolution are recognized as important causes of acquired drug resistance. Glioma is highly heterogeneous; during the process of continuous proliferation and variation, tumors form many cell subclones with genomic and epigenetic diversity. The interaction between subclones and the tumor microenvironment (TME) affects the development of tumors (McGranahan and Swanton, 2017). In the past, it was believed that drug-resistant subclones existed in glioma tissue before the start of treatment. After treatment initiation, TMZ exerts selective pressure on tumor cells, and induces drug-resistant clones to become the dominant cell population, ultimately leading to drug resistance and tumor recurrence (Lee, 2016a; Garnier et al., 2018; Happold et al., 2018). Recent studies based on next-generation sequencing have emphasized the role of epigenetics, whereby some cells can survive initial drug treatment through rapid epigenetic adaptation. Signaling in the TME is transduced into cells, through DNA methylation, histone modifications, chromatin conformational changes, and RNA modifications, resulting in transient or persistent resistance phenotypes in tumors before permanent resistance emerges (Eyler et al., 2020) (Figure 1). This is a flexible and reversible phase in the therapeutic response and is the key to reversing drug resistance in gliomas.
FIGURE 1. The formation of acquired drug resistance in glioma; after treatment with TMZ, some TMZ-sensitive glioma subpopulations die, and pre-existing TMZ-resistant glioma subpopulations have the potential to become dominant. At the same time, some precious TMZ-sensitive tumor cells interact with the tumor microenvironment, which has changed after TMZ exposure, and some of these cells acquire TMZ-resistance via epigenetic and genetic changes.
Among glioma cell lines, A172, U87, U251, and U373 are more sensitive to TMZ, while LN-18, T98G, and U138 are somewhat TMZ-resistant, with the former having a lower IC50 (concentration required to reduce the number of surviving cells by 50%; an indicator of drug efficacy) than the latter. However, sensitive cell lines can also be induced by TMZ to acquire resistance, which to some extent mimics the presence of primary and acquired resistance to TMZ in patients (Lee, 2016b).
RNA modification can act indirectly by regulating MGMT gene expression or directly targeting TMZ-related cell death pathways, including drug metabolism and transport, drug targeting, cell death signaling, DNA damage repair, the TME, immunity responses, and maintenance and differentiation of cancer stem cells, among others (Li et al., 2020; Song et al., 2020). Almost all cells have biochemical modifications in RNA that affect RNA structure, metabolism, and function and are dynamically regulated by RNA modification regulators (Boccaletto et al., 2018). Unsupervised clustering analysis of the expression profiles of RNA modification regulators in glioma tumor samples revealed specific expression profiles with specific immunosuppressive microenvironments and mesenchymal transformation, and this class of patients is highly resistant to standard therapy and has a significant survival disadvantage (Xu et al., 2021). This indicates that RNA modification regulators affect the sensitivity of glioma cells to TMZ and directly affect the prognosis of patients with glioma.
3 Regulation of RNA modification in TMZ resistance in glioma
RNA modifications mainly include N6-methyladenosine (m6A), 5-methylcytidine (m5C), N1-methyladenosine (m1A), N3-methyladenosine (m3A), N7-methylguanosine (m7G), pseudouridine, and adenosine-to-inosine (A-to-I) RNA editing. Each modification is reversible and caused by specific regulators, mainly including methyltransferases, demethylases, and methylation binding proteins. These regulators and their related mechanisms in the regulation of TMZ tolerance are provided in Table 1 and Figure 2. In this section, we focus on RNA modifications of m6A and m5C that play strong roles in TMZ resistance in gliomas and tumor chemoresistance.
FIGURE 2. RNA modification and related regulators involved in TMZ resistance. (A) RNA m6A modification and its “writers”, “erasers” and “readers”, each color represents one individual pathway related to m6A modification regulators, and these pathways may finally cause different effect to TMZ treatment; (B) m5C modifications in rRNAs or lncRNAs may play a potential role in TMZ chemoresistance but the mechanisms still remain unclear; (C) RNA m1A modification related genes expression pattern has been proved to be associated with TMZ resistance; (D) A to I editing may effect TMZ resistance by participating in GSC self-renewal and stemness, and the important role of NF–kB pathway has also been realized.
3.1 m6A
3.1.1 m6A modification and its regulatory proteins
RNA m6A modification occurs in approximately 25% of mRNAs and is a dynamic and reversible process that is coordinated by the activity of methyltransferase writers (METTL3, METTL14, WTAP, and METTL16), demethylase activity erasers (FTO and ALKBH5), and m6A-modified readers (YTHDF2 -3 and YTHDC1) working in concert (Liu et al., 2018; Yang et al., 2018); multiple RNA-binding proteins are also involved in the regulation of this process. These m6A modifier proteins are frequently up- or downregulated in human cancer tissues and are often associated with poor patient outcomes (Lan et al., 2021).
The RNA methyltransferase m6A, which provides the methyl group from S-Adenosyl methionine (SAM) and transfers it to the substrate adenosine, has more than 200 members, but only a few have been identified by specific structural features and functions (Oerum et al., 2021). In mammals, especially in humans, the heterodimers formed by the RNA methyltransferases METTL3 and METTL14 are the predominant m6A RNA methyltransferases, of which WTAP and KIAA1429 are important components (Lence et al., 2016). The METTL3–METTL14 complex functions primarily in mRNA to induce adenosine methylation in a specific conserved domain (DRACH; D = A, G or U; R = A or G; H = A, C, U). This conserved domain is present in many mRNAs, but not all mRNAs are recognized by this conserved domain; only specific mRNAs can be catalyzed by METTL3 to undergo m6A methylation. The mechanisms by which this site and transcript specificity exist are unclear. The majority of methyltransferases and demethylases are fixed in the nucleus, and the process of methylation usually occurs in the nucleus in nascent mRNAs that are bound to chromatin. Transcription factors, histone modifiers, RNA binding proteins, RNA polymerases, and many other molecules can recruit m6A RNA methyltransferase and methylate newly transcribed mRNA (Zaccara et al., 2019). Other molecules are involved in RNA methylation in eukaryotes; for example, METTL16 can bind to many types of RNAs, but only some have been shown to be substrates for methylation modifications (Brown et al., 2016). METTL5:TRMT112 (van Tran et al., 2019) and ZCCHC4 (Ma et al., 2019) are important enzymes regulating rRNA methylation.
m6A-RNA pull-down experiments were the first to identify the YTH structural domain-containing RNA-binding proteins YTHDF2 and YTHDF3 (Dominissini et al., 2012). The YTH domain selectively binds the methylated portion of m6A via an internal tryptophan. RNA reading proteins are mainly divided into three main families, YTHDC1, YTHDC2, and YTHDF. YTHDC1 is the main nucleophile that binds rapidly to mRNA after methylation, regulates mRNA splicing by recruiting pre-mRNA splicing factors, and allows access to the target binding regions of mRNAs (Xiao et al., 2016). FTO and ALKBH5 were successively found to have demethylase activity and are both closely related to the growth and proliferation of glioblastoma stem cells (GSCs). FTO promotes the growth and self-renewal of GSCs, which is required for substantial tumor progression (Cui et al., 2017a). Overexpression of ALKBH5 has also been shown to be essential for GSC proliferation and tumorigenesis and can lead to poor prognosis (Zhang et al., 2017).
3.1.2 m6A modification in glioma
Bioinformatics analysis of the China Glioma Genome Atlas Project (CGGA) microarray and RNA-sequencing database revealed that both m6A writers and readers were significantly increased, while erasers such as FTO were relatively decreased (Dong and Cui, 2020b). Xu et al. (Xu et al., 2020a) revealed that 15 m6A regulators are independent prognostic genes using Kaplan–Meier analysis. Patients with low expression of m6A regulators often have favorable prognoses, whereas patients with low expression of FTO, the eraser of m6A modification, have lower survival time compared to those with high FTO expression. These regulatory proteins are involved in the maintenance of the pathological state of gliomas and are highly related to the progression and development of gliomas; however, how m6A modification regulates TMZ resistance in glioma remains unclear.
A PubMed search for the keywords “m6A” and “drug resistance” showed a rapid increase in relevant literature in the last 3 years. This retrieved a total of 99 articles excluding irrelevant and review articles, of which only four articles were about glioma. Gastrointestinal cancers accounted for 30% of articles, lung cancer and reproductive tumors each for 22%, and leukemia 11%, and the other literature was scattered over the fields of skin cancer, bone tumors, lymphomas, thyroid cancer, and urological tumors. With another search of the “References” of these four articles and “similar articles” provided by PubMed, we found another four articles referring to RNA m6A modification in chemotherapy-resistant gliomas. These studies provide some clues to help understand the role of RNA m6A modification in glioma chemotherapy resistance. One study (Li et al., 2021a) conducted a comparative analysis of m6A modifications in TMZ-resistant and TMZ-sensitive clinical GBM samples using m6A miCLIP-seq and revealed that the m6A methyltransferase METTL3 and the EXH2 methyltransferase H3K27me3 form a regulatory loop that exacerbates resistance to TMZ in patients with GBM. The expression of METTL3 was upregulated in drug-resistant cells, resulting in increased m6A modification in the mRNA of the histone modification regulator EZH2. These modifications protected the mRNA of EZH2 from nonsense-mediated mRNA decay), and EZH2 can increase TMZ resistance in gliomas. Several additional studies have also provided other possible mechanisms. For example, a study in 2017 found that overexpression of METTL3 and treatment with the FTO inhibitor MA2 suppressed the growth and self-renewal of GSCs, which are related to the chemo- and radiotherapy-resistance (Cui et al., 2017b). Another study in 2020 provided further clarification about the molecular mechanisms of FTO’s influence on the reaction of glioma cells to TMZ; this may work increasing the stability and translation efficiency of Myc and then participating in regulating the MYC-miR-155/23a Cluster-MXI1 feedback circuit, which is related to the proliferation and tumorigenesis of gliomas (Xiao et al., 2020). JPX (lncRNA proximal inactive specific transcript) mediates PDK1 mRNA demethylation in glioma drug resistance by combining with the m6A demethylase FTO (Li et al., 2021b). The FTO/m6A/MYC/CEBPA signaling pathway may be related to TMZ resistance in glioma cells (Su et al., 2018). ALKBH5-mediated m6A demethylation can increase the expression of NANOG and SOX2 and promote drug resistance in gliomas (Liu et al., 2020; Ding et al., 2021). Furthermore, one study indicated that knockdown of the m6A reader YTHDF1 caused a reduction of GSCs, which may explain the elevated sensitivity of GBM cells to TMZ (Yarmishyn et al., 2020).
The number of relevant studies is too small to form a systematic map of RNA m6A modifications. In general, recent studies have demonstrated that the RNA m6A writer METTL3 can enhance the resistance of glioma cells to TMZ, and the m6A readers such as YTHDF1 and IGF2BP3 may impair this resistance. Different from the exact roles these regulators play, demethylases like FTO and ALKBH5 are both positively and negatively regulated by different cellular signaling pathways (Table 1). Even the same regulatory proteins can have multiple effects on the TMZ-resistant phenotype, and interactions among different regulators of m6A make this more complex. In addition, lncRNA (Liu et al., 2021a), miRNA, histone modification regulatory proteins, and intracellular small-molecule metabolites are also involved in the regulatory process of RNA m6A modification, which further enables RNA m6A modulation of the regulation of drug resistance in glioma.
3.1.3 The effect of m6A modification on tumor chemotherapy resistance is related to tumor type
RNA m6A-modified proteins inhibit or enhance cancer cell sensitivity and resistance to cancer therapy, depending on the individual m6A-modified protein, cancer type, and specific chemotherapeutic agent. In hepatocellular carcinoma, METTL14 induced HNF3γ mRNA m6A modification, resulting in reduced HNF3γ expression and poor response to treatment with the chemotherapy agent sorafenib (Zhou et al., 2020). METTL3 mediates methylation of FOXO3 mRNA and enhances the stability of this RNA by binding to YTHDF1, while overexpression of FOXO3 led to increased sensitivity to sorafenib (Lin et al., 2020). Furthermore inhibition of METTL3-mediated m6A modification of p53 mRNA through the apoptotic pathway enhanced hepatocellular carcinoma sensitivity to sorafenib (Ke et al., 2022). In addition to mRNA m6A modification, circRNA (Xu et al., 2020b) and lncRNA (Kong et al., 2022) m6A modification in hepatocellular carcinoma can also affect drug resistance through multiple pathways. In ovarian cancer, the demethylase ALKBH5, which mediates JAK2 m6A demethylation to activate the JAK2/STAT3 signaling pathway, enhances resistance to cisplatin treatment (Nie et al., 2021a). In non-small cell lung cancer, METTL3-initiated m6A mRNA methylation increased the m6A/A ratio of YAP mRNA, promoted YAP mRNA translation by recruiting YTHDF1/3 and eIF3b, and increased YAP mRNA stability by regulating the MALAT1-miR-1914-3p-YAP axis; increases in YAP expression and activity can induce drug resistance and metastasis in non-small cell lung cancer (Jin et al., 2019). Similar m6A-mediated pathway changes have been observed in bladder cancer (Wei et al., 2021), colorectal cancer (Chen et al., 2021a), and other tumors.
After summarizing different types of tumor chemotherapy resistance mechanisms, we found that the vast majority of studies focused on m6A. The discovery of modified targets focuses on the regulatory mechanisms and downstream signaling pathways of the target. The effect of this target on tumor chemotherapy resistance was demonstrated through cell experiments and animal experiments. These studies suggested that, in drug-resistant tumor cells, RNA m6A modifications can regulate multiple signaling pathways, constituting complex epitranscriptional regulatory homeostasis. However, these studies rarely explored in depth the biological processes of tumor cell adaptation and survival in the drug environment by m6A modification targets and related signaling pathways, and the interactions between different signaling pathways and their ultimate impact on tumor cell phenotype have been rarely explored.
RNA m6A-modifying regulators are potential therapeutic targets for overcoming tumor chemotherapy resistance, and laboratories and drug researchers have begun to develop small molecule activators or inhibitors of m6A modifying proteins. A study (Yankova et al., 2021) on STM2457, an inhibitor of METTL3–METTL14 catalytic activity, screened 250,000 drug-like compounds, and their pharmacological effects were validated using in vitro experiments, animal leukemia models, and additional structure-based drug development approaches (Bedi et al., 2020; Moroz-Omori et al., 2021). UZH1a, a METTL3 inhibitor with good selectivity and cellular activity, has also been confirmed to significantly reduce the m6A/A ratio in mRNA. Small-molecule activators or inhibitors of METTL3, as well as inhibitors of FTO, ALKBH5, and IGF2BP1, showed considerable anticancer effects in some in vitro and animal models of cancer, suggesting that RNA m6A modification has good application prospects in tumor therapy (Lan et al., 2021). However, the role of epigenetic drugs targeting RNA m6A modification regulators in glioma is unknown and awaits further exploration.
3.2 m5C
3.2.1 m5C modification and regulatory proteins
m5C is a well-known methylation modifier shown in several RNA types, and is catalyzed by NSUN protein family members containing NOL1/NOP2/SUN structural domains or DNA methyltransferase homologs (TRDMT1/DNMT2), with S-Adenosyl methionine as the methyl group donor (Agnihotri et al., 2012). m5C is relatively abundant in tRNA and rRNA, less abundant in mRNA, and relatively conserved across species. The knowledge on m5C regulatory enzymes comes mainly from studies on tRNA cytosine methylation (Trixl and Lusser, 2019a). Similar to m6A, m5C has corresponding writers (m5C methyltransferases), erasers (m5C demethylases), and readers (m5C binding proteins) that dynamically regulate intracellular m5C levels.
Ribosomal RNA in all living organisms shows m5C modifications (Motorin et al., 2010), a process mediated mainly by NSUN1 and NSUN5, while NSUN2, NSUN6, TRDMT1, and DNMT2 (Tuorto et al., 2012) have been shown to act on cytoplasmic tRNA. Studies have also found that NSUN2 has a broad range of targets and targeted RNA types (e.g., mRNA and lncRNA) (Squires et al., 2012). One study of (Yang et al., 2017) knockdowns of NSUN family member genes revealed that, among them, only NSUN2 significantly affected m5C levels in mRNA, while the study characterized the specific binding protein ALYREF using biotin-labeled oligonucleotides with or without m5C for RNA affinity chromatography and mass spectrometry. ALYREF is an mRNA export adapter, suggesting that m5C modification may facilitate the nuclear export of mRNA (Liu et al., 2021b). In addition to ALYREF, the readers of m5C include YBX1 (Chen et al., 2019) and FMRP (Yang et al., 2022b). YBX1 is beneficial for preserving the stability of target mRNA, and FMRP is associated with cellular repair and survival in cancer. m5C erasers are members of the TET family and ALKBH1, but whether these proteins act as demethylases remains to be demonstrated. These dynamic elements regulating m5C are summarized in a detailed report by Chen et al. (Chen et al., 2021b).
The methylation m5C loci in mRNA has been controversial, tied to methodological limitations and a process of continuous refinement. Squires et al. (Squires et al., 2012) first combined RNA bisulfite conversion and next-generation sequencing (RNA-BisSeq) to map over 10,000 m5C sites in 2 million cytosines of the human HeLa cell line transcriptome in 2012. In addition, a variety of immunoprecipitation-based m5C detection methods have also emerged, such as Aza-IP (Khoddami and Cairns, 2013; Khoddami and Cairns, 2014) (direct capture of RNA methyltransferase target RNA), miCLIP (Hussain et al., 2013a), m5C-RIP (using m5C monoclonal antibodies). Although these techniques can corroborate the partially deposited m5C sites in RNA, the degree of overlap between the resulting m5C profiles is low (Hussain et al., 2013b). The shortcomings of the detection methods make the sequencing results inevitably noisy; therefore, the detected m5C-related sites are controversial (Trixl and Lusser, 2019b).
A decade of technological development has updated the picture of m5C distribution in RNA, with ever-optimized bisulfite transformation protocols and more stringent analytical treatments being more often applied to m5C sequencing in different species and tissue classes. The prevailing view is now that m5C is not as widely distributed in mRNA as initially thought (Legrand et al., 2017), as shown using computational tools (Legrand et al., 2017; Huang et al., 2019a). Filtering the noise generated during bisulfite sequencing and eliminating artifacts due to incomplete deamidation or RNA secondary structures are the current directions of improvement for detecting mRNA m5C methylation sites using RNA-BisSeq. This method has confirmed that only a few of the previously identified m5C candidate sites are truly methylated marker sites, with median methylation levels of m5C sites of only 15–18%. In mammals, more than 60% of sites are hypomethylated (<20% methylation) and less than 10% are moderately or highly methylated (>40%) (Huang et al., 2019a).
3.2.2 Effects of RNA m5C modification on chemotherapy resistance and glioma
Due to methodological limitations, the molecular role of m5C in mRNA and its biological links to pathophysiological processes remain unclear. Despite recent advances, a clear consensus on the distribution of m5C sites in mRNA has not been established, and reliable annotation of physiological functions is lacking (García-Vílchez et al., 2019). In non-coding RNAs, the functions of m5C modifications are relatively clear. At the molecular level, cytosine methylation modifications within RNAs can increase their structural stability, which may affect the structure and function of ribosomes and tRNA (Motorin and Helm, 2010; Bohnsack et al., 2019; Gkatza et al., 2019). The m5C methyltransferases function in a variety of cellular pathways, mainly involved in stress responses, development, and cell differentiation (García-Vílchez et al., 2019). Among them, NSUN2 is the target of the proto-oncogene c-Myc, which is most closely linked to tumor development (Okamoto et al., 2012) and most likely to be involved in the formation of chemotherapy drug resistance.
NSUN2 is highly expressed in lung cancer (Pan et al., 2021), gastric cancer (Hu et al., 2021), and bladder cancer (Chen et al., 2019). NSUN2 also plays an important role in chemoresistance of tumors, as inhibition of post-transcriptional RNA cytosine methylation modifications can lead to a general decrease in intracellular protein synthesis by activating stress response pathways, and this inhibition makes tumor stem cells sensitive to cytotoxic stress and unable to regenerate and proliferate under chemotherapy drug stress (Blanco et al., 2016). This idea was verified in vitro experiments, where NSUN2 deletion increased the sensitivity of tumor cells to the anticancer drugs 5-fluorouracil and cisplatin in a mouse tumor model (Blanco et al., 2016). In the HeLa cell line, combined knockdown of NSUN2 and METTL1 enhanced the sensitivity of cells to 5-fluorouracil (Okamoto et al., 2014). In esophageal cell carcinoma, NSUN2 upregulated methylated lncRNA, inhibited cisplatin-induced apoptosis, and increased chemoresistance in esophageal squamous cell carcinoma (Li et al., 2018a).
The number of studies on RNA m5C modifications in tumor drug resistance is very small, and its role in glioma drug resistance has not been reported. In glioma, tumor samples showed changes in the expression of several m5C regulators compared to normal tissue, with significant upregulation of ALYREF and NSUN3, downregulation of NSUN2-5 and DNMT1, and no significant changes in TRDMT (Zhang et al., 2022). Analysis of cancer gene mapping data revealed that differential expression of m5C methyltransferase genes in gliomas was significantly associated with clinicopathological features, malignant progression, and prognosis of gliomas. Wang et al. (2020) constructed a risk signature to predict prognosis in patients with glioma. They found that high expression of NOP2, NSUN4, NSUN5, and NSUN7 was linked to poor prognosis, whereas high expression of NSUN6 was linked to better prognosis (Li and Meng, 2021). The ability of NSUN2 to regulate cell migration by promoting mRNA cytosine methylation has been demonstrated in the brain cancer cell line U87 (Xu et al., 2020c; Wang et al., 2020). This evidence describes the role of RNA m5C modifications in the biological behavior of glioma proliferation, differentiation, migration, and malignancy, suggesting the existence of potential m5C epistatic transcriptional modalities affecting temozolomide resistance in glioma.
Similar to NSUN2, deletion of NSUN5 leads to unmethylated rRNA that can drive an overall reduction in protein synthesis, an adaptive translational program developed by cells to survive stress conditions. Epigenetic inactivation of NSUN5 is a hallmark of long-term survival in patients with glioma, but paradoxically, this adaptive regulation may also allow cancer cells to survive in the face of destructive stress conditions (chemotherapeutic agents) to avoid death and ultimately reduce patients’ therapeutic sensitivity to TMZ (Janin et al., 2019). In addition, different m6A/m5C/m1A/m7G-related lncRNA expression patterns have been reported to correlate with glioma prognosis and drug sensitivity (Shao et al., 2022). The characteristic expression patterns of m5C-related lncRNAs were associated with MGMT promoter methylation, epithelial mesenchymal transition (Zhang et al., 2022), tumor immune response, and cell cycle disruption and hypoxia (Zhou et al., 2022), and it may also alter the sensitivity of tumor cells to TMZ-induced cytotoxic effects (Rajesh et al., 2020; Zhang et al., 2022).
The above studies provide clues that RNA m5C modifications affect drug resistance in glioma, but strong direct evidence is lacking. Current research mainly focuses on the perspectives of the m5C methyltransferase NSUN family and m5C-associated lncRNA expression profiles. Maintenance and differentiation of tumor stem cells, adaptive changes in tumor cells under stressful conditions, and alterations in the tumor immune microenvironment are research directions worthy of further exploration.
3.3 Adenosine-to-inosine RNA editing
3.3.1 A-to-I editing and its regulatory protein
A-to-I editing is one of the most abundant RNA modifications, occurring at specific sites in the structure of double-stranded RNA. Adenine nucleotides, catalyzed by the adenosine deaminase family (adenosine deaminases acting on RNA; ADARs) are deaminated to hypoxanthine nucleotides (inosine), also called inosine. During translation, the I on the mRNA is recognized by the ribosome as G (Guanosine), so A-to-I editing alters the amino acid translation codon in the protein-coding sequence, affecting the structure and function of the protein. This pairing also causes a mismatch between RNA sequences and genetic DNA sequences and is used for genomic localization (Heraud-Farlow and Walkley, 2020). Vertebrates have three ADAR genes, ADAR1 (Kim et al., 1994), ADAR2 (Melcher et al., 1996a) and ADAR3 (Melcher et al., 1996b; Chen et al., 2000). ADARs have common functional domains. The dsRNA-binding domain (∼65 amino acids), which has an α-β-β-β-α configuration, makes direct contact with dsRNA (Stefl et al., 1993). The carboxyl-terminus region contains the deaminase domain that forms the catalytic center of an ADAR. In addition to the three vertebrate ADARs, testis nuclear RNA-binding protein (TENR; also known as ADAD1), which is specifically expressed in testes and is required for spermatogenesis (Schumacher et al., 1995), and TENR-like (TENRL; also known as ADAD2), which is expressed in the brain (McKee et al., 2005), have sequence and domain-structure similarity to ADAR but have no deaminase activity, owing to a lack of amino acid residues that are crucial for the catalytic reaction.
3.3.2 A-to-I editing in glioma
Glutamate receptor-mediated signal transmission is associated with the pathogenesis and invasiveness of glioma. The glutamate receptor is re-edited by A-to-I editing at the glutamine (Q)/arginine (R) site of the glutamate receptor subunit B to a degree of nearly 100% in the central nervous system, which is essential for normal receptor function (Maas et al., 2001). It is widely believed that glutamate receptor subunit B mRNA is insufficiently edited in malignant gliomas. However, there is nearly no difference in the expression of ADAR2 between human glioblastoma cell lines (U87, U251, A172) and normal human glial cells (HA 1800) (Li et al., 2015). An alternative splicing variant of ADAR2 with low enzyme activity was finally proven to be the reason for this, and it is dominantly expressed in gliomas and causes the reduced RNA editing of the GluA2 subunit at the glutamine/arginine site in glioma (Wei et al., 2014; Li et al., 2015). ADAR3, a brain-specific, high-expression adenosine deaminase, has been shown to play a similar role to ADAR2 in gliomas (Oakes et al., 2017), and ADAR3 has shown great clinical significance as a potential therapeutic target and useful prognostic factor (Zhang et al., 2018) In the CGGA dataset, downregulation of ADAR3 expression correlated with glioma progression. Zhang et al. 2018) found that ADAR3 was highly expressed in neural molecular subtype glioma and cases with IDH1/2 mutation. Moreover, high expression of ADAR3 predicted a better prognosis in patients with lower-grade glioma, and multivariate analysis suggested ADAR3 expression was an independent prognostic indicator. Jiang et al. (Jiang et al., 2022) found that ADAR1 and GM2A gene that encodes G (M2) ganglioside activator protein were also elevated in GSCs compared with healthy neural stem cells, which impaired GSC self-renewal and stemness. It has further been widely demonstrated that A-to-I editing plays a role in GBM cell proliferation, migration, and invasion (Cesarini et al., 2018; Patil et al., 2020).
3.3.3 The potential regulatory mechanisms of A-to-I RNA editing in temozolomide resistance
Numerous studies suggest that A-to-I RNA editing is involved in tumor chemotherapy resistance via multiple pathways. Studies have described relevant pathways in breast cancer (Nakano et al., 2017) and lung adenocarcinoma (Hua et al., 2022), but unfortunately there is no sufficient description for glioblastoma. A study in 2022 (Raghava Kurup et al., 2022) showed direct experimental evidence that ADAR3 expression increased temozolomide resistance in GBM cells. In that study, Kurup et al. revealed 641 genes that were differentially expressed between control cells and ADAR3-expressing cells, and most of these genes were regulated by NF-κB signaling. Further, GSCs seem to be an essential breakthrough target for A-to-I RNA editing in GBM. ADAR impair GSC self-renewal and stemness through multiple pathways and then may change the reaction to TMZ in gliomas. For example, ADAR1 binds in the 3′-UTR of GM2A for A-to-I editing and links GSC self-renewal to GM2 ganglioside catabolism. TYK2 inhibitors regulate ADAR1 via the upstream JAK/STAT pathway and provide a potential clinical utility (Jiang et al., 2022).
Usually, A-to-I editing acts directly in resistance-related metabolic pathways. In breast cancer, it upregulates the expression of dihydrofolate reductase (DHFR), a target of the chemotherapeutic drug methotrexate, enhancing the proliferation and chemoresistance of breast cancer cells (Nakano et al., 2017). However, the interaction between A-to-I editing and other regulators seems to make a substantial contribution in gliomas. A-to-I editing is significantly enriched in non-coding RNAs and affects the structure and function of non-coding RNAs, which have been shown to be involved in the regulation of drug resistance (Paul et al., 2017; Anthiya et al., 2018; Li et al., 2018b; Li et al., 2019a; Nie et al., 2021b; Tomar et al., 2021). A-to-I editing can also affect the expression of pharmacokinetic- and pharmacodynamic-related genes by mediating changes in miRNA structure and function that determine drug efficacy and toxicity (Nakano and Nakajima, 2018). In general, A-to-I editing and RNA m6A modification always occur independently and do not compete to bind the same sites due to differences in the sequence or structure of the acting RNAs. However, one study (Xiang et al., 2018) demonstrated that m6A-induced changes in secondary RNA structure may regulate the binding of ADAR to the target RNA, thus affecting A-to-I editing, and there may be mutual interference between the two modifications. A-to-I editing also interacts with non-coding RNAs (Liao et al., 2020).
The above studies provide clues for exploring the biological significance of A-to-I RNA editing in the drug resistance phase of glioma. Future studies could be conducted in the following directions: 1) combining sequencing and existing findings to uncover differences in A-to-I editing in chemoresistant and chemosensitive GBM cells, exploring their upstream and downstream pathways, and uncovering drug therapy related genes; 2) finding and integrating the interactions between A-to-I editing and other regulatory molecules to reveal the complete post transcriptional regulatory map of gene expression in glioma drug resistance, providing information for personalized medicine; 3) validating whether ADAR is a potential therapeutic target for glioma drug resistance, which can be used for adjuvant therapy after TMZ treatment by developing molecules that regulate ADAR expression or activity; 4) verifying whether A-to-I editing in glioma affects the TMZ effects by influencing genes related to alkylated drug metabolism, which may provide new options for the treatment of patients with GBM by adjusting the route of administration or co-administration of metabolic enzyme inhibitors.
3.4 Other RNA modifications and small-molecule regulators
Mao et al. (2022) constructed a signature based on m1A regulators with The Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas (TCGA) datasets, which divided gliomas into different subtypes. Four m1A modification-related patterns were found, with different prognoses. Some characteristics of glioma such as stemness, genomic heterogeneity, TME, and immune cell infiltration were also significantly different between the poor and best prognosis clusters. One of these differentially expressed genes, ABCC3, was related to glioma proliferation and TMZ resistance, which may provide new insights into the function and potential therapeutic role of m1A in glioma. Another study showed the potential role of m7G in gliomas. Chen et al. (2022) found that 17 out of 31 regulators of m7G RNA methylation showed significantly higher expression levels in gliomas. Based on these 17 genes, they developed a risk model involving three m7G methylation regulators. Patients were divided into high- and low-risk groups. Overall survival was significantly lower in the high-risk group than in the low-risk group. This finding provides a new insight into the underlying mechanisms of m7G modification in tumorigenesis and the development of gliomas. However, there is no experimental evidence to support these analyses.
The role of other RNA modifications in GBM resistance has not been well described, but there have been some advances in the progression of GBM occurrence and clinical characteristics of patients with GBM. For example, the enzyme pseudouridine synthase 7 is associated with prognosis in patients with GBM (Hu et al., 2020; Wang et al., 2021b) and is involved in regulating the tumor immune microenvironment of GBM, promoting glioma genesis and growth (Cui et al., 2021). Small-molecule inhibitors targeting pseudouridine synthase 7 prolong the life span of mice with GBM (Zhang et al., 2021b). Dysregulation of m1A regulators is closely related to glioma development and progression, and inhibition of TRMT6, a methyltransferase of m1A, inhibits glioma cell proliferation, migration, and invasion (Wang et al., 2021c). These studies provide the basis for further study on the association between RNA modification and TMZ resistance.
We have also noticed the essential effects of non-coding RNA in post-transcriptional gene regulation in TMZ resistance. As a typical example, miRNA, small non-coding RNAs consisting of 20–24 nucleotides, can change the expression profiles of glioma cells through binding with target mRNAs and causing instability and translation inhibition (Zhang et al., 2010; Ahir and Lakka, 2019). miRNA directly targets TMZ resistance-related signaling pathways or interacts with other regulators indirectly. miR-26a reversed the tumor suppressive role of transcription factor AP-2α, which may enhance TMZ resistance by decreasing the expression of MGMT and programmed death-ligand 1 (Huang et al., 2019b). miR-198 enhanced TMZ sensitivity by reducing levels of MGMT and lncRNAs (H19 and HOXD-AS2), and upregulation by TGF-β1 reversed this phenomenon (Nie et al., 2021b). Overexpression of miR-519a enhanced TMZ sensitivity in glioma by inhibiting the STAT3/Bcl-2 signaling pathways (Li et al., 2018b).
Zhao et al. (2020b) found differential expression of mRNAs, lncRNAs, and circRNAs between TMZ-resistant GBMs and primary GBM tissues. They revealed a complex interaction among regulators by constructing a network of lncRNAs, mRNAs, and circRNAs. In fact, lncRNAs often influence regulation by interacting with miRNA-related pathways or, in other words, blocking them. For example, highly expressed in TMZ-resistant glioma cells, lncRNA SNHG12 completely binds with miR-129-5p, causes the upregulation of MAPK1 and E2F7, and then affects TMZ-induced apoptosis via the MAPK/ERK pathway (Lu et al., 2020). lncRNA SNHG15 enhances TMZ-resistance in glioma by decreasing the expression of miR-627-5p (Li et al., 2019b).
4 Discussion
Over the past three decades, the relative survival rate for malignancies has increased from 23% to 36%, but GBM has seen only a slight increase of 3% (Miller et al., 2021). TMZ, the most widely used and studied chemotherapeutic agent for GBM in clinical practice, has shown improvements in patient survival. However, in long-term standardized chemotherapy, patients have developed resistance to TMZ, which severely limits their survival and is an important reason for the maintenance of the 5-year survival rate in GBM. Therefore, expanding our understanding of the molecular mechanisms of drug resistance in glioma, finding new therapeutic targets to improve treatment regimens for patients with GBM, and regaining TMZ sensitivity in patients are the current focus of our work.
Advances in methodology have led to a better understanding of the regulation of gene expression in glioma drug resistance, and the involvement of RNA modification in the formation of a rapid drug resistance phenotype in glioma is an important finding of recent scientific work. Although epigenetic drugs targeting epigenomic regulators and protein modification-regulated enzymes have won approval as a new type of antitumor drugs, with the first clinical applications, there have been no major satisfactory results (Dawson and Kouzarides, 2012; Park and Turcan, 2019; Uddin et al., 2022). The reversibility of RNA modifications provides a new breakthrough to overcome the challenge of drug resistance in glioma. Studies on chemotherapy resistance in a variety of malignancies have shown that modulators of RNA modification are potential targets for drug therapy. However, the question of how RNA modifications are involved in the formation of GBM drug resistance has not been clearly answered, and there are few relevant studies. Future studies should search more for the biological significance of RNA modifications in TMZ tolerance and uncover relevant drug therapeutic targets, which may provide a theoretical basis for the development of epigenetic drugs targeting RNA modifications.
The interactions between RNA modification regulators and interactions of RNA modification with non-coding RNA regulation, epigenomics, proteomics, and metabolomics, also deserve attention. These different levels or dimensions of regulatory factors constitute a complex and stable network of gene expression regulation and are often responsible for the suboptimal clinical outcomes of new drugs. Mapping the complete epigenetic map of gene expression in glioma chemoresistance may provide clues for developing rational drug treatment regimens. Moreover, applying the detection of RNA modifying regulatory factors to the molecular pathological diagnosis of clinical patients and developing individualized treatment plans to achieve precision medicine is also a possible future development direction.
Author contributions
All authors contributed to the design and writing of the manuscript. XL and JC proposed the ideas and drafted the outlines. YY and SL performed the literature search and completed the manuscript. WW helped revise the manuscript and provided support.
Funding
This work was supported by the National Natural Science Foundation of China (82171326).
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
Agnihotri, S., Gajadhar, A. S., Ternamian, C., Gorlia, T., Diefes, K. L., Mischel, P. S., et al. (2012). Alkylpurine-DNA-N-glycosylase confers resistance to temozolomide in xenograft models of glioblastoma multiforme and is associated with poor survival in patients. J. Clin. Investig. 122 (1), 253–266. doi:10.1172/JCI59334
Ahir, B. K., and Lakka, S. S. (2019). Elucidating the microRNA-203 specific biological processes in glioblastoma cells from comprehensive RNA-sequencing transcriptome profiling. Cell. Signal. 53, 22–38. doi:10.1016/j.cellsig.2018.09.014
Alexander, B. M., and Cloughesy, T. F. (2017). Adult glioblastoma. J. Clin. Oncol. 35 (21), 2402–2409. doi:10.1200/JCO.2017.73.0119
Anthiya, S., Griveau, A., Loussouarn, C., Baril, P., Garnett, M., Issartel, J-P., et al. (2018). MicroRNA-based drugs for brain tumors. Trends Cancer 4 (3), 222–238. doi:10.1016/j.trecan.2017.12.008
Bauchet, L., and Ostrom, Q. T. (2019). Epidemiology and molecular epidemiology. Neurosurg. Clin. N. Am. 30 (1), 1–16. doi:10.1016/j.nec.2018.08.010
Bedi, R. K., Huang, D., Eberle, S. A., Wiedmer, L., Śledź, P., and Caflisch, A. (2020). Small-molecule inhibitors of METTL3, the major human epitranscriptomic writer. ChemMedChem 15 (9), 744–748. doi:10.1002/cmdc.202000011
Belter, A., Barciszewski, J., and Barciszewska, A-M. (2020). Revealing the epigenetic effect of temozolomide on glioblastoma cell lines in therapeutic conditions. PloS one 15 (2), e0229534. doi:10.1371/journal.pone.0229534
Blanco, S., Bandiera, R., Popis, M., Hussain, S., Lombard, P., Aleksic, J., et al. (2016). Stem cell function and stress response are controlled by protein synthesis. Nature 534 (7607), 335–340. doi:10.1038/nature18282
Boccaletto, P., Machnicka, M. A., Purta, E., Piatkowski, P., Baginski, B., Wirecki, T. K., et al. (2018). Modomics: A database of RNA modification pathways. 2017 update. Nucleic Acids Res. 46 (D1), D303-D307–D7. doi:10.1093/nar/gkx1030
Bohnsack, K. E., Höbartner, C., and Bohnsack, M. T. (2019). Eukaryotic 5-methylcytosine (m⁵C) RNA methyltransferases: Mechanisms, cellular functions, and links to disease. Genes. 10 (2), E102. doi:10.3390/genes10020102
Brown, J. A., Kinzig, C. G., DeGregorio, S. J., and Steitz, J. A. (2016). Methyltransferase-like protein 16 binds the 3'-terminal triple helix of MALAT1 long noncoding RNA. Proc. Natl. Acad. Sci. U. S. A. 113 (49), 14013–14018. doi:10.1073/pnas.1614759113
Cancer Genome Atlas Research Network (2008). Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature 455 (7216), 1061–1068. doi:10.1038/nature07385
Cesarini, V., Silvestris, D. A., Tassinari, V., Tomaselli, S., Alon, S., Eisenberg, E., et al. (2018). ADAR2/miR-589-3p axis controls glioblastoma cell migration/invasion. Nucleic Acids Res. 46 (4), 2045–2059. doi:10.1093/nar/gkx1257
Chen, C. X., Cho, D. S., Wang, Q., Lai, F., Carter, K. C., and Nishikura, K. (2000). A third member of the RNA-specific adenosine deaminase gene family, ADAR3, contains both single- and double-stranded RNA binding domains. RNA 6 (5), 755–767. doi:10.1017/s1355838200000170
Chen, H., Yao, J., Bao, R., Dong, Y., Zhang, T., Du, Y., et al. (2021). Cross-talk of four types of RNA modification writers defines tumor microenvironment and pharmacogenomic landscape in colorectal cancer. Mol. Cancer 20 (1), 29. doi:10.1186/s12943-021-01322-w
Chen, X., Li, A., Sun, B-F., Yang, Y., Han, Y-N., Yuan, X., et al. (2019). 5-methylcytosine promotes pathogenesis of bladder cancer through stabilizing mRNAs. Nat. Cell. Biol. 21 (8), 978–990. doi:10.1038/s41556-019-0361-y
Chen, Y-S., Yang, W-L., Zhao, Y-L., and Yang, Y-G. (2021). Dynamic transcriptomic m5 C and its regulatory role in RNA processing. Wiley Interdiscip. Rev. RNA 12 (4), e1639. doi:10.1002/wrna.1639
Chen, Z., Zhang, Z., Ding, W., Zhang, J. H., Tan, Z. L., Mei, Y. R., et al. (2022). Expression and potential biomarkers of regulators for M7G RNA modification in gliomas. Front. Neurol. 13, 886246. doi:10.3389/fneur.2022.886246
Choi, S., Yu, Y., Grimmer, M. R., Wahl, M., Chang, S. M., and Costello, J. F. (2018). Temozolomide-associated hypermutation in gliomas. Neuro. Oncol. 20 (10), 1300–1309. doi:10.1093/neuonc/noy016
Cui, Q., Shi, H., Ye, P., Li, L., Qu, Q., Sun, G., et al. (2017). m(6 A RNA methylation regulates the self-renewal and tumorigenesis of glioblastoma stem cells. Cell. Rep. 18 (11), 2622–2634. doi:10.1016/j.celrep.2017.02.059
Cui, Q., Shi, H., Ye, P., Li, L., Qu, Q., Sun, G., et al. (2017). m6A RNA methylation regulates the self-renewal and tumorigenesis of glioblastoma stem cells. Cell. Rep. 18 (11), 2622–2634. doi:10.1016/j.celrep.2017.02.059
Cui, Q., Yin, K., Zhang, X., Ye, P., Chen, X., Chao, J., et al. (2021). Targeting PUS7 suppresses tRNA pseudouridylation and glioblastoma tumorigenesis. Nat. Cancer 2 (9), 932–949. doi:10.1038/s43018-021-00238-0
Daniel, P., Sabri, S., Chaddad, A., Meehan, B., Jean-Claude, B., Rak, J., et al. (2019). Temozolomide induced hypermutation in glioma: Evolutionary mechanisms and therapeutic opportunities. Front. Oncol. 9, 41. doi:10.3389/fonc.2019.00041
Dawson, M. A., and Kouzarides, T. (2012). Cancer epigenetics: From mechanism to therapy. Cell. 150 (1), 12–27. doi:10.1016/j.cell.2012.06.013
Ding, C., Yi, X., Chen, X., Wu, Z., You, H., Chen, X., et al. (2021). Warburg effect-promoted exosomal circ_0072083 releasing up-regulates NANGO expression through multiple pathways and enhances temozolomide resistance in glioma. J. Exp. Clin. Cancer Res. 40 (1), 164. doi:10.1186/s13046-021-01942-6
Dixit, D., Prager, B. C., Gimple, R. C., Poh, H. X., Wang, Y., Wu, Q., et al. (2021). The RNA m6A reader YTHDF2 maintains oncogene expression and is a targetable dependency in glioblastoma stem cells. Cancer Discov. 11 (2), 480–499. doi:10.1158/2159-8290.CD-20-0331
Dominissini, D., Moshitch-Moshkovitz, S., Schwartz, S., Salmon-Divon, M., Ungar, L., Osenberg, S., et al. (2012). Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq. Nature 485 (7397), 201–206. doi:10.1038/nature11112
Dong, Z., and Cui, H. (2020). The emerging roles of RNA modifications in glioblastoma. Cancers (Basel) 12 (3), E736. doi:10.3390/cancers12030736
Dong, Z., and Cui, H. (2020). The emerging roles of RNA modifications in glioblastoma. Cancers (Basel) 12 (3), E736. doi:10.3390/cancers12030736
Eyler, C. E., Matsunaga, H., Hovestadt, V., Vantine, S. J., van Galen, P., and Bernstein, B. E. (2020). Single-cell lineage analysis reveals genetic and epigenetic interplay in glioblastoma drug resistance. Genome Biol. 21 (1), 174. doi:10.1186/s13059-020-02085-1
Galardi, S., Michienzi, A., and Ciafrè, S. A. (2020). Insights into the regulatory role of m(6)A epitranscriptome in glioblastoma. Int. J. Mol. Sci. 21 (8), E2816. doi:10.3390/ijms21082816
García-Vílchez, R., Sevilla, A., and Blanco, S. (2019). Post-transcriptional regulation by cytosine-5 methylation of RNA. Biochim. Biophys. Acta. Gene Regul. Mech. 1862 (3), 240–252. doi:10.1016/j.bbagrm.2018.12.003
Garnier, D., Meehan, B., Kislinger, T., Daniel, P., Sinha, A., Abdulkarim, B., et al. (2018). Divergent evolution of temozolomide resistance in glioblastoma stem cells is reflected in extracellular vesicles and coupled with radiosensitization. Neuro. Oncol. 20 (2), 236–248. doi:10.1093/neuonc/nox142
Gkatza, N. A., Castro, C., Harvey, R. F., Heiß, M., Popis, M. C., Blanco, S., et al. (2019). Cytosine-5 RNA methylation links protein synthesis to cell metabolism. PLoS Biol. 17 (6), e3000297. doi:10.1371/journal.pbio.3000297
Happold, C., Stojcheva, N., Silginer, M., Weiss, T., Roth, P., Reifenberger, G., et al. (2018). Transcriptional control of O(6) -methylguanine DNA methyltransferase expression and temozolomide resistance in glioblastoma. J. Neurochem. 144 (6), 780–790. doi:10.1111/jnc.14326
Hart, M. G., Grant, R., Garside, R., Rogers, G., Somerville, M., and Stein, K. (2008). Temozolomide for high grade glioma. Cochrane Database Syst. Rev. (4), CD007415. doi:10.1002/14651858.CD007415.pub2
Hegi, M. E., Diserens, A-C., Gorlia, T., Hamou, M-F., de Tribolet, N., Weller, M., et al. (2005). MGMT gene silencing and benefit from temozolomide in glioblastoma. N. Engl. J. Med. 352 (10), 997–1003. doi:10.1056/NEJMoa043331
Heraud-Farlow, J. E., and Walkley, C. R. (2020). What do editors do? Understanding the physiological functions of A-to-I RNA editing by adenosine deaminase acting on RNAs. Open Biol. 10 (7), 200085. doi:10.1098/rsob.200085
Higuchi, F., Nagashima, H., Ning, J., Koerner, M. V. A., Wakimoto, H., and Cahill, D. P. (2020). Restoration of temozolomide sensitivity by PARP inhibitors in mismatch repair deficient glioblastoma is independent of base excision repair. Clin. Cancer Res. 26 (7), 1690–1699. doi:10.1158/1078-0432.CCR-19-2000
Hu, G., Wang, R., Wei, B., Wang, L., Yang, Q., Kong, D., et al. (2020). Prognostic markers identification in glioma by gene expression profile Analysis. J. Comput. Biol. 27 (1), 81–90. doi:10.1089/cmb.2019.0217
Hu, Y., Chen, C., Tong, X., Chen, S., Hu, X., Pan, B., et al. (2021). NSUN2 modified by SUMO-2/3 promotes gastric cancer progression and regulates mRNA m5C methylation. Cell. Death Dis. 12 (9), 842. doi:10.1038/s41419-021-04127-3
Hua, H., Zeng, J., Xing, H., He, Y., Han, L., Zhang, N., et al. (2022). The RNA editing enzyme ADAR modulated by the rs1127317 genetic variant diminishes EGFR-TKIs efficiency in advanced lung adenocarcinoma. Life Sci. 296, 120408. doi:10.1016/j.lfs.2022.120408
Huang, B., Li, X., Li, Y., Zhang, J., Zong, Z., and Zhang, H. (2020). Current immunotherapies for glioblastoma multiforme. Front. Immunol. 11, 603911. doi:10.3389/fimmu.2020.603911
Huang, T., Chen, W., Liu, J., Gu, N., and Zhang, R. (2019). Genome-wide identification of mRNA 5-methylcytosine in mammals. Nat. Struct. Mol. Biol. 26 (5), 380–388. doi:10.1038/s41594-019-0218-x
Huang, W., Zhong, Z., Luo, C., Xiao, Y., Li, L., Zhang, X., et al. (2019). The miR-26a/AP-2α/Nanog signaling axis mediates stem cell self-renewal and temozolomide resistance in glioma. Theranostics 9 (19), 5497–5516. doi:10.7150/thno.33800
Hussain, S., Aleksic, J., Blanco, S., Dietmann, S., and Frye, M. (2013). Characterizing 5-methylcytosine in the mammalian epitranscriptome. Genome Biol. 14 (11), 215. doi:10.1186/gb4143
Hussain, S., Sajini, A. A., Blanco, S., Dietmann, S., Lombard, P., Sugimoto, Y., et al. (2013). NSun2-mediated cytosine-5 methylation of vault noncoding RNA determines its processing into regulatory small RNAs. Cell. Rep. 4 (2), 255–261. doi:10.1016/j.celrep.2013.06.029
Janin, M., Ortiz-Barahona, V., de Moura, M. C., Martínez-Cardús, A., Llinàs-Arias, P., Soler, M., et al. (2019). Epigenetic loss of RNA-methyltransferase NSUN5 in glioma targets ribosomes to drive a stress adaptive translational program. Acta Neuropathol. 138 (6), 1053–1074. doi:10.1007/s00401-019-02062-4
Jiang, G., Jiang, A-J., Xin, Y., Li, L-T., Cheng, Q., and Zheng, J-N. (2014). Progression of O⁶-methylguanine-DNA methyltransferase and temozolomide resistance in cancer research. Mol. Biol. Rep. 41 (10), 6659–6665. doi:10.1007/s11033-014-3549-z
Jiang, L., Hao, Y., Shao, C., Wu, Q., Prager, B. C., Gimple, R. C., et al. (2022). ADAR1-mediated RNA editing links ganglioside catabolism to glioblastoma stem cell maintenance. J. Clin. Investig. 132 (6), e143397. doi:10.1172/JCI143397
Jin, D., Guo, J., Wu, Y., Du, J., Yang, L., Wang, X., et al. (2019). m6A mRNA methylation initiated by METTL3 directly promotes YAP translation and increases YAP activity by regulating the MALAT1-miR-1914-3p-YAP axis to induce NSCLC drug resistance and metastasis. J. Hematol. Oncol. 12 (1), 135. doi:10.1186/s13045-019-0830-6
Ke, W., Zhang, L., Zhao, X., and Lu, Z. (2022). p53 m6A modulation sensitizes hepatocellular carcinoma to apatinib through apoptosis. Apoptosis. 27 (5-6), 426–440. doi:10.1007/s10495-022-01728-x
Khoddami, V., and Cairns, B. R. (2013). Identification of direct targets and modified bases of RNA cytosine methyltransferases. Nat. Biotechnol. 31 (5), 458–464. doi:10.1038/nbt.2566
Khoddami, V., and Cairns, B. R. (2014). Transcriptome-wide target profiling of RNA cytosine methyltransferases using the mechanism-based enrichment procedure Aza-IP. Nat. Protoc. 9 (2), 337–361. doi:10.1038/nprot.2014.014
Kim, U., Wang, Y., Sanford, T., Zeng, Y., and Nishikura, K. (1994). Molecular cloning of cDNA for double-stranded RNA adenosine deaminase, a candidate enzyme for nuclear RNA editing. Proc. Natl. Acad. Sci. U. S. A. 91 (24), 11457–11461. doi:10.1073/pnas.91.24.11457
Kong, H., Sun, J., Zhang, W., Zhang, H., and Li, H. (2022). Long intergenic non-protein coding RNA 1273 confers sorafenib resistance in hepatocellular carcinoma via regulation of methyltransferase 3. Bioengineered 13 (2), 3108–3121. doi:10.1080/21655979.2022.2025701
Lan, Q., Liu, P. Y., Bell, J. L., Wang, J. Y., Hüttelmaier, S., Zhang, X. D., et al. (2021). The emerging roles of RNA m6A methylation and demethylation as critical regulators of tumorigenesis, drug sensitivity, and resistance. Cancer Res. 81 (13), 3431–3440. doi:10.1158/0008-5472.CAN-20-4107
Lee, S., Rauch, J., and Kolch, W. (2020). Targeting MAPK signaling in cancer: Mechanisms of drug resistance and sensitivity. Int. J. Mol. Sci. 21 (3), E1102. doi:10.3390/ijms21031102
Lee, S. Y. (2016). Temozolomide resistance in glioblastoma multiforme. Genes. Dis. 3 (3), 198–210. doi:10.1016/j.gendis.2016.04.007
Lee, S. Y. (2016). Temozolomide resistance in glioblastoma multiforme. Genes. Dis. 3 (3), 198–210. doi:10.1016/j.gendis.2016.04.007
Legrand, C., Tuorto, F., Hartmann, M., Liebers, R., Jacob, D., Helm, M., et al. (2017). Statistically robust methylation calling for whole-transcriptome bisulfite sequencing reveals distinct methylation patterns for mouse RNAs. Genome Res. 27 (9), 1589–1596. doi:10.1101/gr.210666.116
Lence, T., Akhtar, J., Bayer, M., Schmid, K., Spindler, L., Ho, C. H., et al. (2016). m6A modulates neuronal functions and sex determination in Drosophila. Nature 540 (7632), 242–247. doi:10.1038/nature20568
Li, B., Jiang, J., Assaraf, Y. G., Xiao, H., Chen, Z-S., and Huang, C. (2020). Surmounting cancer drug resistance: New insights from the perspective of N6-methyladenosine RNA modification. Drug resist. updat. 53, 100720. doi:10.1016/j.drup.2020.100720
Li, F., Chen, S., Yu, J., Gao, Z., Sun, Z., Yi, Y., et al. (2021). Interplay of m6 A and histone modifications contributes to temozolomide resistance in glioblastoma. Clin. Transl. Med. 11 (9), e553. doi:10.1002/ctm2.553
Li, H., Chen, L., Li, J. J., Zhou, Q., Huang, A., Liu, W. W., et al. (2018). miR-519a enhances chemosensitivity and promotes autophagy in glioblastoma by targeting STAT3/Bcl2 signaling pathway. J. Hematol. Oncol. 11 (1), 70. doi:10.1186/s13045-018-0618-0
Li, X., and Meng, Y. (2021). Comment on "physical performance is associated with visual acuity in uUniversitystudents: Results of a school-based study". Rev. Assoc. Med. Bras. 25 (3), 1383–1393. doi:10.1590/1806-9282.20210141
Li, X. D., Wang, M. J., Zheng, J. L., Wu, Y. H., Wang, X., and Jiang, X. B. (2021). Long noncoding RNA just proximal to X-inactive specific transcript facilitates aerobic glycolysis and temozolomide chemoresistance by promoting stability of PDK1 mRNA in an m6A-dependent manner in glioblastoma multiforme cells. Cancer Sci. 112 (11), 4543–4552. doi:10.1111/cas.15072
Li, Y., Li, J., Luo, M., Zhou, C., Shi, X., Yang, W., et al. (2018). Novel long noncoding RNA NMR promotes tumor progression via NSUN2 and BPTF in esophageal squamous cell carcinoma. Cancer Lett. 430, 57–66. doi:10.1016/j.canlet.2018.05.013
Li, Z., Tian, Y., Tian, N., Zhao, X., Du, C., Han, L., et al. (2015). Aberrant alternative splicing pattern of ADAR2 downregulates adenosine-to-inosine editing in glioma. Oncol. Rep. 33 (6), 2845–2852. doi:10.3892/or.2015.3907
Li, Z., Zhang, J., Zheng, H., Li, C., Xiong, J., Wang, W., et al. (2019). Modulating lncRNA SNHG15/CDK6/miR-627 circuit by palbociclib, overcomes temozolomide resistance and reduces M2-polarization of glioma associated microglia in glioblastoma multiforme. J. Exp. Clin. Cancer Res. 38 (1), 380. doi:10.1186/s13046-019-1371-0
Li, Z., Zhang, J., Zheng, H., Li, C., Xiong, J., Wang, W., et al. (2019). Modulating lncRNA SNHG15/CDK6/miR-627 circuit by palbociclib, overcomes temozolomide resistance and reduces M2-polarization of glioma associated microglia in glioblastoma multiforme. J. Exp. Clin. Cancer Res. 38 (1), 380. doi:10.1186/s13046-019-1371-0
Liao, Y., Jung, S. H., and Kim, T. (2020). A-to-I RNA editing as a tuner of noncoding RNAs in cancer. Cancer Lett. 494, 88–93. doi:10.1016/j.canlet.2020.08.004
Lim, M., Xia, Y., Bettegowda, C., and Weller, M. (2018). Current state of immunotherapy for glioblastoma. Nat. Rev. Clin. Oncol. 15 (7), 422–442. doi:10.1038/s41571-018-0003-5
Lin, Z., Niu, Y., Wan, A., Chen, D., Liang, H., Chen, X., et al. (2020). RNA m6 A methylation regulates sorafenib resistance in liver cancer through FOXO3-mediated autophagy. EMBO J. 39 (12), e103181. doi:10.15252/embj.2019103181
Liu, B., Zhou, J., Wang, C., Chi, Y., Wei, Q., Fu, Z., et al. (2020). LncRNA SOX2OT promotes temozolomide resistance by elevating SOX2 expression via ALKBH5-mediated epigenetic regulation in glioblastoma. Cell. Death Dis. 11 (5), 384. doi:10.1038/s41419-020-2540-y
Liu, T., Hu, J., Han, B., Tan, S., Jia, W., and Xin, Y. (2021). A positive feedback loop of lncRNA-RMRP/ZNRF3 axis and Wnt/β-catenin signaling regulates the progression and temozolomide resistance in glioma. Cell. Death Dis. 12 (11), 952. doi:10.1038/s41419-021-04245-y
Liu, Y., Zhao, Y., Wu, R., Chen, Y., Chen, W., Liu, Y., et al. (2021). mRNA m5C controls adipogenesis by promoting CDKN1A mRNA export and translation. RNA Biol. 18 (2), 711–721. doi:10.1080/15476286.2021.1980694
Liu, Z-X., Li, L-M., Sun, H-L., and Liu, S-M. (2018). Link between m6A modification and cancers. Front. Bioeng. Biotechnol. 6, 89. doi:10.3389/fbioe.2018.00089
Lu, C., Wei, Y., Wang, X., Zhang, Z., Yin, J., Li, W., et al. (2020). DNA-methylation-mediated activating of lncRNA SNHG12 promotes temozolomide resistance in glioblastoma. Mol. Cancer 19 (1), 28. doi:10.1186/s12943-020-1137-5
Ma, H., Wang, X., Cai, J., Dai, Q., Natchiar, S. K., Lv, R., et al. (2019). N6-Methyladenosine methyltransferase ZCCHC4 mediates ribosomal RNA methylation. Nat. Chem. Biol. 15 (1), 88–94. doi:10.1038/s41589-018-0184-3
Maas, S., Patt, S., Schrey, M., and Rich, A. (2001). Underediting of glutamate receptor GluR-B mRNA in malignant gliomas. Proc. Natl. Acad. Sci. U. S. A. 98 (25), 14687–14692. doi:10.1073/pnas.251531398
Mao, M., Chu, Q., Lou, Y., Lv, P., and Wang, L. J. (2022). RNA N1-methyladenosine regulator-mediated methylation modification patterns and heterogeneous signatures in glioma. Front. Immunol. 13, 948630. doi:10.3389/fimmu.2022.948630
McGranahan, N., and Swanton, C. (2017). Clonal heterogeneity and tumor evolution: Past, present, and the future. Cell. 168 (4), 613–628. doi:10.1016/j.cell.2017.01.018
McKee, A. E., Minet, E., Stern, C., Riahi, S., Stiles, C. D., and Silver, P. A. (2005). A genome-wide in situ hybridization map of RNA-binding proteins reveals anatomically restricted expression in the developing mouse brain. BMC Dev. Biol. 5, 14. doi:10.1186/1471-213X-5-14
Melcher, T., Maas, S., Herb, A., Sprengel, R., Higuchi, M., and Seeburg, P. H. (1996). RED2, a brain-specific member of the RNA-specific adenosine deaminase family. J. Biol. Chem. 271 (50), 31795–31798. doi:10.1074/jbc.271.50.31795
Melcher, T., Maas, S., Herb, A., Sprengel, R., Seeburg, P. H., and Higuchi, M. (1996). A mammalian RNA editing enzyme. Nature 379 (6564), 460–464. doi:10.1038/379460a0
Miller, K. D., Ostrom, Q. T., Kruchko, C., Patil, N., Tihan, T., Cioffi, G., et al. (2021). Brain and other central nervous system tumor statistics. Ca. Cancer J. Clin. 71 (5), 381–406. doi:10.3322/caac.21693
Moroz-Omori, E. V., Huang, D., Kumar Bedi, R., Cheriyamkunnel, S. J., Bochenkova, E., Dolbois, A., et al. (2021). METTL3 inhibitors for epitranscriptomic modulation of cellular processes. ChemMedChem 16 (19), 3035–3043. doi:10.1002/cmdc.202100291
Motorin, Y., and Helm, M. (2010). tRNA stabilization by modified nucleotides. Biochemistry 49 (24), 4934–4944. doi:10.1021/bi100408z
Motorin, Y., Lyko, F., and Helm, M. (2010). 5-methylcytosine in RNA: Detection, enzymatic formation and biological functions. Nucleic Acids Res. 38 (5), 1415–1430. doi:10.1093/nar/gkp1117
Nakano, M., Fukami, T., Gotoh, S., and Nakajima, M. (2017). A-to-I RNA editing up-regulates human dihydrofolate reductase in breast cancer. J. Biol. Chem. 292 (12), 4873–4884. doi:10.1074/jbc.M117.775684
Nakano, M., and Nakajima, M. (2018). Significance of A-to-I RNA editing of transcripts modulating pharmacokinetics and pharmacodynamics. Pharmacol. Ther. 181, 13–21. doi:10.1016/j.pharmthera.2017.07.003
Newlands, E. S., Stevens, M. F., Wedge, S. R., Wheelhouse, R. T., and Brock, C. (1997). Temozolomide: A review of its discovery, chemical properties, pre-clinical development and clinical trials. Cancer Treat. Rev. 23 (1), 35–61. doi:10.1016/s0305-7372(97)90019-0
Nie, E., Jin, X., Miao, F., Yu, T., Zhi, T., Shi, Z., et al. (2021). TGF-β1 modulates temozolomide resistance in glioblastoma via altered microRNA processing and elevated MGMT. Neuro. Oncol. 23 (3), 435–446. doi:10.1093/neuonc/noaa198
Nie, S., Zhang, L., Liu, J., Wan, Y., Jiang, Y., Yang, J., et al. (2021). ALKBH5-HOXA10 loop-mediated JAK2 m6A demethylation and cisplatin resistance in epithelial ovarian cancer. J. Exp. Clin. Cancer Res. 40 (1), 284. doi:10.1186/s13046-021-02088-1
Oakes, E., Anderson, A., Cohen-Gadol, A., and Hundley, H. A. (2017). Adenosine deaminase that acts on RNA 3 (ADAR3) binding to glutamate receptor subunit B pre-mRNA inhibits RNA editing in glioblastoma. J. Biol. Chem. 292 (10), 4326–4335. doi:10.1074/jbc.M117.779868
Oerum, S., Meynier, V., Catala, M., and Tisné, C. (2021). A comprehensive review of m6A/m6Am RNA methyltransferase structures. Nucleic Acids Res. 49 (13), 7239–7255. doi:10.1093/nar/gkab378
Okamoto, M., Fujiwara, M., Hori, M., Okada, K., Yazama, F., Konishi, H., et al. (2014). tRNA modifying enzymes, NSUN2 and METTL1, determine sensitivity to 5-fluorouracil in HeLa cells. PLoS Genet. 10 (9), e1004639. doi:10.1371/journal.pgen.1004639
Okamoto, M., Hirata, S., Sato, S., Koga, S., Fujii, M., Qi, G., et al. (2012). Frequent increased gene copy number and high protein expression of tRNA (cytosine-5-)-methyltransferase (NSUN2) in human cancers. DNA Cell. Biol. 31 (5), 660–671. doi:10.1089/dna.2011.1446
Ostrom, Q. T., Cioffi, G., Waite, K., Kruchko, C., and Barnholtz-Sloan, J. S. (2021). CBTRUS statistical report: Primary brain and other central nervous system tumors diagnosed in the United States in 2014-2018. Neuro. Oncol. 23 (2), iii1–iii105. doi:10.1093/neuonc/noab200
Pan, J., Huang, Z., and Xu, Y. (2021). m5C RNA methylation regulators predict prognosis and regulate the immune microenvironment in lung squamous cell carcinoma. Front. Oncol. 11, 657466. doi:10.3389/fonc.2021.657466
Park, J-W., and Turcan, Ş. (2019). Epigenetic reprogramming for targeting IDH-mutant malignant gliomas. Cancers (Basel) 11 (10), E1616. doi:10.3390/cancers11101616
Patil, V., Pal, J., Mahalingam, K., and Somasundaram, K. (2020). Global RNA editome landscape discovers reduced RNA editing in glioma: Loss of editing of gamma-amino butyric acid receptor alpha subunit 3 (GABRA3) favors glioma migration and invasion. PeerJ 8, e9755. doi:10.7717/peerj.9755
Paul, D., Sinha, A. N., Ray, A., Lal, M., Nayak, S., Sharma, A., et al. (2017). A-to-I editing in human miRNAs is enriched in seed sequence, influenced by sequence contexts and significantly hypoedited in glioblastoma multiforme. Sci. Rep. 7 (1), 2466. doi:10.1038/s41598-017-02397-6
Raghava Kurup, R., Oakes, E. K., Vadlamani, P., Nwosu, O., Danthi, P., and Hundley, H. A. (2022). ADAR3 activates NF-κB signaling and promotes glioblastoma cell resistance to temozolomide. Sci. Rep. 12 (1), 13362. doi:10.1038/s41598-022-17559-4
Rajesh, Y., Biswas, A., Kumar, U., Banerjee, I., Das, S., Maji, S., et al. (2020). Lumefantrine, an antimalarial drug, reverses radiation and temozolomide resistance in glioblastoma. Proc. Natl. Acad. Sci. U. S. A. 117 (22), 12324–12331. doi:10.1073/pnas.1921531117
Sanai, N., and Berger, M. S. (2018). Surgical oncology for gliomas: The state of the art. Nat. Rev. Clin. Oncol. 15 (2), 112–125. doi:10.1038/nrclinonc.2017.171
Schumacher, J. M., Lee, K., Edelhoff, S., and Braun, R. E. (1995). Distribution of Tenr, an RNA-binding protein, in a lattice-like network within the spermatid nucleus in the mouse. Biol. Reprod. 52 (6), 1274–1283. doi:10.1095/biolreprod52.6.1274
Shao, D., Li, Y., Wu, J., Zhang, B., Xie, S., Zheng, X., et al. (2022). An m6A/m5C/m1A/m7G-Related long non-coding RNA signature to predict prognosis and immune features of glioma. Front. Genet. 13, 903117. doi:10.3389/fgene.2022.903117
Song, H., Liu, D., Dong, S., Zeng, L., Wu, Z., Zhao, P., et al. (2020). Epitranscriptomics and epiproteomics in cancer drug resistance: Therapeutic implications. Signal Transduct. Target. Ther. 5 (1), 193. doi:10.1038/s41392-020-00300-w
Squires, J. E., Patel, H. R., Nousch, M., Sibbritt, T., Humphreys, D. T., Parker, B. J., et al. (2012). Widespread occurrence of 5-methylcytosine in human coding and non-coding RNA. Nucleic Acids Res. 40 (11), 5023–5033. doi:10.1093/nar/gks144
Stefl, R., Xu, M., Skrisovska, L., Emeson, R. B., and Allain, F. H. (1993). Structure and specific RNA binding of ADAR2 double-stranded RNA binding motifs. Structure 14 (2), 345–355. doi:10.1016/j.str.2005.11.013
Stupp, R., Hegi, M. E., Mason, W. P., van den Bent, M. J., Taphoorn, M. J. B., Janzer, R. C., et al. (2009). Effects of radiotherapy with concomitant and adjuvant temozolomide versus radiotherapy alone on survival in glioblastoma in a randomised phase III study: 5-year analysis of the EORTC-NCIC trial. Lancet. Oncol. 10 (5), 459–466. doi:10.1016/S1470-2045(09)70025-7
Stupp, R., Taillibert, S., Kanner, A., Read, W., Steinberg, D., Lhermitte, B., et al. (2017). Effect of tumor-treating fields plus maintenance temozolomide vs maintenance temozolomide alone on survival in patients with glioblastoma: A randomized clinical trial. JAMA 318 (23), 2306–2316. doi:10.1001/jama.2017.18718
Su, R., Dong, L., Li, C., Nachtergaele, S., Wunderlich, M., Qing, Y., et al. (2018). R-2HG exhibits anti-tumor activity by targeting FTO/m(6)A/MYC/CEBPA signaling. Cell. 172 (1-2), 90–105. e23. doi:10.1016/j.cell.2017.11.031
Tan, A. C., Ashley, D. M., López, G. Y., Malinzak, M., Friedman, H. S., and Khasraw, M. (2020). Management of glioblastoma: State of the art and future directions. Ca. Cancer J. Clin. 70 (4), 299–312. doi:10.3322/caac.21613
Tomar, M. S., Kumar, A., Srivastava, C., and Shrivastava, A. (2021). Elucidating the mechanisms of Temozolomide resistance in gliomas and the strategies to overcome the resistance. Biochim. Biophys. Acta. Rev. Cancer 1876 (2), 188616. doi:10.1016/j.bbcan.2021.188616
Trixl, L., and Lusser, A. (2019). Getting a hold on cytosine methylation in mRNA. Nat. Struct. Mol. Biol. 26 (5), 339–340. doi:10.1038/s41594-019-0217-y
Trixl, L., and Lusser, A. (2019). The dynamic RNA modification 5-methylcytosine and its emerging role as an epitranscriptomic mark. Wiley Interdiscip. Rev. RNA 10 (1), e1510. doi:10.1002/wrna.1510
Tuorto, F., Liebers, R., Musch, T., Schaefer, M., Hofmann, S., Kellner, S., et al. (2012). RNA cytosine methylation by Dnmt2 and NSun2 promotes tRNA stability and protein synthesis. Nat. Struct. Mol. Biol. 19 (9), 900–905. doi:10.1038/nsmb.2357
Uddin, M. S., Mamun, A. A., Alghamdi, B. S., Tewari, D., Jeandet, P., Sarwar, M. S., et al. (2022). Epigenetics of glioblastoma multiforme: From molecular mechanisms to therapeutic approaches. Semin. Cancer Biol. 83, 100–120. doi:10.1016/j.semcancer.2020.12.015
van Tran, N., Ernst, F. G. M., Hawley, B. R., Zorbas, C., Ulryck, N., Hackert, P., et al. (2019). The human 18S rRNA m6A methyltransferase METTL5 is stabilized by TRMT112. Nucleic Acids Res. 47 (15), 7719–7733. doi:10.1093/nar/gkz619
Wang, B., Niu, L., Wang, Z., and Zhao, Z. (2021). RNA m1A methyltransferase TRMT6 predicts poorer prognosis and promotes malignant behavior in glioma. Front. Mol. Biosci. 8, 692130. doi:10.3389/fmolb.2021.692130
Wang, L-J., Lv, P., Lou, Y., and Ye, J. (2021). Gene expression-based predication of RNA pseudouridine modification in tumor microenvironment and prognosis of glioma patients. Front. Cell. Dev. Biol. 9, 727595. doi:10.3389/fcell.2021.727595
Wang, M., Zhang, C., Wang, X., Yu, H., Zhang, H., Xu, J., et al. (2021). Tumor-treating fields (TTFields)-based cocktail therapy: A novel blueprint for glioblastoma treatment. Am. J. Cancer Res. 11 (4), 1069–1086.
Wang, P., Wu, M., Tu, Z., Tao, C., Hu, Q., Li, K., et al. (2020). Identification of RNA: 5-Methylcytosine methyltransferases-related signature for predicting prognosis in glioma. Front. Oncol. 10, 1119. doi:10.3389/fonc.2020.01119
Wei, J., Li, Z., Du, C., Qi, B., Zhao, X., Wang, L., et al. (2014). Abnormal expression of an ADAR2 alternative splicing variant in gliomas downregulates adenosine-to-inosine RNA editing. Acta Neurochir. 156 (6), 1135–1142. doi:10.1007/s00701-014-2004-1
Wei, W., Sun, J., Zhang, H., Xiao, X., Huang, C., Wang, L., et al. (2021). Circ0008399 interaction with WTAP promotes assembly and activity of the m6A methyltransferase complex and promotes cisplatin resistance in bladder cancer. Cancer Res. 81 (24), 6142–6156. doi:10.1158/0008-5472.CAN-21-1518
Wick, W., Platten, M., Meisner, C., Felsberg, J., Tabatabai, G., Simon, M., et al. (2012). Temozolomide chemotherapy alone versus radiotherapy alone for malignant astrocytoma in the elderly: The NOA-08 randomised, phase 3 trial. Lancet. Oncol. 13 (7), 707–715. doi:10.1016/S1470-2045(12)70164-X
Wu, Q., Berglund, A. E., and Etame, A. B. (2021). The impact of epigenetic modifications on adaptive resistance evolution in glioblastoma. Int. J. Mol. Sci. 22 (15), 8324. doi:10.3390/ijms22158324
Xiang, J-F., Yang, Q., Liu, C-X., Wu, M., Chen, L-L., and Yang, L. (2018). N6-Methyladenosines modulate A-to-I RNA editing. Mol. Cell. 69 (1), 126–135. doi:10.1016/j.molcel.2017.12.006
Xiao, L., Li, X., Mu, Z., Zhou, J., Zhou, P., Xie, C., et al. (2020). FTO inhibition enhances the antitumor effect of temozolomide by targeting MYC-miR-155/23a cluster-MXI1 feedback circuit in glioma. Cancer Res. 80 (18), 3945–3958. doi:10.1158/0008-5472.CAN-20-0132
Xiao, W., Adhikari, S., Dahal, U., Chen, Y-S., Hao, Y-J., Sun, B-F., et al. (2016). Nuclear m(6)A reader YTHDC1 regulates mRNA splicing. Mol. Cell. 61 (4), 507–519. doi:10.1016/j.molcel.2016.01.012
Xu, J., Gao, Z., Liu, K., Fan, Y., Zhang, Z., Xue, H., et al. (2021). The non-N(6)-methyladenosine epitranscriptome patterns and characteristics of tumor microenvironment infiltration and mesenchymal transition in glioblastoma. Front. Immunol. 12, 809808. doi:10.3389/fimmu.2021.809808
Xu, J., Wan, Z., Tang, M., Lin, Z., Jiang, S., Ji, L., et al. (2020). N6-methyladenosine-modified CircRNA-SORE sustains sorafenib resistance in hepatocellular carcinoma by regulating β-catenin signaling. Mol. Cancer 19 (1), 163. doi:10.1186/s12943-020-01281-8
Xu, S., Tang, L., Dai, G., Luo, C., and Liu, Z. (2020). Expression of m6A regulators correlated with immune microenvironment predicts therapeutic efficacy and prognosis in gliomas. Front. Cell. Dev. Biol. 8, 594112. doi:10.3389/fcell.2020.594112
Xu, X., Zhang, Y., Zhang, J., and Zhang, X. (2020). NSun2 promotes cell migration through methylating autotaxin mRNA. J. Biol. Chem. 295 (52), 18134–18147. doi:10.1074/jbc.RA119.012009
Yang, H., Wang, Y., Xiang, Y., Yadav, T., Ouyang, J., Phoon, L., et al. (2022). FMRP promotes transcription-coupled homologous recombination via facilitating TET1-mediated m5C RNA modification demethylation. Proc. Natl. Acad. Sci. U. S. A. 119 (12), e2116251119. doi:10.1073/pnas.2116251119
Yang, K., Wu, Z., Zhang, H., Zhang, N., Wu, W., Wang, Z., et al. (2022). Glioma targeted therapy: Insight into future of molecular approaches. Mol. Cancer 21 (1), 39. doi:10.1186/s12943-022-01513-z
Yang, X., Yang, Y., Sun, B-F., Chen, Y-S., Xu, J-W., Lai, W-Y., et al. (2017). 5-methylcytosine promotes mRNA export - NSUN2 as the methyltransferase and ALYREF as an m5C reader. Cell. Res. 27 (5), 606–625. doi:10.1038/cr.2017.55
Yang, Y., Hsu, P. J., Chen, Y-S., and Yang, Y-G. (2018). Dynamic transcriptomic m6A decoration: Writers, erasers, readers and functions in RNA metabolism. Cell. Res. 28 (6), 616–624. doi:10.1038/s41422-018-0040-8
Yankova, E., Blackaby, W., Albertella, M., Rak, J., De Braekeleer, E., Tsagkogeorga, G., et al. (2021). Small-molecule inhibition of METTL3 as a strategy against myeloid leukaemia. Nature 593 (7860), 597–601. doi:10.1038/s41586-021-03536-w
Yarmishyn, A. A., Yang, Y. P., Lu, K. H., Chen, Y. C., Chien, Y., Chou, S. J., et al. (2020). Musashi-1 promotes cancer stem cell properties of glioblastoma cells via upregulation of YTHDF1. Cancer Cell. Int. 20 (1), 597. doi:10.1186/s12935-020-01696-9
Zaccara, S., Ries, R. J., and Jaffrey, S. R. (2019). Reading, writing and erasing mRNA methylation. Nat. Rev. Mol. Cell. Biol. 20 (10), 608–624. doi:10.1038/s41580-019-0168-5
Zhang, C., Han, L., Zhang, A., Yang, W., Zhou, X., Pu, P., et al. (2010). Global changes of mRNA expression reveals an increased activity of the interferon-induced signal transducer and activator of transcription (STAT) pathway by repression of miR-221/222 in glioblastoma U251 cells. Int. J. Oncol. 36 (6), 1503–1512. doi:10.3892/ijo_00000637
Zhang, D. Y., Ming, G. L., and Song, H. (2021). PUS7: A targetable epitranscriptomic regulator of glioblastoma growth. Trends Pharmacol. Sci. 42 (12), 976–978. doi:10.1016/j.tips.2021.10.002
Zhang, J., Chen, C., Li, A., Jing, W., Sun, P., Huang, X., et al. (2021). Immunostimulant hydrogel for the inhibition of malignant glioma relapse post-resection. Nat. Nanotechnol. 16 (5), 538–548. doi:10.1038/s41565-020-00843-7
Zhang, J., Wang, N., Wu, J., Gao, X., Zhao, H., Liu, Z., et al. (2022). 5-Methylcytosine related LncRNAs reveal immune characteristics, predict prognosis and oncology treatment outcome in lower-grade gliomas. Front. Immunol. 13, 844778. doi:10.3389/fimmu.2022.844778
Zhang, S., Zhao, B. S., Zhou, A., Lin, K., Zheng, S., Lu, Z., et al. (2017). m6A demethylase ALKBH5 maintains tumorigenicity of glioblastoma stem-like cells by sustaining FOXM1 expression and cell proliferation program. Cancer Cell. 31 (4), 591–606. doi:10.1016/j.ccell.2017.02.013
Zhang, Y., Wang, K., Zhao, Z., Sun, S., Zhang, K., Huang, R., et al. (2018). ADAR3 expression is an independent prognostic factor in lower-grade diffuse gliomas and positively correlated with the editing level of GRIA2(Q607R). Cancer Cell. Int. 18, 196. doi:10.1186/s12935-018-0695-8
Zhao, C., Gao, Y., Guo, R., Li, H., and Yang, B. (2020). Microarray expression profiles and bioinformatics analysis of mRNAs, lncRNAs, and circRNAs in the secondary temozolomide-resistant glioblastoma. Investig. New Drugs 38 (5), 1227–1235. doi:10.1007/s10637-019-00884-3
Zhao, M., van Straten, D., Broekman, M. L. D., Préat, V., and Schiffelers, R. M. (2020). Nanocarrier-based drug combination therapy for glioblastoma. Theranostics 10 (3), 1355–1372. doi:10.7150/thno.38147
Zhou, H., Meng, M., Wang, Z., Zhang, H., Yang, L., Li, C., et al. (2022). The role of m5C-related lncRNAs in predicting overall prognosis and regulating the lower grade glioma microenvironment. Front. Oncol. 12, 814742. doi:10.3389/fonc.2022.814742
Keywords: epigenetic regulation, RNA modification, glioblastoma, multi-omics interactions, m6A, m5C, tumor-acquired drug resistance
Citation: Yan Y, Wei W, Long S, Ye S, Yang B, Jiang J, Li X and Chen J (2022) The role of RNA modification in the generation of acquired drug resistance in glioma. Front. Genet. 13:1032286. doi: 10.3389/fgene.2022.1032286
Received: 31 August 2022; Accepted: 25 October 2022;
Published: 11 November 2022.
Edited by:
Tsung-I Hsu, Taipei Medical University, TaiwanReviewed by:
Wen-Bin Yang, Taipei Medical University, TaiwanFei-Ting Hsu, China Medical University, Taiwan
Copyright © 2022 Yan, Wei, Long, Ye, Yang, Jiang, Li and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiang Li, bGkueGlhbmdAd2h1LmVkdS5jbg==; Jincao Chen, Y2hlbmppbmNhbzIwMTJAMTYzLmNvbQ==
†These authors have contributed equally to this work and share the first authorship