METTL3 Promotes the Resistance of Glioma to Temozolomide via Increasing MGMT and ANPG in a m6A Dependent Manner

Acquired chemoresistance is a major limiting factor in the clinical treatment of glioblastoma (GBM). However, the mechanism by which GBM acquires therapeutic resistance remains unclear. Here, we aimed to investigate whether METTL3-mediated N6-methyladenosine (m⁶A) modification contributes to the temozolomide (TMZ) resistance in GBM. We demonstrated that METTL3 METTL3-mediated m⁶A modification were significantly elevated in TMZ-resistant GBM cells. Functionally, METTL3 overexpression impaired the TMZ-sensitivity of GBM cells. In contrast, METTL3 silencing or DAA-mediated total methylation inhibition improved the sensitivity of TMZ-resistant GBM cells to TMZ in vitro and in vivo. Furthermore, we found that two critical DNA repair genes (MGMT and APNG) were m⁶A-modified by METTL3, whereas inhibited by METTL3 silencing or DAA-mediated total methylation inhibition, which is crucial for METTL3-improved TMZ resistance in GBM cells. Collectively, METTL3 acts as a critical promoter of TMZ resistance in glioma and extends the current understanding of m⁶A related signaling, thereby providing new insights into the field of glioma treatment.

Introduction

Owing to the introduction of temozolomide (TMZ), an alkylation agent, and the use of radiotherapy in combination with TMZ adjuvant therapy, the median survival of patients with glioblastoma multiforme (GBM) was increased from 12.1 months to 14.6 months (1–5). However, the overall clinical efficacy of this regimen remains disappointing, mainly because of inherent or induced resistance to TMZ treatment (6–11). Mostly, TMZ-resistant cell lines highly expressed O⁶-methylguanine-DNA methyltransferase (MGMT) and alkylpurine–DNA–N-glycosylase (ANPG) (12, 13). TMZ methylated 12 kinds of DNA bases at different sites, of which, O⁶-meG was considered the most toxic lesion (14). MGMT repairs O⁶-meG through a suicidal response, thereby becoming resistant to TMZ. On the other hand, ANPG repairs the cytotoxic lesions N³-methyladenine and N⁷-methylguanine and contributes to TMZ resistance (12). Therefore, the clinical treatment of this deadly tumor urgently requires a more comprehensive understanding of its progression, mechanisms of resistance, and new therapeutic targets.

In the eukaryotic cells, N6-methyladenosine (m⁶A) is the predominant modification of mRNA and long non-coding RNAs (15). m⁶A is a dynamic and reversible RNA modification in mammalian cells that occurs after transcription by the m⁶A methyltransferase complex, which contains the enzyme subunit methyltransferase-like 3 protein (METTL3) and its co-cofactors methyltransferase-like 14 protein (METTL14) and WT1-associated protein (WTAP) (16). With the deepening understanding of RNA methylation, a number of regulatory factors involved in the regulation of mammalian m⁶A have been identified (17).

The m⁶A modification of the methyltransferase-imprinted RNA prioritizes the recognition and delivery of the reader protein and is cleared by RNA demethylase (18). Therefore, three types of regulators dynamically controlling m⁶A are defined as writers, readers, and erasers (19). Under the control of these three regulatory factors, m⁶A methylation epigenetic regulation of a large number of gene expression plays multiple roles in the regulation of biological processes (20). The acquisition of m⁶A reduces the stability of transcription and mediates the attenuation of target mRNA, suggesting that m⁶A modification is a negative regulator of mRNA translation. Instead, m⁶A deficiency increases the abundance and longevity of transcripts, as well as the overall expression of the protein. m⁶A can also change the structure of RNA, promote the binding of protein regulators, affect mRNA maturation, and regulate gene expression.

It has been reported that m⁶A modification plays a variety of regulatory roles in tumor initiation, progression, and radiation resistance (21, 22). In addition, a growing body of evidences suggests that genetic alterations and dysregulation of m⁶A RNA methylation regulators are closely associated with the malignant progression of a variety of cancers (23). In recent years, increasing evidences have shown that METTL3 plays an important role in cancer as an m⁶A methyltransferase, both as an oncogene and as a tumor suppressor gene. In most cases, METTL3 has been reported as an oncogene that promotes the occurrence and progression of a variety of cancers, including hematopoietic malignancies and solid tumors, by depositing m⁶A modifications on key transcripts (24, 25).

