Skip to main content

m6A-mediated upregulation of LINC01003 regulates cell migration by targeting the CAV1/FAK signaling pathway in glioma

Abstract

Background

Long non-coding RNAs (lncRNAs) play important roles in the progression of glioma. Here, we examined the potential functions of a lncRNA, LINC01003, in glioma and characterized the underlying molecular mechanisms.

Methods

The GEIPA2 and Chinese Glioma Genome Atlas (CCGA) databases were employed to analyze gene expression and the overall survival curve in patients with glioma. The functions of LINC01003 in glioma growth and migration were assessed by loss-of-function experiments in vitro and in vivo. RNA sequencing was used to determine the signaling pathways effected by LINC01003. Bioinformatics analysis and RNA immunoprecipitation (RIP) assays were used to explore the mechanism underlying the N6-methyladenine (m6A) modification-dependent upregulation of LINC01003 in glioma.

Results

LINC01003 expression was upregulated in glioma cell lines and tissues. Higher LINC01003 expression predicted shorter overall survival time in glioma patients. Functionally, LINC01003 knockdown inhibited the cell cycle and cell proliferation and migration in glioma cells. Mechanistically, RNA sequencing revealed that LINC01003 mediated the focal adhesion signaling pathway. Furthermore, LINC01003 upregulation is induced by m6A modification regulated by METTL3.

Conclusion

This study characterized LINC01003 as a lncRNA that contributes to tumorigenesis in glioma and demonstrated that the LINC01003-CAV1-FAK axis serves as a potential therapeutic target for glioma.

Introduction

Glioma is one of the most common primary tumors in human central nervous system (CNS), with high rates of occurrence, invasiveness, and recurrence. The overall survival time is short, and the 5-year mortality rate is high [1]. According to the 2021 WHO Classification of Tumors of the Central Nervous System, gliomas are mainly divided into astrocytoma, oligodendrogliomas, glioblastoma, angiocentric glioma, and astroblastoma [2]. The 2022 Cancer Statistics of China indicate that glioma is in the top 10 cancer types by incidence and mortality [3]. Although the treatment strategies for glioma patients include surgery, radiotherapy, chemotherapy, and immunotherapy, the prognosis of glioma patients does not improve significantly upon treatment [4]. Therefore, it is necessary to understand the mechanism underlying the progression of glioma in order to determine clinical biomarkers for timely intervention.

Long non-coding RNAs (lncRNAs) are a group of non-coding RNAs that contain more than 200 nucleotides. Accumulating evidence in recent years has suggested that lncRNAs are essential for diverse cellular and molecular mechanisms [5]. LncRNAs serve as competing endogenous RNAs (ceRNAs) that bind to mircroRNAs (miRNAs) and regulate the expression of mRNAs. LncRNAs also function as scaffolds or guides to modulate protein–protein or protein-DNA interactions [6]. Therefore, lncRNAs are multifunctional regulators associated with DNA replication, transcription, and gene expression. There are several studies showing that lncRNAs are abnormally expressed in glioma and affect multiple cellular processes, including the cell proliferation, migration, apoptosis, and the cell cycle [7, 8]. To date, a large number of glioma-related lncRNAs have been identified. Nevertheless, the molecular mechanism underlying the involvement of lncRNAs in the progression and metastasis of glioma has yet to be clarified.

N6-methyladenine (m6A) RNA modification is one of the most abundant epigenetic processes in eukaryotic mRNA [9, 10]. m6A is also found in lncRNAs and serves an important role in the progression of glioma [11, 12]. m6A is often enriched in stop codons, 3′-untranslated regions (UTRs), and splicing sites and plays a regulatory role in almost every step of the lncRNA life cycle, including transcription, splicing, nuclear export, localization, stability, and degradation [13].

In this study, we demonstrated that lncRNA LINC01003 is highly expressed in glioma and is associated with the poor clinical prognosis of glioma patients. We confirmed that LINC01003 promotes glioma proliferation and migration through in vitro and in vivo functional assays. Furthermore, we showed that LINC01003 can modulate an CAV1/focal adhesion kinase (FAK) signaling pathway in glioma. The RNA m6A methyltransferase METTL3 was shown to epigenetically mediate the upregulation of LINC01003 in glioma.

Methods

Clinical specimens and databases

