重度抑郁障碍（MDD）患者氯胺酮辅助心理治疗（KAP）后纵向外周DNA甲基化重塑及NOTCH1相关信号模块探索性研究

《Epigenetics》：Longitudinal peripheral DNA methylation remodeling following ketamine-assisted psychotherapy in major depressive disorder

　　

氯胺酮可产生快速抗抑郁效应，但其效应持续相关的分子过程在人类中仍界定不足。研究人员在10名接受氯胺酮辅助心理治疗（ketamine-assisted psychotherapy, KAP）的重度抑郁障碍（major depressive disorder, MDD）成人中，使用两次治疗前评估（T1、T2）和一次治疗后评估（T3），检查纵向DNA甲基化变化。参与者在一个结构化心理治疗方案内接受三次鼻内外消旋氯胺酮（intranasal racemic ketamine）治疗。使用Illumina MethylationEPIC v2.0 BeadChip测量颊部DNA甲基化。治疗后抑郁症状显著改善，贝克抑郁量表-II（Beck Depression Inventory–II, BDI-II）平均分从重度范围降至轻度范围（p < .001）。在评价的802,595个CpG位点中，20,770个满足名义T3-versus-T2筛选标准。对甲基化变化幅度以及相对于治疗前波动的移动进行序贯过滤后，识别出962个保留CpG位点，包括517个保留的基因注释位点。补充线性混合效应分析对603个预选基因注释且经幅度过滤的位点中参与者的重复测量进行建模；由于这些位点使用相同对比和队列选择，该分析未被视为独立验证。对这603个位点集的探索性富集分析识别出由NOTCH1、DLL1和ADAM17驱动的NOTCH1相关Reactome条目。在保留特征中，多个核心和邻近Notch相关组件形成连贯信号模块，涵盖配体、受体加工、转录调节和通路交叉对话。这些发现识别出与氯胺酮辅助心理治疗相关的结构化外周甲基化特征。它们不建立脑特异性通路激活或因果关系，但为研究伴随临床改善的调节适应提供纵向基础。

研究背景方面，氯胺酮（ketamine）是一种快速起效抗抑郁药，其临床疗效具有独特时间特征：单次亚麻醉剂量可在数小时内启动治疗效应，并持续超过药物代谢清除时间。其主要作用由N-甲基-D-天冬氨酸（N-methyl-D-aspartate, NMDA）受体拮抗驱动，触发谷氨酸爆发及α-氨基-3-羟基-5-甲基-4-异恶唑丙酸（alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid, AMPA）受体激活，继而涉及真核延伸因子2（eukaryotic elongation factor 2, eEF2）去磷酸化和脑源性神经营养因子（brain-derived neurotrophic factor, BDNF）快速翻译，导致短暂突触生成激增。然而，稳定这些新生连接的分子框架尚不明确。行为证据提示氯胺酮可重开约48小时的社会奖赏学习关键期。DNA甲基化（DNA methylation，指DNA上胞嘧啶碱基的甲基化修饰）等神经表观遗传机制可能将短暂药理触发转化为稳定基因表达改变。现有证据多来自动物模型或死后脑组织，活体人类纵向研究较少，尤其缺乏可区分治疗相关变化与基线波动的受试者内设计。外周组织提供可行窗口，但甲基化具组织特异性；既往研究提示应激生理、内分泌调节和免疫信号相关位点可在中枢与外周组织间呈现协调变化。因此，研究人员开展本研究，旨在刻画MDD患者接受KAP后的纵向DNA甲基化重塑，假设治疗不会广泛改变上游受体或神经递质基因，而会选择性重塑涉及神经可塑性、应激适应和长期细胞调节的下游调节系统。

关键技术方法概括如下：研究纳入10名经Mini International Neuropsychiatric Interview（MINI 7.0.2）诊断的MDD成人和12名健康对照，获University of Ottawa Research Ethics Board批准。MDD参与者接受标准KAP方案，约7周内完成三次鼻内外消旋氯胺酮（180–270 mg）治疗。颊部DNA在T1、T2（T1后2周）和T3（治疗后）采集；健康对照间隔2周提供两次样本，仅用于质量控制，不纳入T3-versus-T2治疗相关分析。全基因组甲基化用Illumina MethylationEPIC v2.0 BeadChip测定，经SeSAMe QCDPB流程处理。主分析以配对t检验比较T2与T3，名义p < .05筛选；随后以Δ_max ≥ .035进行幅度过滤，并以R = |T3 ? T2| / |T2 ? T1|评估稳定性：R ≥ 3为基线稳定治疗反应位点（Baseline Stable Treatment-Responsive loci, BSTR），2 ≤ R < 3为Dynamic Movers，R < 2排除。补充线性混合效应模型（linear mixed-effects model, LMM）对603个基因注释位点建模，g:Profiler用于探索性富集分析。