However, the clinicopathological effects of METTL3-mediated RNA m⁶A modification and the related mechanisms of TMZ resistance in glioma have not been elucidated. In this study, we demonstrated that METTL3 acts as a critical promoter of TMZ resistance in glioma. Based on these findings, we provide new insights into the METTL3-mediated modification of m⁶A. We also explored the molecular mechanisms underlying TMZ resistance of glioma by identifying downstream target genes and signals. Therefore, our work extends the current understanding of m⁶A-related signaling and provides new insights into the field of glioma research.

Materials and Method

Cell Lines and Cell Culture

Human glioblastoma-derived U87-MG and U251 cell lines were obtained from the American Type Culture Collection (ATCC, Manassas, VA, USA). All cell lines were cultured in DMEM (Gibco, Grand Island, NY, USA) supplemented with 10% FBS (Gibco) and 1% PS (Invitrogen, Carlsbad, CA, USA), and maintained at 37°C and 5% CO₂ in a humidified atmosphere. TMZ-resistant cell lines were generated by exposure of U87-MG and U251 cells with 200 μM TMZ for over 6 months. The derived resistant cell lines were designated as U87-MG-TMZ resistant and U251-TMZ resistant, respectively. The cell survival ratio and half maximal inhibitory concentration (IC₅₀) of TMZ for U87-MG and U251 was evaluated using the CCK-8 assay (Supplementary Figure 1).

Real-Time Quantitative PCR (qRT-PCR)

RNA extraction and real-time fluorescent quantitative PCR (qRT-PCR) were performed as previously described. The relative gene expression of mRNA was calculated by 2^-ΔΔCT method. GAPDH was used as an endogenous control to normalize the data.

Plasmid Transfection

