《GeroScience》：Age-associated DNA methylation patterns in mouse blood and tail: feasibility of tissue-specific epigenetic clock development
编辑推荐：
衰老伴随着DNA甲基化（DNA methylation, DNAm）的广泛重塑，这可用于构建表观遗传时钟。在小鼠中，大多数DNAm时钟依赖于内脏器官或终末组织，限制了纵向研究。在此，研究人员比较了匹配的血液和尾部——一种易于获取但尚未充分探索的组织——中与年龄
衰老伴随着DNA甲基化（DNA methylation, DNAm）的广泛重塑，这可用于构建表观遗传时钟。在小鼠中，大多数DNAm时钟依赖于内脏器官或终末组织，限制了纵向研究。在此，研究人员比较了匹配的血液和尾部——一种易于获取但尚未充分探索的组织——中与年龄相关的DNAm模式，并评估了二者作为组织特异性DNAm年龄预测底物的可行性。从7至103周龄的雄性C57BL/6N Per2::Luc小鼠中收集匹配的血液和尾部样本，并使用Infinium小鼠甲基化BeadChip进行图谱分析。经过SeSAMe预处理后，探针水平回归（β ~ 年龄，BH-FDR < 0.05）在血液中鉴定出78,111个与年龄相关的CpG位点，在尾部中鉴定出38,348个。血液主要表现为低甲基化变化和较晚、更突然的年龄转变，而尾部则表现出近乎平衡的高/低甲基化模式并伴有更平滑的转变。片段水平分析（BH-FDR < 0.05，≥ 2个CpG）在血液中鉴定出31,322个与年龄相关的差异甲基化区域（differentially methylated regions, DMRs），在尾部中鉴定出21,883个。方向感知重叠产生了6,229个共享的年龄相关DMR，其特征为CpG岛富集的高甲基化和开放海（Open Sea）富集的低甲基化，并与富含发育、神经、免疫和信号通路的基因相关联。弹性网络（Elastic-net）时钟在训练队列中实现了高准确性（血液：α = 0.8，72个CpG，平均绝对误差（MAE）= 2.30周，R2 = 0.99；尾部：α = 0.5，192个CpG，MAE = 2.75周，R2 = 0.98）。血液时钟推广至一个独立的雌性血液数据集（MAE = 13.48周，R2 = 0.76）。总之，这些发现为使用微创采样进行小鼠组织特异性表观遗传年龄估计提供了概念验证基础，并强调了尾部作为无需终末组织采集的纵向DNAm衰老研究的有前景底物。
**论文解读：小鼠血液与尾部组织特异性表观遗传时钟的可行性研究**
**研究背景、存在问题与研究目的**
衰老是一个时间依赖性的生物学过程，伴随全身性的功能衰退，是全球人口老龄化背景下重大的社会和生物医学挑战。在衰老的生物学标志中，表观遗传改变被视为核心且可调控的组分，其中DNA甲基化（DNA methylation, DNAm）因化学稳定性、可定量测量以及与基因调控的紧密关联，成为生物标志物开发的实用底物。表观遗传时钟通过学习选定CpG位点DNAm水平的加权组合来估算实足年龄或衰老相关表型，已在人类和小鼠中广泛建立。然而，现有的大多数小鼠DNAm时钟依赖于内脏器官或终末组织采集，这限制了在同一动物体内进行纵向衰老研究。血液虽可常规采样并支持稳健的年龄预测，但其他可重复获取的微创组织尚待探索。尾部活检常用于基因分型，可在不处死动物的前提下连续采集，但其与年龄相关的DNAm重塑特征以及尾部特异性表观遗传时钟的可行性尚不清楚。为此，该研究旨在系统性比较匹配的小鼠血液和尾部样本中与年龄相关的DNAm模式，评估各组织作为DNAm年龄预测底物的潜力，并开发组织特异性表观遗传时钟。该论文发表在《GeroScience》。
**主要关键技术与方法**
研究人员使用了来自C57BL/6N Per2::Luc敲入雄性小鼠（年龄范围7–103周）的配对血液和尾部样本（每个年龄点n=2，每个组织共10个样本）。采用Infinium小鼠甲基化BeadChip（MM285平台，约296,000个CpG位点）进行全基因组DNAm图谱分析，并通过SeSAMe软件包进行预处理和质量控制。探针水平的年龄关联分析采用线性回归（β ~ 年龄），并用Benjamini-Hochberg假发现率（BH-FDR）进行多重检验校正。差异甲基化区域（differentially methylated regions, DMRs）通过SeSAMe的片段框架鉴定，要求区域包含至少2个CpG且BH-FDR < 0.05。方向感知的重叠策略用于鉴定血液和尾部共享的年龄相关DMR。组织特异性表观遗传时钟通过弹性网络（Elastic-net）回归（使用glmnet包）在log转换的年龄结果上训练，并通过留一法交叉验证（LOOCV）进行调优。血液时钟在独立雌性小鼠血液数据集（GSE200527）中进行外部验证，并应用于高脂饮食诱导肥胖及饮食逆转干预队列以评估生物学响应性。此外，研究还利用IVIS成像系统监测了Per2::Luc生物发光节律以确认队列的年龄相关功能衰退。
**研究结果**
**组织特异性年龄相关DNA甲基化模式**
通过探针水平线性回归，在血液和尾部中分别鉴定出78,111和38,348个年龄相关CpG（BH-FDR < 0.05）。血液中90.2%的显著CpG呈低甲基化变化，表现出以晚期加剧为主的低甲基化优势模式；尾部则呈现近乎平衡的高/低甲基化分布（50.9%高甲基化 vs 49.1%低甲基化），且年龄转变更平滑。时序距离分析显示，血液的方向性变化峰值出现在50–70周（晚期间隔），而尾部的最大方向性转变出现在30–50周（较早间隔），提示两种组织编码年龄信息的动力学存在差异。
**共享年龄相关DMR揭示保守的区域性DNAm重塑**
片段水平分析在血液和尾部中分别鉴定出31,322和21,883个年龄相关DMR。方向感知重叠产生6,229个共享的年龄相关DMR，其中87.8%为低甲基化。共享高甲基化DMR显著富集于CpG岛（36.8%），而低甲基化DMR主要集中在开放海区域（78.6%），呈现“CpG岛高甲基化/开放海低甲基化”模式。基因特征注释表明共享DMR高度富集于启动子近端TSS200区域。
**共享DMR相关基因的功能富集**
功能富集分析显示，共享高甲基化DMR相关基因富集于发育和神经元系统相关的生物学过程及Reactome通路（如跨化学突触传递），并显著富集ZF5、BEN结构域因子和CPBP/Kid3相关转录因子基序；共享低甲基化DMR相关基因则富集于信号转导、Rho GTPase相关通路及刺激响应过程。这些结果表明保守的发育调控程序与广泛信号网络在跨组织衰老DNAm重塑中共同参与。
**组织特异性表观遗传时钟的构建与验证**
弹性网络回归通过LOOCV选择最优混合参数（血液α=0.8，尾部α=0.5），在完整训练集上重拟合后分别得到包含72个CpG（血液）和192个CpG（尾部）的稀疏时钟，在训练集中表现出高准确性（血液：MAE=2.30周，R2=0.99；尾部：MAE=2.75周，R2=0.98）。血液时钟在独立雌性小鼠血液数据集（GSE200527，n=23，11–117周）中仍与实足年龄显著相关（MAE=13.48周，R2=0.76），支持跨队列迁移性。所有时钟CpG均完全来源于对应组织中年龄相关的CpG集合。
**表观遗传时钟反映热量限制干预的抗衰老效应**
将血液时钟应用于高脂饮食诱导肥胖（diet-induced obesity, DIO）及饮食逆转（DIORe）干预队列（实足年龄69周），结果显示：对照组（CON）预测年龄为70.5周，与实足年龄吻合；DIO组预测年龄显著升高至113.7周，表明持续性高脂饮食诱导的表观遗传年龄加速；DIORe组预测年龄为92.2周，部分逆转了肥胖诱导的年龄加速，提示该时钟可捕获由代谢状态驱动的生物学有意义变异。
**讨论与结论**
讨论部分指出，该研究为概念验证性探索，旨在建立尾部作为微创底物用于表观遗传时钟开发的可行性。血液和尾部在全局甲基化结构和衰老轨迹上可清晰分离，但方向性差异显著：血液以低甲基化主导并呈现晚期阶段性转变，尾部则呈平衡的双向重塑且转变更早、更平滑。这种差异可能源于血液中免疫细胞组成的变化以及尾部作为屏障组织的环境暴露历史。共享DMR的“CpG岛高甲基化/开放海低甲基化”模式与既往研究一致，且富集于发育和信号转导通路，提示保守的衰老表观遗传程序。局限性包括训练样本量较小（每个组织n=10）、仅纳入雄性小鼠、以及MM285阵列覆盖的CpG子集可能遗漏部分年龄相关位点。尽管存在这些限制，研究结果支持血液和尾部作为互补的非致死底物用于小鼠表观遗传衰老研究，并鼓励未来开展更大规模、性别均衡、纵向及干预设计的研究以阐明组织身份和暴露背景如何调节DNAm重塑及时钟性能。
**研究结论**：总之，这些发现为使用微创采样进行小鼠组织特异性表观遗传年龄估计提供了概念验证基础，并强调了尾部作为无需终末组织采集的纵向DNAm衰老研究的有前景底物。
