MicroRNA-1225-5p behaves as a tumor suppressor in human glioblastoma via targeting of IRS1

Abstract

Background

MicroRNAs (miRNAs) play an important role in cancer initiation, progression, and metastasis by directly regulating their target genes.

Materials and methods

In this study, we observed that the miR-1225-5p expression level in glioblastoma tissues was significantly lower as compared with that in normal brain tissues, and its low expression was significantly associated with histopathological grade and poor patient prognosis.

Results

Through establishing a miR-1225-5p overexpression glioblastoma cell line, we found that ectopic overexpression of miR-1225-5p inhibited the proliferation, migration, and invasion of glioblastoma cells in vitro. Moreover, the growth of a glioblastoma xenograft tumor was attenuated by overexpression of miR-1225-5p. Further integrative studies suggested that the insulin receptor substrate 1 (IRS1) was a direct functional target of miR-1225-5p in glioblastoma, and the mRNA and protein levels of IRS1 in six human glioblastoma cell lines (A172, SW1783, U87, LN-229, SW1088, and T98G) were significantly higher as compared with normal human astrocytes.

Conclusion

These results suggest that miR-1225-5p may be a novel candidate for glioblastoma therapy.

Keywords:

• microRNA
• microRNA-1225-5p
• glioblastoma
• insulin receptor substrate 1
• prognosis

Introduction

As one of the most aggressive and common malignancies of the central nervous system in humans, glioblastoma is characterized by rapid infiltrative growth and cellular heterogeneity.^1,2 Patients with high-grade glioblastoma have worse outcomes and shorter survival, with a median overall survival time of only 14.6 months as compared with patients who have low-grade tumors.^3,4 Complete removal of glioblastoma tumor cells infiltrating or migrating to surrounding brain tissues by surgery is almost impossible. Besides surgery, radiotherapy and chemotherapy remain the mainstay of treatment for high-grade astrocytoma. However, due to the high invasiveness of glioblastoma cells and resistance to treatment, tumor recurrence and progression occur in >90% of glioblastoma patients shortly after surgery.^5–8 Therefore, understanding the molecular basis and identifying the key molecule(s) responsible for the malignant phenotype of glioblastoma will help to identify potentially applicable therapeutic agents, and improve the therapeutic efficacy of the management of glioblastoma patients.

MicroRNAs (miRNAs) are a group of evolutionarily conserved, small, endogenous, non-coding RNAs that play important roles in the suppression of translation or induction of mRNA degradation as post-transcriptional gene regulators, by directly pairing and binding with the 3′-untranslated regions (UTRs) of specific target messenger RNAs (mRNAs).^9–11 Accumulating evidence suggests that miRNAs are involved in a number of diverse physiological and pathological processes, particularly cancer pathogenesis, by directly targeting key oncogenes or tumor suppressors.^12–16 It has been reported that miR-1225-5p is downregulated in human gastric carcinoma, and inhibits the proliferation and metastasis of gastric carcinoma cells.¹⁷ However, whether miR-1225-5p plays a role in glioblastoma remains largely unknown.

In the present study, we attempted to determine the role of miR-1225-5p in human glioblastoma. We observed that the downregulation of miR-1225-5p in glioblastoma was significantly associated with histopathological grade and poor patient prognosis. Mechanistic investigations revealed that miR-1225-5p inhibits the proliferation, migration, and invasion of glioblastoma cells in vitro. Moreover, the growth of glioblastoma xenograft tumors was attenuated by overexpression of miR-1225-5p. In addition, the insulin receptor substrate 1 (IRS1) was confirmed to be a direct functional target of miR-1225-5p in glioblastoma.

Materials and methods

Cell culture

Normal human astrocytes and six human glioblastoma cell lines (A172, SW1783, U87, LN-229, SW1088, and T98G) were purchased from the Type Culture Collection of the Chinese Academy of Sciences (Shanghai, People’s Republic of China). Apart from the astrocytes, which were cultured in Prigrow IV medium (ABM, Richmond, BC, Canada), all the other cells were maintained in DMEM (Gibco, Shanghai, People’s Republic of China) supplemented with 10% FBS, 100 U/mL penicillin, and 100 mg/mL streptomycin, and were incubated at 37°C in a humidified incubator supplemented with 5% CO₂.

