Characterizing m6A modification factors and their interactions in colorectal cancer: implications for tumor subtypes and clinical outcomes.

Background: The study aims to comprehensively combine colorectal cancer data cohorts in order to analyze the effects of various DNA methylation-coding genes on colorectal cancer patients. The annual incidence and mortality of colorectal cancer are very high, and there are no effective treatments for advanced colorectal cancer. DNA methylation is a method widely used to regulate epigenetics in the molecular mechanism study of tumors.

Method: Three single-cell cohorts GSE166555, GSE146771, and EMTAB8107, and five transcriptome cohorts GSE17536, GSE39582, GSE72970, and TCGA-CRC (TCGA-COAD and TCGA-READ) were applied in this study. 2 erasers (ALKBH5 and FTO), There are 7 writers (METTL3, METTL14, WTAP, VIRMA, RBM15, RBM15B, and ZC3H13) and 11 readers (YTHDC1, IGF2BP1, IGF2BP2, IGF2BP3, YTHDF1, YTHDF3, YTHDC2, and HNRNPA2B1, YTHDF2, HNRNPC and RBMX), a total of 20 M6A regulators, were used as the basis of the dataset in this study and were applied to the construction of molecular typing and prognostic models. Drugs that are differentially sensitive in methylation-regulated gene-related prognostic models were identified using the ConsensusClusterPlus package, which was also used to identify distinct methylation regulatory expression patterns in colorectal cancer and to model the relationship between tissue gene expression profiles and drug IC50 values. Finally, TISCH2 assessed which immune cells were significantly expressed with M6A scores. The immunosuppression of M6A methylation is spatially explained.

Results: This study used data from 583 CRC patients in the TCGA-CRC cohort. Firstly, the mutation frequency and CNV variation frequency of 20 m6A modification-related factors were analyzed, and the corresponding histogram and heat map were drawn. The study next analyzed the expression variations between mutant and wild forms of the VIRMA gene and explored differences in the expression of these variables in tumor and normal tissues. In addition, the samples were divided into different subgroups by molecular clustering method based on m6A modification, and each subgroup's expression and clinicopathological characteristics were analyzed. Finally, we compared prognostic differences, tumor microenvironment (TME) characteristics, immune cell infiltration, and gene function enrichment among different subpopulations. We also developed a colorectal cancer m6A-associated gene signature and validated its prognostic effects across multiple cohorts. Finally, using single-cell RNA sequencing data, we confirmed that tumor cells show elevated expression of m6A-related gene signatures.

Discussion: This study explored the mutation frequency, expression differences, interactions, molecular clustering, prognostic effect, and association with tumor characteristics of m6A modification-related factors in CRC and validated them at the single-cell level. These results clarify the association between m6A alteration and colorectal cancer (CRC) and offer important insights into the molecular recognition and management of cancer.