# Integrated Single-Cell and Spatial Transcriptomics Coupled with Machine Learning Uncovers MORF4L1 as a Critical Epigenetic Mediator of Radiotherapy Resistance in Colorectal Cancer Liver Metastasis

**Authors:** Yuanyuan Zhang, Xiaoli Wang, Haitao Liu, Yan Xiang, Le Yu

PMC · DOI: 10.3390/biomedicines14020273 · Biomedicines · 2026-01-26

## TL;DR

This study identifies MORF4L1 as a key factor in radiotherapy resistance in colorectal cancer liver metastasis using advanced genomic and machine learning techniques.

## Contribution

The study introduces MORF4L1 as a novel epigenetic mediator of radioresistance in CRLM, supported by integrated single-cell and spatial transcriptomics.

## Key findings

- MORF4L1-high tumor cells show enhanced DNA repair and metabolic pathways linked to radioresistance.
- MORF4L1-expressing tumor regions spatially co-localize with immunosuppressive immune cells like Tregs and TAMs.
- Radiotherapy drives tumor cells toward a resistant state marked by increased MORF4L1 expression.

## Abstract

Background and Objective: Colorectal cancer (CRC) liver metastasis (CRLM) represents a major clinical challenge, and acquired resistance to radiotherapy (RT) significantly limits therapeutic efficacy. A deep and comprehensive understanding of the cellular and molecular mechanisms driving RT resistance is urgently required to develop effective combination strategies. Here, we aimed to dissect the dynamic cellular landscape of the tumor microenvironment (TME) and identify key epigenetic regulators mediating radioresistance in CRLM by integrating cutting-edge single-cell and spatial omics technologies. Methods and Results: We performed integrated single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) on matched pre- and post-radiotherapy tumor tissues collected from three distinct CRLM patients. Employing a robust machine-learning framework on the multi-omics data, we successfully identified MORF4L1 (Mortality Factor 4 Like 1), an epigenetic reader, as a critical epigenetic mediator of acquired radioresistance. High-resolution scRNA-seq analysis of the tumor cell compartment revealed that the MORF4L1-high subpopulation exhibited significant enrichment in DNA damage repair (DDR) pathways, heightened activity of multiple pro-survival metabolic pathways, and robust signatures of immune evasion. Pseudotime trajectory analysis further confirmed that RT exposure drives tumor cells toward a highly resistant state, marked by a distinct increase in MORF4L1 expression. Furthermore, cell–cell communication inference demonstrated a pronounced, systemic upregulation of various immunosuppressive signaling axes within the TME following RT. Crucially, high-resolution ST confirmed these molecular and cellular interactions in their native context, revealing a significant spatial co-localization of MORF4L1-expressing tumor foci with multiple immunosuppressive immune cell types, including regulatory T cells (Tregs) and tumor-associated macrophages (TAMs), thereby underscoring its role in TME-mediated resistance. Conclusions: Our comprehensive spatial and single-cell profiling establishes MORF4L1 as a pivotal epigenetic regulator underlying acquired radioresistance in CRLM. These findings provide a compelling mechanistic rationale for combining radiotherapy with the targeted inhibition of MORF4L1, presenting a promising new therapeutic avenue to overcome treatment failure and improve patient outcomes in CRLM.

## Linked entities

- **Genes:** MORF4L1 (mortality factor 4 like 1) [NCBI Gene 10933]
- **Diseases:** colorectal cancer (MONDO:0005575)

## Full-text entities

