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.     

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.     

Identification of microenvironment-related genes with prognostic value in clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the most common type of kidney tumor. Previous studies have shown that the interaction between tumor cells and microenvironment has an important impact on prognosis. Immune and stromal cells are two vital components of the tumor microenvironment. Our study aimed to better understand and explore the genes involved in immune/stromal cells on prognosis. We used the Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data algorithm to calculate immune/stromal scores. According to the scores, we divided ccRCC patients from The Cancer Genome Atlas database into low and high groups and identified the genes which were differentially expressed and significantly associated with prognosis. The result of functional enrichment analysis and protein-protein interaction networks indicated that these genes mainly were involved in extracellular matrix and regulation of cellular activities. Then another independent cohort from the International Cancer Genome Consortium database was used to validate these genes. Finally, we acquired a list of microenvironment-related genes that can predict prognosis for ccRCC patients.     

利用TCGA数据库构建肾透明细胞癌相关miRNA预后模型

利用TCGA数据库中肾透明细胞癌的miRNA与mRNA数据及临床信息,构建由miRNA组成的预后风险评分模型,并筛选与生存预后相关的miRNA-mRNA调控关系对,为研究提供理论依据。下载并整理TCGA[JP+1]数据库中肾透明细胞癌的miRNA与mRNA数据;对数据进行差异分析,将差异表达的miRNA与临床信息进行合并,利用单因素与多因素Cox回归分析,构建预后模型并进行模型评价;对模型中的miRNA进行靶基因预测,结果与差异表达的mRNA进行取交集,构建miRNA-mRNA调控网络;对网络中的mRNA进行生存分析,筛选生存相关的miRNA-mRNA调控关系对。共得到49个差异表达的miRNA与3 613个差异表达的mRNA;预后模型计算公式为:风险值(risk score)=hsa-miR-21-5p表达量×0.603＋hsa-miR-1251-5p表达量×-0.093;调控网络中共纳入31个miRNA-mRNA调控关系对;对mRNA进行生存分析,共得到7个有价值的关系对。所构建预后模型可有效预测肾透明细胞癌患者生存预后情况,筛选到的miRNA-mRNA调控关系对可为相关研究与治疗提供参考。     

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.     

Expression of CIDE proteins in clear cell renal cell carcinoma and their prognostic significance

Clear cell renal cell carcinoma (ccRCC) is the major and aggressive subtype of renal cell carcinoma. It is known to derive its histologic appearance from accumulation of abundant lipids and glycogens. The cell death-inducing DFF45-like effector (CIDE) family has been characterized as the lipid droplet proteins involved in the metabolism of lipid storage droplets. The purpose of this study was to evaluate the expression of CIDE proteins in ccRCC cells and to investigate their prognostic significance. We examined consecutive patients with sporadic ccRCC, who underwent nephrectomy, to measure their mRNA and protein expression of CIDE proteins. We found that Cidec and ADRP expression were significantly up-regulated in ccRCC, compared with normal kidney tissues. Cideb was down-regulated. We also found that Cideb was expressed more in low-grade ccRCC than in high-grade tumors. To further clarify the relationship between Cideb expression and patient prognosis, we evaluated 57 ccRCC patients followed up for 120 months. Reduced ccRCC Cideb expression was associated with a higher Fuhrman nuclear grade. Patients with high Cideb expression had better overall survival rate than those with low expression (p < 0.05). Cideb expression was an independent predictor of survival (p = 0.001). Although the biologic function of Cideb in ccRCC remains unknown, the expression level of Cideb might be a novel predictor of prognosis in ccRCC.     

Stratification of clear cell renal cell carcinoma by signaling pathway analysis

Investigation of cell signaling pathways in 16 clear cell renal cell carcinomas to identify groups based on commonly shared phosphorylation-driven signaling networks. Using laser capture microdissection and reverse-phase protein arrays, we profiled 75 key nodes spanning signaling pathways important in tumorigenesis. Analysis revealed significantly different (P < 0.05) signaling levels for 27 nodes between two groups of samples, designated A (4 samples; high EGFR, RET, and RASGFR1 levels, converging to activate AKT/mTOR) and B (12 samples; high ERK1/2 and STAT phosphorylation). Group B was further partitioned into groups C (7 samples; elevated expression of LC3B) and D (5 samples; activation of Src and STAT). Network analysis indicated that group A was characterized by signaling pathways related to cell cycle and proliferation, and group B by pathways related to cell death and survival. Homogeneous clear cell renal cell carcinomas could be stratified into at least two major functional groups.     

肾透明细胞癌预后相关miRNAs的生物信息学分析

