LivAge：一种用于小鼠转录组学年龄估计的在线衰老时钟

《Aging Cell》：LivAge: An Online Aging Clock for Murine Transcriptomic Age Estimation

　　

过去一个世纪预期寿命的提升伴随着将年龄识别为包括心血管与神经退行性疾病以及癌症在内广泛病症的主要风险因子。因此，开发能够延缓潜在生物学衰老进程的干预手段，可能对预防或减缓这些疾病具有益处。为评估此类老年保护（geroprotective）干预，通过衰老时钟——尤其是基于转录组生物标志物的时钟——量化生物学损伤已变得高度相关，因其提供具备强生物学可解释性的估计。然而，面向小鼠实验模型的转录组时钟可用性有限，阻碍了这些工具在衰老研究中的实施及对其老年保护效应的评估。研究人员开发了一种准确、可及且即开即用的小鼠转录组时钟，利用肝脏RNA-seq数据为健康小鼠提供稳健的年龄预测。将该时钟应用于加速衰老模型揭示出早衰（progeroid）综合征中转录组年龄升高，证明其捕捉衰老相关生物学进程的能力。最后，鉴于此类工具的主要旨趣，研究人员显示具备老年保护潜力的既有遗传干预（Snell Dwarf、Ames Dwarf及生长激素受体缺陷小鼠）与环境干预（热量、蛋白质或蛋氨酸限制，以及特定情境下的雷帕霉素补充）均降低了该时钟估计的转录组年龄。这些发现表明所提转录组时钟可成为研究衰老生物学及设计评估潜在老年保护干预的宝贵工具。

研究背景与立项依据

衰老被定义为生命累积后果（分子与细胞损伤）导致功能衰退、慢性病及最终死亡风险升高的过程。老龄化推动癌症、心血管疾病与神经退行性疾病负担上升，故延缓衰老标志（hallmarks of aging）而非单病靶点的老年保护策略备受关注。现有衰老时钟以DNA甲基化（epigenetic clocks）为主，但甲基化漂移的因果可解释性弱；人源转录组时钟虽有进展，小鼠（Mus musculus）实验模型专用转录组时钟稀缺，且既有工具面临性能不稳、缺乏外部验证与友好软件之困。为此，研究人员在《Aging Cell》发表LivAge——一种基于C57BL/6小鼠肝脏bulk RNA-seq的在线转录组衰老时钟，用以精准估计转录组年龄（transcriptomic age）并量化老年保护干预效应。

主要技术方法概览

研究人员整合23项独立研究的432例C57BL/6小鼠肝脏RNA-seq样本（1–30月龄，365雄/67雌）构建衰老转录组图谱；采用nf-core RNA-seq流程（STAR+Salmon，GRCm39）统一比对定量，DESeq2以“~age+sex+study”做差异表达分析（likelihood ratio test，FDR<0.05，|log₂FC/月|>0.02），fgsea做Hallmark通路富集，MuSiC基于Tabula Muris做细胞反卷积（deconvolution）。预测端以glmnet弹性网络（Elastic Net）为骨架，经重复5折×10次交叉验证调参，对比Box-Cox（年龄）与YuGene（表达）等数值变换提升泛化；独立外部集（39例→134例）逐级验证，并与TACO、RAPToR基准比对。干预评估以AgeDev（估计龄?实足龄）为因变量做线性回归求β系数，Benjamini-Hochberg校正FDR；平台用Shiny封装，供原始count矩阵直出CPM→YuGene→逆Box-Cox月龄。

研究结果

2.1 小鼠肝脏衰老转录组图谱生成

通过对432例样本差异表达分析获3350个年龄相关基因：上调含Gzmk、Ccl8、Cidea，下调含Mup家族、Ces2b、Ces4a、Hsd3b5。GSEA显示免疫炎症通路（Allograft Rejection NES=1.90；Interferon Gamma Response NES=1.83）正富集，氧化磷酸化（NES=?2.90）、脂肪酸代谢、成脂负富集；细胞反卷积见肝细胞比例随龄下降（斜率?0.440，p=4.33×10^?10），Kupffer（0.246）与B细胞（0.086）上升，印证“inflammaging”。

