Development and verification of a nomogram for prediction of recurrence‐free survival in clear cell renal cell carcinoma

Nowadays, gene expression profiling has been widely used in screening out prognostic biomarkers in numerous kinds of carcinoma. Our studies attempt to construct a clinical nomogram which combines risk gene signature and clinical features for individual recurrent risk assessment and offer personalized managements for clear cell renal cell carcinoma. A total of 580 differentially expressed genes (DEGs) were identified via microarray. Functional analysis revealed that DEGs are of fundamental importance in ccRCC progression and metastasis. In our study, 338 ccRCC patients were retrospectively analysed and a risk gene signature which composed of 5 genes was obtained from a LASSO Cox regression model. Further analysis revealed that identified risk gene signature could usefully distinguish the patients with poor prognosis in training cohort (hazard ratio [HR] = 3.554, 95% confidence interval [CI] 2.261‐7.472, P < .0001, n = 107). Moreover, the prognostic value of this gene‐signature was independent of clinical features (P = .002). The efficacy of risk gene signature was verified in both internal and external cohorts. The area under receiver operating characteristic curve of this signature was 0.770, 0.765 and 0.774 in the training, testing and external validation cohorts, respectively. Finally, a nomogram was developed for clinicians and did well in the calibration plots. This nomogram based on risk gene signature and clinical features might provide a practical way for recurrence prediction and facilitating personalized managements of ccRCC patients after surgery.     

VNN3 is a potential novel biomarker for predicting prognosis in clear cell renal cell carcinoma

Although pathological observations provide approximate prognoses, it is difficult to achieve prognosis in patients with existing prognostic factors. Therefore, it is very important to find appropriate biomarkers to achieve accurate cancer prognosis. Renal cell carcinoma (RCC) has several subtypes, the discrimination of which is crucial for proper treatment. Here, we present a novel biomarker, VNN3, which is used to prognose clear cell renal cell carcinoma (ccRCC), the most common and aggressive subtype of kidney cancer. Patient information analyzed in our study was extracted from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) cohorts. VNN3 expression was considerably higher in stages III and IV than in stages I and II. Moreover, Kaplan–Meier curves associated high VNN3 expression with poor prognoses (TCGA, p?p?=?.00076), confirming that ccRCC prognosis can be predicted via VNN3 expression patterns. Consistent with all patient results, the prognosis of patients with higher VNN3 expression was worse in both low stage (I and II) and high stage (III and IV) (TCGA, p < 0.0001 in stage I and II; ICGC, p = 0.028 in stage I and II; TCGA, p = 0.005 in stage III and IV). Area under the curve and receiver operating characteristic curves supported our results that highlighted VNN3 expression as a suitable ccRCC biomarker. Multivariate analysis also verified the prognostic performance of VNN3 expression (TCGA, p?p?=?.017). Altogether, we suggest that VNN3 is applicable as a new biomarker to establish prognosis in patients with ccRCC.     

Genome-wide screening and cohorts validation identifying novel lncRNAs as prognostic biomarkers for clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the main subtype of renal cell carcinoma with varied prognosis. We aimed to identify and assess the possible prognostic long noncoding RNA (lncRNA) biomarkers. LncRNAs expression data and corresponding clinical information of 619 ccRCC patients were downloaded from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases. Differentially expressed genes analysis, univariate Cox regression, the least absolute shrinkage and selection operator Cox regression model were utilized to identify hub lncRNAs. Multivariate Cox regression was used to establish the risk model. Statistical analysis was performed using R 3.5.3. The expression value of five lncRNAs and the risk-score levels were significantly associated with a survival prognosis of ccRCC patients (all P < .001). In the TCGA validation cohort, the area under the curve (AUC) for the integrated nomogram was 0.905 and 0.91 for 3-, 5-year prediction separately. The AUC reached up to 0.757 in an independent ICGC cohort. Besides, the calibration plots also illustrated well curve-fitting between observation values and predictive values. Weighted gene co-expression network analysis and subsequent pathway analysis revealed that the PI3K-Akt-mTOR and hypoxia-inducible factor signaling crosstalk might function as the most essential mechanisms related to the five-lncRNAs signature. Our study suggested that lncRNA AC009654.1, AC092490.2, LINC00524, LINC01234, and LINC01885 were significantly associated with ccRCC prognosis. The prognostic model based on this five lncRNA may predict the overall survival of ccRCC.     

