- Primary research
- Open Access
A novel prognostic cancer-related lncRNA signature in papillary renal cell carcinoma
Cancer Cell International volume 21, Article number: 545 (2021)
Papillary renal cell carcinoma (pRCC) ranks second in renal cell carcinoma and the prognosis of pRCC remains poor. Here, we aimed to screen and identify a novel prognostic cancer-related lncRNA signature in pRCC.
The RNA-seq profile and clinical feature of pRCC cases were downloaded from TCGA database. Significant cancer-related lncRNAs were obtained from the Immlnc database. Differentially expressed cancer-related lncRNAs (DECRLs) in pRCC were screened for further analysis. Cox regression report was implemented to identify prognostic cancer-related lncRNAs and establish a prognostic risk model, and ROC curve analysis was used to evaluate its precision. The correlation between RP11-63A11.1 and clinical characteristics was further analyzed. Finally, the expression level and role of RP11-63A11.1 were studied in vitro.
A total of 367 DECRLs were finally screened and 26 prognostic cancer-related lncRNAs were identified. Among them, ten lncRNAs (RP11-573D15.8, LINC01317, RNF144A-AS1, TFAP2A-AS1, LINC00702, GAS6-AS1, RP11-400K9.4, LUCAT1, RP11-63A11.1, and RP11-156L14.1) were independently associated with prognosis of pRCC. These ten lncRNAs were incorporated into a prognostic risk model. In accordance with the median value of the riskscore, pRCC cases were separated into high and low risk groups. Survival analysis indicated that there was a significant difference on overall survival (OS) rate between the two groups. The area under curve (AUC) in different years indicated that the model was of high efficiency in prognosis prediction. RP11-63A11.1 was mainly expressed in renal tissues and it correlated with the tumor stage, T, M, N classifications, OS, PFS, and DSS of pRCC patients. Consistent with the expression in pRCC tissue samples, RP11-63A11.1 was also down-regulated in pRCC cells. More importantly, up-regulation of RP11-63A11.1 attenuated cell survival and induced apoptosis.
Ten cancer-related lncRNAs were incorporated into a powerful model for prognosis evaluation. RP11-63A11.1 functioned as a cancer suppressor in pRCC and it might be a potential therapeutic target for treating pRCC.
Renal cell carcinoma (RCC) is a main human malignancy threatening people’s health globally. It is mainly divided into three catalogues, including clear cell RCC (ccRCC), papillary RCC (pRCC) and chromophobe RCC (chRCC), among which pRCC accounts of the second type with 15–20% of cases [1, 2]. PRCC is further divided into two subtypes, but the classification remains unsatisfactory . Neither radiotherapy nor chemotherapy is effective in pRCC. In addition, targeted therapy that widely used in ccRCC is not efficient in pRCC either [3, 4]. Therefore, it is urgent to explore promising biomarker and therapeutic target to improve pRCC patients’ prognosis.
In recent decades, long non-coding RNAs (lncRNAs) are reported to be aberrantly expressed in many malignancies and involved in the oncogenesis and progression of human cancers, and they have great potential to serve as biomarkers for cancer diagnosis and targets for cancer treatment . LncRNAs play roles in multiple cellular processes of cancer, including cell survival, migration, invasion, metastasis, and epithelial–mesenchymal transition (EMT), etc. [6, 7]. Additionally, lncRNAs are crucial regulators in many aspects, such as target mRNA expression regulation and cancer immunity regulation [8, 9]. What’s more, a series of lncRNAs have been identified as cancer-related lncRNAs [10,11,12,13]. For more cancer-related lncRNAs identification, different methods were performed and corresponding database was established [14, 15]. Therefore, lncRNAs, especially cancer-related lncRNAs are promising biomolecules in human cancers.
Immlnc database is a web-based resource to investigate the immune-related function of lncRNAs across cancer types . Moreover, it provides a larger number of cancer-related lncRNAs in seventeen cancer types. In this study, we aimed to combine the Immlnc and The Cancer Genome Altas (TCGA) databases to identify prognostic cancer-related lncRNAs and construct a risk model for prognosis prediction in pRCC. As expected, we finally identified 26 prognostic cancer-related lncRNAs and constructed a prognostic risk model with ten lncRNAs (including RP11-573D15.8, LINC01317, RNF144A-AS1, TFAP2A-AS1, LINC00702, GAS6-AS1, RP11-400K9.4, LUCAT1, RP11-63A11.1, RP11-156L14.1).
More interestingly, RP11-63A11.1 was a renal tissue-specific lncRNA among the ten lncRNAs. It was reported that RP11-63A11.1 was significantly down-regulated in pRCC samples and correlated with clinicopathological features of pRCC patients [17, 18], which was consistent with the analysis of ours. However, the role of RP11-63A11.1 is still not completely understood in pRCC. Thus, we further detected the expression of RP11-63A11.1 in pRCC cells and investigated its role on cell survival and apoptosis, which would add more convincing evidence to previous conclusion.
Acquisition and processing of samples data
The KIRP RNA-seq data (TCGA-KIRP) were acquired from UCSC, and clinical data of pRCC was downloaded from TCGA database. A total of 289 tumor samples and 32 matched non-tumor samples were included based on The Cancer Genome Atlas (TCGA) database. The samples with overall survival (OS) less than 1 month were eliminated, since these patients were likely to die of other reasons, such as serious infection [19, 20]. A total of 275 pRCC samples were screened out and divided into training set with 136 cases and testing set with 139 cases by the “caret” package of R software at random. The testing set and the combination set (the 275 samples) were used to validate the outcomes acquired from the training set.
Identification of differentially expressed mRNAs and cancer-related lncRNAs
The differentially expressed mRNAs (DEMs) between pRCC and its corresponding normal samples were analyzed by the “limma” package of R software. The significant cancer-related lncRNAs in pRCC were obtained from the Immlnc database. And then, their expression data in TCGA database was extracted and the “limma” package of R software was performed to analyze the differentially expressed cancer-related lncRNAs (DECRLs). |log2FC| > 2 and adjust P value < 0.01 were considered as the threshold. The DECRLs were removed if their expression values were equal to 0 in more than 135 samples. Finally, 367 DECRLs were screened out for further investigation.
Identification of prognostic cancer-related lncRNAs and establishment of a signature for prognosis prediction
We firstly implemented the univariate Cox regression analysis on the 367 DECRLs and clinical survival data in training set to determine prognostic cancer-related lncRNAs. And then we further analyzed the significant results using multivariate Cox regression analysis to generate a risk model. Each patient obtained a risk score with this formula: riskscore = β1 × exp(lnc1) + β2 × exp(lnc2)+···+βi × exp(lnci). Explnc is the expression value of DECRLs in the patients, while β refers to the regression co-efficient of cancer-related lncRNAs generated from the risk model in training set. Based on the medial riskscore, patients in all three sets were divided into high and low risk groups. Subsequently, the survival rate of the two groups was evaluated. The area under the ROC curve (AUC) was utilized to assess the precision of the risk model. The testing set and combination set were utilized to validate the above results.
Independence analysis of the cancer-related lncRNA signature in pRCC
Both the univariate and multivariate Cox regression were implemented to analyze the independence of this cancer-related lncRNA signature from the clinical features of pRCC patients, including age, gender, tumor stage and T classification. The data of M and N classifications were excluded in the independence analysis because the number of samples with undetermined classification was far beyond half of the all samples.
Weighted gene co-expression analysis (WGCNA) and gene enrichment analysis
In order to explore correlation between lncRNAs signature and biologic functions of pRCC, we performed WGCNA to construct the gene co-expression modules among differentially expressed mRNAs as previous study mentioned . The modules with the maximal absolute module significance correlated to lncRNAs signature were screened out. Then, gene enrichment analysis of genes in the most lncRNAs signature related module was carried out using DAVID . Gene Ontology (GO) terms with false discovery rate (FDR) less than 0.05 were considered as significant biologic functions.
Aberrant expressed mRNAs correlated with RP11-63A11.1
To explore the potential mechanism of RP11-63A11.1 involved in the progression of pRCC, we analyzed the aberrant expressed mRNAs which were correlated with RP11-63A11.1 expression. We listed aberrant expressed mRNAs in the combination set using Pearson method. The positively related mRNAs were screened out with correlation coefficient over 0.3.
Prediction of subcellular location
As the function of lncRNA is associated with its subcellular localization, the sequence of the RP11-63A11.1 was obtained from the LNCipedia website (https://lncipedia.org/) . The obtained sequence was input into the lncLocator website to predict the subcellular localization of the lncRNA . The results were obtained as a score for each potential subcellular location for RP11-63A11.1, including the cytoplasm, nucleus, ribosome, cytosol and exosome.
The association of RP11-63A11.1 expression and clinicopathological characteristics and prognostic indicators of pRCC
Among the ten lncRNAs, RP11-63A11.1 was primarily expressed in renal tissues and it was the most significant one with the lowest P value for OS in training set. Therefore, we further analyzed the association of RP11-63A11.1 expression and clinicopathological features (including gender, age, tumor stage and TNM classifications) as well as prognostic indicators of pRCC. Furthermore, we also investigated the expression and functions of RP11-63A11.1 in vitro.
Human pRCC SK-RC-39 cell line and normal renal cell line HK-2 were purchased from Cellcook Biotech (Guangzhou, China). The cell mediums (including RPM1-DMEM and RPM1-1640) and Fetal Bovine Serum (FBS) were purchased from Fcmacs Biotech (Nanjing, China). Antibiotic mixture (penicillin and streptomycin) and l-glutamine were purchased from Sangon Biotech (Shanghai, China). 57 mL FBS and 5.7 mL antibiotic mixture as well as 5.7 mL l-glutamine were added into cell mediums to generate complete mediums. SK-RC-39 and HK-2 were fed in complete 1640 and DMEM medium, respectively. Both the two kinds of cell lines were cultured in the 5% CO2 incubator at 37 °C.
Lenti-sgRNA RP11-63A11.1 was constructed by lenti-sgRNA. We also generate the SK-RC-39dCas9+MPH+ cells before they were infected with lenti-sgRNA RP11-63A11.1. We then used the puromycin to screen the cells and picked up clones. The over-expressions of RP11-63A11.1 were subsequently validated in clones. The clones with elevated RP11-63A11.1 were used for further assays.
Following the manufacturer’s protocol, EASYspin Tissue/Cell RNA Rapid Extraction Kit (Chaoyan Biotech, Shanghai, China) were used to extract the total RNA of cells. Then reverse transcription kit (Vazyme, Nanjing, China) was used to synthesize cDNA. Based on the instruction of ChamQ SYBR qPCR Master Mix Kit (Vazyme), qPCR amplification was performed on QuantStudio 5. GAPDH was considered as the internal reference. The expression level of RP11-63A11.1 was analyzed by 2−ΔΔCt method. Primers: RP11-63A11.1, forward, 5′-TCAGCAGGGTTTAGAGCAGC-3′; reverse, 3′-CTGAGGTTTCCATGCTGCTG-5′; GAPDH, forward: 5′-GAGTCAACGGATTTGGTCGT-3′, reverse: 5′-TTGATTTTGGAGGGATCTCG-3′.
Cell proliferation assay
CCK-8 (Vazyme, Nanjing, China) assay was applied for the evaluation of cell proliferation ability. In brief, after counting with the cell count plate, 1 × 104 cells were seeded into each experimental well of a 96-well plate and then cultured at 37 °C for 24 h, 48 h and 72 h. After adding 10% CCK-8 reagent into the corresponding wells, the absorbance value at 450 nm was measured.
Cell apoptosis assay
For cell apoptosis detection, we prepared and adjusted the cell concentration to 1 × 106/mL. Then Annexin V-Alexa Fluor 488/PI Apoptosis detection kit (Fcmacs Biotech, Nanjing, China) was used for cell staining and flow cytometer was used for the assessment of cell apoptosis.
The total protein was extracted using radioimmunoprecipitation assay solution (Beyotime, Shanghai, China) before polypropylene gel electrophoresis. Then the protein was further transferred to polyvinylidene difluoride membranes (Millipore, USA). After blocked by 3% bovine serum albumin (Fcmacs Biotech, Nanjing, China) in room termperature, the membrane was incubated with primary antibody overnight. Sources of primary antibodies (cleaved PARP, cleaved Caspase 3, Bcl-2, Bax, Cytochrome C and GAPDH) were from Cell Signaling (Danvers, MA, USA). RASSF10 antibody was from Abcam (Cambridge, England, UK). After the membrane was incubated with secondary antibody, the proteins were visualized by the gel-imaging system.
R software 4.0.2 and SPSS 25 were used to implement statistical analysis in this study. Univariate and multivariate Cox regression analysis were applied to generate a cancer-related lncRNA signature and identify independent prognostic factors for pRCC. The survival rate of pRCC patients in two risk-groups was appraised by Kaplan–Meier method. The P value < 0.05 was generally considered statistically significant.
Differentially expressed mRNAs and cancer-related lncRNAs in pRCC
A total of 1872 cancer-related lncRNAs with expression data in pRCC were extracted (Additional file 1: Table S1). Through the “limma” package of R software, 2030 DEMs were identified, including 219 up-regulated DEMs and 1811 down-regulated DEMs (Fig. 1A, B), while 579 differentially expressed cancer-related lncRNAs (DECRLs) were found, including 84 up-regulated DECRLs and 495 down-regulated DECRLs (Fig. 1C, D). The DECRLs were removed if their expression values were equal to 0 in more than 135 samples. 367 DECRLs were finally determined for further analysis.
Identification and validation of cancer-related lncRNA signature
The expression profile of 367 cancer-related lncRNAs with OS data was analyzed by the univariate Cox regression to identify the prognostic lncRNAs in pRCC in the training set first. There were 26 lncRNAs identified as prognostic cancer-related lncRNAs in total (Fig. 2). Subsequently, multivariate Cox regression was conducted to generate a risk model of prognosis. 16 genes were screened off by the regression analysis. Then, the remaining 10 lncRNAs (RP11-573D15.8, LINC01317, RNF144A-AS1, TFAP2A-AS1, LINC00702, GAS6-AS1, RP11-400K9.4, LUCAT1, RP11-63A11.1, and RP11-156L14.1) were included in the model (Table 1). Riskscore = 0.7953 × exp(RP11-573D15.8) + 0.4074 × exp(LINC01317) + 0.5932 × exp(RNF144A-AS1) + 0.2636 × exp(TFAP2A-AS1) + 0.5327 × exp(LINC00702) − 0.6131 × exp(GAS6-AS1) − 0.8397 × exp(RP11-400K9.4) + 0.3506 × exp(LUCAT1) − 0.6670 × exp(RP11-63A11.1) − 0.6442 × exp(RP11-156L14.1). According to the riskscore formula, the riskscore of each pRCC patient in all three sets was calculated. On the basis of median riskscore value, pRCC cases were segmented into high and low risk groups. As shown in Fig. 3A, the mortality rate of pRCC patients was progressively increased with the rising of risk score in training-set, testing-set and combination-set. The levels of the ten lncRNAs in two groups were shown in Fig. 3B. The survival curves displayed that patients in high risk group had a shorter survival time than those in low risk group (Fig. 3C). The precision of the 10 cancer-related lncRNA signature in predicting prognosis of pRCC were analyzed through the ROC curves. The results demonstrated that all of the AUC values at 1, 3, and 5 years were over than 0.75 in the training-set, testing-set and the combination-set (Fig. 3D), indicating the ten cancer-related lncRNAs signature had a good sensitivity and specificity in prognosis prediction.
The cancer-related lncRNA signature was an independent factor in prognosis prediction
We further evaluated whether the identified cancer-related lncRNA signature was an independent factor for prognosis prediction in pRCC. Co-variables including age, gender, tumor stage and the risk model in training-set, testing-set and combination-set were analyzed by Cox regression analysis. The outcomes from univariate Cox regression revealed that tumor stage and the risk model were significantly relevant to OS in all three sets. Multivariate Cox regression analysis suggested that tumor stage and risk model were independent factors in the combination set. The results were listed in Table 2. Taken together, the cancer-related lncRNA signature was an independent factor for prognosis prediction in pRCC.
By performing stratification analysis on the clinicopathological features of pRCC patients in the combination cohort, we revealed that the cancer-related lncRNA signature was a good predictor. As shown in Fig. 4A, B, both older (> 60 years) and younger (≤ 60 years) pRCC cases in high risk group exhibited a worse prognosis compared with those cases in low risk group. Of note, similar results of the cancer-related lncRNA signature were found in different genders (Fig. 4C, D) and tumor stages (Fig. 4E, F).
The cancer-related lncRNA signature was associated with biological processes
The cancer-related lncRNA signature was vital in prognosis of pRCC patients; therefore this signature should be closely associated with biological processes of pRCC. 2030 DEMs were used to construct 8 similar expression modules by average linkage clustering. The relevance with P value between each module and cancer-related lncRNA signature was listed in every module (Fig. 5A). The most negative related module (yellow module, R = − 0.13, P = 0.04) and positive related module (blue module, R = 0.27, P = 5e−06) were selected out for gene enrichment analysis. The FDR of GO term of genes in yellow module was over than 0.05. The genes in blue module were mainly involved in cell division, cell proliferation, and cell cycle process, such as G1/S transition of mitotic cell cycle, sister chromatid cohesion (Fig. 5B), etc.
RP11-63A11.1 was associated with many clinical characteristics and prognostic indicators of pRCC
Among all the 10 lncRNAs in the signature, RP11-63A11.1 was the most significant one with the lowest P value (P = 0.0014) by Kaplan-Meier method in the training set (Additional file 2: Fig. S1). In order to further understand the role of RP11-63A11.1 in pRCC, we analyzed the association of RP11-63A11.1 with the clinical features including gender, age, clinical stages, TNM classifications, OS, disease-free survival (DFS), disease-specific survival (DSS), progression-free survival (PFS). In the GEPIA database, we found that RP11-63A11.1 was mainly expressed in kidney, suggesting RP11-63A11.1 may be a renal tissue-specific lncRNA (Fig. 6A). It was indicated that the level of RP11-63A11.1 was significantly higher in lower tumor stages and TNM classifications, while it did not show significant difference in age and gender (Fig. 6B–G). Of note, higher level of RP11-63A11.1 indicates favorable OS (Fig. 7A), PFS (Fig. 7B), as well as DSS (Fig. 7C). Although there is no significant difference of DFS between high expression group and low expression group (Fig. 7D), RP11-63A11.1 could be a very critical prognosis predictor based on the data of OS, DSS, and PFS.
Elevated RP11-63A11.1 inhibited survival and induced apoptosis of pRCC cells
In order to validate the role of RP11-63A11.1, we determined cell survival and apoptosis in pRCC cell line (SK-RC-39). We firstly determined the level of RP11-63A11.1 in pRCC cells. The result showed that the expression of RP11-63A11.1 was decreased in pRCC cells compared with normal renal cells (HK-2) (Fig. 8A). Then, we generated 2 clones with CRISPR/dCas9 system to elevate the level of RP11-63A11.1 in both SK-RC-39 and A498 (Additional file 3: Fig. S2). Subsequently, the cell survival assay indicated that elevated RP11-63A11.1 didn’t inhibit the HK-2 (Fig. 8B); however, it significantly inhibited the cell survival of SK-RC-39 and A498 (Fig. 8C, D). In consistent with that, the elevated RP11-63A11.1 didn’t affect the HK-2 (Fig. 8E) but induced the cell apoptosis of SK-RC-39 (Fig. 8F) and A498 (Fig. 8G), suggesting that RP11-63A11.1 could be a suppressive gene in pRCC. Over-expressed RP11-63A11.1 might be an effective way in therapy. Moreover, five apoptotic markers were investigated. It was shown in Fig. 9 that over-expression of RP11-63A11.1 increased level of cleaved-PARP while uncleaved-PARP significantly decreased (Fig. 9A). Similarly, we found the introduction of RP11-63A11.1 elevated cleaved-Caspase-3 (Fig. 9B), adding more evidence to the apoptosis induced by RP11-63A11.1. Moreover, over-expression of RP11-63A11.1 suppressed the expression of Bcl-2 (Fig. 9C), which was an inhibitor of apoptosis. On the other hand, Bax (Fig. 9D) and Cytochrome C (Fig. 9E) both increased upon introduction of RP11-63A11.1, which was consistent with those data in cell apoptosis assay.
Analysis of the potential mechanism of RP11-63A11.1 in pRCC
How the RP11-63A11.1 regulates the apoptosis progression in pRCC cells is still unknown. To further investigate the possible mechanism, we first predict the subcellular location. The sequence information of RP11-63A11.1 was obtained from LNCipedia (Additional file 4: Table. S2). It is indicated in the LncLocator that RP11-63A11.1 is probably located in the cytoplasm (Table 3), suggesting that RP11-63A11.1 may serve as a ceRNA (competing endogenous RNA). Then, we showed all the genes which were positively correlated with RP11-63A11.1 (Additional file 5: Table S3). Among all the positively correlated genes, we found that RASSF10 was attractive. Finally, we analyzed the expression of RASSF10. Elevated RP11-63A11.1 significantly increased the expression of RASSF10 (Additional file 6: Fig.S3). Thus, we believe that introduction of RP11-63A11.1, as a ceRNA, elevated RASSF10, inducing apoptosis in pRCC cells.
Similar to ccRCC, pRCC patients are usually diagnosed at an advanced stage or metastasis. The major treatment for local pRCC is surgical operation. However, approximately 40% patients reoccur after surgical resection, resulting in a worse prognosis . In recent years, increasing researches of cancers focused on molecular characteristic in early diagnosis and prognosis improvement. Of note, lncRNAs have been revealed to function as main regulators in various biological processes [26, 27]. Notably, it was believed that cancer-related lncRNAs may be helpful for researchers to understand the mechanism of cancer progression and develop efficient measures to improve the prognosis . And cancer-related lncRNAs have been described in some cancer types, such as esophageal adenocarcinoma  and prostate cancer . In the present study, we obtained 1872 caner-related lncRNAs through integrating Immlnc and TCGA databases and screened 367 significant DECRLs in pRCC. Cox regression analysis suggested that 26 significant DECRLs were related to the survival of patients with pRCC, among which 10 lncRNAs (RP11-573D15.8, LINC01317, RNF144A-AS1, TFAP2A-AS1, LINC00702, GAS6-AS1, RP11-400K9.4, LUCAT1, RP11-63A11.1, RP11-156L14.1) were enrolled into the risk model to generate a cancer-related lncRNA signature in pRCC for survival prediction. More surprisingly, the signature was a prognostic risk factor from other clinical features in pRCC.
Based on this signature, survival curves suggested that patients in high risk group had a poorer survival rate than patients in low risk group. Additionally, ROC curves exhibited that the precision of this signature was more than 75% for survival prediction at different years. These results were validated in testing set and combination set and similar outcomes were observed, suggesting the signature was a forceful tool for predicting prognosis. What’s more, this cancer-related lncRNA signature was also a forceful tool for prognosis prediction in different classes of some clinicopathological features, including age, gender, and tumor stage.
Other prognostic signatures in pRCC were previously reported. For instance, Wang et al.  established a prognostic signature with 15 immune-related genes. Although prognosis of patients could be predicted by their signature, the accuracy of their prognostic signature in predicting prognosis in 1, 3 and 5 years was lower than that of ours (Table 4). In addition, Gao et al.  also grouped pRCC patients into training set and testing set. A five mRNAs signature was identified in training set and the accuracy of their signature was 0.82, which was lower than the accuracy in our training set. Although the result of testing set showed that patients could be divided into high and low risk groups by the five mRNAs signature, the accuracy of their signature in testing set was not assessed. Zhang et al.  constructed a risk model by 17 mutant genes. The AUC of their signature was 0.907 in 3 years, which was a considerable result. However, the result was not validated in their study and it was unknown whether the mutant-gene signature was efficient for predicting prognosis in 5 years. Similarly, a signature of four lncRNAs was constructed to predict prognosis . In spite of significant accuracy of this signature displayed, a validation result was not provided in their study. Taken together, our cancer-related lncRNA signature might be more beneficial than the previous signatures mainly because of higher accuracy and validation of the findings.
Among the lncRNAs in our cancer-related signature, RNF144A-AS1, now called GRASLND, was reported to serve as an important regulator in stem cell chondrogenesis . RNF144A-AS1 could also facilitate the migration and invasion of bladder cancer cells . LINC00702 was a newly identified lncRNA and was involved in the progression of several malignancies via tumorigenesis-associated pathways such as Wnt/β-catenin pathway and PTEN/PI3K-AKT pathway [37,38,39,40]. LUCAT1 was the only lncRNA widely investigated in human cancers among the ten lncRNAs. LUCAT1 could promote tumorigenesis and development in various cancers. Besides, it could promote anti-tumor drug resistance in some tumor types such as NSCLC and osteosarcoma [41, 42]. These results suggest that these lncRNAs are reasonable and of importance in human cancers. However, to our knowledge, they were firstly uncovered to be novel prognostic biomarkers in pRCC in this study. In consistent with previous researches [17, 18] in which RP11-63A11.1 was reported, we found that RP11-63A11.1 correlated with clinicopathological features and prognosis of pRCC patients. However, its role in pRCC cells was incompletely understood. We revealed that RP11-63A11.1 was decreased in pRCC cells. Furthermore, increased RP11-63A11.1 inhibited the proliferation and induced apoptosis of pRCC cells, indicating RP11-63A11.1 served as a tumor suppressor in pRCC. Subcellular localization of lncRNA transcripts is of great importance when we investigate its mechanism. Therefore, the subcellular location of the RP11-63A11.1 transcript was predicted using lncLocator based on its sequence acquired from LNCipedia. RP11-63A11.1 was likely to be located in the cytoplasm. LncRNA located in the cytoplasm or cytosol probably serves as a ceRNA, and regulates mRNA stability or translation. Regarding that elevated RP11-63A11.1 significantly increased the expression of RASSF10, we hypothesized that RP11-63A11.1 acted as ceRNA, and regulated the stability or translation of RASSF10. Further studies are indicated to investigate more specific mechanisms of RP11-63A11.1 in pRCC.
The roles of other lncRNAs including RP11-573D15.8, LINC01317, RP11-400K9.4, and RP11-156L14.1 were not reported in previous literatures. These lncRNAs were firstly reported in the present study. When combined with the above-mentioned lncRNAs, these lncRNAs generate a powerful tool for survival prediction in pRCC. To explore the biological functions of the cancer-related lncRNA signature in pRCC, we performed WGCNA to seek gene modules associated with this signature. Two gene modules were finally identified. Gene enrichment analysis showed genes in the blue module (positively correlated) were mainly enriched in biologic processes of cell division, proliferation and cell cycle, indicating that pRCC with high risk score has stronger proliferation capability than pRCC with low risk score.
In summary, we successfully figured out a novel cancer-related lncRNA signature with powerful predictive function for pRCC prognosis. These results of this study may facilitate understanding of molecular mechanisms in pRCC development, and the cancer-related lncRNA signature may become a promising biomarker and therapeutic target in pRCC. Of note, RP11-63A11.1 may serve as a ceRNA and regulate RASSF10, facilitating the understanding of the development in pRCC.
Availability of data and materials
All data generated or analyzed during this study are included in Additional files.
Papillary renal cell carcinoma
Long non-coding RNA
The Cancer Genome Atlas
Differentially expressed cancer-related lncRNAs
Receiver operating characteristic
Area under curve
Competing endogenous RNA
Differentially expressed mRNAs
Moch H, Cubilla AL, Humphrey PA, Reuter VE, Ulbright TM. The 2016 WHO classification of tumours of the urinary system and male genital organs—Part A: renal, penile, and testicular tumours. Eur Urol. 2016;70(1):93–105.
Hsieh JJ, Purdue MP, Signoretti S, Swanton C, Albiges L, Schmidinger M, et al. Renal cell carcinoma. Nat Rev Dis Prim. 2017;3:17009.
Akhtar M, Al-Bozom IA, Al Hussain T. Papillary renal cell carcinoma (PRCC): an update. Adv Anat Pathol. 2019;26(2):124–32.
Krawczyk KM, Nilsson H, Allaoui R, Lindgren D, Arvidsson M, Leandersson K, et al. Papillary renal cell carcinoma-derived chemerin, IL-8, and CXCL16 promote monocyte recruitment and differentiation into foam-cell macrophages. Lab Invest. 2017;97(11):1296–305.
Zhang Y, Tang L. The application of lncRNAs in cancer treatment and diagnosis. Recent Pat Anti-Cancer Drug Discov. 2018;13(3):292–301.
Bach DH, Lee SK. Long noncoding RNAs in cancer cells. Cancer Lett. 2018;419:152–66.
Bolha L, Ravnik-Glavač M, Glavač D. Long noncoding RNAs as biomarkers in cancer. Dis Markers. 2017;2017:7243968.
Rao A, Rajkumar T, Mani S. Perspectives of long non-coding RNAs in cancer. Mol Biol Rep. 2017;44(2):203–18.
Denaro N, Merlano MC, Lo Nigro C. Long noncoding RNAs as regulators of cancer immunity. Mol Oncol. 2019;13(1):61–73.
Li J, Li Z, Leng K, Xu Y, Ji D, Huang L, et al. ZEB1-AS1: a crucial cancer-related long non-coding RNA. Cell Prolif. 2018;51(1):e12423.
Thin KZ, Liu X, Feng X, Raveendran S, Tu JC. LncRNA-DANCR: a valuable cancer related long non-coding RNA for human cancers. Pathol Res Pract. 2018;214(6):801–5.
Yu X, Zheng H, Chan MT, Wu WKK. BANCR: a cancer-related long non-coding RNA. Am J Cancer Res. 2017;7(9):1779–87.
Li J, Li Z, Zheng W, Li X, Wang Z, Cui Y, et al. PANDAR: a pivotal cancer-related long non-coding RNA in human cancers. Mol bioSyst. 2017;13(11):2195–201.
Zhang X, Li T, Wang J, Li J, Chen L, Liu C. Identification of cancer-related long non-coding RNAs using XGBoost with high accuracy. Front Genet. 2019;10:735.
Wang J, Zhang X, Chen W, Li J, Liu C. CRlncRNA: a manually curated database of cancer-related long non-coding RNAs with experimental proof of functions on clinicopathological and molecular features. BMC Med Genom. 2018;11(Suppl 6):114.
Li Y, Jiang T, Zhou W, Li J, Li X, Wang Q, et al. Pan-cancer characterization of immune-related lncRNAs identifies potential oncogenic biomarkers. Nat Commun. 2020;11(1):1000.
Lan H, Zeng J, Chen G, Huang H. Survival prediction of kidney renal papillary cell carcinoma by comprehensive LncRNA characterization. Oncotarget. 2017;8(67):110811–29.
Yang F, Song Y, Ge L, Zhao G, Liu C, Ma L. Long non-coding RNAs as prognostic biomarkers in papillary renal cell carcinoma. Oncol Lett. 2019;18(4):3691–7.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102(43):15545–50.
Wei C, Liang Q, Li X, Li H, Liu Y, Huang X, et al. Bioinformatics profiling utilized a nine immune-related long noncoding RNA signature as a prognostic target for pancreatic cancer. J Cell Biochem. 2019;120(9):14916–27.
Wen Jiang Q, Guo C, Wang Yu, Zhu. A nomogram based on 9-lncRNAs signature for improving prognostic prediction of clear cell renal cell carcinoma. Cancer Cell Int. 2019;19:208.
Huang DW, Brad T, Sherman, Richard A, Lempicki. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Volders P-J, Anckaert J, Verheggen K, Nuytens J, Martens L, Mestdagh P, et al. LNCipedia 5: towards a reference set of human long non-coding RNAs. Nucleic Acids Res. 2019;47(D1):D135–9.
Zhen Cao X, Pan Y, Yang Y, Huang H-B, Shen. The lncLocator: a subcellular localization predictor for long non-coding RNAs based on a stacked ensemble classifier. Bioinformatics. 2018;34(13):2185–94.
Courthod G, Tucci M, Di Maio M, Scagliotti GV. Papillary renal cell carcinoma: a review of the current therapeutic landscape. Crit Rev Oncol Hematol. 2015;96(1):100–12.
Peng WX, Koirala P, Mo YY. LncRNA-mediated regulation of cell signaling in cancer. Oncogene. 2017;36(41):5661–7.
Flippot R, Beinse G, Boilève A, Vibert J, Malouf GG. Long non-coding RNAs in genitourinary malignancies: a whole new world. Nat Rev Urol. 2019;16(8):484–504.
Liu Z, Zhang Y, Han X, Li C, Yang X, Gao J, et al. Identifying cancer-related lncRNAs based on a convolutional neural network. Front Cell Dev Biol. 2020;8:637.
Yu Y, Chen X, Cang S. Cancer-related long noncoding RNAs show aberrant expression profiles and competing endogenous RNA potential in esophageal adenocarcinoma. Oncol Lett. 2019;18(5):4798–808.
Mitobe Y, Takayama KI, Horie-Inoue K, Inoue S. Prostate cancer-associated lncRNAs. Cancer Lett. 2018;418:159–66.
Wang Z, Song Q, Yang Z, Chen J, Shang J, Ju W. Construction of immune-related risk signature for renal papillary cell carcinoma. Cancer Med. 2019;8(1):289–304.
Gao Z, Zhang D, Duan Y, Yan L, Fan Y, Fang Z, et al. A five-gene signature predicts overall survival of patients with papillary renal cell carcinoma. PLoS ONE. 2019;14(3):e0211491.
Zhang C, Zheng Y, Li X, Hu X, Qi F, Luo J. Genome-wide mutation profiling and related risk signature for prognosis of papillary renal cell carcinoma. Ann Transl Med. 2019;7(18):427.
Liu Y, Gou X, Wei Z, Yu H, Zhou X, Li X. Bioinformatics profiling integrating a four immune-related long non-coding RNAs signature as a prognostic model for papillary renal cell carcinoma. Aging. 2020;12(15):15359–73.
Huynh NP, Gloss CC, Lorentz J, Tang R, Brunger JM, McAlinden A, et al. Long non-coding RNA GRASLND enhances chondrogenesis via suppression of the interferon type II signaling pathway. eLife. 2020;9:e49558.
Wang Y, Du L, Yang X, Li J, Li P, Zhao Y, et al. A nomogram combining long non-coding RNA expression profiles and clinical factors predicts survival in patients with bladder cancer. Aging. 2020;12(3):2857–79.
Li T, Ren J, Ma J, Wu J, Zhang R, Yuan H, et al. LINC00702/miR-4652-3p/ZEB1 axis promotes the progression of malignant meningioma through activating Wnt/β-catenin pathway. Biomed Pharmacother = Biomedecine & pharmacotherapi. 2019;113:108718.
Yu W, Li D, Ding X, Sun Y, Liu Y, Cong J, et al. LINC00702 suppresses proliferation and invasion in non-small cell lung cancer through regulating miR-510/PTEN axis. Aging. 2019;11(5):1471–85.
Wang L, Ye TY, Wu H, Chen SY, Weng JR, Xi XW. LINC00702 accelerates the progression of ovarian cancer through interacting with EZH2 to inhibit the transcription of KLF2. Eur Rev Med Pharmacol Sci. 2019;23(3 Suppl):201–8.
Yu D, Wang XY, Jin ZL. Linc00702 inhibits cell growth and metastasis through regulating PTEN in colorectal cancer. Eur Rev Med Pharmacol Sci. 2020;24(7):3624–32.
Wang W, Dong ML, Zhang W, Liu T. Long noncoding LUCAT1 promotes cisplatin resistance of non-small cell lung cancer by promoting IGF-2. Eur Rev Med Pharmacol Sci. 2019;23(12):5229–34.
Han Z, Shi L. Long non-coding RNA LUCAT1 modulates methotrexate resistance in osteosarcoma via miR-200c/ABCB1 axis. Biochem Biophys Res Commun. 2018;495(1):947–53.
This work was supported by grants No. 81402100 from National Natural Science Foundation of China, No. BK20211123 from Natural Science Foundation of Jiangsu Province, No. SH2021033 from the Social Development Foundation of Zhenjiang, and No. 2019JJ40269 from Natural Science Foundation of Hunan Province.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The cancer-related lncRNAs with expression data in pRCC.
Risk modes of every LncRNA in the signature. Among all the 10 lncRNAs in the signature, RP11-63A11.1 was the most significant one with the lowest P value (P = 0.0014).
RP11-63A11.1 was introduced into SK-RC-39 and A498 cells. After introduced into SK-RC-39 and A498 cells, two clones in each cell line were generated and the level of RP11-63A11.1 was significantly increased.
Sequence of RP11-63A11.1 obtained from LNCipedia.
The genes which are positively correlated with RP11-63A11.1
Elevated RP11-63A11.1 significantly increased the expression of RASSF10.
About this article
Cite this article
Chen, B., Dong, D., Yao, Q. et al. A novel prognostic cancer-related lncRNA signature in papillary renal cell carcinoma. Cancer Cell Int 21, 545 (2021). https://doi.org/10.1186/s12935-021-02247-6
- Papillary renal cell carcinoma