17种人类组织中DNA甲基化(DNA methylation, DNAm)衰老特征信号的荟萃分析(Meta-analysis)

《Nature Aging》：Meta-analysis of DNA methylation aging signatures in 17 human tissues

　　

表观遗传改变特别是DNA甲基化（DNA methylation, DNAm，胞嘧啶碱基上加甲基基团）随年龄在不同组织中累积，但这些变化是否遵循跨器官的一致模式仍不清楚。研究人员通过对涵盖17种组织的超过15,000份人甲基化谱进行荟萃分析（meta-analysis），发现衰老产生保守及组织特异性的表观遗传特征信号。研究人员鉴定出甲基化水平的系统性偏移、甲基化变异性的增加及跨组织分子无序度（molecular disorder）的增长。网络分析揭示紧密连接的基因簇不受有益干预修饰，而与NAD+代谢相关的可修饰簇支持NAD+作为衰老潜在治疗靶点。编码细胞黏附蛋白的PCDHGA1基因作为跨组织保守枢纽（hub）浮现，暗示细胞间通讯通路参与多器官衰老。该甲基化图谱（methylation atlas）为解析人类衰老分子基础及鉴定潜在生物标志物与转化疗法提供资源。

研究背景与立项依据

既往表观遗传衰老研究多聚焦血液样本，特定实体组织研究匮乏，关于年龄相关DNA甲基化（DNA methylation, DNAm）特征在多大程度上跨组织共享或分化这一基本问题尚未阐明。明确保守与组织特异性甲基化重塑模式对理解系统性与器官特异性衰老机制、评估外周血作为远端器官生物老化_proxy具有重要意义。本研究通过整合公开与协作数据集构建17种成人终身期人类组织DNAm全景图谱，系统表征年龄相关差异甲基化位点（Differentially Methylated Positions, DMPs）、可变甲基化位点（Variably Methylated Positions, VMPs）及DNAm Shannon熵（Shannon entropy）动态，并借助加权基因共表达网络分析（Weighted Gene Co-expression Network Analysis, WGCNA）、多对比通路富集（mitch）及计算机模拟扰动解析衰老关联关键基因与通路。

主要关键技术方法

研究人员整合来自GEO、ArrayExpress、dbGaP、EGA及合作者的131个独立数据集共15,995例人组织DNA甲基化（DNAm）样本（Illumina 27k/450k/EPIC阵列），涵盖17种组织，排除癌症样本及低质量探针后按标准化ChAMP流程预处理并用ComBat校正批次效应，部分组织用EPISCORE/EpiDISH估计细胞类型比例作协变量。① DMP鉴定：各数据集内以limma拟合年龄多元线性回归（M值作因变量，校性别/BMI/批次），bacon校正后METAL行逆方差固定效应荟萃分析（FDR<0.005）；② VMP鉴定：两步Breusch–Pagan异方差框架检测年龄关联甲基化残差方差变化后荟萃（FDR<0.005）；③ Shannon熵：计算协变量校正后全基因组及各CpG子集熵值并建模年龄依赖；④ 跨组织DMP/VMP Jaccard重叠与方向一致性及置换检验；⑤ WGCNA构建脂肪、血、脑、骨骼肌组织特异性共甲基化模块并做GO富积；⑥ mitch做跨组织基因水平t统计量多对比方向性通路富集；⑦ 基于t统计量构建基因–基因Spearman相关跨组织模块及计算机模拟扰动（in silico perturbation）评估模块脆弱性。

研究结果

A pan-tissue landscape of age-associated DNA methylation

研究人员整合17种组织逾15,000样本建立跨横断面DNAm图谱，全局平均甲基化率38%（宫颈）至63%（视网膜）。在单数据集内鉴定DMP（多元线性回归校协变量）、VMP（Breusch–Pagan异方差框架）及Shannon熵（全基因组及各子集），三者互补刻画位点均值偏移、方差偏移及系统无序度。

Multivariate linear analysis of age reflects distinct tissue-specific aging signatures

