Up-regulation of circular RNA hsa_circ_01844 induces apoptosis and suppresses proliferation and migration of glioblastoma cells.

BACKGROUND: Previous studies have demonstrated that various circular RNAs are involved in the malignant proliferation of cancers, such as liver cancer, lung cancer, breast cancer, and others. The potential role of circular RNAs in glioblastoma, however, is still uncertain. In this study, we aimed to study the potential role of hsa_circ_01844 in glioblastoma.
METHODS: Using reverse transcription-polymerase chain reaction (RT-PCR) method, hsa_circ_01844 expression was measured in five glioblastoma samples and five normal brain samples. To evaluate the potential function of hsa_circ_01844 in glioblastoma, hsa_circ_01844 was overexpressed in glioblastoma cell lines (U251 and U87 cells). Using these two cell lines, in vitro experiments including the flow cytometry assay, 3-(4, 5-dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium bromide assay, Transwell assay, and cell apoptosis assay were performed to investigate the role of hsa_circ_01844 in glioblastoma. Student t test and one-way analysis of variance were used for statistical analysis.
RESULTS: The expression of circular RNA hsa_circ_01844 was lower in glioblastoma tissues when compared with the normal brain tissues by RT-PCR method (0.034 ± 0.036 vs. 1.630 ± 0.891, P < 0.001). Using two glioblastoma cell lines, we found that overexpression of hsa_circ_01844 in glioblastoma cells suppressed their proliferation, colony formation, migration, and increased the apoptotic rate compared with empty vector group and blank control group (all P < 0.05).
CONCLUSION: Hsa_circ_01844 shows decreased expression in glioblastoma and its overexpression induces apoptosis and inhibits proliferation, migration, and invasion of glioblastoma cells.

Introduction

Glioma is the most common primary malignancy of the central nervous system. As the most biologically aggressive subtype, glioblastoma accounts for more than 80% of malignant glioma. Glioblastoma patients have a relatively short median survival of 14 to 17 months after surgery and chemoradiotherapy.[ Only 3% to 5% of patients can survive more than 3 years.[ The aggressive nature of glioblastoma lies in the diversity and complexity of its pathogenesis, which is also the main reason for its poor therapeutic effect.[ Exploring the mechanisms of malignant proliferation in glioblastoma cells and how to inhibit their growth is a hot topic in the field of glioma research.

Circular RNAs (circRNAs) are a novel type of RNA molecule formed by a covalently closed loop. Many studies have indicated that circRNAs play key roles in various physiological and pathological processes.[ Recent studies have demonstrated that circRNAs are expressed at variable levels in solid tumors such as esophageal cancer, liver cancer, colon cancer, and others.[ Similar to long non-coding RNAs, circRNAs can also function as a microRNAs (miRNAs) sponge and interact with related miRNAs. For example, studies have shown that there were 70 miRNA-7 (miR-7) receptor binding sites on circRNA-7. The binding of miR-7 to circRNA-7 inhibits the expression of downstream targets of miR-7 and sequentially promotes the malignant proliferation of cervical cancer.[ Some novel circRNAs have been identified to promote or inhibit the proliferation, migration, and invasion of glioma cells.[ Over-expression of has_circ_067934 promotes proliferation and metastasis of malignant gliomas cells by up-regulating phosphatidylinositol 3 kinase/protein kinase B (PI3K-AKT) pathway.[ Exploring the mechanisms of circRNAs in the etiology and pathogenesis of glioblastoma may provide new ideas for the treatment of this disease.[

A previous study using circRNA microarray analysis showed that a novel circRNA, hsa_circ_01844, had the lowest expression in glioblastoma tissues and its expression level was significantly different between glioblastoma and normal brain tissues.[ In this study, we aimed to study the potential role of hsa_circ_01844 in glioblastoma.

Methods

Ethical approval and acquisition of tissue samples

This study was approved by the Ethics Committee of the 904th Hospital of Joint Logistic Support Force PLA (Wuxi Clinical College of Anhui Medical University) (LLSO20200846). All participants signed informed consent.

From May 2019 to November 2019, 5 patients with primary glioblastoma who received surgical treatment in the 904th Hospital of Joint Logistic Support Force PLA were enrolled for this study. After the glioblastoma tissue was removed from the patients, the specimens were frozen in liquid nitrogen and then stored in the refrigerator at −80°C immediately. Five normal brain tissues, which were taken from trauma patients in the 904th Hospital of Joint Logistic Support Force PLA were used as the control. In these control patients, partial resection of the normal brain was required as decompression treatment for their severe head injuries.

RNA extraction and quantitative real-time reverse transcription polymerase chain reaction (RT-PCR)

Trizol reagent was added to the fresh glioblastoma tissues which were obtained intraoperatively, and the total RNA was extracted following the instructions of the reagent manufacturer. The first strand complementary DNA (cDNA) was synthesized by first strand cDNA synthesis kit (Invitrogen K1622, Carlsbad, California, USA). Quantitative real-time RT-PCR was performed using Roche480II Real-Time PCR System (Roche, Switzerland). Each specimen was tested in triplicate. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as a control. The primer sequences were: GAPDH forward primer, 5′-GCGAGATCCCTCCAAAATCAA-3′ and reverse primer, 5′-GTTCACACCCATGACGAACAT-3′; hsa_circ_01844 forward primer, 5′-AGCCACGTTGCTGTTAAACC-3′ and reverse primer, 5′-TAGTCGAGGGCCTTTTCAAG-3′.

Construction of lentivirus vector overexpressing hsa_circ_01844 in glioblastoma cell lines

The lentivirus vector was linearized by restriction enzyme digestion. Hsa_circ_01844 fragments were prepared by PCR amplification and inserted into the linearized lentivirus vector. The plasmids were transfected into bacteria and agar plates were used for bacterial culture and clone section. The positive clones were sequenced, analyzed, and then expanded to obtain high-quality plasmid for downstream virus packaging.

Cell lines and culture

Human glioblastoma cell lines (U87 and U251) were obtained from the American Type Culture Collection. Both U87 and U251 cells were cultured with 10% fetal bovine serum (FBS) and Dulbecco modified Eagle medium (DMEM) in a humidified atmosphere (5% CO2, 37°C).

Flow cytometric analysis

Logarithmic growth cells were inoculated into six-well plates (Corning, New York, USA) at a density of 5 × 105/mL (triplicate). After 48 h of culture, the cells were washed with phosphate-buffered saline (PBS) and fixed with 75% ethanol. Following the manufacturer's instructions, the cell cycle was measured by flow cytometry (BD Biosciences FACS Calibur, BD Biosciences, Piscataway, NJ, USA).

3-(4, 5-Dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium bromide (MTT) assay

Cells in each group were cultured in DMEM with 10% FBS. Cells in logarithmic growth cycle were inoculated into 96-well plates (5000 cells/well) (triplicate) to grow for 1–4 days. Cell proliferation was analyzed following the protocol of MTT kit. The absorption value of each well was measured at 490 nm using a microplate reader (EL × 800, BIO-TEK, Vermont, USA). Data were collected from three independent assays.

Transwell assay

To measure cell migration, Transwell assay was performed. Transwell compartments with 24 holes and pore size of 8-μm (Corning) were coated with or without 150 μL Matrigel. Cell suspension was added in the upper chamber to postpone the logarithmic growth phase of the cells, and 600 μL medium containing 20% FBS was added in the lower chamber. After 24 h of culture, the disposed cells were fixed with 4% paraformaldehyde, stained with crystal violet solution (Sigma-Aldrich, Saint Louis, MO, USA), washed with PBS, and finally removed with cotton swabs.

Cell apoptosis analysis

Cells were cultured in Roswell Park Memorial Institute-1640 medium containing 10% FBS and 1% double antibody (penicillin-streptomycin mixture) for 24 h at a density of 300,000 cells/well. Cells were then collected 48 h after transfection. After washing with PBC, annexin V-fluorescein isothiocyanate apoptosis detection kit (BD Biosciences) was used to analyze cellular apoptosis by flow cytometry.

Statistical analysis

SPSS 20.0 software (IBM, Chicago, IL, USA) was used for statistical analysis. The data were obtained from at least three independent experiments, expressed as mean ± standard deviation. The Student's t test was used to analyze the differences between two groups. One-way analysis of variance was used to analyze the differences among three groups. A value of P < 0.05 was considered statistically significant.

Results

Hsa_circ_01844 was down-regulated in glioblastoma tissues

The expression level of hsa_circ_01844 was significantly lower in glioblastoma tissues when compared to normal control brain tissues (0.034 ± 0.036 vs. 1.630 ± 0.891, P < 0.001).

Overexpression of hsa_circ_01844 inhibited the cell cycle and proliferation of glioblastoma cells

As shown in Figure 1, the content of G2 phase in hsa_circ_01844 overexpression group (lentivirus group) was significantly lower than that in empty vector (EV) and blank control groups (all P < 0.05). This indicated that hsa_circ_01844 inhibited cell cycle progression and proliferation of U87 and U251 cells.

Figure 1

Flow cytometry was performed to monitor cell cycle in U87 and U251 cells. All data were collected from three independent assays. (A) U87 cells: the content of G2 phase in hsa_circ_01844 overexpression group (hsa_circ_01844 lentivirus) was significantly lower than that in EV and blank control (no transfection) groups. (B) U251 cells: the content of G2 phase in the hsa_circ_01844 overexpression group was significantly lower than that in EV and blank control groups. EV: Empty vector.

Overexpression of hsa_circ_01844 inhibited cell growth and reduced colony formation of glioblastoma cells

To determine the impact of hsa_circ_01844 on cell growth, MTT assay was performed using glioma cells transfected with hsa_circ_01844. Cells transfected with EVs and cells with no transfection were included as controls. As shown in Figure 2A and 2B, hsa_circ_01844 overexpression significantly inhibited cell growth (all P < 0.001). Next, we performed colony formation assay to examine the impact of hsa_circ_01844 overexpression on colony formation. As shown in Figure 2C and 2D, cells with hsa_circ_01844 overexpression showed significantly decreased numbers of colonies (all P < 0.01).

Figure 2

Overexpression of hsa_circ_01844 inhibited the growth and colony formation of U87 and U251 cells. All data were collected from 3 independent assays. (A and B) Cell growth was measured using MTT method. The overexpression of hsa_circ_01844 significantly inhibited the growth of glioma cells. (C and D) Colony formation was measured in hsa_circ_01844 overexpressed cells and control cells. The overexpression of hsa_circ_01844 significantly decreased the numbers of colonies. EV: Empty vector; MTT: 3-(4, 5-Dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium bromide.

Up-regulation of hsa_circ_01844 inhibited the migration of glioblastoma cells

To determine whether hsa_circ_01844 has an impact on the migration of glioblastoma cells, transwell assays were conducted using U251 and U87 cells. As shown in Figure 3, overexpression of hsa_circ_01844 led to a significant decrease in the migration ability of glioblastoma cells when compared with control groups (all P < 0.001).

Figure 3

Hsa_circ_01844 inhibited cell migration in U87 and U251 cell lines in Traswell assay. All data were collected from three independent assays. (A) Overexpression of hsa_circ_01844 inhibited the U87 cells’ migration, and the number of migrating cells was significantly lower than that of control groups. (B) Overexpression of hsa_circ_01844 inhibited the U251 cells’ migration, and the number of migrating cells was significantly lower than that of control groups. Crystal violet staining. Original magnification, ×100.

Overexpression of hsa_circ_01844 promoted tumor cell apoptosis in glioblastoma cell lines

We next analyzed the apoptotic function of hsa_circ_01844 in glioblastoma cell lines by flow cytometry. As shown in Figure 4, cells with hsa_circ_01844 overexpression displayed a significant increased apoptotic rate when compared with control groups in both U87 and U251 cells (all P < 0.001). The result suggested that hsa_circ_01844 might accelerate the apoptosis in glioblastoma cell lines.

Figure 4

Overexpression of hsa_circ_01844 induced glioblastoma cells’ apoptosis. All data were collected from three independent assays. (A) Overexpression of hsa_circ_01844 promoted the apoptotic rates of U87 cells, as compared to the control groups. (B) Overexpression of hsa_circ_01844 promoted the apoptosis rate of U251 cells, as compared to the control groups. EV: Empty vector.

Discussion

With the rapid progress of high-throughput sequencing of non-coding RNAs, more attention turns to circRNAs, which are a special class of endogenous RNA, widely exist in human cells and take a great proportion of transcribing genetic information.[ Recent studies have demonstrated that circRNA plays a crucial role in various physiological and pathological processes, including cancers pathogenesis.[ Previous studies have shown that some novel circRNAs are involved in the proliferation of gliomas, such as circRNA ITCH, circRNA FBXW7, circRNA NFIX, and others.[ In our study, we demonstrated that circRNA hsa_circ_01844 was down-regulated in glioblastoma tissues and it regulated the proliferation, migration, invasion, and apoptosis of glioblastoma cells. These findings have clinical implications for both diagnosis and potential therapeutic purposes.

The main mechanism of circRNAs in the role of tumor's pathogenesis may lie in the competitive adsorption of miRNA receptors or self-binding miRNA.[ Compared with their linear isomers, circRNAs have higher expression levels and more miRNA binding sites, thus are more effective in sequestering miRNAs.[ Therefore, the carcinogenic mechanisms driven by circRNAs may occur through their interaction with miRNAs to regulate subsequent pathways. Increasing evidence has supported this theory. Hsa_circ_005211 promotes the migration and invasion processes of breast cancer cells by acting as a sponge for mir-125a-5p.[ Similarly, hsa_circ_0014359 acts as a miRNA-153 sponge in glioma cells and promotes the malignant proliferation of glioma cells by targeting miRNA-153/PI3K signaling pathway.[

In our study, we showed hsa_circ_01844 can potentially act as an inhibitor in glioblastomas by inducing apoptosis and inhibiting proliferation, migration of U251 and U87 cells. We found that the content of G2 phase in hsa-circ-01844 overexpression group (hsa_circ_01844 lentivirus) was significantly lower than that in control groups, which indicated that hsa_circ_01844 inhibit the U87 and U251 cells’ cell cycle progression. Furthermore, based on bioinformatics prediction, hsa_circ_01844 may interact with several miRNAs, such as miRNA-616 and miRNA-671-5p. MiRNA-616 might act as a tumor promoter and play oncogenic roles in glioma by regulating the Sry related HMG box-7 and wingless/integrated/β-catenin signaling,[ and miR-671-5p can regulate complementarity determining region 1 pathway to inhibit the development of glioblastoma.[ Accordingly, it was proposed that hsa_circ_01844 can sponge tumor promoter miRNAs and then inhibit the expression level of some targeted oncogenic genes, which can promote the glioblastoma progression.

Collectively, the present study investigates the expression level and function of has_circ_01844 in glioblastoma. Has_circ_01844 might serve as a novel therapeutic target for this malignant tumor in the future.

Conflicts of interest

None.