# Identification of immune-mediated aging genes associated with cervical spondylosis through single-cell eQTL Mendelian randomization

**Authors:** Zuozhong Liu, Haitao Xu

PMC · DOI: 10.3389/fimmu.2026.1766215 · Frontiers in Immunology · 2026-02-18

## TL;DR

This study identifies immune-related genes linked to cervical spondylosis and suggests 5-ASA as a potential treatment by analyzing genetic and immune data.

## Contribution

A novel single-cell eQTL Mendelian randomization framework reveals immune aging genes and repurposable drugs for cervical spondylosis.

## Key findings

- 118 genes, especially BMPR2, CHUK, CTNNB1, CTSB, and EZH2, show causal links to cervical spondylosis in CD4+ and CD8+ T cells.
- 5-ASA treatment modulates gene expression and cytokine levels in HUT78 T cells, suggesting therapeutic potential.
- Age-related gene expression changes and cytokine alterations were validated in patient serum samples.

## Abstract

Cervical spondylosis (CS) is an age-related degenerative spinal disorder with increasing evidence linking its pathogenesis to immune aging. The genetic mechanisms connecting immune dysregulation to CS remain largely undefined.

We integrated GWAS data from 484,598 individuals with single-cell eQTL profiles across 14 immune cell types, using Mendelian randomization (MR) and Bayesian colocalization to identify genes whose cell-type–specific expression causally influences CS risk. Entropy analysis quantified lineage specificity of MR effects. Functional enrichment, drug–gene interaction mapping, and protein–protein interaction (PPI) network construction were used for target prioritization. Validation was performed using qPCR and ELISA in CS patient serum and 5-Aminosalicylic Acid (5-ASA) treatment in HUT78 T cells.

A total of 118 genes exhibited suggestive or potential causal associations with CS, primarily enriched in CD4+ and CD8+ T cells. Five immune aging-related genes BMPR2, CHUK, CTNNB1, CTSB, and EZH2 emerged as central regulators based on entropy scores and colocalization evidence, PPI centrality, and druggability. These genes were validated in patient serum samples, showing age-associated expression changes and cytokine alterations. 5-ASA treatment modulated their expression and inflammatory cytokine levels in vitro, supporting its repurposing potential.

This study reveals immune cell–specific genetic regulators linking immune aging to cervical spondylosis and identifies 5-ASA as a candidate therapeutic agent. Our single-cell MR framework offers novel insights into immunogenetic mechanisms driving age-related spinal degeneration and highlights actionable targets for translational research.

## Linked entities

- **Genes:** BMPR2 (bone morphogenetic protein receptor type 2) [NCBI Gene 659], CHUK (component of inhibitor of nuclear factor kappa B kinase complex) [NCBI Gene 1147], CTNNB1 (catenin beta 1) [NCBI Gene 1499], CTSB (cathepsin B) [NCBI Gene 1508], EZH2 (enhancer of zeste 2 polycomb repressive complex 2 subunit) [NCBI Gene 2146]
- **Chemicals:** 5-Aminosalicylic Acid (PubChem CID 4075), 5-ASA (PubChem CID 4075)

## Full-text entities

