PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of oncotargetLink to Publisher's site
 
Oncotarget. 2016 March 15; 7(11): 12937–12950.
Published online 2016 February 12. doi:  10.18632/oncotarget.7346
PMCID: PMC4914333

MicroRNA-138 promotes acquired alkylator resistance in glioblastoma by targeting the Bcl-2-interacting mediator BIM

Abstract

Glioblastoma is the most aggressive brain tumor in adults with a median survival below 12 months in population-based studies. The main reason for tumor recurrence and progression is constitutive or acquired resistance to the standard of care of surgical resection followed by radiotherapy with concomitant and adjuvant temozolomide (TMZ/RT→TMZ). Here, we investigated the role of microRNA (miRNA) alterations as mediators of alkylator resistance in glioblastoma cells. Using microarray-based miRNA expression profiling of parental and TMZ-resistant cultures of three human glioma cell lines, we identified a set of differentially expressed miRNA candidates. From these, we selected miR-138 for further functional analyses as this miRNA was not only upregulated in TMZ-resistant versus parental cells, but also showed increased expression in vivo in recurrent glioblastoma tissue samples after TMZ/RT→TMZ treatment. Transient transfection of miR-138 mimics in glioma cells with low basal miR-138 expression increased glioma cell proliferation. Moreover, miR-138 overexpression increased TMZ resistance in long-term glioblastoma cell lines and glioma initiating cell cultures. The apoptosis regulator BIM was identified as a direct target of miR-138, and its silencing mediated the induced TMZ resistance phenotype. Altered sensitivity to apoptosis played only a minor role in this resistance mechanism. Instead, we identified the induction of autophagy to be regulated downstream of the miR-138/BIM axis and to promote cell survival following TMZ exposure. Our data thus define miR-138 as a glioblastoma cell survival-promoting miRNA associated with resistance to TMZ therapy in vitro and with tumor progression in vivo.

Keywords: chemoresistance, temozolomide, glioblastoma, miR-138, BIM

INTRODUCTION

Glioblastomas are the most aggressive primary brain tumors associated with a rapid and invariably fatal clinical course [1]. The current standard of care includes surgical resection followed by radiotherapy with concomitant and adjuvant chemotherapy with the alkylating agent temozolomide (TMZ/RT→TMZ) [2, 3]. Although TMZ prolongs progression-free and overall survival in glioblastoma patients, recurrence is inevitable because of infiltrative growth and inherent resistance to genotoxic stress. Even glioblastomas with favorable molecular profiles characterized by isocitrate dehydrogenase 1 or 2 mutation or O6-methylguanine DNA methyltransferase (MGMT) promoter methylation are prone to fail therapy.

MicroRNAs (miRNAs) are short, non-coding, single-stranded RNA molecules that function as negative post-transcriptional regulators of gene expression by binding to the 3′-untranslated region (3′UTR) of target mRNA, causing its cleavage or translational repression [4]. MiRNAs thereby regulate important cellular features such as growth, survival and death. Their expression is often deregulated in pathological conditions including cancer, where aberrant miRNA expression contributes to proliferation, invasion and resistance to apoptosis [57]. MiRNAs associated with TMZ resistance in glioblastoma include miR-195 [8], miR-9 [9], miR-125b [10], miR-136 [11], miR-181b [12], and miR-221/222 [13].

We have established a model of acquired resistance to TMZ in three well-characterized long-term glioblastoma cell lines (LTC) and investigated several candidate pathways of resistance to TMZ [14]. Here, we report a microarray-based screen for differential miRNA expression in pairs of parental and TMZ-resistant cultures. We identified a set of differentially expressed miRNA candidates and selected miR-138 for further analysis. miR-138 has been implicated in the pathogenesis of several cancers [1517], including gliomas [18, 19], however, it has not been related to therapy resistance. Here, we demonstrate that miR-138 may mediate TMZ resistance by targeting Bcl-2-like protein 11 (BCL2L11), also known as BIM, a BH3-only apoptosis initiator from the Bcl-2 protein family [20].

RESULTS

miR-138 is upregulated in glioma cells with acquired TMZ resistance and in recurrent glioblastomas in vivo

