Biochemistry, Genetics and Molecular Biology

A Primate-Specific lncRNA LINC01021 Contributes to Cellular and Organismal Aging via DAZAP1-Dependent Destabilization of RBMX.

Zhang Y, Hu L, Dong X, Zeng Q, Zi M, Nisar A, Khan S, Bai R, Liu C, Ge M, Pu S, Li G, He Y. Published July 1, 2026 CC-BY

Aging is characterized by progressive physiological decline and age-related pathologies, yet the molecular determinants underlying lineage- and species-specific aging traits remain poorly understood. Although protein-coding regulators have dominated aging research, the contribution of long non-coding RNAs (lncRNAs), particularly primate-specific lncRNAs, has not been systematically explored. Here, through evolutionary screening and cross-species aging-associated analyses, we identified a set of primate-specific lncRNAs (including LINC01021, CTC-575 l10.1, CTA-150C2.13, and RP11-305F18.1, etc.) associated with human aging, and we functionally characterized LINC01021 as a representative candidate to assess their causal involvement. In human cells, LINC01021 promotes cellular senescence, whereas its silencing attenuates senescence-associated phenotypes. Mechanistically, LINC01021 is predominantly located in the nucleus, where it facilitates DAZAP1-dependent destabilization of RBMX mRNA, leading to activation of the P53 pathway and induction of canonical senescence features. At the organismal level, ectopic expression of human LINC01021 in mice contributes to aging-like phenotypes, including increased frailty and impaired motor coordination. Together, these findings implicate primate-specific lncRNAs in lineage-restricted aging and highlight an evolutionarily recent regulatory layer that may modulate aging trajectories.

Introduction

Aging is a complex, multifaceted biological process marked by progressive physiological decline and increased vulnerability to chronic diseases (Pal and Tyler2016). It is driven by an interconnected set of hallmarks, including genomic instability, telomere attrition, epigenetic alterations, loss of proteostasis, deregulated nutrient sensing, mitochondrial dysfunction, cellular senescence, stem‐cell exhaustion, and disrupted intercellular communication, which together shape tissue dysfunction and age‐related pathologies (López‐Otín et al.2023). These hallmarks reflect the cumulative outcome of dynamic regulatory processes that operate across multiple molecular layers (López‐Otín et al.2023). Although aging research has traditionally emphasized protein‐coding genes and canonical signaling pathways, growing evidence indicates that additional, under‐characterized regulatory layers make important contributions to aging biology (X. Li et al.2023; López‐Otín et al.2023; Wagner et al.2024).

Non‐coding RNAs (ncRNAs) have emerged as a prominent, yet incompletely explored, regulatory layer in aging (Guo et al.2022; Wagner et al.2024; K. Wang et al.2022). NcRNAs are broadly classified into small ncRNAs and long ncRNAs (lncRNAs) (Leng et al.2024; Mattick et al.2023; Statello et al.2021). Small ncRNAs such as microRNAs primarily act post‐transcriptionally to modulate mRNA stability and translation, thereby fine‐tuning pathways involved in cell‐cycle control, stress responses, and metabolic homeostasis (Turko et al.2025). By contrast, lncRNAs, transcripts longer than 200 nucleotides with limited protein‐coding potential, exhibit remarkable structural and functional versatility (Kopp and Mendell2018; Yao et al.2019). LncRNAs regulate gene expression at multiple levels: they recruit chromatin modifiers to shape epigenetic states, scaffold transcriptional regulators to control transcription, and influence post‐transcriptional events through interactions with microRNAs and RNA‐binding proteins; in some cases, they even encode functional micropeptides (Mattick et al.2023; Yao et al.2019). Because of this multilayered regulatory potential, lncRNAs are well positioned to coordinate large gene networks and to integrate diverse aging‐relevant processes such as DNA repair, inflammation, and cell‐cycle regulation (López‐Otín et al.2023; Sherazi et al.2023; Statello et al.2021). Consistent with a systems‐level role, lncRNAs constitute a major fraction of the ncRNA transcriptome and are frequently dysregulated with age across multiple tissues (Marttila et al.2020; Mattick et al.2023; Mercer et al.2009; Sherazi et al.2023).

Explaining how conserved molecular processes produce dramatically different aging trajectories across species remains an unresolved challenge. Maximum lifespan varies dramatically across mammals, spanning from approximately 2–4 years in mice to over 80 years in humans and more than 200 years in bowhead whales (Firsanov et al.2025; Keane et al.2015; López‐Otín et al.2023; Seluanov et al.2018), yet protein‐coding genes are highly conserved across species, suggesting that coding sequence variation alone cannot account for such divergence (Breschi et al.2017; Yue et al.2014). In contrast, lncRNAs display greater evolutionary plasticity, with rapid sequence turnover, frequent lineage‐restricted emergence, and flexible regulatory roles (Mattick et al.2023; Statello et al.2021; Yao et al.2019). Large‐scale comparative genomic studies show that total ncRNA content, driven primarily by lncRNAs, positively correlates with maximum lifespan across mammals, whereas protein‐coding sequence length does not show a similar positive relationship (A. Wang2025). Notably, the human genome contains an approximately 2.7‐fold expansion of lncRNA sequence compared with the mouse, which implicates lncRNA innovation as a plausible contributor to the evolution of extended lifespan and species‐specific aging phenotypes (A. Wang2025). Comparative investigations of exceptionally long‐lived species (e.g., naked mole‐rats, bats, and bowhead whales) further suggest that lineage‐specific regulatory innovations, rather than conserved protein‐coding changes alone, may underlie species‐restricted aging trajectories, with rapidly evolving lncRNAs representing a particularly compelling candidate layer (Cagan et al.2022; Jiang and Kong2020; Seluanov et al.2018; Yu et al.2021).

Primate lineages, and humans in particular, exhibit distinctive aging features, such as pronounced cortical expansion, extended post‐reproductive lifespan, and elevated susceptibility to neurodegenerative diseases, that are thought to arise primarily from lineage‐specific regulatory architectures rather than protein‐coding differences (Field et al.2019; Pollen et al.2023; Rigby Dames et al.2023; Vickery et al.2024). Supporting this view, multiple primate‐specific or primate‐enriched lncRNAs have been linked to aging‐relevant processes. For example, the primate‐specific lncRNAPCAT14, derived from human‐specific HERV‐H elements, is highly expressed and nuclear‐localized in young human endothelial cells but is downregulated during endothelial aging (Drekolia et al.2022). This change impairs migration and angiogenesis as well as increases inflammatory adhesion, all hallmarks of vascular aging that non‐primate models do not fully recapitulate (Drekolia et al.2022). In the developing nervous system, primate‐enriched lncRNAs show tightly regulated, transient expression during human cortical differentiation and actively shape cell type‐specific programs involved in synaptogenesis and neuronal connectivity, processes that are among the earliest and most vulnerable targets of brain aging and neurodegeneration (Field et al.2019). Likewise, primate‐cortex‐specific lncRNAs such asLINC00507display age‐dependent expression changes and associate with neurodegeneration‐related pathways, further linking primate regulatory innovation to both extended lifespan and increased age‐associated disease susceptibility (Mills et al.2016; Yan et al.2020).

Motivated by these observations, we performed an evolutionary screen for aging‐associated transcripts across species and identified a set of primate‐specific lncRNA candidates with robust age‐dependent signatures, includingLINC01021,CTC‐575 l10.1,CTA‐150C2.13,RP11‐305F18.1, andRP11‐700 N1.1. We prioritizedLINC01021as a top candidate for detailed functional and mechanistic dissection. Previous studies have reported thatLINC01021can suppress P53 signaling in certain cancer contexts, primarily through MDM2‐associated mechanisms, thereby promoting cell proliferation and tumorigenesis (Kaller et al.2024; Kaller et al.2017; X. L. Li et al.2017). However, our evolutionary and functional analyses reveal a striking context‐dependent duality: in aging‐associated contexts,LINC01021exhibits a strong positive association with P53 activation. Mechanistically, we show thatLINC01021accelerates cellular senescence by promoting DAZAP1‐dependent destabilization of the RNA‐binding protein RBMX, thereby activating P53 signaling. Importantly, ectopic expression of humanLINC01021in a transgenic mouse model is associated with features of accelerated functional decline, evidenced by increased frailty and impaired motor coordination.

Materials and Methods

LncRNASequence Collection

LncRNA sequences were obtained from the NONCODE database (version 6) for multiple representative species spanning primates, mammals, vertebrates, and invertebrates (Zhao et al.2021). NONCODE annotates lncRNAs based on transcriptomic evidence, transcript structure, length, and computational assessment of coding potential (Zhao et al.2021).LINC01021is annotated in NONCODE as transcript NONHSAT100800.2, comprising 1,112 nt across five exons. Additionally, the transcript exhibits a negative CNCI score (−0.1515520), supporting its low protein‐coding potential and corroborating its annotation as a lncRNA in the NONCODE database. Human tissue‐level RNA‐seq expression data were downloaded from the Genotype‐Tissue Expression (GTEx) project (Consortium2013).

Pairwise Sequence Alignment and Conservation Assessment

To evaluate evolutionary conservation of lncRNAs across species, all species were subjected to exhaustive pairwise sequence comparisons. For each species pair, nucleotide‐level alignments were performed using BLASTN. For each query lncRNA, significant alignments were retained based on sequence similarity, alignment length, and coverage thresholds to ensure biologically meaningful matches. For each query transcript, only the highest‐scoring alignment in the target species was retained. Reciprocal best hit (RBH) analysis was then performed to define putative orthologous lncRNA pairs. A lncRNA pair was considered conserved between two species if each transcript represented the best hit of the other in reciprocal searches. The number of RBH‐supported lncRNA pairs was summarized for all species combinations to characterize the global conservation landscape. LncRNAs were defined as primate‐conserved if they showed reciprocal best‐hit support among primate species but lacked detectable RBH relationships in non‐primate species included in this study. These primate‐specific conserved lncRNAs were subsequently used for downstream aging‐associated expression analyses.

Identification of Aging‐AssociatedlncRNAs

Human tissue RNA‐seq expression data were obtained from the GTEx project. Genomic coordinates of primate‐conserved lncRNAs were intersected with GTEx gene annotations to identify corresponding transcripts, requiring at least 90% reciprocal overlap to ensure high‐confidence mapping. Only lncRNAs satisfying this criterion were retained for downstream analyses. For each tissue, samples with available age information were included. Age was modeled as a continuous variable, and sex was incorporated as a covariate when applicable. Lowly expressed lncRNAs were excluded before modeling to enhance statistical stability. To characterize age‐dependent expression dynamics, generalized additive models (GAMs) were fitted independently for each lncRNA within each tissue according to the following model:

Expressioni=β0+s(Agei)+β1Sexi+ϵi\text{Expression}_{i} = \beta_{0} + s \left(Age_{i}\right) + \beta_{1} Sex_{i} + \epsilon_{i}

where s(Age) represents a penalized regression spline capturing non‐linear age effects, and Sex was included when applicable. The statistical significance of the smooth age term was evaluated for each lncRNA. Multiple testing correction was performed using the Benjamini–Hochberg procedure. LncRNAs with significant age‐dependent effects after correction were defined as aging‐associated lncRNAs.

Cell Culture and Transfection

