MRI引导的个性化tDCS电极优化增强卒中后吞咽运动皮层电场靶向：基于有限元建模的三种电极配置比较研究

《Frontiers in Neurology》：MRI-guided personalized tDCS planning for target engagement of the swallowing motor cortex after stroke: a finite-element modeling study

　　

背景：卒中后吞咽困难（Poststroke Dysphagia, PSD）是卒中后常见且后果严重的并发症，影响高达50%的患者，并可导致多种内科并发症。经颅直流电刺激（transcranial Direct Current Stimulation, tDCS）的治疗效果取决于电极放置方案和个体间解剖差异。使用基于10–20脑电图（Electroencephalography, EEG）坐标的标准C5/C6放置来引导电极定位可能导致与目标位置的偏差。目的：为探究计算建模是否相较标准电极放置方法能够改善电场（Electric Field, EF）参数，以及较小尺寸的电极是否能产生更聚焦的结果，研究人员比较了三种以皮层吞咽区（Cortical Swallowing Area, CSA）为靶点的tDCS电极配置方案的EF值：（i）优化圆形方案（Optimized Circular, Opt-C，半径1.5 cm），（ii）基于EEG的传统圆形方案（Conventional Circular, Conv-C，半径1.5 cm），以及（iii）传统方形方案（Conventional Square, Conv-S，5×5 cm）。方法：分析了50例PSD患者的结构MRI数据。基于tES LAB 4.0平台生成的三维四面体头部模型，采用有限元方法（Finite Element Method, FEM）分析模拟EF分布。在不同电极配置方案间，比较了EF强度（定义为目标感兴趣区内各节点EF的L2范数E_norm的均值）和目标特异性聚焦度（定义为目标感兴趣区内E_norm值不低于全脑灰质和白质E_norm分布第75百分位数的四面体单元总体积）。结果：三种tDCS电极配置方案间EF强度（p<0.001）和目标特异性聚焦度（p<0.001）均存在显著差异。Opt-C方案产生了最高的EF强度（中位数：0.3525 V/m；四分位距：0.306–0.455）和目标特异性聚焦度（中位数：104.69 mm³；四分位距：88.17–120.31）。此外，小圆形电极相较方形电极产生了显著更高的EF强度和目标特异性聚焦度（p<0.001）。结论：MRI引导的个性化tDCS规划改善了模型中传递至皮层吞咽靶点的电场，并揭示了在重度PSD患者中与传统C5/C6放置的显著偏差。这些发现支持以图像引导靶向参与作为PSD精准神经调控的框架，但在推断其治疗效能之前仍需前瞻性的生理学和临床验证。

**论文解读**

**一、研究背景与问题**

卒中后吞咽困难（Poststroke Dysphagia, PSD）是卒中急性期影响约42%患者的常见并发症，常与吸入性肺炎、营养不良、脱水和住院时间延长相关，显著增加死亡率和医疗成本。尽管饮食调整和吞咽康复等传统干预措施仍是标准治疗方案，但部分患者即便接受强化治疗后疗效仍不理想，这凸显了对更有效、更易推广的干预手段的需求。通过无创性脑刺激增强神经可塑性已被证实可改善吞咽功能，其中经颅直流电刺激（transcranial Direct Current Stimulation, tDCS）通过调节皮层兴奋性，在与标准吞咽治疗联合使用时显示出改善预后的潜力，且因其无创性和良好的耐受性，即使是重症患者也可适用。

然而，tDCS面临若干固有局限。其一，与重复经颅磁刺激（repetitive Transcranial Magnetic Stimulation, rTMS）等靶向性更强的技术相比，tDCS通常使用头皮电极传递弱电流，产生的电场（Electric Field, EF）分布广泛、空间聚焦性差，常超出预定靶区范围。其二，常规使用的大尺寸电极（5×5 cm或5×7 cm）产生弥散性EF，聚焦性差。其三，个体间解剖差异——包括颅骨厚度、头皮电导率、脑形态学和病灶特征——显著影响tDCS的有效性，限制了将电流持续传递至运动皮层靶区的能力。

这些问题在PSD中尤为关键。皮层吞咽区（Cortical Swallowing Area, CSA）位于初级运动皮层（Primary Motor Cortex, M1），其位置和半球偏侧化存在显著的个体间差异。利用单脉冲TMS研究显示，CSA定位于外侧中央周围皮层，呈现口–咽–食管躯体分布梯度，对应约1.5–2.0 cm直径的头皮区域。CSA与Broca区及岛叶皮层等功能不同的区域毗邻，传统电极配置可能导致对邻近区域的非预期刺激。因此，依赖传统10–20脑电图坐标引导电极放置可能导致峰值EF指向与吞咽无关的邻近运动区。此外，靠近CSA的颞骨附近的颅骨厚度等解剖特征的个体差异可导致EF强度的显著变异。在诸多解剖因素中，脑脊液（Cerebrospinal Fluid, CSF）的影响尤为突出：由于其电导率远高于颅骨和头皮，CSF作为低电阻通路分流并广泛重新分布电流于皮层表面，使其分布成为决定电流扩散和个体间变异的主要因素。