Clinical samples

A total of 44 glioma specimens and five normal brain samples frozen in liquid nitrogen were obtained from the China-Japan Union Hospital of Jilin University between January 2014 and May 2016. The study protocol was approved by the Ethics Review Committee of China-Japan Union Hospital of Jilin University, and written informed consent was obtained from all patients.

Generation of a cell line stably overexpressing miR-1225-5p (U87-1225)

The lentiviral plasmid pEZXMR03 (GeneCopoeia) expressing miR-1225-5p and Lenti-Pac HIV Expression Packaging mix (GeneCopoeia) were co-transfected into 293 T cells using Lipofectamine 2000 (Thermo Fisher Scientific, Waltham, MA, USA). After transfection for 48 hours, lentiviral particles were harvested and then transduced into the U87 glioblastoma cells, and the stably transfected cells were selected using puromycin and validated by reverse transcription PCR (RT-PCR).

Transfection and siRNA knockdown

SW1088 cells were transfected with siRNA sequence targeting IRS1 (AACCGGUUUCGAAAGAGA) using Lipofectamine 2000 (Invitrogen) according to the manufacturer’s instructions.

RNA extraction and quantitative RT-PCR

Total RNA was extracted from cultured cells or frozen tissues using the miRNeasy Mini kit (Qiagen, Hilden, Germany) and 1 mg RNA was reverse-transcribed using the miScript Reverse Transcription kit (Qiagen) for first-strand complementary DNA synthesis. Detection of miR-1225-5p and U6 utilized their specific primers separately. All the primers were synthesized by Invitrogen. The expression of miR-1225-5p was analyzed by quantitative PCR (qPCR) using the TaqMan microRNA assay kit (Applied Biosystems, Life Technologies).

Western blotting analysis

Western blotting was conducted according to the standard protocol with primary antibodies against IRS1, p-AKT, and AKT, and with β-catenin used as an internal control. All primary antibodies were purchased from Cell Signaling Technology, Inc. (Beverly, MA, USA).

Cell proliferation and colony-formation assay

Cell proliferation was detected using the MTT assay, according to the manufacturer’s instructions. Briefly, the cells were seeded in 96-well culture plates and incubated at 37°C with 5% CO₂ for 3 days. At each designated time point, 10 μL MTT was added to the medium and then incubated at 37°C for 1 hour. Cell numbers were estimated by measuring the absorbance at 450 nm using a 96-well plate reader. For the colony-formation assay, 5,000 cells were plated in a six-well plate for 12 days. Colonies were subjected to methanol/acetone (1:1) fixation and stained with crystal violet.

Cell invasion assays

Cell invasion was assessed using Transwell chambers (24-well inserts; 8 μm-pore size; Millipore, Billerica, MA, USA) that had been pre-coated with Matrigel (200 μg/mL; BD Biosciences, Franklin Lakes, NJ, USA). Briefly, cells were added into the upper chamber and complete DMEM with 20% FBS was placed in the lower chamber. After 24 hours of incubation at 37°C with 5% CO₂, the non-invading cells on the upper chamber were removed by wiping with a cotton swab. The cells that had infiltrated from Matrigel into the pores of the inserted filter were fixed and stained. The cell number in five randomly selected fields of view was visualized using a light microscope and counted.

Luciferase reporter assay

The 3′-UTR sequence of IRS1 was predicted to interact with miR-1225-5p, or a mutated sequence within the predicted target sites was synthesized and inserted into the XbaI and FseI sites of the pGL3 control luciferase reporter vector (Promega, Madison, WI, USA). Luciferase activity was measured 48 hours post-transfection using the Dual-Luciferase Reporter Assay system, according to the manufacturer’s instructions (Promega). For each sample, Renilla luciferase activity was normalized to firefly luciferase expression.

