部分重编程，让什么变年轻了：读数、细胞身份，还是功能？

结论是：分子状态确实发生了变化，撤回因子表达后的细胞在所选标志物上重新接近成纤维细胞；但变化幅度依赖供者、处理时间、算法和标志物选择。现有分析不能把变化分解为“真实年龄恢复”和“去分化”的比例，也不能把培养细胞读数换算成人体年轻了多少年。 功能实验提供了有限而有价值的线索，尚未建立时钟变化与功能改善之间的可靠对应关系。

2026-09-19更正： 本次复核修复了缺失CpG处理、PhenoAge截距、RNA跨矩阵标准化及重复测量单位；撤回原稿“约9年反驳30年”“显著相关证明去分化主导”“共同OSKM疤痕”等判断，并纠正三条错误PMID和一项作者归属。下文数字、六张图及复算包均已更新；原作者的自审保留为历史，本次修正由Codex负责。

01实验测了什么，比较单位是什么

MPTR部分来自三名供者，年龄为38、53、53岁。研究比较未感染载体、经历相应培养和分选程序的阴性对照，与表达因子后按表面标志分出的重编程群体及未进入相同状态的群体。后者简称“失败组”，是分选后的细胞状态分类，并非随机分配的“失败处理”。设计与分选

公开数据中，瞬时重编程部分有96个甲基化样本、95个RNA样本；另有124个Sendai重编程RNA参照样本用于定义共同评分坐标。样本列是培养物或分选群体，不能被当作同等数量的人。我们从GEO的SOFT元数据核对供者、实验批次和天数；例如“iPSC 13”“iPSC 21”是样本名称，元数据中的重编程时间为51天，不能把名称里的数字当成暴露天数。GEO

主比较先在同一供者、同一实验、相同暴露时长和采集阶段内，减去匹配阴性对照，再在同一供者与时长内平均重复实验。各时长分别给供者等权，不把重复培养扩大成独立生物学样本。第一个实验撤药后四周采集，第二个五周采集；这是两个实验的终点，不是同一批细胞从四周到五周的连续追踪。

数据并不平衡：撤药后重编程组在10、13、15、17天处理条件下分别覆盖3、2、1、3名供者。因此，跨处理时长的差异不能直接被叫作最佳剂量窗口。我们不从三名供者报告总体显著性，也不把最有利的一个条件挑作普遍效果。

02时钟读数：有变化，但不是一把统一的尺子

五个时钟都使用作者提供的处理后NOOB甲基化值。修正计算保留每个模型的截距与变换：PhenoAge包含513个CpG和一个60.664的截距，“截距”不是第514个缺失位点。Horvath 2013、Hannum、EpiTOC分别有19、6、31个位点不在当前矩阵中；这次用固定版本Biolearn的外部参考/均值混合插补规则，不能再把全缺失行的空均值称为已完成插补。算法与系数

下表为撤药后重编程组减匹配对照的供者平均差；负值表示模型读数较低，括号是供者差值的最小值到最大值，不是置信区间。

这支持某些培养条件下时钟读数降低，也显示供者间差异和模型敏感性。Skin & blood的391个位点在数据中全部存在，所以它不依赖补齐整个位点；但输入完整不等于已经验证该时钟能测量重编程后的功能年龄。

血液训练时钟也不能概括为“完全不动”。例如PhenoAge在四个时长的配对平均差分别为−17.0、+12.3、+7.3、−9.5年，方向不一致。EpiTOC是另一种分数，不能把它的小数变化与时钟年数直接相比。这些差异说明模型响应不同，不足以挑出一把“真实年龄尺”。

技术核验中，五种时钟×96个样本×两种探针处理共960项结果，与固定版本上游公式一致到浮点精度。额外将检测P值大于0.01的探针值设为缺失再插补，单样本Horvath 2013最大变化约1.27年，Skin & blood约0.015年。这个敏感性分析只检查所用位点，不等于重新完成原始芯片质控，更不验证时钟的生物学解释。

为什么撤回“约9年而非30年”？ 原论文的约30年涉及其转录组估计和甲基化分析；本库的旧数字混用了不同供者、时长和未正确实现的算法。即使数字方向类似，不同预处理、模型、平均方式与比较对象也不能直接构成反驳。本次没有重跑论文完整转录年龄模型，也未证明其时钟实现与本库完全相同。Gill的方法写的是minfi/NOOB，不能额外断言作者做了所谓“黄金标准归一化”。准确的增量是上表透明的匹配比较及其不确定边界。论文方法与图3、4