综上，既往研究存在三个尚未解决的关键问题：第一，基于MRI的个体化计算建模是否能改善EF靶向性，超越传统10–20 EEG坐标引导的电极放置效果；第二，小圆形电极是否能相较传统大方形电极提高刺激精度；第三，以预定义的吞咽运动感兴趣区（Region of Interest, ROI）为基础的目标特异性聚焦度指标，是否比全局性EF扩散指标更具有临床相关性。

**二、研究内容与意义**

研究人员开展了一项多中心、回顾性的计算机模拟（in silico）研究，纳入来自三个医疗中心的50例经视频荧光吞咽造影检查（Videofluoroscopic Swallowing Study, VFSS）确诊为重度PSD且需管饲的患者的结构MRI数据（平均年龄73.6±10.7岁，54.0%为男性，62.0%为左侧半球病变）。利用tES LAB 4.0平台，基于深度神经网络进行卒中病灶分割并生成个体化头部模型，采用有限元方法模拟三种tDCS电极配置方案的EF分布：（i）优化圆形方案（Opt-C，半径1.5 cm），通过网格搜索算法在传统C5/C6位置附近搜索最优阳极位置，以最大化目标ROI内的EF；（ii）基于10–20 EEG系统的传统圆形方案（Conv-C，半径1.5 cm），阳极置于未受损半球的C5或C6位置；（iii）传统方形方案（Conv-S，5×5 cm），电极位置同Conv-C。所有方案中阴极均固定于受损半球的眶上区（AF3或AF4），刺激电流均为2 mA。

研究结果表明，Opt-C方案在目标ROI内产生了显著更高的EF强度（较Conv-C提高8.8%，较Conv-S提高18.2%）和更大的目标特异性聚焦度，且小圆形电极在所有病例中均优于方形电极。电极位移（优化位置与常规位置之间的三维欧氏距离）与EF增益呈显著正相关，表明个性化优化的获益在吞咽相关靶点与标准C5/C6位置解剖错位较大的患者中最为明显。该研究首次为PSD的MRI引导个性化tDCS电极优化提供了计算层面的证据，建立了精准靶向方案的生物物理学理论基础，为后续前瞻性临床验证奠定了基础。

**三、关键技术方法**

本研究使用tES LAB 4.0平台，该平台区别于常规模拟工具的关键特征是能将深度学习驱动的卒中病灶分割直接整合入头部建模流程，从而准确表征受损组织电导率的改变。模拟流程分为三步：（1）MRI分割：利用内置深度神经网络分割模块将T1加权高分辨率三维MRI自动分割为不同组织层；（2）解剖建模：采用Delaunay三角剖分法生成三维四面体网格，在皮肤层密集分布节点以精确表示表面电极位置；（3）EF估算：通过有限元方法离散化tDCS问题，为各组织层（灰质、白质、CSF、颅骨、皮肤及受损组织）赋予各向同性电导率值，在准静态近似下求解Laplace方程，使用共轭梯度法计算电位分布，进而通过一阶Lagrange单元线性插值导出节点级方向性EF分量。目标CSA基于已发表的任务态功能磁共振成像研究中自主吞咽激活区的汇聚证据定位于中央前回，操作化为2 mm半径的球形ROI。EF强度量化为目标ROI内各节点E_norm的均值，目标特异性聚焦度量化为ROI内E_norm不低于全脑灰白质第75百分位数的四面体单元总体积（mm³）。样本为50例来自韩国三家医疗中心（2013年1月至2023年7月）的经VFSS确诊的重度PSD患者。统计方法采用Friedman检验和经Bonferroni校正的Wilcoxon符号秩检验进行配对比较，并采用线性混合效应模型校正病灶体积、侧别和位置的潜在混杂效应。

**四、主要研究结果**

**4.1 患者描述性特征**

50例患者中，27例（54.0%）为男性，平均年龄73.6±10.7岁，31例（62.0%）病灶位于左半球，平均病灶体积54.3±58.8 cm³。病灶部位分布为：皮层19例（38%）、皮层下28例（56%）、混合性3例（6%）。所有患者在MRI采集时均经医疗记录确认存在需管饲的重度吞咽困难。

**4.2 三种电极配置方案的EF强度与目标特异性聚焦度比较**

