Abstract
It remains undetermined regarding the impact of neoadjuvant therapy on immunogenic cell death (ICD) and subsequent tumor microenvironment (TME) remodeling in esophageal squamous cell carcinoma (ESCC). And it is of paramount significance to identify beneficiaries from neoadjuvant therapy in treatment-naïve ESCC. In this study, 88 ESCC samples undergoing neoadjuvant therapy plus surgery (NA+S) or surgery alone (SA) were subjected to bulk-RNA sequencing. A five-gene RINscore incorporating ICD-related signature genes with TME-based hub genes was established to predict clinical outcomes and pharmacological responses, in which SLAMF7 and IL1R1 were selected out as co-expressed genes. The regulatory mechanism of the repressive co-transcription factor BATF of SLAMF7 and IL1R1 was further demonstrated. Our data demonstrated that NA+S led to high abundance in kinds of T helper cells, nature killer T cells and M1-like macrophages with increased CD8+T cells infiltration compared with SA. ICD phenotypes were further characterized in treatment-naïve ESCC to determine their differences in TME and potential benefits from NA. Our findings not only offered novel insights into the distinct TME and ICD profiles of ESCC undergoing different therapeutic modes, but also provided the RINscore, which may aid oncologists in determining individualized (neo)adjuvant immunotherapy regimen.