# Characterising Epigenetic Tipping Points using a Spectral Dimension Reduction Approach

**Authors:** Tomás Alarcón, Javier A Menendez, Josep Sardanyés

PMC · DOI: 10.1007/s11538-026-01602-w · 2026-03-18

## TL;DR

This paper introduces a mathematical framework to predict and analyze tipping points in epigenetic landscapes, which could help understand and prevent cell identity loss in aging and cancer.

## Contribution

A novel spectral dimension reduction method is developed to model and predict epigenetic tipping points under chromatin-modifying enzyme competition.

## Key findings

- A dimension reduction approach accurately predicts global transitions in epigenetic landscapes.
- Metabolic cofactors SAM and acetyl-CoA are identified as potential early warning signals for epigenetic tipping points.
- The framework reveals how chromatin connectivity patterns influence the robustness of epigenetic landscapes.

## Abstract

Epigenetic landscapes (ELs) are defined by the pattern of epigenetic marks (acetylation, methylation, etc.) layed over large chromatin regions. The information contained in the ELs is essential to sustain the patterns of gene expression that shape cell fate and identity. EL maintenance requires the precise regulation of chromatin-modifying enzymes (ChME) and their metabolic cofactors (McF). Competition for ChME or dysregulation of McF abundance can lead to degradation of ELs, triggering large-scale changes in the cell fate information contained in EL. Thus, predicting impending epigenetic tipping points (ETPs) by identifying early warning signals (EWS) may help to anticipate the onset of cell identity loss during aging and cancer. Since ELs are formed (and maintained) by a systems of writer/eraser enzymes that interact both in cis (local) and trans (long-range) modes, their mathematical description involves a high-dimensional dynamical system, where identifying ETPs and characterising the biological mechanisms that control them remains challenging. Here, we develop a general mathematical framework that incorporates different connectivity patterns generated by the 3D chromatin folding structure to analyze competition-induced ETP in large EL. This framework allows us to measure the sensitivity and robustness of ETP to the availability of metabolic cofactors and to identify potential EWS. Using a dimension reduction method, we derived coarse-grained (CG) equations for the collective observables associated with chromatin modifications. Analysis of the CG system allows the prediction of global transitions that shape the large-scale features of EL, accurately reproduce the corresponding microscopic benchmarks, and reveal the existence of tipping points under conditions of ChME competition. We applied the CG method to predict ETP under different connectivity patterns, including heterogeneous profiles such as those found in Hi-C data. Although a robustness measure for stable EL was derived from the CG dynamics in bistable regimes, sensitivity analysis revealed that metabolic cofactors have the greatest impact on EL robustness. In particular, we identified the metabolic cofactors SAM and acetyl-CoA as potential EWS for the catastrophic loss of hyperacetylated EL induced by ChME competition. The ability to predict global ETP can facilitate the discovery of predictive biomarkers and inform metabolic interventions aimed at limiting and reversing pathological cell fate decisions.

The online version contains supplementary material available at 10.1007/s11538-026-01602-w.

## Linked entities

- **Chemicals:** SAM (PubChem CID 34755), acetyl-CoA (PubChem CID 444493)
- **Diseases:** cancer (MONDO:0004992)

## Full-text entities

- **Genes:** Esr1 (estrogen receptor 1 (alpha)) [NCBI Gene 13982] {aka ER, ER-alpha, ERa, ERalpha, ESR, Estr}, Ezh2 (enhancer of zeste 2 polycomb repressive complex 2 subunit) [NCBI Gene 14056] {aka Enx-1, Enx1h, KMT6, mKIAA4065}, Kdm6a (lysine (K)-specific demethylase 6A) [NCBI Gene 22289] {aka Utx}, Cbx2 (chromobox 2) [NCBI Gene 12416] {aka M33, MOD2, pc}, H49 (histocompatibility 49) [NCBI Gene 109816] {aka H(a<t>)}, Ewsr1 (Ewing sarcoma breakpoint region 1) [NCBI Gene 14030] {aka Ews, Ewsh}
- **Diseases:** breast cancer (MESH:D001943), influenza (MESH:D007251), AIDS (MESH:D000163), muscle injury (MESH:D009135), cystic fibrosis (MESH:D003550), cancer (MESH:D009369), fibrosis (MESH:D005355)
- **Chemicals:** Acetyl-CoA (MESH:D000105), S-Adenosyl methionine (MESH:D012436), acetate (MESH:D000085), Serine (MESH:D012694), EL (-), HM (MESH:C100283), methionine (MESH:D008715), carbon (MESH:D002244)
- **Species:** Mus musculus (house mouse, species) [taxon 10090]

## Figures

7 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12999790/full.md

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