《Frontiers in Nutrition》：EWAS in a polyphenol dense, DNA methylation-targeted, controlled diet and lifestyle study
编辑推荐：
背景：饮食和生活方式因素可影响DNA甲基化模式。研究人员先前报告了在健康男性中进行的为期8周的多酚密集型、DNA甲基化靶向饮食和生活方式干预（甲基化饮食与生活方式研究）后，表观遗传年龄的衰减，其中植物营养素/多酚丰富的食物（绿茶、乌龙茶、姜黄素、大蒜和浆果）对
背景：饮食和生活方式因素可影响DNA甲基化模式。研究人员先前报告了在健康男性中进行的为期8周的多酚密集型、DNA甲基化靶向饮食和生活方式干预（甲基化饮食与生活方式研究）后，表观遗传年龄的衰减，其中植物营养素/多酚丰富的食物（绿茶、乌龙茶、姜黄素、大蒜和浆果）对这一效果最具预测性。方法：在此，研究人员在来自甲基化饮食与生活方式研究的38名参与者中开展了一项表观基因组关联研究（EWAS）。干预措施包括一种有意富含甲基化途径底物和辅因子营养素以及已知可改变DNA甲基转移酶（DNMT）酶活性的成分的饮食模式。与先前样本量较小、不太可能发现FDR显著结果的EWAS研究一致，研究人员在探索性分析中使用了预先设定的名义P值阈值（0.001、0.0001）。结果：在P < 0.001（未校正）水平下，干预组中鉴定出676个差异甲基化位点（DML），而对照组为286个。在P < 0.0001（未校正）水平下，干预组中鉴定出50个DML，而对照组为13个。其中15个DML位于基因的转录起始位点近端区域，这些基因涉及锌稳态和营养感知、发育与多能性、蛋白质稳态与基因组稳定性、肿瘤抑制以及突触功能。一项分组-时间交互作用分析在P < 0.0001水平下鉴定出70个干预特异性DML，名义富集包括自噬、mTOR信号传导和染色质重塑通路。一项区域DMR分析鉴定出128个组内和129个交互特异性DMR。DMR功能富集分析揭示与脂质代谢（α-亚麻酸、硫辛酸、不饱和脂肪酸生物合成、PPAR信号传导、胆固醇稳态）、中枢能量代谢（TCA循环、糖酵解/糖异生、磷酸戊糖途径、丙酮酸）和营养感知（PI3K-Akt、mTOR、AMPK、自噬以及其他通路）的趋同名义关联。正如预期，由于队列规模有限和干预时间短，单个CpG发现或富集分析均未通过多重检验校正，因此仅被视为探索性和产生假设的结果。结论：本EWAS在干预组中发现了比对照组更大数量的名义上变化的CpG，以及生物学上一致的甲基化变化。这些发现为先前观察到的表观遗传年龄衰减提供了机制假设。在更大队列、更长干预时间中的重复以及功能验证仍然至关重要。
**论文解读：多酚密集型、DNA甲基化靶向饮食与生活方式干预的表观基因组关联研究**
**研究背景与问题**
饮食和生活方式因素可通过影响DNA甲基化模式，进而调节基因表达与衰老相关疾病风险。表观遗传改变（如DNA甲基化）被视为连接环境暴露与衰老及疾病风险的重要机制，部分变化可通过基于DNA甲基化的“表观遗传时钟”捕捉，用于预测年龄加速相关的发病率和死亡率。然而，现有研究在饮食暴露定义、多因素干预复杂性、个体间基因组或微生物组差异等方面存在精度不足，导致差异甲基化CpG位点缺乏一致性。研究人员此前在一项随机对照试验中，对健康中年男性（50–72岁）实施为期8周的DNA甲基化靶向饮食与生活方式干预，发现干预组唾液表观遗传年龄（基于Horvath 2013时钟）平均减少3.23年，且多酚丰富的食物（绿茶、乌龙茶、姜黄素、大蒜、浆果）是这一效果的最强预测因子。为超越全局表观遗传年龄估计，深入探索干预引起的具体通路、基因集和CpG位点变化，本研究在此基础上开展了二次表观基因组关联研究（EWAS）分析。
**研究内容与结论**
本研究基于原始试验中38名参与者（干预组18人，对照组20人）的配对唾液样本（基线和第8周），利用Illumina Methylation EPIC阵列进行甲基化谱分析，通过三类对比分析（组内干预组、组内对照组、分组-时间交互作用）鉴定差异甲基化位点（DML）和差异甲基化区域（DMR）。主要结论包括：在名义P<0.001阈值下，干预组鉴定出676个DML，对照组为286个；在P<0.0001阈值下，干预组50个DML，对照组13个。其中15个转录起始位点近端DML涉及锌稳态与营养感知（SLC39A7/RXRB、MTF1）、发育与多能性（JARID2、PAX6）、蛋白质稳态与基因组稳定性（UBE2R2、BCLAF1、SSBP2）、肿瘤抑制及突触功能（KCNS3、IGSF11）等关键基因。分组-时间交互作用分析鉴定出70个干预特异性DML，名义富集于自噬、mTOR信号和染色质重塑通路。DMR分析（HMFDR<0.05）确认128个组内和129个交互特异性DMR，功能富集显示与脂质代谢（α-亚麻酸、硫辛酸、PPAR信号）、中枢能量代谢（TCA循环、糖酵解）和营养感知（PI3K-Akt、mTOR、AMPK）通路趋同。然而，所有结果均未通过FDR校正，需视为探索性假设。
**研究意义**
本研究首次在对照饮食干预中系统揭示多酚密集型、甲基化靶向方案对唾液甲基组的广泛影响，为先前观察到的表观遗传年龄衰减提供了机制性假设，提示多酚可能通过调节营养感知、脂质代谢和自噬等保守衰老通路发挥表观遗传调控作用。论文发表在《Frontiers in Nutrition》。
**主要关键技术方法**
研究人员从原始试验（ClinicalTrials.gov NCT03472820）中采集38名参与者（干预组18人，对照组20人）的配对唾液样本（基线、第8周），使用Illumina Methylation EPIC阵列（~866,000 CpG位点）进行甲基化分析，通过minfi包进行ssNoob背景校正和染料偏差归一化。采用limma框架建立多变量线性模型，纳入芯片批次、唾液细胞类型比例（上皮细胞、成纤维细胞、免疫细胞，通过EpiDISH评估）、技术探针主成分和个体固定效应。使用预先指定的名义P值阈值（0.001, 0.0001），通过missMethyl包进行KEGG和MSigDB Hallmark通路富集分析，并利用DMRcate（核带宽1000 bp，C=2，区域FDR通过HMFDR校正）鉴定DMR。交互作用分析类似，并进行了细胞组成敏感性分析（排除5名细胞比例变化最大参与者后，logFC相关性Pearson r=0.98）和交叉研究比较（与Perlmutter等和Dwaraka等EWAS比较CpG重叠）。
**研究结果**
**组内差异甲基化分析**
在名义P<0.001阈值下，干预组鉴定出676个DML（290个低甲基化，386个高甲基化），对照组286个；P<0.0001阈值下，干预组50个DML（27个低甲基化，23个高甲基化），对照组13个。基因组膨胀因子λ分别为0.86（干预组）和1.04（对照组），无异常膨胀或压缩。
**转录起始位点近端DML探索**
在P<0.0001阈值下，50个DML中15个位于转录起始位点近端区域（TSS200、TSS1500），对应15个基因，其中12个低甲基化、3个高甲基化。这些基因包括：锌转运蛋白SLC39A7（与RXRB共享启动子）、金属响应转录因子MTF1、锌指转录因子ZNF503、多梳抑制复合物2（PRC2）辅助因子JARID2、发育转录因子PAX6、肿瘤抑制因子SSBP2、泛素结合酶UBE2R2、内质网/高尔基体转运伴侣BCAP29、SMN复合体组分GEMIN6、DNA损伤/衰老整合因子BCLAF1、SNARE蛋白STX2、钾通道亚基KCNS3、免疫球蛋白超家族黏附分子IGSF11、原钙黏蛋白PCDH19、驱动蛋白轻链KLC2。效应量（logFC）范围为-0.0156至+0.0178，方向与生物学功能一致。
**组内功能富集分析**
对50个DML（P<0.0001）的KEGG富集显示名义显著通路：硫辛酸代谢（P=0.034）、α-亚麻酸代谢（P=0.049）、年轻发病的糖尿病（P=0.051）、TCA循环（P=0.055）、干细胞多能性调控（P=0.057，含PAX6和JARID2）等。Hallmark基因集富集显示胰腺β细胞（P=0.082）和胆固醇稳态（P=0.152）名义富集。对676个DML（P<0.001）的KEGG富集强化了脂质代谢主题：α-亚麻酸代谢（P=0.0031）、醚脂代谢（P=0.038）、脂肪酸代谢（P=0.031）、PPAR信号（P=0.049）、不饱和脂肪酸生物合成（P=0.050）；PI3K-Akt信号通路（18个基因，P=0.041）和叶酸转运与代谢（P=0.042）也名义富集。Hallmark分析中，胆固醇稳态（7个基因，P=0.0051，FDR=0.26）为最接近FDR校正的基因集，其次为UV响应下调（P=0.018）和胰腺β细胞（P=0.022）。糖酵解/糖异生、TCA循环、磷酸戊糖途径、碳代谢等核心能量代谢通路也现名义富集。
**分组-时间交互作用分析**
在P<0.001阈值下，694个CpG呈交互效应；P<0.0001阈值下，70个CpG（36个正交互，34个负交互），其中18个（25.7%）与组内DML重叠。对694个交互CpG的Hallmark分析显示E2F靶标（11个基因，P=0.013）和MYC靶标V2（4个基因，P=0.045）名义富集；ROS通路（P=0.062）和WNT/β-catenin信号（P=0.080）在P<0.1水平名义富集。对70个最显著交互CpG的KEGG富集显示自噬-动物（SQSTM1/p62、TSC2、VPS39，P=0.014）和ATP依赖性染色质重塑（GATAD2A、INO80E、SMARCA4，P=0.050）名义富集；Hallmark分析中PI3K/AKT/mTOR信号（SQSTM1、TSC2，P=0.046）名义显著。mTOR信号和泛素介导的蛋白水解在P<0.1水平额外富集。
**差异甲基化区域（DMR）分析**
DMRcate分析在组内（干预组）鉴定出128个DMR（HMFDR<0.05），包括SLC39A7/RXRB（9 CpG区域）、STX2（5 CpG）、ZNF503（6 CpG）和RBMS1（4 CpG）。交互作用分析鉴定出129个DMR，包括VPS39（15号染色体，9 CpG）、INO80E/HIRIP3（16号染色体，18 CpG）和STX2（12号染色体，3 CpG），与单CpG顶部结果一致。
**细胞组成敏感性分析**
排除5名细胞比例变化最大参与者后（n=13保留），676个DML的logFC与全模型高度相关（Pearson r=0.98），方向100%一致，表明结果不受细胞组成漂移影响。
**交叉研究比较**
与两项独立饮食EWAS（Perlmutter等：90天苦荞多酚补充剂，887 DML；Dwaraka等：8周纯素饮食，607 DML）比较，本研究的694个交互CpG仅与Perlmutter共享1个CpG（cg03757262，低甲基化），与Dwaraka共享1个CpG（cg08127021在LRP5，本研究中高甲基化，Dwaraka中低甲基化方向相反），无三向重叠，表明饮食甲基化信号具有高度干预和背景特异性。
**讨论与结论**
讨论部分指出，本研究作为探索性二次分析，在名义阈值下观察到干预组较对照组更多的甲基化变化，且涉及营养感知、脂质代谢、自噬等生物学一致通路，为先前表观遗传年龄衰减提供了机制假设。但所有结果均未通过FDR校正，需谨慎解释，可能源于样本量小（n=38）、干预时间短（8周）及效应量微弱。单CpG结果的生物学关联性（如SLC39A7、MTF1、JARID2、PAX6等）需通过区域分析、重复验证和功能实验确证。局限性包括：样本量不足、细胞组成偏差（虽经敏感性分析部分缓解）、组织特异性（唾液不能代表肝脏、脂肪等靶组织）、缺乏基因表达和蛋白质功能数据，且多模态干预无法归因于单一成分（如多酚）。与DIRECT PLUS试验（n=256，18个月，绿色地中海饮食）相比，后者在FDR水平鉴定出1573个DMR，提示更大样本和更长干预可能获得更稳健的表观基因组效应。
**研究结论翻译**
在这项探索性二次EWAS（n=38；18例治疗组，20例对照组）中，为期8周的多酚密集型、DNA甲基化支持性饮食和生活方式干预，与对照组相比，与更大数量的名义上个体内甲基化变化相关。在最高排名位点中，转录起始位点近端CpG映射到与衰老生物学具有潜在相关性的基因，包括营养/金属感知和应激响应调控、发育与多能性，以及蛋白质稳态和基因组维护相关通路。一项补充的分组-时间交互作用分析鉴定出70个干预特异性DML。富集分析还提示在脂质、核心能量代谢和营养感知通路中的（名义上）聚类。一项补充的DMR分析鉴定出128个组内和129个交互特异性区域（HMFDR<0.05），提供了与单CpG分析中突出显示的多个位点一致的区域共甲基化证据。与其他多酚丰富干预的EWAS进行交叉研究比较显示，CpG水平重叠极小，表明饮食甲基化信号具有干预和背景特异性，这一发现对饮食EWAS的可重复性和普适性具有重要意义。重要的是，所有位点水平和通路水平发现均未通过FDR校正。因此，这些发现应被解释为产生假设，并需要在更大、有充分效力的队列中进行重复，以及多组学后续研究，包括RNA-seq、蛋白质组学和功能读数。