Friedman非参数检验显示，三种方案间EF强度（χ²(2)=100，p<0.001）和目标特异性聚焦度（χ²(2)=53.47，p<0.001）均存在显著差异。Opt-C方案产生最高的中位EF强度（0.3525 V/m），其次为Conv-C（0.3240 V/m）和Conv-S（0.2982 V/m），所有配对比较均显著（经Bonferroni校正后p<0.001）。Opt-C较Conv-C中位EF提高0.0285 V/m（8.8%），较Conv-S提高0.0543 V/m（18.2%）。目标特异性聚焦度同样以Opt-C最高（中位104.69 mm³），其次为Conv-C（100.80 mm³）和Conv-S（98.31 mm³）。圆形电极在72%（36/50）的患者中表现出优于方形电极的目标特异性聚焦度，而无一例出现方形电极优于圆形电极的情况。在同时校正病灶体积、侧别和位置的线性混合效应模型中，优化方案的优越性保持不变：Opt-C较Conv-C的EF强度校正均值差为0.0447 V/m（95% CI 0.0351–0.0544，p<0.001），目标特异性聚焦度差值为7.44 mm³（95% CI 4.08–10.80，p<0.001）。未发现电极配置方案效应因病灶特征而异的证据。在个体层面，Opt-C在所有皮层下病变（28例）、皮层病变（19例）及混合病变（3例）患者中均产生了高于Conv-C的EF强度。

**4.3 电极位移与刺激参数的相关性分析**

优化阳极位置与传统阳极位置之间的电极位移中位数为20.6 mm（四分位距15.1–27.2 mm），最大值为36.5 mm，相当于约2–3个网格点间距。电极位移与EF强度变化（r_s=0.651，p<0.001）及目标特异性聚焦度变化（r_s=0.526，p<0.001）均呈显著正相关。优化阳极距个体化皮层靶点的距离显著小于传统阳极（中位33.9 mm vs 37.5 mm，p<0.001，配对秩双列相关系数r_rb=?0.545），且优化方案的阳极–阴极间距也更短（中位124.8 mm vs 134.0 mm，p<0.001，r_rb=?0.561）。值得注意的是，在13例优化阳极距靶点更远的患者中，目标EF强度仍全部增加，且EF增益幅度与阳极距靶点更近者无显著差异（0.032 V/m vs 0.035 V/m，p=0.740）。阳极–阴极间距的缩短与EF增益呈中等正相关（r_s=0.444，p=0.0012），而阳极–靶点距离变化与EF增益无显著关联。总体而言，优化方案在全部50例患者中均提高了目标EF强度。病灶体积与EF强度变化（r_s=?0.173）、目标特异性聚焦度变化（r_s=0.096）及电极位移（r_s=?0.122）均无显著相关性。

**五、讨论与结论**

本研究是一项靶向参与建模研究，其目的在于确定MRI引导的电极优化是否能够超越固定坐标放置方案，提高传递至吞咽运动皮层的EF。结果表明，个性化优化带来的获益主要源于整体电流通路的重构（通过改变阳极–阴极几何关系），而非简单地将阳极靠近皮层靶点。电极位移与EF增益之间的正相关性提示，对于那些吞咽相关靶点与标准C5/C6位置解剖错位较大的患者，个性化建模可能最具价值——固定C5/C6方案可能导致这部分患者的靶点剂量不足。

从方法学角度，在卒中人群中实现可靠EF建模的关键前提是对电导率与健康脑组织显著不同的受损组织进行准确分割。tES LAB 4.0平台通过深度学习驱动的病灶分割模块，生成包含受损组织电学特性改变的个体化头部模型，为FEM模拟提供了更具生物物理学真实性的EF分布预测。此外，本研究采用的目标特异性聚焦度指标不同于传统全局性脑区聚焦度测量，而是直接量化预定吞咽运动ROI内高强度EF的集中程度，与预设的基于ROI的优化目标直接对应。

研究的局限性包括：CSA热点被简化为2 mm半径球形ROI，未能涵盖包括辅助运动区、小脑和脑干在内的更广泛吞咽网络；模型未纳入白质病变、脑室周围改变或病灶特异性电导率变异；未比较优化双电极方案与高精度tDCS（High-Definition tDCS, HD-tDCS）的直接效果；EF强度使用未截尾的EF分布计算均值，未对极端值进行敏感性分析；病灶特征仅限于体积、侧别和粗略分类，缺乏病灶–靶点距离和几何关系数据；以及对计算靶向参与的发现尚未经生理学或临床疗效验证。

研究结论：计算建模结果表明，基于个性化MRI的计算优化可显著改善卒中患者CSA的tDCS靶向性，超越基于固定头皮坐标的传统电极配置方案。通过将包含个体化病灶分割的卒中特异性建模框架与小圆形电极和目标特异性聚焦度指标相结合，研究人员为PSD的精准神经调控提供了一种实用策略。这些发现为未来临床试验提供了计算基础，以确定个体化刺激是否能改善吞咽功能恢复并减少卒中后吞咽困难相关并发症。该建模框架的临床转化正在一项多中心、双盲、假对照随机试验（NCT06305949）中进行前瞻性评估，该试验纳入30例PSD患者，采用与本模拟研究相同的MRI引导个体化靶向方案和1.5 cm半径圆形电极，目前已完招募，吞咽结局分析正在进行中。