Comprehensive analysis of cellular senescence-related genes in the prognosis, tumor microenvironment, and immunotherapy/chemotherapy of clear cell renal cell carcinoma.

The transcriptome public database and advances in biological discoveries contributed to significant progresses in identifying the drivers of cancer progression. Cellular senescence (CS) is considered as a leading factor resulting in cancer development. The purpose of this study was to explore the significance of CS-related genes in the molecular classification and survival outcome of clear cell renal cell carcinoma (ccRCC). CS-related genes were obtained from the CellAge database, and patients from TCGA-KIRC dataset and ICGC dataset were clustered by ConsesusClusterPlus. The characteristics of overall survival (OS), genomic variation, and tumor microenvironment (TME) of each cluster were analyzed. Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression analysis was conducted to develop a CS-related risk model to score ccRCC patients and assess the risk scores in predicting patients' response to immunotherapy and chemotherapy. A nomogram based on the risk model was established to improve the risk stratification of patients. CcRCC was divided into three molecular subtypes based on CS-related genes. The three molecular phenotypes showed different OS and clinical manifestations, mutation patterns, and TME states. Five genes were obtained from nine differentially expressed CS-related genes in the three molecular subtypes to develop a risk model. Patients with ccRCC were divided into high- and low-risk subgroups. The former showed an unfavorable OS, with a significantly higher genomic variation rate, TME score, and numerous immune checkpoint expressions when compared to the low-risk subgroup. Risk score reflected the response of patients to axitinib, bortezomib, sorafenib, sunitinib, and temsirolimus. In general, CS-related genes divided ccRCC into three molecular subtypes with distinct OS, mutation patterns, and TME states. The risk model based on the five CS-related genes can predict the prognosis and therapeutic outcome of ccRCC patients, providing a theoretical basis for further study on the molecular mechanism of CS-related ccRCC.

Lu C ，Wang Y ，Nie L ，Chen L ，Li M ，Qing H ，Li S ，Wu S ，Wang Z ... -  《Frontiers in Immunology》

Computational construction of TME-related lncRNAs signature for predicting prognosis and immunotherapy response in clear cell renal cell carcinoma.

The tumor microenvironment (TME) is closely related to clear cell renal cell carcinoma (ccRCC) prognosis, and immunotherapy response. In current study, comprehensive bio-informative analysis was adopted to construct a TME-related lncRNA signature for immune checkpoint inhibitors (ICIs) and targeted drug responses in ccRCC patients. The TME mRNAs were screened following the immune and stromal scores with the data from GSE15641, GSE29609, GSE36895, GSE46699, GSE53757, and The Cancer Genome Atlas (TCGA)-kidney renal clear cell carcinoma (KIRC). And the TME-related lncRNAs were recognized using correlation analysis. The TME-related lncRNAs prognostic model was constructed using the training dataset. Kaplan-Meier analysis, principal-component analysis, and time-dependent receiver operating characteristic were used to evaluate the risk model. The immune cell infiltration in TME was evaluated using the single-sample gene set enrichment analysis (ssGSEA), ESTIMATE, and microenvironment cell populations counter algorithm. The immunophenoscore (IPS) was used to assess the response to immunotherapy with the constructed model. In the current study, 364 TME-related lncRNAs were selected based on the integrated bioinformatical analysis. Six TME-related lncRNAs (LINC00460, LINC01094, AC008870.2, AC068792.1, and AC007637.1) were identified as the prognostic signature in the training dataset and subsequently verified in the testing and entire datasets. Patients in the high-risk group exhibited poor overall survival and disease-free survival than those in the low-risk group. The 1-, 3-, and 5-year areas under the curves of the prognostic signature in the entire dataset were 0.704, 0.683, and 0.750, respectively. The risk score independently predicted ccRCC survival based on univariate and multivariate Cox regression. GSEA analysis suggested that the high-risk group was concentrated on immune-related pathways. The high-risk group were characterized by high immune cell infiltration, high TMB and somatic mutation counters, high IPS-PD-1 + CTLA4 scores, and immune checkpoints expression upregulation, reflecting the higher ICIs response. The half inhibitory concentrations of sunitinib, temsirolimus, and rapamycin were low in the high-risk group. The TME-related lncRNAs signature constructed could reliably predict the prognosis and immunotherapy response and targeted ccRCC patients' therapy.