Xenograft assays in nude mice

All animal experiments were approved by the China-Japan Union Hospital of Jilin University Animal Care and Use Committee and followed up with the United Kingdom Coordinating Committee for Cancer Research (UKCCCR) with regard to the welfare of the animals. A total of 2×10⁶ U87 or U87-1225 cells were injected subcutaneously into the dorsal flank of BALB/c nude mice. Six mice were included in one experimental group. The tumor diameters were measured every 3 days, and tumor volumes were calculated (width²×length×0.5). The mice were sacrificed 22 days after tumor implantation.

Statistical analysis

Statistical analyses were conducted using GraphPad Prism IV software. The Student’s t-test was used to analyze the results expressed as mean ± SD. The survival curves were plotted using Kaplan–Meier analysis. Differences were considered to be statistically significant at P<0.05.

Results

miR-1225-5p is downregulated in glioblastoma cell lines and tissues

We first analyzed the expression of miR-1225-5p in six human glioblastoma cell lines (A172, SW1783, U87, LN-229, SW1088, and T98G). Compared with normal human astrocytes, the expression level of miR-1225-5p in these glioblastoma cell lines was significantly lower (Figure 1A). We further investigated the expression of miR-1225-5p in 44 glioblastoma samples (Table 1) and five normal brain tissue samples. The expression level of miR-1225-5p in the glioblastoma samples was also found to be significantly lower as compared with that in normal brain tissues (Figure 1B). Of note, the expression level of miR-1225-5p in advanced-stage (III and IV) glioblastoma was markedly lower when compared with that in earlier-stage (I and II) samples (Figure 1C).

Table 1 Clinicopathological features of 44 patients with glio-blastoma

Features	WHO I	WHO II	WHO III	WHO IV
No of patients	10	6	13	15
Mean age (years)	45.3	49.5	48.6	51.2
Gender
Male	6	3	7	9
Female	4	3	6	6
Surgery
Gross total resection	10	6	7	6
Partial resection	0	0	4	6
Biopsy	0	0	2	3
Adjuvant treatment
Radiotherapy	0	0	9	9
Chemotherapy	0	0	0	2
Radiotherapy and chemotherapy	0	0	4	4

Figure 1 miR-1225-5p is downregulated in glioblastoma and correlated with patient survival.

Notes: (A) The expression of miR-1225-5p in glioblastoma cell lines (A172, SW1783, U87, LN-229, SW1088, and T98G) and normal human astrocytes was measured by reverse transcription PCR (RT-PCR) and normalized against U6. (B) The expression of miR-1225-5p in five normal brain tissues and 44 glioblastoma samples was measured by RT-PCR and normalized against U6. (C) The expression of miR-1225-5p in stage I–II and III–IV glioblastoma tumors was measured by RT-PCR and normalized against U6. (D) The correlation between the miR-1225-5p level (low or high) with overall patient survival was assessed with Kaplan–Meier analysis. **P<0.01 vs astrocytes; **P<0.01 or ***P<0.001 vs indicated group.

Downregulation of miR-1225-5p is associated with disease progression and poor prognosis of patients with glioblastoma

To further investigate whether downregulation of miR-1225-5p is associated with glioblastoma progression and poor prognosis, the 44 abovementioned patients were divided into two groups – namely, low- and high-level expression of miR-1225-5p, based on the median value of the miR-1225-5p expression level among all patients (Table 2). Although the correlation between the miR-1225-5p expression level and the gender or age of the patients was not significant, the level of miR-1225-5p was negatively correlated with histopathological grade. Furthermore, as shown in Figure 1D as well as Figure S1A and B, the overall survival time of low miR-1225-5p expression patients was significantly shorter compared with that of the high miR-1225-5p expression group. Thus, miR-1225-5p suppresses the proliferation of glioblastoma cells in vitro and in vivo.

Table 2 Correlation between miR-1225-5p expression level and clinicopathological characteristics in 44 patients with glioblastoma

