Bioinformatics analysis
At first, the expression level of lncRNA MIAT was analyzed in common cancers as well as glioma by using the RNAseq data from TCGA (The Cancer Genome Atlas). In this regard, the data of each cancer were downloaded in RAW format (HTSeq-Counts). The cancers include BRCA (Breast invasive carcinoma), KIRC (Kidney renal clear cell carcinoma), KIRP (Kidney renal papillary cell carcinoma), KICH (Kidney Chromophobe), LIHC (Liver hepatocellular carcinoma), LUAD (Lung adenocarcinoma), LUSC (Lung squamous cell carcinoma), PRAD (Prostate adenocarcinoma), UCEC (Uterine Corpus Endometrial Carcinoma), STAD (Stomach adenocarcinoma), COAD (Colon adenocarcinoma), GBM (Glioblastoma multiform), and LGG (Brain Lower Grade Glioma). Using edgeR package (RRID:SCR_012802), gene expressions that were zero or close to zero were removed from the raw data by the counts-per-million (CPM) criterion. In the next step, the data were normalized using limma package (RRID:SCR_010943) based on a trimmed mean of M values (TMM) method, and the data were calculated in a logarithmic mode based on 2. The resulting expression matrix was applied to investigate the expression of lncRNA MIAT in each cancer as well as the correlation test to draw the co-expression network.
Patients
Twenty-seven brain tumor tissues were gathered from the Iran National Tumor Bank, (Tehran, Iran). They were kept at -180oC until being applied for the extraction of total RNA. The clinicopathological features of the participants were accompanied by the samples obtained from the Tumor Bank.
Cell culture
Two human Glioma cell lines, U-87MG (RRID: CVCL_0022) and A172 (RRID: CVCL_0131), were used. The cell lines were taken from the national cell bank of Iran (Pasteur Institute, Tehran, Iran). The U-87MG and A172 cells were cultured in RPMI 1640 (Gibco, CA, USA) and high glucose DMEM (Gibco, Germany) medium, respectively. The mediums were enriched with 10 % fetal bovine serum (FBS), 10 μg/mL streptomycin, and 100 U/mL penicillin, and maintained in a humidified 5% CO2 incubator.
RNA isolation and gene expression analysis
Total RNA was extracted by Trizol (Invitrogen, USA) reagent following the recommended instructions. The isolated total RNA was converted to cDNA by MMLV reverse transcriptase enzyme (200 U/μL) (Fermentas, USA), RNase inhibitor (20 U), dNTP mix (1mM) along with the random hexamer primer after treating with RNase-free DNase (Thermo Scientific, USA) to remove the genomic DNA contamination as described (Scientific). Specific primers were designed using Gene Runner software, version 4.0 (Table S1). QRT-PCR was performed by SYBR Premix Ex TaqTM II (Takara, Japan) with a Rotor-Gene 6000 detection system (Corbett Life Science, Sydney, Australia). The qRT-PCR reaction was carried out with the following condition: initiation at 95 °C for 30 seconds, amplification for 45 cycles, denaturation at 95 °C for 5 seconds, and annealing and extension at 60 °C for 30 seconds (adjusted according to the Tm of the primers). The accuracy of the PCR amplification products was confirmed on a 2% agarose gel. β-Actin was as an internal control to normalize qRT-PCR analyses.
LncRNA MIAT knockdown
LncRNA MIAT SMARTpool siRNA and scramble siRNA(s) were purchased from Dharmacon (siRNA SMARTpool® IKKα, Colorado, USA). The SMARTpool contained four modified siRNAs that guaranteed the efficient knocking down of lncRNA MIAT without off-target. siRNAs were introduced to glioma cells along with Lipofectamine 2000 (Invitrogen, USA) based on the suggested protocol. In brief, 2×104 cancer cells were seeded per well on a 12-well plate in an antibiotic-free growth medium until reach 40-60% confluency. Then, the proper concentration of siRNAs and Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA) according to the recommended protocol was diluted in 250 µL Opti-MEM (Invitrogen, Carlsbad, CA, USA). The resulting compound was incubated for 15 minutes at room temperature, and then mixed with the cultured cancer cells. The transfected cells were further kept in a humidified 5% CO2 atmosphere, for 2 or 3 days. Then, they were used for different molecular and cellular assays. To verify the efficient knocking-down of lncRNA MIAT, its expression was measured 2 and 3 days after transfection.
Clonogenic assay
To assess the influence of lncRNA MIAT knocking-down on the clonogenic capacity of glioma cancer cells, colony formation test was done (Franken et al. 2006). After transfection of cancer cells with siRNAs, one hundred of U-87MG-MIAT-depleted and A172-MIAT-depleted cells were transferred to a new 6-well plate and maintained for 1–2 weeks at 37 °C with 10% FBS. After that, the cells (colonies) were fixed with 4% paraformaldehyde and 0.2 % glutaraldehyde, and after that, stained with crystal violet (0.1%) (Sigma-Aldrich, Vienna, Austria) for 30 minutes. Finally, the colonies were counted in each well by ImageJ software (Version: 1.52v, NIH, USA, RRID: SCR_003070).
Cell cycle assay
Cell cycle assay was done by a nuclear propidium iodide (PI)-staining procedure(Krishan 1975). In brief, the transfected cells were detached, and washed with phosphate-buffered saline (PBS). The suspended cells were stained with PI staining solution (50 μg/mL propidium iodide, 0.1 % sodium citrate and 0.1 % Triton X 100). Then, the cells were kept for 30 minutes at room temperature without light(Fried et al. 1976). The cell cycle distribution was analyzed with a flow cytometer (Partec, Germany). The results were analyzed by flowjo 7.6.1 software (Ashland, OR, USA, RRID:SCR_008520).
Apoptosis assay
The occurrence of apoptosis was detected 48 hours after lncRNA MIAT siRNA and scramble siRNA transfection by staining the cells with Annexin V-FITC and PI apoptosis kit based on the suggested guidelines (Thermo Fisher Scientific, USA). After incubating the cells at room temperature for 20 minutes in dark, the apoptotic cells were detected with a flow cytometer (Partec, Germany). To verify the apoptosis percentage obtained from above staining, we used the CellEventTM Caspase-3/7 Green Flow Cytometry Assay Kit (Invitrogen, USA). The kit detected the activated Caspase-3/7 in apoptotic cells, which were analyzed by flow cytometry. Data analysis was carried out by flowjo 7.6.1 software (Ashland, OR, USA, RRID:SCR_008520).
Senescent cell detection
The activity of Lysosomal β-galactosidase, known as an indicator for detecting the senescent cells, was determined as explained previously (Keshavarz and Asadi 2019). Briefly, 104 cells were fixed with 4% paraformaldehyde and 0.2 % glutaraldehyde. Then, the cells were stained with 5-bromo-4-chloro-3-3indolyl β-D-galactoside (X-gal) that contained 1 mg/mL X-gal, 40 mM citric acid/sodium phosphate, 0.15 M NaCl, 2 mM MgCl2, 5 mM potassium ferrocyanide, and 5 mM potassium ferricyanide (pH 5.9–6.1). The cells were kept without light at 37 °C overnight and then they were visualized under a light microscope, and the stained cells were counted in three random areas.
In vitro scratch assay
In vitro scratch test was done to evaluate the effect of MIAT knockdown on cancer cell migration(Liang et al. 2007). In this regard, the transfected cells were seeded on 12-well plates until they reached approximately 90%–95% confluency. Then, the medium was removed and a cell-free gap was scratched to monolayer cells employing. The speed of wound closure was photographed in different time intervals (12 and 24 hours). To determine the migration percentage, the Image J software (NIH, Bethesda, MD, USA, RRID:SCR_003070) was applied to measure the gap distance.
Immunofluorescence staining
The expression level of E-Cadherin, Snail, and Vimentin was measured by Human EMT 3-Color Immunocytochemistry Kit (R&D Systems, USA) as explained before (Keshavarz and Asadi 2019). In summary, the 3×104 cells were seeded on 8 well cell culture chamber slides. Two days after transfection with lncRNA MIAT siRNA and scramble siRNA, the cells were fixed with 4% paraformaldehyde (PFA). Afterwards, they were kept in 0.3%Triton X-100, 10% normal donkey serum, and 1% BSA. Finally, the cells were incubated with monoclonal fluorochrome-conjugated primary antibodies (NL637-conjugated Goat Anti-Human E-Cadherin, NL557-conjugated Goat Anti-Human Snail, and NL493-conjugated Goat Anti-Human Vimentin). Images were taken under a fluorescence microscope (Zeiss, Germany), and 40x magnification. Quantification of the aforementioned proteins was done by measuring their color intensity using MATLAB software (MathWorks, Natick, MA, USA, RRID:SCR_013499).
Cell lysis and immunoblot analysis
Cell lysis and immunoblot were done as described previously (Alipoor et al. 2018). In summary, 72 hours after transfection of the cells with lncRNA MIAT siRNA and scramble siRNA, 106cells were homogenized in NP40 lysis buffer (Life Technologies, USA) and protease inhibitor cocktail (Sigma-Aldrich, USA). After that, 40 µg of total extracted protein was separated by running it on a 12% SDS-polyacrylamide gel and then it was transferred to a polyvinylidene fluoride membrane (PVDF; Roche). Then, blocking was done with Tris-buffered saline containing 5% nonfat dry milk and 0.5% Tween-20 for 1 hour at the room temperature. In the following step, membranes were incubated with primary antibody β-actin (Santa Cruz Biotechnology Cat# sc-47778 HRP, RRID:AB_2714189), SOX2 (Abcam, Cambridge, UK Cat# ab97959, RRID:AB_2341193), OCT4 (Abcam, Cambridge, UK Cat# ab19857, RRID:AB_445175), Nanog (Abcam, Cambridge, UK Cat# ab80892, RRID:AB_2150114), MAP1LC3 (Abcam, Cambridge, UK Cat# ab48394, RRID:AB_881433), and SQSTM1/p62 (Abcam, Cambridge, UK Cat# ab91526, RRID:AB_2050336) at 4 °C overnight. Next, the membranes probed with secondary horseradish peroxidase-conjugated secondary antibody (1: 2000 GE Healthcare Bio-Sciences, Chicago, IL, USA) at 4 °C overnight. Finally, the visualization of protein bands was done using ECL solution (ab65623, Abcam, Cambridge, UK) and displayed to Lu mi-Film chemiluminescent detection film (Roche, USA). The protein intensity was determined by Image LabTM Software (Biorad, USA). β-actin was considered as the internal control for normalizing the expression of proteins.
Construction of the competitive endogenous RNA (ceRNA) network
The ceRNA network was used to identify the pathways that lncRNA MIAT can act through miRNA sponging. Initially, all miRNA that bind to lncRNA MIAT were extracted using DIANA Tools (diana.e-ce.uth.gr/ LncBase v3.0, RRID:SCR_016510)(Paraskevopoulou et al. 2016). For this purpose, brain tissue, direct validation type, high confidence level and Homo sapiens species criteria were considered in the software, and the score ≥ 0.4 was used as the cut off value for the prediction analysis in the LncBase prediction section. Then, using miRWalk (http://mirwalk.umm.uni-heidelberg.de), the target genes (mRNAs) for the extracted miRNAs in the previous step, were identified. Eventually, the network was constructed using Cytoscape software (Seattle, USA, v3.3.0, RRID: SCR_003032).
Data analysis and statistics
TCGA data were analyzed using R programming language (V 4.5). Cytoscape and GlueGO software (Seattle, USA) were used for enrichment analysis, and drawing the co-expression network. Pearson correlation coefficient test was applied to investigate the correlation between lncRNA MIAT and other genes. All statistical analyses were done by GraphPad Prism 8 software (San Diego, CA, USA, RRID:SCR_002798). A one-way ANOVA followed by a Tukey posttest or Student’s t-test were applied to specify the significance level in the experiments. P-values>0.05 were considered as no significant statistical level. Data are reported as mean ± Standard Deviation (SD) at triplicate independent repeats.