Human embryonic lung fibroblasts (HELFs) were obtained from the Kunming Cell Bank of Type Culture Collection (Kunming, China). Primary human dermal fibroblasts (HDFs) were isolated from circumcised foreskins of healthy male donors (aged 5–20 years), as previously described (Xiao et al.2022). Cells were cultured in Dulbecco's modified Eagle medium (DMEM; C11995500BT, Gibco, MD, USA) supplemented with 10% fetal bovine serum (FBS; RYS‐F22‐05, Royacel Biotechnology, Lanzhou, China) and 1% penicillin–streptomycin (15140–122, Gibco, MD, USA) at 37°C in a humidified incubator with 5% CO2. Routine monitoring for morphology and contamination was performed.

For siRNA or plasmid transfection, cells were seeded in 6‐well plates (3–5 × 105cells per well) and transfected at ~80%–90% confluence using Lipofectamine 2000 (11668019, Thermo Fisher Scientific, MA, USA) according to the manufacturer's protocol. For siRNA transfection, 100 pmol of siRNA was used per well, while 2 μg of plasmid DNA was used for plasmid transfection. Medium was replaced with fresh DMEM 6 h after transfection.

Cells transfected with theLINC01021‐overexpression (LINC01021‐OE) plasmid were harvested 6 d post‐transfection, while other conditions were collected after 3 d. siRNA sequences are listed in TableS1.

Transcriptome Data Processing

RNA sequencing (RNA‐seq) datasets (accession number:GSE63577) were obtained from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/). These datasets include transcriptomic profiles of young and senescent fibroblast cell lines, including human foreskin fibroblasts (HFF), BJ foreskin fibroblasts (BJ), IMR‐90 fetal lung fibroblasts (IMR90), and MRC‐5 lung fibroblasts (MRC5). Differential expression analysis was conducted in RStudio (R 4.4.3) using the DESeq2 package. Genes significantly upregulated in senescent cells were defined as those with an adjustedpvalue < 0.05 and log2Fold Change (log2FC) > 0.

Senescent Cell Induction

Senescent cells were induced by three methods: ionizing radiation (IR), doxorubicin (DOX) treatment, and replicative exhaustion, as previously described (Y. He et al.2020). Low‐passage HELF and HDF fibroblasts (< 25 passages) were used as non‐senescent controls or for senescence induction. IR‐induced senescence was established by exposing cells at ~70% confluence to 15 Gy of ionizing radiation using a small animal radiation research platform (Xstrahl Inc., Camberley, UK). Cells were subcultured at a 1:3 ratio 3 d after irradiation, with full senescence typically observed ~7 d later. DOX‐induced senescence was established by treating cells with 250 nM DOX (HY‐15142, MCE, Shanghai, China) for 24 h. After treatment, the medium was replaced with standard DMEM, and cells were cultured for an additional 7 d. Replicative senescence was induced by serially passaging HELF cells until cell division ceased, typically after approximately 35 passages.

SA‐β‐Gal Staining

SA‐β‐Gal staining was performed using a Senescence β‐Galactosidase Staining Kit (C0602, Beyotime Biotechnology, Shanghai, China) following the manufacturer's instructions. Briefly, cells were seeded in six‐well plates and washed with PBS, fixed at room temperature for 7 min, and then incubated with the staining solution overnight at 37°C. The stained cells were imaged using a microscope (SMZ‐171‐BLED, Motic, Xiamen, China). Blue‐stained cells were considered senescent.

Quantitative Real‐TimePCR(qRT‐PCR)

Total RNA was isolated from cultured cells using TRIzol reagent (15596018, Invitrogen, CA, USA) according to the manufacturer's instructions. RNA concentration and purity were assessed spectrophotometrically. Equal amounts of RNA were reverse‐transcribed into cDNA using the HP All‐in‐One qRT Master Mix II Kit (RT203, YoungGen, Kunming, China). Quantitative PCR (qPCR) was performed using 2 × TSINGKE Master qPCR Mix (SYBR Green I) (TSE201, TSINGKE, Beijing, China) on a CFX Connect Real‐Time PCR System (Bio‐Rad, CA, USA), withGAPDHserving as the internal control. Relative gene expression levels were calculated using the 2^−ΔΔCtmethod. Primer sequences used in this study are provided in TableS2.

Cell Counting

Cell counts were determined manually using a hemocytometer under a light microscope. Cell suspensions were mixed and diluted to an appropriate concentration for counting. An aliquot of the diluted sample was loaded onto a hemocytometer, and cells were counted in duplicate at 100× magnification. The average number of cells per grid was used to calculate the final concentration using the standard hemocytometer formula.

EdUAssay

Cell proliferation was assessed using the Cell‐Light EdU Apollo567 In Vitro Kit (C10310‐1, RiboBio, Guangzhou, China) according to the manufacturer's instructions. Cells were seeded into 96‐well plates at a density of 3000 cells per well and allowed to adhere overnight. The cells were then incubated with the EdU working solution for 2 h, fixed with 4% paraformaldehyde for 15 min at room temperature, and permeabilized with 0.5% Triton X‐100. EdU incorporation was detected by staining with Apollo567, and nuclei were counterstained with DAPI. Fluorescence images were acquired using a fluorescence microscope.

Western Blot

Cells were washed twice with cold 1 × phosphate‐buffered saline (PBS; B540627‐0500, Sangon Biotech, Shanghai, China) and lysed on ice in RIPA buffer (BP‐115DG, Boston BioProducts, MA, USA) supplemented with PMSF (ST507, Beyotime Biotechnology, Shanghai, China) for 30 min. The lysates were sonicated at 35 W for 10 s, followed by centrifugation at 13,000 × g for 15 min at 4°C to remove insoluble debris. Protein concentrations were determined using a BCA Protein Assay Kit (P0010, Beyotime Biotechnology, Shanghai, China). Equal amounts of protein were mixed with 5 × Omni‐Easy Instant SDS‐PAGE Sample Buffer (LT101, EpiZyme, Shanghai, China), denatured at 95°C for 5–10 min, and separated by SDS‐PAGE on 12% polyacrylamide gels. Proteins were transferred onto 0.22‐μm PVDF membranes (Millipore, Bedford, MA, USA) using a wet transfer system (Bio‐Rad, CA, USA). Membranes were blocked with 5% (w/v) skim milk in TBS containing 0.1% Tween‐20 (TBS‐T) for 2 h at room temperature and incubated with primary antibodies overnight at 4°C. After washing with TBS‐T, membranes were incubated with HRP‐conjugated secondary antibodies for 1 h at room temperature. Protein bands were visualized using the SCG‐W3000 Plus chemiluminescence imaging system (Servicebio, Wuhan, China) and quantified with ImageJ software (NIH, MD, USA). All antibodies used in this study are listed in TableS3.

Plasmid Construction

To generate overexpression plasmids, the full‐length sequences ofLINC01021,RPL23A,SKP1,EPB41L2, andDAZAP1were individually cloned into the pCDH‐CMV‐MCS‐EF1‐puro‐3×FLAG‐3×HA vector.

RNAPull‐down and Mass Spectrometry (MS) Analysis

RNA was transcribed in vitro using T7 RNA polymerase (K102A01, Perfect mRNA, Hangzhou, China) and labeled at the 3′‐end with biotin using the Pierce RNA 3′ End Desthiobiotinylation Kit (20163, Thermo Fisher Scientific, MA, USA), according to the manufacturer's protocol. RNA pull‐down assays were conducted as previously described (P. Wang et al.2014). Briefly, 50 pmol of biotin‐labeled RNA was refolded in RNA structure buffer (10 mM Tris, pH 7.0; 100 mM KCl; 10 mM MgCl2) by heating at 85°C for 5 min and snap‐cooling on ice. RNA‐protein complexes were captured using the Pierce Magnetic RNA‐Protein Pull‐Down Kit (20164, Thermo Fisher Scientific, MA, USA), washed, eluted, and denatured in 5 × SDS loading buffer. The recovered proteins were separated by gradient SDS‐PAGE and analyzed by MS for protein identification.

Chromatin Immunoprecipitation (ChIP)‐qPCR

LINC01021‐OE HELFs were crosslinked by incubating with 1% formaldehyde for 5 min at 37°C, followed by quenching with 0.125 M glycine for 10 min at room temperature. Cells were then washed twice with ice‐cold PBS containing PMSF and lysed using SDS lysis buffer from the ChIP Assay Kit (P2078, Beyotime Biotechnology, Shanghai, China). Chromatin was fragmented using a Bioruptor Plus sonicator (Diagenode, Seraing, Belgium) at 40 W for 30 cycles (10 s on/10 s off) at 4°C. The sheared chromatin extracts were precleared with Protein A/G agarose beads pre‐blocked with salmon sperm DNA for 30 min at 4°C. After removing 2% of the input as a control, the lysates were incubated overnight at 4°C with either an anti‐RBMX antibody (F1965, Selleckchem, TX, USA) or a control IgG (A7016, Beyotime Biotechnology, Shanghai, China). Immune complexes were captured with Protein A/G agarose beads for 2 h at 4°C and washed sequentially with low‐salt, high‐salt, LiCl, and TE buffers. DNA–protein complexes were eluted with 1% SDS / 0.1 M NaHCO3and decrosslinked by incubation with 0.2 M NaCl at 65°C for 4 h. The samples were treated with proteinase K at 45°C for 1 h, and the purified DNA was isolated using a commercial DNA purification kit (D6492, OMEGA BIO‐TEK, GA, USA). Quantified DNA was analyzed by qPCR, with primer sequences provided in TableS2.

Protein Half‐Life Assay

Protein stability was assessed as previously described (Q. Zhang et al.2020). Cells transfected with the indicated plasmids were treated with 50 μg/mL cycloheximide (CHX; GC17198, GLPBIO, CA, USA) at the specified time points. At each time point, cells were harvested, lysed in RIPA buffer with protease inhibitors, and denatured in SDS loading buffer at 95°C for 5 min. Protein levels were analyzed by immunoblotting with the specified antibodies.

RNAHalf‐Life Assay

Cells transfected with the indicated plasmids were treated with 10 μg/mL Actinomycin D (Act D; GC16866, GLPBIO, CA, USA) to inhibit transcription. Total RNA was extracted at the indicated time points, and RNA stability was assessed by RT‐qPCR using specific primers. Relative RNA levels at each time point were normalized to the 0 h time point and plotted to evaluate RNA decay.

Protein Degradation Pathway Analysis

HELFs were transfected with theLINC01021‐OE plasmid for 6 d to achieve robust overexpression. Cells were then treated with 1 μM MG132 (HY‐13259, MCE, Shanghai, China), 10 μM chloroquine (HY‐17589AR, MCE, Shanghai, China), or 5 μM Q‐VD‐OPh (HY‐12305, MCE, Shanghai, China) for 12 h before being harvested for Western blot analysis.

Animals

Transgenic mice with a knock‐in (KI) allele forLINC01021‐OE were obtained from Shanghai Model Organisms Center Inc. (Shanghai, China). These mice were generated by CRISPR/Cas9‐mediated homologous recombination, inserting a CAG‐LINC01021‐polyA (PA) cassette into the Rosa26 locus on a C57BL/6J background. The CAG promoter drives constitutive and ubiquitous expression ofLINC01021across tissues, independent of endogenous regulatory elements (Hasegawa et al.2013). Briefly, Cas9 mRNA and guide RNA (gRNA) were synthesized by in vitro transcription, and a donor vector containing a 3.3 kb 5′ homology arm, the CAG‐LINC01021‐PA cassette, and a 3.3 kb 3′ homology arm was constructed via In‐Fusion cloning. Cas9 mRNA, gRNA, and donor vector were microinjected into C57BL/6J zygotes to generate F0 founders, which were confirmed by PCR and Sanger sequencing. These founders were backcrossed to wild‐type C57BL/6J mice to obtain heterozygous F1 progeny. All mice were maintained under specific pathogen‐free (SPF) conditions at the Kunming Institute of Zoology, Chinese Academy of Sciences. Animal experiments were approved by the Animal Ethics Committee of the Kunming Institute of Zoology, Chinese Academy of Sciences (IACUC‐RE‐2025‐08‐008) and conducted following the Guide for the Care and Use of Laboratory Animals.

