LncRNA weighted gene co-expression network analysis reveals novel biomarkers related to prostate cancer metastasis

Most prostate cancer patients die from metastasis and lack accurate efficacious biomarkers to monitor the disease behavior, optimize treatment and assess prognosis. Herein, we aimed to identify meaningful lncRNA biomarkers associated with prostate cancer metastatic progression.

By repurposing microarray probes, 11,624 lncRNAs in prostate cancer were obtained from Gene Expression Omnibus  database (GSE46691, N = 545; GSE29079, N = 235; GSE94767, N = 130). Weighted gene co-expression network analysis was applied to determine the co-expression lncRNA network pertinent to metastasis. Hub lncRNAs were screened. RNA-seq and clinical data from the Cancer Genome Atlas prostate cancer (TCGA-PRAD) cohort (N = 531) were analyzed. Transwell assay and bioinformatic analysis were performed for mechanism research.

The high expression levels of nine hub lncRNAs (FTX, AC005261.1, NORAD, LINC01578, AC004542.2, ZFAS1, EBLN3P, THUMPD3-AS1, GAS5) were significantly associated with Gleason score and increased probability of metastatic progression. Among these lncRNAs, ZFAS1 had the consistent trends of expression in all of the analysis from different cohorts, and the Kaplan-Meier survival analyses showed higher expression of ZFAS1 was associated with shorter relapse free survival. In-vitro studies confirmed that downregulation of ZFAS1 decreased prostate cancer cell migration.

We offered some new insights into discovering lncRNA markers correlated with metastatic progression of prostate cancer using the WGCNA. Some may serve as potential prognostic biomarkers and therapeutic targets for advanced metastatic prostate cancer.