- **Genes:** TGFB1 (transforming growth factor beta 1) [NCBI Gene 7040] {aka CAEND1, CED, DPD1, IBDIMDE, LAP, TGF-beta1}, HLA-C (major histocompatibility complex, class I, C) [NCBI Gene 3107] {aka D6S204, HLA-JY3, HLAC, HLC-C, MHC, PSORS1}, ARID1B (AT-rich interaction domain 1B) [NCBI Gene 57492] {aka 6A3-5, BAF250B, BRIGHT, CSS1, DAN15, ELD/OSA1}, EGR1 (early growth response 1) [NCBI Gene 1958] {aka AT225, G0S30, KROX-24, NGFI-A, TIS8, ZIF-268}, PRKCA (protein kinase C alpha) [NCBI Gene 5578] {aka AAG6, PKC-alpha, PKCA, PKCI+/-, PKCalpha}, SPP1 (secreted phosphoprotein 1) [NCBI Gene 6696] {aka BNSP, BSPI, ETA-1, OPN}, CUX1 (cut like homeobox 1) [NCBI Gene 1523] {aka CASP, CDP, CDP/Cut, CDP1, COY1, CUTL1}, KRAS (KRAS proto-oncogene, GTPase) [NCBI Gene 3845] {aka 'C-K-RAS, C-K-RAS, CFC2, K-RAS2A, K-RAS2B, K-RAS4A}, ZMYND8 (zinc finger MYND-type containing 8) [NCBI Gene 23613] {aka PRKCBP1, PRO2893, RACK7}, STAT5A (signal transducer and activator of transcription 5A) [NCBI Gene 6776] {aka MGF, STAT5}, MT1X (metallothionein 1X) [NCBI Gene 4501] {aka MT-1l, MT1}, IL6 (interleukin 6) [NCBI Gene 3569] {aka BSF-2, BSF2, CDF, HGF, HSF, IFN-beta-2}, Igf2 (insulin-like growth factor 2) [NCBI Gene 16002] {aka Igf-2, Igf-II, M6pr, Mpr, Peg2}, YY1 (YY1 transcription factor) [NCBI Gene 7528] {aka DELTA, GADEVS, INO80S, NF-E1, UCRBP, YIN-YANG-1}, MYC (MYC proto-oncogene, bHLH transcription factor) [NCBI Gene 4609] {aka MRTL, MYCC, bHLHe39, c-Myc}, ELF1 (E74 like ETS transcription factor 1) [NCBI Gene 1997] {aka EFTUD1, RIA1}, STAT3 (signal transducer and activator of transcription 3) [NCBI Gene 6774] {aka ADMIO, ADMIO1, APRF, HIES}, H4C6 (H4 clustered histone 6) [NCBI Gene 8361] {aka H4, H4/c, H4FC, HIST1H4F}, MKI67 (marker of proliferation Ki-67) [NCBI Gene 4288] {aka KIA, MIB-, MIB-1, PPP1R105}, HSPA1A (heat shock protein family A (Hsp70) member 1A) [NCBI Gene 3303] {aka HEL-S-103, HSP70, HSP70-1, HSP70-1A, HSP70-2, HSP70.1}, MORF4L1 (mortality factor 4 like 1) [NCBI Gene 10933] {aka Eaf3, FWP006, HsT17725, MEAF3, MORFRG15, MRG15}, NCOA2 (nuclear receptor coactivator 2) [NCBI Gene 10499] {aka GRIP1, KAT13C, NCoA-2, SRC-2, SRC2, TIF2}, ETV4 (ETS variant transcription factor 4) [NCBI Gene 2118] {aka E1A-F, E1AF, PEA3, PEAS3}, Morf4l1 (mortality factor 4 like 1) [NCBI Gene 21761] {aka MORFRG15, MRG15, TEG-189, Tex189, mKIAA4002}, H2AC18 (H2A clustered histone 18) [NCBI Gene 8337] {aka H2A, H2A.2, H2A/O, H2A/q, H2AFO, H2a-615}, TP53 (tumor protein p53) [NCBI Gene 7157] {aka BCC7, BMFS5, LFS1, P53, TRP53}, CCL4 (C-C motif chemokine ligand 4) [NCBI Gene 6351] {aka ACT2, AT744.1, G-26, HC21, LAG-1, LAG1}, KAT5 (lysine acetyltransferase 5) [NCBI Gene 10524] {aka ESA1, HTATIP, HTATIP1, NEDFASB, PLIP, TIP}, Egr1 (early growth response 1) [NCBI Gene 13653] {aka A530045N19Rik, ETR103, Egr-1, Krox-1, Krox-24, Krox24}, AREG (amphiregulin) [NCBI Gene 374] {aka AR, AREGB, CRDGF, SDGF}, H2ax (H2A.X variant histone) [NCBI Gene 15270] {aka H2A.X, H2afx, Hist5-2ax, gammaH2ax}, BRAF (B-Raf proto-oncogene, serine/threonine kinase) [NCBI Gene 673] {aka B-RAF1, B-raf, BRAF-1, BRAF1, NS7, RAFB1}, HLA-A (major histocompatibility complex, class I, A) [NCBI Gene 3105] {aka HLAA}, HLA-B (major histocompatibility complex, class I, B) [NCBI Gene 3106] {aka AS, B-4901, HLAB}, ORC1 (origin recognition complex subunit 1) [NCBI Gene 4998] {aka HSORC1, ORC1L, PARC1}, IGF2R (insulin like growth factor 2 receptor) [NCBI Gene 3482] {aka CD222, CI-M6PR, CIMPR, M6P-R, M6P/IGF2R, MPR 300}, BANF1 (barrier to autointegration nuclear assembly factor 1) [NCBI Gene 8815] {aka BAF, BCRP1, D14S1460, NGPS}, IGF2 (insulin like growth factor 2) [NCBI Gene 3481] {aka C11orf43, GRDF, IGF-II, PP9974, SRS3}, LAMTOR3 (late endosomal/lysosomal adaptor, MAPK and MTOR activator 3) [NCBI Gene 8649] {aka MAP2K1IP1, MAPBP, MAPKSP1, MP1, PRO0633, Ragulator3}, IL2 (interleukin 2) [NCBI Gene 3558] {aka IL-2, TCGF, lymphokine}
- **Diseases:** Cancer (MESH:D009369), MSI-H cancers (MESH:D053842), cytotoxic (MESH:D064420), injury to (MESH:D014947), inflammation (MESH:D007249), CRLM (MESH:D015179)
- **Chemicals:** oxalate (MESH:D010070), fatty acid (MESH:D005227), H &amp; E (MESH:D006371), Hematoxylin and Eosin (-), eicosanoid (MESH:D015777), ROS (MESH:D017382), glutathione (MESH:D005978), leukotriene D4 (MESH:D017998), Lipid (MESH:D008055)
- **Species:** Mus musculus (house mouse, species) [taxon 10090], 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/PMC12937616/full.md

## Figures

6 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12937616/full.md

## References

24 references — full list in the complete paper: https://tomesphere.com/paper/PMC12937616/full.md

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