Mouse Frailty Index (FI) Assessment

Mouse FI was assessed using a validated clinical deficit accumulation approach, as previously described (Feridooni et al.2015). The FI comprised 23 health deficits spanning multiple physiological systems, including integument, musculoskeletal, auditory, ocular, digestive, urogenital, and body weight (see TableS4for detailed scoring criteria). Each deficit was scored as 0 (absent), 0.5 (mild), or 1 (severe) according to predefined criteria. Body weight was recorded as a continuous numerical variable for each mouse. Forelimb grip strength was measured three times using a grip strength meter, and the mean value was used for analysis. These measurements were reported as raw physiological readouts and were not converted into standard deviation‐based categorical scores. The FI for each mouse was calculated by summing the individual deficit scores and dividing by the total number of assessed items (23), yielding a final FI value ranging from 0 (no deficits) to 1 (all deficits present).

Behavioral Assessments

Beam Walking Test

Motor coordination and balance were assessed using the beam walking test. Mice were trained to traverse a horizontal acrylic beam (1 cm width, 100 cm length) elevated 50 cm above the surface. During the test session, both the latency to traverse the full length of the beam and the number of paw slips were recorded. Each mouse performed three consecutive trials with a 10‐min intertrial interval, and the mean value was used for subsequent analysis.

Rotarod Test

Motor performance and endurance were evaluated using an accelerating rotarod apparatus (SA102, SANS, Nanjing, China). On day 1, mice were habituated to the apparatus at a constant speed of 10 rpm for 10 min. On day 2, mice underwent training at 10 rpm for 5 min, followed by stepwise acceleration at 5 rpm every 2 min. During the test phase, mice were placed on the stationary rod, which accelerated from 4 rpm to 35 rpm over 4 min, including a 1‐min constant‐speed period at 4 rpm before acceleration. Each mouse completed three independent trials with 30‐min rest intervals. Latency to fall was recorded for each trial.

Treadmill Fatigue Test

Exercise capacity was assessed using a motorized treadmill system (SA101, SANS, Nanjing, China). Mice were acclimated for three consecutive days at a speed of 5 m/min for 5 min per day without incline and with a mild electrical stimulus (0.3 mA). On the test day, mice were subjected to a 5‐degree incline and began running at 10 m/min for 10 min. The speed was subsequently increased by 2 m/min every 5 min until exhaustion, defined as remaining on the shock grid for more than 5 s without resuming running.

Pole Test

The pole test was performed to assess bradykinesia and motor coordination. Mice were placed head‐upward on the top of a rough wooden pole (50 cm height, 8 mm diameter, 30‐degree inclination). The time required to turn downward and the total time to descend to the base were recorded. Mice underwent one day of training followed by three test trials (maximum duration of 60 s per trial, with 30‐min intertrial intervals). The average of the three trials was used for statistical analysis.

In Vivo Senescence Imaging

In vivo senescence imaging was performed using the Senescence‐Tracker Near‐Infrared (NIR) probe (C0603, Beyotime Biotechnology, Shanghai, China), as previously described (Hu et al.2023). The probe was dissolved in DMSO to generate a 10 mM stock solution and diluted in 1% DMSO, 2% Tween‐80, and 97% sterile saline to a final concentration of 5 μM. Mice received a single intravenous tail‐vein injection of 100 μL probe solution. Whole‐body fluorescence images were acquired 24 h post‐injection using an IVIS Lumina 3 imaging system (PerkinElmer, CT, USA) under identical acquisition settings. Regions of interest (ROIs) were manually defined over each mouse, and fluorescence intensity was quantified as total radiant efficiency using Living Image Software v4.4. All imaging and analysis parameters were kept constant across experimental groups.

Statistical Analysis

Data are presented as mean ± standard error of the mean (SEM). Comparisons between two groups were performed using an unpaired two‐tailed Student'st‐test, with Welch's correction applied when variances were unequal. Comparisons among more than two groups were conducted using one‐way analysis of variance (ANOVA). Ap‐value < 0.05 was considered statistically significant (*p< 0.05, **p< 0.01, ***p< 0.001; ns, not significant).

Results

Primate‐ConservedlncRNAsExhibit Widespread and Recurrent Associations With Aging Across Human Tissues

LncRNAs constitute a major class of regulatory transcripts in eukaryotic genomes but are characterized by rapid evolutionary turnover and pronounced lineage specificity (Mattick et al.2023; Statello et al.2021). To systematically assess the evolutionary distribution of lncRNAs, we compiled annotated lncRNA catalogs from the NONCODE database across multiple representative species. As shown in Figure1A, the number of annotated lncRNAs varies dramatically among species, ranging from only 52 in yeast to more than 96,000 in humans. Notably, substantial discrepancies were also observed among closely related primates: chimpanzee, gorilla, and orangutan each harbor approximately 6,000–15,000 annotated lncRNAs, markedly fewer than reported in humans. These differences are unlikely to reflect true biological absence, but instead largely arise from variation in sequencing depth, annotation strategies, and research intensity, underscoring the limitations of absolute lncRNA counts for inferring functional relevance across species.

Evolutionary distribution and conservation of lncRNAs across species. (A) Phylogenetic tree of representative species annotated with the number of lncRNAs identified in each genome. (B) Circular visualization of lncRNA conservation across species, with links indicating shared lncRNAs between species pairs. (C) Relative enrichment of conserved human lncRNAs in primates compared with non‐primate mammals. (D) Distribution of primate‐conserved lncRNAs exhibiting age‐associated expression across human tissues. (E) Heatmap showing age‐associated expression changes of primate‐conserved lncRNAs across human tissues under stringent selection criteria (FDR < 0.05 and log2FC > 1). Color scale indicates log2FC with age (red, upregulation; blue, downregulation). The bar plot above summarizes the number of tissues in which each lncRNA exhibits significant age‐associated expression.

Evolutionary distribution and conservation of lncRNAs across species. (A) Phylogenetic tree of representative species annotated with the number of lncRNAs identified in each genome. (B) Circular visualization of lncRNA conservation across species, with links indicating shared lncRNAs between species pairs. (C) Relative enrichment of conserved human lncRNAs in primates compared with non‐primate mammals. (D) Distribution of primate‐conserved lncRNAs exhibiting age‐associated expression across human tissues. (E) Heatmap showing age‐associated expression changes of primate‐conserved lncRNAs across human tissues under stringent selection criteria (FDR < 0.05 and log2FC > 1). Color scale indicates log2FC with age (red, upregulation; blue, downregulation). The bar plot above summarizes the number of tissues in which each lncRNA exhibits significant age‐associated expression.

To move beyond catalog size and directly evaluate evolutionary relationships, we performed pairwise comparisons of lncRNAs across species to assess sequence conservation (Figure1B). Consistent with prior observations, lncRNAs exhibited strong species specificity, with only a small fraction shared between any two species. Importantly, when human lncRNAs were used as the reference, a substantially larger fraction showed detectable conservation within primate species compared with non‐primate mammals. Based on this primate‐skewed conservation pattern, we identified a set of 7098 human lncRNAs that are conserved in at least one other primate species but absent from non‐primate species, hereafter referred to as primate‐conserved lncRNAs (Figure1C). These transcripts combine evolutionary constraints with lineage restriction, suggesting that they may encode regulatory functions preferentially relevant to primate biology.

To investigate whether primate‐conserved lncRNAs are implicated in aging, we analyzed their age‐associated expression patterns across multiple human tissues using GAMs. Across all tissues examined, a subset of primate‐conserved lncRNAs exhibited significant associations with chronological age (Figure1D), with the colon displaying the largest number of aging‐associated candidates. Integrative analysis further revealed that many of these lncRNAs showed age‐dependent expression changes across multiple tissues, with some exhibiting significant associations in more than 16 tissues (FigureS1A), indicating potential roles in systemic or organism‐wide aging regulation.

To define a high‐confidence set of aging‐associated candidates, we applied more stringent criteria requiring both statistical significance (FDR < 0.05) and substantial age‐dependent expression changes (log2FC > 1). This analysis identified a core cohort of primate‐conserved lncRNAs with robust and recurrent age‐associated expression across tissues (Figure1E), such asLINC01021,CTC‐575 l10.1,CTA‐150C2.13,RP11‐305F18.1,RP11‐700 N1.1,RP4‐541C22.5, andRP11‐291B21.2. Among these,LINC01021emerged as the most broadly associated transcript, exhibiting significant age‐related expression changes across seven tissues, including the adrenal gland, cervix uteri, and colon (Figure1E). In addition toLINC01021, several other primate‐conserved lncRNAs, such asCTC‐575 l10.1,CTA‐150C2.13, andRP11‐305F18.1, also displayed consistent age‐associated expression patterns across multiple tissues (Figure1E). Together, these findings positionLINC01021as a representative exemplar of a broader class of aging‐associated primate‐conserved lncRNAs rather than an isolated outlier. Based on its breadth of age association, recurrent induction across independent aging contexts, and suitability for functional interrogation,LINC01021was selected for subsequent mechanistic investigation.

LINC01021Is Robustly Upregulated Across Multiple Senescence Models in Human Fibroblasts

Following the identification of primate‐conserved lncRNAs with recurrent age‐associated expression across human tissues, we next examined whetherLINC01021, as a representative candidate, exhibits consistent induction in established human cellular models of senescence. We first analyzed a publicly available RNA‐seq dataset comprising four human fibroblast strains (HFF, BJ, IMR90, and MRC5) subjected to either replicative or stress‐induced senescence (Marthandan et al.2016). Across these well‐established models, which recapitulate core features of cellular aging,LINC01021emerged as a robust aging‐associated primate‐conserved lncRNA that was reproducibly upregulated across all four fibroblast strains under both senescence paradigms (Figure2A–D).

LINC01021is robustly upregulated across multiple senescence models in human fibroblasts. (A‐D) Volcano plots showing transcriptional changes between young and senescent fibroblasts in four independent human fibroblast strains (HFF, BJ, IMR90, and MRC5), highlighting consistent upregulation ofLINC01021. (E) Representative SA‐β‐Gal staining images of HELF cells subjected to IR‐induced, DOX‐induced, or replicative senescence. (F‐H) Quantification of SA‐β‐Gal+cells corresponding to the conditions shown in (E). (I‐N) qRT‐PCR analysis of senescence markersCDKN2AandCDKN1Ain HELF cells under IR‐induced (I and J), replicative (K and L), or DOX‐induced (M and N) senescence. (O‐Q) qRT‐PCR analysis showing induction ofLINC01021expression in HELF cells across the indicated senescence models. Data are shown as mean ± SEM from at least three independent experiments.

