miR-130a can predict response to temozolomide in patients with glioblastoma multiforme, independently of O6-methylguanine-DNA methyltransferase
© Chen et al.; licensee BioMed Central. 2015
Received: 16 November 2014
Accepted: 11 February 2015
Published: 21 February 2015
Currently, O6-methylguanine-DNA methyltransferase(MGMT) promoter methylation is the most convincing predictive biomarker for temozolomide (TMZ) response in patients with glioblastoma multiforme (GBM). However, technical obstacles prevent this biomarker from being applied widely. On the other hand, microRNAs (miRNAs) are easily investigated in the clinical setting using quantitative real-time polymerase chain reactions. This study aimed to identify miRNAs that could serve as predictive biomarkers for TMZ response.
The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA) databases were used to investigate the significance of associations between miRNA expression and overall survival (OS) in TMZ-treated patients with GBM. Cytotoxicity assays were used to validate the miRNAs’ roles in the response of glioma cells to TMZ. Biological insights concerning the miRNAs were explored using gene set enrichment analysis (GSEA) and gene ontology (GO) analysis.
miR-130a was found to be significantly associated with OS in TMZ-treated patients from TCGA and the CGGA. In contrast, miR-130a appeared to be unassociated with OS in patients who only received radiotherapy. The TMZ cytotoxicity assay showed that miR-130a over-expression could sensitize response to TMZ in glioma cells. GSEA and GO analysis indicated that lower miR-130a could generate a more extensive response to oxidative stress, which in turn could elevate Ape1 and mediate resistance to TMZ. In vitro experiment verified that cells with lower miR-130a express higher Ape1 under oxidative stress.
Our data suggested that miR-130a could be a predictive marker for TMZ response in patients with GBM, independently of the mechanism by which MGMT acts as a biomarker. miR-130a could serve as a guide for treatment strategy selection in cases of GBM.
Glioblastoma multiforme (GBM) is the most malignant and common brain tumor in adults. Temozolomide (TMZ) is currently the first-line chemotherapy for this disease. TMZ can cross the blood–brain barrier and significantly improves the prognosis of patients with GBM . However, survival remains poor among these patients because of resistance to TMZ. Biomarkers have been developed to predict patients’ responses to TMZ.
O6-methylguanine-DNA methyltransferase (MGMT) promoter methylation is currently the most convincing predictive biomarker for TMZ response, based on its correlation with overall survival (OS) in patients with GBM who have received TMZ chemotherapy [2-4]. However, this biomarker remains unsuitable for guiding GBM treatment in widespread application because of technical obstacles .
MicroRNAs (miRNAs) are small non-coding RNAs that are more stable than messenger RNAs. Since quantitative real-time polymerase chain reaction (qRT-PCR) is widely used to quantify RNA expression levels in clinical practice, it would be convenient to have a miRNA that could be used to predict response to TMZ treatment in patients with GBM. Previous studies [5-11] reported several kinds of miRNAs which have predictive value for TMZ response in glioblastoma. These studies mainly focused on one special target of a miRNA. That may not explain the whole function of the miRNA, since one kind of miRNA can regulate hundreds of genes. Bioinformatic analyses, on the other hand, could illustrate the miRNA function from the macroscopic view. In this study, using bioinformatic methods, we sought to identify such a miRNA, finding evidence that miR-130a is correlated with OS in TMZ-treated patients with GBM. Further, the biological process involving miR-130a is independent of MGMT.
Material and methods
The Cancer Genome Atlas (TCGA) database  included miRNA expression profiling of 184 adult patients who had GBM and underwent radiation therapy and TMZ chemotherapy. Data on these patients were used as a discovery set. We downloaded the mRNA expression profiling and clinical information on these patients, as well as the miRNA expression profiling of another 128 glioblastoma patients in TCGA who only received radiation therapy. Additionally, we obtained miRNA expression profiling of 28 adult patients who had glioblastoma and underwent TMZ chemotherapy from the Chinese Glioma Genome Atlas (CGGA) database . Data on these 28 patients were used as a validation set.
Human glioma cell lines U87, LN229, SHG44, and U251 were cultured in DMEM/F-12 culture medium (Gibco, Life Technologies, Carlsbad, CA, USA) supplemented with 10% fetal bovine serum (Gibco). Cells were maintained in a humidified incubator with 5% CO2 at 37°C. The All-in-One miRNA qRT-PCR Detection Kit (FulenGen, Guangzhou, China) was used to evaluate the miR-130a levels in cells.
Regarding miR-130a mimics and control miRNA, inhibitors for miRNA-130a and controls were purchased from FulenGen in China. Transfection was performed using Lipofectamine RNAiMAX Transfection Reagent (Invitrogen, Life Technologies, Carlsbad, CA, USA), following the manufacturer’s instructions. Transfection efficiency was also evaluated using the All-in-One miRNA qRT-PCR Detection Kit.
TMZ cytotoxicity assay
After transfection, cells were seeded at a density of 6000 cells per well on 96 well plates in triplicates. TMZ (Sigma Aldrich, St. Louis, MO, USA) was added at various concentrations ranging from 0 to 1000 μmol/L for 48 hours. The Cell Titer 96 AQueous One Solution Cell Proliferation Assay kit (Promega, Madison, WI, USA) was used to determine the number of viable cells in the TMZ cytotoxicity assay. To test the cleavage of nuclear poly (ADP-ribose) polymerase (PARP), proteins of SHG44 and U251 were extracted after transfection and 96 hours 500 μM TMZ treatment.
Induction of Oxidative Stress with hydrogen peroxide
At 48 hours after transfection, cells were washed twice with PBS, then incubated in complete medium supplemented with 100 μM hydrogen peroxide to induce oxidative stress. The incubation continued for 24 hours before determination of Ape1 expression using western blotting.
Proteins were extracted using RIPA buffer (Cell Signaling Technology) supplemented with protease inhibitor. Proteins were separated by 10% SDS-PAGE and transferred to PVDF membrane. After blocking with 5% skim milk in TBST, the membranes were incubated with primary antibodies overnight at 4°C. The membranes were then washed and incubated with secondary antibodies for 1 hour. Antibody binding was assessed by enhanced chemiluminescence substrate (BioRad). Antibodies to PARP and Ape1 were purchased from Proteintech in China.
For each miRNA, the association between the expression level and OS in the discovery set was calculated using a univariate Cox proportional hazards regression analysis. The 4 miRNAs that were found to be significantly associated with OS in the discovery set were then evaluated in the validation set using the Kaplan-Meier method and 2-sided log-rank tests. Subsequently, the miRNA that was significantly associated with OS in both sets was further analyzed using a Cox multivariate analysis that included O-6-methylguanine-DNA methyltransferase (MGMT) methylation status, age, and gender as covariates. The association between the expression of this miRNA and OS was also investigated in the group of patients from TCGA who had received radiation therapy alone. As previously, the Kaplan-Meier method and 2-sided log-rank tests were also used in this analysis. The cutoff value for this miRNA in TCGA was defined as the median expression level.
The miRNA-associated gene set was identified using a Spearman correlation analysis. The negatively correlated genes were included in a gene ontology (GO) analysis that was performed using DAVID ; (National Institute of Allergy and Infectious Diseases, Bethesda, MD, USA). Additionally, we performed a Gene Set Enrichment Analysis (GSEA) with GSEA software downloaded from the Broad Institute [15,16] (Cambridge, MA, USA).
Paired-sample t-tests were employed to assess the significance of differences between TMZ cytotoxicity assay results.
Statistical analyses of TCGA and CGGA datasets were performed using BRB-Array Tools 4.3.0 Beta_2, which was developed by Dr. Richard Simon and the BRB-Array Tools Development Team (National Cancer Institute, Bethesda, MD, USA). Statistical analyses for in vitro experimental results were performed using SPSS 16.0 (IBM SPSS, Inc., Chicago, IL, USA). Two-sided p-values less than 0.05 were regarded as statistically significant.
miR-130a was correlated with overall survival in TMZ-treated patients with GBM, but not in non-TMZ-treated patients
miRNAs correlated with OS of patients treated with temozolomide in TCGA
2.48 ∙ 10 −5
3.65 ∙ 10 −5
Cox regression analyses of the associations of miR-130a and clinical characteristics with OS in TCGA
Univariate Cox regression
Multivariate Cox regression
Further, the Kaplan-Meier method and 2-sided log-rank tests showed that miR-130a was not significantly associated with OS in the patients who only received radiation therapy (Figure 1F).
Biological insights regarding miR-130a
In addition, GSEA showed enrichment of genes related to response to oxidative stress (ROS) among patients with low miR-130a expression (normalized enrichment score [NES] =1.51, p = 0.037).
The miR-130a levels were relatively high in U251 cells (Figure 3A). This cell line was then transfected with inhibitors (5′-AUGCCCUUUUAACAUUGCACUG-3′) for miR-130a and the scrambled control (Figure 3D). Western blotting showed that U251 transfected with the scrambled control (U251-SC) expressed lower Ape1 than cells transfected with miR-130a inhibitors (U251-I) did under oxidative stress (Figure 2E).
miR-130a over expression sensitized the response of glioma cells to TMZ
The TMZ cytotoxicity assay showed that SHG44-M were more sensitive to TMZ than SHG44-SC (p = 0.027, Figure 3C). Western blotting showed that overexpression of miR-130a increased PARP cleavage in SHG44 induced by TMZ (Figure 3C).
In our univariate Cox analyses of TCGA data, miR-130a, miR-20a, miR-221, and miR-222 showed significant associations with the OS of patients with GBM who received TMZ chemotherapy. However, only the predictive capacity of miR-130a was validated by an analysis of CGGA data on patients with GBM who received TMZ chemotherapy. In a further analysis of TCGA data, we also observed that miR-130a failed to correlate with the OS of patients with GBM who received radiation therapy alone.
The role of miR-130a in the chemo-resistance of GBM cells was also determined in the present study. SHG44 with up-regulated miR-130a showed higher sensitivity than parental cells. On the other hand, U251 with down-regulated miR-130a showed lower sensitivity than parental cells. The cytotoxicity of TMZ was less than 50% in these investigations because of the short duration of treatment. Further in vitro functional research on miR-130a is necessary.
On the whole, higher miR-130a expression was associated with prolonged OS in patients who had GBM and received TMZ chemotherapy. miR-130a could be a candidate predictive biomarker for TMZ response.
miR-130a has been reported to reduce resistance to Gefitinib and TNF-related apoptosis-inducing ligand (TRAIL) in Non-small cell lung cancer by down-regulating MET and miRNA-221/222 [17,18]. To investigate the biological mechanism of miR-130a in GBM, we performed GSEA and GO analysis of miR-130a-associated genes. Oxidation reduction and ROS were showed to be involved in miR-130a-correlated biological processes.
MGMT promotes clinical resistance to chemotherapy through its ability to remove the O6-methylguanine (O6-meG) adduct produced by TMZ [19,20]. In fact, TMZ also generates several other kinds of DNA adducts in cells. N7-methylguanine (7-meG) and N3-methyladenine (3-meA) are the most abundant of these DNA adducts . They are precursors of abasic sites, which can impede DNA replication [21-23]. The base excision repair (BER) pathway is responsible for repairing these lesions. Apurinic/apyrimidinic endonuclease 1 (Ape1) is the key enzyme in BER pathways that remove abasic sites induced by TMZ [24,25]. The function of Ape1, which is independent of MGMT, can lead to cell resistance to TMZ treatment . A small-molecule inhibitor of Ape1 has been found to block proliferation and reduces viability of glioblastoma cells . It has been reported that oxidative stress could elevate Ape1, which in turn mediates chemoresistance by repairing abasic sites produced by TMZ [28,29].
In summary, we have reported that miR-130a could be a predictive marker for TMZ response in patients with GBM, and that this predictive ability is independent of the predictive ability of MGMT. Further, the miR-130a biomarker is more convenient to test than MGMT methylation. Its application could provide guidance when selecting treatment strategies for GBM in the future.
We are grateful to National Natural Science Foundation of China for funding and support. We are grateful to Dr. Jiang Tao for the use of his laboratory resources during the course of this work. We also thank to Drs. Zhang Chuanbao and Wang Yinyan for help in acquisition and interpretation of the data.
- Stupp R, Mason WP, van den Bent MJ, Weller M, Fisher B, Taphoorn MJ. Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. N Engl J Med. 2005;352(10):987–96.View ArticlePubMedGoogle Scholar
- Hegi ME, Liu L, Herman JG, Stupp R, Wick W, Weller M, et al. Correlation of O6-methylguanine methyltransferase (MGMT) promoter methylation with clinical outcomes in glioblastoma and clinical strategies to modulate MGMT activity. J Clin Oncol. 2008;26(25):4189–99.View ArticlePubMedGoogle Scholar
- Weller M, Stupp R, Reifenberger G, Brandes AA, van den Bent MJ, Wick W, et al. MGMT promoter methylation in malignant gliomas: ready for personalized medicine? Nat Rev Neurol. 2010;6:39–51.View ArticlePubMedGoogle Scholar
- von Deimling A, Korshunov A, Hartmann C. The next generation of glioma biomarkers: MGMT methylation, BRAF fusions and IDH1 mutations. Brain Pathol. 2011;21:74–87.View ArticleGoogle Scholar
- Hao W, Qiang L, Tao C, Yu-dan C, Fan L, Zhi-fei W. MiR-136 modulates glioma cell sensitivity to temozolomide by targeting astrocyte elevated gene-1. Diagn Pathol. 2014;9:173.View ArticleGoogle Scholar
- Cristina Q, Davide M, Giuseppina R, Giulia R, Angel DL, Margherita I, et al. miR-221/222 Target the DNA Methyltransferase MGMT in Glioma Cells. Plos One. 2013;8(9):e74466.View ArticleGoogle Scholar
- Zhang W, Zhang J, Hoadley K, Kushwaha D, Ramakrishnan V, Li S, et al. miR-181d: a predictive glioblastoma biomarker that downregulates MGMT expression. Neuro-Oncology. 2012;14(6):712–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Lan J, Xue Y, Chen H, Zhao S, Zhijian W, Fang J, et al. Hypoxia-induced miR-497 decreases glioma cell sensitivity to TMZ by inhibiting apoptosis. FEBS Lett. 2014;588:3333–9.View ArticlePubMedGoogle Scholar
- Deepa K, Valya R, Kimberly N, Tyler S, Thien N, Diahnn F, et al. A genome-wide miRNA screen revealed miR-603 as a MGMT-regulating miRNA in glioblastomas. Oncotarget. 2014;5(12):4026–39.Google Scholar
- Qiang Liu, Ran Zou, Rouxi Zhou, Chaofan Gong, Zhifei Wang, Tao Cai, Chaochao Tan, and Jiasheng Fang. miR-155 Regulates Glioma Cells Invasion and Chemosensitivity by p38 Isforms in vitro. Journal of Cellular Biochemistry 2014, doi: [10.1002/jcb.25073].Google Scholar
- Yan-chang S, Jing W, Cheng-cheng G, Ke S, Jian W, Fu-rong C, et al. miR-181b sensitizes glioma cells to teniposide by targeting MDM2. BMC Cancer. 2014;14:611.View ArticleGoogle Scholar
- The Cancer Genome Atlas database [http://cancergenome.nih.gov]
- Chinese Glioma Genome Atlas database [http://www.cgga.org.cn]
- DAVID [http://david.abcc.ncifcrf.gov/home.jsp]
- Gene Set Enrichment Analysis [http://www.broadinstitute.org/gsea/]
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.View ArticlePubMed CentralPubMedGoogle Scholar
- Acunzo M, Visone R, Romano G, Veronese A, Lovat F, Palmieri D, et al. miR-130a targets MET and induces TRAIL-sensitivity in NSCLC by downregulating miR-221&222. Oncogene. 2012;31(5):634–42.PubMed CentralPubMedGoogle Scholar
- Zhou Y-M, Liu J, Sun W. MiR-130a Overcomes Gefitinib Resistance by Targeting Met in Non-Small Cell Lung Cancer Cell Lines. Asian Pacific J Cancer Prevent. 2014;15(3):1391–6.View ArticleGoogle Scholar
- Mitra S. MGMT: a personal perspective. DNA Repair. 2007;6:1064–70.View ArticlePubMed CentralPubMedGoogle Scholar
- Tubbs JLJL, Pegg AE, Tainer JA. DNA binding, nucleotide flipping, and the helix-turn-helix motif in base repair by O6-alkylguanine-DNA alkyltransferase and its implications for cancer chemotherapy. DNA Repair. 2007;6:1100–15.View ArticlePubMed CentralPubMedGoogle Scholar
- Loeb LA, Preston BD. Mutagenesis by apurinic/apyrimidinic sites. Annu Rev Genet. 1986;20:201–30.View ArticlePubMedGoogle Scholar
- Sutton MD, Walker GC. Managing DNA polymerases: coordinating DNA replication, DNA repair, and DNA recombination. Proc Natl Acad Sci U S A. 2001;98:8342–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Yuan F, Zhang Y, Rajpal DK, Wu X, Guo D, Wang M, et al. Specificity of DNA lesion bypass by the yeast DNA polymerase η. J Biol Chem. 2000;275:8233–9.View ArticlePubMedGoogle Scholar
- Evans AR, Limp-Foster M, Kelley MR. Going APE over ref-1. Mutat Res. 2000;461:83–108.View ArticlePubMedGoogle Scholar
- Wilson III DM, Barsky D. The major human abasic endonuclease: formation, consequences, and repair of abasic sites in DNA. Mutat Res. 2001;485:283–307.View ArticlePubMedGoogle Scholar
- Bobola MS, Blank A, Berger MS, Stevens BA, Silber JR. Apurinic/apyrimidinic endonuclease activity is elevated in human adult gliomas. Clin Cancer Res. 2001;7:3510–8.PubMedGoogle Scholar
- Bapat A, Glass LS, Luo M, Fishel ML, Long EC, Georgiadis MM, et al. Novel small-molecule inhibitor of apurinic/apyrimidinic endonuclease 1 blocks proliferation and reduces viability of glioblastoma cells. J Pharmacol Exp Ther. 2010;334(3):988–98.View ArticlePubMed CentralPubMedGoogle Scholar
- Silber JR, Bobola MS, Blank A, Schoeler KD, Haroldson PD, Huynh MB, et al. The Apurinic/Apyrimidinic Endonuclease Activity of Ape1/Ref-1 Contributes to Human Glioma Cell Resistance to Alkylating Agents and Is Elevated by Oxidative Stress. Clin Cancer Res. 2002;8:3008–18.PubMedGoogle Scholar
- Bobola MS, Kolstoe DD, Blank A, Chamberlain MC, Silber JR. Repair of 3-methyladenine and abasic sites by base excision repair mediates glioblastoma resistance to temozolomide. Frontiers Oncol. 2012;2:176. 1–8.View ArticleGoogle Scholar
- Layton ME, Pazdernik TL. Reactive oxidant species in piriform cortex extracellular fluid during seizures induced by systemic kainic acid in rats. J Mol Neurosci. 1999;13:63–8.View ArticlePubMedGoogle Scholar
- Ziche M, Morbidelli L. Nitric oxide and angiogenesis. J NeuroOncol. 2000;50:139–48.View ArticlePubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.