基于脑形态测量学与类淋巴系统参数的多模态神经影像机器学习模型预测帕金森病STN-DBS术后程控早期运动反应

《CNS Neuroscience & Therapeutics》：Multimodal Neuroimaging-Based Machine Learning Models Leveraging Cerebral Morphometry and Glymphatic Parameters Predict Short-Term Post-Programming STN-DBS Motor Response in Parkinson's Disease

　　

背景与目的：丘脑底核脑深部电刺激（STN-DBS）对药物难治性帕金森病（PD）运动症状有效，但临床反应在症状域间存在差异，尤以震颤与步态为著。准确的术前分层具临床重要性，尤其针对程控早期结局。方法：研究人员回顾性纳入155例接受双侧STN-DBS的PD患者与43名健康对照。术前结构磁共振成像与扩散加权成像用于量化脑形态测量及类淋巴标记，包括沿血管周围空间扩散张量成像（DTI-ALPS）与脉络丛体积（CPV）。所有155例患者评估总运动反应，133例具完整震颤分项者评估震颤反应，66例具配对仪器化步态评估者评估探索性数据驱动步态改善表型。采用折叠内特征选择与超参数调优建立机器学习模型，以五折交叉验证评价。结果：最优三模态模型总运动反应AUC为0.850±0.045（95% CI 0.794–0.906），震颤反应AUC为0.861±0.047（95% CI 0.803–0.919），探索性步态改善表型AUC为0.970±0.019（95% CI 0.946–0.994）。仅形态测量模型仍具相当预测性能，运动、震颤、步态终点最大AUC分别为0.830、0.849、0.955。探索性步态表型中，临床加形态与三模态模型表现相近，提示该类淋巴变量在此亚组增量价值有限。结论：术前脑形态测量辅以选定类淋巴标记与基线临床变量，可辅助分层PD患者STN-DBS程控早期反应。结果支持影像知情DBS结局预测的进一步发展，但临床实施前仍需外部验证与更长随访。

研究背景与问题提出

帕金森病（PD）为第二常见神经退行性疾病，以震颤、僵直、运动迟缓及轴性障碍进行性加重为特征。丘脑底核脑深部电刺激（STN-DBS）可改善药物难治性运动症状，但不同症状域（尤其震颤与步态）获益幅度与模式异质性强。既往DBS疗效预测算法多聚焦临床表型与连接指标，曲线下面积（AUC）多低于0.85。STN-DBS作用于皮质-基底节-丘脑-皮质分布环路，而PD进展伴区域选择性萎缩与网络紊乱可能影响反应性。类淋巴系统负责血管周围液清除，调控α-突触核蛋白稳态与神经炎症，与PD进展及DBS反应相关；动物研究提示STN-DBS可增强间隙液引流并恢复水通道蛋白4（AQP-4）极化。基于此，研究人员整合术前脑形态测量、类淋巴标记（DTI-ALPS、CPV）与基线临床变量，构建多模态机器学习框架，预测STN-DBS程控早期总运动、震颤及探索性步态改善表型反应，并直接比较临床单模、影像单模、双模与三模模型互补性。该研究发表于《CNS Neuroscience & Therapeutics》。

主要关键技术方法

研究为复旦大学附属华山医院单中心回顾性队列（2022年11月至2023年12月），纳入155例双侧STN-DBS的PD患者及43名健康对照；终点特异性队列因缺失分量而异（总运动155例、震颤133例、仪器化步态66例）。术前3T MRI采集高分辨T1-MPRAGE与32方向DWI（b=1000 s/mm²）。FreeSurfer v7.4.1重构皮质厚度、面积、平均曲率、相对皮质/皮下体积；FSL与MRtrix3管线计算DTI-ALPS（投影纤维Dx_proj、Dy_proj与联络纤维Dx_assoc、Dz_assoc比值）；ITK-SNAP手动5 mm³ ROI量化；CPV经颅内体积（ICV）标准化。临床终点：Med-OFF/ON MDS-UPDRS III、震颤分项、LCT改善率、DBS改善率（高反应者总运动≥50%、震颤≥60%）；步态亚组对10项仪器化步态参数变化做UMAP降维与bisecting K-means聚类（K=2）得数据驱动表型。特征选择用RFE+随机森林重要性（≤10特征/终点），六算法（XGBoost、RF、KNN、LDA、NB、SVM）在StratifiedGroupKFold嵌套五折中调优，SHAP解释。

