Key molecular mechanisms of mitochondrial metabolic pathways in specific cell subpopulations of pancreatic cancer based on scRNA-seq and bulk RNA-seq

Abstract

Background Pancreatic cancer (PC) is a highly malignant tumor with increasing incidence, mortality, and a low five-year survival rate. Mitochondrial metabolic reprogramming plays a crucial role in tumor development, but the molecular mechanisms in different cell subpopulations of PC remain unclear. This study aims to integrate single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq to explore mitochondrial metabolism in PC. Methods Transcriptome datasets (GSE183795, GSE16515, GSE197177, and TCGA-PAAD) were downloaded from the GEO and TCGA databases. Differentially expressed genes (DEGs) were identified using the limma package for bulk RNA-seq and the Seurat package for scRNA-seq. Weighted gene co-expression network analysis (WGCNA) was performed to identify key modules and hub genes. Machine learning algorithms screened key genes, and functional enrichment analysis was conducted using clusterProfiler. PPI, ceRNA, and transcription factor networks explored gene regulation. Immune infiltration analysis, drug prediction, and molecular docking were conducted on key genes. The prognostic value of the key genes was evaluated using clinical data from TCGA. Results A total of 1238 bulk RNA-DEGs and 8231 scRNA-DEGs in fibroblasts were identified. Integration of both datasets revealed 536 DEGs. WGCNA identified 6 modules associated with PC. By intersecting DEGs, hub genes, and mitochondrial metabolism-related genes, 18 candidate genes were obtained. These genes were enriched in glucose metabolism and mitochondrial outer membrane pathways. Three key genes (IFI27, PKM, and RSAD2) were selected using machine learning. PPI, ceRNA, and transcription factor networks provided regulatory insights. Immune infiltration analysis showed significant differences in immune cells, particularly in T cells CD4 memory resting and macrophages M2. Drug prediction and molecular docking identified potential drugs for these genes. Survival analysis indicated that high expression of these genes correlated with poor prognosis. Conclusion This study integrates scRNA-seq and bulk RNA-seq data to identify three key genes (IFI27, PKM, and RSAD2) and their immune-related mechanisms in PC. These findings offer new insights into the pathogenesis and potential therapeutic targets for PC.