损伤累积模型识别跨物种的不同衰老模式

《Nature Aging》：A damage accumulation model identifies distinct aging regimes across species

　　

不同物种以相似的方式衰老，但其寿命差异可达数个数量级。研究人员尚不清楚这些相似性与差异性如何从构成衰老基础的损伤累积中产生。长寿究竟源于损伤产生减少、清除增强还是对损伤的抵抗能力增强？研究人员在此应用饱和清除（saturating removal, SR）模型——一种损伤累积与清除的随机模型——并将其拟合到多个经过充分研究的物种的生存数据。研究发现，多个参数具有近乎普适的数值，包括清除速率比、噪声幅度和死亡阈值。其中，最能预测寿命的模型参数是损伤产生率（η），其跨度达七个数量级。研究人员识别出两种截然不同的衰老模式：弹道式衰老（ballistic aging），即损伤产生超过清除，以酵母、线虫、果蝇和小鼠为特征；以及准稳态衰老（quasi-steady-state aging），即损伤追踪由产生与清除平衡所定义的移动设定点，以人类、犬、豚鼠和猫为特征。这些结果为比较衰老的 mechanics-based 模型基础提供了依据，有待实验验证。

不同生物以大致相似的方式衰老，伴随功能丧失及随年龄增长而急剧上升的死亡概率，但其寿命中位数差异巨大。解释生命树上衰老动力学的相似性与差异性，是长寿研究领域的主要挑战之一。寿命（或延迟死亡的机制）似乎是为使生物体繁殖并在死亡前以某种方式造福其亲属而进化出来的，死亡原因包括捕食、感染、意外或饥饿等外部因素。对长寿的分子机制解释源于短寿命与长寿命动物的比较，长寿动物具有增强的DNA修复能力、降低的分子错误率及其他相关过程。

研究人员采用名为饱和清除（SR）模型的数学模型，将损伤累积与生存曲线等人口学模式联系起来。该模型基于损伤可由标量x概括的简化假设，即超过一定x水平则生命无法维持。SR模型通过包含产生、清除和噪声的随机微分方程描述损伤动力学：损伤产生随年龄线性上升（ηt），而清除在高损伤时饱和（βx/(κ+x)），噪声以幅度ε的高斯白噪声建模，死亡发生在损伤超过阈值Xc时。该模型能解释衰老的多种定量模式，包括Gompertz和Weibull风险曲线、人类衰弱指数分布、疾病发病率曲线及寿命遗传性等。

为开展跨物种系统比较，研究人员组装了多个模式生物及人类的高质量生存曲线数据集，采用贝叶斯方法拟合各物种的损伤产生（η）、清除（β）、阈值（Xc）和噪声（ε）参数。数据来源包括：大肠杆菌（n=4,744，饥饿微流控实验）、酵母（酿酒酵母BY4742，n=24,512，基因组缺失筛选的野生型对照）、秀丽隐杆线虫（QZ0株，n=2,908，大规模衰老机器实验对照）、黑腹果蝇（4个不同寿命品系，各约1,000只）、遗传异质性HET3小鼠（国家衰老研究所干预测试项目对照，n~3,000）、英国伴侣犬（4个品种，每品种>1,000只）、英国兽医数据的猫和豚鼠、以及丹麦1890和1900年队列及瑞典1910年队列的人类数据（人类死亡率数据库）。这些生物的寿命中位数跨越四个数量级，从饥饿大肠杆菌的约100小时到人类的约80年。

**主要关键技术方法：**（1）基于SR模型随机微分方程的个体模拟（20,000次/参数集），作为It?过程实现，采用布朗桥校正减少步长敏感性，约束损伤非负；（2）贝叶斯参数估计：通过模拟死亡时间分布与实测数据比较，使用马尔可夫链蒙特卡罗（MCMC）方法（Python Emcee库）扫描参数空间，采用高斯核密度估计生成平滑死亡时间分布和生存函数，基于贝叶斯定律计算后验概率；（3）Shapley可加解释（SHAP）分析：以模拟中位寿命的对数为响应变量，评估各参数对寿命预测的贡献度；（4）风险函数分类：基于贝叶斯信息准则（BIC）区分Weibull型（幂律）和Gompertz型（指数）风险曲线。