- **Genes:** IL1B (interleukin 1 beta) [NCBI Gene 3553] {aka IL-1, IL1-BETA, IL1F2, IL1beta}, CTNNB1 (catenin beta 1) [NCBI Gene 1499] {aka CTNNB, EVR7, MRD19, NEDSDV, armadillo}, PYCARD (PYD and CARD domain containing) [NCBI Gene 29108] {aka ASC, CARD5, TMS, TMS-1, TMS1}, CEACAM1 (CEA cell adhesion molecule 1) [NCBI Gene 634] {aka BGP, BGP1, BGPI}, CHUK (component of inhibitor of nuclear factor kappa B kinase complex) [NCBI Gene 1147] {aka BPS2, IKBKA, IKK-1, IKK-alpha, IKK1, IKKA}, CD8A (CD8 subunit alpha) [NCBI Gene 925] {aka CD8, CD8alpha, IMD116, Leu2, p32}, CASP1 (caspase 1) [NCBI Gene 834] {aka ICE, IL1BC, P45}, TGFB1 (transforming growth factor beta 1) [NCBI Gene 7040] {aka CAEND1, CED, DPD1, IBDIMDE, LAP, TGF-beta1}, IRAK4 (interleukin 1 receptor associated kinase 4) [NCBI Gene 51135] {aka IMD67, IPD1, IRAK-4, NY-REN-64, REN64}, AGO2 (argonaute RISC catalytic component 2) [NCBI Gene 27161] {aka CASC7, EIF2C2, LESKRES, LINC00980, PPD, Q10}, FLT1 (fms related receptor tyrosine kinase 1) [NCBI Gene 2321] {aka FLT, FLT-1, VEGFR-1, VEGFR1}, SIRT2 (sirtuin 2) [NCBI Gene 22933] {aka SIR2, SIR2L, SIR2L2}, NLRP3 (NLR family pyrin domain containing 3) [NCBI Gene 114548] {aka AGTAVPRL, AII, AVP, C1orf7, CIAS1, CLR1.1}, TNF (tumor necrosis factor) [NCBI Gene 7124] {aka DIF, IMD127, TNF-alpha, TNFA, TNFSF2, TNLG1F}, MAPK1 (mitogen-activated protein kinase 1) [NCBI Gene 5594] {aka ERK, ERK-2, ERK2, ERT1, MAPK2, NS13}, MAP2K1 (mitogen-activated protein kinase kinase 1) [NCBI Gene 5604] {aka CFC3, MAPKK1, MEK1, MEL, MKK1, PRKMK1}, TGFBR2 (transforming growth factor beta receptor 2) [NCBI Gene 7048] {aka AAT3, FAA3, LDS1B, LDS2, LDS2B, MFS2}, BMPR2 (bone morphogenetic protein receptor type 2) [NCBI Gene 659] {aka BMPR-II, BMPR3, BMR2, BRK-3, POVD1, PPH1}, CTSB (cathepsin B) [NCBI Gene 1508] {aka APPS, CPSB, KWE, RECEUP}, ATM (ATM serine/threonine kinase) [NCBI Gene 472] {aka AT1, ATA, ATC, ATD, ATDC, ATE}, ATR (ATR checkpoint kinase) [NCBI Gene 545] {aka FCTCS, FRP1, MEC1, SCKL, SCKL1}, CDK6 (cyclin dependent kinase 6) [NCBI Gene 1021] {aka MCPH12, PLSTIRE}, CDK1 (cyclin dependent kinase 1) [NCBI Gene 983] {aka CDC2, CDC28A, P34CDC2}, IL18 (interleukin 18) [NCBI Gene 3606] {aka IGIF, IL-18, IL-1g, IL1F4}, MAP2K2 (mitogen-activated protein kinase kinase 2) [NCBI Gene 5605] {aka CFC4, MAPKK2, MEK2, MKK2, PRKMK2}, SOD2 (superoxide dismutase 2) [NCBI Gene 6648] {aka GC1, GClnc1, IPO-B, IPOB, MNSOD, MVCD6}, TP53 (tumor protein p53) [NCBI Gene 7157] {aka BCC7, BMFS5, LFS1, P53, TRP53}, BMP1 (bone morphogenetic protein 1) [NCBI Gene 649] {aka OI13, PCOLC, PCP, TLD}, JAK2 (Janus kinase 2) [NCBI Gene 3717] {aka JTK10}, CS (citrate synthase) [NCBI Gene 1431], PIK3CA (phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha) [NCBI Gene 5290] {aka CCM4, CLAPO, CLOVE, CWS5, HMH, MCAP}, NFKB1 (nuclear factor kappa B subunit 1) [NCBI Gene 4790] {aka CVID12, EBP-1, KBF1, NF-kB, NF-kB1, NF-kappa-B1}, CD4 (CD4 molecule) [NCBI Gene 920] {aka CD4mut, IMD79, Leu-3, OKT4D, T4}, EZH2 (enhancer of zeste 2 polycomb repressive complex 2 subunit) [NCBI Gene 2146] {aka ENX-1, ENX1, EZH2b, KMT6, KMT6A, WVS}, GAPDH (glyceraldehyde-3-phosphate dehydrogenase) [NCBI Gene 2597] {aka G3PD, GAPD, HEL-S-162eP}
- **Diseases:** CS (MESH:D055009), numbness and (MESH:D006987), MR (MESH:C562757), age-related degeneration (MESH:D008268), age-related diseases (MESH:D010024), dizziness (MESH:D004244), age (MESH:D019588), tinnitus (MESH:D014012), neck, shoulder, and back pain (MESH:D020069), gastrointestinal dysfunction (MESH:D005767), diabetes (MESH:D003920), tumor (MESH:D009369), hypertension (MESH:D006973), blurred vision (MESH:D014786), osteophyte (MESH:D054850), spine (MESH:D016135), chronic degenerative diseases (MESH:D019636), chronic inflammation (MESH:D007249), headache (MESH:D006261), T-cell lymphoma (MESH:D016399), fibrosis (MESH:D005355), ulcerative colitis (MESH:D003093), chronic (MESH:D002908), cervical spine degeneration (MESH:D002575), vomiting (MESH:D014839), immune dysregulation (OMIM:614878), spinal degeneration (MESH:D009410), musculoskeletal, metabolic, or autoimmune disorders (MESH:D009140), nausea (MESH:D009325), inflammatory bowel diseases (MESH:D015212), intervertebral disc degeneration (MESH:D055959)
- **Chemicals:** DB00233 (-), 5-ASA (MESH:D019804)
- **Species:** Homo sapiens (human, species) [taxon 9606]
- **Cell lines:** HUT78 — Homo sapiens (Human), Sezary syndrome, Cancer cell line (CVCL_0337)

## Full text

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

## Figures

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

## References

46 references — full list in the complete paper: https://tomesphere.com/paper/PMC12956539/full.md

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