Zhou L ，Fang H ，Guo F ，Yin M ，Long H ，Weng G ... -  《-》

Novel immune-related signature based on immune cells for predicting prognosis and immunotherapy response in clear cell renal cell carcinoma.

Clear cell renal cell carcinoma (ccRCC) is the most common malignant tumor of the kidney and is characterized by poor prognosis. We sought to build an immune-related prognostic signature and investigate its relationship with immunotherapy response in ccRCC. Immune-related genes were identified by ssGSEA and WGCNA. The prognostic signature was conducted via univariate, least absolute shrinkage and selection operator, and multivariable Cox regression analyses. Kaplan-Meier analysis, PCA, t-SNE, and ROC were used to evaluate the risk model. A total of 119 immune-related genes associated with prognosis were screened out. Six immune-related genes (CSF1, CD5L, AIM2, TIMP3, IRF6, and HHLA2) were applied to construct a prognostic signature for KIRC. Kaplan-Meier analysis showed that patients in high-risk group had a poorer survival outcome than in low-risk group. The 1-, 3- and 5-year AUC of the prognostic signature was 0.754, 0.715, and 0.739, respectively. Univariate and multivariate Cox regression models demonstrated that the risk signature was an independent prognostic factor for KIRC survival. GSEA analysis suggested that the high-risk group was concentrated on immune-related pathways. The high-risk group with more regulatory T-cell infiltration showed a higher expression of immune negative regulation genes. The risk score had positively relationship with TIDE score and negatively with the response of immunotherapy. The IC50 values of axitinib, sunitinib, sorafenib, and temsirolimus were lower in the high-risk group. Our study defined a robust signature that may be promising for predicting clinical outcomes and immunotherapy and targeted therapy response in ccRCC patients.

Zhou L ，Fang H ，Yin M ，Long H ，Weng G ... -  《-》

Comprehensive analysis of LD-related genes signature for predicting prognosis and immunotherapy response in clear cell renal cell carcinoma.

Jia Y ，Dong X ，Yang F ，Zhou L ，Long H ... -  《BMC Nephrology》

Construction of an interferon regulatory factors-related risk model for predicting prognosis, immune microenvironment and immunotherapy in clear cell renal cell carcinoma.

Interferon regulatory factors (IRFs) played complex and essential roles in progression, prognosis, and immune microenvironment in clear cell renal cell carcinoma (ccRCC). The purpose of this study was to construct a novel IRFs-related risk model to predict prognosis, tumor microenvironment (TME) and immunotherapy response in ccRCC. Multi-omics analysis of IRFs in ccRCC was performed based on bulk RNA sequencing and single cell RNA sequencing data. According to the expression profiles of IRFs, the ccRCC samples were clustered by non-negative matrix factorization (NMF) algorithm. Then, least absolute shrinkage and selection operator (LASSO) and Cox regression analyses were applied to construct a risk model to predict prognosis, immune cells infiltration, immunotherapy response and targeted drug sensitivity in ccRCC. Furthermore, a nomogram comprising the risk model and clinical characteristics was established. Two molecular subtypes with different prognosis, clinical characteristics and infiltration levels of immune cells were identified in ccRCC. The IRFs-related risk model was developed as an independent prognostic indicator in the TCGA-KIRC cohort and validated in the E-MTAB-1980 cohort. The overall survival of patients in the low-risk group was better than that in the high-risk group. The risk model was superior to clinical characteristics and the ClearCode34 model in predicting the prognosis. In addition, a nomogram was developed to improve the clinical utility of the risk model. Moreover, the high-risk group had higher infiltration levels of CD8+ T cell, macrophages, T follicular helper cells and T helper (Th1) cells and activity score of type I IFN response but lower infiltration levels of mast cells and activity score of type II IFN response. Cancer immunity cycle showed that the immune activity score of most steps was remarkably higher in the high-risk group. TIDE scores indicated that patients in the low-risk group were more likely responsive to immunotherapy. Patients in different risk groups showed diverse drug sensitivity to axitinib, sorafenib, gefitinib, erlotinib, dasatinib and rapamycin. In brief, a robust and effective risk model was developed to predict prognosis, TME characteristics and responses to immunotherapy and targeted drugs in ccRCC, which might provide new insights into personalized and precise therapeutic strategies.

Pan H ，Lu W ，Zhang M ，Liu C ... -  《Frontiers in Oncology》