LINC01559 is upregulated in LUAD tissues and cell lines
To gain a comprehensive and in-depth understanding of LINC01559 in pan-cancer, we retrieved the expression data of LINC01559 from TCGA and GTEx database. As depicted in Fig. 1A, the expression of LINC01559 was significantly upregulated in various cancers, including LUAD. Among the 58 pairs of LUAD tissues from TCGA, LINC01559 exhibited higher expression in tumor tissues compared to normal tissues (Fig. 1B), and the ROC curve (AUC = 0.81) underscored that LINC01559 hold significant diagnostic value for LUAD (Fig. 1C). Subsequent analysis revealed that heightened expression of LINC01559 was associated with an adverse prognosis among patients (Fig. 1D). To investigate whether LINC01559 was an independent factor affecting the prognosis of LUAD, we performed COX regression analysis on the Biomarker Exploration of Solid Tumors website and identified it as an independent prognostic factor in datasets including GSE26939, GSE13213, GSE72094, and TCGA (Fig. 1E). To delve deeper into the association between LINC01559 expression and M stage, we conducted an analysis of GSE11117 data on the Biomarker Exploration of Solid Tumors website and the result indicated a significant elevation of LINC01559 expression in M1 stage compared to M0 stage (Fig. 1F). To validate these findings, we employed four LUAD cell lines to egauged the relative expression of LINC01559 via qRT-PCR in comparison with the HSF cell lines. The result indicated that LINC01559 was upregulated in the A549 and H1299 cell lines (Fig. 1G).
LINC01559 is correlated with the histological grade of LUAD patients.
To investigate the expression of LINC01559, we employed TMAs consisting of 60 tissues, including 30 LUAD tissues and paired 30 para-cancer normal tissues. The expression level of LINC01559 was evaluated by ISH assay (Supplemental Fig.SA-B). Results indicated that the expression of LINC01559 in LUAD tissues was higher compared to normal tissues (Fig. 2A-B). Further analysis revealed a significant correlation between high level of LINC01559 expression and histological grade among LUAD patients (Fig. 2C).
LINC01559 enhances invasion and migration of LUAD cells
To elucidate the role of LINC01559 in LUAD cells, we employed siRNA-mediated knockdown and overexpression plasmid system respectively to downregulate and upregulate LINC01559 expression in A549 and H1299 cells (Fig. 3A-B). Furthermore, transwell assays were performed to determine the invasive and migratory abilities of LINC01559 in LUAD cells. Our results demonstrated that knockdown of LINC01559 significantly suppressed the invasion and migration of LUAD cells (Fig. 3C-D), while overexpression of LINC01559 had the opposite effect (Fig. 3E-F).
LINC01559 promoted LUAD metastasis in vivo
Furthermore, we performed in vivo assays to assess the effect of LINC01559 on LUAD metastasis. Initially, we established A549 sh-LINC01559 and A549 sh-NC cell lines to guarantee stable knockdown of LINC01559 expression (Fig. 4A). Subsequently, we proceeded with tail vein injection of LINC01559-knockdown A549 cells and control cells to establish tumor metastasis mouse models. We observed the weight of sh-NC group gradually decreased compared to the sh-LINC01559 group and the gap increased over time (Fig. 4B). Notably, the knockdown of LINC01559 led to a visible decrease in the number of lung metastasis lesions (Fig. 4C-D), and we measured that the lungs in the sh-NC group were heavier than those in the sh-LINC01559 group (Fig. 4E). Tumor metastasis was confirmed by HE staining (Fig. 4F). Therefore, both in vitro and in vivo data corroborated the metastatic effect of LINC01559 in LUAD.
LINC01559 directly interacts with VIM
Coding Potential Calculator version 2 (CPC2) is a widely used tool for predicting the coding potential of RNA transcripts, particularly long non-coding ones(27). By utilizing CPC2, we have enhanced our confidence in the biological functionality of LINC01559 based on its non-coding attributes (Fig. 5A). Furthermore, the results of subcellular fractionation assays indicated that LINC01559 was predominantly localized in the cytoplasm, suggesting its involvement in post-transcriptional regulatory processes (Fig. 5B). To further elucidate the molecular mechanism underlying LINC01559-mediated LUAD metastasis, we employed RNA-protein pull down assays to identify LINC01559-binding proteins. Coomassie brilliant blue staining and mass spectrometry revealed that VIM was the primary protein interacting with LINC01559 (Fig. 5C-D, Supplemental Fig.SC). Subsequently, RIP assays confirmed that VIM specifically binds to LINC01559 in LUAD cells (Fig. 5E). These findings suggest that LINC01559 actively contributes to LUAD metastasis by interacting with VIM.
LINC01559 regulates the protein level of VIM
To further investigate the potential regulatory role of LINC01559 on VIM expression, knockdown experiments targeting LINC01559 were conducted in A549 and H1299 cell lines, followed by assessment of mRNA levels for VIM. The results showed that there was no significant change in the mRNA expression level of VIM (Fig. 6A). Similar results were obtained after overexpression of LINC01559 (Fig. 6B). However, we observed alterations in VIM protein levels after knockdown or overexpression of LINC01559 in A549 and H1299 cell lines (Fig. 6C-D).
LINC01559-mediated metastasis promotion is facilitated by VIM.
To explore the potential role of VIM in promoting metastasis of LUAD, we initially evaluated the mRNA expression levels of VIM in LUAD cells and HSF cell lines. Our findings indicated that VIM was significantly upregulated in LUAD cell lines compared to HSF (Fig. 7A). Bioinformatics analysis suggested that higher expression of VIM was positively correlated with the higher stage of LUAD patients (Fig. 7B). Subsequently, transwell assays were conducted to evaluate the invasive and migratory capabilities of A549 and H1299 cell lines with respect to VIM. After knocking down VIM expression using siRNA(Fig. 7C), we observed that the invasion and migration capabilities of LUAD cells were significantly inhibited (Fig. 7D-E). Importantly, our findings suggested that the upregulation of VIM could rescue the reduced migration and invasion caused by knocking down LINC10559 (Fig. 8A-B), indicating that LINC10559-mediated acceleration of migration and invasion was mediated through VIM.
LINC01559 impedes the degradation of VIM by inhibiting its ubiquitination.
Based on the observation that LINC01559 only affects VIM protein levels, we used CHX to treat A549 and H1299 cells transfected with siLINC01559-1 or siNC and monitored the expression of VIM protein at 0, 3, 6, and 9 hours. As depicted in Fig. 9A, knockdown of LINC01559 expedited the degradation of VIM, indicating that LINC01559 protected VIM from degradation in LUAD cells. Moreover, MG132, a proteasome inhibitor, was able to alleviate the decrease of VIM caused by LINC01559 knockdown (Fig. 9B). Based on the aforementioned results, further examination was conducted on the ubiquitination of VIM protein, revealing that LINC01559 knockdown or overexpression led to a increase or reduction in the ubiquitination levels of the VIM protein (Fig. 9C-D). These findings suggest that LINC01559 plays a protective role in preventing proteasomal degradation of VIM.