Abstract
Background: Enhancer RNAs (eRNAs), the transcriptional products of active enhancers, are of great significance in the initial progression of cancers. However, the biological function and bioinformatics profiles of eRNA in gastric cancer remains largely enigmatic. Methods: Firstly, STAD were clustered into three subtypes with the data of eRNA expression from TCeA. Then we explored the difference of the tumor immune microenvironment, transcription levels, and transcription regulation among the three clusters. Finally, samples collected from 12 patients diagnosed with STAD were used to conduct qRT-PCR, verifying the conclusion based on network database. Results: The three clusters were detected to have different tumor microenvironments: Cluster A has an immune "cold" microenvironment. While cluster B features as more infiltration of immune cells, accompanied with higher expression of immune checkpoints such as PDCD1, LAG3, and TIGIT. Besides, Cluster C shows a higher stromal feature with B lineage, neutrophils, and fibroblasts. Further analyses indicated that CpG island methylation level of Cluster B is different from the other two clusters. Meanwhile, Cluster A and B showed significant enrichment of TP53 and KRAS mutation respectively while Cluster C has higher tumor mutation burden (TMB) and microsatellite instability (MSI). With the elaboration of transcriptional regulation of epigenetic clustering, we detected that Cluster A enriched in epithelial phenotype pathways. Cluster B enriched in cell-cell adhesion. Cluster C enriched in fibroblast proliferation. The clinical cohort show that Cluster B patients have lower interstitial cell characteristics and CAF infiltration. Conclusion: We identified three unique epigenetic clusters of STAD through the differential activation of super-enhancers, and identified Cluster B with a higher immune infiltrating and a better prognosis, which provides a novel understanding of eRNAs and potential clinical applicability of eRNA-based molecular subtypes in gastric cancer.