《Epigenetics》：A breast tissue-specific epigenetic clock provides accurate chronological age predictions and reveals de-correlation of age and DNA methylation in tumor-adjacent and tumor samples
编辑推荐：
表观遗传时钟已被广泛用于估算不同组织中的生物学年龄，但它们在乳腺组织中的准确性仍不理想。Horvath和Hannum等经典模型在预测乳腺组织实足年龄方面表现不佳，凸显了对组织特异性方法的需求。在本研究中，研究人员介绍了一种乳腺组织特异性表观遗传时钟（Breas
表观遗传时钟已被广泛用于估算不同组织中的生物学年龄，但它们在乳腺组织中的准确性仍不理想。Horvath和Hannum等经典模型在预测乳腺组织实足年龄方面表现不佳，凸显了对组织特异性方法的需求。在本研究中，研究人员介绍了一种乳腺组织特异性表观遗传时钟（Breast Tissue-specific Epigenetic Clock, BTEC），该时钟利用来自7项不同研究的553份健康乳腺组织样本的DNA甲基化数据开发而成。BTEC显著优于现有时钟，其与实足年龄的相关性更高（r = 0.91），预测误差更低（MAE = 3.0年），且无需进行数据集特异性的回归调整。值得注意的是，BTEC在不同祖先群体间表现一致，而现有时钟则表现出与祖先相关的可变偏倚。研究人员还发现，常用于校正预测偏倚的基于残差的表观遗传年龄加速（epigenetic age acceleration, EAA）高度依赖于参考数据集的选择，并不能提供可推广的校正。BTEC对肿瘤邻近样本的实足年龄预测出现扭曲，平均偏差为?2.3年；在肿瘤样本中这种扭曲更为明显，预测年龄与实足年龄之间的平均差异为?14至?15年。按分子亚型分析时，与HR+肿瘤相比，侵袭性更强的HER2+和三阴性乳腺癌（TNBC）肿瘤中的扭曲更大。总体而言，研究结果表明乳腺肿瘤通常并不表现出加速的表观遗传衰老。BTEC所使用的探针与已知癌基因如JAK1、TP73、PDGFRA、CCND1和STAT5B相关，也与癌症相关通路（包括Wnt、Hedgehog和ErbB信号通路）中的基因相关。重要的是，在探索性分析中，由BTEC衍生的表观遗传年龄加速（EAA）所测量的极端表观遗传年龄偏差与HR+乳腺癌患者的生存相关，凸显了组织特异性表观遗传衰老指标的预后价值。这些发现表明，BTEC不仅改善了乳腺组织中的年龄预测，还捕获了肿瘤表观遗传衰老中具有生物学意义的改变，具有潜在的临床意义。
该研究发表于《Epigenetics》。研究背景：衰老是多种疾病和癌症的重要危险因素，表观遗传变化与衰老密切相关。表观遗传时钟（epigenetic clock，利用DNA甲基化预测实足年龄的数学模型）已被广泛用于估算生物学年龄，但Horvath、Hannum、PhenoAge和AltumAge等泛组织时钟在乳腺组织中的预测准确性低，r值波动于0.35～0.78，部分模型甚至不如常数预测。已有研究将残差化EAA作为校正手段，但残差化回归参数随参考数据集变化显著，不能提供通用校正；同时这些模型存在祖源相关偏倚。因此，需要开发乳腺组织特异性时钟，以准确评估乳腺衰老及肿瘤中的表观遗传年龄改变。
研究人员整合了来自GEO的13项研究、共2023例乳腺组织DNA甲基化样本，包括553例健康供体、362例肿瘤邻近和1108例肿瘤样本，开发了乳腺组织特异性表观遗传时钟（Breast Tissue-specific Epigenetic Clock, BTEC）。在留一队列交叉验证中，BTEC对健康供体样本的年龄预测相关系数r=0.918（SeSAMe）和0.912（BMIQ），中位绝对误差（MAE）<3.2年，优于所有已测试的泛组织模型，且无祖源偏倚。将BTEC应用于肿瘤邻近和肿瘤样本时，预测年龄系统低于实足年龄，平均EAA分别为?2.358年和?13.938年，提示乳腺肿瘤中甲基化衰老信号与年龄发生解偶联，而非加速衰老。论文发表在《Epigenetics》，其意义在于：组织特异性表观遗传时钟不仅改善了年龄预测，还揭示了肿瘤表观遗传衰老扭曲及潜在的预后价值。
主要技术方法：数据来自GEO的13项研究：553例健康供体、362例肿瘤邻近和1108例肿瘤样本；另用TCGA 44例独立样本验证。用Minfi+BMIQ和SeSAMe（QCDPB）并行处理为β值；取450k/EPIC/EPICv2共有且与年龄线性或二次项（β和β
2）相关top 5%的4713个探针，以弹性网回归（alpha=0.0077）训练BTEC，用留一队列交叉验证评估；用EpiDISH估计细胞组分、stemTOC估算有丝分裂年龄、Enrichr做基因集富集分析、Cox比例风险模型和Kaplan-Meier法做生存分析。
研究结果：
乳腺组织DNA甲基化数据：研究人员整合了13项研究的健康供体、肿瘤邻近和肿瘤样本；肿瘤样本按分子亚型分为HR+（256例）、HER2+（82例）和TNBC（362例）。肿瘤纯度中位数为0.548，显著高于肿瘤邻近（0.276）和健康供体（0.284）。
多器官表观遗传时钟的性能：在553例健康样本上，AltumAge预测最佳（r=0.873），Hannum次之，但所有模型RMSE为8.5～22.7年；Horvath和PhenoAge的误差甚至大于常数模型。残差化EAA的回归截距和斜率均随参考数据集显著变化，说明其不能作为通用校正方法。
乳腺组织特异性表观遗传时钟优于泛组织模型：通过弹性网回归训练BTEC，LOCO验证显示r约0.91、MAE<3.2年；BTEC在Black与White群体间无显著误差差异，也不受肥胖状态影响。
表观遗传时钟在肿瘤邻近样本中的性能：BTEC在362例肿瘤邻近样本和44例TCGA独立样本上均优于其他模型，RMSE显著更低；FFPE保存方式对BTEC无显著偏倚。
肿瘤邻近和肿瘤样本中的表观遗传年龄改变：BTEC预测的EAA在肿瘤邻近（平均?2.358年）和肿瘤（平均?13.938年）中显著低于健康样本；HER2+和TNBC比HR+具有更负的EAA。PAM50分类中Basal样肿瘤EAA低于Luminal A/B。细胞组分分析不能解释亚型相关EAA差异。肿瘤分级、分期与EAA无显著关联，但接受化疗或放疗的TNBC样本EAA更低。
肿瘤表观遗传年龄不由复制速率决定：Ki-67高/低分组与BTEC年龄的关系在不同数据集中不一致；stemTOC预测的有丝分裂年龄与实足年龄、BTEC预测和EAA均无显著相关。
BTEC探针甲基化状态在肿瘤邻近和肿瘤样本中与年龄解偶联：BTEC探针在健康样本中与年龄的平均|r|为0.359/0.321，在肿瘤邻近中降至0.130/0.158，在肿瘤中降至0.078/0.074，说明这些探针在肿瘤中丧失了年龄相关甲基化模式。
BTEC探针相关基因突出癌症通路和模型特异性功能：两个BTEC的探针相关基因集有>73%重叠，但仅有少数已知哺乳动物衰老基因；正系数探针富集于蛋白质翻译、G蛋白信号和癌基因通路扩增，负系数探针富集于细胞黏附、Wnt、Hedgehog、ErbB、AMPK等通路；包含JAK1、TP73、PDGFRA、CCND1、STAT5B等癌基因。
与BTEC探针相关的年龄相关基因表达特征：利用GSE102088表达数据，140个基因表达与年龄相关；其中ELF5、MVP、VGLL4等与乳腺分化、耐药或Hippo肿瘤抑制通路相关。
按表观遗传年龄加速分层的肿瘤转录组特征：基于GSE37754，EAA<0的肿瘤过表达细胞周期和有丝分裂调节基因（CDC20、CCNA2、PLK1、CDK1）及免疫功能基因；EAA>0的肿瘤过表达PI3K-Akt信号基因（NTRK2、IGF1R、AREG、FGF1）及预后标志基因（PGR、SCUBE2）。
BTEC的EAA幅度与长期预后的探索性关联：Cox连续模型未见显著关联，但按EAA阈值分层的log-rank分析在HR+中发现最佳阈值约?6.58年，两组生存差异显著，且bootstrap分析支持该阈值的稳定性（中位?6.575年，95% CI ?28.796～1.936年）；TNBC中未发现稳定阈值。Horvath和AltumAge未显示一致结果。
讨论部分总结：现有泛组织时钟在乳腺组织中误差大且存在祖源偏倚，残差化EAA属于临时性调整而非稳健校正。BTEC训练于健康组织，外推到肿瘤时可能超出训练分布，但多项证据支持其肿瘤EAA信号具有生物学意义：BTEC在肿瘤邻近样本仍保持较好性能、探针去相关呈系统性、亚型差异稳健、与独立研究一致，且EAA与生存相关。负向EAA并不表示肿瘤“生物学上更年轻”，而是年龄相关CpG位点的甲基化向与正常衰老相反的方向偏移。Ki-67和stemTOC结果不支持增殖或复制历史解释；细胞组分只能解释小部分EAA方差。BTEC探针相关基因和转录组特征提示肿瘤去分化、干性获得或表观遗传重编程等机制。探索性生存分析提示极度负向EAA可能与HR+不良预后相关，但需独立队列验证。
研究结论（翻译）：总之，该研究表明，乳腺组织特异性表观遗传时钟（BTEC）相比泛组织模型改善了年龄预测性能，但其主要价值在于提供的生物学和临床见解。在一个大型整合数据集中，BTEC和其他已有时钟均未检测到乳腺肿瘤中加速衰老的一致证据；相反，肿瘤表现出稳健的负向表观遗传年龄加速趋势。这表明甲基化衰老信号与实足年龄和增殖历史发生机制性解偶联，可能反映去分化、干性获得或表观遗传重编程等肿瘤特异性过程。除预测准确性外，BTEC还揭示了不同的转录程序、亚型特异性差异，以及先前时钟未能捕获的HR+和TNBC肿瘤预后模式。这些发现凸显了组织特异性表观遗传时钟的更广泛意义：不仅可作为完善风险分层的生物标志物，也可作为揭示癌症相关衰老扭曲生物学基础的工具。
