Dataset Information

A Novel Prognostic Signature Associated With the Tumor Microenvironment in Kidney Renal Clear Cell Carcinoma.

ABSTRACT:

Background

The tumor microenvironment (TME) is a complex and evolving environment, and the tumor immune microenvironment in kidney renal clear cell carcinoma (KIRC) has a strong suppressive profile. This study investigates the potential prognostic role and value of genes of the tumor microenvironment in KIRC.

Methods

The transcriptome sequencing data of 530 cases and 39 cases of KIRC and the corresponding clinical prognosis information were downloaded from TCGA data and GEO data, respectively, and TME-related gene expression profiles were extracted. A prognostic signature was constructed and evaluated using univariate Cox regression analysis and LASSO regression analysis. Gene set enrichment analysis (GSEA) was used to obtain the biological process of gene enrichment in patients with high and low-risk groups.

Results

A prognostic signature consisting of eight TME-related genes (LRFN1, CSF1, UCN, TUBB2B, SERPINF1, ADAM8, ABCB4, CCL22) was constructed. Kaplan-Meier survival analysis yielded significantly lower survival times for patients in the high-risk group than in the low-risk group, and the AUC values for the ROC curves of this prognostic signature were essentially greater than 0.7, and univariate and multifactorial Cox regression analyses indicated that the risk score was independent risk factors for KIRC prognosis. GSEA analysis showed that immune-related biological processes were enriched in the high-risk group and that risk values were strongly associated with multiple immune cell scores and immune checkpoint-related genes (PDCD1, CTLA4).

Conclusions

The prognostic signature can accurately predict the prognosis of KIRC patients, which may provide new ideas for future precision immunotherapy of KIRC.

SUBMITTER: Pei D

PROVIDER: S-EPMC9290677 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background: M7G modification is extremely vital for the development of many cancers, especially tumor immunity. M7G modification is a novel functional regulator of miRNA, and the researches on m7G-related miRNAs in kidney renal clear cell carcinoma (KIRC) are still insufficient. This research aims to establish a risk signature on the foundation of m7G-associated miRNAs, which can precisely forecast the prognosis of KIRC patients. Methods: Transcriptome data and clinical data used in this study come from The Cancer Genome Atlas database. Our team utilized univariable Cox, Lasso and multivariable Cox analyses to construct a m7G-associated miRNAs risk signature that can forecast the prognosis of KIRC patients. Kaplan-Meier method, time-dependent receiver operating characteristic (ROC) curve, and the independent analysis of risk signatures were employed to verify the predictability and accuracy of the risk signature. Subsequently, based on CIBERSORT, ESTIMATE and ssGSEA algorithms, we speculated the potential impact of the proposed risk signature on tumor immune microenvironment. Ultimately, by virtue of the risk signature and tumor immunity, the hub genes affecting the prognosis of KIRC patients were screened out. Results: Our team established and verified a prognostic signature comprising 7 m7G-associated miRNAs (miR-342-3p, miR-221-3p, miR-222-3p, miR-1277-3p, miR-6718-5p, miR-1251-5p, and miR-486-5p). The results of the Kaplan-Meier survival analysis revealed that the prognosis of KIRC sufferers in the high-risk group was often unsatisfactory. The accuracy of the prediction ability of the risk signature was verified by calculating the area under the ROC curve. Univariate-multivariate Cox analyses further showed that this risk signature could be utilized as an independent prognosis-related biomarker for KIRC sufferers. The results of the immune analysis revealed that remarkable diversities existed in immune status and tumor microenvironment between high-risk and low-risk groups. On the foundation of the proposed risk signature and other clinical factors, a nomogram was established to quantitatively forecast the survival of KIRC sufferers at 1, 3 and 5 years. Conclusion: Based on m7G-related miRNAs, a risk signature was successfully constructed, which could precisely forecast the prognosis of sufferers and guide personalized immunotherapy for KIRC patients.

Project description:Background: Necroptosis, a cell death of caspase-independence, plays a pivotal role in cancer biological regulation. Although necroptosis is closely associated with oncogenesis, cancer metastasis, and immunity, there remains a lack of studies determining the role of necroptosis-related genes (NRGs) in the highly immunogenic cancer type, kidney renal clear cell carcinoma (KIRC). Methods: The information of clinicopathology and transcriptome was extracted from TCGA database. Following the division into the train and test cohorts, a three-NRGs (TLR3, FASLG, ZBP1) risk model was identified in train cohort by LASSO regression. The overall survival (OS) comparison was conducted between different risk groups through Kaplan-Meier analysis, which was further validated in test cohort. The Cox proportional hazards regression model was introduced to assess its impact of clinicopathological factors and risk score on survival. ESTIMATE and CIBERSORT algorithms were introduced to evaluate immune microenvironment, while enrichment analysis was conducted to explore the biological significance. Correlation analysis was applied for the correlation assessment between checkpoint gene expression and risk score, between gene expression and therapeutic response. Gene expressions from TCGA were verified by GEO datasets and immunohistochemistry (IHC) analysis. Results: This NRGs-related signature predicted poorer OS in high-risk group, which was also verified in test cohort. Risk score could also independently predict survival outcome of KIRC. Significant changes were also found in immune microenvironment and checkpoint gene expressions between different risk groups, with immune functional enrichment in high-risk group. Interestingly, therapeutic response was correlated with the expressions of NRGs. The expressions of NRGs from TCGA were consistent with those from GEO datasets and IHC analysis. Conclusion: The NRGs-related signature functions as a novel prognostic predictor of immune microenvironment and therapeutic response in KIRC.

