重性抑郁障碍静息态潜在fMRI动力学中短滞自回归延续性降低：多分析配置下的聚类稳健证据

《Psychiatry Research: Neuroimaging》：Reduced short-lag autoregressive carry-over in resting-state latent dynamics in major depressive disorder

　　

重性抑郁障碍（major depressive disorder, MDD）与静息态潜在功能磁共振成像（functional magnetic resonance imaging, fMRI）动力学中短滞时间延续性降低相关。在调整去噪流程（pipeline）、潜在维度K、左右利手（LQ）校正、年龄、性别和平均帧位移（mean framewise displacement, meanFD）后，MDD组的滞后₁自相关（lag-1 autocorrelation）和AR(1)增益φ较低。该效应集中于滞后₁，而非反映滞后1–10自相关函数（autocorrelation function, ACF）谱的均匀降低。短滞效应在多数去噪流程中保留，但在独立成分分析-自动去除运动伪影（ICA-AROMA，AROMA）下减弱。机制比较支持AR增益降低作为候选但非唯一的计算解释。

重性抑郁障碍（major depressive disorder, MDD）与静息态脑组织改变相关，涉及默认模式、额顶控制、显著性及背侧注意等网络的功能协调异常。既有研究多描述空间连接强度、网络耦合或静态差异，较少直接检验自发脑状态的时间延续性。动态功能连接研究提示MDD存在时间变异性、状态转换和时间变化交互，但通常不直接检验潜在脑状态从前一时点到下一时点是否短滞连续性降低。自相关和自回归增益可量化当前状态对后续状态的约束；内在时间尺度被视为神经组织原则。MDD临床特征包括持续情感状态、反刍、灵活性受损和内部思维模式转换改变，提示自发神经活动时间结构可能相关，但方向并非简单。已有研究开始关注抑郁的内在时间尺度和时间梯度，但MDD是否与静息态潜在动力学即时自回归延续性降低仍不清楚。因此，Nobuchika Yamaki与Tenna Churiki基于OpenNeuro ds002748静息态fMRI数据开展二次分析，纳入了51例MDD与20名健康对照，构建控制组定义的紧凑潜在动力学表征，检验MDD是否改变潜在脑状态轨迹的时间一致性。研究人员假设MDD短滞时间延续性降低，具体表现为滞后₁自相关和AR(1)增益φ较低，并在30种分析配置中评估三个主要指标，采用聚类稳健模型调整流程、K、LQ、年龄、性别和meanFD，并进行候选机制比较。结论是MDD在多数配置中显示较低滞后₁自相关和AR(1)增益φ，但三项主要指标经错误发现率（false discovery rate, FDR）校正后q≈0.0562；滞后1–10平均ACF未显著降低，组×滞后交互不显著。最佳表征为短滞自回归延续性较弱，最清晰证据在滞后₁，而非时间持久性均匀降低。该研究提供了受限的候选动力学特征，需更大、独立和纵向数据确认。

关键技术方法概括（不超过250字）：研究使用OpenNeuro ds002748公开静息态fMRI数据集（51例MDD、20名对照；TR=2.5 s，100 volumes；fMRIPrep预处理MNI空间）。采用5种区域时间序列去噪流程（minimal、motion24、acompcor5、acompcor5_gsr、aroma），设潜在维度K∈{2,3,4}与LQ校正开/关，共30种配置。将BOLD体积分为4×4×4网格，定义镜像左右bin，保留至少30个有效体素的配对，按控制组平均方差取前24对，计算双侧差异信号。每个配置仅用控制组拟合主成分分析（principal component analysis, PCA），再投射至MDD。计算滞后₁自相关、AR(1)增益φ、滞后1–10平均ACF及次要指标。统计采用重复配置线性模型、受试者聚类标准误和Benjamini–Hochberg FDR校正；机制比较基于控制组拟合的对角AR(1)模型进行候选变换仿真。

