# Machine learning-based comprehensive analysis of m6A RNA methylation regulators in colorectal cancer: implications for prognosis, immune microenvironment, and immunotherapy response

**Authors:** Feifei Kong, Jiawei Feng, Haixia Shan, Youlong Zhu, Ling-Jun Zhu

PMC · DOI: 10.3389/ebm.2025.10776 · Experimental Biology and Medicine · 2026-01-14

## TL;DR

This study uses machine learning to analyze m6A RNA methylation regulators in colorectal cancer, showing they can predict patient survival and immunotherapy response.

## Contribution

A novel eight-gene m6A regulator signature was developed for CRC prognosis and immunotherapy prediction using machine learning.

## Key findings

- Random Forest achieved best performance with AUC of 0.895 in training and 0.847 in validation.
- Low-risk patients had higher CD8+ T cell infiltration and better predicted immunotherapy response.
- IGF2BP2 and METTL3 were top contributors to risk prediction based on SHAP analysis.

## Abstract

N6-methyladenosine (m6A) RNA methylation regulators have been implicated in colorectal cancer (CRC) progression. However, systematic evaluation using multiple machine learning approaches for prognostic prediction remains limited. This study aimed to develop and validate machine learning models for CRC prognosis based on m6A regulators and assess their potential for immunotherapy response prediction. We analyzed 1,047 CRC patients from TCGA and GEO databases (70% training, 30% validation). Twenty machine learning algorithms were systematically evaluated, with LASSO regression selecting optimal features from 27 m6A regulators. SHAP analysis provided model interpretability. Immune microenvironment characterization and immunotherapy response prediction were performed using established computational methods. LASSO regression selected eight m6A regulators (IGF2BP2, METTL3, HNRNPA2B1, METTL14, YTHDF2, VIRMA, FTO, ALKBH5) for model construction. Among 20 algorithms tested, Random Forest achieved optimal performance (training AUC = 0.895, validation AUC = 0.847). SHAP analysis identified IGF2BP2 (mean |SHAP| = 0.42) and METTL3 (mean |SHAP| = 0.36) as primary contributors to risk prediction. Risk stratification showed significant survival differences (HR = 2.41, 95% CI: 1.73–3.36, p < 0.001). Low-risk patients demonstrated enhanced immune infiltration with higher CD8+ T cells (17.8% vs. 10.2%, p < 0.001) and better predicted immunotherapy response rates (36.5% vs. 20.3%, p = 0.006). Our systematic machine learning analysis demonstrates that m6A regulators can effectively predict CRC prognosis and immunotherapy response. The eight-gene signature provides a practical tool for clinical risk assessment and treatment decision-making.

## Linked entities

- **Genes:** IGF2BP2 (insulin like growth factor 2 mRNA binding protein 2) [NCBI Gene 10644], METTL3 (methyltransferase 3, N6-adenosine-methyltransferase complex catalytic subunit) [NCBI Gene 56339], HNRNPA2B1 (heterogeneous nuclear ribonucleoprotein A2/B1) [NCBI Gene 3181], METTL14 (methyltransferase 14, N6-adenosine-methyltransferase non-catalytic subunit) [NCBI Gene 57721], YTHDF2 (YTH N6-methyladenosine RNA binding protein F2) [NCBI Gene 51441], VIRMA (vir like m6A methyltransferase associated) [NCBI Gene 25962], FTO (FTO alpha-ketoglutarate dependent dioxygenase) [NCBI Gene 79068], ALKBH5 (alkB homolog 5, RNA demethylase) [NCBI Gene 54890]
- **Diseases:** colorectal cancer (MONDO:0005575)

## Full-text entities

- **Genes:** CD8A (CD8 subunit alpha) [NCBI Gene 925] {aka CD8, CD8alpha, IMD116, Leu2, p32}, YTHDF2 (YTH N6-methyladenosine RNA binding protein F2) [NCBI Gene 51441] {aka CAHL, DF2, HGRG8, NY-REN-2}, HNRNPA2B1 (heterogeneous nuclear ribonucleoprotein A2/B1) [NCBI Gene 3181] {aka HNRNPA2, HNRNPB1, HNRPA2, HNRPA2B1, HNRPB1, IBMPFD2}, IGF2BP2 (insulin like growth factor 2 mRNA binding protein 2) [NCBI Gene 10644] {aka IMP-2, IMP2, VICKZ2}, FTO (FTO alpha-ketoglutarate dependent dioxygenase) [NCBI Gene 79068] {aka ALKBH9, BMIQ14, GDFD, IFEX9}, METTL3 (methyltransferase 3, N6-adenosine-methyltransferase complex catalytic subunit) [NCBI Gene 56339] {aka IME4, M6A, MT-A70, Spo8, hMETTL3}, METTL14 (methyltransferase 14, N6-adenosine-methyltransferase non-catalytic subunit) [NCBI Gene 57721] {aka hMETTL14}, ALKBH5 (alkB homolog 5, RNA demethylase) [NCBI Gene 54890] {aka ABH5, OFOXD, OFOXD1}, ITIH2 (inter-alpha-trypsin inhibitor heavy chain 2) [NCBI Gene 3698] {aka H2P, ITI-HC2, SHAP}
- **Diseases:** CRC (MESH:D015179)
- **Chemicals:** N6-methyladenosine (MESH:C010223), m6A (MESH:C005955)
- **Species:** Homo sapiens (human, species) [taxon 9606]

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12847061/full.md

## Figures

9 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12847061/full.md

## References

39 references — full list in the complete paper: https://tomesphere.com/paper/PMC12847061/full.md

---
Source: https://tomesphere.com/paper/PMC12847061