Abstract
N6-methyladenosine (m6A) and 5-methylcytosine (m5C) RNA modifications have garnered significant attention in the field of epigenetic research due to their close association with human cancers. This study we focus on elucidating the expression patterns of m6A/m5C-related long non-coding RNAs (lncRNAs) in esophageal squamous cell carcinoma (ESCC) and assessing their prognostic significance and therapeutic potential. Transcriptomic profiles of ESCC were derived from public resources. m6A/m5C-related lncRNAs were obtained from TCGA using Spearman's correlations analysis. The m6A/m5C-lncRNAs prognostic signature was selected to construct a RiskScore model for survival prediction, and their correlation with the immune microenvironment and immunotherapy response was analyzed. A total of 606 m6A/m5C-lncRNAs were screened, and ESCC cases in the TCGA cohort were stratified into three clusters, which showed significantly distinct in various clinical features and immune landscapes. A RiskScore model comprising ten m6A/m5C-lncRNAs prognostic signature were constructed and displayed good independent prediction ability in validation datasets. Patients in the low-RiskScore group had a better prognosis, a higher abundance of immune cells (CD4 + T cell, CD4 + naive T cell, class-switched memory B cell, and Treg), and enhanced expression of most immune checkpoint genes. Importantly, patients with low-RiskScore were more cline benefit from immune checkpoint inhibitor treatment (P < 0.05). Our findings underscore the potential of RiskScore system comprising ten m6A/m5C-related lncRNAs as effective biomarkers for predicting survival outcomes, characterizing the immune landscape, and assessing response to immunotherapy in ESCC.