- Primary research
- Open Access
CircRNA NALCN acts as an miR-493-3p sponge to regulate PTEN expression and inhibit glioma progression
Cancer Cell International volume 21, Article number: 307 (2021)
An increasing number of studies have shown that circular RNAs (circRNAs) play important roles in the regulation of tumor progression. Therefore, we explored the expression characteristics, function, and related mechanism of the newly identified circNALCN in glioma.
RNA sequencing was used to analyze the expression profiles of circRNAs in brain tissue from five glioma cases and four normal controls. Quantitative real-time polymerase chain reaction was implemented to examine the levels of circNALCN, miR-493-3p, and phosphatase and tensin homolog (PTEN). Cell counting kit 8 assays were performed to analyze cell proliferation, and cell migration was assessed by the wound healing test and Transwell assay. Dual-luciferase reporter, fluorescence in situ hybridization, and RNA pulldown assays were performed to confirm the role of circNALCN as an miR-493-3p sponge, weakening the inhibitory effect of miR-493-3p on target PTEN expression.
The downregulated expression of circNALCN was observed in both glioma tissues and cell lines. CircNALCN expression was negatively correlated with World Health Organization grade and overall survival in patients with glioma. Functionally, the overexpression of circNALCN significantly inhibited the proliferation and migration of glioma cells, whereas miR-493-3p mimics counteracted these effects. The mechanistic analysis demonstrated that circNALCN acted as a competing endogenous RNA for miR-493-3p to relieve the repressive effects of miR-493-3p on its target, PTEN, suppressing glioma tumorigenesis.
CircNALCN inhibits the progression of glioma through the miR-493-3p/PTEN axis, providing a developable biomarker and therapeutic target for glioma patients.
Gliomas account for approximately 30% of all intracranial tumors and are the most common primary brain tumors . Unfortunately, most gliomas (approximately 80%) are malignant in nature, with a worldwide incidence of 3.9 per 100,000 in men and 3.1 per 100,000 in women . Currently, comprehensive treatment methods, including surgery, radiotherapy, chemotherapy, and immunotherapy, are the primary options applied in the clinic. However, due to the unclear pathogenesis of glioma and the lack of specific treatment targets, the overall effect is often unsatisfactory, resulting in a median survival time of only 15–16 months among patients with World Health Organization (WHO) Grade IV glioblastoma . In recent years, gene therapy for glioma has been widely studied, but no specific targets or mechanisms have been clearly identified. Therefore, the treatment and prognosis of glioma are likely to benefit greatly from explorations of the underlying molecular mechanism of glioma.
Circular RNAs (circRNAs) are a newly discovered class of endogenous, non-coding RNA found in eukaryotes. CircRNAs are typically comprised of exons that form covalently closed, continuous loops and contain typical splicing sites [4, 5]. Traditionally, circRNAs were considered to be meaningless by-products of gene rearrangements and mRNA splicing . However, due to advances in high-throughput sequencing and bioinformatics analyses, emerging studies have demonstrated that circRNAs are evolutionarily conserved, stable, and richly expressed, with specific expression patterns across different cells and developmental stages . In addition, circRNAs have been associated with a variety of functions, including acting as a microRNA (miRNA) sponges, binding to RNA binding proteins (RBPs), and participating in protein translation [8, 9].
By analyzing non-coding RNAs expressed in 29 different stages of nerves and tissue, AgnieszkaRybak-Wolf identified the rich and specific expression of circRNAs throughout the brain, which showed a high degree of stability and incompatibility compared with related linear mRNAs , indicating that circRNAs may be involved in neurodifferentiation processes and the development of neurological diseases. However, the abundance of circRNAs found in glioma specimens is lower than that in normal samples, which prompted us to consider the role played by circRNAs in glioma progression and potential clinical applications . An in-depth understanding of the role and mechanism of circRNAs in glioma will be helpful for developing new detection methods and effective treatments.
MicroRNAs (miRNAs) are a type of ubiquitous and conserved non-coding small RNA, which act as negative gene regulators that inhibit the expression of target genes . They are involved in a wide range of biological processes associated with cancer, including proliferation, invasion, differentiation, and apoptosis . Increasing evidence has suggested that circRNAs play important roles in tumorigenesis and cancer development by regulating miRNAs . For example, circHIPK3 can bind to multiple miRNAs, including miR-654 and miR-124-3p, and several studies have demonstrated that circHIPK3 targets miR-654, which targets insulin-like growth factor 2 mRNA binding protein 3 (IGF2BP3), suggesting that circHIPK3 promotes glioma progression by regulating the miR-654/IGF2BP3 axis . Another study showed that circ-HIPK3 promotes the proliferation and invasion of glioma cells by acting as a sponge for miR-124-3p, resulting in the upregulation of signal transducer and activator of transcription 3 (STAT3) . However, there may be other circRNAs that are regulated by the ceRNA network in the occurrence and development of gliomas have not been discovered.
Several studies have shown that phosphatase and tensin homolog (PTEN) is a multifunctional tumor suppressor that contains a catalytic domain and a tensin-like domain [17, 18]. PTEN plays an important role in the regulation of cell proliferation, apoptosis, and tumor invasion. The results of clinical studies examining high-grade gliomas have shown that changes in PTEN gene expression are associated with poor prognosis and may affect the response to specific treatments [19, 20]. One mechanism involved in the regulation of PTEN expression is miRNA-mediated regulation [21, 22].
Based on RNA sequencing (RNA-seq) results, we first identified a novel glioma-related circRNA, termed circNALCN (circBase ID: hsa_circ_0099761), which is derived from exons 2 to 7 of the sodium leak channel, non-selective (NALCN) gene. We found that circNALCN has not yet mentioned its role in tumors (including glioma) and other diseases. To obtain insights into the function and underlying molecular mechanisms through which circNALCN participates in the development and progression of glioma, the clinical significance of circNALCN was explored in glioma tissues. CircNALCN was found to be downregulated in glioma tissues and acted as a sponge for miR-493-3p, affecting the expression of PTEN to modulate glioma tumorigenesis and progression. Collectively, our findings suggest that circNALCN can be used as a prognostic predictor and therapeutic target in glioma.
Sequencing data analysis
High-throughput sequencing was performed on five cases of glioblastoma and four cases of normal brain tissue, and the samples were processed and sequenced by Lc-Bio Technologie (HangZhou, China). The results were processed with the edgeR software package (https://www.Bioconductor.org/). Significant differences in circRNAs and miRNAs expression were screened to identify those with log2 (fold-change) ≥ 2 and P < 0.05. The results showed significant differences between circRNA and miRNA expression levels between glioblastoma and normal cases.
Patient tissue samples
In our study, 76 cases of glioma and 17 normal tissues were collected from August 2018 to May 2020 at Xiangya Hospital. All patients underwent surgical resection of glioma or non-tumor operations, such as craniocerebral trauma, which were performed in accordance with the Helsinki Declaration. After surgical resection, all specimens were frozen in liquid nitrogen and stored at − 80 °C until RNA was extracted. The study was approved by the Ethics Committee of Xiangya Hospital of Central South University, and informed consent was obtained from all patients.
Human glioma SF126 cells were obtained from the Cell Research Institute of Peking Union Medical College (Peking, China). Human U251 and U87 glioma cells were obtained from the American Type Culture Collection (ATCC). The normal glial cell line HEB was purchased from the Cell Bank of the Type Culture Collection of the Chinese Academy of Sciences (Shanghai, China). The cells were cultured in Dulbecco’s modified Eagle medium (DMEM) containing 10% fetal bovine serum (FBS) and 100 U/mL penicillin/streptomycin (all from Gibco, Carlsbad, CA, USA) at 37 °C in a humidified incubator containing 5% CO2.
RNA extraction and quantitative real-time polymerase reaction (qRT-PCR)
Total RNA was extracted from glioma and normal tissues and cell lines using TRIzol Reagent (Takara, Japan), and genomic DNA (gDNA) was isolated with FastPure DNA Isolation (Vazyme, China), according to the manufacturer’s instructions. The concentration and purity of RNA samples were measured using a Nanodrop 2000 (Thermo Fisher Scientific, USA). For circRNA and mRNA, reverse transcriptions were performed using the PrimeScript RT Master Mix (Takara, Japan) with random primers. For miRNA, reverse transcriptions were performed using the PrimeScript RT Reagent Kit (Takara, Japan) with specific stem-loop primers. cDNA amplification was performed using Hieff® qPCR SYBR Green Master Mix (Low Rox Plus, YeSen, Shanghai, China) using an ABI Prism 7500 sequence detection system (Applied Biosystems, USA). PCR amplification conditions were as follows: initial denaturation for 30 s at 95 °C, and 40 cycles of 10 s at 95 °C, followed by 30 s at 60 °C. Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) and U6 were used as internal controls, and each sample was measured three times. The relative quantification of circRNA, miRNA, and mRNA expression levels were compared to the expression level of the internal control and analyzed using the 2−ΔΔCT method. Divergent primers were used to detect the back splice junction of circRNAs, and convergent primers were used to detect linear mRNA. All primers used in this study are listed in Table 1.
Nucleic acid electrophoresis
cDNA and gDNA amplification products were detected by 2% agarose gel electrophoresis in tris-acetate-EDTA (TAE) buffer. DNA was separated by 120 V electrophoresis for 25 min. DNA sizes were compared against Marker L (50–5000 bp) (Sango Biotech, China). The bands were detected by ultraviolet light.
RNase R treatment
The RNAs (1 µg) from U87 and U251 cells were treated with RNase R (2 U/µg, GeneSeed, Guangzhou, China) and incubated at 37 °C for 30 min, followed by inactivation by 10 min at 70 °C. The treated RNAs were reverse transcribed with divergent primers or convergent primers and detected by qRT-PCR assay.
CircNALCN small interfering RNA (si-circNALCN), miR-493-3p mimics, miR-493-3p inhibitor, negative control (NC) mimics, and NC inhibitor were all designed and synthesized by RiboBio (Guangzhou, China). The circNALCN overexpression plasmid (pcDNA3.1-circNALCN), PTEN overexpression plasmid (pcDNA3.1-PTEN), and pcDNA3.1-NC were synthesized by Genscript (Nanjing, China). U87 and U251 cells were inoculated in 6-well plates and cultured to 60–70% confluence, and si-circNALCN, miR-493-3p mimic, miR-493-3p inhibitor and control RNA (Ribio, Guangzhou, China) were transfected into U87 and U251 cells with a final concentration of 20 uM Lipofectamine® 2000 (Invitrogen; Thermo Fisher Scientific, Inc.) for continuous 6 h. Similarly, 2000 ng pcDNA3.1-circNALCN, pcDNA3.1-PTEN and pcDNA3.1-NC were also transfected into U87 and U251 cells. The transfected cells were incubated at 37 °C and 5% CO2 for 24 h, and the transfection efficiency was detected by reverse transcription-polymerase chain reaction (RT-PCR).
Luciferase reporter assay
The sequences of circNALCN and PTEN-3ʹUTR and their corresponding mutation were designed, synthesized and inserted into luciferase reporter vector pmirGLO (Promega, USA), termed circNALCN-WT, circNALCN-MUT, PTEN-WT and PTEN-MUT, respectively. All luciferase reporter vectors were synthesized by GenScript (Nanjing, China). U87 and U251 cells were inoculated in a 96-well plate to achieve 70% confluence. The 20 uM miR-493-3p mimics and NC mimics were cotransfected with 100 ng circNALCN-WT, circNALCN-MUT, PTEN-WT, and PTEN-MUT into U87 and U251 with lipofectamine 2000. Then, the relative luciferase activity was examined using a Dual-Luciferase Assay Kit (Promega, USA), according to the manufacturer’s protocol.
Fluorescence in situ hybridization (FISH)
Fluorescence in situ hybridization (FISH) assay was performed to observe the locations of circNALCN and miR-493-3p expression in U87 and U251 cells. Fluoresceine amidite (FAM)-labeled circNALCN probes and cyanine-3 (Cy3)-labeled miR-493-3p were designed and synthesized by Servicebio (Wuhan, China). Hybridization was performed overnight using the circNALCN and miR-493-3p probes, according to the manufacturer’s instructions. The images were acquired on a Zeiss LSM710 Laser-Scanning Confocal Microscope (Zeiss Instrument Inc., Germany). The sequences of circNALCN and miR-493-3p that were probed during FISH are listed in Additional file 1: Table S1.
RNA pulldown assays
U251 and U87 cells (4 × 107) were collected, fixed with 1% formaldehyde, and add cell lysate (BersinBio, Guangzhou, China) on ice. RiboBio (Guangzhou, China) designed and synthesized the circNALCN probe, PTEN probe, and NC biotin probe. RNA Antisense Purification (RAP) kit was purchased from bersinbio (Guangzhou, China). According to the manufacturer’s instructions, the circNALCN, PTEN and NC biotin probes were added to the cell lysate, hybridized at 37 °C for 30 min, denatured at 50 °C for 5 min, hybridized again at 37 °C for 120 min, and incubated with streptomyces and magnetic beads (BersinBio, Guangzhou, China) at 25 °C for 30 min. RNA complexes bound to magnetic biotin beads were eluted using RNA elution buffer, and qRT-PCR was performed to analyze the relative expression levels of miR-493-3p.
Cell counting kit-8 proliferation assay
U87 and U251 cells were plated at 8 × 103 cells per well of a 96-well plate. A 10 µl volume of the CCK-8 reagent (Yisheng, Shanghai, China) was directly added to the culture medium at specified times (0, 24, 48, 72 and 96 h). The cells were incubated with the CCK-8 reagent for 1 h at 37 °C, and the optional density (OD) was measured at 450 nm using a miniature tablet reader (BioTek Instruments, USA).
Wound healing assays
U87 and U251 cells were cultured in a 6-well plate until a cell monolayer was formed. The monolayer was subsequently scratched using a 10 µl pipette tip. Cell migration images were taken at 0 and 36 h after injury in four high-power fields. After 36 h of measurement, the wound distance was reduced to the standard, normalized against the 0 h control, and expressed as relative mobility.
Transwell migration assays
Transwell analysis was performed as described by the manufacturer. U87 and U251 cells were inoculated in the upper compartment in 100 µl of serum-free medium. The bottom chamber of the Transwell chamber (Corning, USA) was filled with 600 µl medium and 10% FBS as a glioma cell chemoattractant. After incubation for 24 h, the upper chamber was fixed and stained with crystal violet (Kaigen, China) for 15 min. Images were taken from five different fields of vision, and the cells were counted.
Cells were lysed in radioimmunoprecipitation assay lysis buffer (RIPA, Beyotime, China). The protein was prepared and quantified using bicinchoninic acid (BCA) analysis (Beyotime, China). Equal amounts of protein were separated by 10% sodium dodecyl sulfide-polyacrylamide gel electrophoresis (SDS-PAGE) and transferred onto polyvinylidene fluoride (PVDF) membranes (Millipore, MA, USA). Membranes were blocked with 5% skim milk powder and incubated with the primary antibodies anti-PTEN (1:1000, Bioworld Technology Inc.) and anti-GAPDH (1:10,000, Bioworld Technology Inc.) at 4 °C overnight. The membranes were incubated with a secondary antibody (1:10,000, Bioworld Technology Inc.) for 1 h in room temperature. Finally, the proteins were visualized using enhanced chemiluminescence (ECL) reagent (Yisheng, Shanghai, China), and the relevant data were analyzed by ImageLab software.
Statistical analysis was performed using GraphPad Prism 8, and Student’s t-test and one-way analysis of variance (ANOVA) were used to determine significant differences between two or more groups. Pearson’s correlation coefficient was used to perform correlation analyses. The results are expressed as the mean ± standard deviation (mean ± SD). The overall survival (OS) curve was calculated by Kaplan–Meier analysis, and the log-rank test was performed. Differences were considered significant at P < 0.05.
circNALCN (hsa_circ_0099761) is expressed at low levels in gliomas
We analyzed five cases of glioblastoma tissues and four cases of normal brain tissues by high-throughput RNA-seq. A total of 20,204 circRNAs were sequenced using edgeR packets to perform difference analyses of glioma tissue and normal brain tissue. Based on the condition of |log2 (fold change(FC)) > = 2| and P < 0.05, 652 differentially expressed circRNAs (DEcircRNAs), including 67 upregulated DEcircRNAs and 585 downregulated DEcircRNAs were identified. In cancers with high cell proliferation rates, circRNA proliferation has been reported to be diluted and downregulated . Therefore, we focused on the significantly downregulated circRNAs. Using the resulting volcano plot, we identified a significantly downregulated circRNA (hsa_circ_0099761), named circNALCN, which is transcribed from NALCN on human chromosome 13 through the head–tail splicing of exons 2–7 (Fig. 1A, B). And then we found that circNALCN is a novel circRNA, that has not been reported in tumors and other diseases, and its function is not clear, so we studied the function of circNALCN in glioma cells. We verified that the available sequencing data for circNALCN and our RNA-seq data were essentially the same. Compared with normal brain tissue, circNALCN decreased in glioma tissue (Fig. 1C). To further verify whether circNALCN expression was downregulated in glioma tissues and cells, we used qRT-PCR analysis to examine circNALCN expression in 76 cases of glioma tissue and 17 cases of normal brain tissue. We found that the expression of circNALCN was significantly decreased in glioma tissues compared with the expression level in normal tissues (Fig. 1D and Table 2, P < 0.001). We also found that the expression levels of circNALCN in WHO Grade I–II gliomas were higher than that in WHO Grade III–IV gliomas (Fig. 1E, P < 0.001), indicating that circNALCN expression may be related to glioma tumor grade.
Next, we compared circNALCN expression levels between HEB, U87, U251 and SF126 cells lines. Compared with the normal glial cell line (HEB), the expression of circNALCN is lower in glioma cells. No significant differences were observed for the expression levels of circNALCN among the various glioma cell lines (Fig. 1F). Therefore, we selected the U87 and U251 cell lines to study the downstream regulatory pathway associated with circNALCN.
Characteristics and clinical features of circNALCN
CircNALCN originates from NALCN, located on chromosome 13. CircNALCN is formed by the splicing of exons 2–7 (101997616–102051516) from head to tail (Fig. 1B). However, head-to-tail splicing may be the result not only of trans-splicing but also of genome rearrangement. To eliminate these two possibilities, we designed special divergent and convergent primers to differentiate between circNALCN generated from cDNA and gDNA. DNA was extracted from U251 cells and analyzed by PCR and agarose gel electrophoresis. The results showed that circNALCN was only detected as cDNA, whereas no product was detected in the extracted gDNA (Fig. 2A). Stability is considered to be an important feature of circRNA . The stability of circNALCN was verified by ribonuclease R (RNase R) treatment. cDNA extracted from RNase R-treated U251 and U87 cells was analyzed by PCR using specially designed divergent and convergent primers, which revealed that circNALCN could resist restriction endonuclease digestion by RNase R, whereas linear NALCN was digested by RNase R. These results suggested that circNALCN is not the result of genome rearrangement and represents a stable product. The results of the FISH analysis revealed that although circNALCN was expressed in the nucleus and cytoplasm to a certain extent, it was mainly located in the cytoplasm (Fig. 2D).
When we examined the clinical data of the examined patients, we found that the expression level of circNALCN was significantly correlated with the WHO glioma grade but not with age or sex (Fig. 1E, Table 2). In addition, the OS curve, generated using the Kaplan–Meier method, indicated that patients with low expression levels of circNALCN in glioma tissues had significantly reduced survival compared with that of patients with high levels of circNALCN in glioma tissues (P = 0.0009, Fig. 2E). These results suggested that circNALCN is a stable circRNA that may serve as a meaningful prognostic marker in glioma and warrants further study.
CircNALCN inhibits the proliferation and migration of glioma cells
To explore the biological functions of circNALCN in glioma cells, three siRNAs were generated to target circNALCN, and a circNALCN overexpression vector was constructed. The circular transcript expression vector circNALCN was successfully constructed in U87 and U251 cells (Fig. 3A). Among the three siRNAs, si-circNALCN-2 was found to have the highest inhibition efficiency; therefore, it was selected for use in the following experiments (Fig. 3B).
CCK-8 assay revealed that circNALCN downregulation significantly enhanced cell proliferation in U87 and U251 cells, whereas circNALCN upregulation had the opposite effect. These results showed that circNALCN overexpression could inhibit the proliferation of glioma cells (Fig. 3C, D).
The effects of circNALCN on glioma cell migration was determined by the performance of the wound healing assay and the Transwell experiment. The results showed that the downregulation of circNALCN significantly enhanced the migration ability of glioma cells, whereas circNALCN upregulation significantly inhibited the migration ability of glioma cells (Fig. 3E–H). These experiments showed that circNALCN inhibited glioma cells migration.
High expression levels of miR-493-3p in gliomas are related to poor prognosis
According to the theory of competitive endogenous RNAs (ceRNAs), circRNAs act as miRNA sponges to regulate miRNA expression . Because circNALCN is primarily located in the cytoplasm and shows significant stability, we further explored whether circNALCN suppresses the biological behavior of gliomas by absorbing miRNAs. We used miRNA target prediction tools, including circBank, starBase, and sequencing up-regulated miRNA data, to predict potential targets for circNALCN (Fig. 4A). We sequenced a total of 1138 miRNAs and used edgR packet forward to perform differential expression analyses of glioma tissue and normal brain tissue. Based on the conditions of |log2(FC) > = 2| and P < 0.05, 130 differentially expressed miRNAs (DEmiRNAs) were identified, including 49 upregulated DEmiRNAs, and 81 down-regulated DEmiRNAs. The results of three prediction methods were combined with differences in miRNA levels obtained from the sequencing results (Fig. 4B), which identified miR-493-3P as significantly upregulated in glioma tissues. The expression level of miR-493-3P was also found to be upregulated in HEB cells compared with the three glioma cell lines. Therefore, we selected miR-493-3p for further verification (Fig. 4C). We transfected si-circNALCN and circNALCN overexpression vectors into U87 and U251 cells and observed the effects on miR-493-3p expression. The results showed that miR-493-3p expression was significantly increased in the si-circNALCN group and decreased in the circNALCN overexpression group compared with untreated cells (Fig. 4D, E). The expression levels of miR-493-3p in glioma and normal tissues were detected by qRT-PCR, which showed that miR-493-3p expression was significantly increased in glioma tissues compared with that in normal tissues (Fig. 4F). We found a negative correlation between the expression levels of circNALCN and miR-493-3p in both glioma tissues and normal tissues (Fig. 4G). These results suggested that circNALCN might serve as an miR-493-3p sponge. In addition, we analyzed the relationship between miR-493-3p expression levels and the clinical characteristics of patients with gliomas, as shown in Table 2. The expression of miR-493-3p was correlated with WHO tumor grade and age. The Kaplan–Meier survival curve showed that the expression of miR-493-3p was negatively correlated with the OS among glioma patients (Fig. 4H), suggesting that miR-493-3p may serve as a prognostic indicator.
CircNALCN acts as an miR-493-4p sponge in glioma
To study whether circNALCN serves as an miR-493-3p sponge, double-luciferase reporter gene analysis was applied to glioma cells. The full-length circNALCN-WT sequence and a mutant sequence without the miR-493-3p binding site (circNALCN-MUT) were subcloned into the pmirGLO luciferase reporter vector (Fig. 5A). The results indicated that the luciferase activity of the circNALCN-WT group significantly decreased in the miR-493-3p mimic group compared with the NC mimic group in U87 and U251 cells (Fig. 5B, C). No significant difference in luciferase activity was observed between the miR-493-3p mimics and NC mimic group in cells transfected with the circNALCN-MUT vector. In addition, FISH assays were performed in U87 and U251 cells to observe the subcellular localization of circNALCN and miR-493-3p. In our study, circNALCN (green) and miR-493-3p (red) were both primarily located in the cytoplasm, and circNALCN and miR-493-3p colocalized, suggesting that they may interact (Fig. 5D). To further confirm the direct interaction between circNALCN and miR-493-3p, we applied an RNA pulldown assay using a biotin-labeled circNALCN probe. In both U87 and U251 cells, miR-493-3p was significantly enhanced upon circNALCN capture (approximately fourfold increase in U87 cells and a twofold increase in U251 cells; Fig. 5E).
MiR-493-3p reverses the antitumor effects of circNALCN on glioma
Next, we designed a rescue experiment to explore whether circNALCN inhibits the malignant behavior of gliomas by interacting with miR-493-3p. U87 and U251 cells were cotransfected with si-circNALCN or circNALCN overexpression vector and miR-493-3p inhibitor or mimics. CCK-8 experiments showed that miR-493-3p mimics reversed the decrease in proliferation rate caused by the overexpression of circNALCN, whereas the miR-493-3p inhibitor attenuated the increase in the proliferation rate of glioma cells induced by circNALCN depletion (Fig. 5F). In addition, the Transwell migration assay and wound healing analysis showed that the miR-493-3p mimics group reversed the inhibition of glioma cell migration caused by circNALCN overexpression, whereas the miR-493-3p inhibitor group reversed the promotion of glioma cell migration caused by circNALCN depletion (Fig. 5G, H). These findings suggested that circNALCN inhibited the progression of gliomas by sponging miR-493-3p.
The miR-493-3p targeting of PTEN promotes the proliferation and migration of glioma cells in vitro
Few studies have explored the role played by miR-493-3p in glioma. Therefore, we aimed to clarify the mechanism of action and biological functions of miR-493-3p in glioma cells. According to the prediction results from starBase V3, PTEN was suggested as a potential target gene of miR-493-3p. PTEN has been widely reported to act as a tumor suppressor [17, 18] and has been shown to inhibit tumor progression in gliomas . Therefore, we conducted the next experiment to prove that PTEN binds to miR-493-3p. We constructed a PTEN overexpression vector, which was successfully expressed in U87 and U251 cells (Fig. 6A). To explore whether miR-493-3p regulates the expression of PTEN, miR-493-3p mimics and inhibitor were separately transfected into U87 and U251 cells. The results showed that the expression of PTEN decreased in both cell lines expressing the miR-493-3p, whereas PTEN expression increased in the miR-493-3p inhibitor groups (Fig. 6B, C), which suggested that miR-493-3p can affect the expression of PTEN in glioma cells.
PTEN-WT and PTEN-MUT were further cloned into the pmirGLO luciferase reporter vector and cotransfected in U87 and U251 cells, together with miR-493-3p mimics or NC mimics, to determine the interaction between miR-493-3p and PTEN (Fig. 6D). The results showed that the luciferase reporter activity in both cell lines decreased significantly in the PTEN-WT and miR-493-3p mimics group; however, these effects disappeared when the PTEN binding site was mutated. There were showed that PTEN binds to miR-493-3p (Fig. 6E, F). To confirm the predicted interaction, we performed a pulldown analysis using a biotinylated PTEN probe. In U87 and U251 cells, miR-493-3p was significantly enhanced after PTEN capture (in U87 cells, 1000-fold increase; in U251 cells, twofold increase; Fig. 6G).
In addition, we cotransfected U87 and U251 cells with PTEN overexpression vector or NC vector and miR-493-3p mimics to explore the biological effects of the miR-493-3p targeting of PTEN on glioma cells. The results of qRT-PCR and western blot analyses showed that the expression of PTEN in U87 and U251 cells was significantly decreased in the miR-493-3p mimics and NC vector groups, whereas the expression of PTEN in the miR-493-3p mimics and PTEN overexpression vector group was significantly higher than that in the miR-493-3p mimics and NC vector group (Fig. 6H, I).
We further explored the biological function of the miR-493-3p targeting of PTEN in glioma cells. The CCK8 growth curve showed that miR-493-3p mimics and NC vector group showed significantly increased cell proliferation activity, whereas the cell proliferation of the miR-493-3p mimics and PTEN overexpression vector group was lower than that of the miR-493-3p mimics and NC vector group (Fig. 6J, K). The effects of miR-493-3p on the migration of glioma cells were detected by wound healing and Transwell experiments. The results showed that the miR-493-3p mimics and NC vector group were characterized by the enhanced migration abilities of U87 and U251 cells, whereas the migration abilities of U87 and U251 cells decreased in the miR-493-3p mimics and PTEN overexpression group (Fig. 6L–O). These results suggested that miR-493-3p can promote the migration and proliferation of gliomas through the regulation of PTEN.
CircNALCN acts as an miRNA sponge of miR-493-3p to regulate the expression of PTEN
To further explore the interaction among circNALCN, miR-493-3p and PTEN, we performed a series of experiments. First, we detected the expression of PTEN in glioma and normal tissues. The qRT-PCR results showed that the expression of PTEN in glioma tissues was significantly lower than that in normal tissues (Fig. 7A), and the Kaplan–Meier curve showed that the expression level of PTEN was positively correlated with OS among glioma patients (Fig. 7B), which was consistent with the results of previous research . We found a positive correlation between the expression of circNALCN and PTEN in both gliomas and normal tissues (Fig. 7C). In addition, qRT-PCR and western blot assays demonstrated that the overexpression of circNALCN significantly increased the mRNA and protein levels of PTEN, whereas the co-transfection of the circNALCN overexpression vector with miR-493-3p mimics appeared to counteract these effects in U87 and U251 cells (Fig. 7D). Similarly, the knockdown of circNALCN significantly decreased the mRNA and protein levels of PTEN, whereas the co-transfection of si-circNALCN and miR-493-3p inhibitor eliminated these effects in U87 and U251 cells (Fig. 7E). These data indicated that circNALCN acts as a ceRNA for miR-493-3p to regulate the expression of PTEN, inhibiting the occurrence and development of gliomas.
Increasing evidence regarding the biogenesis and functions of circRNAs has provided a new perspective for understanding the pathogenesis of a variety of diseases, including cancer, cardiovascular diseases, and autoimmune diseases [27,28,29]. Some circRNAs are thought to play oncogenic or tumor-suppressive roles in various human cancers, such as colon cancer, lung cancer, hepatocellular carcinoma (HCC), and gastric cancer . Moreover, the mechanisms of certain glioma-related circRNAs have gradually been revealed . However, this study report the role played by the circNALCN/miR-493-3p/PTEN pathway in the occurrence and development of gliomas.
First, we identified a new circular RNA, named circNALCN, through a high-throughput RNA-seq analysis. CircNALCN was found to be significantly downregulated in gliomas and significantly correlated with tumor grade and OS in glioma patients. Second, through circRNA functional tests, we confirmed that circNALCN overexpression could inhibit the proliferation and migration of glioma cells. Third, the binding of circNALCN to miR-493-3p was predicted by bioinformatics analysis and verified by double-luciferase reporting assay, FISH analysis, and pulldown experiments. We also found that miR-493-3p was significantly upregulated in glioma cells and tissues and correlated with OS and tumor grade in patients. Fourth, mechanistic experiments demonstrated that circNALCN acts as an miR-493-3p sponge to counteract the miR-493-3p-mediated inhibition of the target gene PTEN during glioma progression. Together, these data suggested that circNALCN may serve as a new biomarker and therapeutic target, with an important role in glioma progression and development.
Most circRNAs contain miRNA response elements (MREs), and studies have increasingly shown that the ceRNA mechanism represents the primary mechanism through which circRNAs affect biological functions . Using circBank, starBase, and high-throughput sequencing prediction analysis, circNALCN was predicted to contain an miRNA binding sites for miR-493-3p. Therefore, FISH, biotin-labeled probe pulldown, and double-luciferase reporter gene assays were used to confirm that circNALCN and miR-493-3p colocalized in the cytoplasm of glioma cells and bound directly. Subsequent rescue experiments further confirmed that miR-493-3p could reverse the antitumor effects of circNALCN. Therefore, our results provide evidence that circNALCN acts as a sponge for miR-493-3p, providing a developable biomarker and therapeutic target for glioma patients.
PTEN was identified as a target gene of miR-493-3p by bioinformatics analysis, which was confirmed by luciferase reporter gene analysis. In gliomas, miR-492-3p negatively regulates the expression of PTEN, whereas circNALCN positively regulates the expression of PTEN. Previous studies have shown that PTEN acts as a tumor suppressor in most human tumors, such as pancreatic cancer, HCC, and breast cancer [32,33,34]. Decreased PTEN expression is typically associated with aggressive tumors and a worse prognosis. An abnormal decrease in PTEN expression is considered to serve as a predictor of malignant progression in gliomas, and PTEN can be used as a marker to predict the prognosis of patients with gliomas . Previous studies have found that PTEN participates in the negative regulation of phosphatidylinositol-triphosphate (PIP3) and plays a tumor inhibitory role by negatively regulating AKT signal transduction . Recent studies have shown that circFBXW7 plays a tumor inhibitory role by reversing the expression of NEK2 and mammalian target of rapamycin (mTOR), promoting the expression of PTEN in tumors . Our results showed that the expression of PTEN is low in gliomas and is associated with the OS of patients with gliomas. To date, whether and how circNALCN is involved in the progression of gliomas induced by PTEN remains unclear. Our results in this study showed that the interaction between circNALCN and miR-493-3p promotes the expression of PTEN, suggesting a potential mechanism through which the circNALCN/miR-493-3p/PTEN axis regulates the progression of glioma.
We studied the expression, function and clinical significance of circNALCN in glioma and researched on the relationship between miR-493-3p and PTEN. These findings may have implications for the treatment of gliomas. However, the interpretation of our research results has several limitations. Firstly, the glioma tissue used for the RNA-seq analysis in this study was taken from a homologous population at a hospital, and we cannot rule out the possibility that other important downregulated circRNAs may also be involved in the progression of PTEN-induced gliomas. Secondly, our preliminary study showed that the expression of circNALCN in glioma tissue was significantly downregulated, and up-regulation of circNALCN can inhibit the proliferation and migration of glioma cells. However, we only performed functional tests of circNALCN in vitro. Additional animal experiments are being performed to further explore the tumor inhibitory effect of circNALCN in vivo. Thirdly, our study confirmed the ability of circNALCN to combine with miR-493-3p. However, other miRNAs may also bind circNALCN to regulate the progress and development of gliomas. Fourth, whether circNALCN regulates the progression of gliomas through other mechanisms, such as interactions with RBPs or the sponge absorption of trans-acting elements, requires further study. Therefore, a more in-depth understanding of the therapeutic potential of circNALCN in gliomas remains necessary.
In summary, our results provide evidence to support the downregulation of circNALCN expression in gliomas, which may play a role as a tumor suppressor and prognostic biomarker in gliomas. In addition, we also indicated that circNALCN inhibits the proliferation and migration of glioma cells by downregulating miR-493-3p and upregulating the expression of PTEN. Our results demonstrate the potential mechanism through which circNALCN may regulate the progression of glioma cells through the circNALCN/miR-493-3p/PTEN axis, providing a new potential therapeutic target for glioma patients.
Availability of data and materials
The data sets used and analyzed during the current study are available from the corresponding author on reasonable request.
World Health Organization
Phosphatase and tensin homolog
Fetal bovine serum
Quantitative real-time polymerase reaction
Glyceraldehyde 3-phosphate dehydrogenase
CircNALCN small interfering RNA
RNA binding proteins
Insulin-like growth factor 2 mRNA binding protein 3
Signal transducer and activator of transcription 3
Dulbecco’s modified Eagle medium
Fluorescence in situ hybridization
RNA antisense purification
Differentially expressed circRNAs
Dolecek TA, Propp JM, Stroup NE, Kruchko C. CBTRUS statistical report: primary brain and central nervous system tumors diagnosed in the United States in 2005–2009. Neuro Oncol. 2012;14(Suppl 5):v1-49.
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Ahmed R, Oborski MJ, Hwang M, Lieberman FS, Mountz JM. Malignant gliomas: current perspectives in diagnosis, treatment, and early response assessment using advanced quantitative imaging methods. Cancer Manag Res. 2014;6:149–70.
Li X, Yang L, Chen LL. The biogenesis, functions, and challenges of circular RNAs. Mol Cell. 2018;71(3):428–42.
Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, Maier L, Mackowiak SD, Gregersen LH, Munschauer M, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495(7441):333–8.
Liu J, Liu T, Wang X, He A. Circles reshaping the RNA world: from waste to treasure. Mol Cancer. 2017;16(1):58.
Jeck WR, Sharpless NE. Detecting and characterizing circular RNAs. Nat Biotechnol. 2014;32(5):453–61.
Salzman J. Circular RNA expression: its potential regulation and function. Trends Genet. 2016;32(5):309–16.
Zhou WY, Cai ZR, Liu J, Wang DS, Ju HQ, Xu RH. Circular RNA: metabolism, functions and interactions with proteins. Mol Cancer. 2020;19(1):172.
Rybak-Wolf A, Stottmeister C, Glazar P, Jens M, Pino N, Giusti S, Hanan M, Behm M, Bartok O, Ashwal-Fluss R, et al. Circular RNAs in the mammalian brain are highly abundant, conserved, and dynamically expressed. Mol Cell. 2015;58(5):870–85.
Song X, Zhang N, Han P, Moon BS, Lai RK, Wang K, Lu W. Circular RNA profile in gliomas revealed by identification tool UROBORUS. Nucleic Acids Res. 2016;44(9):e87.
Achkar NP, Cambiagno DA, Manavella PA. miRNA biogenesis: a dynamic pathway. Trends Plant Sci. 2016;21(12):1034–44.
Hayes J, Peruzzi PP, Lawler S. MicroRNAs in cancer: biomarkers, functions and therapy. Trends Mol Med. 2014;20(8):460–9.
Kristensen LS, Hansen TB, Veno MT, Kjems J. Circular RNAs in cancer: opportunities and challenges in the field. Oncogene. 2018;37(5):555–65.
Jin P, Huang Y, Zhu P, Zou Y, Shao T, Wang O. CircRNA circHIPK3 serves as a prognostic marker to promote glioma progression by regulating miR-654/IGF2BP3 signaling. Biochem Biophys Res Commun. 2018;503(3):1570–4.
Hu D, Zhang Y. Circular RNA HIPK3 promotes glioma progression by binding to miR-124-3p. Gene. 2019;690:81–9.
Li J, Yen C, Liaw D, Podsypanina K, Bose S, Wang SI, Puc J, Miliaresis C, Rodgers L, McCombie R, et al. PTEN, a putative protein tyrosine phosphatase gene mutated in human brain, breast, and prostate cancer. Science. 1997;275(5308):1943–7.
Helseth R, Helseth E, Johannesen TB, Langberg CW, Lote K, Ronning P, Scheie D, Vik A, Meling TR. Overall survival, prognostic factors, and repeated surgery in a consecutive series of 516 patients with glioblastoma multiforme. Acta Neurol Scand. 2010;122(3):159–67.
Endersby R, Baker SJ. PTEN signaling in brain: neuropathology and tumorigenesis. Oncogene. 2008;27(41):5416–30.
Ma J, Benitez JA, Li J, Miki S, PontedeAlbuquerque C, Galatro T, Orellana L, Zanca C, Reed R, Boyer A, et al. Inhibition of nuclear PTEN tyrosine phosphorylation enhances glioma radiation sensitivity through attenuated DNA repair. Cancer Cell. 2019;35(3):504-518.e507.
Guo X, Qiu W, Liu Q, Qian M, Wang S, Zhang Z, Gao X, Chen Z, Xue H, Li G. Immunosuppressive effects of hypoxia-induced glioma exosomes through myeloid-derived suppressor cells via the miR-10a/Rora and miR-21/Pten pathways. Oncogene. 2018;37(31):4239–59.
Wang ZF, Liao F, Wu H, Dai J. Glioma stem cells-derived exosomal miR-26a promotes angiogenesis of microvessel endothelial cells in glioma. J Exp Clin Cancer Res. 2019;38(1):201.
Bachmayr-Heyda A, Reiner AT, Auer K, Sukhbaatar N, Aust S, Bachleitner-Hofmann T, Mesteri I, Grunt TW, Zeillinger R, Pils D. Correlation of circular RNA abundance with proliferation–exemplified with colorectal and ovarian cancer, idiopathic lung fibrosis, and normal human tissues. Sci Rep. 2015;5:8057.
Vo JN, Cieslik M, Zhang Y, Shukla S, Xiao L, Zhang Y, Wu YM, Dhanasekaran SM, Engelke CG, Cao X, et al. The landscape of circular RNA in cancer. Cell. 2019;176(4):869-881.e813.
Zhong Y, Du Y, Yang X, Mo Y, Fan C, Xiong F, Ren D, Ye X, Li C, Wang Y, et al. Circular RNAs function as ceRNAs to regulate and control human cancer progression. Mol Cancer. 2018;17(1):79.
Khan IN, Ullah N, Hussein D, Saini KS. Current and emerging biomarkers in tumors of the central nervous system: possible diagnostic, prognostic and therapeutic applications. Semin Cancer Biol. 2018;52(Pt 1):85–102.
Guarnerio J, Bezzi M, Jeong JC, Paffenholz SV, Berry K, Naldini MM, Lo-Coco F, Tay Y, Beck AH, Pandolfi PP. Oncogenic role of fusion-circRNAs derived from cancer-associated chromosomal translocations. Cell. 2016;166(4):1055–6.
Wang K, Long B, Liu F, Wang JX, Liu CY, Zhao B, Zhou LY, Sun T, Wang M, Yu T, et al. A circular RNA protects the heart from pathological hypertrophy and heart failure by targeting miR-223. Eur Heart J. 2016;37(33):2602–11.
Liu CX, Li X, Nan F, Jiang S, Gao X, Guo SK, Xue W, Cui Y, Dong K, Ding H, et al. Structure and degradation of circular RNAs regulate PKR activation in innate immunity. Cell. 2019;177(4):865-880.e821.
Meng S, Zhou H, Feng Z, Xu Z, Tang Y, Li P, Wu M. CircRNA: functions and properties of a novel potential biomarker for cancer. Mol Cancer. 2017;16(1):94.
Sun J, Li B, Shu C, Ma Q, Wang J. Functions and clinical significance of circular RNAs in glioma. Mol Cancer. 2020;19(1):34.
Zhang H, Liu A, Feng X, Tian L, Bo W, Wang H, Hu Y. MiR-132 promotes the proliferation, invasion and migration of human pancreatic carcinoma by inhibition of the tumor suppressor gene PTEN. Prog Biophys Mol Biol. 2019;148:65–72.
Wu Q, Yu S, Chen J, Li Y, Gao Y. Downregulation of STRAP promotes tumor growth and metastasis in hepatocellular carcinoma via reducing PTEN level. IUBMB Life. 2018;70(2):120–8.
Csolle MP, Ooms LM, Papa A, Mitchell CA. PTEN and other PtdIns(3,4,5)P3 lipid phosphatases in breast cancer. Int J Mol Sci. 2020;21(23):9189.
Choe G, Horvath S, Cloughesy TF, Crosby K, Seligson D, Palotie A, Inge L, Smith BL, Sawyers CL, Mischel PS. Analysis of the phosphatidylinositol 3’-kinase signaling pathway in glioblastoma patients in vivo. Cancer Res. 2003;63(11):2742–6.
Lu H, Yao B, Wen X, Jia B. FBXW7 circular RNA regulates proliferation, migration and invasion of colorectal carcinoma through NEK2, mTOR, and PTEN signaling pathways in vitro and in vivo. BMC Cancer. 2019;19(1):918.
This study was supported by the National Natural Science Foundation of China (No. 81372183) and the National Natural Science Foundation of China Youth Fund Project (No. 81801908, No. 81702490).
Ethics approval and consent to participate
This study was approved by the Medical Ethics Committee of Xiangya Hospital of Central South University, and the ethics approval number is 201703741. All patients signed the informed consent before surgery.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Liu, Y., Chen, S., Peng, G. et al. CircRNA NALCN acts as an miR-493-3p sponge to regulate PTEN expression and inhibit glioma progression. Cancer Cell Int 21, 307 (2021). https://doi.org/10.1186/s12935-021-02001-y