Prognostic value of a gene signature in clear cell renal cell carcinoma

Renal cancer is a common urogenital system malignance. Novel biomarkers could provide more and more critical information on tumor features and patients’ prognosis. Here, we performed an integrated analysis on the discovery set and established a three-gene signature to predict the prognosis for clear cell renal cell carcinoma (ccRCC). By constructing a LASSO Cox regression model, a 3-messenger RNA (3-mRNA) signature was identified. Based on the 3-mRNA signature, we divided patients into high- and low-risk groups, and validated this by using three other data sets. In the discovery set, this signature could successfully distinguish between the high- and low-risk patients (hazard ratio (HR), 2.152; 95% confidence interval (CI),1.509–3.069; p < 0.0001). Analysis of internal and two external validation sets yielded consistent results (internal: HR, 2.824; 95% CI, 1.601–4.98; p < 0.001; GSE29609: HR, 3.002; 95% CI, 1.113–8.094; p = 0.031; E-MTAB-3267: HR, 2.357; 95% CI, 1.243–4.468; p = 0.006). Time-dependent receiver operating characteristic (ROC) analysis indicated that the area under the ROC curve at 5 years was 0.66 both in the discovery and internal validation set, while the two external validation sets also suggested good performance of the 3-mRNA signature. Besides that, a nomogram was built and the calibration plots and decision curve analysis indicated the good performance and clinical utility of the nomogram. In conclusion, this 3-mRNA classifier proved to be a useful tool for prognostic evaluation and could facilitate personalized management of ccRCC patients.     

Serum exosomal miR-210 as a potential biomarker for clear cell renal cell carcinoma

Exosomal microRNAs (miRNAs) are suggested to reflect molecular changes occurring in their cells of origin and are potential indicators in the early detection of cancers. This study aimed to determine whether certain exosomal miRNAs from tumor tissue can be used as noninvasive biomarkers for clear cell renal cell carcinoma (ccRCC). Based on ccRCC miRNA expression profiles and the literature, we selected six miRNAs (miR-210, miR-224, miR-452, miR-155, miR-21, and miR-34a) and analyzed their expression in tissues, sera, and serum exosomes through quantitative real-time polymerase chain reaction in hypoxia-induced (with CoCl₂) renal cell lines. miR-210, miR-224, miR-452, miR-155, and miR-21 were upregulated in tumor tissues compared with normal tissues. Serum miR-210 and miR-155 levels were higher in patients with ccRCC than in healthy controls (HCs). Furthermore, only exosomal miR-210 was significantly upregulated in patients with ccRCC than in HCs. Moreover, receiver operating characteristic (ROC) curve analysis revealed an area under the ROC curve of 0.8779 (95% confidence interval, 0.7987-0.9571) and a sensitivity and specificity of 82.5% and 80.0%, respectively. Moreover, exosomal miR-210 was upregulated at an advanced stage, and Fuhrman grade and metastasis decreased significantly one month after surgery. Acute hypoxia exposure activates miR-210 and release of exosomes with upregulated miR-210 in both normal and tumor RCC cell lines and interferes with vacuole membrane protein 1 mRNA expression, especially in the metastatic ccRCC cell line. In conclusion, Serum exosomal miR-210 originating from tumor tissue has potential as a novel noninvasive biomarker for the detection and prognosis of ccRCC.     

Identification of a three-miRNA signature as a novel potential prognostic biomarker in patients with clear cell renal cell carcinoma