LINC01021is robustly upregulated across multiple senescence models in human fibroblasts. (A‐D) Volcano plots showing transcriptional changes between young and senescent fibroblasts in four independent human fibroblast strains (HFF, BJ, IMR90, and MRC5), highlighting consistent upregulation ofLINC01021. (E) Representative SA‐β‐Gal staining images of HELF cells subjected to IR‐induced, DOX‐induced, or replicative senescence. (F‐H) Quantification of SA‐β‐Gal+cells corresponding to the conditions shown in (E). (I‐N) qRT‐PCR analysis of senescence markersCDKN2AandCDKN1Ain HELF cells under IR‐induced (I and J), replicative (K and L), or DOX‐induced (M and N) senescence. (O‐Q) qRT‐PCR analysis showing induction ofLINC01021expression in HELF cells across the indicated senescence models. Data are shown as mean ± SEM from at least three independent experiments.

To experimentally validate these observations, we established three independent senescence models in HELFs, including IR‐induced senescence, DOX‐induced senescence, and replicative senescence induced by serial passaging. Successful induction of senescence was confirmed by a marked increase in SA‐β‐Gal positive cells (Figure2E–H) and significant upregulation of canonical senescence regulatorsCDKN2AandCDKN1Aat the mRNA level (Figure2I–N). Consistent with the bioinformatic analyses, qRT‐PCR revealed a robust and highly reproducible increase inLINC01021expression across all three senescence models (Figure2O–Q). Notably, the magnitude ofLINC01021induction was comparable to, and in some conditions exceeded that ofCDKN2AandCDKN1A, further supporting its tight association with the senescent state.

LINC01021Modulates Cellular Senescence in Human Fibroblasts

Given the consistent and reproducible induction ofLINC01021across multiple cellular senescence models, we next investigated whetherLINC01021functionally contributes to the establishment of senescent phenotypes in human fibroblasts. Gain‐of‐function experiments were first performed in HELFs, in which overexpression ofLINC01021was confirmed by qRT‐PCR (Figure3A). Compared with control cells,LINC01021‐OE HELFs exhibited a significantly increased proportion of SA‐β‐Gal positive cells, a reduced fraction of EdU‐positive proliferating cells, and an overall impairment of proliferative capacity, as reflected by slower cell growth rates (Figure3B–F). Consistent phenotypic effects were also observed in HDFs, in whichLINC01021‐OE similarly increased SA‐β‐Gal positivity and reduced EdU incorporation (FigureS2A–D), indicating that the pro‐senescent effects ofLINC01021are not restricted to a single fibroblast lineage.

LINC01021modulates cellular senescence in human fibroblasts. (A) qPCR analysis confirmingLINC01021overexpression in HELF cells. (B, C) Representative SA‐β‐Gal staining and quantification showing increased senescence uponLINC01021overexpression. (D, E) Representative EdU staining and quantification indicating reduced DNA synthesis inLINC01021‐OE HELF cells. (F) Proliferation curves comparingLINC01021‐OE HELF cells with vector controls. (G‐K) qRT‐PCR ofCDKN2A,CDKN1A, and SASP‐associated genes inLINC01021‐OE HELF cells. (L‐O) Western blot analysis and quantification of P16, P21, and P53 inLINC01021‐OE HELF cells. (P‐S) qRT‐PCR analysis ofLINC01021,CDKN2A,CDKN1A, and SASP‐associated genes in replicative senescent HELF cells followingLINC01021knockdown. (T‐W) Western blot analysis and quantification of P21, P16, and P53 in replicative senescent HELF cells afterLINC01021knockdown. Data are shown as mean ± SEM from at least three independent experiments. SASP, senescence‐associated secretory phenotype.

LINC01021modulates cellular senescence in human fibroblasts. (A) qPCR analysis confirmingLINC01021overexpression in HELF cells. (B, C) Representative SA‐β‐Gal staining and quantification showing increased senescence uponLINC01021overexpression. (D, E) Representative EdU staining and quantification indicating reduced DNA synthesis inLINC01021‐OE HELF cells. (F) Proliferation curves comparingLINC01021‐OE HELF cells with vector controls. (G‐K) qRT‐PCR ofCDKN2A,CDKN1A, and SASP‐associated genes inLINC01021‐OE HELF cells. (L‐O) Western blot analysis and quantification of P16, P21, and P53 inLINC01021‐OE HELF cells. (P‐S) qRT‐PCR analysis ofLINC01021,CDKN2A,CDKN1A, and SASP‐associated genes in replicative senescent HELF cells followingLINC01021knockdown. (T‐W) Western blot analysis and quantification of P21, P16, and P53 in replicative senescent HELF cells afterLINC01021knockdown. Data are shown as mean ± SEM from at least three independent experiments. SASP, senescence‐associated secretory phenotype.

At the molecular level,LINC01021overexpression induced the expression of core senescence effectors, includingCDKN2AandCDKN1A, and promoted the transcription of multiple senescence‐associated secretory phenotype (SASP) factors, such asIL1β,IL6, andCXCL2(Figure3G–K). Concordantly, protein levels of P16 and P21 were markedly increased followingLINC01021overexpression (Figure3L–N). Notably, P53 protein abundance was also elevated under these conditions (Figure3L,O), a pattern that contrasts with previous reports describing a suppressive role ofLINC01021on P53 activity in cancer cells, highlighting potential context‐dependent regulation (Han et al.2021).

To complement these gain‐of‐function analyses, loss‐of‐function experiments were performed across three independent senescence models. Knockdown ofLINC01021consistently reduced the expression ofCDKN2AandCDKN1Aand attenuated the induction of SASP‐associated genes (Figures3P–SandS2E–R). At the protein level, P16 and P21 were correspondingly decreased across all senescence models uponLINC01021depletion, whereas P53 protein levels decreased upon knockdown and increased upon overexpression, remaining positively correlated withLINC01021expression (Figures3T–WandS2S–Z). Together, these gain‐ and loss‐of‐function analyses demonstrate thatLINC01021actively modulates key cellular senescence features in human fibroblasts and reveal a context‐dependent relationship betweenLINC01021and P53 regulation that differs from its previously reported role in cancer cells (Kaller et al.2024; Li et al.2017).

RBMXMediatesLINC01021‐Dependent Senescence Regulation

To explore the molecular mechanism underlyingLINC01021‐mediated senescence regulation, particularly its impact on P53‐associated pathways, we first examined the subcellular localization ofLINC01021. Subcellular fractionation analyses revealed thatLINC01021is predominantly enriched in the nucleus (Figure4A), prompting us to focus on nuclear‐associated proteins in subsequent mechanistic investigations. We therefore performed RNA pull‐down coupled with MS to identifyLINC01021‐interacting proteins, which yielded 8 candidate binding partners with annotated nuclear localization (Figure4B,C).

RBMX is a key mediator ofLINC01021‐dependent senescence regulation. (A) Subcellular fractionation showing predominant nuclear localization ofLINC01021. (B) Schematic overview of the workflow for identifying and validatingLINC01021‐interacting proteins, including RNA pull‐down coupled with mass spectrometry (MS) and subsequent chromatin immunoprecipitation (ChIP) assays. (C) Candidate nuclear proteins identified byLINC01021RNA pull‐down and MS. (D) qRT‐PCR analysis of selected candidate mRNAs inLINC01021‐OE HELF cells. (E) qRT‐PCR analysis ofRBMX,CDKN2A,CDKN1A, and SASP‐associated genes in si‐RBMXHELF cells. (F, G) Representative SA‐β‐Gal staining and quantification afterRBMXknockdown. (H‐L) Western blot analysis and quantification of RBMX, P16, P21, and P53 in si‐RBMXHELF cells. Data are shown as mean ± SEM from at least three independent experiments.

RBMX is a key mediator ofLINC01021‐dependent senescence regulation. (A) Subcellular fractionation showing predominant nuclear localization ofLINC01021. (B) Schematic overview of the workflow for identifying and validatingLINC01021‐interacting proteins, including RNA pull‐down coupled with mass spectrometry (MS) and subsequent chromatin immunoprecipitation (ChIP) assays. (C) Candidate nuclear proteins identified byLINC01021RNA pull‐down and MS. (D) qRT‐PCR analysis of selected candidate mRNAs inLINC01021‐OE HELF cells. (E) qRT‐PCR analysis ofRBMX,CDKN2A,CDKN1A, and SASP‐associated genes in si‐RBMXHELF cells. (F, G) Representative SA‐β‐Gal staining and quantification afterRBMXknockdown. (H‐L) Western blot analysis and quantification of RBMX, P16, P21, and P53 in si‐RBMXHELF cells. Data are shown as mean ± SEM from at least three independent experiments.

To assess the functional relevance of these candidates inLINC01021‐associated senescence, we examined their expression changes inLINC01021‐OE HELFs. Among the 8 candidates, the mRNA levels of genes encoding 3 candidate proteins (NUDT21, DAZAP1, and RBMX) were negatively correlated withLINC01021expression, 3 (RPL23A, EPB41L2, and SKP1) were positively correlated, and 2 (FBL and DKC1) showed no significant changes (Figure4D). To functionally interrogate their roles, negatively correlated candidates were individually depleted, whereas positively correlated candidates were ectopically expressed in HELFs. Notably, knockdown ofRBMX, but notNUDT21orDAZAP1, resulted in robust induction of the senescence regulatorsCDKN2AandCDKN1A, along with multiple SASP genes, includingCXCL1,CXCL3,IGFBP5,MMP1, andCCL2(Figures4EandS3A,B). In contrast, overexpression ofRPL23A,EPB41L2, orSKP1failed to elicit comparable senescence‐associated transcriptional responses (FigureS3C–E) indicating RBMX as the primary functional mediator among the identifiedLINC01021‐interacting candidates.

Consistent with these transcriptional changes,RBMXknockdown markedly increased the proportion of SA‐β‐Gal‐positive HELFs and elevated protein levels of the senescence effectors P16 and P21 (Figure4F–K). Notably,RBMXknockdown also led to an increase in P53 protein (Figure4H,L), closely mirroring the molecular phenotype observed uponLINC01021overexpression. To assess P53 dependency, we knocked downp53via siRNA inLINC01021‐overexpressing cells and examined P53/P21 pathway activity.p53depletion markedly attenuatedLINC01021‐induced P21 activation (FigureS3F,G) indicating thatLINC01021promotes senescence at least in part through P53‐dependent signaling downstream of RBMX. Together, these findings position RBMX as a key effector downstream ofLINC01021that mediates senescence regulation via P53, and underscore the essential role of P53 in this pathway.

LINC01021RegulatesRBMXAbundance Predominantly ThroughDAZAP1‐Dependent Post‐Transcriptional Mechanisms

Given that RBMX functions as a key downstream mediator ofLINC01021‐induced cellular senescence, we next investigated howLINC01021regulates RBMX abundance. AlphaFold‐based structural modeling predicted a putative three‐dimensional association betweenLINC01021and RBMX (Figure5A); however, ChIP‐qPCR using an anti‐RBMX antibody revealed only low but reproducible enrichment ofLINC01021(Figures4Band5B), suggesting an indirect chromatin‐associated relationship rather than a stable RNA–protein interaction. Consistent with a regulatory effect,LINC01021overexpression in HELFs dramatically reduced RBMX protein levels (Figure5C,D) and shortened the RBMX protein half‐life (Figure5E), indicating thatLINC01021negatively regulates RBMX abundance.

