《Brain Stimulation》：Brain-State-Resolved Consistency of Corticospinal Responses to TMS
编辑推荐：
刺激前感觉运动mu动态塑造运动诱发电位（MEPs）的试次间变异性。高mu功率与跨参与者MEP变异性降低相关。mu相位在组水平调节MEP波幅，但不调节MEP变异性。研究结果可能为脑状态依赖性经颅磁刺激（TMS）降低反应变异性提供基础。
刺激前感觉运动mu动态塑造运动诱发电位（MEPs）的试次间变异性。高mu功率与跨参与者MEP变异性降低相关。mu相位在组水平调节MEP波幅，但不调节MEP变异性。研究结果可能为脑状态依赖性经颅磁刺激（TMS）降低反应变异性提供基础。
研究背景与问题：经颅磁刺激（transcranial magnetic stimulation, TMS）作用于初级运动皮层（M1）可在对侧目标肌肉诱发运动诱发电位（motor-evoked potentials, MEPs），用于读取皮质脊髓兴奋性。TMS研究的关键挑战是诱发反应存在明显变异性，削弱基于TMS的兴奋性估计的可靠性与临床效用。MEP试次间变异性由方法学（外源性）与生理性（内源性）来源共同构成；即使方法学参数严格控制，经脑电图（electroencephalography, EEG）捕捉的短暂电生理脑状态波动仍贡献反应变异性。既有脑状态依赖性TMS研究多关注组水平不同脑状态间平均MEP波幅差异，通常在状态定义的分箱内对试次取平均以提高信噪比，并将扣除脑状态后的MEP离散度视为随机、不可控的生理噪声。然而，状态内MEP变异性本身可能是感觉运动系统兴奋性状态的有信息特征。MEP变异性与波幅同皮质脊髓输入-输出（input-output, I/O）曲线工作点相关；外部操纵如提高刺激强度或随意肌收缩可移动工作点并可能降低变异性。内源性脑状态是否同样影响MEP变异性尚不清楚。因此，研究人员开展本研究，检验刺激前mu动态是否不仅影响MEP幅度，也影响MEP一致性；以状态内变异性而非状态间平均波幅差异为指标，聚焦mu功率、mu相位及半球间M1–M1功能连接。
研究开展与结论意义：主数据集来自University of Tübingen的29名健康右利手参与者，完成左侧M1单次TMS–EEG；验证数据集来自Aalto University（Finland）的22名健康右利手参与者；另10名参与者完成两日相同会话用于重测。刺激为左侧M1手部运动热点1200个单脉冲，强度110%静息运动阈值（resting motor threshold, RMT）。结果显示，高mu功率与MEP变异性降低相关；mu相位在组水平不调节MEP变异性；半球间M1–M1连接与变异性关系较不稳健。跨站点验证支持mu功率效应。研究人员结论为皮质脊髓变异性并非纯随机，而包含结构化的脑状态依赖性成分；这为发展脑状态依赖性TMS以减少反应变异性、提高神经调控效应一致性提供基础。
关键技术方法（≤250字）：主队列为Tübingen 29名健康右利手参与者，验证队列为Aalto University 22名健康右利手参与者，重测亚组10名。左侧M1手部运动热点给予1200个单脉冲TMS（110% RMT），同步128通道EEG与右侧目标肌肌电图（electromyography, EMG）。以C3表面Laplacian提取刺激前8–13 Hz mu特征；FFT估计单试次mu功率，PHASTIMATE估计mu相位，C3–C4的相位锁定值虚部（imaginary part of phase-locking value, iPLV）量化半球间M1–M1连接。按mu功率/连接八分位或四分位及mu相位四窗分层，计算状态内归一化变异系数（coefficient of variation, CV），结合bootstrap、滑窗、线性回归、重复测量方差分析、错误发现率（false discovery rate, FDR）校正、一致性相关系数（concordance correlation coefficient, CCC）重测与跨站点验证。
研究结果：
3.1 高刺激前mu功率与MEP变异性降低相关：通过按mu功率八分位分层并对每名参与者拟合线性回归，研究人员发现MEP变异性随mu功率升高而降低（斜率–1.62±0.41，t(28)=–3.96，p<.001），75.9%参与者呈负关系；四分位重复测量方差分析显示高功率分箱变异性显著低于低功率分箱（t(28)=3.39，p=.006）且低于随机抽样（t(28)=?2.69，p=.012）；滑窗分析显示该关系在全实验内稳定，替代变异性指标与补充控制分析支持mu频段相关效应。平均MEP波幅随mu功率升高而增加（p<.001）。
3.2 mu相位不调节MEP变异性：四相位窗重复测量方差分析未发现组水平显著效应（F(2.93,82.04)=1.77，p=.160）；mu相位显著调节平均MEP波幅（p=.002）。联合mu功率与相位分层时，mu功率主效应显著（F(1,28)=12.42，p=.001），相位主效应（F(3,84)=0.37，p=.772）及交互（F(3,84)=0.12，p=.951）不显著；替代变异性指标与相位估计参数结果一致。
3.3 半球间M1–M1功能连接对MEP变异性的调节：按mu频段iPLV八分位分层回归，连接增强与MEP变异性降低相关（斜率–0.82±0.25，t(28)=–3.29，p=.003），72.4%参与者呈负关系；高连接分箱变异性低于低连接分箱（t(28)=3.67，p=.002）及随机分箱（t(28)=–4.35，p<.001），滑窗显示全实验稳定；连接对平均MEP波幅无显著影响（p=.055）。联合mu功率与连接分层时，mu功率主效应显著（F(1,25)=5.91，p=.023），连接主效应（F(1,25)=2.92，p=.10）及交互（F(1,25)=0.01，p=.92）不显著。
3.4 mu功率分辨的MEP变异性的跨会话重测信度：10名参与者两日相同方案；重复测量方差分析显示状态主效应显著（F(1.18,10.59)=6.22，p=.027），会话（F(1,9)=0.02，p=.900）及交互（F(1.84,16.52)=0.14，p=.853）不显著。Lin一致性相关系数（CCC）在高功率分箱为0.42（95% CI [–0.08,0.75]）、低功率分箱为0.26（95% CI [–0.32,0.61]），随机分组为–0.42（95% CI [–0.74,0.10]）；置信区间均跨零。
3.5 脑状态分辨MEP变异性的跨站点泛化性：Aalto验证数据中mu功率与MEP变异性仍呈负线性关系（斜率–1.42±0.59，t(21)=–2.40，p=.026），合并数据集为–1.53±0.34（t(50)=–4.49，p<.001）；mu相位在验证集（F(2.45,48.95)=0.54，p=.624）及合并数据（F(2.81,137.77)=0.28，p=.825）均无显著效应。M1–M1连接负关系在验证集不显著（斜率–0.52±0.31，t(21)=–1.68，p=.108），合并数据显著（斜率–0.69±0.19，t(50)=–3.56，p<.001）。
讨论总结：讨论指出，既往TMS–EEG研究已显示刺激前mu功率与相位关联MEP波幅的试次间波动，本研究复现组水平效应并进一步显示mu动态也可调节反应一致性。mu功率与MEP变异性负相关提示内源性皮层动态参与稳定皮质脊髓反应；mu相位未显示组水平一致性效应，提示相位更多影响反应幅度而非一致性，且相位效应短暂、效应量较弱并存在个体差异。机制上，110% RMT位于I/O曲线上升斜率附近、处于非饱和工作区，对状态依赖性变异敏感；高mu功率可能作为内源性调节因素使皮质脊髓系统工作点向更饱和募集区移动，使更多皮质运动神经元被一致驱动至阈值附近，降低内在背景波动对输出的相对影响；这更可能是相对信噪比现象而非绝对抑制神经噪声。较强半球间M1–M1连接可能通过更协调的同步突触输入降低变异性，但该关联对指标及验证集不稳健，且在联合mu功率后减弱。讨论还强调，随机试次分组重测一致性为负，而按脑状态条件化后一致性改善但置信区间宽，提示个体内可重复性仍不确定；状态分辨一致性可能作为反映皮质脊髓募集稳定性的生物标志物。局限包括CV可能残留均值-方差耦合，状态分箱内仍有未测量脑状态及方法学变异，不能将状态分辨变异性完全归因于该脑状态，且变异性估计对预处理与分析选择敏感。
研究结论（原文翻译）：研究人员量化了刺激前感觉运动mu节律动态如何塑造皮质脊髓输出的试次间变异性。MEP变异性主要受mu功率状态调节，而mu相位和半球间M1–M1功能连接在跨参与者中并非MEP变异性的足够强或一致的决定因素。这些发现提示皮质脊髓变异性并非纯随机，而包含结构化的脑状态依赖性成分。总体而言，结果推进了对TMS诱发反应试次间变异性的机制理解，并可能为未来旨在减少未解释反应变异性、提高神经调控效应一致性的TMS方法提供基础。