研究结果如下。3.1 样本与分析配置（Sample and analysis configurations）：所有71名受试者纳入，人口学与运动特征无显著组间差异；30种配置覆盖5种流程、3个K值和2种LQ设置；潜在空间解释双侧差异方差中位数为0.174，故结果描述紧凑动力学投影而非全脑BOLD动力学。3.2 MDD中滞后₁时间一致性降低（Reduced lag-1 temporal consistency in MDD）：30种配置中24种MDD滞后₁自相关较低（中位d=?0.363，中位置换p=0.168）；整合聚类稳健模型β=?0.0321，p=0.0364，三项主要指标BH校正后q=0.0562；结论为方向一致且未校正整合检验达常规阈值，但校正后略高于0.05。3.3 MDD中AR(1)增益φ降低（Reduced AR(1) gain φ in MDD）：24/30配置较低（中位d=?0.364，中位置换p=0.173）；整合模型β=?0.0316，p=0.0374，q=0.0562；参与者内平均φ对照0.0510、MDD 0.0196（Δ=?0.0314，d=?0.510，Welch p=0.0454）；结论为与滞后₁自相关一致，提示即时自回归延续性较弱，但主要家族FDR校正后略高于0.05。3.4 自相关差异的滞后-wise profile（Lag-wise profile of the autocorrelation difference）：滞后1–10平均ACF第三主要指标无清晰组降低（14/30方向符合，中位d=0.019，中位置换p=0.704；β=?0.0099，p=0.2192，q=0.2192）；联合重复测量模型组×滞后交互F(9,70)=1.340，p=0.233；滞后₁与后续滞后直接对比BH后最小q=0.097；结论为滞后₁证据最清晰，但未统计确立滞后₁特异性。3.5 流程敏感性（Pipeline sensitivity）：滞后₁ ACF与AR(1)增益在minimal、motion24、acompcor5、acompcor5_gsr下中位效应为负，最大在minimal与motion24；AROMA下略呈正值（d=+0.065）；结论为短滞组差异非预处理不变，但在多数流程存在。3.6 次要指标（Secondary metrics）：迹协方差在30/30配置MDD较低（β=?0.4379，p=0.0047，q=0.0189）；低频分数26/30较低（β=?0.0171，p=0.0382，q=0.0765）；平均潜在速度方向较低（β=?0.0600，p=0.0915，q=0.1220），内在时间尺度无清晰组差（β=?0.2979，p=0.2750，q=0.2750）；结论为不支持所有时间或动力学属性均匀改变，仅迹协方差通过次要家族FDR。3.7 机制比较（Mechanism comparison）：18种配置中，增益修改模型在8/18优于最佳非增益替代；AR增益降低最佳6种，噪声增加5种，振幅缩放4种，组合及基线各1种；结论为AR增益降低是候选计算解释之一，但非唯一。3.8 总结（Summary）：MDD在24/30配置中滞后₁自相关和AR(1)增益较低；未校正整合模型达显著，三项主要指标校正后略高于q=0.05；更广ACF无显著总体降低，组×滞后未确立滞后₁特异；最清晰次要结果为迹协方差降低；机制比较将AR增益降低列为竞争性候选而非唯一生成机制。

讨论部分总结：主要发现是MDD短滞时间延续性较低，由滞后₁自相关和平均AR(1)增益φ两个收敛指标表达，不应视作独立现象。效应在多数配置方向一致并在整合模型达常规显著，但三项主要指标FDR校正后未低于0.05；未均匀扩展至更广ACF，组×滞后未显著，故不能称“更慢”或全局持久性降低，应表述为控制组定义紧凑潜在表征中短滞自回归延续性较弱。该结果将焦点从空间连接转向潜在脑状态局部时间连续性；与内在时间尺度框架一致，但不表明全脑内在时间尺度普遍降低。临床持续情绪和反刍与神经短滞延续性属不同尺度，不能简化成“过度持续”或“过度不稳定”。滞后1–10平均ACF无稳健组差，说明不支持时间持久性均匀降低，滞后依赖未解。次要指标只有迹协方差通过校正，低频分数未通过，速度和内在时间尺度无清晰组差，不支持所有持久性、方差或速度指标一致偏移。五流程分析是优势：四个流程下呈负效应，AROMA下缺失并略呈正值，提示对预处理敏感，但meanFD调整和多数流程存在反对单纯运动解释。机制比较混合：AR增益降低最频繁但未主导，噪声与振幅接近，仅少于半数配置优于最佳非增益，故只能视为候选计算解释。方法学上，控制组定义PCA避免组混合轴，但解释方差中位0.174，仅为紧凑投影；双侧左右差异信号非图谱网络；30种配置为重复分析视图非独立重复，采用受试者聚类标准误；证据应看方向、效应量和敏感性而非二元显著。局限包括样本量较小尤其对照20名，横断面不能因果或状态标记，未建模药物、病程、症状维度、睡眠、焦虑等，潜在表征解释方差有限且解剖解释受限，AROMA不支持，主要指标FDR略高且滞后₁特异性未确立，机制比较仅18种配置。未来需独立MDD数据、其他精神疾病特异性、纵向症状与治疗变化、图谱与网络定位、更丰富生成模型及临床意义检验。

研究结论翻译：MDD与多数分析配置中较低的滞后₁自相关和较低AR(1)增益φ相关，两个效应在整合模型中达到常规显著性，但在三项主要指标校正后略高于q=0.05。滞后1–10平均自相关未显著降低，直接组×滞后检验未确立滞后₁特异性。因此，经验模式最好描述为较弱的短滞自回归延续性，最清晰证据在滞后₁，而非时间持久性的均匀降低。AR增益降低是一种合理计算变换，但未唯一优于噪声或振幅替代。这些发现提供MDD静息态活动的受约束候选动力学表征，需在更大、独立和纵向数据集中确认。