研究结果

2.1 受试者基线：PD组155例（62.50±8.75岁，60%男）。高运动反应者（n=59）LCT总改善率更高、Med-OFF/ON轴性分项更低；高震颤反应者（n=72）僵直分项更高。

2.2 类淋巴与形态差异：PD vs 健康对照全脑及右半球DTI-ALPS降低（p=0.014/0.018），右投射纤维Dy_proj、右联络纤维Dz_assoc升高，左右绝对/相对CPV升高。高运动与高震颤反应者绝对及相对CPV显著低于低反应者；DTI-ALPS各组间无差。形态学在运动/震颤应答亚组间多个皮质厚度、面积、曲率、体积特征差异显著。

2.3 相关分析：CPV随龄正相关，与LCT及DBS运动改善率负相关；左投射纤维Dx_proj、Dy_proj及右联络纤维Dx_assoc与震颤DBS改善率正相关。形态特征中楔前叶面积、距状沟面积、颞横回面积、角回体积、颞横回体积等与总运动DBS改善率正相关；左侧伏隔核体积、脑室体积、CSF体积与总运动改善率负相关。

2.4 特征选择与模型验证

2.4.1 总运动预测：临床精选MMSE、Med-OFF MDS-UPDRS III、震颤分项等；类淋巴选右Dx_assoc、左相对CPV；形态选扣带回、距状沟、颞枕区体积/面积/厚度等。单模最大AUC：临床0.689、类淋巴0.653、形态0.830；三模SVM最优AUC 0.850±0.045，ACC 0.787，特异0.696，灵敏0.862。SHAP示左相对CPV、右Dx_proj、Med-OFF严重程度、震颤分项及扣带/距状/颞/枕形态贡献最大。

2.4.2 震颤预测：临床选LEDD、Med-ON震颤分项、LCT改善率等；类淋巴选右DTI-ALPS_r、右Dy_proj、右相对CPV；形态选眶额、颞、膝下扣带、嗅区等。单模最大AUC临床0.682、类淋巴0.704、形态0.849；三模LDA最优AUC 0.861±0.047，ACC 0.751，特异0.737，灵敏0.765。

2.4.3 探索性步态表型预测：UMAP+bisecting K-means将66例分两簇（簇I步长/摆期/步速/步频/跨步/步周期改善更大）。三模SVM AUC 0.970±0.019，ACC 0.834，特异0.700，灵敏0.925；临床+形态SVM AUC 0.971，提示类淋巴增量有限。SHAP主导为额盖、颞、顶枕、后中央、次中央形态及右支撑相、左步长、CPV、DTI-ALPS。

讨论总结

研究人员指出三方面进展：影像特征在常规临床预测外具增量价值，形态单模已近0.83–0.95，但临床变量（Med-OFF MDS-UPDRS III、震颤分项、LCT、MMSE）仍不可替；类淋巴标记呈终点特异性，CPV负向关联运动/震颤改善且入选重要特征，支持其作补充间接生物标记而非主决定因子，DTI-ALPS与总运动改善无直接关联且在步态模型无增量；偏侧化特征反映疾病负荷不对称与感觉运动/视空间网络半球特异，但回顾性设计不能确立因果偏侧。局限含单中心回顾、无独立外验、仅程控早期（约1月）结局、步态亚组小且标签为无监督聚类非临床阈值、未校校准度、高维形态多共线性与扫描仪变异。

结论部分翻译

通过整合基线临床变量与术前形态测量及补充性类淋巴生物标记，本研究提供一影像知情框架，用于预测PD双侧STN-DBS术后程控早期总运动反应、震颤反应及探索性步态相关表型。结果支持多模态术前分层之潜在价值，同时强调在临床实施前仍须外部验证、前瞻性确认与更长期随访。