Factor	Low-level miR-1225-5p expression (n=25)	High-level miR-1225-5p expression (n=18)	P-value
Gender			0.63
Male	13	12
Female	8	11
Age (median), in years			0.55
>48	10	11
<48	12	11
Histopathologic grade			0.023
Low grade (I+II)	5	11
High grade (III+IV)	20	8

To investigate the role of miR-1225-5p in glioblastoma, we generated a glioblastoma cell line (U87-1225) stably overexpressing miR-1225-5p using retroviral vectors. As shown in Figure 2A, the miR-1225-5p expression level in U87-1225 cells was seven-fold higher when compared with that in parental U87 cells. The MTT assay demonstrated that overexpression of miR-1225-5p significantly suppressed the proliferation of U87 glioblastoma cells (Figure 2B). Furthermore, after being cultured for 14 days, the number of colonies in the U87-1225 group was significantly lower when compared with that in parental U87 cells (Figure 2C and D). Similar results were found in another glioblastoma cell line SW1088 (Figure S2A and D). The overexpression of miR-1225-5p significantly suppressed proliferation and colony formation of SW1088 cells. Moreover, miR-1225-5p significantly inhibited xenograft tumor growth in nude mice (Figure 2E and F).

Figure 2 Overexpression of miR-1225-5p suppresses the proliferation of glioblastoma cells in vitro and in vivo.

Notes: (A) The expression of miR-1225-5p in U87 and U87-1225 (stably expressing miR-1225-5p) cells was measured by reverse transcription PCR (RT-PCR) and normalized against U6. (B) U87 and U87-1225 cells were subjected to an MTT assay to assess cell proliferation at each time point. (C) Representative image of crystal violet-stained cell colonies analyzed by colony-formation assay at 14 days. (D) Quantification plots of (C). (E) Tumor-growth curves demonstrated that miR-1225-5p inhibited xenograft tumor growth in nude mice. (F) Relative tumor weight at 22 days after inoculation. **P<0.01 or ***P<0.001 vs U87 group.

miR-1225-5p suppresses the migration and invasion of glioblastoma cells

To further investigate the role of miR-1225-5p in glioblastoma, we next examined the effect of miR-1225-5p on the migration and invasion of U87 and SW1088 glioblastoma cell. In the scratch assay, overexpression of miR-1225-5p significantly attenuated wound healing (Figure 3A and B, Figure S3A and B). Moreover, we observed that overexpression of miR-1225-5p significantly abrogated the invasive ability of glioblastoma cells (Figure 3C and D, Figure S3C and D).

Figure 3 Overexpression of miR-1225-5p suppresses migration and invasion of glioblastoma cells.

Notes: (A) Representative picture of U87 and U87-1225 cell migration was measured by scratch wound assay. (B) The wound-healing percentage is plotted in (A). (C) Representative image of U87 and U87-1225 cell invasion using Matrigel-based invasion assay. (D) The average number of invading cells of five randomly selected fields from each sample was calculated. **P<0.01 or ***P<0.001 vs U87 group.

IRS1 is a direct target of miR-1225-5p in glioblastoma

To investigate the mechanism underlying the inhibitory effect of miR-1225-5p on the proliferation, migration, and invasion of glioblastoma cells, the activation of the PI3K/AKT pathway in U87 cells was evaluated by Western blotting. We found that the overexpression of miR-1225-5p reduced both the expression and phosphorylation level of IRS1 and the phosphorylation of AKT, whereas it did not affect the total AKT expression (Figure 4A). Further, RT-PCR analysis demonstrated that overexpression of miR-1225-5p decreased the IRS1 mRNA level (Figure 4B). It was previously suggested that IRS1 is a direct target of miR-1225-5p in gastric carcinoma.¹⁷ As shown in Figure 4C, the relative luciferase activity of the pGL3-IRS1-3′-UTR-wt reporter was significantly inhibited by overexpression of miR-1225-5p; however, the luciferase activity of the pGL3-IRS1-3′-UTR-mutant was barely affected. Furthermore, the mRNA and protein levels of IRS1 in six human glioblastoma cell lines (A172, SW1783, U87, LN-229, SW1088, and T98G) was significantly higher when compared with that in normal human astrocytes (Figure 4D). Next, we conducted an expression correlation analysis. We identified 63 glioblastoma samples in The Cancer Genome Atlas. We determined the correlation of miR-1225-5p with IRS1. Our results showed that IRS1 negatively correlated miR-1225-5p (Figure 4E). To confirm miR-1225-5p targeting IRS1 in glioblastoma cell, we knocked-down IRS1 in SW1088 cells (Figure S4A). As shown in Figure S4B–F, when compared to si con cells, knockdown IRS1 suppressed proliferation, migration, and invasion in SW1088 cells.

