Identification of colon cancer subtypes based on multi-omics data—construction of methylation markers for immunotherapy

Background: Being the most widely used biomarker for immunotherapy, the microsatellite status has limitations in identifying all patients who benefit in clinical practice. It is essential to identify additional biomarkers to guide immunotherapy. Aberrant DNA methylation is consistently associated with changes in the anti-tumor immune response, which can promote tumor progression. This study aims to explore immunotherapy biomarkers for colon cancers from the perspective of DNA methylation. Methods: The related data (RNA sequencing data and DNA methylation data) were obtained from The Cancer Genome Atlas (TCGA) and UCSC XENA database. Methylation-driven genes (MDGs) were identified through the Pearson correlation analysis. Unsupervised consensus clustering was conducted using these MDGs to identify distinct clusters of colon cancers. Subsequently, we evaluated the immune status and predicted the efficacy of immunotherapy by tumor immune dysfunction and exclusion (Tide) score. Finally, The Quantitative Differentially Methylated Regions (QDMR) software was used to identify the specific DNA methylation markers within particular clusters. Results: A total of 282 MDGs were identified by integrating the DNA methylation and RNA-seq data. Consensus clustering using the K-means algorithm revealed that the optimal number of clusters was 4. It was revealed that the composition of the tumor immune microenvironment (TIME) in Cluster 1 was significantly different from others, and it exhibited a higher level of tumor mutation burdens (TMB) and stronger anti-tumor immune activity. Furthermore, we identified three specific hypermethylation genes that defined Cluster 1 (PCDH20, APCDD1, COCH). Receiver operating characteristic (ROC) curves demonstrated that these specific markers could effectively distinguish Cluster 1 from other clusters, with an AUC of 0.947 (95% CI 0.903-0.990). Finally, we selected clinical samples for immunohistochemical validation. Conclusion: In conclusion, through the analysis of DNA methylation, consensus clustering of colon cancer could effectively identify the cluster that benefit from immunotherapy along with specific methylation biomarkers.