Glioma tissues (n = 66) and normal brain tissues (n = 14) were obtained from patients who had undergone surgery in the Department of Neurosurgery of the First Affiliated Hospital of Wannan Medical College. The fresh tissue samples were immediately preserved in liquid nitrogen and stored at − 80 °C until RNA was extracted. This project was reviewed and approved by the Ethics Committee of the First Affiliated Hospital of Wannan Medical College. All participants gave written informed consent to allow their tissue to be used for scientific research. The transcriptional data and overall survival in glioma were downloaded from GEPIA2 (Gene Expression Profiling Interactive Analysis 2, http://gepia2.cancer-pku.cn/) [14] and the CGGA database (the Chinese Glioma Genome Atlas, http://www.cgga.org.cn/) [15].

Cell lines and culture

The human glioma cell lines (U87MG, U251, U373, and T98G) and normal human astrocyte cell line (HEB) were obtained from American Type Culture Collection (ATCC, Manassas, USA) and cultured in DMEM (Hyclone, GE Healthcare Life Sciences, USA) medium with 10% fetal bovine serum (FBS; Gibco, Life Technologies, USA). All cell lines were grown in a humidified incubator at 37 °C containing 5% CO2.

Transfection with siRNA

The lncRNA smart silencer negative control (NC) siRNA, LINC01003 smart silencer siRNA, negative control (NC) siRNA, and METTL3 siRNA were purchased from RiboBio (Guangzhou, China). The sequence of the siRNA targeting METTL3 was CAAGTATGTTCACTATGAA. A 6-well plate was seeded with glioma cells at a density of 15 × 104 cells per well. Twenty-four hours later, the smart silencer and siRNA (100 nM) were transfected in cells using the riboFECT™ CP Transfection Kit (RiboBio) following the manufacturer's instructions.

Infection with lentivirus

Recombinant lentiviruses for knocking down LINC01003 (sh-LINC01003) and the corresponding control viruses (sh-NC) were purchased from HANBIO (Shanghai, China). Glioma cell U87MG (15 × 104 cells/well) were seeded into a 6-well plate and transfected with the lentivirus according to the manufacturer’s instructions. Infected cells were treated with 5 μg/ml puromycin for more than 7 days.

Real-time qPCR (RT-qPCR) analysis

Total RNA from tissue or cell lines was isolated using TRIzol reagent (Ambion, Life Technologies, USA). The RNA was then reverse transcribed to complementary DNAs (cDNA) by RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific, Lithuania, USA). RT-qPCR was performed using a QuantiNova™ SYBR® Green PCR Kit (Qiagen, Hilden, Germany) with a Bio-Rad CXF96 PCR system (Hercules, CA, USA) according to the manufacturer's instructions. GAPDH was used as reference gene for mRNA and lncRNA detection. Relative mRNA and lncRNA expression was determined by the 2−ΔΔCt method. The qRT-PCR primers were purchased from Sangon Biotech (Shanghai, China). LINC01003: F-5′-TCTCTGCGTCGCTCCTTCTCGC-3′, R-5′-CCTTCCGCGCCAAGTAGTGC-3′; CAV1: F-5′-GACTCGGAGGGACATCT CTACACC-3, R-5′-GTCGTACACTTGCTTCTCGCTCAG-3′.

Cell proliferation and cycle assay

EdU staining assay (Cell-LightTM Apollo 488 Stain Kit, RiboBio) was used to detect cell proliferation following the manufacturer’s protocol. Briefly, transfected cells were seeded in 96-well plates and incubated with DMEM medium containing EdU for 2 h. Following cell fixation (4% formaldehyde) and cell membrane permeation (0.5% Triton X-100), the cells were stained with Apollo dye solution and Hoechst33342. The images were photographed with a Nikon Ti-U fluorescent microscope (Tokyo, Japan). The EdU-positive cells were counted using Photoshop CS6 (Adobe Systems Incorporated, San Jose, CA, USA).

For the cell cycle assay, transfected cells were collected and fixed with 75% pre-chilled ethanol for 24 h at 4 °C. Then the cells were stained with RNase A/PI (propidium iodide) mixture and incubated in the dark at 37 °C for 30 min. The cells were detected by a CytoFLEX flow cytometer (Beckman Coulter, CA, USA). The cell cycle was analyzed with Modfit LT32 software (Verity Software House, Topsham, ME).

Cell migration assay

Cell migration was determined by a transwell assay using a transwell chamber with an 8.0-mm pore size (Corning, NY, USA). Then 2 × 104 transfected cells, which were suspended in 100 μl of serum-free BMDM medium, were seeded in the upper chamber and 600 μl of BMDM medium containing 10% FBS was added to the lower chamber. Then the cells were incubated for 24 h, and the bottom side of the upper chamber was fixed with 4% formaldehyde for 10 min. After staining with 0.1% crystal violet, the cells were captured by an inverted microscope (Nikon Ti-U). The migrated cells were counted using Photoshop CS6 (Adobe Systems Incorporated).

The migration ability of the cells was also determined by a wound-healing assay. The transfected cells were cultured to 100% confluence in a 6-well plate and scratched with a 200-μl pipette tip. After washing with PBS, the cells were incubated in DMEM medium without FBS. The wound closures were captured using an inverted microscope (Nikon Ti-U) at h 0, 24, and 48.

Western blots

The protein expression of p-FAK, FAK, E-Cadherin, N-Cadherin (Rabbit mAb; CST), matrix metallopeptidase 2 (MMP2), MMP9, and CAV1 (Rabbit Polyclonal; Proteintech) was determined as described previously [16, 17]. The total protein of cells was extracted using Laemmli 2 × concentrate sample buffer (Sigma, MO, USA). The proteins were separated by 10% polyacrylamide gel electrophoresis and transferred to 0.22-nm nitrocellulose membranes (Pall Corporation, East Hills, USA). Then the membranes were blocked with 5% BSA (bovine serum albumin) and incubated with primary and secondary antibodies. After incubation with an enhanced chemiluminescence kit (ECL, Millipore, Burlington, USA), the protein bands were captured by a ChemiDoc Touch Imaging System (Bio-Rad, Hercules, USA).

RNA sequencing (RNA-seq) analysis

Glioma cells were transfected with lncRNA smart silencer NC (n = 2) or LINC01003 smart silencer (n = 2) siRNA for 48 h. The RNA-seq analysis was performed by Novogene Co., Ltd (Beijing, China). RNA integrity was determined by the RNA Nano 6000 Assay Kit using the Bioanalyzer 2100 system (Agilent Technologies, CA, USA). The clustering of the index-coded samples was conducted by TruSeq PE Cluster Kit v3-cBot-HS (Illumia) on a cBot Cluster Generation System. Following cluster generation, the library was sequenced on the Illumina Novaseq platform (Illumina, CA, USA). Genes with alog2FoldChange ≥ 0 and adjusted P value ≤ 0.05 between NC and LINC01003 smart silencer cells were assigned as significantly differentially expressed genes (DEGs). The Cluster Profiler R package was used to analyze the statistical enrichment of DEGs in KEGG pathways.

