# The Immune Cell Atlas of “Longevity Molecular Tag”: Identification of Principal Immune Cell Subsets and Their Underlying Molecular Regulatory Mechanisms

**Authors:** Zhiling Zhang, Huabin Su, Shihui Fu, Fansen Ji, Liuguan Liang, Wanlu Song, Caiyou Hu, Liuxiang Wei, Erping Long, Yang Lin, Xiaolin Ni

PMC · DOI: 10.1111/acel.70431 · 2026-03-05

## TL;DR

This study identifies immune cell subsets and molecular mechanisms linked to longevity in centenarians, offering insights into healthy aging.

## Contribution

The study introduces an immune cell atlas of 'Longevity Molecular Tags' using the Scissor algorithm to uncover novel immune mechanisms in centenarians.

## Key findings

- Scissor+ cells include NK cells, CD8+ T cells, and γδ T cells with enhanced cytotoxic functions.
- Scissor− cells involve CD4+ T cells, B cells, and DCs linked to inflammatory pathways.
- Five eQTL-colocalized events were identified, suggesting genetic links to longevity.

## Abstract

Immunosenescence represents a critical aspect of the aging process. Centenarians, serving as a nature model of “healthy aging,” demonstrate a distinctive immune “compensatory adaptation” mechanism that contributes to the maintenance of immune homeostasis. However, the specific immune cell subsets involved and the molecular mechanisms underlying these phenotypic traits remain incompletely understood. In this study, we integrated single‐cell RNA sequencing data spanning the entire lifespan of East Asian populations with bulk transcriptomic data from a centenarian cohort in Guangxi. Utilizing the Scissor algorithm, we identified immune cell subpopulations positively (Scissor+) and negatively (Scissor−) associated with longevity phenotypes, thereby constructing an immune cell atlas of “Longevity Molecular Tag.” Our findings indicate that Scissor+ cells predominantly comprise natural killer (NK) cells, CD8+ T cells, and γδ T cells, characterized by enhanced cytotoxic and immunomodulatory functions. Conversely, Scissor− cells mainly include CD4+ T cells, B cells, and dendritic cells (DCs), which are linked to inflammatory signaling pathways and Th17/Th1 differentiation. Trajectory analysis elucidated the differentiation pathways of NK, CD8+ T cells, CD4+ T cells, and B cells. Differentially expressed genes were enriched in pathways such as NF‐κB signaling, T cell receptor signaling, and NK cell cytotoxicity. Furthermore, co‐localization analysis revealed five eQTL‐colocalized events (rs3793537–GLIPR2/CD72/TLN1 and rs8019902–TRDV2/TRDC) associated with longevity. Collectively, these results suggest that centenarians achieve immune equilibrium by remodeling cytotoxic immune lineages and finely tuning inflammatory responses, thereby promoting health span and longevity. This study offers novel insights into potential strategies for modulating immunosenescence.

This study aims to develop an immune cell atlas of the “Longevity Molecular Tags” by employing the Scissor algorithm to identify key immune cell subsets associated with longevity phenotypes. Subsequently, the underlying immune molecular mechanisms contributing to the longevity‐associated immune cell subsets were investigated through differential gene expression analysis, pseudotime analysis, and co‐localization analysis.

## Linked entities

- **Genes:** GLIPR2 (GLI pathogenesis related 2) [NCBI Gene 152007], CD72 (CD72 molecule) [NCBI Gene 971], TLN1 (talin 1) [NCBI Gene 7094], TRDV2 (T cell receptor delta variable 2) [NCBI Gene 28517], TRDC (T cell receptor delta constant) [NCBI Gene 28526]

## Full-text entities