目的寻找可作为肾透明细胞癌（ccRCC）生物标志物的miRNA,以及ccRCC与正常组织间miRNA差异表达情况。 方法利用TCGA数据库下载ccRCC中miRNA表达数据,分析肿瘤与正常组织间差异表达miRNA。使用Kaplan-Meier曲线对患者进行生存分析,筛选出表达情况与临床预后相关的miRNA。通过生物信息学对miRNA的靶基因进行预测,然后运用FunRich软件和ClueGO对靶基因进行GO和KEGG富集分析。 结果通过TCGA数据库分析发现,ccRCC较正常组织差异表达miRNA共54个,其中上调33个,下调21个。通过生存分析发现hsa-miR-21和hsa-miR-155与患者预后相关,P≤0.05。进一步通过Perl软件在Targetscan、miRDB、miRTarBase、miRPath这四个数据库中预测miRNA靶基因并将结果取交集,共发现129个靶基因。GO和KEGG分析结果表明,这些靶基因主要与转录因子活性、信号转导以及FoxO、TNF等信号通路密切相关。 结论通过生物信息学分析发现了ccRCC与正常组织的差异表达miRNA;其中hsa-miR-21和hsa-miR-155与患者总体生存率相关,并通过调控靶基因参与相关的信号通路进而影响ccRCC的发生发展进程,提示hsa-miR-21和hsa-miR-155可能是ccRCC潜在的生物标志物。     

Single-cell analysis reveals metastatic cell heterogeneity in clear cell renal cell carcinoma

Renal cell carcinoma (RCC) is one of the leading causes of cancer-related death worldwide. Tumour metastasis and heterogeneity lead to poor survival outcomes and drug resistance in patients with metastatic RCC (mRCC). In this study, we aimed to assess intratumoural heterogeneity (ITH) in mRCC cells by performing a combined analysis of bulk data and single-cell RNA-sequencing data, and develop novel biomarkers for prognosis prediction on the basis of the potential molecular mechanisms underlying tumorigenesis. Eligible single-cell cohorts related to mRCC were acquired using the Gene Expression Omnibus (GEO) dataset to identify potential mRCC subpopulations. We then performed gene set variation analysis to understand the differential function in primary RCC and mRCC samples. Subsequently, we applied weighted correlation network analysis to identify coexpressing gene modules that were related to the external trait of metastasis. Protein-protein interactions were used to screen hub subpopulation-difference (sub-dif) markers (ACTG1, IL6, CASP3, ACTB and RAP1B) that might be involved in the regulation of RCC metastasis and progression. Cox regression analysis revealed that ACTG1 was a protective factor (HR < 1), whereas the other four genes (IL6, CASP3, ACTB and RAP1B) were risk factors (HR > 1). Kaplan-Meier survival analysis suggested the potential prognostic value of these sub-dif markers. The expression of sub-dif markers in mRCC was further evaluated in clinical samples by immunohistochemistry (IHC). Additionally, the genetic features of sub-dif marker expression patterns, such as genetic variation profiles, correlations with tumour-infiltrating lymphocytes (TILs), and targeted signalling pathway activities, were assessed in bulk RNA-seq datasets. In conclusion, we established novel subpopulation markers as key prognostic factors affecting EMT-related signalling pathway activation in mRCC, which could facilitate the implementation of a treatment for mRCC patients.     

Familial non-syndromic clear cell renal cell carcinoma

The diagnosis of familial non-syndromic clear cell renal cell carcinoma is one of exclusion. In families presenting with clear cell RCC a germline VHL mutation and a constitutional translocation of chromosome 3 must be excluded before familial non-syndromic clear cell RCC can be diagnosed. Large familial non-syndromic clear cell RCC kindreds are uncommon and a predisposing gene has not been identified. However inheritance is autosomal dominant in most cases and age at onset is earlier than in sporadic cases. Recognition and appropriate screening of familial non-syndromic clear cell RCC cases will reduce morbidity and mortality. Large scale collaborative linkage studies may provide a basis for the identification of familial non-syndromic clear cell RCC susceptibility gene(s).     

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.     

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.     

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.     

Integrated glycoproteomic characterization of clear cell renal cell carcinoma

1. Download : Download high-res image (269KB)
2. Download : Download full-size image

     

Construction and validation of a novel ten miRNA-pair based signature for the prognosis of clear cell renal cell carcinoma

BackgroundClear cell renal cell carcinoma (ccRCC) is the most predominate pathological subtype of renal cell carcinoma, causing a recurrence or metastasis rate as high as 20% to 40% after operation, for which effective prognostic signature is urgently needed.MethodsThe mRNA and miRNA profiles of ccRCC specimens were collected from the Cancer Genome Atlas. MiRNA-pair risk score (miPRS) for each miRNA pair was generated as a signature and validated by univariate and multivariate Cox proportional hazards regression analysis. Functional enrichment was performed, and immune cells infiltration, as well as tumor mutation burden (TMB), and immunophenoscore (IPS) were evaluated between high and low miPRS groups. Target gene-prediction and differentially expressed gene-analysis were performed based on databases of miRDB, miRTarBase, and TargetScan. Multivariate Cox proportional hazards regression analysis was adopted to establish the prognostic model and Kaplan-Meier survival analysis was performed.FindingsA novel 10 miRNA-pair based signature was established. Area under the time-dependent receiver operating curve proved the performance of the signature in the training, validation, and testing cohorts. Higher TMB, as well as the higher CTLA4-negative PD1-negative IPS, were discovered in high miPRS patients. A prognostic model was built based on miPRS (1 year-, 5 year-, 10 year- ROC-AUC=0.92, 0.84, 0.82, respectively).InterpretationThe model based on miPRS is a novel and valid tool for predicting the prognosis of ccRCC.FundingThis study was supported by research grants from the China National Natural Scientific Foundation (81903972, 82002018, and 82170752) and Shanghai Sailing Program (19YF1406700 and 20YF1406000).