Comprehensive profiling of enhancer RNA in stage II/III colorectal cancer defines two prognostic subtypes with implications for immunotherapy.

Abstract

Recently, enhancer RNAs (eRNAs) have garnered attention as pivotal biomarkers for the onset and progression of cancer. However, the landscape of eRNAs and the implications of eRNA-based molecular subtypes in stage II/III colorectal cancer (CRC) remain largely unexplored. Comprehensive profiling of eRNAs was conducted on a public stage II/III CRC cohort with total RNA-seq data. We used unsupervised clustering of prognostic eRNAs to establish an eRNA-based subtyping system. Further evaluations included molecular characteristics, immune infiltration, clinical outcomes, and drug responses. Finally, we validated the eRNA-based subtyping system in The Cancer Genome Atlas (TCGA) CRC cohort. We identified a total of 6453 expressed eRNAs, among which 237 were prognostic. A global upregulation of eRNAs was observed in microsatellite-stable (MSS) CRCs when compared to microsatellite instability-high (MSI-H) CRCs. Through consensus clustering, two novel molecular subtypes, termed Cluster 1(C1) and Cluster 2(C2), were further identified. C1, associated with the activation of epithelial-mesenchymal transition (EMT), hypoxia, and KRAS signaling pathways, showed poorer prognosis. C2, correlated with the canonical CRC subtype, exhibited superior survival outcomes. In addition, C1 showed enrichment with immune infiltration and more sensitivity to immune checkpoint inhibitors. Our study unravels the molecular heterogeneity of stage II/III CRC at the eRNA level and highlights the potential applications of the novel eRNA-based subtyping system in predicting prognosis and guiding immunotherapy.