Abstract
BACKGROUND: Pancreatic neuroendocrine tumor (pNET) is a heterogeneous tumor originating from pancreatic endocrine cells. Emerging evidence suggests that oxidative stress plays a crucial role in pNET pathogenesis, yet the precise molecular mechanisms and their interplay with the tumor microenvironment remain unclear. This study aims to systematically elucidate how oxidative stress-related pathways drive pNET progression through an integrated multi-omics approach. METHODS: We designed a three-tier analytical strategy to address interconnected scientific questions. First, to identify which oxidative stress-related genes are dysregulated in pNET, we performed differential expression analysis and weighted gene co-expression network analysis (WGCNA) on the GSE73338 dataset (63 pNET samples, 5 controls), intersecting the. results with oxidative stress gene sets to obtain 71 candidate genes. Second, to understand the functional implications of these genes, we conducted GO/KEGG enrichment analysis and constructed protein-protein interaction (PPI) networks, from which we identified BCL2L1 and PHGDH as key hub genes using three independent algorithms. We then assessed their diagnostic value through ROC analysis and built a prognostic nomogram model. Third, to explore how these key genes influence the tumor microenvironment, we performed immune infiltration analysis using CIBERSORTx. Fourth, to reveal upstream regulatory mechanisms, we constructed ceRNA networks and predicted transcription factors. Fifth, to identify potential therapeutic interventions, we conducted drug prediction and molecular docking analyses. Finally, to validate our findings at cellular resolution and understand cellular heterogeneity, we analyzed single-cell RNA sequencing data from GSE256136 (20 samples), identifying cell types, quantifying cell-cell communications, and confirming key gene expression patterns across different cell populations. RESULTS: Our systematic analysis revealed that oxidative stress-related genes in pNET were significantly enriched in the PI3K-Akt signaling pathway, cysteine and methionine metabolism, and HIF-1 signaling pathway. BCL2L1 and PHGDH emerged as central regulators with excellent diagnostic performance (AUC > 0.9). Immune infiltration analysis demonstrated significant alterations in activated dendritic cells, memory B cells, and resting NK cells, which correlated strongly with BCL2L1 and PHGDH expression, suggesting these genes link oxidative stress to immune dysfunction. The ceRNA network centered on KCNQ1OT1 and hsa-miR-15a-5p revealed multi-layered transcriptional and post-transcriptional regulation. Drug prediction identified sertindole and cabozantinib as promising therapeutic candidates. Single-cell analysis identified 11 cell types and confirmed that endocrine cells are the primary site of BCL2L1 and PHGDH dysregulation, with extensive crosstalk between endocrine cells and T cells potentially mediating immune evasion. CONCLUSION: Through integrated multi-omics analysis, we established that oxidative stress pathways may drive pNET progression through a coordinated mechanism involving metabolic reprogramming (via BCL2L1 and PHGDH downregulation), immune microenvironment remodeling (through altered dendritic cell and NK cell function), and complex regulatory networks. BCL2L1 and PHGDH represent potential diagnostic biomarkers and candidate therapeutic targets that require experimental validation, providing new directions for precision medicine in pNET.