Microarray-based miRNA expression profiling of parental and TMZ-resistant (TMZR) LN-18, LN-229 and LN-308 cells revealed several differentially expressed MiRNAs (Figure (Figure1A),1A), including several miRNAs previously implicated in TMZ resistance, such as mir-125b [10], miR-181 [12] or miR-221/222 [13]. We focused on MiRNAs upregulated in LN-308 cells since the mechanism of TMZ resistance in this cell line remained unclear [14]. We selected miR-138 for further investigations which showed increased expression in LN-18-resistant (LN-18_R) and LN-308-resistant (LN-308_R) cells (Figure (Figure1B),1B), and confirmed upregulation in these cells versus parental cells by real-time qPCR (Figure (Figure1C).1C). Next, we screened the baseline expression of miR-138 in a panel of 9 LTC and 5 GIC (glioma-initiating cells) lines (Figure (Figure1D).1D). There was no significant correlation of baseline miR-138 expression with sensitivity to TMZ in clonogenic survival assays either for all cell lines pooled (r = 0.25, p = 0.39) (Figure S1A), LTC alone (r = 0.03, p = 0.95) (Figure S1B), or GIC alone (r = 0.3. p = 0.68) (Figure S1C). We confirmed that the expression of miR-138 was increased in 9 of 10 paired tissues from primary and recurrent glioblastomas following TMZ/RT→TMZ, corroborating the potential significance of this miRNA in human glioblastoma in vivo (Figure (Figure1E1E).

Figure 1
mir-138 is significantly upregulated in TMZ-resistant glioma cell lines and in recurrent glioblastoma

miR-138 overexpression increases glioma cell proliferation

Next we assessed the effects of miR-138 overexpression in LN-308 and LN-319 cells selected for their low baseline expression of miR-138 (Figure (Figure1D,1D, Figure Figure2A).2A). miR-138 induced proliferation of LN-308 glioma cells as indicated by BrdU incorporation assay (Figure (Figure2B),2B), confirmed by trypan blue dye exclusion assays over time in both cell lines (Figure (Figure2C).2C). Cell cycle analysis by flow cytometry showed a reduced G2/M fraction in the miR-138 overexpressing cells at day 6 post transfection (Figure (Figure2D2D).

Figure 2
miR-138 induces proliferation of glioma cells in vitro

Overexpression of miR-138 induces resistance to TMZ and CCNU in glioma cells

In LN-308 and ZH-305 cells, chosen for their low baseline expression of miR-138, overexpression of miR-138 conferred protection from TMZ. The EC50 shifted from 40 to 120 μM in LN-308 and from 8 to 16 μM in ZH-305 (Figure (Figure3A).3A). In contrast, no such effect was seen in LN-18 or LN-428 cells (Figure S2A), which have higher baseline expression of miR-138 (Figure (Figure1D).1D). When miR-138-overexpressing LN-308 and ZH-305 cells were exposed to increasing concentrations of another alkylating agent, CCNU, a gain of resistance was also achieved (Figure (Figure3B).3B). No cross-resistance was observed to irradiation, to the topoisomerase II inhibitor etoposide, or to the vinca alkaloid vincristine (Figure S2B). To address whether the TMZ resistant phenotype is reversed by downregulation of endogeneous miR-138, we transiently transfected miR-138-expressing LN-18 and LN-229 cells with a specific miR-138 inhibitor or a non-targeting negative control molecule. Transfection efficacy was assessed by real-time qPCR (Figure S3A) and by fluorescent light microscopy (Figure S3B). Immunoblot analysis for expression of the miR-138 target SOX4 [17] in LN-229 cells confirmed the biological activity of miR-138 inhibition (Figure S3C). Conversely, the overexpression of miR-138 decreased SOX4 levels by 40% in LN-308 cells selected for their low baseline miR-138 expression. However, the inhibitor did not increase sensitivity to TMZ in LN-18, LN-18_R, LN-229 or LN-229_R cells (Figure S3D). We assume that upon miR-138 inhibitor application, residual endogenous levels of miR-138 are still high enough to preserve the primary resistant phenotype.

Figure 3
Overexpression of miR-138 confers TMZ and CCNU resistance in glioma cells

Predicted miR-138 targets in mediating TMZ sensitivity – no role for activated leukocyte cell adhesion molecule (ALCAM)

ALCAM, an in silico predicted target of miR-138 [21], emerged as a candidate decreased in LN-18_R and LN-229_R cells from the microarray-based mRNA expression profiling (data not shown). According to The Cancer Genome Atlas (TCGA) online dataset [22], low expression levels of ALCAM correlate with shorter overall survival of glioblastoma patients (n = 504; average cut-off = 590.3; p = 5.2e−03) (Figure S4A). We confirmed a downregulation at protein level by immunoblot in LN-18_R and LN-229_R cells (Figure S4B). However, siRNA-mediated gene silencing of ALCAM (Figure S4C) did not confer TMZ resistance (Figure S4D), which indicates that ALCAM downregulation alone is not sufficient to mediate TMZ resistance in glioma cells.

BIM is a direct target of miR-138 and modulates TMZ resistance

We further investigated additional predicted miR-138 targets focusing on those that were found to be downregulated at the mRNA level by Affymetrix gene chip analyses of the three parental and resistant glioma cell lines. The predicted target BCL2L11 (or BIM) was found to be downregulated in all three resistant cell lines in the Affymetrix array. Accordingly, in the TCGA online dataset, low expression levels of BCL2L11/BIM correlate with shorter overall survival of glioblastoma patients (Figure S5). Down-regulation of BIM protein was confirmed by immunoblot analysis in LN-229_R and LN-308_R cells (Figure (Figure4A).4A). BIM has three predominant isoforms generated by alternative splicing, BIMEL (extra-large), BIML (large) and BIMs (small). These isoforms are ubiquitously expressed with a tissue-specific variation, where BIMEL is the most abundant isoform. Only BIMEL was detected in our glioma cells. This case of preferential expression of one BIM isoform in different cell types/tissues has been previously described [23, 24]. BIM levels were reduced or increased, respectively, in glioma cells transfected with miR-138 mimic or miR-138 inhibitor (Figure (Figure4B).4B). Luciferase reporter assays assessing miR-138 binding to the BIM-3′UTR region confirmed BIM as a direct target of miR-138. Additionally, control vectors carrying either each at a time or both mutated miR-138-binding sites in the BIM-3′UTR region revealed that the binding site 1223-1249 nt is more potently regulated by miR-138 (Figure (Figure4C;4C; Figure S6A). SiRNA-mediated BIM gene silencing (Figure (Figure4D,4D, Figure S6B) increased TMZ resistance in LN-308 cells (Figure (Figure4E).4E). To confirm the role of the BIM-related apoptotic pathway, we employed sub-lethal concentrations of ABT-737 (Figure S6C), a small-molecule BH3 mimetic that directly inhibits anti-apoptotic proteins Bcl-2, Bcl-xL and Bcl-w [25]. The resistant phenotype induced by BIM gene silencing was indeed reversed by ABT-737 (Figure (Figure4F4F).

Figure 4
BIM downregulation via miR-138 confers TMZ resistance to glioma cells

TMZ treatment does not induce acute apoptosis, but modulates expression of Bcl-2 family proteins

BIM belongs to the BH-3 only group of pro-apoptotic Bcl-2 family members and is a direct regulator of the intrinsic cell death pathway. To assess whether TMZ protection of glioma cells by miR-138 over-expression or by BIM silencing involves cell death regulation, we performed flow cytometry cell death analysis by Annexin V/PI staining at days 1 and 3 after TMZ treatment. The overall amount of cell death induced by TMZ was low (up to 10% on day 3). Still, miR-138 mimic or siBim -transfected cells showed reduced cell death upon TMZ stimulation at day 3 (Figure (Figure5A),5A), suggesting that regulation of apoptosis via miR-138 and BIM is only partially involved in acquired TMZ resistance. In order to assess induction of apoptosis by a complementary approach, we checked for caspase 3 cleavage on protein level in LN-308 cells treated with TMZ, using staurosporin as a positive control. We failed to detect caspase 3 cleavage after a short period of TMZ treatment (1–2 days), as well as after an extended period of 4–8 days (Figure (Figure5B).5B). Next we investigated the protein levels of the pro-apoptotic protein BIM and of the anti-apoptotic Bcl-2 family members Bcl-2, Bcl-xL and Mcl-1 after TMZ treatment. Interestingly, protein levels of Bcl-2 and Mcl-1 increased as early as 1 day after TMZ treatment, and in the case of Bcl-2, the protein levels were further elevated at later time points (day 2–8). BIM levels, however, decreased dramatically already after one day of TMZ treatment and remained constantly down up to 8 days after stimulation (Figure (Figure5B),5B), indicating an induction of a pro-survival signaling pathway mediated by the anti-apoptotic protein stabilization as well as BIM degradation rather than apoptotic cell death.

Figure 5
TMZ and miR-138 induce little apoptosis but cytoprotective autophagy

Autophagy is induced by TMZ and miR-138, negatively correlates with BIM protein levels, and acts as a pro-survival pathway

Different pathways have been reported to promote survival or bypass cell death under stress conditions, including autophagy [2629]. Accordingly, LN-308 were exposed to TMZ and activation of autophagy was measured by conversion of the light chain 3-I (LC3-I) to LC3-II, a marker for autophagy [30]. Indeed, TMZ treatment led to a strong increase in LC3-II levels over time (Figure (Figure5C).5C). Since BIM protein levels are negatively regulated by miR-138, and miR-138 expression in turn is associated with TMZ-induced resistance, we assessed the correlation between miR-138-dependent BIM regulation and autophagy in response to TMZ in glioma cells. Indeed, we recapitulated induction of autophagy in LN-308_R cells, where BIM levels were reduced (Figure (Figure5C).5C). More importantly, knockdown of BIM or transfection of miR-138 mimic led to enhanced autophagy as assessed by increased LC3-II production (Figure (Figure5C).5C). Band intensity quantification confirmed a reduction of BIM protein levels in the miR-138 mimic-transfected cells. To assess whether the induced autophagy in LN-308 cells plays a chemoprotective role, we performed clonogenic survival assays with TMZ in the presence or absence of the autophagy inhibitor 3-methyl adenine (3-MA). When autophagy was blocked, cells were sensitized to TMZ (Figure (Figure5D),5D), suggesting that autophagy is a pro-survival signal mediating TMZ resistance in glioma cells.

DISCUSSION

Increasing evidence suggests that miRNAs play important roles in proliferation, invasion and migration as well as spontaneous and induced cell death in multiple cancer entities, including glioblastoma [31, 32]. Moreover, numerous studies have implicated altered miRNA expression in constitutive or acquired resistance to specific cancer treatments, such as irradiation or chemotherapy [33, 34]. Here we used microarray-based miRNA expression profiling of paired parental and TMZ-resistant glioma cell lines [14] to identify novel miRNA candidates potentially mediating resistance of glioma cells to TMZ, which together with radiotherapy constitutes the current standard of care in glioblastoma patients [3]. We provide evidence that miR-138 mediates acquired TMZ resistance in glioblastoma by regulating the apoptosis inducer BIM and autophagy.

miR-138 has been implicated in the promotion of growth and survival of glioma stem cells [18], however, it appears to have an opposite role as a tumor suppressor in other cancers [16, 35]. Upregulation of miR-138 in recurrent glioblastomas of patients progressing after radiotherapy and TMZ confirmed the potential in vivo relevance of our in vitro profiling approach (Figure (Figure1E).1E). We found no clear-cut association between baseline miR-138 expression levels and TMZ sensitivity of various glioma cell lines (Figure S1). This finding may partially be due to the overwhelming role of MGMT promoter methylation and expression in TMZ sensitivity [36], whose strong impact on alkylator resistance might override the effects of miR-138-mediated TMZ resistance in these cell lines. Moreover, miR-138 may be of particular importance in acquired rather than constitutive TMZ resistance, consistent with the observation that DNA damage induces expression of different microRNAs via transcriptional or post-transcriptional regulation [37].

We confirmed as previously described [18] that overexpression of miR-138 promotes cell cycle progression and drives glioma cell proliferation (Figure 2B–2D), but this is at odds with the report that miR-138 overexpression in glioma reduced cell proliferation in vitro and tumorigenicity in vivo through inducing a G1/S cell cycle arrest by inhibition of a EZH2-CDK4/6-pRb-E2F1 signaling loop [19].

Consistent with our observation in primary tumors in vivo (Figure (Figure1E),1E), transfection of glioma cell lines or glioma stem-like cells with miR-138 mimics increased TMZ resistance (Figure (Figure3A);3A); yet, this effect was only seen in models with low baseline miR-138 expression. Consistently, there was cross-resistance to another alkylating drug, CCNU, but not to irradiation or other chemotherapeutic agents such as vincristine or etoposide, confirming a specific role of this miRNA in the resistance pathway to alkylator-induced DNA damage (Figure (Figure3B;3B; Figure S2B). One challenge that we faced while studying the role of miR-138 was that miR-138 depletion by an inhibitor approach did not show the opposite effect of miR-138 over-expression by adding mimic (Figure S3). We may have encountered this problem due to a technical limitation of the used miRNA silencing approach: miRNA inhibitors are designed as microRNA antisense oligonucleotides, which are supposed to saturate endogenous levels of miRNA by binding them, thus preventing the binding to the target mRNA molecules. However, if the number of antisense oligonucleotides is not in several-fold excess of the miRNA copies per cell, the miRNA saturation is not complete, hence the residual levels that we observe upon miR-138 silencing (Figure S3A) might be sufficient to maintain the original phenotype [38, 39]. We also investigated potential mechanisms that could underlie the induction of TMZ resistance by miR-138. For this purpose, we studied two predicted miR-138 targets for their involvement in TMZ resistance of glioma cells, namely CD166/ALCAM and BCL2L11/BIM. Expression of ALCAM on glioma stem cells has been previously linked to the regulation of glioblastoma invasion [40]. Moreover, RNAi-mediated silencing of ALCAM induced chemotherapy resistance in pancreatic cancer cells [41]. However, while we found decreased ALCAM expression levels in TMZ-resistant cells, our functional analyses did not confirm an effect of ALCAM downregulation on glioma cell sensitivity to TMZ (Figure S4D). In contrast, BIM, a pro-apoptotic BH3-only protein of the Bcl-2 family emerged as an interesting target to study. We succeeded to confirm the in silico predictions that BIM is a direct target of miR-138 by luciferase reporter assays (Figure (Figure4C).4C). We further show that transient knockdown of BIM increased TMZ resistance of glioma cells (Figure (Figure4E),4E), thus mimicking the miR-138 over-expression phenotype. Therefore, we postulate that the resistance phenotype mediated by miR-138 is partly dependent on BIM downregulation. Two other microRNAs, miR-17 and miR-20a, have been previously shown to target BIM, and their suppression sensitized chemoresistant leukemia cells to DNA damage-induced apoptosis [42]. MiR-148a has also been shown to target BIM in glioma stem cells and thereby control cell proliferation and apoptosis [43]. However, the contribution of miR-mediated BIM regulation to the response of glioma cells to TMZ has not been addressed before. We observed that BIM silencing led to increased TMZ resistance, which was reversed by adding sub-lethal concentrations of the BH3-mimetic ABT-737 (Figure (Figure4F).4F). All experiments assessing TMZ resistance employed clonogenic survival assays, which measure the potential of single cells to form colonies over long periods of time (up to 20 days) upon different stimulations. Since BIM regulated resistance to TMZ, we measured TMZ-induced apoptosis by AnnexinV/PI (AV/PI) assay and detection of cleaved caspase 3. Our flow cytometry AV/PI measurements (Figure (Figure5A)5A) revealed a minor induction of apoptosis 72 h after TMZ treatment (up to 10%) as described before [27, 44]. However, these studies did not investigate caspase 3 cleavage after TMZ exposure of glioma cells, but only after alkylating MNNG treatment, reporting cleavage at 96 h post treatment. Yet, LN-308 cells stimulated with TMZ did not exhibit caspase 3 cleavage at early or later time-points (Figure (Figure5B).5B). Literature reports suggest that TMZ can induce autophagy or senescence rather than apoptosis [27, 4547]. These studies offer conflicting conclusions whether autophagy is a pro- or anti-survival mechanism. Indeed, autophagy as a mechanism for maintenance of cellular homeostasis can have tumor-suppressive properties in healthy cells. In neoplastic cells, however, autophagic responses rather protect from intracellular and environmental stress, thus promoting tumor progression and resistance to therapy [29]. We indeed saw induction of autophagy after TMZ treatment at day 1 up to day 8 (Figure (Figure5C).5C). This consequence of TMZ exposure was accompanied by reduction of the pro-apoptotic protein BIM and increase of anti-apoptotic Bcl-2 family proteins (Figure (Figure5B).5B). Interestingly, BIM has previously been described to inhibit autophagy by binding and mislocalizing Beclin 1, an autophagy regulator [48]. Additionally, increased autophagy was seen in TMZ-resistant cells, and upon miR-138 over-expression or BIM silencing to a lesser extent (Figure (Figure5C).5C). Finally, blocking autophagy by the inhibitor 3-MA led to increased TMZ sensitivity (Figure (Figure5D).5D). It is tempting to speculate that TMZ treatment, leading to increased miR-138 expression, which in turn reduces BIM levels, leads to autophagy that results in cell survival, thus TMZ resistance. However, probably additional TMZ-induced mechanisms further regulate autophagy (Figure (Figure66).

Figure 6
Autophagy mediates TMZ resistance via the miR-138/BIM axis

In conclusion, we demonstrate a role for miR-138 in acquired TMZ resistance in glioma cells. Moreover, we establish a functional link between miR-138-mediated resistance and BIM, and confirm that down-regulation of BIM by miR-138 promotes TMZ resistance. Additionally, we propose that the miR-138/BIM axis mediates resistance to TMZ by autophagy. Hence, its targeting may represent a novel strategy to overcome acquired TMZ resistance in glioblastoma.

MATERIALS AND METHODS

Cells and reagents

The human long-term glioma cell lines (LTC) LN-18, LN-229, LN-308, LN-319, LN-428, D247MG, A172, U87MG and T98G were kindly provided by Dr. N. de Tribolet (Lausanne, Switzerland). TMZR sublines (LN-18_R, LN-229_R and LN-308_R) have been characterized [14]. The following glioma-initiating cells (GIC) grown as sphere cultures were used: ZH-161, ZH-305, T-325, T-269 and S-24 [49]. Stock solutions of TMZ (Schering-Plough, Kenilworth, NJ) and lomustine (CCNU) (Dr. Reddy's Laboratories Ltd., Hyderabad, India) at 100 mM were prepared in dimethylsulfoxide (DMSO); aliquots were stored at −20°C. Etoposide (Sigma-Aldrich, St. Louis, MO) stock solution was prepared in DMSO at 50 mM, ABT-737 (ChemieTek, Indianapolis, IN) was prepared as 50 mM stock solution in DMSO. Vincristine was purchased from Tocris Bioscience (Bristol, UK) and prepared as 48 mM stock solution in DMSO. 3-methyl adenine (3-MA) (Sigma solutions) were freshly prepared in water at 100 mM. Irradiation experiments were performed by irradiation of the cells in a Co-radiation source (Gebrüder Sulzer, Thermische Energiesysteme, 60-Co, Winterthur, Switzerland).

MiRNA and mRNA extraction from glioblastoma tissue samples and glioma cells

Deep-frozen tissue samples from primary and recurrent glioblastomas of 10 patients treated with TMZ/RT→TMZ after surgery were selected from the tumor tissue bank of the Department of Neuropathology, Heinrich Heine University, Düsseldorf, Germany, and investigated as approved by the institutional review board (study number 3825). Tissue samples were snap frozen immediately after resection and stored at −80°C. Each tissue specimen was histologically evaluated for sufficient tumor cell content. miRNA was extracted from frozen tissue samples or cell Qiagen, Hilden, Germany). Total RNA was prepared using the NucleoSpin System (Macherey-Nagel AG, Oensingen SO, Switzerland). RNA quantity was determined with the NanoDrop ND-1000 spectrophotometer (PeqLab, Erlangen, Germany). miRNA and mRNA quality assessment was done with an Agilent 2100 bioanalyzer (Agilent Technologies, Santa Clara, CA). Only RNAs with a RNA integrity number of ≥ 7 were used for molecular analyses.

Microarray-based miRNA and mRNA expression profiling

For microarray-based miRNA expression analysis 400 ng of high-molecular weight RNA was biotin-labelled using the Affymetrix® FlashTag™ Biotin HSR RNA Labeling Kit (Affymetrix, Santa Clara, CA) and hybridized to GeneChip® miRNA 2.0 arrays (Affymetrix). Transcriptome-wide changes in gene expression were determined by hybridization to Gene Chip® Human Genome U133 Plus 2.0 arrays (Affymetrix) [50]. Sample preparation for mRNA expression analysis was done with 2.5 μg of total RNA using the One Cycle Target Labelling and Controls kit (Affychip hybridizations were performed at the Center for Biological and Medical Research (BMFZ) at Heinrich Heine University Düsseldorf using a GeneChip Hybridization Oven 645 and a GeneChip® Scanner 3000 7G. Three independent batches from each parental and TMZR cell line were subjected to the Affymetrix chip analyses.

Bioinformatic and statistical analyses of microarray data

All steps of the microarray analysis pipeline were set up within the R framework for statistical computing [51]. The microarray dataset was normalized and log2-transformed using the robust multi-array average (RMA) method implemented in the affy package [52]. Differential expression of miRNA was assessed by computing the moderated t-statistics using the LIMMA package [53]. The resulting p-values were corrected for multiple testing by controlling the false discovery rate [54]. A miRNA was considered to be differentially expressed if the corrected p-value was below 0.05. For in silico predictions of potential miRNA targets, the miRWalk online database on predicted and validated miRNA targets was used [21].

Transfection reagents, miRNA mimics and inhibitors, and siRNA molecules

The miRIDIAN miR-138 mimics and miR-138 miRIDIAN hairpin inhibitors and respective control molecules were purchased from Dharmacon (Dharmacon GE Healthcare, Lafayette, CO). For ALCAM knockdown ON-TARGET plus siRNA-SMART pool (L-004574) was used (Dharmacon). BIM gene silencing was performed using FlexiTube siRNA (siBim1-SI02655359, siBim2 – SI03056466, siBim3 – SI03062080, siBim4 – SI04951961; Qiagen, Hilden, Germany). Transient transfection of glioma cells with mimics, inhibitors or siRNA was performed using Lipofectamine RNAiMAX (Life Technologies, Carlsbad, CA, USA). Lipofectamine 2000 was used for co-transfection of luciferase reporter vector and miRNA mimics or inhibitor molecules (Life Technologies).

Transient transfection

LTC were seeded in six-well plates at 250,000 to 600,000 cells/well. Cells were transfected with specific miRNA mimics, miRNA inhibitors, siRNA, or respective control molecules 24 h later. For clonogenic assays and mRNA or miRNA expression analyses, cells were harvested 24 h post transfection, while for preparation of protein lysates, medium was changed 24 h post transfection and the cells were harvested 48 h later.

GICs were dissociated using Accutase cell detachment solution (PAA Laboratories; Pasching, Austria). Two million cells were transfected using the Neon Transfection System electroporation device and kit (Invitrogen, Life Technologies) with 100 nM mimic, inhibitor, siRNA, or control. After transfection, cells were seeded in 6-well plates in neurobasal conditions for recovery. One day later cells were processed for clonogenic survival assay, mRNA or microRNA isolation, or 72 h later for protein lysate preparation.

Quantitative reverse transcriptase PCR

Complementary DNA transcription from miR-138 and RNU48 and subsequent real-time PCR-FAM detection was performed using specific TaqMan™ MicroRNA Assays (hsa-miR-138 assay ID: #002284; RNU48 assay ID: #001006) (Applied Biosystems, Life Technologies). For mRNA expression analyses, total RNA was transcribed into cDNA using the iScript cDNA Synthesis Kit (Bio-Rad Laboratories; Hercules, CA). cDNA amplification was monitored using SYBRGreen chemistry on the 7300 Real time PCR System (Applied Biosystems). Conditions for PCR reactions were: 40 cycles, 95°C/15 sec, 60°C/1 min, using the following specific primers: GAPDH fwd: 5′-CTCTCTGCTCCTCCTGTTCGAC-3′ (NM_002046, 114-135 nt), GAPDH rv: 5′-TGAGCGATGTGGCTCGGCT-3′ (NM_002046, 182-164 nt), ALCAM fwd: 5′-TTCTTAGCACCTGGCGTTTC-3′ (NM_001627, 293-312 nt), ALCAM rv: 5′-CGGTTCTTTTCGCTGGTATC-3′ (NM_001627, 388-369 nt), BIM fwd: 5′-TTGATTCTTGCAGCCACCCT-3′ (NM_138621, 120-139 nt), BIM rv: 5′-GAACGCAGCGAACCGAATAC-3′ (NM_138621, 193-174 nt). Data analysis was done using the ΔCT method for relative quantification, GAPDH and RNU48 were used as control genes for mRNA or miRNA expression [55].

miRNA expression analyses in glioma tissue specimens

Reverse transcription of mature miR-138 in primary and recurrent glioblastomas was carried out with 100 ng of total RNA as template using stem-loop primers specific for miR-138 and Superscript® II reverse transcriptase (Life Technologies). Real-time PCR on the StepOnePlus™ system (Life Technologies) was done with TaqMan® miRNA assays (Applied Biosystems, Life Technologies) for miR-138 (Assay ID: #002284) and U6 snRNA (Assay ID: #001973) as reference. Expression changes relative to a calibrator sample (Universal Human Reference RNA, Stratagene, La Jolla, CA) were calculated with the 2−ΔΔ Ct method [56].

Immunoblotting

Immunoblot analysis was performed as described previously [57]. The following primary antibodies were used: anti-β-actin (csc-1616) (Santa Cruz Biotechnology, Dallas, TX), anti-SOX4 antibody (ab41891) (Abcam, Cambridge, UK), anti-ALCAM antibody (MAB656) (R & D Systems, Minneapolis, MN) anti-Cleaved Caspase-3 (Asp175) (#9661), anti-Bcl-xL (#2764) and anti-LC3 A/B (#4108) (Cell Signaling Technology, Danvers, MA), and anti-BIM antibody (B7929) (Sigma-Aldrich).

Proliferation, viability, clonogenicity and flow cytometry

BrdU cell proliferation assays were performed using a Cell Signaling Technology kit (Danvers, MA). Briefly, 24 h after miRNA mimic or control transfection, the cells were seeded at 5,000 cells/well in 96-well plates. BrdU reagent was added 24 h later in serum-free DMEM medium, and cells were grown for another 24 h. BrdU detection was performed according to the user's manual. Trypan blue dye exclusion assay was also used to monitor cell viability. Transfected cells were seeded at a density of 5,000 cells/well in 12-well plates. Medium was changed to serum-free DMEM 24 h later for specified experiments; after 24 h floating and attached cells were harvested, and living and dead cells were counted in a Neubauer chamber using 0.4% trypan blue. Clonogenic assays for adherent LTC and GIC cultures were performed as described previously [14, 50]. For cell cycle and cell death analysis, transfected LN-308 cells were treated as indicated, harvested at the indicated time points after transfection (cell cycle analysis) or post-treatment (cell death analysis) and analyzed [14].

3′UTR luciferase gene reporter assay

A 1281-bp fragment of the BIM 3′UTR sequence containing two predicted binding sites for miR-138 was subcloned into the psiCHECK™-2 dual luciferase vector (Promega, Madison, WI). For this purpose cDNA from LN-308 cells was used and the following primers for amplification of the insert: fwd: 5′-CTCGAGATTATGCAGCCAGCGGTTCT-3′, rev: 5′-GCGGCCGCAGCAGGCACAGAGAAAGAGC-3′. Restriction enzymes used for cloning the insert 3′ to the Renilla luciferase were XhoI and NotI (New England Biolabs, Ipswich, MA). For the control vectors the same fragment of BIM 3′UTR sequence was synthesized de novo, bearing mutations in the two miR-138 recognition sites: psiCheck2_mut1 (mutations in site 161-185 nt); psiCheck2_mut2 (mutations in site 1223-1249 nt); psiCheck2_mut3 (mutations in both recognition sites) (Eurofins Genomics, Ebersberg, Germany). Dual luciferase reporter assays were performed seeding 20,000 LN-308 cells per well in a 96-well plate. Following attachment overnight, cells were co-transfected with 100 ng of the respective reporter construct and 50 nM of miR-138 mimic or control molecules. After 48 h, cells were lysed, luciferase activity was measured, and Renilla values normalized to Firefly luciferase.

The cancer genome atlas (TCGA) database

Data was obtained from the glioblastoma data set of the Cancer Genome Atlas network (http://cancergenome.nih.gov/) [22]. The affymetrix probesets used in the TCGA database were 201952_at for ALCAM and 208536_s_at for BCL2L11. Survival analysis within the glioblastoma data set of the TCGA database was performed using the Kaplan-Meier analysis module of the R2 microarray analysis and visualization platform (http://r2.amc.nl).

Statistical analyses

Data are expressed as means and standard deviations. The experiments shown were independently repeated a minimum of three times with similar results. If not otherwise indicated, results of representative experiments are shown. Significance was assessed using two-sided unpaired Student's t-test, paired Student's t-test for comparison of miR-138 expression in primary and recurrent glioblastomas, one-sample Student's t-test for analysis of luciferase reporter results or two-way ANOVA with post hoc Bonferroni test, and correlation was calculated by assessing the Spearman's coefficient (r) (GraphPad Prism 5, La Jolla, CA) (*p < 0.05, **p < 0.01, ***p < 0.001).

SUPPLEMENTARY MATERIALS FIGURES

Acknowledgments

T. Wachtmeister, Center for Biological and Medical Research (BMFZ), Heinrich Heine University Düsseldorf, is greatly acknowledged for expert help with Affymetrix chip hybridization and scanning experiments.

Footnotes

FUNDING

Supported by a grant from the Wilhelm Sander Foundation for Cancer Research to GR and MiW (2012.088.1). GS was supported by MD student grants from the Research Commission of the Medical Faculty and the Research Training Group 1033 at Heinrich Heine University Düsseldorf.

CONFLICTS OF INTEREST

N. Stojcheva, G. Schechtmann, S. Sass, N. Müller, F. Theis, A. Stefanski, K. Stühler, A.-M. Florea and M. Wolter declare no conflicts of interest. C. Happold has received compensation from MSD as a member of the advisory board. G. Reifenberger has received support for research projects from Roche and Merck, and honoraria for lectures or advisory boards from Amgen and Celldex. P. Roth has received honoraria from Novartis and compensation from Roche, MSD and Molecular Partners as a member of advisory boards. M. Weller has received honoraria from Celldex, Isarna Therapeutics, Magforce, MSD Oncology, Novartis, Novocure, Pfizer, Roche Pharma AG, and Teva, and compensation from Acceleron Pharma, Alpinia, Bayer Schering Pharma, Isarna Therapeutics, MSD Oncology, Merck Serono, Novocure, Piqur and Roche/Genentech as a member of advisory boards.

REFERENCES

1. Johnson DR, O'Neill BP. Glioblastoma survival in the United States before and during the temozolomide era. J Neurooncol. 2012;107:359–364. [PubMed]
2. Weller M, van den Bent M, Hopkins K, Tonn JC, Stupp R, Falini A, Cohen-Jonathan-Moyal E, Frappaz D, Henriksson R, Balana C, Chinot O, Ram Z, Reifenberger G, et al. European Association for Neuro-Oncology Task Force on Malignant G EANO guideline for the diagnosis and treatment of anaplastic gliomas and glioblastoma. Lancet Oncol. 2014;15:e395–403. [PubMed]
3. Stupp R, Mason WP, van den Bent MJ, Weller M, Fisher B, Taphoorn MJ, Belanger K, Brandes AA, Marosi C, Bogdahn U, Curschmann J, Janzer RC, Ludwin SK, et al. Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. N Engl J Med. 2005;352:987–996. [PubMed]
4. Shen J, Hung MC. Signaling-Mediated Regulation of MicroRNA Processing. Cancer Res. 2015 [PMC free article] [PubMed]
5. Di Leva G, Garofalo M, Croce CM. MicroRNAs in cancer. Annu Rev Pathol. 2014;9:287–314. [PMC free article] [PubMed]
6. Floyd D, Purow B. Micro-masters of glioblastoma biology and therapy: increasingly recognized roles for microRNAs. Neuro Oncol. 2014;16:622–627. [PMC free article] [PubMed]
7. Garofalo M, Croce CM. MicroRNAs as therapeutic targets in chemoresistance. Drug Resist Updat. 2013;16:47–59. [PMC free article] [PubMed]
8. Ujifuku K, Mitsutake N, Takakura S, Matsuse M, Saenko V, Suzuki K, Hayashi K, Matsuo T, Kamada K, Nagata I, Yamashita S. miR-195, miR-455-3p and miR-10a(*) are implicated in acquired temozolomide resistance in glioblastoma multiforme cells. Cancer Lett. 2010;296:241–248. [PubMed]
9. Munoz JL, Rodriguez-Cruz V, Ramkissoon SH, Ligon KL, Greco SJ, Rameshwar P. Temozolomide resistance in glioblastoma occurs by miRNA-9-targeted PTCH1, independent of sonic hedgehog level. Oncotarget. 2015;6:1190–1201. doi: 10.18632/oncotarget.2778. [PMC free article] [PubMed] [Cross Ref]
10. Shi L, Zhang S, Feng K, Wu F, Wan Y, Wang Z, Zhang J, Wang Y, Yan W, Fu Z, You Y. MicroRNA-125b-2 confers human glioblastoma stem cells resistance to temozolomide through the mitochondrial pathway of apoptosis. Int J Oncol. 2012;40:119–129. [PubMed]
11. Wu H, Liu Q, Cai T, Chen YD, Liao F, Wang ZF. MiR-136 modulates glioma cell sensitivity to temozolomide by targeting astrocyte elevated gene-1. Diagn Pathol. 2014;9:173. [PMC free article] [PubMed]
12. Li P, Lu X, Wang Y, Sun L, Qian C, Yan W, Liu N, You Y, Fu Z. MiR-181b suppresses proliferation of and reduces chemoresistance to temozolomide in U87 glioma stem cells. J Biomed Res. 2010;24:436–443. [PMC free article] [PubMed]
13. Quintavalle C, Mangani D, Roscigno G, Romano G, Diaz-Lagares A, Iaboni M, Donnarumma E, Fiore D, De Marinis P, Soini Y, Esteller M, Condorelli G. MiR-221/222 target the DNA methyltransferase MGMT in glioma cells. PLoS One. 2013;8:e74466. [PMC free article] [PubMed]
14. Happold C, Roth P, Wick W, Schmidt N, Florea AM, Silginer M, Reifenberger G, Weller M. Distinct molecular mechanisms of acquired resistance to temozolomide in glioblastoma cells. J Neurochem. 2012;122:444–455. [PubMed]
15. Girgis AH, Iakovlev VV, Beheshti B, Bayani J, Squire JA, Bui A, Mankaruos M, Youssef Y, Khalil B, Khella H, Pasic M, Yousef GM. Multilevel whole-genome analysis reveals candidate biomarkers in clear cell renal cell carcinoma. Cancer Res. 2012;72:5273–5284. [PubMed]
16. Chakrabarti M, Banik NL, Ray SK. miR-138 overexpression is more powerful than hTERT knockdown to potentiate apigenin for apoptosis in neuroblastoma in vitro and in vivo. Exp Cell Res. 2013;319:1575–1585. [PMC free article] [PubMed]
17. Yeh YM, Chuang CM, Chao KC, Wang LH. MicroRNA-138 suppresses ovarian cancer cell invasion and metastasis by targeting SOX4 and HIF-1alpha. Int J Cancer. 2013;133:867–878. [PubMed]
18. Chan XH, Nama S, Gopal F, Rizk P, Ramasamy S, Sundaram G, Ow GS, Ivshina AV, Tanavde V, Haybaeck J, Kuznetsov V, Sampath P. Targeting glioma stem cells by functional inhibition of a prosurvival oncomiR-138 in malignant gliomas. Cell Rep. 2012;2:591–602. [PubMed]
19. Qiu S, Huang D, Yin D, Li F, Li X, Kung HF, Peng Y. Suppression of tumorigenicity by microRNA-138 through inhibition of EZH2-CDK4/6-pRb-E2F1 signal loop in glioblastoma multiforme. Biochim Biophys Acta. 2013;1832:1697–1707. [PubMed]
20. Strasser A, Cory S, Adams JM. Deciphering the rules of programmed cell death to improve therapy of cancer and other diseases. EMBO J. 2011;30:3667–3683. [PubMed]
21. Dweep H, Sticht C, Pandey P, Gretz N. miRWalk--database: prediction of possible miRNA binding sites by “walking” the genes of three genomes. J Biomed Inform. 2011;44:839–847. [PubMed]
22. Cancer Genome Atlas Research Network Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature. 2008;455:1061–1068. [PMC free article] [PubMed]
23. O'Reilly LA, Cullen L, Visvader J, Lindeman GJ, Print C, Bath ML, Huang DC, Strasser A. The proapoptotic BH3-only protein bim is expressed in hematopoietic, epithelial, neuronal, and germ cells. Am J Pathol. 2000;157:449–461. [PubMed]
24. Chen D, Zhou Q. Caspase cleavage of BimEL triggers a positive feedback amplification of apoptotic signaling. Proc Natl Acad Sci U S A. 2004;101:1235–1240. [PubMed]
25. Oltersdorf T, Elmore SW, Shoemaker AR, Armstrong RC, Augeri DJ, Belli BA, Bruncko M, Deckwerth TL, Dinges J, Hajduk PJ, Joseph MK, Kitada S, Korsmeyer SJ, et al. An inhibitor of Bcl-2 family proteins induces regression of solid tumours. Nature. 2005;435:677–681. [PubMed]
26. Singh K, Matsuyama S, Drazba JA, Almasan A. Autophagy-dependent senescence in response to DNA damage and chronic apoptotic stress. Autophagy. 2012;8:236–251. [PMC free article] [PubMed]
27. Knizhnik AV, Roos WP, Nikolova T, Quiros S, Tomaszowski KH, Christmann M, Kaina B. Survival and death strategies in glioma cells: autophagy, senescence and apoptosis triggered by a single type of temozolomide-induced DNA damage. PLoS One. 2013;8:e55665. [PMC free article] [PubMed]
28. Harder LM, Bunkenborg J, Andersen JS. Inducing autophagy: a comparative phosphoproteomic study of the cellular response to ammonia and rapamycin. Autophagy. 2014;10:339–355. [PubMed]
29. Galluzzi L, Pietrocola F, Bravo-San Pedro JM, Amaravadi RK, Baehrecke EH, Cecconi F, Codogno P, Debnath J, Gewirtz DA, Karantza V, Kimmelman A, Kumar S, Levine B, et al. Autophagy in malignant transformation and cancer progression. EMBO J. 2015;34:856–880. [PubMed]
30. Kabeya Y, Mizushima N, Ueno T, Yamamoto A, Kirisako T, Noda T, Kominami E, Ohsumi Y, Yoshimori T. LC3, a mammalian homologue of yeast Apg8p, is localized in autophagosome membranes after processing. EMBO J. 2000;19:5720–5728. [PubMed]
31. Calin GA, Croce CM. MicroRNA signatures in human cancers. Nat Rev Cancer. 2006;6:857–866. [PubMed]
32. Moller HG, Rasmussen AP, Andersen HH, Johnsen KB, Henriksen M, Duroux M. A systematic review of microRNA in glioblastoma multiforme: micro-modulators in the mesenchymal mode of migration and invasion. Mol Neurobiol. 2013;47:131–144. [PMC free article] [PubMed]
33. Raza U, Zhang JD, Sahin O. MicroRNAs: master regulators of drug resistance, stemness, and metastasis. J Mol Med (Berl) 2014;92:321–336. [PubMed]
34. Low SY, Ho YK, Too HP, Yap CT, Ng WH. MicroRNA as potential modulators in chemoresistant high-grade gliomas. J Clin Neurosci. 2014;21:395–400. [PubMed]
35. Jin Y, Chen D, Cabay RJ, Wang A, Crowe DL, Zhou X. Role of microRNA-138 as a potential tumor suppressor in head and neck squamous cell carcinoma. Int Rev Cell Mol Biol. 2013;303:357–385. [PMC free article] [PubMed]
36. Hermisson M, Klumpp A, Wick W, Wischhusen J, Nagel G, Roos W, Kaina B, Weller M. O6-methylguanine DNA methyltransferase and p53 status predict temozolomide sensitivity in human malignant glioma cells. J Neurochem. 2006;96:766–776. [PubMed]
37. Wan G, Mathur R, Hu X, Zhang X, Lu X. miRNA response to DNA damage. Trends Biochem Sci. 2011;36:478–484. [PMC free article] [PubMed]
38. Gentner B, Schira G, Giustacchini A, Amendola M, Brown BD, Ponzoni M, Naldini L. Stable knockdown of microRNA in vivo by lentiviral vectors. Nat Methods. 2009;6:63–66. [PubMed]
39. Brown BD, Naldini L. Exploiting and antagonizing microRNA regulation for therapeutic and experimental applications. Nat Rev Genet. 2009;10:578–585. [PubMed]
40. Kijima N, Hosen N, Kagawa N, Hashimoto N, Nakano A, Fujimoto Y, Kinoshita M, Sugiyama H, Yoshimine T. CD166/activated leukocyte cell adhesion molecule is expressed on glioblastoma progenitor cells and involved in the regulation of tumor cell invasion. Neuro Oncol. 2012;14:1254–1264. [PMC free article] [PubMed]
41. Hong X, Michalski CW, Kong B, Zhang W, Raggi MC, Sauliunaite D, De Oliveira T, Friess H, Kleeff J. ALCAM is associated with chemoresistance and tumor cell adhesion in pancreatic cancer. J Surg Oncol. 2010;101:564–569. [PubMed]
42. Weng H, Huang H, Dong B, Zhao P, Zhou H, Qu L. Inhibition of miR-17 and miR-20a by oridonin triggers apoptosis and reverses chemoresistance by derepressing BIM-S. Cancer Res. 2014;74:4409–4419. [PubMed]
43. Kim J, Zhang Y, Skalski M, Hayes J, Kefas B, Schiff D, Purow B, Parsons S, Lawler S, Abounader R. microRNA-148a is a prognostic oncomiR that targets MIG6 and BIM to regulate EGFR and apoptosis in glioblastoma. Cancer Res. 2014;74:1541–1553. [PMC free article] [PubMed]
44. Roos WP, Batista LF, Naumann SC, Wick W, Weller M, Menck CF, Kaina B. Apoptosis in malignant glioma cells triggered by the temozolomide-induced DNA lesion O6-methylguanine. Oncogene. 2007;26:186–197. [PubMed]
45. Kanzawa T, Germano IM, Komata T, Ito H, Kondo Y, Kondo S. Role of autophagy in temozolomide-induced cytotoxicity for malignant glioma cells. Cell Death Differ. 2004;11:448–457. [PubMed]
46. Natsumeda M, Aoki H, Miyahara H, Yajima N, Uzuka T, Toyoshima Y, Kakita A, Takahashi H, Fujii Y. Induction of autophagy in temozolomide treated malignant gliomas. Neuropathology. 2011;31:486–493. [PubMed]
47. Lin CJ, Lee CC, Shih YL, Lin CH, Wang SH, Chen TH, Shih CM. Inhibition of mitochondria- and endoplasmic reticulum stress-mediated autophagy augments temozolomide-induced apoptosis in glioma cells. PLoS One. 2012;7:e38706. [PMC free article] [PubMed]
48. Luo S, Garcia-Arencibia M, Zhao R, Puri C, Toh PP, Sadiq O, Rubinsztein DC. Bim inhibits autophagy by recruiting Beclin 1 to microtubules. Mol Cell. 2012;47:359–370. [PMC free article] [PubMed]
49. Lemke D, Weiler M, Blaes J, Wiestler B, Jestaedt L, Klein AC, Low S, Eisele G, Radlwimmer B, Capper D, Schmieder K, Mittelbronn M, Combs SE, et al. Primary glioblastoma cultures: can profiling of stem cell markers predict radiotherapy sensitivity? J Neurochem. 2014 [PubMed]
50. Happold C, Roth P, Silginer M, Florea AM, Lamszus K, Frei K, Deenen R, Reifenberger G, Weller M. Interferon-beta induces loss of spherogenicity and overcomes therapy resistance of glioblastoma stem cells. Mol Cancer Ther. 2014;13:948–961. [PubMed]
51. R Core Team . R: A language and environment for statistical computing. R Foundation for Statistical Computing; Vienna, Austria: 2015.
52. Gautier L, Cope L, Bolstad BM, Irizarry RA. affy—analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004;20:307–315. [PubMed]
53. Smyth GK, Michaud J, Scott HS. Use of within-array replicate spots for assessing differential expression in microarray experiments. Bioinformatics. 2005;21:2067–2075. [PubMed]
54. Benjamini Y, Hochberg Y. Controlling the False Discovery Rate-a Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B-Methodological. 1995;57:289–300.
55. Codo P, Weller M, Meister G, Szabo E, Steinle A, Wolter M, Reifenberger G, Roth P. MicroRNA-mediated down-regulation of NKG2D ligands contributes to glioma immune escape. Oncotarget. 2014;5:7651–7662. doi: 10.18632/oncotarget.2287. [PMC free article] [PubMed] [Cross Ref]
56. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−ΔΔC(T)) Method. Methods. 2001;25:402–408. [PubMed]
57. Silginer M, Weller M, Ziegler U, Roth P. Integrin inhibition promotes atypical anoikis in glioma cells. Cell Death Dis. 2014;5:e1012. [PMC free article] [PubMed]

Articles from Oncotarget are provided here courtesy of Impact Journals, LLC