- Open Access
miR-130a can predict response to temozolomide in patients with glioblastoma multiforme, independently of O6-methylguanine-DNA methyltransferase
Journal of Translational Medicinevolume 13, Article number: 69 (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
Cox univariate analyses showed that miR-130a, miR-20a, miR-221, and miR-222 were correlated with OS in TMZ-treated patients with GBM from TCGA (Table 1, Figure 1A,B,C and D). These 4 miRNAs were then evaluated in the validation dataset using the Kaplan-Meier method and 2-sided log-rank tests. miR-130a was found to be significantly correlated with OS in both of the datasets (Figure 1E).
MGMT methylation status, age, and gender were analyzed using univariate Cox analyses. MGMT methylation status and age were found to be significantly associated with OS (Table 2). miR-130a was then analyzed in the group of TMZ-treated GBM patients from TCGA using a multivariate Cox analysis that included MGMT methylation status and age as covariates (Table 2). miR-130a remained significantly associated with OS in this analysis.
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
GO analysis and GSEA were conducted to provide biological insights concerning the role of miR-130a. The mRNA expression profiling data were used to find miR-130a-correlated genes, as assessed using Spearman’s correlation coefficient. The eighty-six genes (Figure 2A) that were negatively correlated with miR-130a were then included in GO analyses. The top 10 GO terms are presented in Figure 2B. As shown in the figure, miR-130a-related genes had a relatively tight association with oxidation reduction.
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 low in SHG44 cells (Figure 3A). This cell line was then transfected with miR-130a mimics and the scrambled control (Figure 3B). Western blotting showed that SHG44 transfected with the scrambled control (SHG44-SC) expressed higher Ape1 than cells transfected with miR-130a mimics (SHG44-M) did under oxidative stress (Figure 2D).
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].
GO analysis and GSEA showed a negative correlation between response to oxidation stress and miR-130a. Patients with lower miR-130a expression possess a higher ability to respond to oxidative stress. In the present study, we found that cells with lower miR-130a could express higher Ape1 under oxidative stress. Several factors could induce oxidative stress in glioma, including seizures , angiogenesis associated nitric oxide , surgery associated inflammation, and radiotherapy-delivered reactive oxygen species. We hypothesize that, in the presence of these factors, patients with lower miR-130a could generate a more extensive response to oxidative stress, which in turn could elevate Ape1 and the repair of abasic sites, finally mediating resistance to TMZ (Figure 4). Because Ape1-repaired DNA adducts are different from O6-mG, this process is independent of MGMT. It appears that miR-130a was an even better predictive marker than MGMT methylation in the present study, since the former was significantly correlated with OS in a multivariate Cox analysis, while the latter was not (Table 2).
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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].
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.
The Cancer Genome Atlas database [http://cancergenome.nih.gov]
Chinese Glioma Genome Atlas database [http://www.cgga.org.cn]
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.
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.
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.
Mitra S. MGMT: a personal perspective. DNA Repair. 2007;6:1064–70.
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.
Loeb LA, Preston BD. Mutagenesis by apurinic/apyrimidinic sites. Annu Rev Genet. 1986;20:201–30.
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.
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.
Evans AR, Limp-Foster M, Kelley MR. Going APE over ref-1. Mutat Res. 2000;461:83–108.
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.
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.
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.
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.
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.
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.
Ziche M, Morbidelli L. Nitric oxide and angiogenesis. J NeuroOncol. 2000;50:139–48.
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.
The authors declare that they have no competing interests.
HC carried out the statistical analysis and the experiments, and drafted the manuscript. XL contributed in manuscript preparation. WL and HZ designed the research and revised the manuscript. All authors read and approved the final manuscript.