脑、肝、肺DMP数最多（脑>180,000），胰腺、视网膜、前列腺极少（FDR<0.005）；除骨骼肌与肺外多数组织年龄关联DMP以高甲基化（hypermethylation）为主，青年低/中甲基化CpG随龄向高甲基化漂移；脑DMP高/低甲基化均衡分布。DMP检出数与年龄跨度显著正相关（R=0.66, P=0.0072），受样本量及效应量影响，少数组织未检出主因统计功效不足。

VMPs are scarce and highly tissue-specific

仅脑、颊黏膜、宫颈、肺、皮肤有较多年龄关联VMP，其余组织极少或无；跨组织VMP重叠极少（任一对<100），经置换检验仅脑–颊、颊–宫颈等少数组对超出随机期望且多为同向变异性改变。VMP检出与年龄范围/样本量无显著相关，支持其生物学组织特异性而非技术假象。

Shannon entropy and epigenetic disorder reveal tissue-specific aging patterns

Shannon熵在多种组织随龄升高，脂肪/乳腺/肾等熵升高优先发生于DMP位点，提示老化关联位点表观遗传保真度丧失；脑/颊/肝/肺/皮肤熵升高出现于DMP–VMP重叠位点，反映方向性与变异性均紊乱的更随机表观遗传衰老模式。

WGCNA highlights tissue-specific hubs and conserved aging modules

WGCNA鉴定组织特异性共甲基化模块，与年龄正/负相关模块分别富积发育黏附通路与突触/NF-κB/代谢调控通路。原钙粘蛋白γ（protocadherin gamma, PCDHG）基因家族在脂肪、血、脑、骨骼肌中反复作为模块影响力基因（module-influential genes，高模块内连通性kME且强年龄关联）出现。跨组织保守富集同型细胞黏附（homophilic cell adhesion via plasma membrane adhesion molecules）及小GTP酶信号转导（small GTPase-mediated signal transduction）。

Cross-tissue overlap reveals core and divergent epigenetic aging programs

DMP两两组织重叠Jaccard指数0%–25%，但重叠位点甲基化方向多数一致（>59%），仅脑–视网膜方向一致性低（28%）。ELOVL2位点cg16867657在15种组织显著。经置换检验多组织DMP重叠显著超越随机。≥9种组织共有DMP富积转录调控/染色质组织/基因调控通路；mitch整合分析确认转录调控超甲基化及线粒体和大分子代谢相关低甲基化为跨组织共性。

Cross-tissue network level analysis through module enrichment identifies key pathways and genes as module-influential genes in aging

基于mitch t统计量建基因–基因相关矩阵并分层聚类得跨组织共老化行为基因模块，仅模块5（module 5）显示全组织一致高甲基化，其枢纽为PCDHG家族PCDHGA1，提示该家族参与各组织衰老程序。多数模块呈组织特异性老化甲基化方向，部分模块在不同组织方向相反。

In silico robustness analysis reveals functionally fragile genes driving module architecture

计算机模拟上调/下调单个基因甲基化t统计量评估模块结构扰动，PCDHG、MEST、HDAC4及HOX家族具最高跨模块扰动评分，多数模块扰动恶化衰老关联特征（脆弱），仅富集NAD⁺补救合成（nicotinamide riboside salvage pathway）的模块可被正向调节（韧性），支持NAD⁺代谢为潜在干预节点。

讨论与结论翻译

本荟萃分析涵盖17种人组织逾15,000甲基化组的三个互补层面（DMP、VMP、Shannon熵），表明衰老关联甲基化组改变兼具组织特异与跨组织保守效应，反映随机变异与结构化组织依赖重塑并存而非单纯随机漂移。跨组织与计算机模拟分析鉴定衰老关联通路模块层级脆弱性，并突出以NAD⁺生物合成为代表的韧性模块具治疗相关性。总结数据与结果公开于 https://bioinformatics.erc.monash.edu/apps/human-aging-atlas/。局限含组织间样本量不均、批量组织掩盖细胞类型特异信号、横断面设计无法直接捕捉个体内时序轨迹、缺染色质可及性或表达等多组学验证，需未来纵向与大样本细化。