**研究结果：**

**贝叶斯方法估计SR模型参数：** SR模型将损伤x的动力学描述为产生减去清除再加噪声的过程。方程为dx/dt = ηt ? βx/(κ+x) + √(2ε)ξ，其中线性上升的产生项ηt随年龄增长，饱和清除项βx/(κ+x)在高损伤时饱和，噪声项为幅度ε的高斯白噪声。死亡为x首次超过阈值Xc的首达过程。此外还引入外部死亡率m_ex（Makeham项）描述外部死因。通过20,000次个体模拟生成死亡时间分布，与实测数据比较后使用MCMC调整参数，获得各参数的后验分布。模型普遍能较好描述数据，仅线虫老年期死亡率下降未被捕捉。

**损伤产生参数是跨物种寿命差异的最佳单一预测因子：** 产生参数η在物种间变化达七个数量级，寿命越长则η越低。通过SHAP分析量化各参数对寿命的贡献，η的平均SHAP重要性为2.85，与次重要的阈值X_c（1.53）和清除β（1.38）之和相当，噪声ε最低（0.33）。方差分析（ANOVA）和单参数替换分析均确认η为最强预测因子。其他参数随寿命变化趋势较温和。

**物种呈现两种不同的衰老模式——弹道式与准稳态：** 根据产生与清除的相对关系，识别出两种模式。第一，弹道式衰老：酵母、线虫及多数果蝇品系中，损伤产生率ηt在整个寿命期内远大于清除，平均损伤近似二次方上升（x ~ ηt²/2），导致η与寿命呈反平方关系（η ~ 1/L²）。此类生物的风险曲线更偏Weibull型——风险先急剧上升后在中位寿命前转为幂律增长。第二，准稳态衰老：较长寿命哺乳动物及一种果蝇品系中，清除与产生在大部分寿命期内平衡，损伤追踪缓慢上升的移动设定点，仅在极老年时生产才超过清除，寿命粗略由L ~ β/η决定，η与寿命呈反线性关系（η ~ 1 adultery ~ 1/L）。其风险曲线呈Gompertz指数上升，在远超中位寿命的极老年才减缓。各生物处于两种模式的连续谱上， female小鼠约前2年处于准稳态，之后转为弹道式。

**哺乳动物具有近乎不变的噪声和清除参数组合：** 发现某些参数组合在哺乳动物中高度保守。噪声首达死亡阈值的时间T_ε = X_c²/ε约2年，在鼠、犬、豚鼠、猫和人中相差不超过3倍；清除与噪声时间尺度比βX_c/ε ~ 10，相差约2倍。这种守恒性提示噪声可能源于清除过程的系统波动，如免疫系统清除损伤/衰老细胞的昼夜节律波动。

**酵母复制性衰老具有可忽略的清除：** 酵母风险曲线在前约10代（约半数中位寿命）急剧上升后转为多项式增长。最佳拟合清除参数β极小，可忽略而不显著损害拟合质量，是弹道式衰老的极端情况。分析近似解显示晚期风险渐近于h(t) = (3/16)(η/ε)²t²，与观测曲线合理吻合。这提示酵母母细胞中存在未被有效清除的损伤累积，如滞留母细胞的受损线粒体、表观遗传状态或碱性液泡等。