结果部分，基线稳定性：临床量表显示治疗前抑郁症状和人格域高度稳定，提示KAP前无自发临床漂移。抑郁症状：重复测量ANOVA显示BDI-II显著下降；因Mauchly检验提示球形性违反，采用Greenhouse–Geisser校正，F(1.17, 10.50) = 55.67，p < .001，效应量大，η² = .861；两次基线间无显著差异（p = .753），基线与治疗后均显著改善（均p < .001）；BDI-II均值从T1的36.7、T2的35.2降至T3的14.8，由重度范围转为轻度范围。全基因组DNA甲基化与表观基因组范围变化概览：802,595个CpG位点进入分析，20,770个满足T3-versus-T2名义p < .05；T3-versus-T2甲基化变化分布中心近零，部分位点尾部变化较大；T2-versus-T1描述治疗前波动，与T3-versus-T2比较显示治疗后部分位点移动更大。治疗相关甲基化信号识别：幅度过滤Δ_max ≥ .035后，20,770个位点降至1,180个；稳定性比值分类得BSTR 692个、Dynamic Movers 270个、R < 2排除218个，最终保留962个CpG位点。1,180个幅度过滤位点中603个有基因注释，包括359个BSTR、158个Dynamic Movers和86个排除位点，最终基因注释保留位点为517个。3.5百分点阈值用于区分生物学上可解释的治疗相关信号与技术或随机测量变异；该阈值低于强环境暴露相关有时观察到的>10%差异，同时覆盖阵列研究中常见约1–7%的干预和疾病相关甲基化信号；幅度过滤也限制了R在分子分母均很小时出现虚假增大的不稳定性。线性混合效应敏感性分析：对603个基因注释位点在10名MDD参与者的T1、T2、T3数据建模，时间点为固定效应，参与者为随机截距；所有603个位点T3-versus-T2对比名义p < .05；由于这些位点用同一对比和队列预选，该分析不作为独立显著性验证，仅用于检查在显式建模参与者内依赖后变化方向和幅度是否仍明显。功能富集分析：g:Profiler对603个基因注释幅度过滤位点映射的唯一基因符号进行分析，识别出两个显著且密切相关的Reactome条目，涉及NOTCH1 t(7;9) translocation mutant背景下的NOTCH1信号，校正p = .0249，均由NOTCH1、DLL1和ADAM17重叠驱动；还观察到GO生物过程术语multicellular organismal process和若干转录因子motif显著富集。因86个位点未满足最终稳定性标准，这些发现被解释为刻画更广泛的幅度过滤甲基化集，而非最终517个位点保留签名。Notch相关信号模块识别：保留的基因注释集包含多个经典Notch信号组件，跨越胞外信号、受体加工、胞内信号整合和转录调节。核心Notch组件包括NOTCH1受体、DLL1配体和转录效应因子RBPJ；邻近调节因子包括ADAM17、BAMBI和SKI。模块还包含DYRK2、TSC2、STK32B、ACTR1B、BCL2L1、HSF1、ERN1等胞内信号调节因子，以及核内RBPJ、SKI、NR5A2和CRHBP。DYRK2在模块内显示最大治疗相关甲基化移动，因治疗后变化仍超过基线波动两倍以上而被保留为Dynamic Mover。数据库注释支持NOTCH1相关信号机制在该集合中的表现，而非参与者存在translocation或致癌性NOTCH1激活的证据。Notch模块内纵向甲基化动态：模块内位点显示治疗相关甲基化变化超过基线波动；ACTR1B和SKI相关保留位点符合BSTR标准；DYRK2治疗相关变化最大但归类为Dynamic Mover。T3与T2比较的Δβ显示模块内甲基化变化呈双向，并非统一升高或降低；所示位点均满足幅度和稳定性保留标准。

讨论部分总结：研究采用受试者内双基线设计，观察到KAP后抑郁症状显著降低，治疗后甲基化变化超过两次治疗前评估的组水平波动，符合治疗相关生物重塑，也与氯胺酮短暂增强可塑性后由下游调节过程稳定的模型一致。中心发现是定量过滤后保留962个CpG位点，其中517个有基因注释并纳入生物学解释；这些位点分布于转录调节、信号整合、应激适应、结构重塑、免疫调节、突触组织和电生理功能等多个生物学域，提示KAP相关表观遗传变化具有选择性和组织性，而非弥漫性漂移。探索性富集分析在更广泛603位点集中识别NOTCH1相关信号条目，保留签名中包含多个相互作用的Notch组件，涵盖配体、受体加工、转录调节和通路交叉对话，因此Notch相关模块被选作假设生成性解释；其余保留位点未在本次分析中分配正式功能类别。DYRK2为双特异性激酶，实验研究提示其可通过磷酸化依赖性蛋白酶体降解调节NOTCH1胞内结构域稳定性，并影响Hedgehog信号、神经元细胞骨架组织和突起生长，但这些发现未确立DYRK2在氯胺酮反应或抑郁中的作用。应激和内分泌相关基因在保留签名中较突出，涉及HPA轴调节、代谢控制和细胞应激反应，与成功抗抑郁反应伴随应激相关生物系统协调正常化及神经可塑性改善的模型一致。外周甲基化不应解释为脑特异性表观遗传状态的直接代理，颊部与脑组织甲基化一致性高度依赖位点。局限性包括样本量小（n = 10）、颊部组织而非脑组织、甲基化变化不能假定反映基因表达、蛋白丰度或通路活性、仅有一个治疗后时间点，尚不清楚观察到的变化是否持续、正常化或继续演化。结论翻译：氯胺酮辅助心理治疗与显著临床改善及纵向DNA甲基化重塑的协调模式相关。应用定量过滤框架识别出962个保留CpG位点，包括517个保留的基因注释位点，涉及转录调节、信号整合、应激适应和神经可塑性相关过程。对更广泛的603个基因注释且经幅度过滤位点集的探索性富集分析识别出NOTCH1相关信号机制的表现。在保留特征内，一个由NOTCH1、DLL1、RBPJ、ADAM17、BAMBI、SKI及相关调节基因组成的Notch相关模块作为突出生物学特征出现。虽然仅凭甲基化数据不能得出因果推断，这些发现为未来研究氯胺酮辅助心理治疗后临床反应相关的调节系统提供了基础。