帕金森病运动症状的共享和症状特异性局部场电位频谱特征

《Neurobiology of Disease》：Shared and symptom-specific local field potential spectral signatures of motor symptoms in Parkinson's disease

　　

自适应深部脑刺激（adaptive deep brain stimulation, aDBS）用于帕金森病（Parkinson's disease, PD）旨在通过基于反映患者临床状态的神经生理反馈信号调整刺激振幅来提高治疗效果。尽管来自局部场电位（local field potential, LFP）记录的β频带（± 13–35?Hz）频谱功率是迄今为止特征最明确的生理标志物，但其与运动症状严重程度的相关性仍然有限。为了探究结合多个频率的频谱信息是否能提高运动症状严重程度的解释方差，研究人员将典型相关分析（canonical correlation analysis, CCA）应用于来自67名接受丘脑底核DBS的PD患者的LFP记录的全功率谱密度（1–100?Hz），并结合运动障碍学会统一帕金森病评定量表第三部分（Movement Disorders Society Unified Parkinson's Disease Rating Scale part III, UPDRS-III）的聚类项目和单个项目。CCA冗余指数量化表明，UPDRS-III评分中超过20%的方差可由跨频率频谱功率的线性组合解释，而仅β频带功率解释约10%。β频率对总UPDRS-III评分贡献最强且为正相关，其中低β（约13–25?Hz）范围的贡献略高于高β（约26–35?Hz）范围。低频（<10?Hz）和刺激频率一半处的精细调谐伽马活动（62–63?Hz）贡献为负。症状特异性分析显示，与运动迟缓相关的清晰双β峰，而僵直主要与低β功率相关。震颤表现出较不明显的频谱模式。CCA结果在对侧和同侧症状中高度相似。总体而言，CCA结果揭示了PD运动症状背后既有共享也有症状特异性的频谱模式，并证明了更广泛频谱生理标志物的附加价值。

**论文解读：帕金森病运动症状的共享与症状特异性局部场电位频谱特征**

**研究背景与问题**

帕金森病（Parkinson's disease, PD）是一种进行性神经退行性疾病，以运动迟缓、震颤和/或僵直为主要运动症状。深部脑刺激（deep brain stimulation, DBS）是治疗中晚期PD的有效手段，但传统DBS采用固定参数刺激，无法根据患者症状波动进行实时调整。自适应DBS（adaptive DBS, aDBS）通过神经电生理反馈信号（如局部场电位（local field potential, LFP）中的β频带功率）调整刺激幅度，有望改善疗效。然而，β频带（约13–35?Hz）功率与运动症状严重程度的相关性仅能解释约16–17%的个体间变异，尤其对震颤相关性弱。此外，不同频率频带（如θ频带、γ频带）与不同症状存在关联，提示单一频带难以完整反映运动症状的复杂性。因此，需要探索全频谱信息能否提高症状严重程度的解释力，并识别各运动症状的特征性频谱模式。

**研究内容与结论**

研究人员对来自67名接受丘脑底核（subthalamic nucleus, STN）DBS的PD患者（共131个半球）的839次LFP静息态记录和同步的UPDRS-III评分进行分析，采用典型相关分析（canonical correlation analysis, CCA）将全功率谱密度（1–100?Hz）与UPDRS-III总分、症状子分数（运动迟缓、僵直、震颤）、对侧/同侧半身分数及所有单项分数进行关联。研究发现，CCA冗余指数（RI）显示全频谱线性组合能解释UPDRS-III总分超过20%的方差，而仅β频带功率仅解释约10%。β频率贡献最强且为正，其中低β（~13–25?Hz）略强于高β（~26–35?Hz）；低频（<10?Hz）和精细调谐伽马（62–63?Hz，刺激频率一半处）贡献为负。症状特异性分析显示，运动迟缓具有清晰的双β峰（低β和高β均强贡献），僵直主要与低β相关，震颤则无明确频谱模式。对侧和同侧症状的频谱特征高度相似。该研究论文发表在《Neurobiology of Disease》。

**关键技术与方法**

本研究为前瞻性观察性队列研究（Amsterdam University Medical Center），纳入2022年9月至2024年7月期间接受双侧STN-DBS手术的PD患者（植入Medtronic Percept? PC神经刺激器和SenSight? B33005电极）。主要技术方法包括：（1）采用BrainSense? Streaming技术以250?Hz采样频率记录静息态STN-LFP，并使用奇异值分解自动去除心电图伪迹，再选取最长连续无伪迹段（≥10秒）进行Welch法功率谱密度估计（1–100?Hz），相对化处理（1–49 Hz和51–99 Hz总功率归一化）；（2）临床评估采用UPDRS-III评分，通过视频回顾确定（僵直项由认证医师现场评估），包括总分、症状子分数（运动迟缓、僵直、震颤）、对侧/同侧半身分数及18个聚类；（3）统计分析：先计算β频带（13–35 Hz）功率与各UPDRS-III评分的Pearson相关（及Bonferroni校正），再应用CCA对全频谱密度与UPDRS-III评分进行多变量分析，计算典型负载、典型相关系数、冗余指数等。

