ABSTRACT Objectives To screen programmed cell death (PCD)-related genes in esophageal squamous cell carcinoma (ESCC) based on transcriptomic data and to explore its clinical value. Methods Differentially expressed PCD genes (DEPCDGs) were screened from ESCC transcriptome and clinical data in TCGA database. Univariate COX and LASSO COX were performed on prognostically DEPCDGs in ESCC to develop prognostic model. Differences in immune cell infiltration in different RiskScore groups were determined by ssGSEA and CIBERSORT. The role of RiskScore in immunotherapy response was explored using Tumor Immune Dysfunction and Exclusion (TIDE) and IMvigor210 cohorts. Results Fourteen DEPCDGs associated with prognosis were tapped in ESCC. These DEPCDGs form a RiskScore with good predictive performance for prognosis. RiskScore demonstrated excellent prediction accuracy in three data sets. The abundance of M2 macrophages and Tregs was higher in the high RiskScore group, and the abundance of M1 macrophages was higher in the low RiskScore group. The RiskScore also showed good immunotherapy sensitivity. RT-qPCR analysis showed that AUP1, BCAP31, DYRK2, TAF9 and UBQLN2 were higher expression in KYSE-150 cells. Knockdown BCAP31 inhibited migration and invasion. Conclusion A prognostic risk model can predict prognosis of ESCC and may be a useful biomarker for risk stratification and immunotherapy assessment.