44岁和60岁，身体真的会突然变老吗？

01先说我们得到的答案

目前最有依据的表述是：分子随年龄的变化不是匀速的，这一点有多个来源的支持；但"44岁和60岁发生两次群体分子剧变"这一具体时刻表，其统计证据来自一个会把平滑变化、甚至纯噪声都渲染成双峰的分析管线。在正确构造的零假设下，同样的波照常出现。 这既不等于原论文"全是错误"，也不等于非线性衰老研究没有意义——被削弱的是那两个精确到岁的数字，以及围绕它们的确定性叙述。

本篇的五项计算各回答一个问题：

范围限定为这项研究的群体计数曲线及其统计识别。不评价一切"非线性衰老"文献，不提供按年龄择时的健康建议，也不把分子计数峰外推为疾病、功能或寿命结局。

02这个说法从何而来，现在处于什么状态

2024 年 8 月，Shen 等在《Nature Aging》发表一项纵向多组学研究：108 名 25–75 岁参与者、中位随访 1.7 年，每人多次采样后取平均，得到横断化的 11,305 个分子变量。论文摘要报告"约 44 岁和 60 岁两个主要时期发生显著分子失调"，这句话随后以"衰老断崖"的形式广泛传播，本站旧版也曾按表面结论报道。原论文

2026 年 6 月，Carbonneau 等发布方法学预印本，系统论证该管线的三个伪影来源；值得注意的是，原论文的第一作者与通讯作者（Shen、Snyder）是该预印本的共同作者。质疑预印本

2026 年 7 月 15 日，期刊在论文页面挂出编辑说明（原文照录）：

Readers are alerted that, following further review by the authors and independent researchers, the reliability of some of the analyses used in this article to detect peaks of changes in molecular aging markers around certain ages is currently in question. … Appropriate editorial action will be taken in due course.

大意：经作者与独立研究者进一步审查，本文用于检测特定年龄分子衰老标志物峰值的若干分析，其可靠性目前存疑；编辑行动将在适当时候采取。截至本篇检索日（2026-09-17），编辑行动尚未落地，论文既未更正也未撤稿。

所以本篇要回答的不是"这项研究红不红"，而是一个更窄的问题：那条给出 44/60 的计数曲线，在合理的分析选择下能否稳定地区分真实信号与管线行为。

03一条计数曲线是怎么做出来的

作者的管线分三步（代码已固定提交存档）。作者代码库

第一步，横断化：把每人多次采样取平均，每人只剩一个向量。我们恢复的转录组横断矩阵为 8,556 个基因 × 97 人（108 名参与者中进入转录组层的人数），年龄 25.9–75.2 岁、中位 56.2 岁，其中 25–40 岁只有 7 人。

第二步，平滑插值（LOESS）：对每个变量沿年龄做局部回归，在 26–75 岁间每 0.5 岁插出一个值。从此数据不再是 97 个人，而是 91 个等距网格点。

第三步，滑窗计数（DE-SWAN）：对每个中心年龄 c，比较窗口 [c−10, c) 与 [c, c+10) 内的"样本"，逐变量做 Wilcoxon 检验并做 BH 校正，数出显著变量个数，连成曲线。

两个设计后果值得记住。其一，平滑让相邻网格点高度相关，窗口比较里的"样本量"有相当部分是插值造出来的；其二，真实参与者在各年龄分布不均，平滑与窗口检验的功效随年龄变化。这两点后面都会回来。

04先复现：曲线一模一样，但基线高得惊人

我们先验证自己的实现：在作者的合并组学网格矩阵上重写滑窗检验，26 个中心 × 11,305 个变量共 293,930 对 p 值与作者保存的输出零差异，计数曲线逐值一致；把窗宽从 20 年换成 15、25、30 年，依然逐值一致。实现是忠实的。复算包（含逐值对照说明）

复现之后看曲线的真实形状：每个中心都有 60%–82% 的变量显著。峰在 44（9,312 个）与 60（8,276 个），谷在 54（7,117 个）——所谓"两次剧变"，是在"几乎所有变量在所有窗口都显著"的基线上，显著比例从约六成涨到八成再回落，峰谷比约 1.3。论文自己的另一项分析也称 81.03% 的变量至少在某个年龄段显著变化。原论文

这不是挑刺：计数 = 真实效应 × 检验功效 × 变量间相关结构的混合体。把这条曲线直接读成"分子变化速率"，本身就跨过了一步未加论证的翻译。

05峰依赖哪一步：把平滑去掉之后

我们接着在恢复的 97 人转录组上完整重走管线：逐基因 LOESS（span 按作者的交叉验证协议选取）→ 网格 → 滑窗，得到的计数曲线与作者的转录组曲线几乎重合：26 个中心最大相差 74 个（约 6,400 的 1.1%），相关 0.9993，峰、谷、第二峰的位置一致。作者曲线复算

然后只做一处改动：去掉 LOESS，把 97 人的原始值直接送进同一个滑窗检验。结果：26 个中心、8,556 个基因，FDR 0.05 下显著基因数全部为零。质疑方在论文数据上做了同一实验，同样全零。质疑预印本 逐中心结果

