突变积累得慢，是否就能活得久？

复核后的答案是：肠道隐窝中的负相关相当稳健，不能据此认定突变率决定寿命。细胞分裂、修复、组织更新与生活史可能共同演化；干预一条路径还可能同时改变其他病理过程。DREAM研究就是一个具体提醒：脑组织突变较少的遗传模型，其来源研究却报告了系统性淀粉样变性和更短生存。

2026-09-20更正： 原版系统发育协方差写反，并用过大的对角修补和不完整似然计算，造成“λ=1时关联消失”的错误。现已取得作者原树，用Python和R分别重算；撤回深部支系主导及与作者结论冲突的推断。同时更正汇总权重、缺失数据、偏相关自由度与物种名称，补入细胞分裂及DREAM原模型的反证。

01先弄清楚数据中的一个点是什么

Cagan等研究的是结直肠隐窝。一个隐窝可由一支干细胞谱系占据，因而能从其细胞群中识别累积的体细胞突变。论文涉及16种哺乳动物，寿命关联分析使用15种；港湾鼠海豚缺少年龄信息，不能计算年率。原稿将其误译成了海豹。s1

公开逐隐窝表有208行、53个不同的“Individual”编号。两行鼠海豚缺年龄；另两行环尾狐猴有年龄但替换年率记为缺失。因此现成可分析的年率有204个，来自52个表内编号、15种物种，不能说四行都是年龄缺失。跨物种关联的信息单位主要是物种，并不是204次独立的寿命比较。s2s4

原文代码先在个体内平均，再在物种内平均；原稿直接平均隐窝，给隐窝较多的个体更大权重。修正后14种物种的均值与作者表几乎精确一致。人类的差异也找到了算术来源：按表内6个Individual编号平均为46.395/年；把PD36813再按ac、ad、x三组样本前缀分别计权，就精确得到作者的47.125/年。这三组也对应不同的匹配参照样本。s10s11s14

但算术对齐不能证明它们是三名独立的人。物种汇总表写8名人类、2只狐猴，逐样本表却分别只有6个和1个Individual编号；这些“个体”的生物学对应仍未完全解释。下面以作者发表的物种均值核对原模型，同时把按表内编号等权重算作为单独的敏感性分析，保留这个界限。s12

02负相关在正确的系统发育模型里保留下来

先用作者15个物种均值，对log₁₀年突变率与log₁₀寿命作回归：斜率−0.862，95%区间−1.078至−0.645，R²=0.851。这与“寿命越长，年率越低”一致；斜率区间包含−1，但并不证明恰好遵循严格的反比例定律。

近缘物种共享演化历史，需要在模型中反映这种非独立性。原代码恰好做反了：用“共同祖先距今年数”直接作协方差，使人和鼠比狮和虎更相似。它生成了不合法的协方差矩阵，又加了远大于舍入误差的对角项，不能再解释为原假设下的系统发育模型。

本次从作者Zenodo存档取到了真正的TimeTree树，以从根到共同祖先的共享支长构造协方差，不作任意对角修补。结果如下：s13

λ描述残差中的系统发育相关强度，不是“关联由演化解释的百分比”。估计λ时还必须同时估计残差方差；原代码没有完成这一步。表内区间仍以当前树和λ为条件，没有合并树、λ及寿命测量的全部不确定性。

还有一个原稿忽略的比较口径：作者发表的PGLS是原尺度上“年率对寿命倒数、零截距”的模型，而不是上表的自由截距log–log模型。按原代码公式和原树，系数k由普通回归的2439变为PGLS的2000；去掉鼠和大鼠时分别为3326和3078。它们属于同量级的估计，但不能用“k显著大于零”证明反比例律；零截距和倒数形式本身就是模型设定。s10

Python的树解析和GLS计算，与独立R实现ape/caper的225项协方差和90项参数核对通过。此前的“校正后关联消失”来自我们的实现错误，不能当作反驳原研究的证据。

03哪些稳健性结果，仍然有条件

五种可分析物种只有一个表内Individual编号。我们在物种内重抽这些编号，得到斜率约−0.907至−0.824的95%自助区间；这个很窄的区间只反映已观察到的物种内波动。单编号物种重抽多少次都不会产生新的个体差异，它没有覆盖跨物种抽样、系统发育或未知生物个体对应的全部不确定性。

寿命本身也不是一个无歧义的数。主分析的Lifespan80是达到成年后、80%已经死亡的年龄；动物主要来自圈养记录，人类来自三个国家1900年出生队列，狗的记录还合并了犬与狼。它不是平均寿命，也不是所有物种都用同一种来源测量。s1s3

换用50%、90%死亡年龄或AnAge最大寿命，负相关仍在。但若直接替换Lifespan80置信界限或90%寿命，因缺失会少掉不同物种；“六种口径”并不都是同一15种。固定为六套数据都齐全的11种后，斜率仍介于−0.906和−0.816。这是口径敏感性，不是寿命误差已经被统计消除。

