- Research
- Open Access
- Published:
Glutaminase isoforms expression switches microRNA levels and oxidative status in glioblastoma cells
Journal of Biomedical Science volume 28, Article number: 14 (2021)
Abstract
Background
Glutaminase isoenzymes GLS and GLS2 play apparently opposing roles in cancer: GLS acts as an oncoprotein, while GLS2 (GAB isoform) has context specific tumour suppressive activity. Some microRNAs (miRNAs) have been implicated in progression of tumours, including gliomas. The aim was to investigate the effect of GLS and GAB expression on both miRNAs and oxidative status in glioblastoma cells.
Methods
Microarray profiling of miRNA was performed in GLS-silenced LN229 and GAB-transfected T98G human glioblastoma cells and their wild-type counterparts. Results were validated by real-time quantitative RT-PCR. Oxidative status and antioxidant enzymes were determined by spectrophotometric or fluorescence assays in GLS-silenced LN229 and T98G, and GAB-transfected LN229 and T98G.
Results
MiRNA-146a-5p, miRNA-140-3p, miRNA-21-5p, miRNA-1260a, and miRNA-92a-3p were downregulated, and miRNA-1246 was upregulated when GLS was knocked down. MiRNA-140-3p, miRNA-1246, miRNA-1260a, miRNA-21-5p, and miRNA-146a-5p were upregulated when GAB was overexpressed. Oxidative status (lipid peroxidation, protein carbonylation, total antioxidant capacity, and glutathione levels), as well as antioxidant enzymes (catalase, superoxide dismutase, and glutathione reductase) of silenced GLS glioblastoma cells and overexpressed GAB glioblastoma cells significantly changed versus their respective control glioblastoma cells. MiRNA-1246, miRNA-1260a, miRNA-146a-5p, and miRNA-21-5p have been characterized as strong biomarkers of glioblastoma proliferation linked to both GLS silencing and GAB overexpression. Total glutathione is a reliable biomarker of glioblastoma oxidative status steadily associated to both GLS silencing and GAB overexpression.
Conclusions
Glutaminase isoenzymes are related to the expression of some miRNAs and may contribute to either tumour progression or suppression through certain miRNA-mediated pathways, proving to be a key tool to switch cancer proliferation and redox status leading to a less malignant phenotype. Accordingly, GLS and GAB expression are especially involved in glutathione-dependent antioxidant defence.
Background
Although metabolic profiles of tumours depend on both the genotype and tissue of origin, cancer cells have a very different metabolic profile compared to normal cells [1]. Glutaminase (GA; EC 3.5.1.2) catalyses the first step of glutaminolysis transforming glutamine into glutamate. Glutamate is one of the precursor amino acids for glutathione, which is present within the cell in a reduced (GSH) and an oxidized form (GSSG) [2]. The ratio of GSH/GSSG is the most important regulator of the cellular redox status, which together with metabolic stress switched by glutamine and GA can be a key regulatory tool to determine cell fate, leading to normal or tumour proliferation [3]. In humans, the GA family consists of two genes, GLS and GLS2, located in chromosome 2 and 12, respectively [2]. The GLS gene encodes two isoforms, known as kidney (K-type) glutaminase or KGA, and a shorter spliced form named glutaminase C or GAC [4]. These two isoenzymes are usually referred to as GLS [3]. On the other hand, the GLS2 gene codes for the liver (L-type) isozymes, named LGA, as well as for a longer isoform named GAB originally described in breast cancer cells [2]. Both isoenzymes are collectively designated as GLS2 [3].
MiRNAs are a family of 19- to 25-nucleotides noncoding small RNAs that primarily function as gene regulators suppressing mRNA translation. Aberrant miRNA expression has been described for several human malignancies, and this new class of small regulatory RNAs can have both oncogenic and tumour suppressor functions [5]. MiRNAs may interact with the 3’-UTRs of target mRNA molecules through partial sequence homology. They form an important class of regulators that participate in diverse biological functions including cell proliferation, differentiation and apoptosis. In fact, recent studies have shown that some miRNAs regulate the expression of several tumour suppressors and oncogenes and their dysregulation contributes to carcinogenesis [6]. It is now evident that changes in miRNA are involved in cancer progression, but the mechanisms of transcriptional regulation of miRNAs remain unknown [5]. To test if GLS isoenzyme expression is essential for miRNA regulation in glioblastoma (GBM) cells, an independent GLS-targeted shRNA was expressed in LN229 cells. We chose this cell line because it uses glutamine as the preferred anaplerotic precursor when given access to a physiological mixture of nutrients [7]. Furthermore, this cell line does not express appreciable amounts of the GLS2 isoforms LGA or GAB, increasing the likelihood of reducing cell growth by targeting GLS [7]. Additionally, stable transfection of GBM T98G cells with a vector carrying human full-length GAB sequence, T98G-GAB(+) [3], was carried out to ascertain the role of GAB in GBM cells, which has been linked to p53 pathway and antioxidant function [2, 3]. Hence, GLS2 increased cellular levels of GSH and NADH and decreases reactive oxygen species (ROS) levels in hepatocarcinoma cells [8]. Consistently, GLS2 protected these cells from DNA oxidation and ROS-sensitive apoptosis [9]. On the other hand, miRNA-200c induced ROS generation in ischemic cardiomyocytes through GA as GLS has been characterized as a direct target of miR-200c in this cell model [10]. Interestingly, long non-coding ribonucleic acid (lncRNA) urothelial carcinoma-associated 1 (UCA1) was highly expressed in bladder cancer cells. UCA1 regulates the expression of GLS2 by interfering with miR-16, and blocked ROS formation in bladder cancer [11]. In this study, we aim to elucidate whether some miRNAs can be modulated either by inhibition of GLS or by overexpression of GAB, impacting the redox state of cancer cells and contributing to the alterations in the markers of oxidative stress.
Methods
Cell lines, culture conditions and stable transfections
GLS-silenced LN229 and its control LN229 cells were established and characterized as described [7]. Both GLS isoforms KGA and GAC have been silenced as previously published by authors [3, 7]. Those cells were cultured in Dulbecco’s Modified Eagle Medium (DMEM) supplemented with 10 % foetal bovine serum (FBS), 100 I.U./mL penicillin, 100 µg/mL streptomycin, and 6 mM L-glutamine as previously described by authors [7]. In all stable knockdown experiments, very few detached cells were noted in the culture, and these were not included in growth and viability counts. T98G human GBM cell line was purchased from American Type Culture Collection (ATCC, Rockville, MD, USA). T98G-GAB(+) and control cells T98G-pcDNA were maintained in minimum essential medium supplemented with 10 % FBS, 1 % non-essential amino acids, and 100 I.U./mL penicillin and 100 µg/mL streptomycin. LN229-GAB(+) and control LN229-pcDNA human GBM cell lines were kindly provided by Monika Szeliga, Department of Neurotoxicology, Mossakowski Medical Research Centre, Polish Academy of Sciences, Warsaw, Poland. LN229-GAB(+) and LN229-pcDNA cell lines were maintained in DMEM as previously described for other LN229 cell lines [7]. GLS-silenced T98G and its control T98G were generated by transient silencing, using GLS-targeted small interferent RNAs (siRNAs) and a negative control siRNA transfection, respectively. Cells were cultured in DMEM as previously described for other T98G models. All cultures were maintained at 37 °C in a humidified atmosphere with 95 % air and 5 % CO2. T98G-GAB(+), T98G-pcDNA, LN229-GAB(+) and LN229-pcDNA cell lines were obtained by stable transfection of T98G or LN229 cells with the full-length cDNA sequence encoding human GAB isoform or empty pcDNA3 vector, respectively, as described previously [3, 12]. The culture medium for the polyclonal populations of T98G-GAB(+) and its control T98G-pcDNA, and for LN229-GAB(+) and its control LN229-pcDNA cells containing the neomycin-resistance gene was supplemented with 0.5 and 0.75 mg/mL geneticin, respectively. All cell lines were authenticated by both microscopic cellular morphology and growth curve analysis. All cell lines were also checked to be free of mycoplasma contamination.
Protein expression
KGA and GAC downregulation was confirmed in the newly prepared cell line, T98G-GLS(−) by Western-blot (Additional file 1: Fig. S1). Cells were washed in PBS and harvested in lysis buffer (1 % NP-40, 150 mM NaCl, 50 mM Tris-HCl pH 8.0, supplemented with cOmplete™ Mini Protease Inhibitor, Roche, Indianapolis, IN, USA). Lysates were centrifuged at 12,000 g, 15 min at 4 °C. Supernatants were collected and stored at -80 °C. Protein quantification was performed using Pierce™ BCA Protein Assay Kit (ThermoFisher Scientific, Waltham, MA, USA). After heating 5 min at 95 °C, twenty micrograms of protein was resolved on 10 % polyacrylamide-SDS gels and then electrotransferred to nitrocellulose membranes. Membranes were blocked for 1 h at RT with 5 % non-fat milk in TBST. For GLS detection, two primary antibodies targeting specifically KGA or GAC isoforms were used, diluted in 5 % BSA in TBST and incubated overnight at 4 °C in agitation. These two GLS-targeted antibodies were kindly donated by Prof. Norman P. Curthoys (Colorado State University, CO, USA), both were generated in Rabbit again peptides containing the 20 C-terminal aminoacids of human KGA or GAC, respectively. Both antibodies were demonstrated to react specifically with its GLS isoform (KGA or GAC), and not to react with other glutaminase isoforms [3]. After incubation with respective primary antibodies, the membranes were washed in TBST and incubated 1 h RT with a commercial secondary HRP-conjugated Goat anti rabbit antibody (A0545, Merck) in 5 % BSA in TBST, and visualized in a ChemiDoc™ Gel Imaging System (Bio-Rad, Hercules, CA, USA), using an enhanced chemiluminiscence (ECL) detection system (SuperSignal™ West Pico, ThermoFisher Scientific). A commercial antibody targeting β-Actin (MAB1501, Merck) was used as a loading control, following the same procedure described before for the GLS antibodies, but using a Rabbit anti Mouse secondary HRP-conjugated antibody (A9044, Merck) in 5 % BSA in TBST (Additional file 1: Fig. S1a, c). For data quantification, Image Lab Software (Bio-Rad) was used, following the manufacturer’s instructions. In brief, each band was detected and selected, and background signal was descarted. After that, net values corresponding to the signal of each KGA or GAC band was normalized to β-Actin signal. For a more clear representation, the values obtained for T98G-GLS(−) and its T98G control were referred to the T98G control values, and the resulting ratio was represented (Additional file 1 : Fig. S1b, d).
RNA interference and miRNA array experiments
All RNA interference (RNAi) experiments used pools of cells. Vectors for RNAi, lentiviral particles and details to obtain GLS knockdown LN229 cells have been described previously by authors [7]. Stably infected pools with adequate silencing were maintained with 1 µg/mL puromycin. LN229-GLS (−) express an shRNA targeting GLS, while its control LN229 cells contain the same construct but expressing green fluorescent protein (GFP) (LN229-GFP) instead of the GLS-targeted shRNA. GLS silencing in T98G cells were carried out by transfection using Lipofectamine™ 3000 reagent (Invitrogen, Grand Island, NY, USA) with GLS-targeted siRNAs, using a 1:1 mixture of two commercial siRNAs targeting GLS (NM_014905, sequence start around 1072 and 780, Merck KGaA, Darmstadt, Germany).
For miRNA expression analysis of different GBM cell lines, RNA was extracted using standard RNA extraction protocols (Trizol™, Invitrogen, Grand Island, NY, USA) according to manufacturer’s instructions. Four vials containing the GBM cell samples were shipped to Miltenyi Biotec, Bergisch Gladbach (Germany), on dry ice. RNA quality was assessed with the Agilent 2100 Bioanalyser (Agilent Technologies, Palo Alto, CA). Samples included in the present analysis had RNA Integrity Numbers (RIN) between 9.8 and 10.0 (Additional file 2: Fig. S2a, b). Labelling and hybridization were performed according to user manuals of the miRXplore™ instrument (Miltenyi Biotec). In brief, 1.2 µg/sample total RNA were labelled with the red fluorescent Hy5 using the miRNA/LNA labelling Exiqon kit. A pool of synthetic miRNAs in equimolar concentrations was designed by Miltenyi based on sequences of miRBase 9.2 and labelled with Hy3 (Additional file 3: Table S1). Subsequently, the labelled material was hybridized overnight to miRXplore™ microarrays using the a-Hyb™ Hybridization Station (Miltenyi Biotec). Fluorescence signals of the hybridized miRXplore™ microarrays were detected using ImaGene® software (Biodiscovery, Hawthorne, CA, USA) and raw data were acquired with PIQOR™ analyser software (Miltenyi Biotec), as detailed in Additional file 2: Fig. S2c, d.
RNA isolation and real‐time quantitative reverse transcription-PCR
Total RNA was isolated from cultured cells, using miRCURY RNA isolation kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Twelve miRNAs were selected from the miRNA microarray analysis (see Table 1). Only the RNA samples with 260/280 ratios of 1.8–2.0 were used for further investigation. cDNA was generated from 1 µg total RNA using the miRCURY LNA Universal RT microRNA PCR, Universal cDNA Synthesis Kit II (Qiagen). Real-time quantitative reverse transcriptase-polymerase chain reaction (real-time qRT-PCR) was performed using the miRCURY LNA Universal RT microRNA PCR, ExiLENT SYBR® Green PCR Kit (Qiagen) according to the manufacturer’s instructions, using CFX96 Real-Time System C1000 Thermal Cycler (Bio-Rad). Oligonucleotides sequences are shown in Table 1. Samples that lacked either a template or reverse transcriptase were used as controls. To normalize the data for quantification of miRNAs, U6 small nuclear RNA was selected as the best endogenous control. Control LN229 and T98G cell lines were used as the miRNA expression level references. Real-time qRT-PCR amplifications for the miRNAs and U6 snRNA were performed after an initial denaturation step of 10 min at 95 °C followed by 40 cycles of 10-sec denaturation at 95 °C and 1-min annealing/extension at 60 °C. The relative expression of each miRNA was calculated according to the comparative 2−ΔΔCt quantification method where ΔCt=Ct(sample) − Ct(normaliser) and ΔΔCt = ΔCt(sample)−ΔCt(calibrator). The control U6 reactions were implemented with 5′-CTCGCTTCGGCAGCACA-3′ as the forward primer and 5′-AACGCTTCACGAATTTGCGT-3′ as the reverse primer. Standard deviation (SD) was incorporated into fold change by expressing the fold change as a range. This is done by incorporating the SD into the formula: 2 − ΔΔCt + SD and 2 − ΔΔCt – SD.
GSH assays
The level of total GSH in every cell line was determined by using glutathione assay kit obtained from Cell Biolabs INC, San Diego, California, USA (#STA-312-T), according to the manufacturer’s instructions. Briefly, approximately 3.5 × 106 cells were resolved. Supernatant samples were transferred to a microplate and added 25 µL of glutathione reductase (GR), 25 µL of NADPH solution, 100 µL of cell lysate, and 50 µL of assay buffer with 5,5′-dithiobis-(2-nitrobenzoic acid) (DTNB), reading every 60 s, during 10 min, at 405 nm. Absorbance was determined using a microplate spectrophotometer Eon, Biotek Instruments, Inc. (Winooski, VT, USA). All samples were run in duplicate. Total GSH amounts were normalized by the protein content (Bradford) in each sample, and a ratio was calculated by dividing the experimental values by the control values for each independent experiment.
GR activity was measured using the Glutathione Reductase Assay kit from Cell Biolabs, INC. (#STA-812), as described by suppliers. In brief, approximately 3.5 × 106 cells were used for each experiment. Then, 25 µL of NADPH solution was added to a well in the 96 well plate and then 100 µL of cell lysate, 25 µL of GSSG, and 50 µL of assay buffer containing DTNB were added to each well. The absorbance was read at 405 nm every 1 min for 10 min in a microplate spectrophotometer Eon, Biotek Instruments, Inc. to calculate the activity. All samples were run in duplicate. GR activities were normalized by the protein content (Bradford) in each sample, and a ratio was calculated by dividing the experimental values by the control values for each independent experiment.
Measurement of antioxidant activities
Catalase (CAT) activity was measured using the Catalase Assay kit from Cell Biolabs, INC (#STA-341-T), according to the manufacturer’s instructions. In brief, approximately 3.5 × 106 cells were used for each experiment. Then, 20 µL of cytosolic lysate was mixed with 50 µL of 12 mM H2O2 solution and incubated for 1 min at room temperature. The reaction was stopped by adding 50 µL of stop solution (containing sodium azide). Then, 5 µL out of the 120 µL reaction mixture was mixed with 250 µL of the colour reagent in a new microplate. After 50 min of incubation for colour development, the absorbance was measured at 520 nm in a microplate spectrophotometer Eon, Biotek Instruments, Inc. All samples were run in duplicate. CAT activities were normalized by the protein content (Bradford) in each sample, and a ratio was calculated by dividing the experimental values by the control values for each independent experiment.
Superoxide dismutase (SOD) activity was measured using a SOD assay kit (Cell Biolabs, INC, #STA-340-T), as described by suppliers. Approximately 3.5 × 106 cells were used for each experiment. SOD activity was measured by assessing the degree of inhibition of the reduction of nitroblue tetrazolium by xanthine oxidase. One unit of activity was considered as the amount of protein that gives 50 % inhibition of the reaction. Briefly, the cell lysates were incubated in a 96-well plate and the absorbance at 490 nm was measured kinetically at 37 °C for 30 min. Absorbance was determined using a microplate spectrophotometer Eon, Biotek Instruments, Inc. All samples were run in duplicate. SOD activities were normalized by the protein content (Bradford) in each sample, and a ratio was calculated by dividing the experimental values by the control values for each independent experiment.
Oxidative stress assessment
Total antioxidant capacity (TAC) was determined using Cell Biolabs, INC kit (#STA-360-T), evaluating cupric reducing antioxidant capacity, based on the complex equilibria between Cu(II)–bathocuproinedisulfonic acid and Cu(I)–bathocuproinedisulfonic acid. Briefly, 3.5 × 106 cells were used for each experiment. Then, 20 µL of cell lysate was loaded into a well in the 96 well plate and added 180 µL of reaction mix. Absorbance was quantified at 490 nm before adding 50 µL of Cu(II) solution to each well. Samples were incubated 5 min at RT; then, redox reactions were halted with 50 µL of a stop solution and final absorbance was measured. Absorbance was determined using a microplate spectrophotometer Eon, Biotek Instruments, Inc. All samples were run in duplicate. TAC values were calculated as nmol Copper Reducing Equivalents (CRE) and normalized by the protein content (Bradford) in each sample, and then a ratio was calculated by dividing the experimental values by the control values for each independent experiment.
Lipid peroxidation was evaluated by thiobarbituric acid reactive substances (TBARS) method, using Cell Biolabs, INC kit (#STA-330-T), which is based on malondialdehyde (MDA) reaction with TBARS, producing a pink complex with a peak absorbance at 532 nm. To ensure that no lipid oxidation occurs during the assay, butylated hydroxytoluene was added to the sample prior to lysis. Approximately 3.5 × 106 cells were used for each experiment. Cell lysates were used to measure lipid peroxidation level according to the manufacturer’s instruction. Briefly, 100 µL of cell lysate were treated with reagents during 1 h at 95 °C, cooled on ice for 5 min, centrifuged at 1000 g for 15 min and used for measurement in a 96-well plate. Absorbance was determined at 532 nm using a microplate spectrophotometer Eon, Biotek Instruments, Inc. All samples were run in duplicate. TBARS values were normalized by the protein content (Bradford) in each sample and a ratio was calculated by dividing the experimental values by the control values for each independent experiment.
Protein carbonylation was performed as described by Cell Biolabs, INC kit (#STA-307-T) using a fluorescein 5-thiosemicarbazide (FTC)-based method. Approximately 7 × 106 cells were used for each experiment. In brief, 50 µL of cell lysate was mixed with 50 µL of FCT solution. Proteins were precipitated by the addition of 4-volumes of cold 20 % trichloroacetic acid. Following 10 min incubation on ice, tubes were centrifuged at 9000 g for 10 min at 4 °C. Supernatants were decanted; precipitates were washed 3 times by vortexing with 1 mL of acetone. Finally, acetone supernatant was carefully decanted out and protein precipitates were air dried, solubilized with 50 µL of 6 M guanidine hydrochloride and diluted by the addition of 450 µL of pH 7.0 buffer. Protein concentration in each of these samples was measured by Bradford assay. The samples were aliquoted 100 µL per well in triplicate, and fluorescence measured in a microplate fluorescence reader FL600, Biotek Instruments, Inc., with excitation at 480 nm and emission at 530 nm. Fluorescence readings from 6 wells for each sample (3 for each of the duplicates) were averaged and nanomoles of FTC-reacted carbonyls were calculated using a standard curve generated from the readings of various concentrations of FTC prepared in a medium similar to that of the samples. Values were normalized to the protein content and expressed as pmol carbonyl/mg protein, and finally a ratio was calculated by dividing the experimental values by the control values for each independent experiment.
Statistical analysis
For array and real-time qRT–PCR experiments, miRNAs with > 10 % undetectable expression values were excluded from the study. Statistical analysis was performed using software Statgraphics. Values were expressed as means ± SEM. Differences between groups were calculated with the Student’s t-test or Mann–Whitney U test. In addition, Spearman’s rank-order correlation test was performed. All results were considered significantly as p < 0.05. Experiments were performed at least three times.
Results
Up- and downregulation of miRNAs following GLS silencing or GAB overexpression
The miRXplore™ chips enabled examination of 993 human miRNAs (Additional file 4: Table S2). To investigate the role of GLS and GAB isoforms in miRNAs regulation, we first selected the 23 miRNAs that strongly changed their expression levels when GLS was silenced or GAB was overexpressed (Table 2). Then, we analysed these profiles to find eight miRNAs (miRNA-21-5p, miRNA-92a-3p, miRNA-140-3p, miRNA-146a-5p, miRNA-762, miRNA-1246, miRNA-1260a, and miRNA-1274) which were highly under- or overexpressed in both models (GLS silenced cells or cells overexpressing GAB). In LN229-GLS(−) vs. LN229-GFP cells, thirteen miRNAs demonstrated enhanced overexpression, highlighting miRNA-1246 that showed 2.83-fold augmented expression. Among five downregulated miRNAs, miRNA-146a showed the highest reduced expression: 0.21-fold in these GLS silenced cells when compared to control LN229-GFP cells (Table 2). In T98G-GAB(+) compared to T98G-pcDNA control cells, 13 miRNAs showed higher expression, including miRNA-140-3p, miRNA-92a-3p, miRNA-1260, miRNA-21, miRNA-1246, and miRNA-146a-5p (Table 2). Only miRNA-762, miRNA-1246, and miRNA-1274 were similarly expressed in both models. The results from the microarray analysis were validated by real-time quantitative RT-PCR, as described in "Methods" section. From the twelve selected miRNAs (Table 1), six were analysed in two GBM cell models: LN229 and T98G cells. Both models were studied when GLS was silenced and when GAB was overexpressed (Fig. 1): LN229-GLS(−) cells vs. control LN229-GFP, T98G-GLS(−) cells vs. control (scramble siRNA-transfected) T98G, LN229-GAB(+) cells vs. control LN229-pcDNA, and T98G-GAB(+) cells vs. control T98G-pcDNA. Five miRNAs (miRNA-146a-5p, miRNA-140-3p, miRNA-21-5p, miRNA-1260a, and miRNA-92a-3p) were significantly repressed in LN229-GLS(−) vs. LN229-GFP cells, whilst miR-1246 was 2.3-fold upregulated in LN229-GLS(−) vs. LN229-GFP cells (Fig. 1a). Similar results were found comparing T98G-GLS(−) vs. control T98G cells, but miR-1246 was 1.5-fold upregulated in T98G-GLS(−) vs. control T98G cells (Fig. 1b). All six miRNAs assayed were overexpressed in LN229-GAB(+) cells (miRNA-21-5p, miRNA-92a-3p, miRNA-140-3p, miRNA-146a-5p, miRNA-1246 and miRNA-1260a) as stated in Fig. 1c. Five miRNAs: miRNA-140-3p, miRNA-1246, miRNA-1260a, miRNA-21-5p, and miRNA-146a-5p, were all significantly overexpressed in T98G-GAB(+) vs. control T98G-pcDNA cells. No significant differences were found for miRNA-92a-3p in this cell model (Fig. 1d).
Down- and up-regulated miRNAs. Results are represented for LN229 (LN229-GFP) and LN229-GLS(−) cells (a), T98G (scrambled-transfected T98G) and T98G-GLS(−) (b), LN229 (LN229-pcDNA) and LN229-GAB(+) (c) and for T98G (T98G-pcDNA) and T98G-GAB(+) cells (d). The relative expression of each miRNA was calculated according to ΔΔCt method. Vertical coordinates represent the average fold change
As stated in "Background" section, miRNAs over- or under-expression has been associated with GLS [10] or GLS2 [11] expression levels, which are related with antioxidant capacity [2, 8, 9]. Following, we have measured a number of oxidative stress biomarkers in our GBM cell models that up- or down-express GA isoenzymes.
Glutathione-dependent antioxidant capacity is an associated-biomarker to GLS and GAB expression
When GLS was silenced total glutathione levels were reduced by 45 % in LN229-GLS(−) vs. LN229-GFP cells. Simultaneously, GR activity was increased by 50 % (Fig. 2a). Almost identical results were found comparing T98G-GLS(−) vs. T98G cells (Fig. 2b). When GAB was overexpressed in LN229 cells, glutathione levels were augmented 1.5-fold and GR activity was also significantly increased (Fig. 2c). Similar significant results were found for both GSH-dependent oxidative biomarkers when GAB was overexpressed in T98G cells (Fig. 2d).
Total glutathione levels and GR activity. Results are represented for LN229 (LN229-GFP) and LN229-GLS(−) cells (a), T98G (scrambled-transfected T98G) and T98G-GLS(−) (b), LN229 (LN229-pcDNA) and LN229-GAB(+) (c) and for T98G (T98G-pcDNA) and T98G-GAB(+) cells (d). The results were calculated as nmoles GSH/mg protein, and mU/mg protein, respectively, and a ratio was calculated by dividing the values of the experimental sample by those of its control for each biological replicate. Data are expressed as mean ± S.E.M. (N = 3) of the ratios calculated for at least three independent experiments. *p < 0.05, **p < 0.01, ***p < 0.001, compared to the respective control group. GR glutathione reductase, tGSH total glutathione
Catalase and superoxide dismutase activities change following both GLS silencing and GAB overexpression
When GLS was silenced both CAT and SOD activities were significantly reduced in LN229-GLS(−) vs. LN229-GFP cells (Fig. 3a). However, CAT and SOD activities were both significantly enhanced when GLS was silenced in T98G cells (Fig. 3b). Decreased CAT and SOD activities were found when GAB was overexpressed in LN229 cells (Fig. 3c). SOD activity was significantly increased when GAB was overexpressed in T98G cells; however, CAT activity was slightly reduced (Fig. 3d).
Antioxidant enzymatic activities. SOD and CAT were determined, for LN229 (LN229-GFP) and LN229-GLS(−) cells (a), T98G (scrambled-transfected T98G) and T98G-GLS(−) (b), LN229 (LN229-pcDNA) and LN229-GAB(+) (c) and for T98G (T98G-pcDNA) and T98G-GAB(+) cells (d). The results are expressed as U/mg protein and a ratio was calculated by dividing the values of the experimental sample by those of its control for each biological replicate. Data are expressed as mean ± S.E.M. (N = 3) of the ratios calculated for at least three independent experiments. **p < 0.01, ***p < 0.001, compared to the respective control group. CAT catalase, SOD superoxide dismutase
Lipid peroxidation was attenuated by both GLS silencing and GAB overexpression. TAC decreased by GLS silencing and increased when GAB was overexpressed
Lipid peroxidation fell 10-fold in LN229-GLS(−) vs. LN229-GFP cells. In this model, TAC decreased to two thirds after GLS silencing (Fig. 4a). TAC trend was identical after GLS silencing in T98G model. Lipid peroxidation dropped almost to half in in T98G-GLS(−) vs. control T98G cells. (Fig. 4b). TAC was significantly increased and lipid oxidative damage, measured as TBARS, was significantly decreased when GAB was overexpressed in LN229 cells (Fig. 4c). Similar results were found when GAB was overexpressed in T98G cells (Fig. 4d). In addition, when GAB was overexpressed in T98G cells, protein damage declined to half. This oxidative stress parameter was 25 % reduced when GLS was silenced in LN229 cells (Fig. 4e).
Oxidative status values. Lipid peroxidation and TAC are represented for LN229 (LN229-GFP) and LN229-GLS(−) (a), T98G (scrambled-transfected T98G) and T98G-GLS(−) (b), LN229 (LN229-pcDNA) and LN229-GAB(+) (c) and for T98G (T98G-pcDNA) and T98G-GAB(+) cells (d). Protein damage are depicted for LN229 (LN229-GFP) and LN229-GLS(−) cells and for T98G (T98G-pcDNA) and T98G-GAB(+) cells (e). The results are expressed as nmol CRE/mg protein, pmol MDA/mg protein, and pmol carbonyl/mg protein, respectively, and a ratio was calculated by dividing the values of the experimental sample by those of its control for each biological replicate. Data are expressed as mean ± S.E.M. (N = 3) of the ratios calculated for at least three independent experiments. *p < 0.05, **p < 0.01, ***p < 0.001, compared to the respective control group. CRE Copper Reducing Equivalents, MDA malondialdehyde, TAC total antioxidant capacity
Discussion
Because cancer signalling pathways often exist in atypical forms, our understanding of cancer networks is currently not sufficient to allow a priori predictions of the cellular response [4]. However, combination therapies and sequential treatments to gain anticancer power are profusely recommended [7]. In GBMs, a thorough understanding of the vast genetic aberrations and signalling pathways involved in gliomagenesis as well as heterogeneous clinicopathological presentations remains elusive. The discovery of miRNAs and their capability of simultaneously regulating multiple downstream genes may play a key role in explaining the complex mechanisms underlying GBM formation [13]. Recent efforts against GBM malignancies have focused on therapies which target key intracellular apoptotic pathways which may confer tumour resistance, such as p53, and more recently some miRNA sequences [14]. Of particular interest in this research, direct repression of miRNA-23a and miRNA-23b by c-Myc is responsible for increased expression of GLS protein in human Burkitt lymphoma and prostate cancer [15]. In addition, miRNA-137 inhibits growth of malignant melanoma by targeting GLS [16]. On the other hand, expression of GLS2 by bladder cancer cells was regulated through interfering with miRNA-16 [11]. Furthermore, the capacity of GLS2 to alter miRNA patterns in cancer cells would be explained by recent findings showing that GLS2 is involved in miRNA regulation through Dicer stabilization [17]. As showed in Table 1, some of the sequences of the miRNAs up- or downregulated after changes in GLS or GLS2 expression are present in some genes that have been involved in glioma proliferation and inhibition of apoptosis as NGB [18], WWP2 [19] and PTTG1 [20]. Of interest, miRNA-1260a, miRNA-140-3p, miRNA-21-5p, and miRNA-146a-5p are among the validated by real-time qRT-PCR in each model with targeted GA expression, LN229-GLS(−), T98G-GLS(−), LN229-GAB(+) and T98G-GAB(+) (Fig. 1). On the other hand, miRNA-1246, which is the only that hybridize with both GLS(−) and GAB(+) transcripts (Table 1), was also the only miRNA overexpressed in all models, as shown in Fig. 1a–d.
Among miRNAs having notable expression differences between T98G-GAB(+) vs. T98G-pcDNA cells, and LN229-GAB(+) vs. LN229-pcDNA cells, miRNA-140-3p were increased 8.8- and 3.5-fold, respectively, as assessed by qRT-PCR. Although studies dealing with this particular miRNA in GBM or other brain tumours have not been yet published, our results showing a marked increase in the less-malignant and more differentiated GAB(+) cell line are in agreement with the downregulation of miRNA-140-3p found in lung cancer [21]. MiRNA-1246 huge expression (8.2- and 6.8-fold, respectively) in T98G and LN229 cells overexpressing GAB supports that this miRNA is present in less malignant human GBM cells [22]. MiRNA-1246 has been described as a novel interactor of p53 and its analogous p63 and p73, which are associated with apoptosis [23] and inhibition of growth in hepatocellular carcinoma (HCC) cells [23, 24]. Accordingly, miRNA-1246 targets p53, phosphatidylinositol-3-kinase/protein kinase B (PI3K/AKT), mitogen-activated protein kinase (MAPK), and mammalian target of rapamycin (mTOR) signalling pathways, leading to apoptosis in paclitaxel-treated HCC cells [23]. Of note, GAB has shown context-dependent tumour suppression activity as a target of the p53 family of tumour suppressors [2], and a recent work linked this GAB function with the negative regulation of PI3K/AKT pathway in hepatocellular carcinoma [25]. Thus, induction of miRNA-1246 after GAB overexpression may be one key mechanism to inhibit cancer progression through p53 and PI3K/AKT dependent pathways.
T98G and LN229 GAB overexpressing cells showed respectively a 6.5- and 21.4-fold overexpression of miRNA-1260a. Correspondingly, miRNA-1260a was up-regulated in the less proliferative atypical Spitz lesions relative to Spitz skin tumours [26]. MiRNA-21-5p was 3.4- and 3.2-fold expressed respectively in T98G and LN229 GBM cells overexpressing GAB. This miRNA targets the gene for p-21 protein, which is correlated with the anti-proliferative and cell death effects in T98G cells [27]. The tumour suppressing function of this miRNA was also characterized in thyroid papillary carcinoma via downregulation of Bcl-2 [28]. On the other hand, miRNA-146a-5p was 2.2- and 1.7-fold overexpressed in T98-GAB(+) and LN229-GAB(+) models, respectively. This is in accordance with its function as tumour suppressor in GBM by negatively regulating epidermal growth factor receptor (EGFR) and platelet-derived growth factor receptor beta (PDGFRB) genes [29]. MiRNA-146a-5p has been described as a native molecular brake for oncogenesis in aggressive and deadly brain tumours via NOTCH1 signalling [30]. It also functions as a tumour-suppressive miRNA targeting NOTCH2, inhibiting the epithelial-mesenchymal transition (EMT) progression of oesophageal squamous cell carcinoma (ESCC) [31]. Similarly, functional studies showed that this miRNA is a negative regulator of tumourigenic gene expression in microglia via its target SMAD4 [32]. On the other hand, miRNA-146a-5p inhibits cervical cancer cell proliferation, migration and invasion, repressing the expression of histone demethylase [33]. Of note, miRNA-146a-5p has been found downregulated in prostate cancer cells treated with the anticancer agent metformin [34]. Finally, concerning miRNA-92a-3p, although no differences were found in T98G-GAB(+) vs. T98G-pcDNA cells, it was 4.1-fold overexpressed in LN229-GAB(+) vs. LN229-pcDNA cells. This data agrees with a very recent publication showing that miRNA-92a-3p overexpression is a potential biomarker of favourable outcome of chronic lymphocytic leukaemia patients [35]. Importantly, miRNA-92a-3p profile was markedly different (even or upregulated) from that downregulation described below for both LN229-GLS(−) vs. LN229-GFP cells and T98G-GLS(−) vs. control T98G cells (0.75-fold and 0.79-fold, respectively).
GBM LN229-GLS(−) cells dramatically repressed expression of miRNA-146a-5p. Its expression in T98G-GLS(−) was also decreased (0.6-fold). Although it is currently unknown whether miRNA-146a-5p possesses oncogenic properties in GBM cells, results obtained in breast cancer tissue and MCF-7 cell line support this role through the BRCA1 pathways [36]. These findings are in agreement with the role of miRNA-146a as an inflammatory miRNA that regulates the NF-κB pathway, transforming growth factor-β (TGF-β), and induces angiogenesis [37]. Noteworthy, miRNA-146a-5p has been recognized as a positive biomarker in lung cancer [38] and as an oncogenic factor in oral squamous cell carcinoma (OSCC). Diverse investigations have pointed out that miRNA-146a-5p affects proliferation and apoptosis in a cellular context-dependent manner, defining to this miRNA as a pleiotropic regulator of carcinogenesis [39].
LN229-GLS(−) showed 0.05-fold expression of miRNA-140-3p vs. LN229-GFP cells. However, in T98G-GLS(−) vs. control T98G cells its expression was only decreased 0.85-fold. Although not previous results have been published in GBM cells, miRNA-140-3p was up-regulated in other malignancies [40]. Versus LN229 and T98G, in both LN229-GLS(−) and T98G-GLS(−) miRNA-21-5p was underexpressed vs. their controls (0.44-fold and 0.60-fold, respectively). The oncogenic function of miRNA-21 in GBM cells was associated to the p53 circuit, targeting TGF-β, modulating cell growth, and repressing p53-mediated apoptosis in response to anticancer drugs like doxorubicin [41]. MiRNA-21-5p is also linked to tumour development in colorectal cancer, modulating Wnt circuit, RAS-MAPK-PI3K-AKT signalling pathways, as well as TGF-β and p53 components [42]. MiRNA-21 regulated cell proliferation and apoptosis in ESCC through the AKT signalling pathway via targeting cell adhesion molecules [43]. MiRNA-21 has oncogenic properties that play important roles in many cancers, including leukaemia [44]. Multiple data show that miRNA-21 may be an independent prognostic marker for glioma proliferation and invasion, especially those gliomas with high pathological grades, and thus could also be a potential therapeutic target for molecular glioma therapy [45]. MiRNA-21 target sites include a number of genes involved in apoptosis, i.e.: expression of programmed cell death four gene (PDCD4) correlates inversely with expression of miRNA-21 in a number of human GBM cell lines such as T98G. Inhibition of miRNA-21 increases endogenous levels of PDCD4 in cell line T98G and overexpression miRNA-21 inhibits PDCD4-dependent apoptosis, functioning as an anti-apoptotic factor [46]. MiRNA-21 exerts its oncogenic function in gliomas through its post-transcriptional targets and downstream signal pathways: PDCD4, phosphatase and tensin homolog (PTEN), tumour protein 63 (TP63), tropomyosin 1 (TPM1), and tissue inhibitor of metalloproteases 3 (TIMP3) [47]. MiRNA-1260a was 0.58-fold and 0.65-fold underexpressed in LN229-GLS(−) and T98G-GLS(−) GBM cells. This data agrees with significant increase in the expression of miRNA-1260a in metastatic tumours as compared to the non-metastatic neuroblastoma in vivo [48]. MiRNA-92a was 0.74-fold and 0.80-fold underexpressed when GLS was silenced in LN229 and T98G cells. Properly, depleting endogenous miRNA-92a-3p inhibits the rate of cell viability, cell migration, and metastasis of U87 and U251 glioma cells in vitro while diminishing tumour volume and weight of xenograft in vivo [49]. Interestingly, miRNA-92a is overexpressed in human GBM cells, including 6-fold in LN229 cells. This research showed that Bim is a direct functional target of miRNA-92a and that high grade gliomas show augmented miRNA-92a levels [50]. MiRNA-92a-3p has showed oncogenic properties, i.e.: increasing cell proliferation, migration and invasion and inhibited cell apoptosis in vitro, as well as augmenting xenograft tumour formation in vivo, in a variety of cancers including leukaemia [44] and gastric [51] cancer. MiRNA-1246 was 2.2-fold and 1.5-fold overexpressed in GLS silenced LN229 and T98G cells, respectively, when compared with control LN229 and T98G cells. These results agree with findings in T98G-GAB(+) vs. T98G-pcDNA cells. In addition, upregulation of miRNA-1246 in prostate cancer significantly decreased growth, augmented apoptosis and inhibited proliferation, invasiveness, and migration in vitro and in vivo [52]. Similar conclusions were obtained in lung cancer [53]. Conversely, miRNA-1246 was overexpressed in MDA-MB-231 breast cancer cells [54].
Concerning antioxidant defence systems, although not similar studies have been found in literature, discordant results have been described, i.e.: antioxidant phenethyl isothiocyanate decreases SOD activity and GSH levels, inducing oxidative stress and repressing LN229 cell growth [55]; but apoptosis induced by cisplatin treatment diminished SOD to impair cancer growth [56]. It has been also reported differences between SOD and CAT responses to oxidative stress and apoptosis in LN229 cells, increasing SOD but decreasing CAT activity [57]. Accordingly, not concluding results were found in our models with reference to SOD and CAT activities (Additional file 5: Fig. S3c, d). Although ROS and antioxidant enzymes show a complex regulation [2], why a different tendency in SOD activity upon GAB overexpression in LN229 as compared to T98G cells is an open question. Impressive, total GSH was established as a clear biomarker of GA-mediated redox regulation: when GLS was silenced GSH levels declined and when GAB was overexpressed GSH content was increased (Additional file 5: Fig. S3a). Similarly, GR activity was another positive biomarker because rose following both GLS silencing and GAB overexpression (Additional file 5: Fig. S3b). Interestingly, GR has been associated to drug resistance in GBM cells, being a potential target for ameliorating GBM treatment [58]. In all models studied in this research, a reduction in both lipid peroxidation and protein damage was found (Additional file 5: Fig. S3e, g), supporting the antioxidant effect of altering GLS and GAB expression and subsequently increasing GR activity (Additional file 5: Fig. S3b). Results depicted in Additional file 5: Fig. S3f strengthen these conclusions: when GLS is silenced in both LN229 and T98G cells TAC is significantly dipped but when GAB was overexpressed in one and the other, LN229 and T98G cells, TAC climbed significantly. Therefore, in spite of scarce studies in GBM cells, changes in oxidative status and antioxidant capacity linked to differential GA expression [3] emerge as promising support against GBM aggressiveness.
In several GBM cell lines, including T98G and LN229, GAB increased the sensitivity towards ROS throughout downregulation of the PI3K/AKT cascade [12]. Changes in oxidative stress through PI3K/AKT linked to miRNA expression levels support our results [59]. Of note, upregulation of miRNA-34a, which is induced by GLS2 through Dicer [17], suppresses antioxidant circuits (SIRT1/PGC-1α/NRF2 pathway) and enhances susceptibility of wild-type p53 cancer cells to ROS [60]. Sirtuin 1 (SIRT1) is involved in the redox balance and carcinogenesis, including the proliferation of GBM cells [61]. Peroxisome proliferator-activated receptor γ coactivator1α (PGC-1α) regulates energy metabolism and oxidative stress, evoking a more aggressive behaviour of GBM cells [62]. Nuclear erythroid factor 2-related factor 2 (NRF2) is a redox-sensitive transcription factor regulating the expression of antioxidant enzymes such as SOD [60].
On the other hand, GLS2 displays a novel non-metabolic mechanism interacting and stabilizing Dicer protein to thrust miRNA maturation and inhibit Snail expression. These results indicate that GLS2 can act as a suppressor for migration and invasion of cancer cells, repressing EMT [17]. Appropriately, GLS and GLS2 target the expression of Snail through p53 regulation linking GAs to some microRNA expression, which translationally can repress Snail and oncogenic phenotype (Fig. 5). In fact, Snail induces a glycolytic switch to Warburg effect (higher glucose uptake and lactate production) by suppressing mitochondrial oxidative phosphorylation that may contribute to metastasis and tumour progression in breast, colon and ovary cancer [63]. Thus, GA can regulate the expression of some miRNAs, uncovering a previously unsuspected link between glutamine metabolism, miRNAs, ROS homeostasis and cancer, that complement data of other miRNAs regulating GA isoenzymes [11, 15, 16]. Future studies might involve investigating the possibility that GLS silencing and/or GAB overexpression can show synergic effects in combination with miRNA targeting as a new therapeutic strategy against cancer.
Some metabolic and signalling circuits that intersect with glutaminase isoenzymes. This network schematic representation shows both glutaminase isoenzymes, GLS and GLS2, executing a dual function. First, glutaminase is the limiting enzyme in glutaminolysis to generate glutamate from glutamine, that can be used for GSH synthesis, the most important intracellular antioxidant. Second, both isoenzymes can interact with other proteins and signalling pathways. GLS is regulated by the oncogenic transcription factor c-Myc and activates mTORC1, which is a major positive regulator of the Warburg effect. Tumour suppressor factor p53 is hampered by GLS, but stimulates GLS2. GLS2 blocks PI3K/AKT signalling pathway while interacts and stabilizes Dicer to mature several miRNAs which repress Snail favouring Warburg effect of tumour cells, and impacting oxidative status. P53 connects GLS, GLS2, microRNAs, ROS homeostasis and cancer. GLS glutaminase isoenzyme, GLS2 glutaminase 2 isoenzyme, GSH glutathione, LRH-1 nuclear receptor liver receptor homolog 1, mTORC1 mammalian target of rapamycin complex 1, PI3K/AKT phosphatidylinositol-3-kinase/protein kinase B, ROS reactive oxygen specie
Conclusions
Importantly, this study is the first to identify GA-dependent miRNA expression profiles in glutamine-addicted GBM cells with a consistent pattern of GA expression: high levels of GLS isoforms but only traces or lack of GLS2 transcripts [64]. Noteworthy, miRNA-153 was found to inhibit GBM cell growth in U87MG cell line, modulating GLS expression by directly targeting the 3′-UTR of GLS mRNA [65]. Although some miRNAs have been shown to regulate GA expression in certain types of cancers [66], it is currently unknown whether similar miRNA mechanisms would contribute to GLS activation and GLS2 repression in gliomas [64]. Of interest, a combination therapy which incorporates chemotherapy, radiotherapy, as well as GA- and miRNA-based therapy is likely to be the future direction for the treatment of gliomas [67]. Nonetheless, miRNAs studied in this article will require further functional analysis to evaluate whether they can play a direct role in the clinical outcome of patients suffering from GBM. Second, the marked change in glutathione and oxidative status following GLS silencing or GAB overexpression are a promising tool against cancer progression [68], including glioma cells [69]. Combinatorial and synergic strategies including glutaminase inhibitors are gaining important battles against cancer [70,71,72,73]. Isoenzyme specific GA expression [74] and associated miRNAs [75] can contribute to improve the efficiency of new drugs [76] and to block medicament resistance [77], targeting key tumour pathways [78], inducing cancer metabolic reprogramming [79], and helping to win this war in XXI century.
Availability of data and materials
All data generated in this study are available from the corresponding author on reasonable request.
Abbreviations
- AKT:
-
Protein kinase B
- ATCC:
-
American Type Culture Collection
- CAT:
-
Catalase
- CRE:
-
Copper Reducing Equivalents
- DMEM:
-
Dulbecco's Modified Eagle Medium
- DTNB:
-
5,5′-dithiobis-(2-nitrobenzoic acid)
- EGFR:
-
Epidermal growth factor receptor
- EMT:
-
Epithelial-mesenchymal transition
- ESCC:
-
Oesophageal squamous cell carcinoma
- FBS:
-
Foetal bovine serum
- FTC:
-
Fluorescein 5-thiosemicarbazide
- GA:
-
Glutaminase
- GAB:
-
Longer GLS2 transcript
- GAC:
-
Shorter GLS transcript
- GFP:
-
Green fluorescent protein
- GBM:
-
Glioblastoma
- GLS:
-
Glutaminase isoenzyme
- GLS2:
-
Glutaminase isoenzyme 2
- GR:
-
Glutathione reductase
- GSSG:
-
Oxidized glutathione
- GSH:
-
Reduced glutathione
- HCC:
-
Hepatocellular carcinoma
- KGA:
-
Longer GLS transcript
- LGA:
-
Shorter GLS2 transcript
- lncRNA:
-
Long non-coding ribonucleic acid
- MAPK:
-
Mitogen-activated protein kinase
- MDA:
-
Malondialdehyde
- miRNA:
-
microRNA
- mTOR:
-
Mammalian target of rapamycin
- NRF2:
-
Nuclear erythroid factor 2-related factor 2
- OSCC:
-
Oral squamous cell carcinoma
- PDCD4:
-
Programmed cell death four gene
- PDGFRB:
-
Platelet-derived growth factor receptor beta
- PGC-1α:
-
Peroxisome proliferator-activated receptor γ coactivator1α
- PI3K:
-
Phosphatidylinositol-3-kinase
- PTEN:
-
Phosphatase and tensin homolog
- qRT-PCR:
-
Quantitative reverse transcription-polymerase chain reaction
- RIN:
-
RNA integrity numbers
- RNAi:
-
RNA interference
- ROS:
-
Reactive oxygen species
- SD:
-
Standard deviation
- siRNA:
-
Small interferent RNA
- SIRT1:
-
Sirtuin 1
- SOD:
-
Superoxide dismutase
- TAC:
-
Total antioxidant capacity
- TBARS:
-
Thiobarbituric acid reactive substances
- tGSH:
-
Total glutathione
- TGF-β:
-
Transforming growth factor-β
- TIMP3:
-
Tissue inhibitor of metalloproteases 3
- TP63:
-
Tumour protein 63
- TPM1:
-
Tropomyosin 1
- UCA1:
-
Urothelial carcinoma-associated 1
References
Yuneva MO, Fan TW, Allen TD, Higashi RM, Ferraris DV, Tsukamoto T, et al. The metabolic profile of tumors depends on both the responsible genetic lesion and tissue type. Cell Metab. 2012;15:157–70.
Matés JM, Segura JA, Martín-Rufián M, Campos-Sandoval JA, Alonso FJ, Márquez J. Glutaminase isoenzymes as key regulators in metabolic and oxidative stress against cancer. Curr Mol Med. 2013;4:514–34.
Martín-Rufián M, Nascimento-Gomes R, Higuero A, Crisma AR, Campos-Sandoval JA, Gómez-García MC, et al. Both GLS silencing and GLS2 overexpression synergize with oxidative stress against proliferation of glioma cells. J Mol Med (Berl). 2014;92:277–90.
Márquez J, Tosina M, de la Rosa V, Segura JA, Alonso FJ, Matés JM, et al. New insights into brain glutaminases: beyond their role on glutamatergic transmission. Neurochem Int. 2009;55:64–70.
Kashyap D, Tuli HS, Garg VK, Goel N, Bishayee A. Oncogenic and tumor-suppressive roles of microRNAs with special reference to apoptosis: molecular mechanisms and therapeutic potential. Mol Diagn Ther. 2018;22:179–201.
Beyer S, Fleming J, Meng W, Singh R, Haque SJ, Chakravarti A. The role of miRNAs in angiogenesis, invasion and metabolism and their therapeutic implications in gliomas. Cancers (Basel). 2017;9:pii:E85.
Cheng T, Sudderth J, Yang C, Mullen AR, Jin ES, Matés JM, et al. Pyruvate carboxylase is required for glutamine-independent growth of tumor cells. Proc Natl Acad Sci U S A. 2011;108:8674–9.
Hu W, Zhang C, Wu R, Sun Y, Levine A, Feng Z. Glutaminase 2, a novel p53 target gene regulating energy metabolism and antioxidant function. Proc Natl Acad Sci U S A. 2010;107:7455–60.
Suzuki S, Tanaka T, Poyurovsky MV, Nagano H, Mayama T, Ohkubo S, et al. Phosphate-activated glutaminase (GLS2), a p53-inducible regulator of glutamine metabolism and reactive oxygen species. Proc Natl Acad Sci U S A. 2010;107:7461–6.
Liu F, Li Y, Liu G. MicroRNA-200c exacerbates the ischemia/reperfusion injury of heart through targeting the glutaminase (GLS)-mediated glutamine metabolism. Eur Rev Med Pharmacol Sci. 2017;21:3282–9.
Li HJ, Li X, Pang H, Pan JJ, Xie XJ, Chen W. Long non-coding RNA UCA1 promotes glutamine metabolism by targeting miR-16 in human bladder cancer. Jpn J Clin Oncol. 2015;45:1055–63.
Majewska E, Márquez J, Albrecht J, Szeliga M. Transfection with GLS2 glutaminase (GAB) sensitizes human glioblastoma cell lines to oxidative stress by a common mechanism involving suppression of the PI3K/AKT pathway. Cancers (Basel). 2019;11:115.
Mercatelli N, Galardi S, Ciafrè SA. MicroRNAs as multifaceted players in glioblastoma multiforme. Int Rev Cell Mol Biol. 2017;333:269–323.
Liu J, Albrecht AM, Ni X, Yang J, Li M. Glioblastoma tumor initiating cells: therapeutic strategies targeting apoptosis and microRNA pathways. Curr Mol Med. 2013;13:352–7.
Gao P, Tchernyshyov I, Chang TC, Lee YS, Kita K, Ochi T, et al. c-Myc suppression of miR-23a/b enhances mitochondrial glutaminase expression and glutamine metabolism. Nature. 2009;458:762–5.
Luan W, Zhou Z, Zhu Y, Xia Y, Wang J, Xu B. miR-137 inhibits glutamine catabolism and growth of malignant melanoma by targeting glutaminase. Biochem Biophys Res Commun. 2018;495:46–52.
Kuo TC, Chen CK, Hua KT, Yu P, Lee WJ, Chen MW, et al. Glutaminase 2 stabilizes Dicer to repress Snail and metastasis in hepatocellular carcinoma cells. Cancer Lett. 2016;383:282–94.
Zhang B, Liu Y, Li Y, Zhe X, Zhang S, Zhang L. Neuroglobin promotes the proliferation and suppresses the apoptosis of glioma cells by activating the PI3K/AKT pathway. Mol Med Rep. 2018;17:2757–63.
Liang J, Qi WF, Xie S, Wang WF, Zhang XL, Zhou XP, et al. Expression of WW domain-containing protein 2 is correlated with pathological grade and recurrence of glioma. J Cancer Res Ther. 2017;13:1032–7.
Su X, Chen J, Ni L, Shi W, Shi J, Liu X, et al. Inhibition of PTTG1 expression by microRNA suppresses proliferation and induces apoptosis of malignant glioma cells. Oncol Lett. 2016;12:3463–71.
Dong W, Yao C, Teng X, Chai J, Yang X, Li B. MiR-140-3p suppressed cell growth and invasion by downregulating the expression of ATP8A1 in non-small cell lung cancer. Tumour Biol. 2016;37:2973–85.
Ames HM, Yuan M, Vizcaíno MA, Yu W, Rodriguez FJ. MicroRNA profiling of low-grade glial and glioneuronal tumors shows an independent role for cluster 14q32.31 member miR-487b. Mod Pathol. 2017;30:204–16.
Yan H, Wang S, Yu H, Zhu J, Chen C. Molecular pathways and functional analysis of miRNA expression associated with paclitaxel-induced apoptosis in hepatocellular carcinoma cells. Pharmacology. 2013;92:167–74.
Zhang Q, Cao LY, Cheng SJ, Zhang AM, Jin XS, Li Y. p53-induced microRNA-1246 inhibits the cell growth of human hepatocellular carcinoma cells by targeting NFIB. Oncol Rep. 2015;33:1335–41.
Liu J, Zhang C, Lin M, Zhu W, Liang Y, Hong X, et al. Glutaminase 2 negatively regulates the PI3K/AKT signaling and shows tumor suppression activity in human hepatocellular carcinoma. Oncotarget. 2014;5:2635–47.
Latchana N, Regan K, Howard JH, Aldrink JH, Ranalli MA, Peters SB, et al. Global microRNA profiling for diagnostic appraisal of melanocytic Spitz tumors. J Surg Res. 2016;205:350–8.
Paolini A, Curti V, Pasi F, Mazzini G, Nano R, Capelli E. Gallic acid exerts a protective or an anti-proliferative effect on glioma T98G cells via dose-dependent epigenetic regulation mediated by miRNAs. Int J Oncol. 2015;46:1491–7.
Zhang H, Cai Y, Zheng L, Zhang Z, Lin X, Jiang N. LncRNA BISPR promotes the progression of thyroid papillary carcinoma by regulating miR-21-5p. Int J Immunopathol Pharmacol. 2018;32:2058738418772652.
Xu G, Li JY. Differential expression of PDGFRB and EGFR in microvascular proliferation in glioblastoma. Tumour Biol. 2016;37:10577–86.
Hu HQ, Sun LG, Guo WJ. Decreased miRNA-146a in glioblastoma multiforme and regulation of cell proliferation and apoptosis by target Notch1. Int J Biol Markers. 2016;31:e270–5.
Wang C, Zhang W, Zhang L, Chen X, Liu F, Zhang J, et al. miR-146a-5p mediates epithelial-mesenchymal transition of oesophageal squamous cell carcinoma via targeting Notch2. Br J Cancer. 2016;115:1548–54.
Karthikeyan A, Gupta N, Tang C, Mallilankaraman K, Silambarasan M, Shi M, et al. Microglial SMAD4 regulated by microRNA-146a promotes migration of microglia which support tumor progression in a glioma environment. Oncotarget. 2018;9:24950–69.
Peta E, Sinigaglia A, Masi G, Di Camillo B, Grassi A, Trevisan M, et al. HPV16 E6 and E7 upregulate the histone lysine demethylase KDM2B through the c-MYC/miR-146a-5p axys. Oncogene. 2018;37:1654–68.
Avci CB, Harman E, Dodurga Y, Susluer SY, Gunduz C. Therapeutic potential of an anti-diabetic drug, metformin: alteration of miRNA expression in prostate cancer cells. Asian Pac J Cancer Prev. 2013;14:765–8.
Papageorgiou SG, Diamantopoulos MA, Kontos CK, Bouchla A, Vasilatou D, Bazani E, et al. MicroRNA-92a-3p overexpression in peripheral blood mononuclear cells is an independent predictor of prolonged overall survival of patients with chronic lymphocytic leukemia. Leuk Lymphoma. 2019;60:658–67.
Gao W, Hua J, Jia Z, Ding J, Han Z, Dong Y, et al. Expression of miR-146a-5p in breast cancer and its role in proliferation of breast cancer cells. Oncol Lett. 2018;15:9884–8.
Simanovich E, Brod V, Rahat MM, Rahat MA. Function of miR-146a-5p in tumor cells as a regulatory switch between cell death and angiogenesis: macrophage therapy revisited. Front Immunol. 2018;8:1931.
Lv S, Xue J, Wu C, Wang L, Wu J, Xu S, et al. Identification of a panel of serum microRNAs as biomarkers for early detection of lung adenocarcinoma. J Cancer. 2017;8:48–56.
Min SK, Jung SY, Kang HK, Park SA, Lee JH, Kim MJ, et al. Functional diversity of miR-146a-5p and TRAF6 in normal and oral cancer cells. Int J Oncol. 2017;51:1541–52.
Asano M, Umezu T, Katagiri S, Kobayashi C, Tauchi T, Gotoh M, et al. Up-regulated exosomal miRNA-140-3p in CML patients with musculoskeletal pain associated with discontinuation of tyrosine kinase inhibitors. Int J Hematol. 2017;105:419–22.
Papagiannakopoulos T, Shapiro A, Kosik KS. MicroRNA-21 targets a network of key tumor-suppressive pathways in glioblastoma cells. Cancer Res. 2008;68:8164–72.
Falzone L, Scola L, Zanghì A, Biondi A, Di Cataldo A, Libra M, et al. Integrated analysis of colorectal cancer microRNA datasets: identification of microRNAs associated with tumor development. Aging. 2018;10:1000–14.
Li X, Chen D, Li M, Gao X, Shi G, Zhao H. The CADM2/Akt pathway is involved in the inhibitory effect of miR-21-5p downregulation on proliferation and apoptosis in esophageal squamous cell carcinoma cells. Chem Biol Interact. 2018;288:76–82.
Cui J, Li Q, Luo M, Zhong Z, Zhou S, Jiang L, et al. Leukemia cell-derived microvesicles induce T cell exhaustion via miRNA delivery. Oncoimmunology. 2018;7:e1448330.
Wu L, Li G, Feng D, Qin H, Gong L, Zhang J, et al. MicroRNA-21 expression is associated with overall survival in patients with glioma. Diagn Pathol. 2013;8:200.
Chen Y, Liu W, Chao T, Zhang Y, Yan X, Gong Y, et al. MicroRNA-21 down-regulates the expression of tumor suppressor PDCD4 in human glioblastoma cell T98G. Cancer Lett. 2008;272:197–205.
Li C, Sun J, Xiang Q, Liang Y, Zhao N, Zhang Z, et al. Prognostic role of microRNA-21 expression in gliomas: a meta-analysis. J Neurooncol. 2016;130:11–7.
Khan FH, Pandian V, Ramraj S, Aravindan S, Herman TS, Aravindan N. Reorganization of metastamiRs in the evolution of metastatic aggressive neuroblastoma cells. BMC Genom. 2015;16:501.
Song H, Zhang Y, Liu N, Zhao S, Kong Y, Yuan L. miR-92a-3p exerts various effects in glioma and glioma stem-like cells specifically targeting CDH1/β-catenin and notch-1/akt signaling pathways. Int J Mol Sci. 2016;17:pii: E1799.
Niu H, Wang K, Zhang A, Yang S, Song Z, Wang W, et al. miRNA-92a is a critical regulator of the apoptosis pathway in glioblastoma with inverse expression of BCL2L11. Oncol Rep. 2012;28:1771–7.
Zhang G, Li S, Lu J, Ge Y, Wang Q, Ma G, et al. LncRNA MT1JP functions as a ceRNA in regulating FBXW7 through competitively binding to miR-92a-3p in gastric cancer. Mol Cancer. 2018;17:87.
Bhagirath D, Yang TL, Bucay N, Sekhon K, Majid S, Shahryari V, et al. microRNA-1246 is an exosomal biomarker for aggressive prostate cancer. Cancer Res. 2018;78:1833–44.
Xu X, Cao L, Zhang Y, Lian H, Sun Z, Cui Y. MicroRNA-1246 inhibits cell invasion and epithelial mesenchymal transition process by targeting CXCR4 in lung cancer cells. Cancer Biomark. 2018;21:251–60.
Li XJ, Ren ZJ, Tang JH, Yu Q. Exosomal microRNA MiR-1246 promotes cell proliferation, invasion and drug resistance by targeting CCNG2 in breast cancer. Cell Physiol Biochem. 2017;44:1741–8.
Su JC, Lin K, Wang Y, Sui SH, Gao ZY, Wang ZG. In vitro studies of phenethyl isothiocyanate against the growth of LN229 human glioma cells. Int J Clin Exp Pathol. 2015;8:4269–76.
Duan L, Aoyagi M, Tamaki M, Yoshino Y, Morimoto T, Wakimoto H, et al. Impairment of both apoptotic and cytoprotective signalings in glioma cells resistant to the combined use of cisplatin and tumor necrosis factor alpha. Clin Cancer Res. 2004;10:234–43.
Kusaczuk M, Krętowski R, Naumowicz M, Stypułkowska A, Cechowska-Pasko M. Silica nanoparticle-induced oxidative stress and mitochondrial damage is followed by activation of intrinsic apoptosis pathway in glioblastoma cells. Int J Nanomedicine. 2018;13:2279–94.
Zhu Z, Du S, Du Y, Ren J, Ying G, Yan Z. Glutathione reductase mediates drug resistance in glioblastoma cells by regulating redox homeostasis. J Neurochem. 2018;144:93–104.
Docrat TF, Nagiah S, Krishnan A, Naidoo DB, Chuturgoon AA. Atorvastatin induces MicroRNA-145 expression in HEPG2 cells via regulation of the PI3K/AKT signalling pathway. Chem Biol Interact. 2018;287:32–40.
Do MT, Kim HG, Choi JH, Jeong HG. Metformin induces microRNA-34a to downregulate the Sirt1/Pgc-1α/Nrf2 pathway, leading to increased susceptibility of wild-type p53 cancer cells to oxidative stress and therapeutic agents. Free Radic Biol Med. 2014;74:21–34.
Li Q, Wang C, Cai L, Lu J, Zhu Z, Wang C, et al. miR–34a derived from mesenchymal stem cells stimulates senescence in glioma cells by inducing DNA damage. Mol Med Rep. 2019;19:1849–57.
Bruns I, Sauer B, Burger MC, Eriksson J, Hofmann U, Braun Y, et al. Disruption of peroxisome proliferator-activated receptor γ coactivator (PGC)-1α reverts key features of the neoplastic phenotype of glioma cells. J Biol Chem. 2019;294:3037–50.
Lee SY, Jeon HM, Ju MK, Jeong EK, Kim CH, Park HG, et al. Dlx-2 and glutaminase upregulate epithelial-mesenchymal transition and glycolytic switch. Oncotarget. 2016;7:7925–39.
Márquez J, Alonso FJ, Matés JM, Segura JA, Martín-Rufián M, Campos-Sandoval JA. Glutamine addiction in gliomas. Neurochem Res. 2017;42:1735–46.
Liu Z, Wang J, Li Y, Fan J, Chen L, Xu R. MicroRNA-153 regulates glutamine metabolism in glioblastoma through targeting glutaminase. Tumour Biol. 2017;39:1010428317691429.
Ortiz-Pedraza Y, Muñoz-Bello JO, Olmedo-Nieva L, Contreras-Paredes A, Martínez-Ramírez I, Langley E, et al. Non-coding RNAs as key regulators of glutaminolysis in cancer. Int J Mol Sci. 2020;21:2872.
McBrayer SK, Mayers JR, DiNatale GJ, Shi DD, Khanal J, Chakraborty AA, et al. Transaminase inhibition by 2-hydroxyglutarate impairs glutamate biosynthesis and redox homeostasis in glioma. Cell. 2018;175:101–16.
Matés JM, Campos-Sandoval JA, de los Santos-Jiménez J, Márquez J. Glutaminases regulate glutathione and oxidative stress in cancer. Arch Toxicol. 2020;94:2603–23.
Meng G, Xia M, Xu C, Yuan D, Schnurr M, Wei J. Multifunctional antitumor molecule 5’-triphosphate siRNA combining glutaminase silencing and RIG-I activation. Int J Cancer. 2014;134:1958–71.
Song M, Kim SH, Im CY, Hwang HJ. Recent development of small molecule glutaminase inhibitors. Curr Top Med Chem. 2018;18:432–43.
Akins NS, Nielson TC, Le HV. Inhibition of glycolysis and glutaminolysis: an emerging drug discovery approach to combat cancer. Curr Top Med Chem. 2018;18:494–504.
Wu C, Jin S, Chen L, Li H. Glutaminase inhibitors: a patent review. Expert Opin Ther Pat. 2018;28:823–35.
Matés JM, Campos-Sandoval JA, de los Santos-Jiménez J, Segura JA, Alonso FJ, Márquez J. Metabolic reprogramming of cancer by chemicals that target glutaminase isoenzymes. Curr Med Chem. 2020;27:5317–39.
Márquez J, Matés JM, Campos-Sandoval JA. Glutaminases Adv Neurobiol. 2016;13:133–71.
Xia HL, Lv Y, Xu CW, Fu MC, Zhang T, Yan XM, et al. MiR-513c suppresses neuroblastoma cell migration, invasion, and proliferation through direct targeting glutaminase (GLS). Cancer Biomark. 2017;20:589–96.
Xu X, Meng Y, Li L, Xu P, Wang J, Li Z, et al. Overview of the development of glutaminase inhibitors: achievements and future directions. J Med Chem. 2019;62:1096–115.
Shukla SK, Purohit V, Mehla K, Gunda V, Chaika NV, Vernucci E, et al. MUC1 and HIF-1alpha signaling crosstalk induces anabolic glucose metabolism to impart gemcitabine resistance to pancreatic cancer. Cancer Cell. 2017;32:71–87.
Herranz D, Ambesi-Impiombato A, Sudderth J, Sánchez-Martín M, Belver L, Tosello V, et al. Metabolic reprogramming induces resistance to anti-NOTCH1 therapies in T cell acute lymphoblastic leukemia. Nat Med. 2015;21:1182–9.
Matés JM, Campos-Sandoval JA, Márquez J. Glutaminase isoenzymes in the metabolic therapy of cancer. Biochim Biophys Acta Rev Cancer. 2018;1870:158–64.
Acknowledgements
Thanks are due to Monika Szeliga, Department of Neurotoxicology, Mossakowski Medical Research Centre, Polish Academy of Sciences, Warsaw, Poland. She kindly provided LN229-GAB(+) GBM cell model.
Funding
This work was financially supported by Ministerio de Ciencia y Tecnología of Spain, RTI2018-096866-B-I00 (to JMM and JM) and Ministerio de Educación of Spain PHB2010-0014-PC (to JMM). RJD is supported by the Howard Hughes Medical Institute, the National Cancer Institute, the Cancer Prevention and Research Institute of Texas, and the Moody Foundation. Thanks are also due to CAPES/DGU 250/11, Brazil (to FCS, RC and AC). JDSJ is granted by FPU17/04084, Ministerio de Ciencia, Innovación y Universidades.
Author information
Authors and Affiliations
Contributions
JMM designed research and wrote the paper. JDSJ, JACS, CMT, NUP, DB, JMC and FCS: performed experiments. JDSJ, MMR, AP, MCGG, CC, LC: cultures of cancer cells and supervised experiments. TC and RJDB: provided LN229 and LN229-GLS(−) cells. CL, FJA and JAS: supervised microarray experiments. JDSJ, JACS, JM and JMM: analysed data. JDSJ and JACS prepared the figures. JDSJ, JACS, FJA, JAS, RC, AC, RJD and JM: advised the design and reviewed the paper. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable since no patients subjects were used in this study.
Consent for publication
Not applicable. All authors read and approved the manuscript.
Competing interests
RJD is an advisor for Agios Pharmaceuticals.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Figure S1.
Analysis of GLS protein isoforms (KGA and GAC) downregulation in T98G-GLS(−) and its control (scramble siRNA-transfected) T98G.
Additional file 2: Figure S2.
MicroRNA expression analysis of T98G(T98G-pcDNA), T98G-GAB(+), LN229 (LN229-GFP), and LN229-GLS(−) cell samples, using miRXplore™ microarrays.
Additional file 3:
TableS1. Steep-loop sequences, mature sequences and probe sequences of the one thousand microRNAs used for microarray experiments.
Additional file 4:
TableS2. The microRNAs ratio list was determined for T98G-GAB(+) vs. T98G(T98G-pcDNA), depicted at column F, and for LN229-GLS(−) vs. LN229 (LN229-GFP) cells, depicted at column G.
Additional file 5: Figure S3.
Summary of all oxidative stress parameters in the four cell models studied LN229-GLS(−)/LN229 (LN229-GFP), T98G-GLS(−)/T98G (scrambled siRNA-transfected T98G), LN229-GAB(+)/LN229 (LN229-pcDNA) and T98G-GAB(+)/T98G(T98G-pcDNA), with the exception of Protein Damage parameter (g), which shows LN229-GLS(−)/LN229 (LN229-GFP) and T98G-GAB(+)/T98G (T98G-pcDNA).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
de los Santos-Jiménez, J., Campos-Sandoval, J.A., Márquez-Torres, C. et al. Glutaminase isoforms expression switches microRNA levels and oxidative status in glioblastoma cells. J Biomed Sci 28, 14 (2021). https://doi.org/10.1186/s12929-021-00712-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12929-021-00712-y
Keywords
- Antioxidant enzymes
- Cancer
- Glioblastoma
- Glutaminase
- microRNA
- Oxidative stress