LINC01021regulates RBMX abundance predominantly through DAZAP1‐dependent post‐transcriptional mechanisms. (A) AlphaFold‐based modeling predicting a potential structural interface betweenLINC01021and RBMX. (B) ChIP‐qPCR analysis showing low but reproducible chromatin‐associated enrichment ofLINC01021in RBMX immunoprecipitates. (C, D) Western blot analysis and quantification of RBMX protein levels inLINC01021‐OE HELF cells. (E) Cycloheximide chase analysis assessing RBMX protein stability uponLINC01021overexpression. (F, G) MG132 treatment partially restores RBMX protein levels inLINC01021‐OE HELF cells. (H‐J) Inhibition of lysosomal degradation (chloroquine), caspase‐dependent proteolysis (Q‐VD‐OPh), or SCF E3 ligase component SKP1 does not rescue RBMX protein levels. (K) Actinomycin D chase analysis showing acceleratedRBMXmRNA decay inLINC01021‐OE HELF cells. (L) Western blot analysis of DAZAP1 protein levels inLINC01021‐OE HELF cells. (M) RBMX protein levels followingDAZAP1knockdown. (N, O) qRT‐PCR and mRNA stability analyses demonstrating DAZAP1‐dependent regulation ofRBMXtranscript abundance. (P) RBMX protein levels followingDAZAP1overexpression. Data are shown as mean ± SEM from at least three independent experiments.

LINC01021regulates RBMX abundance predominantly through DAZAP1‐dependent post‐transcriptional mechanisms. (A) AlphaFold‐based modeling predicting a potential structural interface betweenLINC01021and RBMX. (B) ChIP‐qPCR analysis showing low but reproducible chromatin‐associated enrichment ofLINC01021in RBMX immunoprecipitates. (C, D) Western blot analysis and quantification of RBMX protein levels inLINC01021‐OE HELF cells. (E) Cycloheximide chase analysis assessing RBMX protein stability uponLINC01021overexpression. (F, G) MG132 treatment partially restores RBMX protein levels inLINC01021‐OE HELF cells. (H‐J) Inhibition of lysosomal degradation (chloroquine), caspase‐dependent proteolysis (Q‐VD‐OPh), or SCF E3 ligase component SKP1 does not rescue RBMX protein levels. (K) Actinomycin D chase analysis showing acceleratedRBMXmRNA decay inLINC01021‐OE HELF cells. (L) Western blot analysis of DAZAP1 protein levels inLINC01021‐OE HELF cells. (M) RBMX protein levels followingDAZAP1knockdown. (N, O) qRT‐PCR and mRNA stability analyses demonstrating DAZAP1‐dependent regulation ofRBMXtranscript abundance. (P) RBMX protein levels followingDAZAP1overexpression. Data are shown as mean ± SEM from at least three independent experiments.

To assess whether RBMX destabilization involves proteolytic pathways, we examined the major mechanisms of protein degradation. Pharmacological inhibition of the proteasome using MG132 marginally restored RBMX protein levels inLINC01021‐OE cells (Figure5F,G), whereas inhibition of lysosomal degradation or caspase‐dependent proteolysis had no detectable effect (Figures5H,IandS4A,B). SKP1, a core adaptor component of SCF (SKP1‐CUL1‐F‐box) E3 ubiquitin ligase complexes that mediate proteasomal protein turnover (Zeng et al.2025), was identified among candidateLINC01021‐interacting proteins in the RNA pull‐down assay; however,SKP1knockdown failed to rescue RBMX protein levels (Figures5JandS4C), indicating that proteasome‐mediated degradation contributes weakly and in an SKP1‐independent manner. In parallel, qPCR analysis revealed thatLINC01021overexpression significantly reducedRBMXmRNA levels (Figure4D), prompting us to examine post‐transcriptional regulation.

Act D chase experiments demonstrated thatLINC01021markedly acceleratedRBMXmRNA decay (Figure5K), indicating a substantial contribution of mRNA destabilization toRBMXdownregulation. Among proteins identified in the RNA pull‐down screen, the RNA‐binding protein DAZAP1, previously implicated in mRNA stability control (Qiu, Li, et al.2025; Wang et al.2021; Zhang et al.2024), was significantly reduced at both mRNA and protein levels uponLINC01021overexpression (Figure4D; Figures5LandS4D). Functionally,DAZAP1depletion decreased RBMX protein and mRNA levels and accelerated RBMX transcript decay (Figures5M–OandS4E,F), whereasDAZAP1overexpression robustly increased RBMX protein abundance (Figures5PandS4G,H). Together, these findings indicate thatLINC01021suppresses RBMX expression predominantly through a DAZAP1‐dependent post‐transcriptional mechanism governingRBMXmRNA stability, with a limited contribution from proteasome‐mediated protein turnover.

LINC01021Overexpression Promotes Aging‐Like Phenotypes in Humanized Transgenic Mice

To assess whether ectopic expression of the primate‐specific lncRNALINC01021is sufficient to elicit aging‐associated phenotypes in vivo, we generated a humanizedLINC01021KI mouse model using CRISPR‐Cas9. Correct transgene integration was confirmed by PCR genotyping, andLINC01021expression in the spleen was verified by qRT‐PCR (Figures6A,BandS5A), as rodents lack an endogenous ortholog of this primate‐specific lncRNA. Body weight measurements revealed no significant differences between KI and WT mice at the examined ages, indicating that overall growth was not significantly affected (Figures6CandS5B). To further quantify senescence burden in vivo, we applied near‐infrared senescence imaging using the activatable near‐infrared β‐galactosidase probe, a β‐galactosidase‐responsive fluorescent sensor that emits in the NIR spectrum upon enzymatic cleavage by SA‐β‐Gal, enabling real‐time detection of senescent cells in deep tissues (Hu et al.2023). This analysis revealed a significantly elevated senescence signal in young adult KI/KI mice compared with WT controls, whereas aged KI mice exhibited a trend toward increased signal intensity that did not reach statistical significance (Figures6D,EandS5C,D), suggesting thatLINC01021expression may contribute to the onset of cellular senescence in vivo.

LINC01021overexpression promotes aging‐like phenotypes in humanized transgenic mice. (A) Schematic of the humanizedLINC01021‐KI mouse model. (B) qRT‐PCR analysis of humanLINC01021expression in spleen from young WT, KI/WT, and KI/KI mice. (C) Body weight of young mice. (D, E) Representative NIR fluorescence imaging using a SA‐β‐Gal responsive probe and quantification of senescence‐associated signal intensity in young mice. (F) Frailty index in young mice. (G‐I) Western blot analysis and quantification of P21 and P16 protein levels in spleen from young mice. (J, K) Representative images and quantification of fur depigmentation in 10‐month‐old male mice. (L) Forelimb grip strength in aged mice. (M‐O) Balance and coordination performance in aged mice assessed by balance beam traversal time (M), pause frequency (N), and foot slips (O). (P–T) Hematological parameters in aged mice, including white blood cells (WBC), lymphocytes, granulocytes, monocytes, and red blood cells (RBC).

LINC01021overexpression promotes aging‐like phenotypes in humanized transgenic mice. (A) Schematic of the humanizedLINC01021‐KI mouse model. (B) qRT‐PCR analysis of humanLINC01021expression in spleen from young WT, KI/WT, and KI/KI mice. (C) Body weight of young mice. (D, E) Representative NIR fluorescence imaging using a SA‐β‐Gal responsive probe and quantification of senescence‐associated signal intensity in young mice. (F) Frailty index in young mice. (G‐I) Western blot analysis and quantification of P21 and P16 protein levels in spleen from young mice. (J, K) Representative images and quantification of fur depigmentation in 10‐month‐old male mice. (L) Forelimb grip strength in aged mice. (M‐O) Balance and coordination performance in aged mice assessed by balance beam traversal time (M), pause frequency (N), and foot slips (O). (P–T) Hematological parameters in aged mice, including white blood cells (WBC), lymphocytes, granulocytes, monocytes, and red blood cells (RBC).

Consistent with increased senescence burden, comprehensive phenotypic analyses uncovered signs of accelerated functional decline inLINC01021‐expressing mice. Frailty indices increased more steeply in young KI/KI animals relative to WT littermates, while aged transgenic mice displayed a similar upward trend without statistical significance, likely reflecting increased inter‐individual variability (Figures6FandS5E). In parallel, both young and aged KI mice showed a tendency toward reduced whisker number compared with WT controls (FigureS5F,G). At the molecular level, both heterozygous (KI/WT) and homozygous (KI/KI) young mice displayed marked upregulation of the senescence markers P16 and P21 in the spleen (Figure6G–I), indicating that ectopicLINC01021expression enhances senescence‐associated gene expression in vivo, thereby linking organismal phenotypes to underlying cellular senescence processes.

In addition, male transgenic mice exhibited early‐onset coat depigmentation beginning around 8 months of age, with homozygous KI/KI animals developing pronounced whitening by 10 months (Figure6J,K), further supporting premature aging‐associated phenotypes.

Functional assessments further revealed aging‐related impairments in motor performance and systemic immune homeostasis. From 10 months of age onward, transgenic mice exhibited reduced grip strength compared with WT controls, with the most pronounced decline observed in KI/KI mice compared to their WT counterparts, while younger KI mice showed a similar downward trend (Figures6LandS5H). In beam‐walking assays, aged transgenic mice required more time to traverse the beam and paused more frequently, accompanied by a trend toward increased foot slips (Figures6M–OandS5I–K). In contrast, rotarod and treadmill performance did not differ significantly between WT and KI mice (FigureS5L–S), suggesting that motor deficits associated withLINC01021expression are task‐specific rather than globally impaired. Hematological analyses revealed systemic aging‐associated immune alterations, including increased white blood cell and lymphocyte counts in both KI genotypes (Figures6P,QandS5T,U), with additional elevations in granulocytes and monocytes in KI/KI mice (Figures6R,SandS5V,W), while platelet and red blood cell counts remained unchanged (Figures6TandS5X–Z).

Collectively, these findings demonstrate that ectopic expression of humanLINC01021contributes to multiple aging‐like phenotypes in vivo, linking primate‐specificLINC01021activity to organismal senescence and functional decline.

Discussion

Aging is shaped by conserved molecular hallmarks yet unfolds along distinct trajectories across species, implying the existence of lineage‐restricted regulatory layers that modulate canonical aging programs (López‐Otín et al.2023). Here, through comparative evolutionary analyses integrated with cross‐species aging‐associated transcriptomic profiling, we identify a set of primate‐specific lncRNAs that exhibit recurrent age‐dependent expression patterns. Functional interrogation across multiple senescence paradigms and a humanized KI mouse model demonstrates that lineage‐specific lncRNAs can engage conserved senescence machinery and influence aging phenotypes. UsingLINC01021as a representative example, we further delineate a mechanism by which this primate‐specific lncRNA, through interaction with the RNA‐binding protein DAZAP1, destabilizes RBMX, ultimately activating the P53 pathway and translating lineage‐specific regulatory information into measurable cellular and organismal aging outcomes. Together, these findings provide causal evidence that non‐coding regulatory novelty contributes to lineage‐restricted aging phenotypes.