**研究结果**

**3.1 人口学与临床特征**：最终分析纳入67名参与者（24女/43男），年龄中位65岁，病程中位8年，共839次观察（131半球）。**3.2 β频带功率与UPDRS-III评分的相关性**：Pearson相关显示，β频带功率与运动迟缓+僵直子分数（r=0.3233, R2=10.45%）、运动迟缓子分数（r=0.3143, R2=9.88%）、UPDRS-III总分（r=0.3101, R2=9.62%）相关性最强；与僵直相关性较弱（r=0.2212, R2=4.89%）；与震颤无显著相关（r=0.0934）。对侧症状相关性略高于同侧。**3.2.1 CCA：全频谱与UPDRS-III总分**：典型函数显著（λ=0.77, p<0.001），典型相关系数r=0.4806，冗余指数RI=23.10%（即全频谱线性组合解释23.10%的总分方差），远高于β频带功率的9.62%。负载显示β频带正向贡献最强（低β略强于高β），而低频（<10 Hz）和62–63 Hz（精细调谐伽马）呈负向贡献。**3.2.2 CCA：全频谱与症状子分数**：运动迟缓（RI=24.14%）和僵直（RI=23.70%）均与β频带强相关，但运动迟缓涉及更宽β范围（14–34 Hz）且具有双β峰，僵直主要限于低β（16–23 Hz）；震颤（RI=22.94%）无明显频谱模式，负载均<0.3。**3.2.3 CCA：全频谱与对侧/同侧半身分数**：对侧和同侧半身分数的频谱模式高度相似（RI分别为21.87%和21.53%），联合分析仅第一典型函数显著（λ=0.691, p<0.001），负载体现β正向和低频、精细调谐伽马负向贡献。**3.2.4 CCA：全频谱与所有UPDRS-III单项**：前9个典型函数显著（p<0.05），累计解释UPDRS-III单项12.6%的方差。函数1显示运动迟缓（尤其对侧）与低频和低β相关；函数2显示腿部运动迟缓和步态与双β峰（低β和高β）相关；函数3显示僵直贡献但无特异性频谱；函数4显示震颤与低频相关但伴有僵直项目。结果提示症状类别存在部分聚类及特征性频谱。

**讨论与结论**

本研究通过CCA证明，全频谱（1–100?Hz）线性组合比单独β频带功率能解释更多PD运动症状方差（>20% vs ~10%），支持更宽频谱生理标志物的价值。β频带中低β和高β均有贡献，但运动迟缓存在双β峰，僵直主要与低β相关，提示低β与高β可能具有不同功能角色（低β更多反映病理状态，高β与生理运动活动相关）。低频负向贡献与过往研究一致（尤其在左旋多巴开启状态），而精细调谐伽马（刺激频率一半）的负向贡献可能反映运动促进状态与运动障碍的反向关系，尽管其负载未达强贡献阈值（可能因刺激幅度中位仅1.1 mA）。震颤缺乏明显频谱模式，可能因STN-LFP未完全捕获震颤相关的小脑-丘脑-皮层环路活动，且仅37%的记录存在震颤。CCA结果受限于线性假设、非正态数据、刺激和药物混杂效应及纵向稳定性未直接评估。结论：考虑全功率谱密度（1–100?Hz）信息能显著提高对PD运动症状变异的解释力，不同运动症状存在共享和特异性的频谱模式，这些发现有助于发展更个性化的aDBS治疗。研究结论部分翻译：综上所述，研究人员的发现表明，与单独的平均β频带功率相比，考虑全功率谱密度（1–100?Hz）的信息可以解释更多PD运动症状的变异。低频和低β频带功率与运动迟缓和僵直严重程度相关，而运动迟缓还额外涉及高β频带。震颤缺乏独特的频谱模式，但在震颤和僵直中均观察到刺激频率一半处的精细调谐伽马活动的负向贡献。总体而言，不同运动症状的频谱模式大部分重叠，提示共享的神经机制。然而，微妙但独特的症状特异性频谱模式可能有助于识别个体运动症状的神经特征，从而实现更个性化的PD aDBS治疗。