含义要读准：计数双峰依赖平滑插值这一步骤，它不是原始数据经过同一检验直接呈现的属性。反过来，全零也不证明"这些基因与年龄无关"——97 人、窗内人数有限，功效本来就低。它证明的是：这条曲线的形状归管线所有，不归数据所有。

06为什么论文自己的置换检验"通过"了

论文用置换检验论证峰不是巧合。但看其实现顺序：先完成 LOESS 插值，再打乱 91 个网格点的年龄标签，重新计数。打乱发生在平滑之后，平滑制造的相关结构被一并打碎，于是零假设下几乎长不出波，观测曲线自然显得"显著"。我们按作者方案复跑 200 次：168 次所有中心的计数不超过 100，最大单中心计数 2,979——而观测值在 5,000–7,000，对照之下显得"绝不可能是巧合"。

统计上正确的做法是先打乱 97 名参与者的年龄，再走完整管线（质疑方称为 Algorithm 1）。这样零假设代表的是"没有真实年龄结构、但分析行为原样保留"的世界。结果：波照常出现。质疑方跑了 80 次，各中心计数 5,955–7,738；我们用恢复的矩阵独立跑了 100 次，各中心零假设中位数 5,317–6,995，95% 带上限约 6,800–7,900。观测曲线 26 个中心全部落在零假设带内；就连"波形有多明显"（峰谷摆幅）也是观测 3,317、零假设中位数 3,952——63% 的零假设运行，摆幅比真实数据还大。质疑方置换结果 本文置换结果

一句话：作者方案的对照组是"没有平滑结构的世界"，而不是"没有年龄效应的世界"。用错了对照，检验自然通过。

07纯噪声、线性真值，都能长出同样的双峰

如果峰是管线的产物，那么连信号都不需要提供。我们生成 2,000 个纯噪声"基因"，按 97 名参与者的真实年龄走同一管线，重复 100 次：每一次都长出双峰——峰 1 全部落在 43–46 岁，峰 2 全部落在 60–65 岁；其中 66 次的双峰同时落在 43–46 与 60–62 岁，与观测位置几乎重合。把年龄换成同范围均匀分布再跑 100 次：峰位散开（峰 1 多在 40–42 岁，峰 2 散布 57–63 岁），波形摆幅中位数从 443 降到 150。年龄分布的不均——25–40 岁仅 7 人——和平滑在数据边缘的行为，共同决定了峰落在哪。模拟结果 均匀年龄对照

再进一步：让每个基因的真值严格线性（只有斜率和噪声，100 次模拟）。这时同一滑窗检验（不加平滑）的计数曲线依然大幅起伏——各中心显著基因中位数从 310 到 1,177，约 3.8 倍。真值明明匀速，计数却不平：因为窗口里的人数随年龄变化，功效随之起伏，计数跟踪的是功效，不是效应。把平滑加回去，计数曲线变平，但整体抬升到每中心约一半基因显著——插值伪样本让"显著"变得廉价，这正是观测曲线那 60%–82% 高基线的来历。线性真值模拟

08不数个数，直接看一眼效应量

计数可以被功效牵着走，效应量不会。我们换一个做法：不数"显著"的个数，直接计算每个窗口两侧人群表达差异的效应量（Cliff's δ，可理解为"右边的人比左边大"的概率与 50% 的距离），再看全部 8,556 个基因的中位 |δ|。

结果有两层。第一，效应量普遍很小：所有窗口的中位 |δ| 在 0.08–0.17 之间，换算过来，组间方向差异离纯随机最多约 4–9 个百分点——与"六成到八成分子剧变"的计数叙事完全不在一个量级。第二，变化确实不均匀：把纯噪声基因按同一年龄分布跑 200 次画出 95% 噪声地板带，观测曲线在 45–50 岁与 56–59 岁窗口小幅超出带上沿；在 40–43、51–55 与 60–65 岁则落在带内甚至略低于带——低于带与测量值中的大量并列有关（74% 的基因在 97 个样本中存在重复取值，会把 δ 压向 0）。效应量曲线 噪声基线

这是本文一个与"全是伪影"不同的发现：窗口效应在四十多岁末段与五十多岁末段有小幅真实富集——非匀速是存在的，只是量级是几个百分点的概率偏移，而不是断崖。 同时必须标注：该矩阵未见混杂校正痕迹，队列构成（BMI、性别、胰岛素抵抗、族裔）随龄变化也可能贡献部分窗口差异；这小幅富集需要在独立队列确认后才能谈归属。

09"44岁"和"60岁"能精确到岁吗

即便暂且接受这条管线，两个数字的精确性也经不起三方面的推敲。

其一，重抽样。对 97 名参与者做有放回重抽样、重走管线 100 次：峰 1 的位置在 40–50 岁之间漂移（95% 范围，中位 44 岁），峰 2 在 57–65 岁之间（中位 61 岁）；波的摆幅同样不稳，95% 范围从约 1,200 到 9,000 个基因。自助结果