Current studies suggest that some microRNAs (miRNAs) are associated with prognosis in clear cell renal cell carcinoma (ccRCC). In this paper, we aimed to identify a miRNAs signature to improve prognostic prediction for ccRCC patients. Using ccRCC RNA-Seq data of The Cancer Genome Atlas (TCGA) database, we identified 177 differentially expressed miRNAs between ccRCC and paracancerous tissue. Then all the ccRCC tumor samples were divided into training set and validation set randomly. Three-miRNA signature including miR130b, miR-18a, and miR-223 were constructed by the least absolute shrinkage and selection operator (LASSO) Cox regression model in training set. According to optimal cut-off value of three-miRNA signature risk score, all the patients could be classified into high-risk group and low-risk group significantly. Survival of patients was significantly different between two groups (hazard ratio, 5.58, 95% confidence interval, 3.17-9.80; P < 0.0001), and three-miRNA signature performed favorably prognostic and predictive accuracy. The results were further validated in the validation set and total set. Multivariate Cox regression analyses and subgroup analyses showed that three-miRNA signature was an independent prognostic factor. Two nomograms that integrated three-miRNA signature and three clinicopathological risk factors were constructed to predict overall survival and disease-free survival after surgery for ccRCC patients. Functional enrichment analysis showed the possible roles of three-miRNA signature in some cancer-associated biological processes and pathways. In conclusion, we developed a novel three-miRNA signature that performed reliable prognostic for patient survival with ccRCC, it might facilitate ccRCC patients counseling and individualize management.     

Urinary exosome miR‐30c‐5p as a biomarker of clear cell renal cell carcinoma that inhibits progression by targeting HSPA5

Exosome‐derived miRNAs are regarded as biomarkers for the diagnosis and prognosis of many human cancers. However, its function in clear cell renal cell carcinoma (ccRCC) remains unclear. In this study, differentially expressed miRNAs from urinal exosomes were identified using next‐generation sequencing (NGS) and verified using urine samples of ccRCC patients and healthy donors. Then, the exosomes were analysed in early‐stage ccRCC patients, healthy individuals and patients suffering from other urinary system cancers. Thereafter, the target gene of the miRNA was detected. Its biological function was investigated in vitro and in vivo. The results showed that miR‐30c‐5p could be amplified in a stable manner. Its expression pattern was significantly different only between ccRCC patients and healthy control individuals, but not compared with that of other urinary system cancers, which indicated its specificity for ccRCC. Additionally, the overexpression of miR‐30c‐5p inhibited ccRCC progression in vitro and in vivo. Heat‐shock protein 5 (HSPA5) was found to be a direct target gene of miR‐30c‐5p. The depletion of HSPA5 caused by miR‐30c‐5p inhibition reversed the promoting effect of ccRCC growth. In conclusion, urinary exosomal miR‐30c‐5p acts as a potential diagnostic biomarker of early‐stage ccRCC and may be able to modulate the expression of HSPA5, which is correlated with the progression of ccRCC.     

Identification of a glycolysis-related lncRNA prognostic signature for clear cell renal cell carcinoma

Background: The present study investigated the independent prognostic value of glycolysis-related long noncoding (lnc)RNAs in clear cell renal cell carcinoma (ccRCC).Methods: A coexpression analysis of glycolysis-related mRNAs–long noncoding RNAs (lncRNAs) in ccRCC from The Cancer Genome Atlas (TCGA) was carried out. Clinical samples were randomly divided into training and validation sets. Univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression analyses were performed to establish a glycolysis risk model with prognostic value for ccRCC, which was validated in the training and validation sets and in the whole cohort by Kaplan–Meier, univariate and multivariate Cox regression, and receiver operating characteristic (ROC) curve analyses. Principal component analysis (PCA) and functional annotation by gene set enrichment analysis (GSEA) were performed to evaluate the risk model.Results: We identified 297 glycolysis-associated lncRNAs in ccRCC; of these, 7 were found to have prognostic value in ccRCC patients by Kaplan–Meier, univariate and multivariate Cox regression, and ROC curve analyses. The results of the GSEA suggested a close association between the 7-lncRNA signature and glycolysis-related biological processes and pathways.Conclusion: The seven identified glycolysis-related lncRNAs constitute an lncRNA signature with prognostic value for ccRCC and provide potential therapeutic targets for the treatment of ccRCC patients.     