**部分品系和物种的衰老模式与人类相似而其他则不同：** 以X_c为损伤单位、L为时间单位，三个无量纲参数（βL/X_c、εL/X_c²、ηL²/X_c）呈现对数线性关系。以噪声与产生时间尺度比s = η^0.5X_c^1.5/ε（反映过程确定性）和弹道-准稳态参数P₁ = (β/η)/L（反映清除与产生可比较的寿命比例）为坐标，最接近人类的物种为犬（除德国牧羊犬外）、猫、豚鼠、饥饿大肠杆菌和果蝇853品系，具有Gompertz型风险曲线；小鼠、线虫、其他果蝇品系和酵母则更接近弹道模式。

**讨论部分总结：** 研究人员通过损伤累积模型比较了不同生物的衰老，从高质量生存曲线中确定参数。变化最大的参数是损伤产生率η，其与寿命成反比；其他参数在物种间更相似。哺乳动物中噪声与清除时间尺度比守恒，提示噪声可能代表清除过程的波动。识别的两种衰老模式为：短寿命物种中生产占主导的弹道式衰老，以及哺乳动物中损伤追踪产生-清除平衡移动设定点的准稳态衰老。

长期以来，延长寿命究竟源于减少损伤产生、增强清除、降低噪声还是提高损伤耐受（死亡阈值）尚不明确。研究发现，解释物种间寿命差异的最佳单一参数是损伤产生率——酵母比人类高七个数量级。该模型可应用于不同生物中不同形式的损伤x；在哺乳动物中，x的候选者为衰老细胞。研究人员建议区分修复（repair，细胞内过程，影响产生参数η）与清除（removal，如免疫系统清除衰老细胞）。长寿哺乳动物的增强修复和降低错误率已在体细胞突变、表观遗传失调、转录和翻译错误等方面得到证实。

两种衰老模式的发现有助于解决风险曲线形状的争议：弹道式衰老呈Weibull型，准稳态衰老呈Gompertz型。并非所有短寿命生物均为弹道式——饥饿大肠杆菌和部分果蝇品系具有准稳态参数和类Gompertz风险曲线。女性小鼠处于模式过渡区。所有生物的死亡阈值X_c至少为清除饱和点κ的10倍，这种饱和在年轻阶段即发生，选择压力可能使死亡发生在远超清除饱和之后。

研究提出关于哺乳动物衰老过程中生物噪声起源的假说：守恒的T_ε ~ 10T_β关系将噪声与清除联系起来，可能源于免疫系统昼夜节律的随机波动，这与昼夜节律在衰老中作用的研究领域相关联。

基于比较分析的一个结论是，进化可能主要通过改变损伤产生率（通过增强分子修复）来调节物种间寿命数量级的差异，而非清除率或死亡阈值。其他参数的相关变化提示长寿是产生、清除、阈值和噪声的协调改变。因此，旨在大幅延长寿命的未来干预需聚焦于产生参数，而非仅关注清除或阈值。近期分析表明，仅降低产生率会同等比例延长寿命和疾病期（sickspan）；因此要延长寿命并压缩疾病期，需同时降低产生、提高阈值、减少噪声和/或增强清除的联合干预。

尽管存在差异，各生物在三维无量纲参数空间中呈对数线性关系，提示从人类（大参数极端）到酵母（小参数极端）的策略连续谱。与人类参数相近的物种可能最适合作为人类衰老动态研究的模型，包括特定果蝇品系、猫、多数犬种和豚鼠。

本研究的局限性在于物种寿命和SR参数均从相同生存曲线估计，因此推断的参数是数学构建而非直接测量的生物学量，这些基于数学的假说有待未来实验验证。

**研究结论翻译：** 研究人员比较了物种生存曲线以确定SR损伤累积模型的参数。发现两种衰老模式：生产占主导的弹道式衰老，以及生产与清除平衡直至极高龄的准稳态衰老。寿命的最佳单一预测因子是跨越七个数量级的产生率；其他参数也协调变化。这揭示了衰老中随机噪声的起源，表明某些衰老参数在类群间的近普适性。进化可能主要通过改变损伤产生率（通过增强分子修复）而非清除率或死亡阈值来调节寿命，暗示旨在延长最大寿命的干预的潜在靶点。