其二，论文内部就在漂移：摘要写 44/60，结果一节写 45/65，图 4a 图注写 44/60，蛋白质组一节写 40/60，讨论定稿写 40/60，轨迹聚类的转折在"约 60"。同一篇论文给出过至少五组数字。原论文

其三，跨研究并不重合：Lehallier 等 2019 年在 4,263 人的血浆蛋白组上用同族方法报告的三个峰在 34、60、78 岁——第一个峰与 44 差了十岁。Lehallier 2019

所以合理的读法是"中年段"与"60 岁前后"两个宽区域，而不是两个精确到岁的时刻。

10那什么还站得住？

把一张时刻表放回抽屉，不等于把整片领域否定掉。按现有证据，以下三点仍然有依据：

分子变化不是匀速的。 我们自己的效应量重算就看到 45–50 与 56–59 岁窗口的小幅真实富集（超出纯噪声地板，但只是几个百分点的概率偏移）；论文的轨迹聚类分析（与计数曲线不同的另一条证据线）显示多组学轨迹存在非线性形态；Lehallier 的三峰用的是没有 LOESS 前置的原始 DE-SWAN（线性模型加 F 检验），其受试者级置换在他们的管线里是有效的，且置换后峰消失、窗宽稳健。Lehallier 2019 Pyrkov 等用完全另一类方法（纵向血常规的恢复力动态）发现，以自相关时间度量的恢复时间从 40 岁到 90 岁单调上升约四倍（约 2 周到 8 周以上）——"恢复力丧失"是平滑过程，不需要断崖。Pyrkov 2021

60 岁前后是一个反复出现的候选区域。 Lehallier 的 60、论文聚类的"约 60"、以色列 HPP 队列 10,019 人预印本的 60–65 宽峰，都指向同一地带。注意 HPP 预印本沿用滑窗家族方法且未经同行评议，同族方法的相互复现不能排除共同的方法学伪影，所以它只是"候选"，不是"证实"。HPP 预印本

部分波动可能有生命阶段的解释。 HPP 预印本中，加速衰老子群的女性在 50–55 岁出现尖峰，与绝经平均年龄重合，男性峰延后 10–15 年；Pyrkov 也提到女性中段拐点或与绝经相关。激素转变是真实存在的生理事件，它可以在群体曲线上留下痕迹，而不需要"所有人同岁断崖"。HPP 预印本

11四种解释，现在怎么判

需要诚实说明：平滑累积、个体年龄不一与管线伪影三者可以并存，这份横断数据不足以在它们之间完全裁决。被推翻的是"精确时刻表"这个版本，不是"衰老非线性"这个方向。

12我们知道什么，还不知道什么

知道：作者的关键输出可以完全复现；计数双峰的统计显著性来自一个无效的置换对照；该信号无法与管线自身行为区分；"44/60"应改读为两个宽区域；编辑审查正在进行。

不知道：个体层面是否存在转变、发生在什么年龄（需要数十年纵向追踪）；转录组以外的组学层在独立重算下表现如何（本次只独立复算了占合并变量 76% 的转录组）；期刊将采取何种编辑行动；其他人群中同法表现。

13研究说明、复算与更新

这是一次针对单一论文主结论的复算与压力测试，不是系统综述，也没有做新的湿实验。检索截止 2026-09-17，覆盖原论文及其编辑说明、方法学质疑预印本、双方代码库与同领域对照研究，检索式与覆盖缺口见研究记录。

下载复算包：包含全部复算程序、固定输入与哈希、分析计划（含先登记后执行的判定逻辑与工程偏差记录）及输出说明；作者代码库按 MIT 许可使用，参与者级数据来自作者仓库的 git 历史对象，不随包分发论文全文。分析计划（ANALYSIS_PLAN.md）与全部输出说明（README.md）随包附载。

审计中发现的另外三件事实，如实记录、不下结论：论文内部峰位数字至少五组且互相不一致；作者仓库中转录组峰 1 与峰 2 的分子清单文件逐字节相同（明显的复制错误）；方法段声明"所有分析前校正 BMI、性别、胰岛素抵抗与族裔"，而转录组滑窗的代码路径中未见该校正步骤，矩阵数值亦无残差化痕迹——作者声明与代码可见性之间的差异，已列入待澄清事项。作者代码库

计算分开标注：A 系列是对作者保存输出的逐值复算；B 系列是从恢复的 97 人矩阵出发的管线重跑、置换、自助与受控模拟；C 系列为外部证据的论文层面核对。固定随机种子、先登记的判定逻辑、与作者输出的逐值对照用于核查实现，不构成新的实验验证。

本篇更新长寿派原《自然·衰老》断崖报道：撤回既定化表述与行动建议，保留观察性谨慎提示，补入复算、置换检验与 2026 年编辑说明后的证据状态。首次公开此版本：2026-09-17；下次计划检索复核：2027-03-17；期刊编辑行动、质疑预印本同行评议版本或作者正式回应时提前更新。

作者、计算、译文核对及编辑均由同一 Kimi AI agent 完成，没有独立 agent、人类或临床专业审阅。逐条来源的实际读取范围、研究利益关系与剩余限制见下方记录。