Innovative prognostic modeling in ESCC: leveraging scRNA-seq and bulk-RNA for dendritic cell heterogeneity analysis

Globally, esophageal squamous cell carcinoma (ESCC) stands out as a common cancer type, characterized by its notably high rates of occurrence and mortality. Recent advancements in treatment methods, including immunotherapy, have shown promise, yet the prognosis remains poor. In the context of tumor development and treatment outcomes, the tumor microenvironment (TME), especially the function of dendritic cells (DCs), is significantly influential. Our study aims to delve deeper into the heterogeneity of DCs in ESCC using single-cell RNA sequencing (scRNA-seq) and bulk RNA analysis.

In the scRNA-seq analysis, we utilized the SCP package for result visualization and functional enrichment analysis of cell subpopulations. CellChat was employed to identify potential oncogenic mechanisms in DCs, while Monocle 2 traced the evolutionary trajectory of the three DC subtypes. CopyKAT assessed the benign or malignant nature of cells, and SCENIC conducted transcription factor regulatory network analysis, offering a preliminary exploration of DC heterogeneity. In Bulk-RNA analysis, we constructed a prognostic model for ESCC prognosis and immunotherapy response, based on DC marker genes. This model was validated through quantitative PCR (qPCR) and immunohistochemistry (IHC), confirming the gene expression levels.

In this study, through intercellular communication analysis, we identified GALECTIN and MHC-I signaling pathways as potential oncogenic mechanisms within dendritic cells. We categorized DCs into three subtypes: plasmacytoid (pDC), conventional (cDC), and tolerogenic (tDC). Our findings revealed that pDCs exhibited an increased proportion of cells in the G2/M and S phases, indicating enhanced cellular activity. Pseudotime trajectory analysis demonstrated that cDCs were in early stages of differentiation, whereas tDCs were in more advanced stages, with pDCs distributed across both early and late differentiation phases. Prognostic analysis highlighted a significant correlation between pDCs and tDCs with the prognosis of ESCC (P< 0.05), while no significant correlation was observed between cDCs and ESCC prognosis (P = 0.31). The analysis of cell malignancy showed the lowest proportion of malignant cells in cDCs (17%), followed by pDCs (29%), and the highest in tDCs (48%), with these results being statistically significant (P< 0.05). We developed a robust ESCC prognostic model based on marker genes of pDCs and tDCs in the GSE53624 cohort (n = 119), which was validated in the TCGA-ESCC cohort (n = 139) and the IMvigor210 immunotherapy cohort (n = 298) (P< 0.05). Additionally, we supplemented the study with a novel nomogram that integrates clinical features and risk assessments. Finally, the expression levels of genes involved in the model were validated using qPCR (n = 8) and IHC (n = 16), thereby confirming the accuracy of our analysis.

This study enhances the understanding of dendritic cell heterogeneity in ESCC and its impact on patient prognosis. The insights gained from scRNA-seq and Bulk-RNA analysis contribute to the development of novel biomarkers and therapeutic targets. Our prognostic models based on DC-related gene signatures hold promise for improving ESCC patient stratification and guiding treatment decisions.