2.2 小鼠转录组衰老时钟开发

由差异基因初筛（FDR<0.25）入模，Elastic Net交叉验证后，Box-Cox转年龄+YuGene转表达组合最优。第三外部集134例（雄50/雌84）得雄R²=0.930、Pearson r=0.964、MAE=1.53月；雌R²=0.918、MAE=2.47月；显著优于同期TACO与RAPToR。最终模型倚重268个预测基因；ILR转细胞比例仅解释估计龄方差35.2%，说明时钟捕获超越细胞组成改变的复杂转录重编程。

2.3 早衰小鼠模型呈现加速转录组年龄

将时钟用于Zmpste24^?/?（β=0.881，FDR=0.004）、Lmna^G609G/G609G（β=0.943，FDR=0.004）、Terc^?/?（β=1.296，FDR=0.008）、Sprtn^Y118C/Y118C（β=1.524，FDR=0.008）肝病样/基因组不稳定早衰模型，AgeDev均显著为正，证时钟反映生物学加速而非实足时间。

2.4 部分遗传长寿干预降低转录组年龄

Ames Dwarf（Prop1缺陷，β=?4.101，FDR=0.004）、Snell Dwarf（Pit1缺陷，β=?2.519，FDR<0.001）及他莫昔芬诱导型全局GHRKO（12月诱导至24月，β=?2.412，FDR=0.006）呈显著负AgeDev；组成型GHRKO仅趋势不显著（β=?1.337，FDR=0.114），提示品系/诱导时机影响。

2.5 饮食与药物策略减速估计年龄

8月40%热量限制（caloric restriction，CR）趋显著（β=?1.428，FDR=0.125），2月CR无效；低脂 vs 高脂（high-fat，HF）β=?3.224（FDR=0.029），标准料 vs HF β=?1.562；HF下CR再降β=?3.917，低fat（low-fat，LF）下CR β=?3.066；蛋白限制β=?0.819（FDR=0.015），蛋氨酸限制（methionine restriction，MR）β=?3.917（FDR=0.004）。药理端，HF背景下雷帕霉素2月β=?4.574（FDR<0.001）；UM-HET3群中二甲双胍、8月雷帕霉素仅非显著趋势，余化合物无效应。

讨论部分总结

讨论指出，LivAge以肝脏为窗，克服表观时钟因果模糊、小鼠转录时钟荒缺之弊；268基因签名兼顾代谢下行与炎症上行，且细胞组成仅释35%方差，余为实质转录程序。早衰四类模型全阳、遗传/营养干预多阴向AgeDev，与寿命操纵已知表型一致；但组成型GHRKO不显著、UM-HET3药理响应弱，归因为非C57BL/6背景、亚株差异、batch效应及药物作用过窄未撼动全局转录态。局限含：雌样本少致MAE偏高；bulk RNA-seq不能纵向（肝取样致死）；部分预测基因功能未明；绝对龄受技术噪声影响。作者主张单细转录时钟、蛋白/代谢多组学融合为补充方向。

结论部分翻译

本研究开发人员提出LivAge，一种准确、可及且即开即用的转录组时钟，基于肝脏RNA-seq提供稳健转录组年龄估计。汇432例1–30月小鼠肝样，融合差异表达分析、年龄Box-Cox变换与表达YuGene归一化训练Elastic Net，外部基准达金标准表观时钟量级。LivAge在早衰小鼠捕获加速转录衰老，反之在遗传与环境老年保护干预下检出显著年龄减速。为便于抗衰疗法评估与社区采用，配套Shiny网页端（https://gencanage.uniovi.es/resources/liv_age ），实现标准化、可重复、用户友好的转录组年龄评定。