A key conceptual contribution of this study is the establishment of an evolutionary framework for prioritizing aging‐associated lncRNAs, in which lineage restriction is leveraged as an informative biological signal rather than treated as a technical limitation. By systematically integrating primate‐specific conservation with age‐dependent expression across multiple human tissues, we uncovered a substantial subset of lncRNAs exhibiting robust and recurrent associations with aging. The broad tissue distribution of these associations, spanning reproductive, immune, muscular, gastrointestinal, and endocrine systems, suggests that primate‐specific lncRNAs are not confined to tissue‐specific contexts but may participate in higher‐order, systemic aging programs that operate across physiological axes most vulnerable to age‐related decline (Marttila et al.2020; Sherazi et al.2023). Importantly, the absence of orthologs for many of these lncRNAs in rodents highlights an inherent gap in conventional mammalian aging models and raises the possibility that certain dimensions of human aging are governed by regulatory layers that emerged relatively late in evolution (Mattick et al.2023; Yue et al.2014). In this context, primate‐conserved lncRNAs may serve as molecular interfaces through which conserved stress‐response pathways are fine‐tuned in a lineage‐dependent manner, providing a plausible explanation for how shared aging hallmarks give rise to divergent organismal aging trajectories across species (López‐Otín et al.2023; Mattick et al.2023; Sherazi et al.2023).

Extending this evolutionary framework into the context of cellular aging, our findings identifyLINC01021as a representative primate‐specific lncRNA that directly interfaces with the core senescence regulatory network in human fibroblasts. Its consistent upregulation across replicative, IR‐induced, and DOX‐induced senescence models indicates thatLINC01021responds to convergent aging signals rather than to a single category of cellular insult. This pattern is in line with independent reports linkingLINC01021to senescence‐associated states in human cells (Frediani et al.2025; Rossi et al.2023; Wang et al.2025), supporting its relevance as a bona fide regulator of cellular aging. Functionally, elevatedLINC01021expression reinforces growth arrest and amplifies senescence‐associated outputs. Notably, the role ofLINC01021is highly context‐dependent. In proliferative cancer cells, it exerts oncogenic effects primarily by attenuating P53 responses via MDM2‐related pathways, thereby promoting cell cycle progression and tumor development (Kaller et al.2024; Kaller et al.2017; Li et al.2017; Malakar2024). In contrast, in senescent fibroblasts,LINC01021activates rather than suppresses the P53 pathway through the DAZAP1‐RBMX axis, thereby reinforcing cellular senescence. This switch from P53 inhibition in cancer to P53 activation in senescence exemplifies context‐dependent lncRNA regulation, which may arise from differences in cellular state, RNA secondary structure, or protein and epigenetic interactions. Further studies are warranted to elucidate the molecular basis of this switch, which could reflect an evolutionary adaptation in primates for fine‐tuning regulatory pathways according to distinct physiological demands, such as proliferation versus long‐term tissue maintenance.

Mechanistically, our study identifies RBMX as a critical regulatory node through which primate‐specific lncRNAs interface with conserved senescence pathways. RBMX is a multifunctional RNA‐binding protein implicated in RNA processing, genome stability, and transcriptional restraint of P53‐responsive programs (Cai et al.2021; He et al.2024), positioning it as a gatekeeper of stress‐induced growth arrest. We show that age‐associated upregulation ofLINC01021leads to a dramatic reduction in RBMX abundance, thereby releasing RBMX‐mediated constraint on P53 signaling and lowering the threshold for senescence induction. Notably, this regulation does not occur through enhanced proteasomal degradation but instead through accelerated decay ofRBMXmRNA, highlighting post‐transcriptional control as a key interface between lncRNAs and aging pathways (Li et al.2023; Mattick et al.2023). Our data further indicate thatLINC01021disrupts DAZAP1‐mediated stabilization ofRBMXtranscripts, providing a mechanistic link between lncRNA scaffolding functions and RNA‐binding protein networks. Although DAZAP1 itself has not been extensively characterized as an aging regulator, its established role in RNA metabolism places it within a class of RNA‐binding proteins that are increasingly recognized as vulnerable nodes in aging‐associated decline (Choudhury et al.2014; Conlon and Manley2017; Di Fraia et al.2025; Yang et al.2009). This RNA‐centric and modular mode of regulation may therefore represent a common strategy through which evolutionary innovation shapes lineage‐specific aging phenotypes.

Importantly, the regulatory logic uncovered here aligns with a broader emerging framework in which lncRNAs act as integrative modulators of aging across cellular and tissue contexts. While many aging‐associated lncRNAs are evolutionarily conserved, comparing their regulatory magnitude and functional scope with primate‐specific lncRNAs provides important context for assessing lineage‐specific regulatory innovation. In skeletal muscle, age‐associated decline of the conserved lncRNAMALAT1enhances TGF‐β1 signaling and disrupts myoblast proliferation and cell‐cycle progression through P53‐dependent mechanisms, contributing to impaired regenerative capacity (Ruan et al.2021). Conversely, in vascular smooth muscle cells undergoing replicative aging,NEAT1is robustly upregulated and reinforces senescence‐associated secretory programs, including IL8 expression, via P53‐centered feedback loops (Rossi et al.2023). Notably, within these aging‐relevant contexts, the primate‐specificLINC01021exhibits even stronger induction than conserved lncRNAs, highlighting the potential for lineage‐restricted regulatory layers to exert disproportionate influence on aging pathways. Beyond cellular models, dysregulation of circulating lncRNAs such asH19,NEAT1, andPVT1correlates with sarcopenia, frailty, cardiovascular dysfunction, and neurodegenerative progression, suggesting that lncRNA‐mediated regulatory shifts are reflected at the organismal level (Aparicio et al.2025). Together, these observations position lncRNAs as molecular conduits linking cell‐intrinsic senescence programs to tissue dysfunction and age‐related disease phenotypes, providing a broader context for understanding how primate‐specific lncRNAs such asLINC01021may shape human aging trajectories.

In a broader evolutionary context, lineage‐restricted regulatory innovations have increasingly been recognized as important modulators of species‐specific aging trajectories. Rather than altering conserved core pathways, such innovations appear to fine‐tune stress response, genome maintenance, and inflammatory programs in a lineage‐dependent manner (Firsanov et al.2025; Huang et al.2020; Jiang and Kong2020; Seluanov et al.2018). Our findings extend this framework by demonstrating that primate‐specific lncRNAs can interface directly with conserved senescence effectors, thereby providing a mechanistic link between evolutionary novelty and canonical aging pathways. Within this perspective,LINC01021exemplifies recently emerged non‐coding regulators that recalibrate the threshold and dynamics of P53 activation to generate species‐restricted aging phenotypes.

LINC01021acts as an aging accelerator by lowering the threshold for senescence induction. In the humanized KI mice, constitutive expression ofLINC01021from the Rosa26 locus under the CAG promoter resulted in increased frailty, depigmentation, and motor coordination deficits, indicating that this primate‐specific lncRNA modulates aging‐related processes in vivo. Immune remodeling emerged more clearly in aged cohorts, and differences in older mice were often trends rather than robust. This pattern supports a model whereLINC01021advances the onset of functional decline by sensitizing tissues to senescence signals (Di Micco et al.2021; Sun et al.2022). Behavioral impairments and immune changes resembling human immunosenescence indicate preferential impact on specific systems (Huang et al.2021; Sherazi et al.2023; Zayoud et al.2024). Partial penetrance and age‐dependent expressivity reflect cross‐species differences (Mattick et al.2023). Notably, these phenotypes were relatively modest compared with the pronounced senescence responses observed in cultured human fibroblasts. This difference may reflect the gradual accumulation of senescent cells during aging, tissue‐ and organism‐level compensatory mechanisms, and the challenges of modeling a primate‐specific lncRNA within a murine genetic background (Karin and Alon2021; Ogrodnik2021; Pal and Tyler2016). In particular, transgene insertion at the Rosa26 locus may not fully recapitulate the native regulatory environment of the humanLINC01021locus (Chen et al.2011), while divergence in downstream regulatory networks could further attenuate its biological effects (Ghanam et al.2022; Mattick et al.2023). Consistent with this interpretation, phenotypic alterations were most evident in young adults and became less pronounced in aged cohorts, suggesting thatLINC01021may accelerate the onset of functional decline rather than drive severe aging phenotypes. Mechanistically, our cellular studies support a model wherebyLINC01021promotes senescence through the DAZAP1‐RBMX‐P53 pathway. Although direct validation of this mechanism in vivo remains an important direction for future investigation, both DAZAP1 and RBMX have established roles in regulating cell proliferation, development, and tissue homeostasis. DAZAP1‐deficient mice display growth retardation, reproductive defects, and impaired tissue development, indicating an essential role in cellular growth and survival (Hsu et al.2008). Similarly, RBMX has been implicated in RNA processing, genome stability, and proliferative control, and its disruption leads to defects in development and tissue maintenance (Qiu, Pu, et al.2025). These observations are consistent with the cellular functions affected byLINC01021and suggest that components of the DAZAP1‐RBMX pathway likely contribute to the aging‐like phenotypes observed inLINC01021KI mice.

Despite these advances, several considerations highlight the boundaries and future directions of this work. First, the constitutive expression of a human lncRNA in a species lacking endogenous orthologs may amplify certain phenotypes, underscoring the need for tissue‐specific and inducible models to disentangle developmental effects from aging‐associated processes. Second, although our data support a DAZAP1‐RBMX‐P53 regulatory axis, the precise molecular interfaces governing these interactions remain to be defined. WhetherLINC01021directly scaffolds RNA‐binding proteins or instead reshapes post‐transcriptional regulatory complexes through indirect mechanisms will require detailed structural and biochemical interrogation. Moreover, whileLINC01021was selected as a representative exemplar, our evolutionary screening uncovered multiple primate‐specific lncRNAs with robust age‐dependent expression patterns. Determining whether these candidates converge on shared senescence pathways or engage distinct regulatory circuits will be essential for delineating the broader architecture of lineage‐restricted aging control.

Altogether, this study proposes that primate‐specific lncRNAs constitute an evolutionarily encoded regulatory layer that interfaces with conserved senescence machinery and may contribute to species‐restricted aging trajectories. By demonstrating thatLINC01021modulates RBMX stability to influence P53 signaling, thereby promoting cellular senescence and aging‐like phenotypes, our work provides causal evidence that non‐coding regulatory novelty can actively reprogram core aging pathways. More broadly, these findings caution against exclusive reliance on rodent models for dissecting aging mechanisms and highlight the necessity of incorporating lineage‐specific regulators into aging research frameworks. Elucidating such primate‐specific regulatory axes may ultimately enable more precise strategies for modulating senescence and age‐related pathologies within the molecular context of human evolution.

Author Contributions

Conceptualization: Y.H. Data curation: Y.Z., L.H., and X.D. Formal analysis: Y.Z., X.D., M.Z., C.L., M.G., G.L., and S.P. Funding acquisition: X.D., R.B., G.L., and Y.H. Methodology: Y.Z., L.H., and X.D. Project administration: Y.H., and G.L. Supervision: Y.H. Validation: Y.Z. and L.H. Visualization: Y.Z. and X.D. Writing – original draft: Y.Z. and X.D. Writing – review and editing: A.N., S.K., and Y.H. Investigation: Y.Z., L.H., X.D., and Q.Z.

Funding

This work was supported by the National Key R&D Program of China (2023YFC3603300 to YHH), National Natural Science Foundation of China (82471599 to YHH, 32500660 to RXB), Yunnan Fundamental Research Projects (202602AS100006 to YHH, 202305AH340006 to YHH, 202401CF070064 to XD), CAS “Light of West China” Program (xbzg‐zdsys‐202312 to YHH), The Postdoctoral Fellowship Program of China Postdoctoral Science Foundation (CPSF) (GZC20232763 to XD), State Key Laboratory Conservation and Utilization of Bio‐Resources in Yunnan (2023KF008 to XD), China Postdoctoral Science Foundation (2025M772550 to RXB), and Yunnan Revitalization Talent Support Program Young Talent Project (GHL). YHH is supported by the Pioneer Hundred Talents Program of the Chinese Academy of Sciences. In addition, YHH is supported by the Yunnan Revitalization Talent Support Program Young Talent Project.