控制体重的偏相关约−0.907，按一个协变量正确扣除自由度后p≈7.8×10⁻⁶。它表明这份数据中的关联不能被一个log体重协变量消除，但没有排除其他共同演化因素，更不等于体重的因果作用已被分解。

04“一生的突变差不多”放在什么位置

作者物种年率乘Lifespan80，得到的预测终生负担为1828—5379个替换/基因组，约2.94倍的范围。

这个收敛现象值得解释，却不是另一组独立验证：当年率近似随寿命倒数缩放，两者的乘积自然更接近。供者通常也不是活到各物种Lifespan80才被测序。由此既可提出“损伤累积参与限制寿命”，也可提出“较长的生活史推动更强维护、更慢更新”，还可能是两者共同作用。它不能直接识别一个普遍的致死突变阈值。

05细胞分裂这个竞争解释，并没有被排除

2022年原文使用少数物种、不同方法的分裂时间估计，提示细胞分裂不足以独自解释全部差距。原稿进一步写成“解释不了主体”，超出了证据。

2023年Montazid等直接测量小鼠和裸鼹鼠的肠道干细胞动力学，并结合人类资料，发现隐窝底部干细胞分裂率也随寿命变化。其结肠估计为小鼠181次/年、人类50次/年：按这些数值归一化，鼠/人的年突变率比16.90降为每次分裂归一化比4.67。分裂值得作为重要竞争解释保留；该研究只有三个物种的跨物种对照，也不能独自确定普遍规律。s16

这仍不是“有多少百分比由分裂介导”的因果拆分。若年突变率r=d×μ+ν，其中d是分裂频率、μ是每次分裂相关贡献、ν是其他随时间累积的贡献，那么r/d=μ+ν/d。将突变数除以分裂次数，并不能自动把非分裂来源去掉。

组织边界同样具体。Kapadia2025测得小鼠造血干/祖细胞的年龄斜率约45.3个SBS/年（42.2—48.4），在人类造血细胞中相关研究约为14—17个。该研究还在同一批老鼠中比较血液与肠道隐窝，验证了组织差异。45.3的年龄斜率与Cagan的累计负担除年龄汇总并不完全同义，不能把二者比值直接称为普适的“18倍组织差”。隐窝结果也不能直接推广为全身统一的突变时钟。s6s7 该文勘误把图4e上方重复的5×10⁻⁵刻度改为5×10⁻⁴，不改变这里引用的年斜率。s8

06DREAM扰动：改变突变，不等于延寿

Koch2026把表达推算的DREAM活性与多个尺度联系起来：跨物种寿命、细胞突变及人类脑病理。表达—突变分析部分来自同一RNA数据；人类疾病分析还使用死后脑表达去关联较早发生的诊断，不能据此确立前瞻性保护或因果方向。s9

其干预部分更直接：8周龄诱导破坏DREAM装配，最后分析4只功能缺失鼠和5只对照的存档脑组织，报告协变量调整后的SBS减少4.2%、indel减少19.6%。起初两组各5只，剔除一只离群样本后得到上述人数。脑组织在死亡时取得，模型包含年龄、测序覆盖和读段家族等项；小样本、模型选择及终点取材都需要考虑，不能当作同龄纵向轨迹的实测减速。

更关键的是，这些组织来自Perampalam2021的模型。那项研究的功能缺失组30只、对照37只，报告约16%的生存缩短（log-rank p=0.0236），并出现心、肝、脾、肾的淀粉样沉积与器官损伤。其生存终点包括达到人道终点后安乐死；2026稿称脑在自然死亡后取得，但本次没有逐鼠对应表，不能假定两篇的取材定义完全相同。s15

这并不证明降低突变本身有害，也不意味着所有选择性、暂时性DREAM调节都会复制这种广泛遗传破坏。它说明：一项干预同时影响多个过程时，较低的某组织突变负担不足以证明净健康获益。真正打通证据链，需要在明确的干预、组织和时间窗口下，把突变、功能、安全和生存一起检验。

07可以保留的结论

这份数据中，肠道隐窝年突变率与物种寿命的负相关，经正确系统发育模型、汇总口径与寿命定义检查后仍然存在。解释空间包括修复、分裂、选择和生活史的共同变化，尚不能由相关性决定因果方向。终生负担收敛是同一关系的派生描述；DREAM为突变可被干预提供线索，同时原模型的疾病与生存结果阻止我们把它直接写成延寿证据。

本次没有从原始测序重新调用变异，也未完整重跑原论文的层级混合及贝叶斯模型；复核的是公开表、作者代码对应的物种模型及新增敏感性。编号与生物个体对应的缺口、跨研究估计对象差异、树与λ的不确定性均保留。证据核查至2026-09-20；Codex完成返修、自审及编辑，不构成人类专业审阅或干预建议。

复算包包含固定来源表、原始树、分析代码、独立R验证与输出；第三方公开数据按其CC BY 4.0许可注明来源，原论文图未再分发。