《Psychiatry and Clinical Neurosciences》：Genome-wide DNA methylation study uncovers epigenetic markers associated with early life stress in individuals with major depressive disorder
编辑推荐：
目的：早期生活压力（early life stress, ELS）是重性抑郁障碍（major depressive disorder, MDD）的重要风险因素；然而，该关联背后的表观遗传机制仍知之甚少，部分原因在于既往研究缺乏按创伤分层设立的对照组。
方法：
目的：早期生活压力（early life stress, ELS）是重性抑郁障碍（major depressive disorder, MDD）的重要风险因素；然而，该关联背后的表观遗传机制仍知之甚少，部分原因在于既往研究缺乏按创伤分层设立的对照组。
方法：研究人员采用简化代表性亚硫酸氢盐测序（Reduced-Representative Bisulfite Sequencing, RRBS）技术，对90名成人全血样本进行了全基因组DNA甲基化谱分析。受试者分为四组：伴ELS的MDD（n=30）、不伴ELS的MDD（n=30）、伴ELS的对照组（n=15）及不伴ELS的对照组（n=15）。差异甲基化区域（differentially methylated regions, DMRs）通过MethylKit和DMRcate软件识别（≥5个CpG位点，错误发现率FDR<0.05），并利用FUMA GENE2FUNC及clusterProfiler进行功能注释。
结果：在MDD vs. 对照、MDD+ELS vs. 对照+ELS、以及MDD非ELS vs. 对照非ELS三组比较中，分析分别揭示了551、581和479个差异甲基化区域。核心基因SELENOO、FASN和THOC7在所有组间均呈现一致的差异甲基化，反映了共享的线粒体、代谢和RNA加工机制。MDD+ELS组表现出更高的DMR负荷，尤其在活性启动子染色质相关区域；ELS特异性基因HIVEP3和NRL展现出与难治性抑郁和注意缺陷多动障碍（attention deficit hyperactivity disorder, ADHD）相关的全基因组关联研究（genome-wide association study, GWAS）证据。相反，MDD非ELS组的DMR富集于抑制性染色质状态，其中CAMK4和CALM1位点的超甲基化提示钙依赖性神经元信号通路的组成性沉默。上述发现强调了ELS对MDD相关特异性表观遗传改变的显著贡献，特别是在神经发育和应激反应通路中。相比之下，不伴ELS的MDD表现出影响突触及突触可塑性的独特甲基化谱。
结论：因此，该研究强调了在精神科研究中采用创伤知情（trauma-informed）方法的必要性。
**研究背景与问题**
早期生活压力（early life stress, ELS）与重性抑郁障碍（major depressive disorder, MDD）之间的关系已成为当代精神病学研究的重要焦点。MDD影响全球超过20%的人口，给现代社会带来沉重负担。尽管对MDD病理和治疗的理解已有显著进展，但其病因学仍然复杂。ELS涵盖虐待和忽视等事件，是日后发生MDD及其他精神障碍的公认风险因素。暴露于ELS的个体常表现出独特的临床和生物学特征，因此探究ELS增加情感障碍易感性的精确机制成为关键研究领域。
表观遗传修饰日益被认为是环境暴露（如ELS）与个体遗传易感性相互作用从而塑造抑郁表型的关键机制。表观遗传指不改变DNA序列而改变基因功能或其他基因组过程的分子修饰，包括DNA甲基化、组蛋白修饰和非编码RNA调控。这些改变可影响神经内分泌反应、神经可塑性和神经传递，均与MDD病理生理相关。既往的全表观基因组关联研究（epigenome-wide association studies, EWAS）已鉴定出若干候选CpG位点和通路，如BDNF、KITLG、NR3C1和FKBP5等基因的甲基化改变与神经发育和应激反应相关，但存在重复性不一致的问题。然而，当前文献的一个显著局限是包含创伤分层对照组的研究相对匮乏。许多研究未能区分有无ELS的个体，从而将创伤暴露效应与MDD诊断本身的影响相混淆。同时纳入ELS阳性和ELS阴性对照，对于厘清早期逆境特有的表观遗传效应与单纯MDD相关的效应至关重要。
为弥补这一空白，研究人员开展了全基因组DNA甲基化研究，对按ELS暴露分层的MDD患者及创伤匹配的对照组的全血样本进行分析，旨在鉴定与ELS相关抑郁特异性关联的差异甲基化区域（differentially methylated regions, DMRs），并通过功能注释和通路富集阐明其生物学意义及潜在干预靶点。
**研究设计与结果概述**
研究发现，ELS与MDD患者中独特的表观遗传改变相关，提示早期逆境可能在更广泛的抑郁景观中塑造独特的分子轨迹。与无ELS者相比，同时患有MDD和ELS的个体表现出更广泛且潜在功能相关的DNA甲基化改变，尤其在神经发育、应激调节和细胞信号传导相关基因中。相比之下，无ELS史的MDD个体表现出更有限且系统性分布的甲基化改变。这些结果突显了早期应激暴露持久的生物学影响。
**主要关键技术方法**
研究人员从美国阿拉巴马大学伯明翰分校（University of Alabama at Birmingham, UAB）卫生系统招募住院及门诊MDD患者，并从伯明翰社区招募志愿者对照，共90名受试者，经DSM-5标准及MINI国际神经精神访谈确诊，使用童年创伤问卷（Childhood Trauma Questionnaire, CTQ）评估ELS。对全血样本采用简化代表性亚硫酸氢盐测序（Reduced-Representative Bisulfite Sequencing, RRBS）进行全基因组DNA甲基化谱分析。使用Bismark比对至人类参考基因组（Hg19），MethylKit和DSS软件识别差异甲基化CpG位点（differentially methylated CpGs, DMCs），DMRcate鉴定DMRs（≥5个CpGs, FDR<0.05, 5%甲基化差异）。采用ANNOVAR进行组学注释，FUMA GENE2FUNC和clusterProfiler进行功能富集分析，UCSC基因组浏览器进行局部可视化。
**研究结果**
**参与者特征**：四组分别为MDD伴ELS（n=30）、MDD不伴ELS（n=30）、对照伴ELS（n=15）、对照不伴ELS（n=15）。MDD组与对照组在汉密尔顿抑郁量表（HAMD-21）、哥伦比亚自杀严重程度评定量表（C-SSRS）、贝克抑郁量表第二版及CTQ各分量表上存在显著差异（P<0.05），组间种族和性别分布无显著差异。
**三组比较的整体甲基化差异**：全基因组分析在MDD vs. 对照、MDD+ELS vs. 对照+ELS、MDD非ELS vs. 对照非ELS三组比较中分别鉴定出29,094、34,091和37,218个DMCs；对应DMRs数量分别为551、581和479个。MDD vs. 对照比较中显著差异甲基化位点映射至NRXN2、CHTF8、DLEU1等基因；MDD+ELS组突出信号包括CNNM3、SPPL3、CHTF8、PBX1等；MDD非ELS组关键位点包括LINC02352、SOBP、KIF26B等。DMRs主要位于内含子、外显子和上游区域。染色质状态注释显示，MDD vs. 对照和MDD+ELS vs. 对照+ELS的DMRs主要富集于活性启动子（分别为169和194个），而MDD非ELS组的DMRs在异染色质和抑制区域富集更多。
**协变量调整敏感性分析**：在校正年龄、性别、种族、吸烟状态、BMI及估计的血细胞组成后，MDD vs. 对照中143个DMRs、MDD+ELS组126/127个（99.2%）DMRs、MDD非ELS组46/49个（93.9%）DMRs保持与原始分析一致的作用方向。
**基于CNS性状关联和组学注释的DMR优先化**：通过复合评分（最高4分）和甲基化差异对DMRs进行优先化排序，分别从三组比较中筛选出17、29和12个最显著DMRs。MDD vs. 对照中评分最高的基因包括SELENOO、NRL、SFN、ARID2、IRF2BPL、STK11等；MDD+ELS组包括HIVEP3、SFN、NRL、MEX3B、CHSY1、STK11、USF2、LSR等；MDD非ELS组包括CAMK4、IRF2BPL、SELENOO、TMEM201、NRDE2和CALM1。
**ELS与非ELS比较组的DNA甲基化差异**：三组比较的DMR关联基因Venn图显示，SELENOO、FASN和THOC7在所有组间一致改变，代表MDD不分ELS状态的共同表观遗传特征。MDD+ELS特异性基因包括HIVEP3、MEX3B、CHSY1、USF2、LSR、CRTAP、SPPL3、GATAD2A、PIF1等；MDD非ELS特异性基因包括CAMK4、TMEM201、NRDE2、CALM1、ZNF783、PAK1、SAC3D1、KIF26B等。Circos图显示ELS组比较中检测到更强的甲基化差异信号，且HIVEP3、SFN、NRL等基因的DMRs与难治性抑郁、一般认知能力、注意缺陷多动障碍、阿尔茨海默病精神病性症状及MDD等CNS性状的GWAS显著关联。
**差异甲基化CNS疾病基因的组学景观**：ELS相关MDD的LONRF1基因DMR区域呈现高水平H3K27Ac修饰、CpG岛、增强子和启动子元件富集，且存在MIR3926与LONRF1的潜在增强子调控元件互作。SFN和NRL基因的DMR区同样具有H3K27Ac标记、CpG岛和ATAC-seq信号。非ELS组DMR基因IRF2BPL、CAMK4和TMEM201也显示H3K27Ac、CpG岛、增强子及启动子元件富集。共享基因SELENOO、FASN和THOC7的DMRs均定位于CpG岛富集区域，并重叠活性染色质标记、注释增强子/启动子及转录因子结合位点。
**基因集富集分析**：FUMA GENE2FUNC分析显示，MDD vs. 对照的DMR基因在脑组织中高表达，包括皮层、额叶皮层（BA9）、前扣带皮层（BA24）、海马、下丘脑、杏仁核和基底节。KEGG通路富集于GnRH信号、脂联素信号、Wnt信号、谷氨酸能突触、Hippo信号和神经活性配体-受体相互作用。MDD+ELS组的DMR基因表达富集于额叶皮层、前扣带皮层、海马、杏仁核和基底节；KEGG通路富集于Hippo信号、AMPK信号、MAPK信号、Hedgehog信号、Rap1信号和凋亡通路；GO分析富集于细胞命运决定、肌肉分化、激酶活性调控等。MDD非ELS组的DMR基因中CALM1、IRF2BPL和FASN在海马、前扣带皮层和额叶皮层高表达；KEGG通路富集于Hedgehog、AMPK、Hippo、mTOR和GnRH信号；GO分析富集于脊髓细胞分化、上皮管形态发生、中枢神经系统神经元分化及谷氨酸能突触、突触后致密膜、树突棘等突触相关结构。
**讨论与结论**
本研究的系统评估表明，ELS暴露显著影响MDD患者的全基因组DNA甲基化谱，在ELS相关MDD与非ELS相关MDD之间鉴定出截然不同的甲基化模式，凸显了MDD通过复杂基因-环境相互作用形成的多面性病因。SELENOO、FASN和THOC7作为跨所有比较组一致超甲基化的核心基因，分别代表线粒体氧化还原调控、脂肪酸代谢和mRNA输出机制，提示这些通路可能是MDD不分ELS状态的基础表观遗传特征。MDD+ELS组表现出最高的DMR负荷和以活性启动子富集为特征的独特甲基化谱，关键基因共定位H3K27Ac标记、CpG岛及ATAC-seq信号，且富集于Hippo、AMPK、MAPK、Hedgehog和Rap1等信号通路，其优先化DMR基因与多种CNS性状的GWAS信号收敛，支持ELS介导的表观遗传重塑赋予跨神经精神结局的多基因风险负荷假说。相反，MDD非ELS组的DMRs富集于异染色质和抑制状态，CAMK4和CALM1的超甲基化提示钙依赖性神经元信号通路的组成性沉默，mTOR等通路的富集则关联突触可塑性和抗抑郁反应机制。
**研究结论**：该研究在伴与不伴ELS暴露的MDD个体中鉴定出独特的外周血DNA甲基化特征，提示ELS对MDD表观遗传变异具有显著贡献。然而，这些关联尚不能确立ELS、DNA甲基化改变与疾病发生之间的因果关系，需要纵向研究阐明其时间关系；且全血甲基化谱可能无法直接反映脑区特异性改变，样本量有限、缺乏匹配的基因表达数据也限制了对转录后果的直接评估。未来需更大规模独立队列、纵向采样、甲基化与转录组整合分析及外周和脑组织的功能验证，以确立这些发现的生物学意义并评估其作为MDD和ELS暴露生物标志物的潜在价值。综上，该研究强调了在精神疾病分子基础研究中纳入创伤史的重要性，创伤知情方法对理解MDD异质性和指导更精准的个体化干预具有关键意义。
