《Human Brain Mapping》：Epigenetic and Brain Structure Covariation Patterns in Typical and Atypical Development
编辑推荐：
通过研究不改变DNA序列的基因表达变化，表观遗传学为揭示生物因素如何与环境相互作用以塑造脑发育提供了有前景的线索。然而，既往发育研究多聚焦于候选脑区或候选表观遗传位点，限制了对儿童期全脑与全表观基因组关联的理解。本研究在144名伴或不伴注意缺陷多动障碍（ADH
通过研究不改变DNA序列的基因表达变化，表观遗传学为揭示生物因素如何与环境相互作用以塑造脑发育提供了有前景的线索。然而，既往发育研究多聚焦于候选脑区或候选表观遗传位点，限制了对儿童期全脑与全表观基因组关联的理解。本研究在144名伴或不伴注意缺陷多动障碍（ADHD）的儿童中，联合全脑结构影像（T1加权、扩散加权）与全表观基因组DNA甲基化（DNAm）数据，采用平行独立成分分析（pICA）识别关联模式。研究人员识别出5个显著关联的脑-表观基因组成分，其中甲基化解释了关联脑成分中显著比例的方差，并进一步通过多元随机森林分析在最强表观遗传-脑关联中定位908个CpG位点。这些CpG注释基因富集于神经元发育及精神/神经发育障碍相关通路，凸显功能相关性。探索性分析提示ADHD诊断与人口学变量可能调节这些脑-表观基因组关联。本研究将脑解剖与参与神经元发育及精神障碍的功能相关基因DNAm模式相联系，揭示遗传与环境因素影响脑发育的潜在生物学洞察，并提供表观遗传至脑簇映射及配套提取脚本，为未来神经发育研究提供资源。
研究背景方面，儿童期人脑发育依赖遗传易感与环境线索的高度调控互作，双生与全基因组关联研究（GWAS）证实脑特征具显著遗传率，但遗传率随脑区与发育阶段而异，童年环境影响力上升。DNA甲基化（DNAm）作为核心表观修饰，在神经发生、迁移与突触形成中调控基因表达，并与典型发育及ADHD等神经发育条件相关。既往影像-表观研究受高维数据限制多聚焦先验脑区或候选标记，缺乏全脑-全表观基因组整合。为此，研究人员在《Human Brain Mapping》发表该研究，以数据驱动方式绘制典型与非典型发育中脑结构（体积、纤维密度FD、纤维横截面FC）与DNAm的宽域共变地图，并解析其功能意义与临床调节效应。
主要关键技术方法上，研究基于澳大利亚墨尔本NICAP纵向亚队列，纳入144名9至11岁儿童（73典型发育、71 ADHD），同期采集唾液DNAm（Illumina EPIC 850K阵列）与3T MRI（T1加权及高角分辨率扩散成像）。影像经FreeSurfer与MRtrix处理得体积、FD、FC体素图；甲基化经minfi质控、EpiDISH校正上皮细胞比例与批次效应后取M值。采用Fusion ICA工具箱行平行独立成分分析（pICA）获脑-甲基化共变成分，Bonferroni校正筛5组成对显著关联；成分内取绝对因子载荷>3的脑簇（弃<20体素）与CpG，再借randomForestSRC包“变量狩猎”以最小深度选最强CpG-脑簇关联（908 CpG/748基因）；功能注释用FUMA、SynGO、DisGeNET，临床调节用limma与线性混合模型。
研究结果部分，3.2节平行ICA共变总模式中，研究人员获5组显著配对：1组FC-DNAm（R2=?0.295）、1组FD-DNAm（R2=0.309）、3组体积-DNAm（R2=?0.360、0.335、?0.277），筛选93脑簇与19119唯一CpG。3.3节各成分脑簇描述显示，FC成分负关前额顶叶联络通路（上纵束SLF、钩束、下额枕束IFOF）；FD成分正关长程整合通路（SLF、IFOF、弓状束、丘脑辐射、胼胝体压部）；体积成分1正关广泛整合网络（额中回、颞中回、楔前叶、小脑VI），成分2负关感觉运动与视觉区（中央后回、舌回、旁扣带回），成分3负关局限颞极-楔前叶-枕外侧。3.4节核心表观位点中，cg03743597（RNU11）、cg04016767（ZNF571）高甲基化总体伴更高结构指标，cg09860529（WDR25）呈反向；三探针跨全部5成分出现。3.5节随机森林最强关联获908 CpG注748基因，右外侧枕叶体积预测子最密（46），cg09507771（NBPF1）等跨8脑簇；3.5.1节功能映射示基因脑高表达，富集模式规范、器官形态发生、染色质、树突、突触后特化（SynGO），DisGeNET示685项精神障、477项神经系统病关联。3.6节人口临床变异中，ADHD组FD成分因子载荷更低（T=?3.981，padj=0.006），左SLF的FD较低但未过全协变量校正；性别差异示女性额中回与左枕外侧体积更大，9个性别差异甲基化位点；年龄无直接效应。3.7节调节分析示ADHD诊断调节胼胝体压部与左SLF的FD-甲基化（RTCD1、HTR6）关联，仅在对照组显著；性别调节右SLF FD-ZNF25（女）及左旁扣带回体积-cg01298254（男）；年龄增强HTR6、ERG相关关联。
讨论与结论翻译部分，研究人员指出五组跨FC、FD、体积模态的显著共变表明DNAm与脑宏观缩放及微结构细化多样关联，而非孤立；白质SLF、IFOF、皮质脊髓束映射_fronto-parietal及半球间成熟，灰质体积成分1涵盖默认网络节点楔前叶与小脑VI，契合晚期童年皮质-小脑整合“微调”；RNU11失活致小鼠小头、WDR25涉早期胚胎发育、ZNF571属C2H2锌指转录家族调控神经干细胞，佐证跨成分一致性。908 CpG注释基因富集神经元发育与精神/神经发育疾障（双相、精神分裂、ASD、智力障碍），共享遗传架构提示典型与非典型发育通路重叠。综上，本研究以pICA与随机森林在144名9至11岁ADHD/典型儿童中建立全脑-全表观基因组共变框架，锁定5成分与908 CpG（748基因），阐明甲基化解释脑结构方差的比例，功能注释验证神经元发育与精神障碍相关，探索性示ADHD诊断、性别、年龄调节关联强度，并提供表观-脑映射数据库与查询脚本，为后续神经发育机制发现提供非侵入、可复用资源。
