- Research
- Open access
- Published:
LncRNA weighted gene co-expression network analysis reveals novel biomarkers related to prostate cancer metastasis
BMC Medical Genomics volume 15, Article number: 256 (2022)
Abstract
Background
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.
Methods
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.
Results
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.
Conclusion
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.
Introduction
According to the Encyclopedia of DNA Elements (ENCODE) project, less than 2% of the human genome is translated into proteins. More than 80% is transcribed into a versatile group of RNA. Those RNA transcripts without protein-coding capability, known as non-coding RNA (ncRNA), can be divided into “housekeeping ncRNA” such as mRNA and “regulatory ncRNA” such as small ncRNAs or long ncRNAs (lncRNAs) [1]. However, the functional roles of lncRNAs (> 200 nucleotides) are still not nearly as well-known as small ncRNAs, particularly miRNA. It gradually becomes clear that lncRNAs play important roles in transcription, translation, splicing, nuclear/cytoplasmic trafficking and interact with major pathways controlling proliferation, differentiation, or survival [2]. Accumulating evidence in hematological and solid malignancies have shown that lncRNAs can regulate tumor initiation, progression, or metastasis due to alterations of gene expressions [3]. In addition, for their tissue and tumor specificity, lncRNAs can be used as promising candidates for biomarkers and therapeutic targets.
Some lncRNAs have been applied successfully in clinical practice. Prostate cancer antigen 3 (PCA3), the highly prostate-specific lncRNA overexpressed in more than 95% of primary prostate tumors, was developed as a urine test for the supplementary diagnosis of prostate cancer (PCa) [4]. Furthermore, the Progensa PCA3 test was successfully approved by the US Food and Drug Administration (FDA) in 2012 [5]. The lncRNA Second Chromosome Locus Associated with Prostate 1 (SChLAP1) has also been identified as a molecular driver and predictive biomarker of aggressive PCa in multiple clinical studies [6]. PCa is the foremost non-cutaneous malignancy in men worldwide, with increasing morbidity and mortality. The leading cause of prostate cancer death is metastasis. The five-year survival rate for metastatic prostate cancer (mPCa) dramatically decrease to only 31% while for localized PCa is almost 100%1. Worse, some of them exhibit rapid resistance to standard androgen deprivation therapy (ADT) combined with/without microtubule-targeted taxane chemotherapy and lead to fatal disease. It is important to study clonal diversity and cancer heterogeneity alongside the development and impact of the tumour microenvironment (TME), which includes the extracellular matrix (ECM),blood vessels, immune cells and fibroblasts2. However, the present clinical predictors (such as serum prostate-specific antigen, biopsy Gleason score) of biological outcomes for PCa still remain imprecise and imperfect [7]. Meanwhile, molecular biomarkers, especially non-coding genes, set a new phenotypic frontier in prostate cancer metastasis and resistance3. Therefore, we aimed to find the clinically relevant lncRNA network and hub lncRNAs tightly associated with PCa metastasis by applying systems biology.
The functions of lncRNAs are greatly influenced by their transcript abundance due to not encoding proteins. However, the overall expression abundance of lncRNA is lower than mRNA, which leads to higher requirements for the amount of sequencing data for the detection of lncRNA. In comparison, chips have higher reliability in detecting low-abundance RNA, particularly in lncRNA with overlapping regions of mRNA. Chip-based probes are shorter and can be designed according to specific sequences, which are more suitable for high-throughput lncRNA detection than RNA sequencing. Therefore, we developed a computational pipeline to re-annotate probes that are uniquely mapped to lncRNA.
Moreover, the Affymetrix Human Exon 1.0 ST Array platform was used to query the expression of lncRNA in PCa samples. This platform retains the most annotated lncRNA data currently [8]. Based on the expression data, weighted gene co-expression network analysis (WGCNA), which allows a systems-level investigation of correlations between genotype and phenotype, was chosen to examine how lncRNAs jointly affect prostate cancer metastasis. WGCNA approach can identify the highly correlated genes that share quite similar expression profiles on different samples and construct a “co-expression network” using these genes. WGCNA can cluster highly correlated genes into a module, a subregion of the network, and might be linked to specific functions for the members involved in well-defined pathways. The central genes in a module are defined as “hub genes” that are regarded as the most suitable candidates to imply some physiological or pathological processes [9, 10]. Comparing co-expression lncRNAs networks according to the metastatic and non-metastatic associated expression data can recognize the distinguished modules and the most relevant hub lncRNAs. Thus, some potential diagnostic or prognostic lncRNAs that are highly associated with PCa metastasis could be found.
Materials and methods
Data recruitment and preprocessing
PCa related gene expression data were searched in the Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/) of National Center of Biotechnology Information (NCBI), and were downloaded according to the accession numbers GSE46691(N = 545), GSE29079 (N = 235), GSE94767 (N = 130). The microarray platform used in these above datasets was Affymetrix Human Exon 1.0 ST array (HuEx-1_0-st). GSE46691 consisting of 545 PCa patients, 212 (39%) developed metastatic progression after long-term follow-up (mean 13.4 years), was used to construct the co-expression gene module and discover hub lncRNAs. GSE29079 and GSE94767 consisting of the tumor stage information, were exploited as independent validation datasets. These datasets were all labeled with their corresponding accession numbers. For the TCGA prostate adenocarcinoma (PRAD) cohort (N = 531), the gene expression profiles of the prostate cancer tumor and adjacent normal tissue and the most recent clinicopathological information of these patients were extracted from Genomic Data Commons Data Portal (https://portal.gdc.cancer.gov).
Re-purposing data to interrogate lncRNA expression
We collected lncRNA annotations from the Ensembl database (Homo sapiens GRCh38, release 96) [11] and re-annotated probe sets of the Human Exon 1.0 ST array (GPL5188 or GPL5175) for lncRNAs through mapping all probes to the human genome with SeqMap [12], just as the computational pipeline developed by Liu et al. [8]. Only the probes with no mismatch and mapped uniquely to the genome were kept. Based on the Ensembl database annotations (http://www.ensembl.org/), probes mapped to pseudogene transcripts or protein-coding transcripts were all removed. At last, we obtained 11,624 corresponding lncRNA genes with no less than 4 probes. The raw intensity of the array probes was preprocessed identically using the R package oligo for normalization and background correction. The lncRNA expression level was measured as the average value of the background-corrected intensity with all probes corresponding to this lncRNA.
WGCNA for lncRNA and hub lncRNAs definition
Four co-expression networks were constructed using the WGCNA package in R. LncRNA expression profiles of the GSE46691 dataset were used to conduct the co-expression networks of the selected lncRNAs. An unsupervised co-expression relationship was initially built on the adjacency matrix of connection strengths using Pearson’s correlation coefficients for all pair-wise lncRNAs. The parameter of the matrix emphasized the strong connections between lncRNAs, whereas it penalized weaker connections (β = 6). Then, the adjacency matrix has been transformed into a topological overlap matrix (TOM), which represented the connectivity of each lncRNA in the module. Hierarchical cluster analysis dendrogram was also performed to identify co-expression clusters with the gene profiler with a minimum gene number of 30.
The module eigengenes (MEs) were calculated following principal components analysis, representing the mean expression level of a module. The association between tumor metastasis (metastasis = 1; non-metastasis = 0) and MEs were tested with a logistic regression model. For single-lncRNA association analysis, each lncRNA expression was dichotomized around its median expression, and then logistic regression analysis was performed. Meanwhile, we test the association between the expression of each lncRNA and metastasis. Gene significance (GS) of each lncRNA, which represented the biological significance of the lncRNA for metastasis, was calculated as minus log 10 of the p values [13]. In our study, hub genes were selected by the following norm, top 10 genes with the largest k.in the module and GS larger than 1.3 (the p-value of the association test < 0.05). The significance of differences between groups was assessed by Student’s t-test or one-way ANOVA. P values of statistical significance are represented as ∗p < 0 05, ∗∗p < 0 01.
Survival analyses
To appraise the prognostic values of the hub lncRNAs, we explored the relationship between hub lncRNAs and survival by using the ‘survival’ R package on the TCGA PRAD dataset. Biochemical recurrence-free survival (RFS) was used as the survival endpoints respectively. Patients were divided into the metastatic group and the non-metastatic group according to the TNM stages defined by the AJCC 8th edition system. Each lncRNA was dichotomized into high and low expression subgroups around its median value in each group. RFS was determined and compared using the Kaplan–Meier and the log-rank (Mantel-Cox) test.
Cell lines and cell culture
The human prostate cancer cell lines PC-3 and LNCaP were purchased from the Cell Resource Center of Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences (Shanghai, PR China). The cells were cultured according to the instructions from American Type Culture Collection (ATCC). PC-3 and LNCaP cells were cultured in F-12 K medium (Biological Industries) and RPMI-1640 medium (Biological Industries), respectively, supplemented with 10% (v/v) fetal bovine serum (Biological Industries), 100 units/mL penicillin, and 100 mg/L streptomycin (Invitrogen).
Plasmid and RNA interference
The siRNA targeting ZFAS1 (si-ZFAS1) was obtained from Ribobio (Guangzhou, PR China), while as control, the mammalian non-targeting siRNA (si-NC) was used. Appropriate siRNA oligos for ZFAS1 were transfected into PC-3 or LNCaP cells by using Lipofectamine RNAiMAX (Invitrogen) according to the manufacturer’s protocol. Cells were harvested and examined for PCR analysis to ensure the efficiency of siRNA treatment after 24 h transfection or harvested for transwell assay after 48 h transfection. Detailed information about the sequences of siRNA could be found in Additional file 1: Table S1.
RNA extraction and real-time PCR analysis
Real-time quantitive RT-PCR was carried out to confirm the inhibition of lncRNA expression. Briefly, total RNA was isolated with the TRIzol reagent (Invitrogen) and reverse transcribed using PrimeScript RT reagent Kit with gDNA Eraser (TaKaRa) according to the manufacturer’s instructions. Then the PCR reaction was conducted at 95 °C for 2 min followed by 40 cycles at 95 °C for 15 s and 61 °C for 20 s using the SYBR® Premix Ex TaqTM II real-time PCR system (TaKaRa). All gene expression levels in different treatments were represented relative to their relevant control (Ct) and normalized to GAPDH gene levels (ΔΔCt). All assays were repeated at least three times. The final data were presented as mean ± SEM, and p < 0.05 was considered a statistical difference. The sequences of primers were listed in Additional file 1: Table S2.
Transwell assay
For transwell assay, 2 × 104 cells of each group were re-suspended in 200ml serum-free medium and planted in the upper chamber (8.0 mm pore, BD, USA), and 600ml of appropriate medium with 10% FBS was added to the lower chamber. After incubating for 24 h at 37 ℃, the medium in the upper chamber was removed. The attached cells in the lower section were fixed with 4% ice-cold polymethanol for 15 min before stained with 0.1% crystal violet for 30 min. The cells that migrated to the reverse side of the upper chambers were photographed. The migration rate was calculated by counting the migration cells in five random fields. All assays were repeated at least three times. The final data were presented as mean ± SEM, and p < 0.05 was considered a statistical difference.
LncRNA–miRNA–mRNA network prediction
For predicting potential lncRNA–target regulations, we used LncTarD (http://bio-bigdata.hrbmu.edu.cn/LncTarD/), a manually curated database that offered a comprehensive resource of experimentally supported lncRNA-target regulations and integrated pan-cancer transcriptome data from TCGA[14]. The experimental module from DIANA-LncBase v2 (www.microrna.gr/LncBase), which provided a broad compendium of experimentally supported microRNA-lncRNA interactions from high-throughput and low yield methodologies extracted from published literature and analyzed raw sequencing data[15], was used to predict lncRNA-microRNA interactions (validation type was Direct). MiRcode (http://www.mircode.org)[16], which provides putative miRNA target sites based on the complete GENCODE gene annotation, and ctcRbase (http://www.origin-gene.cn/database/ctcRbase/), a database collected expression data of circulating tumor cells/microemboli (CTCs/CTM) in all kinds of cancer types derived from RNA-seq data analyses [17], were used to predict miRNA–mRNA target relationships. Cytoscape_3.7.2 software was used to visualize networks.
Results
LncRNA modules association with PCa metastasis
Totally, 4 lncRNA modules were identified, labeled with blue, green, turquoise, or yellow (Fig. 1 A, B), which indicated a high degree of independence among the 4 modules and relative independence of lncRNAs expression in each module. Similar results were shown by the heatmap plotted according to the correlation between the 4 modules (Fig. 1 C). The yellow module (OR = 1.44, 95% CI = 1.02–2.04, p = 3.83 × 10− 2) was significantly associated with PCa metastasis (Table 1). The MEs in the yellow module was calculated and revealed a positively correlation with tumor metastasis status (p = 0.0039) as well as the Gleason score(p = 1.9 × 10− 4) (Fig. 1D, E).
Hub lncRNAs definition and correlation with metastasis
9 lncRNAs (FTX, AC005261.1, NORAD, LINC01578, AC004542.2, ZFAS1, EBLN3P, THUMPD3-AS1, GAS5) in the yellow module met the criteria (Table 2). The Results indicated that the expression of all the 9 hub lncRNAs was up-regulated in metastatic tumor group relative to non-metastatic tumor group (Fig. 2).
Altered expression of hub lncRNAs throughout prostate tumor progression
According to GSE29079 dataset, the samples were divided into the neoplasm and benign groups, T-test analysis showed that 5 lncRNAs (AC005261.1 p = 2.16 × 10− 13, AC004542.2 p = 6.40 × 10− 17, ZFAS1 p = 1.06 × 10− 6, EBLN3P p = 8.09 × 10− 8, GAS5 p = 0.014) had been up-regulated in the tumor samples, while other 4 lncRNAs were not significant (Fig. 3). Furthermore, the one-way ANOVA analysis showed seven lncRNAs increased with tumor progression (NORAD p = 0.025, LINC01578 p = 0.023, AC004542.2 p = 0.025, ZFAS1 p = 1.17 × 10− 8, EBLN3P p = 0.002, THUMPD3-AS1 p = 0.044, GAS5 p = 2.26 × 10− 7), which implied the expression levels of these lncRNAs had a highly consistent variation tendency with the grades of PCa (Fig. 4). Taken together, the lncRNAs in common (EBLN3P, AC004542.2, ZFAS1, and GAS5) seemed the most potentially predictive or prognostic genes for PCa that were closely related to the tumor progression.
Hub lncRNAs validation using TCGA data
The results showed that AC005261.1, LINC01578, ZFAS1, EBLN3P, THUMPD3-AS1 and GAS5 were highly expressed in tumor than adjacent normal samples (Additional file 2: Fig. S1). Differential expressions of FTX and NORAD were not significant. However, when the patients were stratified into metastatic and non-metastatic groups according to TNM classification (N > 0 or M > 0), the expressions of FTX, AC005261.1, NORAD, ZFAS1, EBLN3P and THUMPD3-AS1 were significantly higher in metastatic group (Fig. 5 A). Interestingly, when the patients were classified into two categories determined by tumor extra-prostatic extension or Gleason score > 8, Stage I-IIIA and Stage IIIB-IV, only LINC01578 and GAS5 were not significant higher expressions in high-grade tumors (Additional file 3: Fig. 2 A).
The curves revealed that higher expression of almost each hub lncRNA in non-metastatic group was connected to shorter RFS. Conversely, lower expression of almost each of them in metastatic group was closely correlated with worse survival (Fig. 5B). But unexpectedly, the expressions of these hub lncRNAs except FTX seem unable to predict the RFS in advanced patients when the patients were classified into Stage I-IIIA and Stage IIIB-IV group (Additional file 3: Fig. 2B). These results implied the 8 hub lncRNAs might be involved in the pathogenesis of tumor metastasis, while their functions were likely to change once the tumor had metastasized.
Downregulation of ZFAS1 inhibited the migration of prostate cancer cells
In view of above findings, the expression of ZFAS1 was consistent in all different independent cohorts. To further verify the its effect on prostate cancer metastasis, the mobility of cells was measured by the transwell experiment. Downregulation of ZFAS1 showed a decreased migratory capacity of prostate cancer cells (Fig. 6). Compared with negative control cells, there was a significant 55%, or 54% decrease in LNCaP cells or PC-3 cells transfected with si-ZFAS1, respectively. The changes in migratory cellular behavior impled that downregulation of ZFAS1 in prostate cancer might cause a less aggressive phenotype consistent with the clinical data analysis.
lncRNA–miRNA–mRNA ceRNA predicted regulatory networks
Based on the LncTarD records, we obtained the target networks of 4 hub lncRNAs (GAS, ZFAS1, NORAD, and FTX) (Fig. 7 A). It was noted that ZEB2 appeared to be not only co-regulated by NORAD and ZFAS1 but also involved in cell migration. The expression correlations of ZFAS1-ZEB2 and NORAD-ZEB2 across 33 TCGA cancer types were checked (Fig. 7B). Interestingly, ZFAS1 was negatively correlated with ZEB2, while NORAD was positively correlated with ZEB2 in PRAD cancer type. For mechanism research, the relevant lncRNA–miRNA–mRNA networks were constructed depending on the data from DIANA-LncBase, miRcode, and ctcRbase (Fig. 7 C). Among the potential miRNA targets, ZFAS1, NORAD, and ZEB2 (has-miR-144-3p, has-miR-101-3p, has-miR-26a-5p, has-miR-26b-5p) might be worthy of further verification.
Discussion
LncRNAs have emerged as new important regulators in the pathogenesis of human urologic cancers over the last decade. Herein, we undertook an unbiased study to identify lncRNA biomarkers associated with prostate cancer metastatic progression. Microarray/RNA-seq and clinical data of prostate cancer were obtained from the GEO and TCGA-PRAD cohort. The WGCNA was used to identify lncRNA networks in GSE46691. One lncRNA module was confirmed to be associated with cancer metastasis, and nine hub lncRNAs, namely FTX, AC005261.1, NORAD, LINC01578, AC004542.2, ZFAS1, EBLN3P, THUMPD3-AS1, and GAS5, were discovered. These nine lncRNAs’ expressions were all up-regulated in metastasis samples and significantly associated with the Gleason score. Among them, ZFAS1was not only high expressed in tumor groups compared with adjacent normal/benign groups, but significantly increased with deterioration by verifying in all the cohorts. More potentially, survival analysis suggested that early PCa patients with high expression of ZFAS1 might have a worse prognosis. Taken together, ZFAS1 appeared to be involved in the development of tumor metastasis. The in vitro analysis further confirmed that down-regulation of ZFAS1 expression decreased prostate cancer cell migration. In addition, the hub lncRNA-target networks were constructed, which implied that ZFAS1 and NORAD might participate in the process of cell migration through their common target ZEB2.
Metastatic PCa remains a fatal condition when cancer cells spread to the pelvic, bladder, rectum, bone, retroperitoneal lymph nodes, spinal cord, and other body areas, despite localized PCa has a high long-term survival [18]. Metastatic PCa may fall into two broad classifications, disease treated with no androgen deprivation and disease resistant to such therapy. The standard of care for metastatic patients remains androgen-deprivation therapy, which uses gonadotropin-releasing hormone (GnRH) analogs as an alternative to surgical castration [19]. The diagnosis of the metastatic disease depends on imaging, including radionuclide bone scanning, computed tomography (CT), positron emission tomography (PET) with choline or prostate-specific membrane antigen (PSMA), whole-body magnetic resonance imaging (MRI). However, the pathological confirmation from imaging studies is incomplete and not sufficiently sensitive [20]. Second, While androgen deprivation therapy (ADT) combined with microtubule-targeted taxane chemotherapy can provide a survival benefit in recurrent or metastatic disease, treatment resistance invariably develops, leading to fatal disease. Acquired drug resistance is the result of genetic/epigenetic changes that confer a drug-resistant phenotype in cancer cells4. LncRNA mediated mechanisms have been associated with epigenetic changes in prostate cancer. According to the results of the COSMIC-021, combination of the cabozantinib, a kind of antiangiogenic drugs, with the immune checkpoint inhibitor atezolizumab achieved encouraging activity in patients with metastatic castration-resistant PC (mCRPC), while multiple previous phase III clinical trials5,6 in refractory CRPC have demonstrated increased toxicity with no clinical benefit. Tumor angiogenesis is regulated through interaction between anti-angiogenic factors and pro-angiogenic within the tumor microenvironment. A growing number of studies have shown that different lncRNAs affect tumor angiogenesis, for example, lncRNA H19 targets miR-199a-5p and negatively regulates VEGFA 7. The vast spectrum of cancer traits warrants better molecular biomarkers to identify the more aggressive and clinically significant tumor subtypes. Accumulating recent evidence has shown that lncRNAs play key roles in cancer initiation or progression and appear to have many features for diagnostic markers or therapeutic targets [5].
Similar to PCA3, the clinically-relevant lncRNAs in PCa continue to be discovered. Studies in two separate independent cohorts validated prostate cancer-associated transcript-14 (PCAT-14), an androgen-regulated lncRNA, might be used as a novel PCa diagnosis biomarker for its high expression prostate tumors and low expression linked to poor outcomes [21, 22]. Studies found that prostate cancer gene expression marker 1 (PCGEM1) with highly prostate-specificity was also related to androgen receptor (AR) signaling and was overexpressed in therapy-resistant PCa [23, 24]. Metastasis-associated lung adenocarcinoma transcript-1 (MALAT1) reported a close connection between high expression and poor prognosis indicators[25, 26]. Similar findings had been made with the second chromosome locus associated with prostate-1 (SCHLAP1) [27]. Other lncRNAs clinically relevant to PCa were such as urothelial carcinoma-associated 1 (UCA1) [28, 29], nuclear enriched abundant transcript 1 (NEAT1) [30], HLA complex group 11 (HCG11) [31]. Nevertheless, only a few lncRNAs related to PCa metastasis have been identified so far. Thus, we performed WGCNA by using publicly available data to characterize the metastasis-associated lncRNAs.
WGCNA, also known as the weighted correlation network analysis, is widely used in bioinformatics applications based on pairwise correlations amongst variables [9]. It has been extensively used in various genomic applications. WGCNA can also be used as a data reduction technique, clustering method, feature selection method, framework for integrating genomic data, and data exploratory technique [32]. According to the theory, genes involved in similar pathways or closely related functions may have similar expression profiles. The WGCNA approach as a systems biology strategy transforms gene expression data across samples into a co-expression module that can provide perception in signaling networks. The results from the GSE46691 dataset showed that nine hub lncRNAs had noticeable associations for metastasis in the yellow module by the WGCNA package tool. Though the original array design did not target lncRNA measurement, microarray probes could be re-annotated for interrogating the expression of lncRNA. Especially comparing with RNA-seq data, the array-based expression data had better detection sensitivity and lower technical variation due to the low abundance of lncRNA [32]. Expression analyses from the other two independent datasets (GSE29079, GSE94767) showed four (AC004542.2, ZFAS1, EBLN3P, GAS5) out of nine lncRNAs were significantly overexpressed in tumor. Moreover, the TCGA data also showed again that ZFAS1, EBLN3P, or GAS5 was overexpressed in PCa tumor compared to adjacent normal samples. Besides, ZFAS1 showed the potential to be a predictor of poor prognosis.
For further mechanism research, the hub lncRNA-target networks based on experimentally supported functional regulations and expression associations in human diseases from LncTarD were constructed. ZEB2 was found to be a common potential target for both ZFAS1 and NORAD. The expression of ZEB2 in the TCGA-PRAD dataset was negatively correlated with ZFAS1 and positively correlated with NORAD, which indicated that the “co-highly expressed” ZFAS1 and NORAD had more complicated molecular mechanisms for regulating cell metastasis. ZEB2, a key activator of epithelial-mesenchymal transition (EMT), has been found to promote EMT processes in various cancer metastases [33]. Moreover, NORAD and ZFAS1 have been reported to promote cell invasion and migration in cervical cancer and bladder cancer, respectively, by up-regulating ZEB2 [34, 35]. Detaching from the tumor mass and becoming motile are necessary for metastasizing cancer cells. Once these cells intravasate into the lymphatic or blood circulation, the survivors among them will attach to the target organ endothelium, extravasate into the organ parenchyma, and proliferate. The abnormal activation of an EMT contributes to cancer metastasis, so treatments targeting the EMT pathway have become an attractive strategy [36]. In order to investigate the regulatory mechanism among NORAD, ZFAS1 and ZEB2, DIANA-LncBase v2, miRcode, and ctcRbase databases were used to discover their respective potential target miRNAs. The results showed 69 common potential miRNA targets between ZFAS1 and NORAD and 28 common potential miRNA targets between NORAD and ZEB2. The four miRNAs of has-miR-144-3p, has-miR-101-3p, has-miR-26a-5p, and has-miR-26b-5p were potential targets shared by ZFAS1, NORAD, and ZEB2. There have been no reports about the involvement of these four miRNAs in tumor migration or invasion. As we know, multiple targets could be regulated by one miRNA, while one target could also be co-regulated by multiple miRNAs. Similarly, one lncRNA could target more than one miRNA, and some lncRNAs could form a complex lncRNA-miRNA-mRNA regulatory network [37]. In addition, the regulation of multiple lncRNAs on the same target gene may have differences in their response time.
Indeed, ZFAS1 was up-regulated in metastasis samples and related to poor RFS [38]. High expression of ZFAS1 was closely associated with worse disease-free survival for PCa or gliomas patients [39, 40]. Recent studies disclosed that miR-135a-5 and miR-150-5p could be used as binding targets of ZFAS1 to participate in the regulation of tumor proliferation, invasion, and migration [41, 42]. As an important mechanism of lncRNA function, some lncRNAs with miRNA-binding sites can serve as ‘sponges’ to seclude endogenous miRNAs, which disturbs the modulation of gene expression mediated by miRNA. This phenomenon called competing endogenous RNAs (ceRNAs) belongs to the genome-wide finetuning regulatory molecules network [43,44,45]. PTEN pseudogene (PTENP1) is one of the ceRNAs correlated with PCa, which is linked to PTEN phosphatase as a tumor-suppressive pseudogene. Because of the 3’ UTR of PTENP1 RNA binding the same region of regulatory miRNA sequences, PTENP1 could reduce the downregulation of PTEN mRNA [46]. At present, the research on the mechanism of ZFAS1 as ceRNA in migration is still in progress.
NORAD, a newly identified lncRNA, is unique for its high conservatism and the key role in maintaining chromosomal stability by modulating the Pumilio (PUM) proteins’ activity. Moreover, genome instability drives metastasis via a cytosolic DNA response [47]. Recent studies indicated that NORAD was highly expressed in numerous human cancers, including breast cancer, esophageal squamous cell carcinoma, pancreatic cancer, colorectal cancer, bladder cancer, and cervical cancer, with poor overall survival [34, 48,49,50]. However, in hepatocellular carcinoma, NORAD acts contradictorily as a tumor suppressor [51]. In the present research, our results suggested that NORAD might involve cell migration and be highly expressed in the metastasis community of PCa. Significantly, regulation of NORAD affects the PUM proteins and changes the extracellular vesicle (EV) proteins that participate in communication between cells in the tumor microenvironment [52, 53]. The effects of NORAD are implicated in almost all aspects of tumors, including carcinogenesis proliferation, apoptosis, invasion, and metastasis [54]. However, the precise molecular mechanisms of NORAD are still at a preliminary stage that requires further systematic investigation.
In summary, our analyses focused on the metastasis of PCa and uncovered some relevant hub lncRNAs. Although other lncRNAs did not show significant differences in all cohort like ZFAS1 that might be related to the different stages of the disease and the different types of samples we used for analysis. The other 8 hub lncRNAs still shown some potential as a biomarker correlated with metastatic progression of PCa. For ZFAS1, further molecular mechanism investigations are needed in order to understand its function adequately.
Conclusion
To conclude, the yellow module derived from the microarray-based dataset via WGCNA was significantly associated with metastasis of PCa. The nine hub lncRNAs within yellow were successfully validated in different clinical cohots, especially ZFAS1. These results suggest some new lncRNAs as potential prognostic biomarkers and therapeutic targets for prostate cancer development and progression worthy of further investigation.
Availability of data and materials
Publicly available datasets were analyzed in this study. All datasets for this study are included in the article. The authors acknowledge the efforts of all of the researchers who have contributed the data to the public databases of GEO (https://www.ncbi.nlm.nih.gov/geo/) with accession numbers GSE46691; GSE29079 and GSE94767, and the cancer genome database TCGA (https://www.cancer.gov/about-nci/organization/ccg/research/structural-genomics/tcga). The interpretation and reporting of these data are the sole responsibility of the authors.
References
Jiang C, Li X, Zhao H, Liu H. Long non-coding RNAs: potential new biomarkers for predicting tumor invasion and metastasis. Mol Cancer. 2016;15(1):62.
Guttman M, Amit I, Garber M, French C, Lin MF, Feldser D, Huarte M, Zuk O, Carey BW, Cassady JP, et al. Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals. Nature. 2009;458(7235):223–7.
Mouraviev V, Lee B, Patel V, Albala D, Johansen TE, Partin A, Ross A, Perera RJ. Clinical prospects of long noncoding RNAs as novel biomarkers and therapeutic targets in prostate cancer. Prostate Cancer Prostatic Dis. 2016;19(1):14–20.
Ghafouri-Fard S, Khoshbakht T, Hussen BM, Baniahmad A, Taheri M, Rashnoo F. A review on the role of PCA3 lncRNA in carcinogenesis with an especial focus on prostate cancer. Pathol Res Pract. 2022;231:153800.
Bolton EM, Tuzova AV, Walsh AL, Lynch T, Perry AS. Noncoding RNAs in prostate cancer: the long and the short of it. Clin Cancer Res. 2014;20(1):35–43.
Falese JP, McFadden EJ, Hargrove AE. Structure-function relationships for the lncRNA SChLAP1 in aggressive prostate cancer. FASEB J. 2022. https://doi.org/10.1096/fasebj.2022.36.S1.0R251.
Prensner JR, Zhao S, Erho N, Schipper M, Iyer MK, Dhanasekaran SM, Magi-Galluzzi C, Mehra R, Sahu A, Siddiqui J, et al.RNA biomarkers associated with metastatic progression in prostate cancer: a multi institutional high throughput analysis of SChLAP1. Lancet Oncol. 2014;15(13):1469–80.
Du Z, Fei T, Verhaak RG, Su Z, Zhang Y, Brown M, Chen Y, Liu XS. Integrative genomic analyses reveal clinically relevant long noncoding RNAs in human cancer. Nat Struct Mol Biol. 2013;20(7):908–13.
Langfelder P, Horvath S. WGCNA. An R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.
Giulietti M, Righetti A, Principato G, Piva F. LncRNA co-expression network analysis reveals novel biomarkers for pancreatic cancer. Carcinogenesis. 2018;39(8):1016–25.
Zerbino DR, Achuthan P, Akanni W, Amode MR, Barrell D, Bhai J, Billis K, Cummins C, Gall A, Girón CG, et al. Ensembl 2018. Nucleic Acids Res. 2018;46(D1):D754-61.
Jiang H, Wong WH. SeqMap: mapping massive amount of oligonucleotides to the genome. Bioinformatics. 2008;24(20):2395–6.
Chang X, Shi L, Gao F, Russin J, Zeng L, He S, Chen TC, Giannotta SL, Weisenberger DJ, Zada G, et al. Genomic and transcriptome analysis revealing an oncogenic functional module in meningiomas. Neurosurg Focus. 2013;35(6):E3.
Zhao H, Shi J, Zhang Y, Xie A, Yu L, Zhang C, Lei J, Xu H, Leng Z, Li T, et al. LncTarD. A manually-curated database of experimentally-supported functional lncRNA-target regulations in human diseases. Nucleic Acids Res. 2020;48(D1):D118-26.
Paraskevopoulou MD, Vlachos IS, Karagkouni D, Georgakilas G, Kanellos I, Vergoulis T, Zagganas K, Tsanakas P, Floros E, Dalamagas T, et al. DIANA-LncBase v2: indexing microRNA targets on non-coding transcripts. Nucleic Acids Res. 2016;44(D1):D231-8.
Jeggari A, Marks DS, Larsson E. miRcode: a map of putative microRNA target sites in the long non-coding transcriptome. Bioinformatics. 2012;28(15):2062–3.
Zhao L, Wu X, Li T, Luo J, Dong D. ctcRbase: the gene expression database of circulating tumor cells and microemboli. Database. 2020;2020:baaa020. https://doi.org/10.1093/database/baaa020.
Schatten H. Brief overview of prostate cancer statistics, grading, diagnosis and treatment strategies. Adv Exp Med Biol. 2018;1095:1–14.
Sun M, Choueiri TK, Hamnvik OP, Preston MA, De Velasco G, Jiang W, Loeb S, Nguyen PL, Trinh QD. Comparison of gonadotropin-releasing hormone agonists and orchiectomy: effects of androgen-deprivation therapy. JAMA Oncol. 2016;2(4):500–7.
Sartor O, de Bono JS. Metastatic prostate cancer. New Engl J Med. 2018;378(17):1653–4.
White NM, Zhao SG, Zhang J, Rozycki EB, Dang HX, McFadden SD, Eteleeb AM, Alshalalfa M, Vergara IA, Erho N, et al. Eur Urol. 2017;71(2):257–66.
Shukla S, Zhang X, Niknafs YS, Xiao L, Mehra R, Cieslik M, Ross A, Schaeffer E, Malik B, Guo S, et al. Identification and validation of PCAT14 as prognostic biomarker in prostate cancer. Neoplasia. 2016;18(8):489–99.
Srikantan V, Zou Z, Petrovics G, Xu L, Augustus M, Davis L, Livezey JR, Connell T, Sesterhenn IA, Yoshino K, et al. PCGEM1, a prostate-specific gene, is overexpressed in prostate cancer. P Natl Acad Sci 2000;97(22):12216–21.
Petrovics G, Zhang W, Makarem M, Street JP, Connelly R, Sun L, Sesterhenn IA, Srikantan V, Moul JW, Srivastava S. Elevated expression of PCGEM1, a prostate-specific gene with cell growth-promoting function, is associated with high-risk prostate cancer patients. Oncogene. 2004;23(2):605–11.
Wang R, Sun Y, Li L, Niu Y, Lin W, Lin C, Antonarakis ES, Luo J, Yeh S, Chang C. Preclinical study using malat1 small interfering RNA or androgen receptor splicing variant 7 degradation enhancer ASC-J9((R)) to suppress enzalutamide-resistant prostate cancer progression. Eur Urol. 2017;72(5):835–44.
Wang F, Ren S, Chen R, Lu J, Shi X, Zhu Y, Zhang W, Jing T, Zhang C, Shen J, et al. Development and prospective multicenter evaluation of the long noncoding RNA MALAT-1 as a diagnostic urinary biomarker for prostate cancer. Oncotarget. 2014;5(22):11091–102.
Prensner JR, Iyer MK, Sahu A, Asangani IA, Cao Q, Patel L, Vergara IA, Davicioni E, Erho N, Ghadessi M, et al. The long noncoding RNA SChLAP1 promotes aggressive prostate cancer and antagonizes the SWI/SNF complex. Nat Genet. 2013;45(11):1392–8.
Yu Y, Gao F, He Q, Li G, Ding G. lncRNA UCA1 Functions as a ceRNA to promote prostate cancer progression via sponging miR143. Mol Ther Nucleic Acids. 2020;19:751–8.
Fotouhi Ghiam A, Taeb S, Huang X, Huang V, Ray J, Scarcello S, Hoey C, Jahangiri S, Fokas E, Loblaw A, et al. Long non-coding RNA urothelial carcinoma associated 1 (UCA1) mediates radiation response in prostate cancer. Oncotarget. 2017;8(3):4668–89.
Chakravarty D, Sboner A, Nair SS, Giannopoulou E, Li R, Hennig S, Mosquera JM, Pauwels J, Park K, Kossai M, et al. The oestrogen receptor alpha-regulated lncRNA NEAT1 is a critical modulator of prostate cancer. Nat Commun. 2014;5:5383.
Zhang Y, Zhang P, Wan X, Su X, Kong Z, Zhai Q, Xiang X, Li L, Li Y. Downregulation of long non-coding RNA HCG11 predicts a poor prognosis in prostate cancer. Biomed Pharmacother. 2016;83:936–41.
Raghavachari N, Barb J, Yang Y, Liu P, Woodhouse K, Levy D, O’Donnell CJ, Munson PJ, Kato GJ. A systematic comparison and evaluation of high density exon arrays and RNA-seq technology used to unravel the peripheral blood transcriptome of sickle cell disease. BMC Med Genom. 2012;5:28.
Sathyanarayanan A, Chandrasekaran KS, Karunagaran D. microRNA-145 modulates epithelial-mesenchymal transition and suppresses proliferation, migration and invasion by targeting SIP1 in human cervical cancer cells. Cell Oncol. 2017;40(2):119–31.
Huo H, Tian J, Wang R, Li Y, Qu C, Wang N. Long non-coding RNA NORAD upregulate SIP1 expression to promote cell proliferation and invasion in cervical cancer. Biomed Pharmacother. 2018;106:1454–60.
Yang H, Li G, Cheng B, Jiang R. ZFAS1 functions as an oncogenic long non-coding RNA in bladder cancer. Biosci Rep. 2018;38(3):BSR20180475. https://doi.org/10.1042/BSR20180475.
Mittal V. Epithelial mesenchymal transition in tumor metastasis. Annu Rev Pathol. 2018;13:395–412.
Liang H, Pan Z, Zhao X, Liu L, Sun J, Su X, Xu C, Zhou Y, Zhao D, Xu B, et al. LncRNA PFL contributes to cardiac fibrosis by acting as a competing endogenous RNA of let-7d. Theranostics. 2018;8(4):1180–94.
Ghafouri-Fard S, Kamali MJ, Abak A, Shoorei H, Taheri M. LncRNA ZFAS1: role in tumorigenesis and other diseases. Biomed Pharmacother. 2021;142:111999.
Chen X, Yang C, Xie S, Cheung E. Long non-coding RNA GAS5 and ZFAS1 are prognostic markers involved in translation targeted by miR-940 in prostate cancer. Oncotarget. 2018;9(1):1048–62.
Lv QL, Chen SH, Zhang X, Sun B, Hu L, Qu Q, Huang YT, Wang GH, Liu YL, Zhang YY, et al. Upregulation of long noncoding RNA zinc finger antisense 1 enhances epithelial-mesenchymal transition in vitro and predicts poor prognosis in glioma. Tumour Biol. 2017;39(3):1010428317695022. https://doi.org/10.1177/1010428317695022.
Xia B, Hou Y, Chen H, Yang S, Liu T, Lin M, Lou G. Long non-coding RNA ZFAS1 interacts with mir-150-5p to regulate Sp1 expression and ovarian cancer cell malignancy. Oncotarget. 2017;8(12):19534–46.
Pan J, Xu X, Wang G. lncRNA ZFAS1 is involved in the Proliferation, Invasion and Metastasis of prostate Cancer cells through competitively binding to miR-135a-5p. Cancer Manag Res. 2020;12:1135–49.
Salmena L, Poliseno L, Tay Y, Kats L. A ceRNA hypothesis: the rosetta stone of a hidden RNA language? Cell. 2011;146(3):353–8.
Su Y, Hou W, Zhang C, Ji P, Hu R, Zhang Q, Wang Y, Li P, Zhang H, Chen Y, et al. Long non-coding RNA ZFAS1 regulates cell proliferation and invasion in cervical cancer via the miR-190a-3p/KLF6 axis. Bioengineered. 2022;13(2):3840–51.
Wang H, Chen Y, Liu Y, Li Q, Luo J, Wang L, Chen Y, Sang C, Zhang W, Ge X, et al. The lncRNA ZFAS1 regulates lipogenesis in colorectal cancer by binding polyadenylate-binding protein 2 to stabilize SREBP1 mRNA. Mol Ther Nucleic Acids. 2022;27:363–74.
Poliseno L, Salmena L, Zhang J, Carver B, Haveman WJ, Pandolfi PP. A coding-independent function of gene and pseudogene mRNAs regulates tumour biology. Nature. 2010;465(7301):1033–8.
Bakhoum SF, Ngo B, Laughney AM, Cavallo JA, Murphy CJ, Ly P, Shah P, Sriram RK, Watkins TBK, Taunk NK, et al. Chromosomal instability drives metastasis through a cytosolic DNA response. Nature. 2018;553(7689):467–72.
Wu X, Lim ZF, Li Z, Gu L, Ma W, Zhou Q, Su H, Wang X, Yang X, Zhang Z. NORAD expression is associated with adverse prognosis in esophageal squamous cell carcinoma. Oncol Res Treat. 2017;40(6):370–4.
Li H, Wang X, Wen C, Huo Z, Wang W, Zhan Q, Cheng D, Chen H, Deng X, Peng C, et al. Long noncoding RNA NORAD, a novel competing endogenous RNA, enhances the hypoxia-induced epithelial-mesenchymal transition to promote metastasis in pancreatic cancer. Mol Cancer. 2017;16(1):169.
Zhang J, Li XY, Hu P, Ding YS. LncRNA NORAD contributes to colorectal cancer progression by inhibition of miR-202-5p. Oncol Res. 2018;26(9):1411–8.
Hu B, Cai H, Zheng R, Yang S, Zhou Z, Tu J. Long non-coding RNA 657 suppresses hepatocellular carcinoma cell growth by acting as a molecular sponge of miR-106a-5p to regulate PTEN expression. Int J Biochem Cell B. 2017;92:34–42.
Spiniello M, Knoener RA, Steinbrink MI, Yang B, Cesnik AJ, Buxton KE, Scalf M, Jarrard DF, Smith LM. HyPR-MS for multiplexed discovery of MALAT1, NEAT1, and NORAD lncRNA protein interactomes. J Proteome Res. 2018;17(9):3022–38.
Xu R, Rai A, Chen M, Suwakulsiri W, Greening DW, Simpson RJ. Extracellular vesicles in cancer - implications for future improvements in cancer care. Nat Rev Clin Oncol. 2018;15(10):617–38.
Fletcher CE, Deng L, Orafidiya F, Yuan W, Lorentzen M, Cyran OW, Varela-Carver A, Constantin TA, Leach DA, Dobbs FM, et al. A non-coding RNA balancing act: miR-346-induced DNA damage is limited by the long non-coding RNA NORAD in prostate cancer. Mol Cancer. 2022;21(1):82. https://doi.org/10.1186/s12943-022-01540-w.
Acknowledgements
Not applicable.
Funding
This work was supported by the National Key Research and Development Program of China (No. 2021YFA1301203), the National Natural Scientific Foundation of China (Nos. 82073945, 81874329, 81472802 and 82074000), the Hunan Provincial Natural Science Foundation of China (2020JJ4878), the Scientific Research Project of Department of Education of Hunan Province (20K136).
Author information
Authors and Affiliations
Contributions
RL, WZ and ML were responsible for study conception and design of the study. ML and WHH were responsible for drafting and revision of the manuscript. MYC, ML, QL LW, NY and JMH were responsible for data acquisition and analysis. ML and WZ were responsible for providing critical advice and interpretation of the results. HWH and WZ was responsible for funding support and final manuscript approval. All authors are read and approved by the final manuscript.
Corresponding authors
Ethics declarations
Ethical approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
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.
siRNA and shRNA sequences. Table S2. Primer sequences of primers used in PCR analysis
Additional file: 2 Fig. S1
Hub lncRNAs validation in TCGA-PRAD data. The expression levels of the 8 hub lncRNAs in adjacent normal and tumor tissues of PCa.
Additional file : 3 Fig. S2
Hub lncRNAs validation in tumor tissues from TCGA-PRAD data AThe expression levels of the 8 hub lncRNAs in I-IIIA stage and IIIB-IV stage group of PCa. B RFS curve of the 8 hub lncRNAs in PCa based on Kaplan–Meier analysis and log-rank test. Patients were divided into the high expression level and the low expression level based on the median value in I-IIIA stage or IIIB-IV stage group.
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.
About this article
Cite this article
Liu, M., Chen, MY., Huang, JM. et al. LncRNA weighted gene co-expression network analysis reveals novel biomarkers related to prostate cancer metastasis. BMC Med Genomics 15, 256 (2022). https://doi.org/10.1186/s12920-022-01410-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12920-022-01410-w