Stable overexpression of METTL3 was achieved by constructing a lentiviral vector (Biospec Technology, Shanghai). In addition, we synthesized shRNA-targeting genes. Transfection of the expression plasmid in glioma cells was performed using Lipofectamine 3000 kit (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions.

Western Blot

The cells were directly lysed in 1× SDS-PAGE loading buffer. Protein bands were detected sequentially with primary and HRP-bound secondary antibodies, visualized using a Chemiluminescence Detection Kit (Servicebio, Wuhan, China), and detected with an imaging system (Bio-Rad, USA). Antibodies against METTL3 (AB195352, 1:2000) were obtained from Abcam. GAPDH (60004, 1:5000) antibodies were purchased from Proteintech.

Total RNA m⁶A Quantification

The total level of m⁶A in the treated glioma cells was determined using the EpiQuik™ m⁶A RNA Methylation Quantitative Kit (Epigentek, USA). Briefly, 200 ng of RNA was added to each well, followed by a mixture of capture and detection antibodies. After several weeks of incubation, the m⁶A content was quantified at 450 nm and calculated according to the standard curve.

Dot Blot

The mRNA samples were dissolved in a 3-fold volume of RNA incubation buffer, denatured at 65°C for 5 min, and loaded onto an Amersham Hybond-N+membrane (GE Healthcare, USA) mounted on a Bio-Dot device (Bio-Rad, USA). After blocking the membrane with 5% skimmed milk, the specific m⁶A antibody (1:1000, Abcam) was incubated overnight at 4°C. Mouse immunoglobulin G (IgG) was incubated with HRP-conjugated immunoglobulin G (IgG) for 1 h, and imaging was performed using an imaging system (Bio-Rad, USA).

Methylated RNA Immunoprecipitation (Me-RIP)

Total RNA or poly(A)+mRNA was isolated using the above methods. The purified mRNAs and magnetic bead-antibody complexes were then added to IP buffers and incubated overnight at 4°C, followed by elution with eluent and purification. MGMT and ANPG in RNA were extracted using RT-qPCR.

Cell Viability Assay

Cell viability was measured after treatment with different concentrations of TMZ (Selleck Chemicals, Houston, TX, USA). After 4 h of normal culture, 10 µL CCK-8 reagent (Dojindo) was added and absorbance at 450 nm was detected using an ultra-multifunctional microplate analyzer (Tecan, Durham, NC, USA). Using GraphPad Prism 9.0 (GraphPad Software, San Diego, CA, USA), the “log (inhibitor) vs normalized slope of response variable” method was used to calculate the 50% inhibition concentration (IC₅₀) of TMZ.

Colony Formation Experiments

Glioma cells were seeded in a 6-well culture plate containing 500 cells per well for 14 d. The colonies were washed with PBS and fixed with 4% paraformaldehyde. Photographs were taken using a microscope (Olympus, Ishikawa, Japan).

Subcutaneous Glioma Xenograft Model

All experiments involving mice were conducted in accordance with the ethical standards of the animal care and use committee of the third hospital affiliated to Soochow University and the NIH guidelines for the care and use of laboratory animals. To establish the xenograft model of glioma in mice, 1×10⁷ human U87-MG-TMZ cells (sh-Con, sh-METTL3, or normal U87-MG-TMZ) were subcutaneously inoculated into the right posterior limb of BALB/c nude mice (6-week-old, female) in 80 μL PBS. Tumor volume was measured with calipers every 5 d. After approximately 30 d, all mice were euthanized, and the tumor masses were removed, weighed, and embedded for further pathological study.

Statistical Analyses

SPSS 21.0 statistical software (IBM Corp. Armonk, NY) was used for statistical analyses, and statistical significance was set at P<0.05. Data are expressed as mean ± standard deviation. Multiple sets of data were evaluated using one-way analysis of variance (ANOVA), and multiple comparisons were performed using Tukey’s post-hoc test. Time-based multiple comparisons were tested by repeated analysis of variance and the Bonferroni post facto test.

Results

METTL3 Mediated m⁶A Is Elevated in the TMZ-Resistant GBM Cells

Previously, elevated METTL3 levels have been associated with malignant characteristics of cancer cells (21), but its role in TMZ resistance in GBM has not been fully understood. Here, upon comparing the METTL3 levels between the TMZ-sensitive cells and the resistant cells, we found that the mRNA level of METTL3 was significantly higher (about 4.78-fold in U87-MG-TMZ and 4.48-fold in U251-TMZ) in the TMZ-resistant group than in the sensitive group (Figure 1A), which was further confirmed by western blot analysis (Figure 1B). We then examined m⁶A levels in the total RNAs from TMZ-sensitive cells and resistant cells using the colorimetric m⁶A quantification strategy, revealing significantly increased m⁶A levels in TMZ-resistant cells (Figure 1C) compared with TMZ-sensitive cells (approximately 4.16-fold in U87-MG-TMZ and 5.92-fold in U251-TMZ), confirmed by dot blot analysis (Figure 1D). These results suggest that METTL3 mediated m⁶A may contribute to TMZ-resistant GBM cells.

FIGURE 1

Figure 1 METTL3-mediated m⁶A is elevated in the TMZ-resistant GBM cells. (A) The mRNA level of METTL3 in TMZ sensitive and resistant U87-MG/U251 cells was analyzed by real-time PCR. (B) The protein level of METTL3 in TMZ sensitive and resistant U87-MG/U251 cells was analyzed by western blot. (C) The colorimetric m⁶A quantification assay was used to examine the total m⁶A levels in the TMZ sensitive cells and the resistant U87-MG/U251 cells. (D) The dot blot was used to confirm the total m⁶A levels in the TMZ sensitive cells and the resistant U87-MG/U251 cells. *P < 0.05, **P < 0.01, and ***P< 0.001 versus normal U87-MG/U251 cells.

METTL3 Contributes to the TMZ Resistance in GBM Cells

To further study the functional role of METTL3 in the regulation of TMZ resistance, we established METTL3-stable overexpression and knockdown U87-MG-TMZ and U251-TMZ cell lines. The efficiency of overexpression and knockdown on the mRNA and protein levels of METTL3 was verified by qRT-PCR (Figure 2A) and western blot (Figure 2B), respectively. Consistently, the m⁶A levels were significantly increased in METTL3 overexpressed U87-MG-TMZ and U251-TMZ cells, whereas decreased in METTL3 knockdown U87-MG-TMZ and U251-TMZ cells (Figure 2C). Compared with parental control, METTL3 knockdown GBM cells had a significantly lower ability to form colonies (Figure 2D), while TMZ-Resistant cells overexpressing METTL3 had no effect (Figure 2D). More importantly, METTL3 level was positively correlated with TMZ sensitivity. When METTL3 was knocked down, the IC₅₀ value decreased from approximately 268.9 μM to 95.6 μM in U87-MG-TMZ cells and 296.0 μM to 110.6 μM in U251-TMZ, whereas the IC₅₀ value remained unchanged in METTL3 overexpressing cells (Figure 2E). These results suggest that METTL3 silencing caused TMZ-resistant cells more sensitive to TMZ.

FIGURE 2

Figure 2 METTL3 contributes to the TMZ resistance in GBM cells. (A) The efficiency of overexpression and knockdown on the mRNA levels of METTL3 were analyzed by qRT-PCR. (B) The efficiency of overexpression and knockdown on the protein levels of METTL3 were analyzed by western blot. (C) The effect of METTL3 overexpression and knockdown on the total m⁶A RNA level was analyzed by the colorimetric m⁶A quantification assay. (D) The effect of METTL3 overexpression and knockdown on the cell proliferation was analyzed by the colony formation assay. (E) The effect of METTL3 overexpression and knockdown on the sensitivity to TMZ was analyzed by CCK-8 assay. *P < 0.05, **P < 0.01, and ***P < 0.001 versus indicated control U87-MG/U251 cells.

METTL3 Contributes to the TMZ Resistance via m⁶A Modification

To further study the functional role of METTL3-mediated m⁶A modification in the regulation of TMZ resistance, we inhibited methylation with a methylation inhibitor, 3-deazaadenosine (DAA, 100 µM). Consistent with our hypothesis, treating U87-MG/U251-TMZ cells with DAA led to a remarkable reduction in total m⁶A level (Figure 3A), which was verified by dot blot (Figure 3B). Moreover, compared with the parental control, DAA-treated GBM cells had a significantly lower ability to form colonies (Figure 3C). Furthermore, the IC₅₀ value decreased from approximately 275.4 μM to 98.6 μM in U87-MG-TMZ cells and 288.2 μM to 108.3 μM in U251-TMZ (Figure 3D). The major repair enzymes, O6-methylguanine–DNA methyltransferase (MGMT) and alkylpurine–DNA–N-glycosylase (APNG), repairs the most cytotoxic lesions generated by TMZ. To analyze the underlying mechanism of METTL3-mediated m⁶A modification in the regulation of TMZ resistance, we screened a series of TMZ-resistant genes (ANPG, CBX5, MGMT, MSH2, MSH6, MLH1, MPG, XRCC3, and XPC), revealing that METTL3 overexpression significantly increased the MGMT and ANPG expression in GBM cells (Figure 4A).

FIGURE 3

Figure 3 METTL3-mediated m⁶A modification contributes to the TMZ resistance. (A) The colorimetric m⁶A quantification assay was used to examine the total m⁶A levels in the control or DAA-treated U87-MG/U251-TMZ cells. (B) The dot blot was used to confirm the total m⁶A levels in the control or DAA-treated U87-MG/U251-TMZ cells. (C) The effect of DAA treatment on the cell proliferation was analyzed by the colony formation assay. (D) The effect of DAA treatment on the sensitivity to TMZ was analyzed by CCK-8 assay. *P < 0.05 and **P < 0.01 versus blank U87-MG/U251 cells.

FIGURE 4

Figure 4 METTL3 contributes to the TMZ resistance via m⁶A modification of MGMT and ANPG mRNAs. (A) The expression mRNA level of TMZ resistant genes (ANPG, CBX5, MGMT, MSH2, MSH6, MLH1, MPG, XRCC3, and XPC) in normal and TMZ-resistant U87-MG/U251 cells were analyzed by real-time PCR. (B) The m⁶A-methylated level of MGMT and ANPG in normal and TMZ resistant U87-MG/U251 cells were analyzed by Me-RIP-real-time PCR. (C) The effect of METTL3 overexpression and knockdown on the m⁶A methylated level of MGMT and ANPG were analyzed by Me-RIP-real-time PCR. (D) The effect of DAA treatment on the m⁶A-methylated level of MGMT and ANPG were analyzed by Me-RIP-real-time PCR. *P < 0.05 and **P < 0.01 versus indicated control U87-MG/U251 cells.

Furthermore, the m⁶A methylation level (Figure 4B) of MGMT and ANPG were significantly increased in TMZ-resistant GBM cells. Notably, the m⁶A methylation level (Figure 4C) of MGMT and ANPG was significantly increased by METTL3 overexpression, which decreased by METTL3 knockdown or DAA treatment (Figure 4D). Collectively, these results demonstrate that METTL3 contributes to TMZ resistance via m⁶A modification.

METTL3-Mediated m⁶A Modification Contributes to the TMZ Resistance In Vivo

To investigate whether METTL3-mediated m⁶A modification was TMZ-resistant in vivo, we subcutaneously injected shMETTL3 or shNC-expressing U87-MG-TMZ cells into BALB/c NOD mice. After confirmation of GBM implantation, mice were treated with TMZ (66 mg/kg/d, 5 d per week, for 3 cycles). The tumor volume (Figures 5A, B) and weight (Figure 5C) of mice injected with shMETTL3 were significantly lower than those of xenografts expressing shNC. In contrast, mice treated with DAA (50 mg/kg/d, 5 d per week, for 3 cycles) and TMZ also resulted in a smaller tumor volume (Figures 5A, B) and weight (Figure 5C) than the blank group. IHC staining was performed to verify the expression of cleaved caspase-3. TMZ-treated xenografts with shMETTL3 expressing or DAA treatment had significantly increased level of cleaved caspase-3 compared with shNC or blank xenografts (Figure 5D). Taken together, these results demonstrate that METTL3-mediated m⁶A modification contributes to TMZ resistance in vivo.

FIGURE 5

Figure 5 METTL3-mediated m⁶A modification contributes to the TMZ resistance in vivo. (A) Representative images of four groups of U87-MG-TMZ cells (shNC, shMETTL3, blank, and DAA)-derived subcutaneous tumors in the presence of TMZ. (B) Growth curve of tumor xenografts originated from four groups of U87-MG-TMZ cells (shNC, shMETTL3, blank, and DAA) in the presence of TMZ, (C) Weight of tumor xenografts originated from four groups of U87-MG-TMZ cells (shNC, shMETTL3, blank, and DAA) in the absence of TMZ. (D) IHC analysis of cleaved caspase-3 in tumor xenografts originated from four groups of U87-MG-TMZ cells (shNC, shMETTL3, blank, and DAA) in the presence of TMZ. *P < 0.05 versus indicated control U87-MG/U251 cells.

Discussion

GBM is one of the most aggressive types of cancer, for which no effective way of treatment is available (1). Despite advances in the development of chemotherapeutic agents, including targeted therapies, the overall survival after diagnosis is usually less than two years (26, 27). Recent experimental and clinical studies have shown that epigenetic regulation of GBM also plays an important role in promoting tumorigenesis and the development of drug resistance (2, 28, 29). m⁶A RNA methylation is an important RNA modification that has been shown to play an important role in the genesis and development of glioblastoma (30). In this study, we investigated the potential role of m⁶A methylation modification in the regulation of TMZ resistance and the feasibility of using the M⁶A inhibitor DAA as a therapeutic candidate.

In this study, we first analyzed the level of m⁶A RNA methylation in TMZ-sensitive and TMZ-resistant GBM cells and critical role of a major m⁶A methyltransferase METTL3 in TMZ resistance. METTL3 is an effective therapeutic target for various cancers, including pancreatic cancer (31), melanoma (32), colorectal cancer (33), and lung adenocarcinoma (16). METTL3 is overexpressed in hepatic cancer cells (HCC), and is associated with poor prognosis (34). METTL3 is highly expressed in ovarian cancer, significantly correlating with ovarian cancer grade, PT status, PN/PM status, and FIGO staging (35). These studies suggest that METTL3 is a potential oncogene. METTL3 enhances the m⁶A methylation by improving the stability of SOX2 in GBM, thereby promoting the stemness of glioma stem cells (GSCs) (25). Controversially, another team found that downregulation of METTL3 significantly promoted GSC self-renewal and tumorigenesis (36). In addition, ALKBH5 reduces m⁶A modification in GSCs and plays an important role in tumorigenesis in the progression of GBM by regulating FOXM1 expression (37). These findings highlight the importance of modifying m⁶A methylation in GBM progression. However, its role in TMZ resistance in GBM remains unclear. We found no significant difference in the METTL3 expression between normal and GBM tissues, and no association was observed between its expression level and the prognosis in GBM patients (GEPIA, data not shown). However, its expression is significantly elevated in TMZ-resistant GBM cells, compared to its parent TMZ-sensitive cells. Moreover, we verified the critical role of METTL3-mediated m⁶A modification in TMZ resistance in GBM cells. Both METTL3 silencing or total methylation inhibition with DAA increased the sensitivity of GBM cells to TMZ in vitro and in vivo. Meanwhile, we discovered that METTL3 overexpression dramatically increased the m⁶A methylation of MGMT and APNG, but did not affect the level of METTL14 (Supplementary Figure 2). However, METTL3 overexpression showed no effect on the colony formation of TMZ-resistant GBM cells, suggesting that a highly expressed and super-functional role of METTL3 in TMZ-resistant GBM cells, thus further overexpression of METTLE3 increased the total m⁶A methylated mRNAs, but did not enhance the cell proliferation ability of TMZ-resistant GBM cells.

Considering the molecular mechanism underlying the resistance of glioma cells to TMZ, a DNA alkylation agent, is currently the only chemotherapeutic drug having some efficacy against GBM, accompanied by surgery and radiation therapy (28). In in vitro and animal models, TMZ resistance can be mediated by MGMT, a DNA repair protein that removes the methyl group produced by TMZ from the O6 site of guanine, which represents the most cytotoxic damage (13, 38). GBM patients with methylated MGMT promoter had an increased overall survival compared with radiotherapy alone, and responded better in combination with TMZ and radiotherapy (14). However, 50% of GBM patients with MGMT methylation promoters do not survive for 2 years, and therefore receive only moderate benefits from TMZ treatment, suggesting additional resistance factors. Similarly, GBM patients with unmethylated MGMT also showed some response to TMZ, strongly suggesting that MGMT promoter methylation was not the only predictor of response to TMZ (39). In this study, we demonstrated that the expression of MGMT mRNA is also regulated by METTL3-mediated m⁶A modification, which contributes to TMZ resistance. Moreover, our investigation into other DNA repair modulating systems, including GATA4-mediated TMZ sensitivity (40), showed increased levels of APNG in METTL3 over-expressed GBM cells. The m⁶A or total mRNA levels of MGMT and APNG were elevated by METTL3 overexpression, whereas decreased by METTL3-silencing or DAA treatment. In summary, we have demonstrated that METTL3 promotes the TMZ resistance of glioma cells by increasing MGMT and ANPG in an m⁶A-dependent manner.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Ethics Statement

The animal study was reviewed and approved by The Ethics Committee of The Second Affiliated Hospital of Soochow University.

Author Contributions

JS and GC were major contributors to the molecular experiments and manuscript writing. XD and HL performed statistical analyses of the experimental data. SL and SC cultured the cells. YL and LW performed animal experiments. JY and ZQ reviewed the manuscript. JD designed and supervised the experiments. All authors contributed to the article and approved the submitted version.

Funding

1. The Natural Science Foundation of Jiangsu Province (grant no. BK20201172). 2. Clinical Special Disease Diagnosis and Treatment Technology in Suzhou (grant no. LCZX201807). 3. Key project of the Jiangsu Health Commission (ZDB2020016). 4. Young Talent Development Plan of Changzhou Health Commission (grant no. 2020-233-CZQM2020013).

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.702983/full#supplementary-material

Supplementary Figure 1 | Generation of TMZ-resistant U87-MG and U251 cell lines. TMZ-resistant cell lines were generated by exposure of U87-MG and U251 cells to 200 μM of TMZ for over 6 months. The derived resistant cell lines were designated as U87-MG-TMZ resistant and U251-TMZ resistant, respectively. The cell survival ratio and half maximal inhibitory concentration (IC₅₀) of TMZ for U87-MG and U251 was evaluated by CCK-8 assay.

Supplementary Figure 2 | METTL3 overexpression or knockdown did not affect the level of METTL14 in the TMZ-resistance GBM cells. (A) The mRNA level of METTL14 in the METTL3 overexpression or knockdown TMZ-resistance GBM cells were analyzed by qRT-PCR. (B) The protein level of METTL14 in the METTL3 overexpression or knockdown TMZ-resistance GBM cells were analyzed by western blot.

References

1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Phillips RE, Soshnev AA, Allis CD. Epigenomic Reprogramming as a Driver of Malignant Glioma. Cancer Cell (2020) 38(5):647–60. doi: 10.1016/j.ccell.2020.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Sidaway P. Low-Grade Glioma Subtypes Revealed. Nat Rev Clin Oncol (2020) 17(6):335. doi: 10.1038/s41571-020-0380-4

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Tan AC, Ashley DM, Lopez GY, Malinzak M, Friedman HS, Khasraw M. Management of Glioblastoma: State of the Art and Future Directions. CA Cancer J Clin (2020) 70(4):299–312. doi: 10.3322/caac.21613

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Tuleasca C, Knisely J, Leroy HA, Hottinger AF, Peciu-Florianu I, Levivier M, et al. Glioma Patient-Reported Outcome Assessment in Clinical Care. Lancet Oncol (2020) 21(5):e230. doi: 10.1016/S1470-2045(20)30141-8

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Chen X, Zhang M, Gan H, Wang H, Lee JH, Fang D, et al. A Novel Enhancer Regulates MGMT Expression and Promotes Temozolomide Resistance in Glioblastoma. Nat Commun (2018) 9(1):2949. doi: 10.1038/s41467-018-05373-4

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Guo G, Gong K, Puliyappadamba VT, Panchani N, Pan E, Mukherjee B, et al. Efficacy of EGFR Plus TNF Inhibition in a Preclinical Model of Temozolomide-Resistant Glioblastoma. Neuro Oncol (2019) 21(12):1529–39. doi: 10.1093/neuonc/noz127

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Lu C, Wei Y, Wang X, Zhang Z, Yin J, Li W, et al. DNA-Methylation-Mediated Activating of lncRNA SNHG12 Promotes Temozolomide Resistance in Glioblastoma. Mol Cancer (2020) 19(1):28. doi: 10.1186/s12943-020-1137-5

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Meng X, Zhao Y, Han B, Zha C, Zhang Y, Li Z, et al. Dual Functionalized Brain-Targeting Nanoinhibitors Restrain Temozolomide-Resistant Glioma via Attenuating EGFR and MET Signaling Pathways. Nat Commun (2020) 11(1):594. doi: 10.1038/s41467-019-14036-x

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Wu L, Bernal GM, Cahill KE, Pytel P, Fitzpatrick CA, Mashek H, et al. BCL3 Expression Promotes Resistance to Alkylating Chemotherapy in Gliomas. Sci Transl Med (2018) 10(448):eaar2238. doi: 10.1126/scitranslmed.aar2238

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Yi GZ, Huang G, Guo M, Zhang X, Wang H, Deng S, et al. Acquired Temozolomide Resistance in MGMT-Deficient Glioblastoma Cells Is Associated With Regulation of DNA Repair by DHC2. Brain (2019) 142(8):2352–66. doi: 10.1093/brain/awz202

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Agnihotri S, Gajadhar AS, Ternamian C, Gorlia T, Diefes KL, Mischel PS, et al. Alkylpurine-DNA-N-Glycosylase Confers Resistance to Temozolomide in Xenograft Models of Glioblastoma Multiforme and Is Associated With Poor Survival in Patients. J Clin Invest (2012) 122(1):253–66. doi: 10.1172/JCI59334

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Bell EH, Zhang P, Fisher BJ, Macdonald DR, McElroy JP, Lesser GJ, et al. Association of MGMT Promoter Methylation Status With Survival Outcomes in Patients With High-Risk Glioma Treated With Radiotherapy and Temozolomide: An Analysis From the NRG Oncology/RTOG 0424 Trial. JAMA Oncol (2018) 4(10):1405–9. doi: 10.1001/jamaoncol.2018.1977

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Lee SY. Temozolomide Resistance in Glioblastoma Multiforme. Genes Dis (2016) 3(3):198–210. doi: 10.1016/j.gendis.2016.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Frye M, Harada BT, Behm M, He C. RNA Modifications Modulate Gene Expression During Development. Science (2018) 361(6409):1346–9. doi: 10.1126/science.aau1646

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Lin S, Choe J, Du P, Triboulet R, Gregory RI. The M(6)A Methyltransferase METTL3 Promotes Translation in Human Cancer Cells. Mol Cell (2016) 62(3):335–45. doi: 10.1016/j.molcel.2016.03.021

PubMed Abstract | CrossRef Full Text | Google Scholar

17. 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(7397):201–6. doi: 10.1038/nature11112

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Akichika S, Hirano S, Shichino Y, Suzuki T, Nishimasu H, Ishitani R, et al. Cap-Specific Terminal N (6)-Methylation of RNA by an RNA Polymerase II-Associated Methyltransferase. Science (2019) 363(6423):eaav0080. doi: 10.1126/science.aav0080

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Wang H, Hu X, Huang M, Liu J, Gu Y, Ma L, et al. Mettl3-Mediated mRNA M(6)A Methylation Promotes Dendritic Cell Activation. Nat Commun (2019) 10(1):1898. doi: 10.1038/s41467-019-09903-6

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Winkler R, Gillis E, Lasman L, Safra M, Geula S, Soyris C, et al. M(6)A Modification Controls the Innate Immune Response to Infection by Targeting Type I Interferons. Nat Immunol (2019) 20(2):173–82. doi: 10.1038/s41590-018-0275-z

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Liu J, Harada BT, He C. Regulation of Gene Expression by N(6)-Methyladenosine in Cancer. Trends Cell Biol (2019) 29(6):487–99. doi: 10.1016/j.tcb.2019.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Pan Y, Ma P, Liu Y, Li W, Shu Y. Multiple Functions of M(6)A RNA Methylation in Cancer. J Hematol Oncol (2018) 11(1):48. doi: 10.1186/s13045-018-0590-8

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Wang S, Chai P, Jia R, Jia R. Novel Insights on M(6)A RNA Methylation in Tumorigenesis: A Double-Edged Sword. Mol Cancer (2018) 17(1):101. doi: 10.1186/s12943-018-0847-4

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Visvanathan A, Patil V, Abdulla S, Hoheisel JD, Somasundaram K. N(6)-Methyladenosine Landscape of Glioma Stem-Like Cells: METTL3 Is Essential for the Expression of Actively Transcribed Genes and Sustenance of the Oncogenic Signaling. Genes (Basel) (2019) 10(2):141. doi: 10.3390/genes10020141

CrossRef Full Text | Google Scholar

25. Visvanathan A, Patil V, Arora A, Hegde AS, Arivazhagan A, Santosh V, et al. Essential Role of METTL3-Mediated M(6)A Modification in Glioma Stem-Like Cells Maintenance and Radioresistance. Oncogene (2018) 37(4):522–33. doi: 10.1038/onc.2017.351

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Cheng X, Geng F, Pan M, Wu X, Zhong Y, Wang C, et al. Targeting DGAT1 Ameliorates Glioblastoma by Increasing Fat Catabolism and Oxidative Stress. Cell Metab (2020) 32(2):229–42.e8. doi: 10.1016/j.cmet.2020.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Gao X, Li S, Ding F, Liu X, Wu Y, Li J, et al. A Virus-Mimicking Nucleic Acid Nanogel Reprograms Microglia and Macrophages for Glioblastoma Therapy. Adv Mater (2021) 33(9):e2002116. doi: 10.1002/adma.202006116

CrossRef Full Text | Google Scholar

28. Kamson DO, Grossman SA. The Role of Temozolomide in Patients With Newly Diagnosed Wild-Type IDH, Unmethylated MGMTp Glioblastoma During the COVID-19 Pandemic. JAMA Oncol (2021) 7(5):675–6. doi: 10.1001/jamaoncol.2020.6732

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Luo X, Weiss WA. Utility of Human-Derived Models for Glioblastoma. Cancer Discov (2020) 10(7):907–9. doi: 10.1158/2159-8290.CD-20-0493

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Dixit D, Prager BC, Gimple RC, Poh HX, Wang Y, Wu Q, et al. The RNA M6a Reader YTHDF2 Maintains Oncogene Expression and Is a Targetable Dependency in Glioblastoma Stem Cells. Cancer Discov (2021) 11(2):480–99. doi: 10.1158/2159-8290.CD-20-0331

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Zhang J, Bai R, Li M, Ye H, Wu C, Wang C, et al. Excessive miR-25-3p Maturation via N(6)-Methyladenosine Stimulated by Cigarette Smoke Promotes Pancreatic Cancer Progression. Nat Commun (2019) 10(1):1858. doi: 10.1038/s41467-019-09712-x

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Wang L, Hui H, Agrawal K, Kang Y, Li N, Tang R, et al. M(6) A RNA Methyltransferases METTL3/14 Regulate Immune Responses to Anti-PD-1 Therapy. EMBO J (2020) 39(20):e104514. doi: 10.15252/embj.2020104514

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Li T, Hu PS, Zuo Z, Lin JF, Li X, Wu QN, et al. METTL3 Facilitates Tumor Progression via an M(6)A-IGF2BP2-Dependent Mechanism in Colorectal Carcinoma. Mol Cancer (2019) 18(1):112. doi: 10.1186/s12943-019-1038-7

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Zuo X, Chen Z, Gao W, Zhang Y, Wang J, Wang J, et al. M6A-Mediated Upregulation of LINC00958 Increases Lipogenesis and Acts as a Nanotherapeutic Target in Hepatocellular Carcinoma. J Hematol Oncol (2020) 13(1):5. doi: 10.1186/s13045-019-0839-x

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Liang S, Guan H, Lin X, Li N, Geng F, Li J. METTL3 Serves an Oncogenic Role in Human Ovarian Cancer Cells Partially via the AKT Signaling Pathway. Oncol Lett (2020) 19(4):3197–204. doi: 10.3892/ol.2020.11425

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Cui Q, Shi H, Ye P, Li L, Qu Q, Sun G, et al. M(6)A RNA Methylation Regulates the Self-Renewal and Tumorigenesis of Glioblastoma Stem Cells. Cell Rep (2017) 18(11):2622–34. doi: 10.1016/j.celrep.2017.02.059

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Zhang S, Zhao BS, Zhou A, Lin K, Zheng S, Lu Z, et al. M(6)A Demethylase ALKBH5 Maintains Tumorigenicity of Glioblastoma Stem-Like Cells by Sustaining FOXM1 Expression and Cell Proliferation Program. Cancer Cell (2017) 31(4):591–606.e6. doi: 10.1016/j.ccell.2017.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Lamba N, Chukwueke UN, Smith TR, Ligon KL, Aizer A, Reardon DA, et al. Socioeconomic Disparities Associated With MGMT Promoter Methylation Testing for Patients With Glioblastoma. JAMA Oncol (2020) 6(12):1972–4. doi: 10.1001/jamaoncol.2020.4937

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Killock D. Lomustine-Temozolomide Combination Efficacious in Newly Diagnosed Glioblastoma. Nat Rev Clin Oncol (2019) 16(5):273. doi: 10.1038/s41571-019-0192-6

CrossRef Full Text | Google Scholar

40. Fosmark S, Hellwege S, Dahlrot RH, Jensen KL, Derand H, Lohse J, et al. APNG as a Prognostic Marker in Patients With Glioblastoma. PloS One (2017) 12(6):e0178693. doi: 10.1371/journal.pone.0178693

PubMed Abstract | CrossRef Full Text | Google Scholar