In vivo xenograft model

BALB/c nude mice (5 weeks old) were acquired from the Experimental Animal Center of Qinglongshan (Nanjing, China) and housed in a specific pathogen-free (SPF) mouse colony. For in vivo tumorigenesis experiments, U87MG cells (1 × 107) were suspended in 100 μl PBS and implanted subcutaneously into the mice. The width (W) and length (L) of tumors were monitored every 5 days with a caliper. The tumor volume was calculated using the following formula: W2 × L × 0.5. After four weeks, the mice were sacrificed, and the tumors were isolated and weighed. The experiments using the nude mice were approved by the ethics committee of the First Affiliated Hospital of Wannan Medical College.

m 6 A modification prediction

The online tool RMBase v2.0 (http://rna.sysu.edu.cn/rmbase/) [18], which integrates epitranscriptome sequencing data, was used to explore the m6A modification of LINC01003. The m6A2Target (http://m6a2target.canceromics.org) database was used to predict the relationship between m6A enzymes with LINC01003 expression.

RNA immunoprecipitation

The m6A-modified LINC01003 levels were determined using the Magna RIP RNA-binding Protein Immunoprecipitation Kit (Millipore) combined with RT-qPCR. The RIP assay for RNA binding with anti-m6A (Mouse monoclonal, Abcam, MA, USA) or anti-IgG (CST) was carried out in accordance with the manufacturer’s instructions. Then RT-qPCR was performed to detect the expression level of LINC01003 m6A modification.

Statistical analyses

All experimental data are presented as the mean ± SD from three independent replicates. All statistical analyses were conducted using GraphPad Prism 5 (La Jolla, CA, USA). Student’s t-test was used to compare the data differences between two groups. One-way ANOVA followed by post hoc test was used to compare the data differences between three groups. *P < 0.05, **P < 0.01, and ***P < 0.001.

Results

LncRNA LINC01003 is highly expressed in glioma patients and cells

Previous research has shown that LncRNA LINC01003 expression in glioblastoma stem-like cells is significantly higher than in conventional glioma cells [5]. From this report we speculated that LINC01003 would play a critical role in glioma tumorigenesis. The GEPIA2 database showed that LINC01003 expression in GBM (glioblastoma multiforme) or LGG (lower grade glioma) is significantly higher than in normal brain tissues (Fig. 1A). Importantly, the Kaplan–Meier survival analysis from the GEPIA2 database indicated that GBM and LGG patients with a higher LINC01003 expression level were associated with shorter overall survival than those with lower LINC01003 expression level (Fig. 1B). We then performed RT-qPCR to analyze the expression of LINC01003 in brain tissues from 66 glioma patients and 14 normal patients of Yijishan Hospital. We observed that LINC01003 expression was significantly upregulated in glioma patients compared with normal controls (Fig. 1C). In addition, LINC01003 expression in glioma cell lines was significantly higher than in normal HEB cells (Fig. 1D). As LINC01003 expression was highest in U87MG and U251 cells, LINC01003 smart silence was subsequently transfected into U87MG and U251 cells. RT-qPCR assays showed that LINC01003 smart silence system significantly knockdown LINC01003 expression (Fig. 1E). Taken together, these results suggested that LINC01003 expression was upregulated in glioma and positively associated with poor prognosis.

Fig. 1
figure 1

LINC01003 is upregulated in glioma tissues and cells. A The expression of LINC01003 in glioblastoma multiforme (GBM), lower grade glioma (LGG), and normal brain tissues from the GEPIA2 database. B Kaplan–Meier curves from the GEPIA2 database determined the relevance between LINC01003 expression and overall survival in GBM and/or LGG patients. C RT-qPCR was used to detect the relative expression of LINC01003 in 14 normal brain tissues and 66 glioma tissues. D RT-qPCR was used to detect the relative expression of LINC01003 in HEB, U87MG, U251, U373, and T98G cell lines. E Real-time qPCR (RT-qPCR) was used to detect the relative expression of LINC01003 in U87MG and U251 cells lines, where the cells were transfected with lncRNA smart silencer negative control (NC) and LINC01003 smart silencer siRNA; *P < 0.05, **P < 0.01, and ***P < 0.001

Effects of LINC01003 on glioma cell cycle and cell proliferation

To further verify the function of LINC01003 in glioma, we determined the effect of LINC01003 on cell proliferation and the cell cycle by silencing LINC01003 in U87MG and U251 glioma cells using lncRNA smart silence. EdU assays indicated that LINC01003 knockdown significantly inhibited DNA replication in glioma cells (Fig. 2A–B). Uncontrolled cell cycle progression is considered to be a key marker to promote malignant proliferation of tumor cells. Thus, inhibition of cell cycle is an effective strategy in the clinical treatment of tumors. We then employed flow cytometry assays to determine the effect of LINC01003 on the cell cycle. We found that LINC01003-silenced glioma cells showed changes in the distribution of the cell cycle, with increased arrest at S phase and decreased proportion of G0/G1 phase (Fig. 2C–D). These results suggested that LINC01003 positively mediates glioma cell proliferation in vitro.

Fig. 2
figure 2

Effects of LINC01003 on the proliferation and cell cycle of glioma cells. AB EdU (5-Ethynyl-2′-deoxyuridine) staining assay was performed to determine the proliferation of U87MG and U251 cells with LINC01003 knockdown. The percentage of EdU-positive cells were counted; scale bar, 200 μm. CD Flow cytometry with propidium iodide (PI) staining was used to detect the cell cycle. A histogram showing the proportion of glioma cells in G0/G1, S, and G2/M phase. Data represent the mean ± SD from three independent experiments; *P < 0.05 and **P < 0.01

Effects of LINC01003 on glioma cell migration

Next, we examined the effects of LINC01003 on the migration of glioma cells by transwell and wound-healing assays. In the transwell assays, LINC01003 knockdown decreased the number of migrated cells (Fig. 3A). Wound-healing assays revealed that LINC01003 depletion decreased the speed of cell migration (Fig. 3B). In order to verify that LINC01003 knockdown reduced cell migration, we detected the protein expression of key regulators required for cell migration. MMP2 and MMP9 are two matrix metalloproteinases that play a critical role in tumor cell invasion and metastasis [19, 20]. Epithelial–mesenchymal transition (EMT), which is involved in the early stage of metastasis, is characteristic by N-cadherin (N-Cad) upregulation and E-cadherin (E-Cad) downregulation [21, 22]. LINC01003 knockdown significantly inhibited the expression of MMP2 and MMP9, whereas LINC01003 silencing increased the expression of E-Cad and decreased the expression of N-Cad (Fig. 3C). Collectively, this result indicated that LINC01003 knockdown inhibits the migration of glioma cells.

Fig. 3
figure 3

Effects of LINC01003 on the cell migration of glioma cells. A A transwell assay was employed to determine the migration of U87MG and U251 cells with LINC01003 knockdown. The glioma cells that migrated were counted; scale bars, 200 μm. B Wound healing assays was also employed to determine the migration of U87MG and U251 cells with LINC01003 knockdown. The scratch was measured 0 and 48 h after the scratch was performed; scale bars, 200 μm. C Western blot was used to determine the expression of matrix metallopeptidase 2 (MMP2), MMP9, E-Cadherin (E-cad) and N-Cadherin (N-cad) in U87MG and U251 cells with LINC01003 knockdown. *P < 0.05, **P < 0.01, and ***P < 0.001

LINC01003 regulates the focal adhesion kinase (FAK) signaling pathway

To investigate the mechanism of LINC01003 action in glioma, RNA transcriptome sequencing was performed in U251 glioma cells with or without LINC01003 knockdown. LINC01003 suppression upregulated 368 genes and downregulated 557 genes (log2FoldChange ≥ 0, adjusted P value ≤ 0.05) (Fig. 4A–B and Additional file 1: Table S1). KEGG pathway enrichment analysis revealed that the 557 downregulated genes were enriched in influenza A, c-type lectin receptor signaling pathway, TNF signaling pathway, measles, focal adhesion, herpes simplex infection, and FoxO signaling pathway (Fig. 4C and Additional file 2: Table S2). Focal adhesion ranked fifth in the KEGG enrichment category and included 12 genes (PIK3R3, PGF, VCL, LAMB3, THBS2, CAV1, KDR, PRKCA, LAMC2, MYL5, TNC, LAMA2, PIK3R1, BIRC3, ITGB6, and MYLPF). Our sequencing data suggested that LINC01003 mediates the expression of focal adhesion-related genes, indicating that LINC01003 might regulate focal adhesion in glioma cells. Focal adhesion kinase (FAK) is a key component of focal adhesions and crucial for cancer cell migration [23,24,25]. We found that LINC01003 depletion significantly reduced FAK phosphorylation in glioma cells (Fig. 4D).

Fig. 4
figure 4

Knockdown of LINC01003 inhibited the CAV1-FAK signaling pathway. A The heatmap reflected the differently expressed genes (DEGs) regulated by LINC01003 knockdown. B A volcano plot of DEGs from (A). C Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was performed for the 557 downregulated genes regulated by LINC01003 knockdown from (B). D Western blot was used to detect the expression levels of p-FAK and CAV1 in U87MG and U251 cells with LINC01003 knockdown. Up: upregulation; Down: downregulation; No: no change

CAV1 (Caveolin-1), a membrane protein of caveolae, was reported to be involved in invasion via the FAK-mediated adhesion signaling pathway [26]. The GEPIA2 database showed that CAV1 expression in GBM or LGG was significantly higher than in normal brain tissues (Additional file 3: Fig. S1A). The CGGA database showed that CAV1 expression was increased in glioma tissues with advanced tumor stage (Additional file 3: Fig. S1B). The GEPIA2 and CGGA databases demonstrated that glioma patients with higher CAV1 expression level were associated with shorter overall survival than those with a lower CAV1 expression level (Additional file 3: Fig. S1C and D). Next, western blot and RT-qPCR were performed to verify whether LINC01003 could regulate CAV1 expression. Consistent with the transcriptome sequencing data, CAV1 expression was decreased in the LINC01003 knockdown cells (Fig. 4D and Additional file 3: Fig. S1E). Collectively, these data demonstrated that LINC01003 may regulate FAK phosphorylation through CAV1.

LINC01003 knockdown inhibited glioma proliferation in vivo

We further employed a xenograft tumor model to assess the tumor-suppressive effects of LINC01003 in vivo. The tumor growth in nude mice injected with sh-LINC01003 cells was slower as compared with those injected with sh-NC cells (Fig. 5A). LINC01003 knockdown inhibited glioma tumor volume and weight in nude mice (Fig. 5B–C). RT-qPCR assays showed that LINC01003 expression was decreased in the sh-LINC01003 xenograft tumors (Fig. 5D). The expression of CAV1 and p-FAK was then examined by western blot in xenograft tumors. As shown in Fig. 5E, the expression of CAV1 and p-FAK was significantly downregulated in the LINC01003-silenced group. In conclusion, these results suggested that knockdown of LINC01003 inhibits glioma growth in vivo.

Fig. 5
figure 5

LINC01003 knockdown represses glioma growth in vivo. A A subcutaneous xenograft tumor model was generated by transfecting U87MG cells with shNC or shLINC01003. Tumors were removed from the nude mice and photographed at the end of experiment. B Once the mice were injected with the cells, he tumor volume in nude mice was measured every 5 days. The volume curves were plotted from day 5 to day 30. C Quantification of tumor weight in nude mice from (A). D RT-qPCR was used to determine the expression of LINC01003 in the tumors from nude mice. E Western blot was used to detect the expression levels of p-FAK and CAV1 in the tumors from nude mice. Data represent the mean ± SD; *P < 0.05, **P < 0.01, and ***P < 0.001

To determine the role of the LINC01003/CAV1 axis in glioma, we employed functional rescue experiments in vitro and in vivo. Genetic addition of CAV1 alleviated the inhibitory effect of LINC01003 deficiency on the migration and p-FAK levels of glioma cells (Additional file 4: Fig. S2A and B). The xenografts result in nude mice further indicated that overexpression of CAV1 partially relieved the inhibition of tumor formation induced by knockdown of LINC01003 in vivo (Additional file 4: Fig. S2C). These findings suggested that the LINC01003/CAV1 axis regulates glioma tumorigenesis in vitro and in vivo.

m 6 A modification is associated with LINC01003 expression in glioma cells

Recent advances have suggested that N6-methyladenine (m6A) modifications play a critical role in the regulation of the posttranscriptional level of lncRNAs. Thus, we speculated that m6A modifications were responsible for the upregulation of LINC01003 in glioma. The m6A modification analysis online software, RMBase v2.0 database, predicted 15 potential m6A modification sites in the exon of LINC01003 transcript (Additional file 5: Table S3). Then the m6A modification level of LINC01003 in glioma cells was determined by RNA immunoprecipitation and RT-qPCR assays. Compared with the IgG control group, the relative enrichment level of m6A on LINC01003 in glioma cells was significantly increased (Fig. 6A).

Fig. 6
figure 6

N6-methyladenosine (m6A) modification regulates LINC01003 expression in glioma cells. A A RNA immunoprecipitation (RIP) assay combined with RT-qPCR was used to determine the relative m6A enrichment at LINC01003 in U87MG and U251 cells. B RT-qPCR was used to determine the expression of LINC01003 in glioma cells with NC or METTL3 siRNA. C A RIP assay combined with RT-qPCR was used to determine the relative m6A enrichment at LINC01003 in glioma cells with NC or METTL3 siRNA; ***P < 0.001

The online tool, m6A2Target, was used to explore the relationship between m6A enzymes with LINC01003. We found 4 “writer” enzymes (RBM15B, WTAP, METTL3, HAKAI) and 1 “eraser” enzyme (ALKBH5) that may perturb the expression of LINC01003 (Additional file 6: Fig. S3A). The GEPIA2 database showed that METTL3 had the highest positive correlation with LINC01003 in glioma tissues (Additional file 6: Fig. S3B–F). METTL3 (termed “writers”), the key component of the m6A methyltransferase complex, has been reported to regulate the malignant progression of glioma [27]. Previous research has shown that METTL3 siRNA treatment significantly knocked down the expression of METTL3 in U87MG and U251 cells [16]. We then examined the regulatory effect of METTL3 on LINC01003. In glioma cells, knockdown of METTL3 significantly reduced the expression of LINC01003 (Fig. 6B). An anti-m6A RIP assay combined with RT-qPCR indicated that the m6A modification of LINC01003 was remarkably decreased in LINC01003-silenced glioma cells (Fig. 6C). Collectively, these data suggested that METTL3-mediated m6A modification is associated with the upregulation of LINC01003 in glioma.

Discussion

Glioma is the most common malignant tumor in the central nervous system [28]. Recent advances have found that numerous functional lncRNAs are abnormally expressed in glioma tissues. LncRNAs have been considered as central regulators in glioma progression [29]. A recent study demonstrated that LINC01003 acts as a suppressor in multiple myeloma [30]. In glioma, LINC01003 is highly expressed in glioma stem cells [5]. In this study, we identified a new lncRNA, LINC01003, which showed a significant trend of increased expression in glioma tissues. Importantly, we found that higher LINC01003 expression level is associated with poor prognosis in glioma patients. Furthermore, knockdown of LINC01003 significantly suppressed the proliferation of glioma cells in vitro and tumorigenesis in vivo. Based on these results, LINC01003 may represent a predictor of prognosis and potential therapeutic target for glioma therapy.

EMT is an important process that bestows cells an increased invasive or migrated phenotype. The main reason for the recurrence of glioma is largely attributed to the increased invasion and migration ability induced by EMT [31, 32]. Cadherins, including N-cad and E-cad, matrix metalloproteinase 2 (MMP2), and MMP9 were demonstrated to be EMT-associated proteins [33,34,35]. A loss-of-function assay indicated that LINC01003 can trigger the migration of glioma cells in vitro. This was confirmed by the observed upregulation of E-Cad and downregulation of N-Cad, MMP2, and MMP9. Our findings suggested that LINC01003 may be a crucial effector molecule that induces EMT in glioma.

To undercover the mechanism underlying LINC01003 in the tumorigenesis of glioma, RNA-seq analysis were performed. KEGG enrichment indicated that LINC01003 may be involved in a focal adhesion signaling pathway that is closely associated with cell proliferation and migration. Focal adhesion kinase (FAK), a central mediator of cell adhesion, was shown to play important roles in glioma cell proliferation and migration [36,37,38]. However, glioma-associated lncRNAs targeting focal adhesion signaling pathways have not been fully determined yet. In our study, we found LINC01003 knockdown repressed the phosphorylation of FAK. Interestingly, we identified a new signaling axis, LINC01003/CAV1/FAK, through which LINC01003 regulates focal adhesion in glioma. Notably, expression of CAV1 was significantly upregulated in high-grade glioma patients, and high expression of CAV1 is associated with a poor prognosis in glioma. Previous studies have demonstrated that CAV1 promotes glioma cell proliferation and vasculogenic mimicry and regulates focal adhesion by effecting FAK phosphorylation [26, 39]. Therefore, our study indicated that LINC01003 is a positive regulator in the CAV1/FAK signaling pathway in glioma cells.

To determine the upstream mechanism of LINC01003 overexpression in glioma, we investigated whether m6A methylation regulates LINC01003 expression. In our study, we combined a RIP assay with bioinformatics analysis and found that LINC01003 is significantly enriched with m6A methylation in glioma cells. The m6A methylation is the most common RNA modification in eukaryotic cells. Studies have demonstrated that m6A methylation exists in many lncRNAs and modulates almost all stages of the lncRNA life cycle, including transcription, nuclear export, splicing, and degradation. Therefore, the m6A methylation is related to cancer progression, including cell proliferation and migration [40,41,42]. However, the functions of m6A modification in lncRNAs remain largely unknown. METTL3 is a pivotal m6A methyltransferase and has recently been reported to be involved in lncRNA stability in tumor progression [43]. We found that METTL3 mediates the expression and m6A modification of LINC01003, suggesting that METTL3-mediated m6A modification may be responsible for the upregulation of LINC01003 in glioma.

Conclusions

Our results demonstrated the pro-metastatic function of LINC01003 in the diffusion of glioma cells. Additionally, we uncovered a new signaling axis, LINC01003/CAV11/FAK, involved in the modulation of cell migration. Moreover, METTL3-mediated m6A modification may increase LINC01003 expression and function in glioma cells. Taken together, our present study provides a potential target that may serve as a promising prognostic biomarker and potential therapeutic target for glioma patients. Thus, clarifying the functions and mechanism of LINC01003 will provide anti-metastatic therapies in glioma.

Availability of data and materials

The processed data that support the findings of this study are included in the manuscript and the supplementary materials.

Abbreviations

CNS:

Central nervous system

WHO:

World Health Organization

lncRNAs:

Long non-coding RNAs

ceRNAs:

Competing endogenous RNAs

m6A:

N6-Methyladenosine

UTRs:

Untranslated regions

FAK:

Focal adhesion kinase

GEPIA2:

Gene expression profiling interactive analysis 2

CGGA:

Chinese glioma genome atlas

HEB:

Human astrocyte cell line

ATCC:

American type culture collection

DMEM:

Dulbecco's modified eagle medium

FBS:

Fetal bovine serum

NC:

Negative control

GBM:

Glioblastoma multiforme

LGG:

Lower grade glioma

EMT:

Epithelial–mesenchymal transition

KD:

Knockdown

KEGG:

Kyoto encyclopedia of genes and genomes

DEGs:

Differentially expressed genes

ECL:

Enhanced Chemiluminescence

RIP:

RNA imunoprecipitation

References

  1. Xu S, Tang L, Li X, Fan F, Liu Z. Immunotherapy for glioma: current management and future application. Cancer Lett. 2020;476:1–12.

    Article  CAS  PubMed  Google Scholar 

  2. Louis DN, Perry A, Wesseling P, et al. The 2021 WHO classification of tumors of the central nervous system: a summary. Neuro Oncol. 2021;23(8):1231–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Xia C, Dong X, Li H, et al. Cancer statistics in China and United States, 2022: profiles, trends, and determinants. Chin Med J (Engl). 2022;135(5):584–90.

    Article  PubMed  Google Scholar 

  4. Lapointe S, Perry A, Butowski NA. Primary brain tumours in adults. Lancet. 2018;392(10145):432–46.

    Article  PubMed  Google Scholar 

  5. Xiong Z, Wang L, Wang Q, Yuan Y. LncRNA MALAT1/miR-129 axis promotes glioma tumorigenesis by targeting SOX2. J Cell Mol Med. 2018;22(8):3929–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Li Y, Liu X, Cui X, et al. LncRNA PRADX-mediated recruitment of PRC2/DDX5 complex suppresses UBXN1 expression and activates NF-κB activity, promoting tumorigenesis. Theranostics. 2021;11(9):4516–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Qi ZY, Wang LL, Qu XL. lncRNA LINC00355 acts as a novel biomarker and promotes glioma biological activities via the regulation of miR-1225/FNDC3B. Dis Markers. 2021;2021:1683129.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Li W, Dong X, He C, et al. LncRNA SNHG1 contributes to sorafenib resistance by activating the Akt pathway and is positively regulated by miR-21 in hepatocellular carcinoma cells [published correction appears in J Exp Clin Cancer Res 2021 Dec 1;40(1):377]. J Exp Clin Cancer Res. 2019;38(1):183.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Chang YZ, Chai RC, Pang B, et al. METTL3 enhances the stability of MALAT1 with the assistance of HuR via m6A modification and activates NF-κB to promote the malignant progression of IDH-wildtype glioma. Cancer Lett. 2021;511:36–46.

    Article  CAS  PubMed  Google Scholar 

  10. Qian B, Wang P, Zhang D, Wu L. m6A modification promotes miR-133a repression during cardiac development and hypertrophy via IGF2BP2. Cell Death Discov. 2021;7(1):157.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Cui Y, Zhang C, Ma S, et al. RNA m6A demethylase FTO-mediated epigenetic up-regulation of LINC00022 promotes tumorigenesis in esophageal squamous cell carcinoma. J Exp Clin Cancer Res. 2021;40(1):294.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Liu H, Xu Y, Yao B, Sui T, Lai L, Li Z. A novel N6-methyladenosine (m6A)-dependent fate decision for the lncRNA THOR. Cell Death Dis. 2020;11(8):613.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Qian X, Yang J, Qiu Q, et al. LCAT3, a novel m6A-regulated long non-coding RNA, plays an oncogenic role in lung cancer via binding with FUBP1 to activate c-MYC. J Hematol Oncol. 2021;14(1):112.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Tang Z, Kang B, Li C, Chen T, Zhang Z. GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 2019;47(W1):W556–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Zhao Z, Zhang KN, Wang Q, et al. Chinese Glioma Genome Atlas (CGGA): a comprehensive resource with functional genomic data from Chinese Glioma patients. Genom Proteom Bioinform. 2021;19(1):1–12.

    Article  CAS  Google Scholar 

  16. Zhu X, Yang H, Zhang M, et al. YTHDC1-mediated VPS25 regulates cell cycle by targeting JAK-STAT signaling in human glioma cells. Cancer Cell Int. 2021;21(1):645.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Zhu X, Chen T, Yang H, Lv K. Lactate induced up-regulation of KLHDC8A (Kelch domain-containing 8A) contributes to the proliferation, migration and apoptosis of human glioma cells. J Cell Mol Med. 2020;24(20):11691–702.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Xuan JJ, Sun WJ, Lin PH, et al. RMBase v2.0: deciphering the map of RNA modifications from epitranscriptome sequencing data. Nucleic Acids Res. 2018;46(D1):D327–34.

    Article  CAS  PubMed  Google Scholar 

  19. Banday MZ, Sameer AS, Mir AH, Mokhdomi TA, Chowdri NA, Haq E. Matrix metalloproteinase (MMP) -2, -7 and -9 promoter polymorphisms in colorectal cancer in ethnic Kashmiri population—a case-control study and a mini review. Gene. 2016;589(1):81–9.

    Article  CAS  PubMed  Google Scholar 

  20. Das S, Amin SA, Jha T. Inhibitors of gelatinases (MMP-2 and MMP-9) for the management of hematological malignancies. Eur J Med Chem. 2021;223: 113623.

    Article  CAS  PubMed  Google Scholar 

  21. Giannelli G, Koudelkova P, Dituri F, Mikulits W. Role of epithelial to mesenchymal transition in hepatocellular carcinoma. J Hepatol. 2016;65(4):798–808.

    Article  CAS  PubMed  Google Scholar 

  22. Loh CY, Chai JY, Tang TF, et al. The E-Cadherin and N-Cadherin switch in epithelial-to-mesenchymal transition: signaling, therapeutic implications, and challenges. Cells. 2019;8(10):1118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Shen J, Cao B, Wang Y, et al. Hippo component YAP promotes focal adhesion and tumour aggressiveness via transcriptionally activating THBS1/FAK signalling in breast cancer. J Exp Clin Cancer Res. 2018;37(1):175.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Jiang K, Yao G, Hu L, et al. MOB2 suppresses GBM cell migration and invasion via regulation of FAK/Akt and cAMP/PKA signaling. Cell Death Dis. 2020;11(4):230.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. de Semir D, Bezrookove V, Nosrati M, et al. PHIP drives glioblastoma motility and invasion by regulating the focal adhesion complex. Proc Natl Acad Sci U S A. 2020;117(16):9064–73.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Trimmer C, Whitaker-Menezes D, Bonuccelli G, et al. CAV1 inhibits metastatic potential in melanomas through suppression of the integrin/Src/FAK signaling pathway. Cancer Res. 2010;70(19):7489–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Visvanathan A, Patil V, Arora A, et al. Essential role of METTL3-mediated m6A modification in glioma stem-like cells maintenance and radioresistance. Oncogene. 2018;37(4):522–33.

    Article  CAS  PubMed  Google Scholar 

  28. Chao M, Liu N, Sun Z, et al. TGF-β signaling promotes glioma progression through stabilizing Sox9. Front Immunol. 2021;11: 592080.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Peng Z, Liu C, Wu M. New insights into long noncoding RNAs and their roles in glioma. Mol Cancer. 2018;17(1):61.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Wu L, Xia L, Chen X, Ruan M, Li L, Xia R. Long non-coding RNA LINC01003 suppresses the development of multiple myeloma by targeting miR-33a-5p/PIM1 axis. Leuk Res. 2021;106: 106565.

    Article  CAS  PubMed  Google Scholar 

  31. Noronha C, Ribeiro AS, Taipa R, et al. Cadherin expression and EMT: a focus on gliomas. Biomedicines. 2021;9(10):1328.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Wang W, Hao Y, Zhang A, et al. miR-19a/b promote EMT and proliferation in glioma cells via SEPT7-AKT-NF-κB pathway. Mol Ther Oncolytics. 2021;20:290–305.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Xu A, Wang X, Luo J, et al. Overexpressed P75CUX1 promotes EMT in glioma infiltration by activating β-catenin. Cell Death Dis. 2021;12(2):157.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Zhang J, Cai H, Sun L, et al. LGR5, a novel functional glioma stem cell marker, promotes EMT by activating the Wnt/β-catenin pathway and predicts poor survival of glioma patients. J Exp Clin Cancer Res. 2018;37(1):225.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Peng L, Wen L, Shi QF, et al. Scutellarin ameliorates pulmonary fibrosis through inhibiting NF-κB/NLRP3-mediated epithelial-mesenchymal transition and inflammation. Cell Death Dis. 2020;11(11):978.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Lin X, Zhuang S, Chen X, et al. lncRNA ITGB8-AS1 functions as a ceRNA to promote colorectal cancer growth and migration through integrin-mediated focal adhesion signaling. Mol Ther. 2022;30(2):688–702.

    Article  CAS  PubMed  Google Scholar 

  37. Tapial Martínez P, López Navajas P, Lietha D. FAK structure and regulation by membrane interactions and force in focal adhesions. Biomolecules. 2020;10(2):179.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Chen X, Guo ZQ, Cao D, Chen Y, Chen J. MYC-mediated upregulation of PNO1 promotes glioma tumorigenesis by activating THBS1/FAK/Akt signaling. Cell Death Dis. 2021;12(3):244.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Chen W, Cheng X, Wang X, Hu W, Wang J, Liao C. Caveolin-1 promotes tumor cell proliferation and vasculogenic mimicry formation in human glioma. Braz J Med Biol Res. 2021;54(10): e10653.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Huang H, Weng H, Chen J. m6A modification in coding and non-coding RNAs: roles and therapeutic implications in cancer. Cancer Cell. 2020;37(3):270–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Lan Y, Liu B, Guo H. The role of M6A modification in the regulation of tumor-related lncRNAs. Mol Ther Nucleic Acids. 2021;24:768–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Yi YC, Chen XY, Zhang J, Zhu JS. Novel insights into the interplay between m6A modification and noncoding RNAs in cancer. Mol Cancer. 2020;19(1):121.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Liu HT, Zou YX, Zhu WJ, et al. lncRNA THAP7-AS1, transcriptionally activated by SP1 and post-transcriptionally stabilized by METTL3-mediated m6A modification, exerts oncogenic properties by improving CUL4B entry into the nucleus. Cell Death Differ. 2022;29(3):627–41.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by grants from the National Natural Science Foundation of China (grant nos. 81870017, 82170368, 82100019, 82072370), Natural Science Foundation of Anhui Province (grant nos. 2008085QH367, 2108085J44, 2108085QH307, 1908085QH380), Key Projects of Natural Science Research of Universities in Anhui Province (grant nos. KJ2020A0607), Funding of “Peak” Training Program for Scientific Research of Yijishan Hospital, Wannan Medical College (grant nos.GF2019G10, GF2019G15), Key Scientific Research Project Training Fund of Wannan Medical College (grant no. WK2021ZF05), Youth Project Training Fund of Wannan Medical College (grant no. WK2021F08).

Author information

Authors and Affiliations

Authors

Contributions

KL and XZ designed the experiments. XZ and XW performed the functional examination of cells, HY and XL performed the western blot examination, QX and MZ performed the animal examination. XZ analyzed all of the data. XZ wrote the manuscript. KL and XZ revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Kun Lv.

Ethics declarations

Ethics approval and consent to participate

This study was reviewed and approved by the Ethics Committee of the Yijishan Hospital of Wannan Medical College. The informed consent was obtained from all participants.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Table S1. The DEGs between NC and LINC01003 smart silencer siRNA.

Additional file 2

: Table S2. Top 20 KEGG enrichment pathways in the downregulated mRNAs.

Additional file 3

: Fig. S1. CAV1 is upregulated in glioma tissues and cells.The expression of CAV1 in GBM, LGG, and normal brain tissues from the Gene Expression Profiling Interactive Analysis 2database.The expression of CAV1 in World Health Organizationgrade II–IV tumors from the Chinese Glioma Genome Atlasdatabase. Kaplan–Meier curves from theGEPIA2 orCGGA database determined the relevance between CAV1 expression and overall survival in each subtype or each grade of glioma patients.RT-qPCR was used to detect the expression level of CAV1 in U87MG and U251 cells with LINC01003 knockdown. * P < 0.05, ** P < 0.01, and ***P < 0.001.

Additional file 4

: Fig. S2. The LINC01003/CAV1 axis regulates glioma in vitro and in vivo.A transwell assay was employed to determine the migration ability of U87MG cells. The migrating glioma cells were counted.Western blot was used to detect the expression levels of p-FAK in U87MG cells.An xenograft tumor model was generated by U87MG cells. The tumor volume and weight in nude mice were measured. Data represent the mean ± SD; * P < 0.05, ** P < 0.01, and ***P < 0.001.

Additional file 5

: Table S3. The m6A sites on LINC01003 were predicted by the GEIPA2 database.

Additional file 6

: Fig. S3. The relationship between m6A enzymes with LINC01003.m6A enzymes may perturb the expression of LINC01003.The correlation between m6A enzymes with LINC01003 in glioma tissues from the GEPIA2 database.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhu, X., Wu, X., Yang, H. et al. m6A-mediated upregulation of LINC01003 regulates cell migration by targeting the CAV1/FAK signaling pathway in glioma. Biol Direct 18, 27 (2023). https://doi.org/10.1186/s13062-023-00386-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13062-023-00386-6

Keywords