Combination of bulk RNA sequencing and scRNA sequencing uncover the molecular characteristics of MAPK signaling in kidney renal clear cell carcinoma

Abstract

The MAPK signaling pathway significantly impacts cancer progression and resistance; however, its functions remain incompletely assessed across various cancers, particularly in kidney renal clear cell carcinoma (KIRC). Therefore, there is an urgent need for comprehensive pan-cancer investigations of MAPK signaling, particularly within the context of KIRC. In this research, we obtained TCGA pan-cancer multi-omics data and conducted a comprehensive analysis of the genomic and transcriptomic characteristics of the MAPK signaling pathway. For in-depth investigation in KIRC, status of MAPK pathway was quantitatively estimated by ssGSEA and Ward algorithm was utilized for cluster analysis. Molecular characteristics and clinical prognoses of KIRC patients with distinct MAPK activities were comprehensively explored using a series of bioinformatics algorithms. Subsequently, a combination of LASSO and COX regression analyses were utilized sequentially to construct a MAPK-related signature to help identify the risk level of each sample. Patients in the C1 subtype exhibited relatively higher levels of MAPK signaling activity, which were associated with abundant immune cell infiltration and favorable clinical outcomes. Single-cell RNA sequencing (scRNA-seq) analysis of KIRC samples identified seven distinct cell types, and endothelial cells in tumor tissues had obviously higher MAPK scores than normal tissues. The immunohistochemistry results indicated the reduced expression levels of PAPSS1, MAP3K11, and SPRED1 in KIRC samples. In conclusion, our study represents the first integration of bulk RNA sequencing and single-cell RNA sequencing to elucidate the molecular characteristics of MAPK signaling in KIRC, providing a solid foundation for precision oncology.