- **Genes:** 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}, MS4A1 (membrane spanning 4-domains A1) [NCBI Gene 931] {aka B1, Bp35, CD20, CVID5, FMC7, LEU-16}, Stat5a (signal transducer and activator of transcription 5A) [NCBI Gene 20850] {aka STAT5}, Ifng (interferon gamma) [NCBI Gene 15978] {aka IFN-g, If2f, Ifg}, TRDV2 (T cell receptor delta variable 2) [NCBI Gene 28517] {aka hDV102S1}, ZNF683 (zinc finger protein 683) [NCBI Gene 257101] {aka Hobit}, ANXA1 (annexin A1) [NCBI Gene 301] {aka ANX1, LPC1}, MFT2 (Trichoepithelioma, multiple familial, 2) [NCBI Gene 100188881] {aka TEM}, Il7r (interleukin 7 receptor) [NCBI Gene 16197] {aka CD127, IL-7Ralpha}, TRDC (T cell receptor delta constant) [NCBI Gene 28526] {aka TCRD}, CCL4 (C-C motif chemokine ligand 4) [NCBI Gene 6351] {aka ACT2, AT744.1, G-26, HC21, LAG-1, LAG1}, GZMB (granzyme B) [NCBI Gene 3002] {aka C11, CCPI, CGL-1, CGL1, CSP-B, CSPB}, SLC4A10 (solute carrier family 4 member 10) [NCBI Gene 57282] {aka NBCn2, NCBE}, LTB (lymphotoxin beta) [NCBI Gene 4050] {aka TNFC, TNFSF3, TNLG1C, p33}, CCL5 (C-C motif chemokine ligand 5) [NCBI Gene 6352] {aka D17S136E, RANTES, SCYA5, SIS-delta, SISd, TCP228}, Il7 (interleukin 7) [NCBI Gene 16196] {aka A630026I06Rik, Il-7, hlb368}, IGHM (immunoglobulin heavy constant mu) [NCBI Gene 3507] {aka AGM1, MU, VH}, CCR7 (C-C motif chemokine receptor 7) [NCBI Gene 1236] {aka BLR2, CC-CKR-7, CCR-7, CD197, CDw197, CMKBR7}, NKG7 (natural killer cell granule protein 7) [NCBI Gene 4818] {aka GIG1, GMP-17, p15-TIA-1}, Tnf (tumor necrosis factor) [NCBI Gene 21926] {aka DIF, TNF-a, TNF-alpha, TNFSF2, TNFalpha, Tnfa}, HLA-DRA (major histocompatibility complex, class II, DR alpha) [NCBI Gene 3122] {aka HLA-DRA1}, PRF1 (perforin 1) [NCBI Gene 5551] {aka HPLH2, P1, PFP}, AQP3 (aquaporin 3 (Gill blood group)) [NCBI Gene 360] {aka AQP-3, GIL}, CD72 (CD72 molecule) [NCBI Gene 971] {aka CD72b, LYB2}, FGFBP2 (fibroblast growth factor binding protein 2) [NCBI Gene 83888] {aka HBP17RP, KSP37}, TCL1A (TCL1 family AKT coactivator A) [NCBI Gene 8115] {aka TCL1}, JCHAIN (joining chain of multimeric IgA and IgM) [NCBI Gene 3512] {aka IGCJ, IGJ, JCH}, TLN1 (talin 1) [NCBI Gene 7094] {aka ILWEQ, TLN, talin-1}, KLRK1 (killer cell lectin like receptor K1) [NCBI Gene 22914] {aka CD314, D12S2489E, KLR, NKG2-D, NKG2D}, TNF (tumor necrosis factor) [NCBI Gene 7124] {aka DIF, IMD127, TNF-alpha, TNFA, TNFSF2, TNLG1F}, CD79B (CD79b molecule) [NCBI Gene 974] {aka AGM6, B29, IGB, Igbeta}, CD8A (CD8 subunit alpha) [NCBI Gene 925] {aka CD8, CD8alpha, IMD116, Leu2, p32}, CST7 (cystatin F) [NCBI Gene 8530] {aka CMAP}, NCAM1 (neural cell adhesion molecule 1) [NCBI Gene 4684] {aka CD56, MSK39, NCAM}, GZMA (granzyme A) [NCBI Gene 3001] {aka CTLA3, HFSP}, TGFB1 (transforming growth factor beta 1) [NCBI Gene 7040] {aka CAEND1, CED, DPD1, IBDIMDE, LAP, TGF-beta1}, GNLY (granulysin) [NCBI Gene 10578] {aka D2S69E, LAG-2, LAG2, NKG5, TLA519}, GLIPR2 (GLI pathogenesis related 2) [NCBI Gene 152007] {aka C9orf19, GAPR-1, GAPR1, gliPR 2}, CD79A (CD79a molecule) [NCBI Gene 973] {aka IGA, IGAlpha, MB-1, MB1}, GZMH (granzyme H) [NCBI Gene 2999] {aka CCP-X, CGL-2, CSP-C, CTLA1, CTSGL2}, Bcl2 (B cell leukemia/lymphoma 2) [NCBI Gene 12043] {aka Bcl-2, C430015F12Rik, D630044D05Rik, D830018M01Rik}
- **Diseases:** type 2 diabetes mellitus (MESH:D003924), AADs (MESH:C564653), stroke (MESH:D020521), osteoarthritis (MESH:D010003), infectious diseases (MESH:D003141), chronic inflammation (MESH:D007249), neurodegenerative disorders (MESH:D019636), died (MESH:D003643), atherosclerosis (MESH:D050197), cardiovascular diseases (MESH:D002318), myocardial infarction (MESH:D009203), Alzheimer's disease (MESH:D000544), cancers (MESH:D009369), endothelial dysfunction (MESH:D014652), cytotoxic (MESH:D064420), insulin resistance (MESH:D007333), neuroinflammatory (MESH:D000090862)
- **Chemicals:** lipid (MESH:D008055), agarose (MESH:D012685), TRIzol (MESH:C411644), glucose (MESH:D005947), AMPure XP (-), EDTA (MESH:D004492)
- **Species:** Cytomegalovirus (genus) [taxon 10358], Mus musculus (house mouse, species) [taxon 10090], Homo sapiens (human, species) [taxon 9606]
- **Mutations:** rs3793537, rs8019902
- **Cell lines:** S2 — Drosophila melanogaster (Fruit fly), Spontaneously immortalized cell line (CVCL_Z232)

## Figures

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

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