Figure 4 Insulin receptor substrate 1 (IRS1) is a direct target of miR-1225-5p in glioblastoma.

Notes: (A) Protein expression in U87 and U87-1225 cells was measured by Western blotting. (B) The expression of IRS1 in U87 and U87-1225 (stably expressing miR-1225-5p) cells was measured by reverse transcription PCR (RT-PCR) and normalized against U6. (C) Top, wild-type (wt) miR-1225-5p target sequences of IRS1 3′-untranslated region (UTR) and mutant (mut) miR-1225-5p target sequences of IRS1 3′-UTR. Bottom, relative luciferase activity of wt or mut pGL3-IRS1-3′-UTR in U87 and U87-1225 cells. (D) The mRNA (top) and protein (bottom) expression level of IRS1 in glioblastoma cell lines (A172, SW1783, U87, LN-229, SW1088, and T98G) and normal human astrocytes was measured by Western blotting or RT-PCR. (E) Scatter plot of the RNA Sequence data (TCGA) of 63 samples from glioblastoma patients. The red line indicates an inverse correlation of expression between the samples for miR-1225-5p and IRS1 genes (R=−0.5132; P=2.248e-05). **P<0.01.

Discussion

Based on the efforts of research groups worldwide, the abnormal expression and related role of miRNAs have been well established in certain types of cancer, acting as either oncogenes or tumor suppressors.^18–25 It is crucial to identify cancer-specific miRNAs and further elucidate their role in tumor cell proliferation, invasion, and metastasis through specific targets, in order to identify novel therapeutic targets.^26,27 miR-1225-5p was originally identified as a differentially expressed miRNA in adult renal progenitor cells, and was suggested to regulate the stemness and repair capacity of progenitor cells by targeting CD133, PAX2, and TLR2.²⁸ It was recently reported that frequent downregulation of miR-1225-5p in gastric carcinoma was strongly correlated with an aggressive phenotype and poor prognosis in patients with gastric carcinoma.¹⁷ However, the expression pattern and clinical relevance of miR-1225-5p in human glioblastoma remains largely unknown.

In the present study, we demonstrated that miR-1225-5p was significantly downregulated in six human glioblastoma cell lines and 44 glioblastoma tissue samples, as compared with normal brain tissues. We further analyzed the correlation between miRNA expression and clinicopathological characteristics. A robust negative association between low expression levels of miR-1225-5p and advanced clinical stage was confirmed in patients with glioblastoma. Moreover, we demonstrated that low miR-1225-5p expression in glioblas-toma was associated with poorer overall survival. To the best of our knowledge, this is the first study to provide evidence that miR-1225-5p may play a crucial role in the initiation and development of glioblastoma, and potentially acts as a prognostic biomarker in glioblastoma.

High proliferative, migratory, and invasive abilities are common characteristics of tumor cells.²⁹ To characterize the function of miR-1225-5p in human glioblastoma, we induced stable overexpression of miR-1225-5p in U87 human glioblastoma cells, in which a relatively low expression level of miR-1225-5p was originally found. Our results from the MTT and colony-formation assays demonstrated that ectopic miR-1225-5p overexpression impaired short-term as well as long-term cell proliferation. Further, an in vivo xenograft assay confirmed the inhibitory role of miR-1225-5p in glioblastoma cells. Moreover, overexpression of miR-1225-5p was found to impair the migration and invasion of glioblastoma cells. These data suggest that miR-1225-5p plays a tumor-suppressing role in glioblastoma, similar to its function in gastric carcinoma.¹⁷