Systematic analysis of alternative splicing signature unveils prognostic predictor for kidney renal clear cell carcinoma

There is growing evidence that alternative splicing (AS) plays an important role in cancer development. However, a comprehensive analysis of AS signatures in kidney renal clear cell carcinoma (KIRC) is lacking and urgently needed. It remains unclear whether AS acts as diagnostic biomarkers in predicting the prognosis of KIRC patients. In the work, gene expression and clinical data of KIRC were obtained from The Cancer Genome Atlas (TCGA), and profiles of AS events were downloaded from the SpliceSeq database. The RNA sequence/AS data and clinical information were integrated, and we conducted the Cox regression analysis to screen survival-related AS events and messenger RNAs (mRNAs). Correlation between prognostic AS events and gene expression were analyzed using the Pearson correlation coefficient. Protein-protein interaction analysis was conducted for the prognostic AS-related genes, and a potential regulatory network was built using Cytoscape (version 3.6.1). Meanwhile, functional enrichment analysis was conducted. A prognostic risk score model is then established based on seven hub genes (KRT222, LENG8, APOB, SLC3A1, SCD5, AQP1, and ADRA1A) that have high performance in the risk classification of KIRC patients. A total 46,415 AS events including 10,601 genes in 537 patients with KIRC were identified. In univariate Cox regression analysis, 13,362 survival associated AS events and 8,694 survival-specific mRNAs were detected. Common 3,105 genes were screen by overlapping 13,362 survival associated AS events and 8,694 survival-specific mRNAs. The Pearson correlation analysis suggested that 13 genes were significantly correlated with AS events (Pearson correlation coefficient >0.8 or <−0.8). Then, We conducted multivariate Cox regression analyses to select the potential prognostic AS genes. Seven genes were identified to be significantly related to OS. A prognostic model based on seven genes was constructed. The area under the ROC curve was 0.767. In the current study, a robust prognostic prediction model was constructed for KIRC patients, and the findings revealed that the AS events could act as potential prognostic biomarkers for KIRC.     

Plasmatic carbonic anhydrase IX as a diagnostic marker for clear cell renal cell carcinoma

Carbonic anhydrase (CA, EC 4.2.1.1) IX is regarded as a tumour hypoxia marker and CA inhibitors have been proposed as a new class of antitumor agents, with one such agent in Phase II clinical trials. The expression of some CAs, in particular the isoforms CA IX and CA XII, has been correlated with tumour aggressiveness and progression in several cancers. The aim of this study was to evaluate the possibility that CA IX could represent a marker related to clear cell Renal Cell Carcinoma (ccRCC). Bcl-2 and Bax, and the activity of caspase-3, evaluated in tissue biopsies from patients, were congruent with resistance to apoptosis in ccRCCs with respect to healthy controls, respectively. In the same samples, the CA IX and pro-angiogenic factor VEGF expressions revealed that both these hypoxia responsive proteins were strongly increased in ccRCC with respect to controls. CA IX plasma concentration and CA activity were assessed in healthy volunteers and patients with benign kidney tumours and ccRCCs. CA IX expression levels were found strongly increased only in plasma from ccRCC subjects, whereas, CA activity was found similarly increased both in plasma from ccRCC and benign tumour patients, compared to healthy volunteers. These results show that the plasmatic level of CA IX, but not the CA total activity, can be considered a diagnostic marker of ccRCCs. Furthermore, as many reports exist relating CA IX inhibition to a better outcome to anticancer therapy in ccRCC, plasma levels of CA IX could be also predictive for response to therapy.     

