Comprehensive Single-Cell Transcriptomic Profiling Identifies Molecular Subtypes and Prognostic Biomarkers in Esophageal Squamous Cell Carcinoma: Implications for Targeted Therapy

Abstract Background: Esophageal squamous cell carcinoma (ESCC) is a genetically heterogeneous disease with poor clinical outcomes. Identification of biomarkers linked to DNA replication stress may enable improved prognostic risk stratification and guide therapeutic decision making. We performed integrated single-cell RNA sequencing and computational analyses to define the molecular determinants and subtypes underlying ESCC heterogeneity. Methods: Single-cell RNA sequencing was performed on ESCC samples and analyzed using Seurat. Differential gene expression analysis was used to identify esophageal cell phenotypes. DNA replication stress-related genes were intersected with single-cell differential expression data to identify potential prognostic genes, which were used to generate a DNA replication stress (DRS) score. This score and associated genes were evaluated in survival analysis. Putative prognostic biomarkers were evaluated by Cox regression and consensus clustering. Mendelian randomization analyses assessed the causal role of PRKCB. Results: High DRS score associated with poor survival. Four genes (CDKN2A, NUP155, PPP2R2A, PRKCB) displayed prognostic utility. Three molecular subtypes were identified with discrete survival and immune properties. A 12-gene signature displayed robust prognostic performance. PRKCB was overexpressed in ESCC, while PRKCB knockdown reduced ESCC cell migration. Conclusions: This integrated single-cell sequencing analysis provides new insights into the molecular heterogeneity and prognostic determinants underlying ESCC. The findings identify potential prognostic biomarkers and a gene expression signature that may enable improved patient risk stratification in ESCC. Experimental validation of the role of PRKCB substantiates the potential clinical utility of our results.