编辑推荐：
背景/目标: 在众多提出的衰老理论中，越来越多的研究指出表观遗传改变，特别是DNA甲基化模式的变化，是生物衰退的关键因素。DNA甲基化可以被靶向和调控，以诱导抗衰老和再生效应，这种新兴的治疗方法被称为部分重编程，由Yamanaka因子的不同组合诱导。方法: 在
背景/目标: 在众多提出的衰老理论中，越来越多的研究指出表观遗传改变，特别是DNA甲基化模式的变化，是生物衰退的关键因素。DNA甲基化可以被靶向和调控，以诱导抗衰老和再生效应，这种新兴的治疗方法被称为部分重编程，由Yamanaka因子的不同组合诱导。方法: 在本研究中，研究人员对来自多个哺乳动物物种、组织、细胞类型和转录因子组合的部分重编程实验的五个公开可用的DNA甲基化数据集（n = 189个样本）进行了整合二次分析。结果: 在对929,009个CpG位点进行差异甲基化分析后，研究人员识别出14个显示出可重复且显著变化的CpG位点，其中两个CpG在所有数据集中一致存在。结论: 这些发现表明，部分重编程可能通过选择性调控保守的DNA甲基化位点，部分发挥其抗衰老和再生效应。
**论文解读：保守的部分重编程对甲基化组的影响**
**研究背景与问题**
衰老是一个复杂的生物学过程，其特征包括功能衰退、组织稳态丧失和分子损伤积累。在众多衰老标志物中，表观遗传修饰，尤其是DNA甲基化模式的变化，已成为生物年龄的可靠指标。近年来，一种称为“部分重编程”的新兴治疗策略引起了广泛关注。该策略通过表达不同组合的Yamanaka因子（OCT4、SOX2、KLF4、c-MYC、LIN-28和NANOG，统称为OSKMLN），在不诱导完全去分化的情况下，促进年龄相关表型的再生，并逆转表观遗传时钟，甚至使人类真皮成纤维细胞（HDFs）的表观遗传年龄降低多达30年。然而，不同研究之间在物种、组织、细胞类型、因子组合（OSK、OSKM、OSKMLN）以及处理方案（如短期表达、长期循环诱导）上存在显著异质性，导致部分重编程的效果在不同模型中表现不一。这种异质性提出了一个核心问题：是否存在一些特定的甲基化组区域，能够不受物种、组织类型或处理方案的影响，对部分重编程产生一致的反应？识别这类保守的甲基化标记将有助于理解部分重编程是通过全局表观遗传调控还是选择性调控关键调控位点来实现其抗衰老效应。截至目前，尚未有研究系统性地回答这一问题。因此，研究人员开展了这项跨数据集的整合分析，旨在鉴定部分重编程过程中保守且可重复的甲基化特征。
**研究内容与结论**
为探究保守的甲基化变化，研究人员对来自三个哺乳动物物种（人类、小鼠、大鼠）、多种组织与细胞类型、以及不同部分重编程方案（OSK、OSKM、OSKMLN）的五个独立DNA甲基化数据集（n=189个样本，对照91个，处理98个）进行了整合二次分析。通过严格的过滤、数据协调和差异甲基化分析，研究人员在929,009个CpG位点中识别出14个在部分重编程后显著差异甲基化的CpG位点（DMCs）。其中，两个CpG位点（cg20011459和cg01397139）在所有数据集中均一致存在，另外两个CpG位点存在于所有基因组和处理条件下但未在Olova等人的数据集中检出，其余十个CpG仅存在于大鼠和小鼠基因组中。功能富集分析表明，这些共享差异甲基化位点相关的基因富集于细胞周期调控、端粒维持、发育模式和脂质代谢等生物学过程。注释显示，这些位点位于启动子、CpG岛、基因体、外显子、内含子和基因间区，涉及多个基因，如HOXA10、Irx3、Ptch1、ACSL6、Glcci1、Tmem26、FBXL7、Rnf19b、NABP2、Ptprm等。进化保守性分析发现，58.3%的注释CpG位点位于中等至高度保守的基因组区域。转录因子结合富集分析未发现与重编程相关因子（如KLF4、MYC、NANOG）的显著富集。该研究表明，部分重编程诱导了有限但可重复的甲基化变化，提示其抗衰老和再生效应可能通过选择性靶向保守的调控位点实现。该论文发表在《Epigenomes》。
**主要技术方法**
研究人员使用了R语言（v.4.5.2）进行数据分析。从Gene Expression Omnibus（GEO）数据库检索了五个公开DNA甲基化数据集（GSE252401、GSE142439、GSE147436、GSE190665、GSE54848），这些数据集采用Illumina HumanMethylation450 BeadChip (HM450K)、MethylationEPIC BeadChip (EPIC)或Horvath's pan-mammalian methylation array (Mammal40k)平台。数据预处理包括将β值转换为M值，通过ε因子调整极端值，并采用基于数据集来源的协变量进行批次校正。差异甲基化分析使用limma框架的线性模型和经验贝叶斯调节，分别采用并集（保留在所有数据集中出现的CpG，剔除>20%缺失值的位点，最终37,554个CpG）和交集（仅保留五个数据集共有的CpG，无缺失值的5,306个CpG）两种策略。基因组注释基于相应平台注释包，进化保守性分析使用UCSC Genome Browser的PhyloP分数，转录因子富集分析采用ChIP-Atlas工具，针对共享CpG位点±1 kb区域进行。
**研究结果**
**2.1 数据集特征**
五个数据集覆盖人类、小鼠和大鼠样本，包括真皮成纤维细胞、内皮细胞、神经元、海马组织等，处理方案从短期OSK表达（5天）到长期循环OSKM（7-10个月）不等。经筛选后保留189个样本用于分析。
**2.2 跨平台整合二次分析识别出有限的差异甲基化**
通过limma差异甲基化分析，在并集分析中识别出14个显著CpG位点（adj.p.Val < 0.05; |log2FC| > 0.5），其中8个高甲基化，6个低甲基化。交集分析仅保留两个CpG（cg20011459和cg01397139）仍显著。热图和箱线图显示处理组与对照组之间清晰的甲基化模式分离，且方向一致。
**2.3 基因组分布、功能注释与进化保守性**
14个CpG位点分布于多条染色体，不聚集于单一区域。注释显示位于启动子、CpG岛、基因体等区域。功能富集分析表明共享位点相关基因富集于细胞周期、端粒维持、发育模式和脂质代谢。进化保守性分析显示58.3%的CpG位点位于中等至高度保守区域（PhyloP > 1）。
**2.4 转录因子结合富集分析**
ChIP-Atlas分析未发现与重编程相关因子（如KLF4、MYC、NANOG）的显著富集，仅观察到孤立重叠，且未达到统计学显著性。
**讨论与结论**
本研究通过跨物种整合分析，发现部分重编程诱导了有限的甲基化变化，但其中保守的CpG位点可能代表了重编程过程的核心表观遗传特征。这些位点相关的基因涉及发育可塑性、代谢稳态、DNA修复和蛋白稳态等过程，提示部分重编程可能通过调控这些保守位点来发挥抗衰老和再生效应。尽管进化保守性和功能注释支持这些位点的潜在生物学相关性，但转录因子富集分析未提供直接证据表明这些位点位于重编程相关因子的结合区域内。研究局限性包括：显著CpG位点数量较少，可能受效应量异质性、跨数据集变异和平台间CpG覆盖不完全等因素影响。研究结论为：部分重编程通过选择性靶向保守的调节位点（而非全局甲基化重塑）发挥作用，这些位点可作为未来系统性抗衰老治疗中监测部分重编程效果的生物标志物，并可能指导优化因子组合和处理方案，以提高疗效并降低风险。然而，这些发现仍需进一步的实验验证。