Project description:BackgroundRecently, an increasing number of studies have uncovered the aberrant expression of methyltransferase-like family (METTL) plays an important role in tumorigenesis, such as METTL3 (an m6A writer). In our recent work, we discovered METTL24 expression was highly associated with the hazard ratio (HR) of kidney renal clear cell carcinoma (KIRC) compared to other tumors, implying a special function of METTL24 in KIRC carcinogenesis. Until now, the functions and mechanisms of METTL24 in KIRC have remained mostly unknown.MethodsThe mRNA expression of METTL24 in KIRC was analyzed using the TIMER 2.0, GEPIA, and UALCAN databases. The immunohistochemical assay was performed to validate METTL24 expression in our self-built Chinese cohort (n tumor = 88, n normal = 85). The gene set enrichment analysis (GSEA) was used to investigate the biological processes in which METTL24 might be engaged. The Spearman analysis was used to evaluate the expression correlations between METTL24 and a range of immunological variables, and the effects of METTL24 on the infiltration levels of multiple immune cells were explored using TCGA data. The upstream transcription factors of METTL24 were screened through a multi-omics analysis.ResultsMETTL24 expression in KIRC tissues was significantly decreased compared to normal adjacent kidney tissues, which was associated with the lower survival rate of KIRC patients. METTL24 potentially participated in the immune-relevant biological processes such as cytokine binding, NF-kappa B binding, MHC protein complex, and interleukin-12 action. Besides, METTL24 expression was linked to a number of immune checkpoints, cytokines, chemokines, and chemokine receptors, and also correlated with the infiltration levels of 10 types of immune cells in KIRC. Meanwhile, METTL24 expression differently affected the overall survival rates (OS) of KIRC patients with high or low levels of immune infiltration. Finally, CTCF and EP300 were discovered to be the probable transcription factors of METTL24 in KIRC.ConclusionThis study revealed that METTL24 might serve as a prognostic marker in KIRC and as one immune-relevant target for clinical treatment.

Project description:Abnormal autophagic levels have been implicated in the pathogenesis of multiple cancers, however, its role in tumors is complex and has not yet been explored clearly. Hence, we aimed to explore the prognostic values of autophagy-related genes (ARGs) for kidney renal clear cell carcinoma (KIRC). Differentially expressed ARGs and transcription factors (TFs) were identified in KIRC patients obtaining from the The Cancer Genome Atlas (TCGA) database. Then, networks between TFs and ARGs, gene ontology functional annotations and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis were conducted. Next, we performed consensus clustering, COX regression analysis and Lasso regression analysis to identify the prognostic ARGs. Finally, an individual prognostic index (PI, riskScore) was established. Based on TCGA cohort and ArrayExpress cohort, Survival analysis, ROC curve, independent prognostic analysis, and clinical correlation analysis were also performed to evaluate this PI. Based on differentially expressed ARGs, KIRC patients were successfully divided into two clusters (P = 5.916e-04). AS for PI, it was constructed based on 11 ARGs and significantly classified KIRC patients into high-risk group and low-risk group in terms of OS (P = 4.885e-15 for TCGA cohort, P = 6.366e-03 for ArrayExpress cohort). AUC of its ROC curve reached 0.747 for TCGA cohort and 0.779 for ArrayExpress cohort. What's more, this PI was proven to be a valuable independent prognostic factor in both univariate and multivariate COX regression analysis (P < .001). Prognostic nomograms were also performed to visualize the relationship between individual predictors and survival rates in patients with KIRC. By means of connectivity map database, emetine, cephaeline and co-dergocrine mesilate related to ARGs were found to be negatively correlated with KIRC. This study provided an effective PI for KIRC and also displayed networks between TFs and ARGs. KIRC patients were successfully divided into two clusters based on differentially expressed ARGs. Besides, small molecule drugs related to ARGs were also identified for KIRC.

Project description:BackgroundRNA N6-methyladenosine (m6A) has been found to have a critical impact on clear cell renal cell carcinoma (ccRCC) by affecting the tumor microenvironment (TME) and immune cell (IC) infiltration and is related to the treatment and survival rate of patients with ccRCC. However, the mechanism of m6A in TME and IC infiltration remained unclear.MethodsNonnegative Matrix Factorization (NMF) clustering was performed on 650 ccRCC cases from the Cancer Genome Atlas (TCGA) and the Gene-Expression Omnibus (GEO) datasets. The immune infiltration was generated by the single-sample gene-set enrichment analysis (ssGSEA) algorithm. Survival analyses were performed using the Kaplan-Meier method, and the significance of the differences was determined using the log-rank test. The m6A score was constructed based on the expression of m6A regulators to quantify m6A modification. The package "survminer R" was employed to layer patients' low and high scores groups and predict the immunotherapy response.ResultsThree different patterns of m6A modification were established, and significant differences in TME and IC infiltration features were found in these three patterns. Survival analysis demonstrated that m6A cluster A and m6A gene cluster A experienced a longer survival time. Evaluation of m6A modification patterns in individual tumors was initiated by the m6A score. The low m6A score subtype was characterized by increased tumor mutation burden (TMB) and immune infiltration, whereas a high m6A score with a lack of immune cell infiltration showed significantly better overall survival. m6A score was also associated with the expression of programmed cell death protein 1 (PD-L1) and cytotoxic T lymphocyte antigen 4 (CTLA-4). Patients in the high m6A score group had high PD-L1 expression and low CTLA-4 expression. Significant differences in prognosis were identified among types of different TMB and m6A scores, where low TMB and high m6A score had longer survival time.ConclusionsThis research indicated that m6A modification greatly affected TME and IC infiltration. Physicians can develop practical immunotherapy strategies for patients with ccRCC by evaluating m6A-associated genes.