3.1 LINC00857 is upregulated in PC and associated with poor prognosis.
To identify lncRNAs involved in the progression and metastasis of PC, we compared the expression of lncRNAs in PC and normal samples from TCGA (The Cancer Genome Atlas, https://cancergenome.nih.gov/) and GTEx (Genotype-Tissue Expression, https://commonfund.nih.gov/GTex) databases. The differential gene expression between them is shown in volcano plots and a heatmap (|log2FoldChange| >1, P < 0.05; Fig. 1A, Fig. S1A). One representative lncRNA, LINC00857, was considered to have abnormal upregulation.
Through Coding Potential Calculator analysis, LINC00857 had a very low protein-coding ability (Supplementary Fig. S1B-E), which means that LINC00857 was indeed a long noncoding RNA. After extracting and comparing the expression data of LINC00857 from the TCGA and GTEx databases, we found that LINC00857 was significantly overexpressed in PC tissues (Fig. 1B), and this result was further confirmed in our centre (Fig. 1C). Moreover, TCGA data showed that LINC00857 expression levels were significantly associated with overall survival (OS) and progression-free survival (PFS) in PC patients (Fig. 1D, E). To verify the accuracy of LINC00857 expression for predicting prognosis, a ROC curve was constructed, and the results suggested that the area under the curve (AUC) of LINC00857 was 0.6238, which showed a good ability to predict prognosis. Subsequently, through PCR verification after nuclear-cytoplasmic separation, we found that LINC00857 was localized in both the nucleolus and cytoplasm (Fig. 1G, H). In addition, combined immunofluorescence and FISH assays confirmed this phenomenon (Fig. 1I). Finally, we compared the expression of LINC00857 in a normal pancreatic ductal epithelial cell line (hTERT-HPNE) and a PC cell line, and the results showed that LINC00857 expression was significantly increased in the PC cell line (Fig. 1J). Taken together, these results implied that LINC00857 was highly expressed in PC and associated with poor clinical prognosis.
3.2 LINC00857 promotes PC cell migration, invasion and epithelial mesenchymal transformation (EMT).
Considering the clinical association between LINC00857 and PC, we explored the biological significance of LINC00857 during the metastasis of PC cells in vitro. We first generated multiple PC cell lines with either silencing (by shRNAs) or overexpression of LINC00857, and the qRT–PCR results confirmed that LINC00857 expression was significantly decreased in Panc-1 and MIA Paca-2 cells stably expressing the shRNA (Fig. 2A). Additionally, LINC00857 overexpression was observed in PC cells stably expressing the pLVX-LINC00857 construct (Fig. 2B).
Notably, EMT is known as the key factor affecting PC cell metastasis. Next, we used western blotting to explore the effect of LINC00857 on modulating EMT-related molecules, and we found that LINC00857 knockdown with shRNAs inhibited the expression of N-cadherin/Vimentin and elevated the expression of E-cadherin, while forced expression of LINC00857 showed the logical inverse regulatory effects (Fig. 2C). In addition, we used wound healing and Transwell assays to evaluate the migration and invasion abilities of PC cells in different contexts. The wound healing assay showed that LINC00857 knockdown suppressed the migration of Panc-1 and MIA PaCa-2 cells (Fig. 2D), and the difference in the healed area was statistically significant (Fig. 2E), while overexpression of LINC00857 promoted migration (Fig. 2F, G). In addition, the Transwell assays demonstrated that lower LINC00857 expression significantly reduced the migration and invasion abilities of PC cells (Fig. 2H, I), while higher LINC00857 expression enhanced these abilities (Fig. 2J, K). These results revealed that LINC00857 facilitated the migration and invasion of PC cells.
3.3 LINC00857 promotes metastasis by enhancing FOXM1 expression in PC both in vitro and in vivo
Since the above results confirmed the role of LINC00857 in promoting the metastasis of PC cells, we expect to further explore the mechanism. Through website prediction, we found that LINC00857 had a high interaction probability with FOXM1, and the prediction scores were 0.9 and 0.97 with the RF classifier and SVM classifier, respectively (Table S1). Moreover. FOXM1 has been shown to promote the metastasis of PC in previous studies. Hence, we wondered whether LINC00857 could influence PC metastasis by regulating FXOM1. As Fig. 3A shows, FOXM1 expression was detected in Panc-1 and MIA Paca-2 cells with different LINC00857 expression levels. LINC00857 knockdown significantly inhibited FOXM1 protein expression, while overexpression of LINC00857 increased FOXM1 protein expression but not its mRNA level (Fig. S2A).
Then, we wanted to determine whether LINC00857's promotion of PC metastasis was dependent on FOXM1. si-FOXM1 was transfected into LINC00857-overexpressing cells, and EMT-related indicators were detected by western blotting. The results showed that N-cadherin, Vimentin and E-cadherin expression induced by LINC00857 overexpression could be significantly inhibited by si-FOXM1 (Fig. 3B). Moreover, FOXM1 knockdown partially rescued the LINC00857-induced increase in the PC cell wound healing ability (Fig. 3C, D). Subsequently, changes in PC cell migration and invasion were further examined, and the results suggested that LINC00857 promoted PC cell migration and invasion and that silencing FOXM1 evidently rescued this phenomenon in Panc-1 and MIA Paca-2 cells (Fig. 3E, F).
Next, we analysed the effect of the LINC00857/FOXM1 signalling pathway on tumour metastasis in vivo. Here, a xenograft metastasis model was established by injecting luciferase-expressing cells into BALB/c nude mice via the tail vein. The mice were randomly divided into three groups (A, B, C), and LINC00857-overexpressing KPC cells were subsequently injected into the mice in Groups B and C, while the mice in Group A was injected with control Panc-1 cells. Beginning on the third day after successful injection, the mice in Group C were given the FOXM1 inhibitor thiostrepton (17 mg/kg) by intraperitoneal injection two times a week, while the mice in the other groups were treated with saline (Fig. 3G). After three weeks of observation, small animal imaging technology was used to evaluate metastasis in the mice. By comparing Group A with Group B, we found that LINC00857-overexpressing PC cell injection resulted in higher fluorescence in lung tissues, and this fluorescence was reduced by administration of the FOXM1 inhibitor thiostrepton (Fig. 3H, I). In addition, mice injected with LINC00857-overexpressing cells (Group B) developed more metastatic lung nodules than control mice (Group A). Similarly, we found that the number of pulmonary metastatic nodules was also significantly reduced after thiostrepton treatment in the comparison between Groups B and C (Fig. J, K). Finally, the metastatic lung tissue samples were fixed and sectioned for staining. The HE staining results showed that the visual field area of lung metastasis was most extensive in Group B, followed by Groups A and C. IHC staining also confirmed that the percentage of FOXM1-positive cells in Group B was significantly higher than that in the other groups (Fig. 3L, M). Collectively, these results suggested that overexpression of LINC00857 enhanced the metastasis of PC, and the effect was inhibited by FOXM1 inhibitors both in vitro and in vivo.
3.4 LINC00857 stabilizes FOXM1 via OTUB1-mediated deubiquitination
Given that LINC00857 affected the FOXM1 protein level but not its mRNA level, we hypothesized that LINC00857 regulates FOXM1 through protein degradation. Common protein degradation methods always include the ubiquitin–proteasome pathway and the lysosomal pathway. Next, cycloheximide (CHX, an inhibitor of protein synthesis), MG132 (a proteasome inhibitor), and chloroquine (CQ, a lysosomal inhibitor) were applied in further experiments. We divided the PC cell interventions into three groups: CHX alone, CHX combined with MG132, and CHX combined with CQ. FOXM1 protein levels were analysed at 0, 1, 2 and 4 hours after intervention. The results showed that FOXM1 protein expression after MG132 treatment was higher than that after CQ treatment, which implied that FOXM1 was degraded mainly through the ubiquitin–proteasome pathway (Fig. S2C). Hence, we induced FOXM1 protein degradation using CHX in PC cells with different expression levels, and the results suggested that LINC00857 knockdown accelerated the degradation of the FOXM1 protein (Fig. 4A, B), while LINC00857 overexpression resulted in slower FOXM1 degradation (Fig. 4C, D). We then examined the ubiquitination level of FOXM1 and found after immunoprecipitation of endogenous FOXM1 in Panc-1 and MIA Paca-2 cells, obviously increased ubiquitin signals were detected in cells with stable LINC00857 silencing compared to the corresponding control cells. Consistent with this finding, FOXM1 ubiquitination was lower in cells overexpressing LINC00857 than in control cells (Fig. 4E).
Next, we explored how LINC00857 inhibits FOXM1 ubiquitination. Accumulating evidence has demonstrated that lncRNAs may function as scaffolds for binding proteins; thus, we hypothesized that LINC00857 recruits a deubiquitinase to bind FOXM1. To test this hypothesis, we performed RNA pull-down assays with PC cells and observed that multiple target proteins were pulled down by LINC00857, including the deubiquitinase OTUB1 (Fig. S3A). OTUB1 has been reported to block FOXM1 ubiquitination; thus, we further verified this finding. First, we found that OTUB1 knockdown reduced the FOXM1 protein level (Fig. 4F) but did not affect its mRNA expression level (Fig. S3B). Other experiments also confirmed that OTUB1 slowed FOXM1 degradation and reduced the FOXM1 ubiquitination level (Fig. S3C, S3D). Additionally, the western bloting results showed that LINC00857 knockdown did not affect the OTUB1 protein level (Fig. S3E), which means that LINC00857 does not regulate FOXM1 by altering OTUB1 expression but instead via another mechanism. Finally, we simultaneously knocked down OTUB1 in LINC00857-overexpressing cells, and the results showed that FOXM1 upregulation induced by LINC00857 was partially inhibited (Fig. 4E), and the decrease in FOXM1 ubiquitination induced by LINC00857 was also reversed by siOTUB1 (Fig. 4F). Collectively, the results suggested that LINC00857 may reduce ubiquitination-mediated degradation by recruiting the deubiquitinase OTUB1.
3.5 LINC00857 serves as a protein scaffold that promotes the interaction between FOXM1 and OTUB1
Accumulating evidence has demonstrated that lncRNAs may function in different ways during cancer development, including as scaffolds for protein interactions. Therefore, we hypothesized that LINC00857 may provide a scaffold for the interaction between OTUB1 and FOXM1, which increases FOXM1 deubiquitination. According to the previous RNA pull-down/mass spectrometry results, OTUB1 can be pulled down by LINC00857 (Fig. S3A). To verify this, we performed independent RNA pull-down assays in both Panc-1 and MIA PaCa-2 cells and found that OTUB1 was successfully pulled down by LINC00857 compared with its antisense RNA as the negative control (Fig. 5A). Furthermore, RIP assays suggested that LINC00857 was enriched in RNA–protein complexes precipitated with anti-OTUB1 antibody in PC cells (Fig. 5B). Combined immunofluorescence and FISH analysis showed that the colocalization of LINC00857 and OTUB1 was mainly in the nucleolus and partly in the cytoplasm (Fig. 5C). Similarly, we found that FOXM1 was pulled down by LINC00857 (Fig. 5D), and qRT–PCR results indicated obvious enrichment of LINC00857 mRNA with an anti-FOXM1 antibody (Fig. 5E). Moreover, we also confirmed the colocalization of LINC00857 and FOXM1 (Fig. 5F).
The above results indicated that LINC00857 can bind separately to FOXM1 and OTUB1. We then set out to demonstrate the effect of this binding on the interaction between FOXM1 and OTUB1. To this end, we first revealed the interaction between OTUB1 and FOXM1. Co-IP assays showed that FOXM1 could be precipitated with OTUB1 and that endogenous OTUB1 could be precipitated with FOXM1 in Panc-1 and MIA Paca-2 cells (Fig. 5G). Moreover, immunofluorescence assays demonstrated that FOXM1 and OTUB1 were colocalized with each other in PC cells (Fig. 5H). As shown in Fig. 5I, we then knocked down LINC00857, and the results suggested that less FOXM1/OTUB1 protein was precipitated with anti-OTUB1/FOXM1 antibodies from Panc-1 cells in comparison with the corresponding control cells. In contrast, more FOXM1/OTUB1 protein was precipitated with anti-OTUB1/FOXM1 antibodies in cells with stable forced expression of LINC00857 compared to the corresponding control cells (Fig. 5J), and similar results were confirmed in MIA Paca-2 cells (Fig. 5K).
Since we attempted to demonstrate the role of LINC00857 as a scaffold, we were eager to clarify the possible binding site. Then, we applied ChIRP, in which we used 20 biotin-labelled probes for LINC00857 segments to pull down FOXM1 and OTUB1separately (Fig. 5L). The results showed that compared with the control probe, probes 2 and 7 could pull down FOXM1, while OTUB1 could be pulled down by probes 3, 6 and 10 (Fig. 5M), which means that the binding of LINC00857 to FOXM1 may be localized between nt 422–441 or 1253–1272 and that the binding of LINC00857 to OTUB1 may be localized between nt 51–70, 401–420 or 721–740 (Fig. 5N). Collectively, these results suggested that LINC00857 may function as a scaffold to facilitate interactions between FOXM1 and OTUB1, which could decelerate FOXM1 degradation.
3.6 LINC00857 is transcriptionally regulated by mutant p53
As p53 mutation is one of the predisposing factors for the development of PC, we interestingly observed that the mutant p53 group expressed a higher level of LINC00857 than the wild-type p53 group in the TCGA database (Fig. 6A). Thus, we speculated that LINC00857 could be regulated by mutant p53. Panc-1 and MIA Paca-2, two p53 mutant cell lines, were convenient for further study. We knocked down and overexpressed mutant p53 in these two cell lines and verified the knockdown and overexpression efficiencies; the results showed that mutant p53 expression was successfully changed at the mRNA (Fig. 6B, C) and protein (Fig. 6D) levels. Indeed, depletion of mutant p53 significantly abolished LINC00857 expression (Fig. 6E). In contrast, the expression of LINC00857 was increased with mutant p53 overexpression (Fig. 6F). To confirm the possible binding site, we searched the JASPAR and PROMO databases. Sequence analysis showed that the LINC00857 promoter contains 5 putative binding sites for mutant p53 (Fig. 6G). Based on these binding sites, we designed 5 pairs of primers (P1-P5) and conducted ChIP–qPCR experiments, and the results revealed obvious enrichment of P1 and P5 but not P2-P4 by the anti-p53 antibody (Fig. 6H); this result was also confirmed by DNA gel electrophoresis in PC cells (Fig. 6I). Then, dual-luciferase reporter plasmids containing the wild-type (WT) or mutant (MUT) promoter sequence were designed and transfected into Panc-1 and MIA Paca-2 cells (Fig. 6J). The results showed that overexpression of mutant p53 significantly increased the luciferase activity in the WT2 group, while no obvious changes were observed in the MUT2, WT1 and MUT2 groups (Fig. 6K), which means that mutant p53 binding at LINC00857 promoter P4 promotes transcription. Therefore, these data strongly indicated that LINC00857 was transcriptionally regulated by mutant p53 in PC cells.
3.7 Mutant p53-LINC00857-mediated metastasis of PC was inhibited by atorvastatin
Statins are among the most commonly used lipid-lowering drugs in clinical practice, and recent studies have reported that statins can degrade mutant p53. Here, we applied one of the commonly used drugs in the clinic, ATOR (Fig. 7A), to treat cells and observed its regulatory effect on mutant p53 and LINC00857. Using a CCK-8 assay, we evaluated the viability of Panc-1, MIAPaCa-2 and BxPC-3 pancreatic cancer cells treated with different doses of ATOR. ATOR significantly reduced the growth of pancreatic cancer cells in a dose-dependent manner, and its half-maximal inhibitory concentrations (IC50) in Panc-1 and MIA PaCa-2 cells were 43.04 µM and 27.63 µM, respectively (Fig. 7B, 7C). Then, ATOR was applied to PC cells, and we found that ATOR significantly downregulated the protein expression of mutant p53 and the mRNA expression of LINC00857 (Fig. 7D, 7E). Previous data suggested that LINC00857-FOXM1 can promote metastasis in PC. Therefore, we further verified whether ATOR can inhibit PC metastasis, and the results indicated that ATOR elevated the expression of E-cadherin and reduced the expression of N-cadherin/Vimentin (Fig. 7F), which seemed to influence the EMT behaviour of PC cells. Moreover, through Transwell assays, we found that ATOR suppressed the migration and invasion of Panc-1 cells in a dose-dependent manner (Fig. 7G, 7H), and a similar effect was confirmed by a Transwell assay in MIA Paca-2 cells (Fig. 7I, 7J). Taken together, our data suggested that ATOR could inhibit p53/LINC00857 axis mutation-mediated pancreatic cancer cell metastasis.