Ethics Statement

All animal procedures were reviewed and approved by the Animal Ethics Committee of the Kunming Institute of Zoology, Chinese Academy of Sciences (IACUC‐RE‐2025‐08‐008), and were carried out in compliance with institutional guidelines for laboratory animal care and use.

Consent

The authors have nothing to report.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Aparicio, P. , T. López‐Royo, D. Navarrete‐Villanueva, et al. 2025. “Serum lncRNAs NEAT1, PVT1 and H19 as Novel Biomarkers for Sarcopenia Diagnosis and Treatment Response. ”Non‐coding RNA Research14: 166–176. . doi.org/10.1016/j.ncrna.2025.06.003
  2. Breschi, A. , T. R. Gingeras, andR. Guigó. 2017. “Comparative Transcriptomics in Human and Mouse. ”Nature Reviews Genetics18, no. 7: 425–440. . doi.org/10.1038/nrg.2017.19
  3. Cagan, A. , A. Baez‐Ortega, N. Brzozowska, et al. 2022. “Somatic Mutation Rates Scale With Lifespan Across Mammals. ”Nature604, no. 7906: 517–524. . doi.org/10.1038/s41586-022-04618-z
  4. Cai, T. , J. K. Cinkornpumin, Z. Yu, O. D. Villarreal, W. A. Pastor, andS. Richard. 2021. “Deletion of RBMX RGG/RG Motif in Shashi‐XLID Syndrome Leads to Aberrant p53 Activation and Neuronal Differentiation Defects. ”Cell Reports36, no. 2: 109337. . doi.org/10.1016/j.celrep.2021.109337
  5. Chen, C. M. , J. Krohn, S. Bhattacharya, andB. Davies. 2011. “A Comparison of Exogenous Promoter Activity at the ROSA26 Locus Using a ΦiC31 Integrase Mediated Cassette Exchange Approach in Mouse ES Cells. ”PLoS One6, no. 8: e23376. . doi.org/10.1371/journal.pone.0023376
  6. Choudhury, R. , S. G. Roy, Y. S. Tsai, A. Tripathy, L. M. Graves, andZ. Wang. 2014. “The Splicing Activator DAZAP1 Integrates Splicing Control Into MEK/Erk‐Regulated Cell Proliferation and Migration. ”Nature Communications5: 3078. . doi.org/10.1038/ncomms4078
  7. Conlon, E. G. , andJ. L. Manley. 2017. “RNA‐Binding Proteins in Neurodegeneration: Mechanisms in Aggregate. ”Genes & Development31, no. 15: 1509–1528. . doi.org/10.1101/gad.304055.117
  8. Consortium, G. 2013. “The Genotype‐Tissue Expression (GTEx) Project. ”Nature Genetics45, no. 6: 580–585. . doi.org/10.1038/ng.2653
  9. Di Fraia, D. , A. Marino, J. H. Lee, et al. 2025. “Altered Translation Elongation Contributes to Key Hallmarks of Aging in the Killifish Brain. ”Science389, no. 6759: eadk3079. . doi.org/10.1126/science.adk3079
  10. Di Micco, R. , V. Krizhanovsky, D. Baker, andF. d'Adda di Fagagna. 2021. “Cellular Senescence in Ageing: From Mechanisms to Therapeutic Opportunities. ”Nature Reviews. Molecular Cell Biology22, no. 2: 75–95. . doi.org/10.1038/s41580-020-00314-w
  11. Drekolia, M. K. , S. Talyan, R. Cordellini Emídio, et al. 2022. “Unravelling the Impact of Aging on the Human Endothelial lncRNA Transcriptome. ”Frontiers in Genetics13: 1035380. . doi.org/10.3389/fgene.2022.1035380
  12. Feridooni, H. A. , M. H. Sun, K. Rockwood, andS. E. Howlett. 2015. “Reliability of a Frailty Index Based on the Clinical Assessment of Health Deficits in Male C57BL/6J Mice. ”Journals of Gerontology. Series A, Biological Sciences and Medical Sciences70, no. 6: 686–693. . doi.org/10.1093/gerona/glu161
  13. Field, A. R. , F. M. J. Jacobs, I. T. Fiddes, et al. 2019. “Structurally Conserved Primate lncRNAs Are Transiently Expressed During Human Cortical Differentiation and Influence Cell‐Type‐Specific Genes. ”Stem Cell Reports12, no. 2: 245–257. . doi.org/10.1016/j.stemcr.2018.12.006
  14. Firsanov, D. , M. Zacher, X. Tian, et al. 2025. “Evidence for Improved DNA Repair in the Long‐Lived Bowhead Whale. ”Nature648, no. 8094: 717–725. . doi.org/10.1038/s41586-025-09694-5
  15. Frediani, E. , C. Anceschi, J. Ruzzolini, et al. 2025. “Divergent Regulation of Long Non‐Coding RNAs H19 and PURPL Affects Cell Senescence in Human Dermal Fibroblasts. ”Geroscience47, no. 2: 2079–2097. . doi.org/10.1007/s11357-024-01399-3
  16. Ghanam, A. R. , W. B. Bryant, andJ. A. ‐O. Miano. 2022. “Of Mice and Human‐Specific Long Noncoding RNAs. ”Mammalian Genome: Official Journal of the International Mammalian Genome Society33, no. 2: 281–292. . doi.org/10.1007/s00335-022-09943-2
  17. Guo, J. , X. Huang, L. Dou, et al. 2022. “Aging and Aging‐Related Diseases: From Molecular Mechanisms to Interventions and Treatments. ”Signal Transduction and Targeted Therapy7, no. 1: 391. . doi.org/10.1038/s41392-022-01251-0
  18. Han, S. , X. Li, K. Wang, et al. 2021. “PURPL Represses Autophagic Cell Death to Promote Cutaneous Melanoma by Modulating ULK1 Phosphorylation. ”Cell Death & Disease12, no. 11: 1070. . doi.org/10.1038/s41419-021-04362-8
  19. Hasegawa, Y. , Y. Daitoku, K. Sekiguchi, et al. 2013. “Novel ROSA26 Cre‐Reporter Knock‐In C57BL/6N Mice Exhibiting Green Emission Before and Red Emission After Cre‐Mediated Recombination. ”Experimental Animals62, no. 4: 295–304. . doi.org/10.1538/expanim.62.295
  20. He, D. , X. He, D. Shen, et al. 2024. “Loss‐Of‐Function Variants in RNA Binding Motif Protein X‐Linked Induce Neuronal Defects Contributing to Amyotrophic Lateral Sclerosis Pathogenesis. ”MedComm5, no. 9: e712. . doi.org/10.1002/mco2.712
  21. He, Y. , X. Zhang, J. Chang, et al. 2020. “Using Proteolysis‐Targeting Chimera Technology to Reduce Navitoclax Platelet Toxicity and Improve Its Senolytic Activity. ”Nature Communications11, no. 1: 1996. . doi.org/10.1038/s41467-020-15838-0
  22. Hsu, L. C. L. , H. Y. Chen, Y. W. Lin, et al. 2008. “DAZAP1, an hnRNP Protein, Is Required for Normal Growth and Spermatogenesis in Mice. ”RNA14, no. 9: 1814–1822. . doi.org/10.1261/rna.1152808
  23. Hu, L. , C. Dong, Z. Wang, et al. 2023. “A Rationally Designed Fluorescence Probe Achieves Highly Specific and Long‐Term Detection of Senescence in Vitro and in Vivo. ”Aging Cell22, no. 8: e13896. . doi.org/10.1111/acel.13896
  24. Huang, Z. , B. Chen, X. Liu, et al. 2021. “Effects of Sex and Aging on the Immune Cell Landscape as Assessed by Single‐Cell Transcriptomic Analysis. ”Proceedings of the National Academy of Sciences of the United States of America118, no. 33: e2023216118. . doi.org/10.1073/pnas.2023216118
  25. Huang, Z. , C. V. Whelan, D. Dechmann, andE. C. Teeling. 2020. “Genetic Variation Between Long‐Lived Versus Short‐Lived Bats Illuminates the Molecular Signatures of Longevity. ”Aging (Albany NY)12, no. 16: 15962–15977. . doi.org/10.18632/aging.103725
  26. Jiang, J. J. , andQ. P. Kong. 2020. “Comparative Analysis of Long Noncoding RNAs in Long‐Lived Mammals Provides Insights Into Natural Cancer‐Resistance. ”RNA Biology17, no. 11: 1657–1665. . doi.org/10.1080/15476286.2020.1792116
  27. Kaller, M. , I. Forné, A. Imhof, andH. Hermeking. 2024. “LINC01021 Attenuates Expression and Affects Alternative Splicing of a Subset of p53‐Regulated Genes. ”Cancers16, no. 9: 1639. . doi.org/10.3390/cancers16091639
  28. Kaller, M. , U. Götz, andH. Hermeking. 2017. “Loss of p53‐Inducible Long Non‐Coding RNA LINC01021 Increases Chemosensitivity. ”Oncotarget8, no. 61: 102783–102800. . doi.org/10.18632/oncotarget.22245
  29. Karin, O. , andU. Alon. 2021. “Senescent Cell Accumulation Mechanisms Inferred From Parabiosis. ”Geroscience43, no. 1: 329–341. . doi.org/10.1007/s11357-020-00286-x
  30. Keane, M. , J. Semeiks, A. E. Webb, et al. 2015. “Insights Into the Evolution of Longevity From the Bowhead Whale Genome. ”Cell Reports10, no. 1: 112–122. . doi.org/10.1016/j.celrep.2014.12.008
  31. Kopp, F. , andJ. T. Mendell. 2018. “Functional Classification and Experimental Dissection of Long Noncoding RNAs. ”Cell172, no. 3: 393–407. . doi.org/10.1016/j.cell.2018.01.011
  32. Leng, X. , M. Zhang, Y. Xu, et al. 2024. “Non‐Coding RNAs as Therapeutic Targets in Cancer and Its Clinical Application. ”Journal of Pharmaceutical Analysis14, no. 7: 100947. . doi.org/10.1016/j.jpha.2024.02.001
  33. Li, X. , Q. Liu, andJ. Liu. 2023. “Long Non‐Coding RNAs: Discoveries, Mechanisms, and Research Strategies in Seeds. ”Genes14, no. 12: 2214. . doi.org/10.3390/genes14122214
  34. Li, X. L. , M. Subramanian, M. F. Jones, et al. 2017. “Long Noncoding RNA PURPL Suppresses Basal p53 Levels and Promotes Tumorigenicity in Colorectal Cancer. ”Cell Reports20, no. 10: 2408–2423. . doi.org/10.1016/j.celrep.2017.08.041
  35. López‐Otín, C. , M. A. Blasco, L. Partridge, M. Serrano, andG. Kroemer. 2023. “Hallmarks of Aging: An Expanding Universe. ”Cell186, no. 2: 243–278. . doi.org/10.1016/j.cell.2022.11.001
  36. Malakar, P. 2024. “Long Noncoding RNA, PURPL Is Associated With Aneuploidy and Its Magnitude of Expression Level Is Dependent on P53 Status. ”Frontiers in Cell and Developmental Biology12: 1410308. . doi.org/10.3389/fcell.2024.1410308
  37. Marthandan, S. , M. Baumgart, S. Priebe, et al. 2016. “Conserved Senescence Associated Genes and Pathways in Primary Human Fibroblasts Detected by RNA‐Seq. ”PLoS One11, no. 5: e0154531. . doi.org/10.1371/journal.pone.0154531
  38. Marttila, S. , K. Chatsirisupachai, D. Palmer, andJ. P. de Magalhães. 2020. “Ageing‐Associated Changes in the Expression of lncRNAs in Human Tissues Reflect a Transcriptional Modulation in Ageing Pathways. ”Mechanisms of Ageing and Development185: 111177. . doi.org/10.1016/j.mad.2019.111177
  39. Mattick, J. S. , P. P. Amaral, P. Carninci, et al. 2023. “Long Non‐Coding RNAs: Definitions, Functions, Challenges and Recommendations. ”Nature Reviews. Molecular Cell Biology24, no. 6: 430–447. . doi.org/10.1038/s41580-022-00566-8
  40. Mercer, T. R. , M. E. Dinger, andJ. S. Mattick. 2009. “Long Non‐Coding RNAs: Insights Into Functions. ”Nature Reviews. Genetics10, no. 3: 155–159. . doi.org/10.1038/nrg2521
  41. Mills, J. D. , M. Ward, B. J. Chen, A. M. Iyer, E. Aronica, andM. Janitz. 2016. “LINC00507 Is Specifically Expressed in the Primate Cortex and Has Age‐Dependent Expression Patterns. ”Journal of Molecular Neuroscience: MN59, no. 4: 431–439. . doi.org/10.1007/s12031-016-0745-4
  42. Ogrodnik, M. A. ‐O. 2021. “Cellular Aging Beyond Cellular Senescence: Markers of Senescence Prior to Cell Cycle Arrest in Vitro and in Vivo. ”Aging Cell20, no. 4: e13338. . doi.org/10.1111/acel.13338
  43. Pal, S. , andJ. K. Tyler. 2016. “Epigenetics and Aging. ”Science Advances2, no. 7: e1600584. . doi.org/10.1126/sciadv.1600584
  44. Pollen, A. A. , U. Kilik, C. B. Lowe, andJ. G. Camp. 2023. “Human‐Specific Genetics: New Tools to Explore the Molecular and Cellular Basis of Human Evolution. ”Nature Reviews Genetics24, no. 10: 687–711. . doi.org/10.1038/s41576-022-00568-4
  45. Qiu, L. , W. Li, L. Zhang, et al. 2025. “p52‐ZER6/DAZAP1 Axis Promotes Ferroptosis Resistance and Colorectal Cancer Progression via Regulating SLC7A11 mRNA Stabilization. ”Acta Pharmaceutica Sinica B15, no. 4: 2039–2058. . doi.org/10.1016/j.apsb.2025.02.013
  46. Qiu, Y. , C. Pu, C. Wang, andZ. Quan. 2025. “The RNA‐Binding Protein RBMX Mediates the Immunosuppressive Microenvironment of Osteosarcoma by Regulating CD8(+)T Cells. ”Cancers17, no. 17: 2928. . doi.org/10.3390/cancers17172928
  47. Rigby Dames, B. A. , H. Kilili, C. J. Charvet, et al. 2023. “Evolutionary and Genomic Perspectives of Brain Aging and Neurodegenerative Diseases. ”Progress in Brain Research275: 165–215. . doi.org/10.1016/bs.pbr.2022.10.004
  48. Rossi, C. , M. Venturin, J. Gubala, et al. 2023. “PURPL and NEAT1 Long Non‐Coding RNAs Are Modulated in Vascular Smooth Muscle Cell Replicative Senescence. ”Biomedicine11, no. 12: 3228. . doi.org/10.3390/biomedicines11123228
  49. Ruan, L. , B. Mendhe, E. Parker, et al. 2021. “Long Non‐Coding RNA MALAT1 Is Depleted With Age in Skeletal Muscle in Vivo and MALAT1 Silencing Increases Expression of TGF‐β1 in Vitro. ”Frontiers in Physiology12: 742004. . doi.org/10.3389/fphys.2021.742004
  50. Seluanov, A. , V. N. Gladyshev, J. Vijg, andV. Gorbunova. 2018. “Mechanisms of Cancer Resistance in Long‐Lived Mammals. ”Nature Reviews. Cancer18, no. 7: 433–441. . doi.org/10.1038/s41568-018-0004-9
  51. Sherazi, S. A. M. , A. Abbasi, A. Jamil, et al. 2023. “Molecular Hallmarks of Long Non‐Coding RNAs in Aging and Its Significant Effect on Aging‐Associated Diseases. ”Neural Regeneration Research18, no. 5: 959–968. . doi.org/10.4103/1673-5374.355751
  52. Statello, L. , C. J. Guo, L. L. Chen, andM. Huarte. 2021. “Gene Regulation by Long Non‐Coding RNAs and Its Biological Functions. ”Nature Reviews. Molecular Cell Biology22, no. 2: 96–118. . doi.org/10.1038/s41580-020-00315-9
  53. Sun, Y. , Q. Li, andJ. L. Kirkland. 2022. “Targeting Senescent Cells for a Healthier Longevity: The Roadmap for an Era of Global Aging. ”Lifestyle Medicine1, no. 2: 103–119. . doi.org/10.1093/lifemedi/lnac030
  54. Turko, R. , A. Hajja, A. M. Magableh, et al. 2025. “The Emerging Role of miRNAs in Biological Aging and Age‐Related Diseases. ”Non‐coding RNA Research13, no. 131: 131–152. . doi.org/10.1016/j.ncrna.2025.05.002
  55. Vickery, S. , K. R. Patil, R. Dahnke, et al. 2024. “The Uniqueness of Human Vulnerability to Brain Aging in Great Ape Evolution. ”Science Advances10, no. 35: eado2733. . doi.org/10.1126/sciadv.ado2733
  56. Wagner, V. , F. Kern, O. Hahn, et al. 2024. “Characterizing Expression Changes in Noncoding RNAs During Aging and Heterochronic Parabiosis Across Mouse Tissues. ”Nature Biotechnology42, no. 1: 109–118. . doi.org/10.1038/s41587-023-01751-6
  57. Wang, A. 2025. “Noncoding RNAs Evolutionarily Extend Animal Lifespan. ”Global Medical Genetics12, no. 2: 100034. . doi.org/10.1016/j.gmg.2024.100034
  58. Wang, J. , X. Yang, X. Su, et al. 2025. “Targeting PURPL RNA Enabled Rejuvenation of Senescence Cells via Epigenetic Reprogramming. ”Journal of Translational Medicine23, no. 1: 1127. . doi.org/10.1186/s12967-025-07208-5
  59. Wang, K. , H. Liu, Q. Hu, et al. 2022. “Epigenetic Regulation of Aging: Implications for Interventions of Aging and Diseases. ”Signal Transduction and Targeted Therapy7, no. 1: 374. . doi.org/10.1038/s41392-022-01211-8
  60. Wang, P. , Y. Xue, Y. Han, et al. 2014. “The STAT3‐Binding Long Noncoding RNA Lnc‐DC Controls Human Dendritic Cell Differentiation. ”Science344, no. 6181: 310–313. . doi.org/10.1126/science.1251456
  61. Wang, Q. , Y. Guo, W. Wang, et al. 2021. “RNA Binding Protein DAZAP1 Promotes HCC Progression and Regulates Ferroptosis by Interacting With SLC7A11 mRNA. ”Experimental Cell Research399, no. 1: 112453. . doi.org/10.1016/j.yexcr.2020.112453
  62. Xiao, F. ‐H. , Q. Yu, Z. ‐L. Deng, et al. 2022. “ETS1 Acts as a Regulator of Human Healthy Aging via Decreasing Ribosomal Activity. ”Science Advances8, no. 17: eabf2017. . doi.org/10.1126/sciadv.abf2017
  63. Yan, Y. , H. Yan, Y. Teng, et al. 2020. “Long Non‐Coding RNA 00507/miRNA‐181c‐5p/TTBK1/MAPT Axis Regulates Tau Hyperphosphorylation in Alzheimer's Disease. ”Journal of Gene Medicine22, no. 12: e3268. . doi.org/10.1002/jgm.3268
  64. Yang, H. T. , M. Peggie, P. Cohen, andS. Rousseau. 2009. “DAZAP1 Interacts via Its RNA‐Recognition Motifs With the C‐Termini of Other RNA‐Binding Proteins. ”Biochemical and Biophysical Research Communications380, no. 3: 705–709. . doi.org/10.1016/j.bbrc.2009.01.166
  65. Yao, R. W. , Y. Wang, andL. L. Chen. 2019. “Cellular Functions of Long Noncoding RNAs. ”Nature Cell Biology21, no. 5: 542–551. . doi.org/10.1038/s41556-019-0311-8
  66. Yu, Z. , I. Seim, M. Yin, et al. 2021. “Comparative Analyses of Aging‐Related Genes in Long‐Lived Mammals Provide Insights Into Natural Longevity. ”Innovation (Camb)2, no. 2: 100108. . doi.org/10.1016/j.xinn.2021.100108
  67. Yue, F. , Y. Cheng, A. Breschi, et al. 2014. “A Comparative Encyclopedia of DNA Elements in the Mouse Genome. ”Nature515, no. 7527: 355–364. . doi.org/10.1038/nature13992
  68. Zayoud, K. , A. Chikhaoui, I. Kraoua, et al. 2024. “Immunity in the Progeroid Model of Cockayne Syndrome: Biomarkers of Pathological Aging. ”Cells13, no. 5: 402. . doi.org/10.3390/cells13050402
  69. Zeng, X. , J. Cao, J. Xu, et al. 2025. “SKP1‐CUL1‐F‐Box: Key Molecular Targets Affecting Disease Progression. ”FASEB Journal: Official Publication of the Federation of American Societies for Experimental Biology39, no. 2: e70326. . doi.org/10.1096/fj.202402816RR
  70. Zhang, J. , Z. Ni, Y. Zhang, et al. 2024. “DAZAP1 Phase Separation Regulates Mitochondrial Metabolism to Facilitate Invasion and Metastasis of Oral Squamous Cell Carcinoma. ”Cancer Research84, no. 22: 3818–3833. . doi.org/10.1158/0008-5472.Can-24-0067
  71. Zhang, Q. , X. Li, K. Cui, et al. 2020. “The MAP3K13‐TRIM25‐FBXW7α Axis Affects c‐Myc Protein Stability and Tumor Development. ”Cell Death and Differentiation27, no. 2: 420–433. . doi.org/10.1038/s41418-019-0363-0
  72. Zhao, L. , J. Wang, Y. Li, et al. 2021. “NONCODEV6: An Updated Database Dedicated to Long Non‐Coding RNA Annotation in Both Animals and Plants. ”Nucleic Acids Research49, no. D1: D165–d171. . doi.org/10.1093/nar/gkaa1046

Republished from the open web under CC-BY. Authors: Zhang Y, Hu L, Dong X, Zeng Q, Zi M, Nisar A, Khan S, Bai R, Liu C, Ge M, Pu S, Li G, He Y. Read the original.