A large-scale transcriptome analysis identified ELANE and PRTN3 as novel methylation prognostic signatures for clear cell renal cell carcinoma

DNA methylation was involved in the progress of many types of cancer including clear cell renal cell carcinomas (ccRCCs). This study aimed to identify the prognostic DNA methylation biomarkers for the ccRCCs by a large-scale RNA-seq analysis. The DNA methylation data and the corresponding clinical information of the patients with ccRCCs were extracted from TCGA database and randomly divided into the training group and the validation group. The differentially expressed CpG sites and the survival-related CpG sites were further identified, which was combined into CpG sites pair and followed by screening the survival-related pairs. The C-index and the forward search algorithms were constructed to identify the prognostic signatures for the patients with ccRCCs. The prognostic signatures were verified by the validation dataset and the protein–protein interactions (PPI) network analysis was performed on the CPG sites of the signature. A total of 9,861 differentially expressed CPG sites were identified and 567 CpG sites were found to relate to the overall survival (OS) of the patients with ccRCCs. Besides, 1,146 CPG sites pairs were found to be related to the OS of the ccRCCs samples and the signature composed of seven CpG sites pairs were obtained to predict the prognosis of patients with ccRCCs and the results were verified in the validation dataset. Besides, the PPI network analysis showed that ELANE and PRTN3 gene may be associated with the invasion and metastasis of ccRCCs and could function as potential prognostic and therapeutic signatures for ccRCCs.     

IDUA,NDST1, SAP30L,CRYBA4, and SI as novel prognostic signatures clear cell renal cell carcinoma

Carcinoma of the kidney is one of the most prevalent carcinoma worldwide. The majority types of carcinoma are clear cell renal cell carcinoma (CCRCC), which consist more than 80% of the cases. As a genetically diverse disease, identification of prognosis-related genes has utmost importance in the early diagnosis and prognosis of the CCRCC. In this study, we performed gene expression profiling to identify prognosis-related genes for CCRCC. In addition, we developed and validated a gene signature-based risk score to comprehensively assess the prognostic function of differentially expressed genes. Furthermore, we performed a ROC analysis to identify the optimal cut-off point for classification risk level of the patients. Univariate Cox regression models were used to assess the association between differentially expressed genes in relation to the prognosis of patients with different stages of CCRCC. Five genes were identified significantly differentially expressed in CCRCC and associated with their survival time, namely: IDUA, NDST1, SAP30L, CRYBA4, and SI. A 5-gene signature-based risk score was developed based on the Cox coefficient of the individual genes. The prognostic value of this risk score was validated in an internal testing data set. In summary, a gene-based risk score was identified and validated, which can predict CCRCC patient survival. The potential functions of this gene expression signature and individual differentially expressed gene as prognostic targets of CCRCC were revealed by this study. Furthermore, these findings may have important implications in the understanding of the potential therapeutic method for the CCRCC patients.     

A novel seven-lncRNA signature for prognosis prediction in hepatocellular carcinoma

Liver cancer is still one of the leading causes of cancer-related death worldwide. This study is dedicated to developing a multi–long noncoding RNA (lncRNA) model for risk stratification and prognosis prediction on patients with hepatocellular carcinoma (HCC). We first downloaded lncRNA expression profiles and corresponding clinical information of patients with liver cancer from The Cancer Genome Atlas database. Differentially expressed (DE) lncRNAs between HCC samples and normal samples were identified. In total, 308 patients with HCC were randomly divided into a training group (n = 154) and a testing group (n = 154). Univariate Cox regression and least absolute shrinkage and selection operator Cox regression analyses were performed to select the best survival-related candidates from these DE lncRNAs in the training set. Seven lncRNAs (AC009005.2, RP11-363N22.3, RP11-932O9.10, RP11-572O6.1, RP11-190C22.8, RP11-388C12.8, and ZFPM2-AS1) were finally identified and used to construct a seven-lncRNA signature. The signature could classify patients into high-risk and low-risk groups with significantly different overall survival. The area under the curve of receiver operating characteristic curve for the signature to predict 5-year survival reached more than 0.75. Besides, the prognostic value of the seven-lncRNA signature was independent of conventional clinical factors. The predictive performance of the signature was further validated in the testing set and the whole set. Functional enrichment analysis indicated that the seven prognostic lncRNAs may be involved in several essential biological processes and pathways. The current study demonstrated the potential clinical implications of the seven-lncRNA signature for survival prediction of patients with HCC.     