The major principle of miRNA function is direct binding to the 3′-UTR of the downstream target of miRNAs, and then inhibition of target gene expression at both the mRNA and protein levels.^9,30 Consistently with previous reports, we found that overexpression of miR-1225-5p was significantly correlated with decreased expression of IRS1, as well as phosphorylation of AKT. As a docking protein, IRS1 plays a key role in transmitting signals from the insulin and insulin-like growth factor-1 (IGF-1) receptors to intracellular pathways, such as the PI3K/Akt pathway.³¹ Previous studies have shown that IRS1 may act as an oncogene; is involved in cancer cell growth, proliferation, migration, invasion, and differentiation;³² and is highly expressed in pancreatic cancer,³³ breast cancer,³⁴ colorectal cancer,³⁵ and glioblastoma.³⁶ Further, our results demonstrated the opposite pattern of expression when compared with miR-1225-5p, as IRS1 expression was commonly elevated in the six human glioblastoma cell lines. More importantly, the effect of miR-1225-5p overexpression on the luciferase activity of pGL3-IRS1-3′-UTR-wt, but not the mutant construct, confirmed that IRS1 is a direct target of miR-1225-5p in glioblastoma. Because it was previously demonstrated that miR-1225-5p inhibits the proliferation and metastasis of gastric carcinoma through IRS1, we consider IRS1 to be the major target of miR-1225-5p in cancer initiation and progression. Given the multiple targets of miR-1225-5p, future research should focus on investigating other potential regulators of glioblastoma progression.

In summary, the present study provided evidence that miR-1225-5p is downregulated in glioblastoma tissues and cell lines, and its low expression is correlated with advanced WHO grade and shorter overall survival of glioblastoma patients. The tumor suppressor effect of miR-1225-5p overexpression in vitro and in vivo makes it interesting to investigate for further development of miR-1225-5p-based therapeutic strategies.

Supplementary materials

Figure S1 (A) The correlation between the miR-1225-5p level (low or high) with overall patient survival was assessed with Kaplan–Meier analysis in patients with grade I/II tumors. (B) The correlation between the miR-1225-5p level (low or high) with overall patient survival was assessed with Kaplan–Meier analysis in patients with grade III/IV tumors.

Figure S2 (A) The expression of miR-1225-5p in SW1088 and SW1088-1225 (stably expressing miR-1225-5p) cells was measured by reverse transcription PCR (RT-PCR) and normalized against U6. (B) SW1088 and SW1088-1225 cells were subjected to an MTT assay to assess cell proliferation at each time point. (C) Representative image of crystal violet-stained cell colonies analyzed by colony-formation assay at 14 days. (D) Quantification plots of (C). **P<0.01; ***P<0.001.

Figure S3 (A) Representative picture of SW1088 and SW1088-1225 cell migration was measured by scratch wound assay. (B) The wound healing percentage was plotted in (A). (C) Representative image of SW1088 and SW1088-1225 cell invasion using Matrigel-based invasion assay. (D) The average number of invading cells of five randomly selected fields from each sample was calculated. **P<0.01 vs SW1088 group.

Figure S4 (A) Western blotting of IRS1 in SW1088 cell transfected with si IRS1. (B) SW1088 si con and SW1088 si IRS1 cells were subjected to an MTT assay to assess cell proliferation at each time point. (C) Representative picture of SW1088 si con and SW1088 si IRS1 cell migration was measured by scratch wound assay. (D) The wound-healing percentage was plotted in (C). (E) Representative image of SW1088 si con and SW1088 si IRS1 cell invasion using Matrigel-based invasion assay. (F) The average number of invading cells of five randomly selected fields from each sample was calculated. **P<0.01 vs SW1088 si con group.

Disclosure

The authors report no conflicts of interest in this work.