03身份评分：用同一个坐标系，结论才可比较

我们预选13个成纤维细胞标志和12个多能性标志，将每个基因按Sendai参照矩阵中的均值与标准差变换，再把同一组固定变换应用到MPTR矩阵。两类标志的差值以参照起始成纤维细胞均值为0、参照iPSC端点均值为1。这是探索性表达评分；0.3不代表30%的细胞已重编程，也不是细胞纯度。

修正后，MPTR重编程群体的轴中位数从暴露期约0.26回到撤药后约−0.23，MPTR的iPSC参照约1.02。阴性对照也落在负值区域：两个系列之间仍有培养和批次差异，不能强迫所有成纤维细胞都等于0。方向上，这与论文观察到的成纤维细胞样特征重获相容；不能据此证明每个细胞完全回到原身份。原文身份分析

跨模态只连接相同培养条件标签的群体读数，不是逐细胞轨迹。分别看三名供者，暴露期Skin & blood时钟与身份轴的描述性相关系数约为−0.90、−0.83、−0.55；撤药后约为0.06、0.04、−0.18。这里的点仍有条件与时间结构，相关性也混合了对照和不同分选状态。它们提示两个读数在暴露期共同变化，不能确定哪个过程主导，也不能计算去分化解释了多少年龄变化。原先把68个组格视为独立样本得到的P值已撤回。

04“多能性残留”对标志物选择有多敏感

原稿把撤药后多能性分数抬高称为“共同OSKM疤痕”。但12个标志中含POU5F1和SOX2，它们本身属于被外源表达的因子；当前基因级数据不能可靠区分内源和载体来源。

同供者配对后，重编程组四个时长的12标志平均分差约为+0.17至+0.42；去掉上述两项后，10标志的差降为约+0.001至+0.096。失败组也对这个选择敏感。单位是基于参照标准化后的基因分数平均值，不是已经标定的生物效应SD。

这不是证明差异全由载体造成，也不是证明剩余小差异无意义；它表明“持续多能性程序”的结论依赖评分定义。由于成功、失败群体由处理后状态分选，失败组同样不能给“选择效应”画出因果上限。

四种解释可以同时发生，当前设计不足以替它们分配比例。

05功能实验存在，但不能由时钟代替

Gill论文报告了胶原I/IV表达与免疫荧光结果，以及划痕迁移实验。胶原I转录恢复未达其统计阈值，蛋白读数有向年轻方向改变；迁移比较中部分配对培养群体改善，部分没有明显改善。作者也指出迁移变化没有表现出与其他年龄指标的一致对应。图3及补充图、方法和审稿回复

这些是功能与表型线索，不能因GEO组学矩阵不包含它们就写成“没有功能证据”。反过来，小样本未显示相关也不证明时钟与功能已经脱钩。本次没有取得可与组学条件可靠连接的功能原始数据，故没有重新估计时钟—迁移相关，更不能把伤口培养皿结果当成人体组织修复或健康寿命。

其他方案也报告过培养细胞的年龄相关表型改善，例如Sarkar等的非整合瞬时表达研究；不同细胞类型、因子组合和测量流程之间，没有直接共同对照就不能排名“年轻化多少年”。Sarkar 2020

06安全与体内获益：必须保留具体条件

撤药后四至五周采集到某种状态，不证明它在此前每一天连续稳定，也不能回答长期肿瘤或器官风险。以下小鼠证据提供边界，不是对Gill培养细胞结果的直接重复。

组织损伤和细胞衰老本身还会影响重编程环境，进一步限制简单的全身外推。Mosteiro 2016 这些研究的物种、因子、递送和观察期不同；本篇没有系统合并其效应，也不提供重编程操作或用药方案。

07本次分析的范围与复算

本次使用作者处理后的甲基化与表达矩阵，未重跑原始IDAT或FASTQ，也未复现整套转录年龄模型。三个MPTR供者、分选后的bulk样本、不同条件的缺样、参照批次和探索性标志评分，均限制推广。即使程序计算正确，也不能自动把模型输出解释为生物学年龄的真实改变。

复算包包含修正代码、输入来源及哈希、配对输出、标志物参考和图表；原始大矩阵按清单从GEO取得。旧程序与结论留有内部修订快照，当前附件不再保留失效的显著性表。来源核查更新至2026-09-19；正文及英文版由本次返修agent自审，这不是独立临床或人类专业审核。