Systematic analysis of survival-associated alternative splicing signatures in clear cell renal cell carcinoma

Alternative splicing (AS) constitutes a major reason for messenger RNA (mRNA) and protein diversity. Increasing studies have shown a link to splicing dysfunction associated with malignant neoplasia. Systematic analysis of AS events in kidney cancer remains poorly reported. Therefore, we generated AS profiles in 533 kidney renal clear cell carcinoma (KIRC) patients in The Cancer Genome Atlas (TCGA) database using RNA-seq data. Then, prognostic models were developed in a primary cohort (N = 351) and validated in a validation cohort (N = 182). In addition, splicing networks were built by integrating bioinformatics analyses. A total of 11 268 and 8083 AS variants were significantly associated with patient overall survival time in the primary and validation KIRC cohorts, respectively, including STAT1, DAZAP1, IDS, NUDT7, and KLHDC4. The AS events in the primary KIRC cohorts served as candidate AS events to screen the independent risk factors associated with survival in the primary cohort and to develop prognostic models. The area under the curve of the receiver-operator characteristic curve for prognostic prediction in the primary and validation KIRC cohorts was 0.84 and 0.82 at 2500 days of overall survival, respectively. In addition, splicing correlation networks revealed key splicing factors (SFs) in KIRC, such as HNRNPH1, HNRNPU, KHDBS1, KHDBS3, SRSF9, RBMX, SFQ, SRP54, HNRNPA0, and SRSF6. In this study, we analyzed the AS landscape in the TCGA KIRC cohort and detected predictors (prognostic) based on AS variants with high performance for risk stratification of the KIRC cohort and revealed key SFs in splicing networks, which could act as underlying mechanisms.     

Identification of a four immune-related genes signature based on an immunogenomic landscape analysis of clear cell renal cell carcinoma

Renal clear cell carcinoma (ccRCC) is the most common type of renal cell carcinoma, which has strong immunogenicity. A comprehensive study of the role of immune-related genes (IRGs) in ccRCC is of great significance in finding ccRCC treatment targets and improving patient prognosis. In this study, we comprehensively analyzed the expression of IRGs in ccRCC based on The Cancer Genome Atlas datasets. The mechanism of differentially expressed IRGs in ccRCC was analyzed by bioinformatics. In addition, Cox regression analysis was used to screen prognostic related IRGs from differentially expressed IRGs. We also identified a four IRGs signature consisting of four IRGs (CXCL2, SEMA3G, PDGFD, and UCN) through lasso regression and multivariate Cox regression analysis. Further analysis results showed that the four IRGs signature could effectively predict the prognosis of patients with ccRCC, and its predictive power is independent of other clinical factors. In addition, the correlation analysis of immune cell infiltration showed that this four IRGs signature could effectively reflect the level of immune cell infiltration of ccRCC. We also found that the expression of immune checkpoint genes CTLA-4, LAG3, and PD-1 in the high-risk group was higher than that in the low-risk group. Our research revealed the role of IRGs in ccRCC, and developed a four IRGs signature that could be used to evaluate the prognosis of patients with ccRCC, which will help to develop personalized treatment strategies for patients with ccRCC and improve their prognosis. In addition, these four IRGs may be effective therapeutic targets for ccRCC.