《Neurotherapeutics》：Electroconvulsive therapy shifting brain dynamics toward more adaptive and integrated states in treatment-resistant depression based on a hidden markov model
编辑推荐：
中文翻译：难治性抑郁症（treatment-resistant depression, TRD）以症状持续存在且对常规治疗反应不佳为特征，但其动态神经机制及治疗相关变化仍未完全阐明。本研究应用隐马尔可夫模型（hidden Markov model, HMM）结
中文翻译：难治性抑郁症（treatment-resistant depression, TRD）以症状持续存在且对常规治疗反应不佳为特征，但其动态神经机制及治疗相关变化仍未完全阐明。本研究应用隐马尔可夫模型（hidden Markov model, HMM）结合神经递质图谱映射，探究TRD的全脑时空动态及其经电休克疗法（electroconvulsive therapy, ECT）的调节作用。研究人员收集了72例TRD患者（52例完成随访）和63例健康对照的静息态功能磁共振成像（resting-state fMRI）数据。HMM识别出六个由不同功能连接和平均功能活动模式定义的反复出现的脑状态。基线时，TRD患者表现出全局协调状态的占用率和稳定性降低，而皮层下主导状态的参与增加，提示全局稳定性受损及皮层–皮层下平衡紊乱。ECT后，临床显著改善，伴随默认网络（default mode network, DMN）和额顶网络（frontoparietal network, FPN）内活动较高的状态以及高度整合的全脑状态占用率增加。这些变化提示整体脑协调和皮层参与的改善，而非基线异常的完全正常化。相比之下，皮层下主导状态在治疗后持续存在，与TRD的类特质特征一致。神经递质映射显示，治疗响应状态与血清素能、多巴胺能和谷氨酸能系统相关，而持续存在的皮层下状态与去甲肾上腺素能和胆碱能标志物相关。总体而言，这些发现提示ECT通过将脑动态转向更具适应性和整合性的状态来改善症状，而非完全恢复基线功能。这一动态视角为TRD及其治疗的神经机制提供了新见解。
论文解读
一、研究背景与问题
难治性抑郁症（treatment-resistant depression, TRD）指对至少两种足量足疗程抗抑郁药治疗反应不佳，影响高达55%的重性抑郁障碍（major depressive disorder, MDD）患者，且缓解率低于15%。TRD不仅伴随显著的生活质量下降，还与高自杀风险相关，造成沉重的社会和经济负担。电休克疗法（electroconvulsive therapy, ECT）是TRD实现快速缓解最有效的治疗手段，但其神经机制尚未完全阐明。既往神经影像研究多依赖静态功能连接（functional connectivity, FC）分析，发现ECT可调节边缘系统和边缘–前额叶网络，并恢复默认网络（default mode network, DMN）内的连接。然而，静态FC假设脑活动随时间保持稳定，可能忽略与抑郁症状相关的短暂且反复出现的脑状态。脑是持续变化的动态系统，越来越多的证据提示ECT通过改变这些随时间变化的脑活动模式发挥作用。既往动态研究多采用滑动窗口方法，依赖任意的窗口长度，难以捕捉离散脑状态之间的快速转换。隐马尔可夫模型（hidden Markov model, HMM）无需时间平滑即可识别不同脑状态，并计算每个时间点状态激活的精确概率，为精神疾病中的微状态转换提供更敏感的视角，但将其应用于纵向治疗研究仍有限。此外，理解这些动态变化背后的神经化学基础也是重要挑战。因此，本研究将HMM与神经递质图谱映射相结合，旨在从动态角度阐明TRD的神经机制及ECT的治疗作用。
二、研究内容与结论
研究人员对74例TRD患者和66例健康对照（healthy controls, HCs）进行静息态功能磁共振成像（resting-state fMRI）采集，排除头动过大者后，最终纳入72例TRD基线患者、63例HCs和52例完成ECT随访的患者。ECT采用spECTrum 5000Q系统，标准疗程6–12次。HMM分析识别出六个脑状态，并计算分数占用率（fractional occupancy, FO）、平均寿命（averaged lifetime, ALT）和转换率（switching rate, SR）三个动态指标。研究发现，基线时TRD患者在全局协调状态（状态1）的FO和ALT降低，在皮层下主导状态（状态5）的FO升高；ECT后，患者在DMN和FPN活动较高的状态（状态2）及高度整合状态（状态3）的FO增加，但状态1和状态5的异常持续存在。神经递质映射显示，状态2的活动与5-HT2a、D1、GABAa、M1和NMDA受体分布正相关，状态3与5-HT2a、GABAa和M1负相关，状态5与NET和VAChT正相关。研究结论认为，ECT通过将脑动态转向更具适应性和整合性的状态来改善症状，而非完全恢复基线功能，为TRD的治疗机制提供了动态视角。
三、主要技术方法
本研究采用高斯HMM，使用hmmlearn库在Python中实现，通过期望最大化算法估计参数，状态数从5到12进行测试，最终基于模型拟合和可解释性选择k=6。每个状态由平均向量（代表平均功能活动，AFA）和协方差矩阵（代表区域间FC）定义。动态指标包括FO、ALT和SR。神经递质关联分析使用JuSpace工具箱，将各状态AFA图谱与来自1200余名健康个体的19种神经递质受体/转运体正电子发射断层扫描（positron emission tomography, PET）图谱进行Spearman相关分析，并采用错误发现率（false discovery rate, FDR）校正。样本来源于内蒙古自治区精神卫生中心，TRD患者于2018年6月至2023年12月间招募。影像数据在3.0T Skyra扫描仪上采集，使用Schaefer-100皮层图谱和Tian 2020皮层下图谱提取116个感兴趣区的时序信号。统计采用线性混合效应模型（linear mixed-effects model, LMM），以条件为固定效应、受试者ID为随机效应，年龄、性别和教育年限为协变量，事后比较采用Wald检验和Bonferroni校正。
四、研究结果
（一）人口学和临床特征
最终样本包括72例TRD基线患者、63例HCs和52例完成随访的TRD患者。各组在年龄、性别和教育程度上无显著差异。随访完成者与未完成者、完整TRD队列与随访子集之间基线特征无差异。ECT后，患者HAMD、HAMA、SDS、SAS、ATQ和BPRS评分均显著降低（均p<0.001），表明临床显著改善。
（二）时空状态动态和转换模式
HMM识别出六个脑状态。状态1在所有网络内和网络间FC较高（仅次于状态3），DAN和皮层下区域AFA降低；状态2 FC中等，DMN和FPN的AFA升高；状态3 FC最高，DMN、FPN和视觉网络（visual network, VN）AFA降低；状态4 DAN内FC中等，DMN和边缘网络（limbic network, LN）AFA升高，SMN AFA降低；状态5皮层下FC较强，LN和VN AFA降低；状态6各网络FC普遍较弱，DAN、SMN和VN AFA升高。状态转换呈非随机、非顺序性，HCs、TRD基线和TRD治疗后三组转换模式不同。
（三）时间动态的组间差异
FO方面，HCs在状态1高于TRD基线和治疗后（p=0.010和p<0.001）；治疗后TRD在状态2高于HCs（p=0.004）；HCs在状态3高于两个TRD组（p<0.001和p=0.003），但ECT后较基线显著增加（p=0.003）；TRD在状态5两个时间点均高于HCs（p=0.016和p<0.028）。ALT方面，治疗后TRD在状态1的持续时间短于HCs和基线（p=0.009）。SR无显著组间差异。
（四）与临床指标的相关性
状态6的FO与SAS评分负相关（r=?0.283, p=0.037）。状态4 FO的减少和状态6 ALT的减少分别与SDS和SAS的降低率相关（r=?0.308, p=0.031；r=?0.328, p=0.021）。基线预测分析显示，HAMA降低率与状态2基线ALT正相关（r=0.414, p=0.017），与状态6基线ALT（r=?0.340, p=0.034）和状态6基线FO（r=?0.321, p=0.024）负相关，但均为未校正的名义相关。
（五）与神经递质的空间相关性
状态2的AFA与5-HT2a（r=0.363, p<0.001）、D1（r=0.376, p<0.001）、GABAa（r=0.360, p<0.001）、M1（r=0.328, p<0.001）和NMDA（r=0.321, p<0.001）受体密度正相关。状态3的AFA与5-HT2a（r=?0.603, p<0.001）、GABAa（r=?0.618, p<0.001）和M1（r=?0.519, p<0.001）负相关。状态5的AFA与NET（r=0.620, p=0.002）和VAChT（r=0.512, p=0.002）分布正相关。
五、讨论与结论
讨论部分指出，TRD在特质层面表现为全局稳定性和皮层–皮层下平衡的缺陷。状态1的参与和稳定性降低，反映DAN功能障碍和皮层下情绪系统调节能力减弱；状态5的持续参与提示皮层对核心情绪系统的控制减弱。这些异常在ECT后持续存在，属于类特质特征，可能增加TRD复发风险。临床改善主要与治疗相关的皮层参与和全局协调变化有关。ECT后状态2占用率增加，反映DMN和FPN相关皮层活动恢复和自上而下调控增强；状态3占用率增加提示全脑高效信息处理和协调能力改善。神经递质映射进一步阐明，状态2与单胺能和谷氨酸能系统相关，状态3与抑制性系统失衡相关，状态5与去甲肾上腺素能和胆碱能系统失调相关。研究提出ECT诱导脑动态的功能性“重置”模型，即通过增加适应性脑状态的可及性来补偿特质性缺陷，而非完全纠正潜在功能障碍。局限性包括患者持续服药无法完全排除药物影响、样本流失、使用群体水平受体图谱而非个体化成像、对照组仅基线扫描缺乏重复测量，以及HMM状态数选择的方法学挑战。
结论部分翻译：本研究基于HMM为TRD及其ECT治疗提供了动态视角。TRD以基线稳定性和皮层–皮层下调节的持续性缺陷为特征，而临床改善与皮层参与增加和全脑协调改善相关。这些发现强调将脑功能视为动态过程的重要性，并提示有效治疗可能通过重塑而非完全正常化脑活动模式来发挥作用。
