A08 · Motor cortical representation of speed and direction during reaching
Daniel W. Moran & Andrew B. Schwartz (1999) Journal of Neurophysiology, 82, 2676–2692. doi:10.1152/jn.1999.82.5.2676
二维伸手中,含速度项及速度—方向交互项的模型能够描述相当一部分运动皮层时变活动。单细胞拟合R²中位数为0.68,群体平均可达0.99;后者不能代替单试次预测性能,也不能把约145 ms的最佳相关领先量当作固定神经传导时间。
核对版本:Zotero motor_control_interception_64 集合本地全文(PDF 全文提取,图注与正文互相校订,2026-09-18;图注所标页码均为 PDF 文件页序,非期刊印刷页码);Zotero 条目。
背景
按整个动作求平均能揭示方向调谐,却忽略活动在动作过程中的变化。本文在方向大致保持的点对点运动中,研究速度及其与方向的交互能解释多少时变活动。
研究思路
先对八方向响应求平均,提取方向共同成分,再分析相对偏好方向的差异及试次间速度变化。最后把速度、方向交互与时间偏移纳入同一模型,并构建群体向量与运动轨迹比较。
方法
四只恒河猴在竖直触摸屏上从中心伸向6 cm外的八个目标。1039个符合稳定记录标准的单元进入分析,其中897个来自M1,142个来自一个半球的背侧前运动皮层。每单元八方向、各5次重复,共41560次试次;这不是1039个同步单元的记录。
运动期归一为10个时间箱,并加入15个运动前时间箱。主要分析使用平滑、平方根变换及基线处理后的活动;原始放电也作对照。模型包含常数项、方向无关速度项,以及速度与方向正弦/余弦的交互项;时间偏移通过相关或拟合选择。速度差异主要来自任务允许范围内的自然试次波动。
主要结果
M1单元方向共同成分与速度的最佳相关和时滞具有明显差异。将897个单元平均后,共同成分与平均速度的R²为0.99,最佳领先量为145 ms;EMG平均成分对应R²为0.96、领先65 ms。它们是不同汇总信号的最佳匹配,不直接测量逐单元的传导时间,也不证明每个M1细胞都同样精确地表示速度。
在所测速度范围内,M1群体的速度影响随相对偏好方向变化,支持在描述模型中加入速度—方向交互。所采肌肉的方向分布及速度效应与皮层样本不同,但肌肉取样和竖直运动的重力需求限制了比较范围。
逐单元选择时滞后,平方根变换模型的R²中位数为0.68;不变换时为0.65,统一145 ms时滞时为0.54。按相对偏好方向平均的群体拟合达到0.99。M1群体向量与速度向量的最大场相关为0.97,向量长度与速度的相关为0.94。上述指标对应不同分析层级,不能互换,也不是独立留出数据的单试次解码成绩。
逐图解读
原文 Fig. 1 · PDF 第 4 页
图1 · 记录位置。 汇总四只猴八个半球的穿刺,区分M1与PMd取样。样本量较大,但仍是特定任务和区域的取样,不能称作全运动皮层的普查。
原文 Fig. 2 · PDF 第 5 页
图2 · 平均运动轨迹与速度。 A为八方向轨迹及离散范围,B为速度时程。平均路径接近直线、速度轮廓相似,支持本任务中近似固定方向的分析;平均图不能保证每个试次完全相同。
原文 Fig. 3 · PDF 第 5 页
图3 · 单细胞方向共同成分。 外圈为八方向放电,中心比较预处理后跨方向平均的神经成分与速度。示例在选择最佳时滞后较相似,不是一次试次中的无误差预测。
原文 Fig. 4 · PDF 第 6 页
图4 · 相关程度与时滞的分布。 A/B为M1,C/D为PMd。每个单元选择使相关最大的时滞,因此需结合完整分布判断异质性,不能由群体峰值推论每个细胞具有同一领先量。
原文 Fig. 5 · PDF 第 6 页
图5 · 汇总信号的时间轮廓。 PMd、M1、EMG与手速度分别作平均及对齐比较。不同领先量描述相关时序,不是直接测量PMd依次驱动M1和肌肉的传递链。
原文 Fig. 6 · PDF 第 6 页
图6 · 单细胞方向调谐。 周围为八方向栅格,中心为平均响应与余弦拟合。中心使用平方根变换后的活动,不能把纵轴误写为弧度乘以脉冲频率。
原文 Fig. 7 · PDF 第 6 页
图7 · 方向拟合与偏好方向分布。 A汇总M1与PMd单元的拟合R²,B为偏好方向。未拒绝均匀性只说明该检验未检测到明显偏向,不是严格均匀分布的证明。
原文 Fig. 8 · PDF 第 7 页
图8 · 肌肉的方向和时间响应。 多块肌肉在多个方向参与,其时序和偏好不同。肌肉解剖功能可帮助解释,但单凭EMG不能确认各肌的净关节力矩或唯一协作机制。
原文 Fig. 9 · PDF 第 8 页
图9 · 肌肉调谐及偏好方向。 A为肌电拟合,B/C为两只猴的肌肉方向。多数偏向上方,与皮层样本分布不同;作者将重力需求作为解释,所取肌肉并非所有执行肌的完整集合。
原文 Fig. 10 · PDF 第 9 页
图10 · 按相对偏好方向平均的活动。 A为M1、B为PMd,比较方向共同成分和不同方向差的残余。它展示经过分解及跨单元平均的轮廓,不是每个单元都具有相同的形状。
原文 Fig. 11 · PDF 第 10 页
图11 · 不同速度对应的拟合调谐曲线。 由表中回归结果构建25与35 cm/s的曲线;A/B为皮层,C/D为肌肉示例。两条速度曲线不是额外安排的两组速度指令实验,也不能据此将有限范围关系外推到任意速度。
原文 Fig. 12 · PDF 第 10 页
图12 · 模型拟合的单元分布。 A为变换活动及逐单元时滞,B为原始活动,C统一时滞;R²中位数分别为0.68、0.65、0.54。这是拟合表现,选择时滞本身也使用了数据。
原文 Fig. 13 · PDF 第 11 页
图13 · 时滞与偏好方向估计。 A/B显示模型时滞,C/D比较动态模型和时间平均模型所得方向。两种方法的一致性基于同一批记录,不能当作独立数据集的复现。
原文 Fig. 14 · PDF 第 11 页
图14 · 群体向量与重建轨迹。 外圈是反应期及运动期的向量序列,中心为选定时滞后积分的轨迹。它说明群体平均可形成合理运动学读出,没有直接检验在线控制或未见条件的泛化。
原文 Fig. 15 · PDF 第 12 页
图15 · 向量长度与速度。 A/B为M1和PMd的相关,C/D比较向量长度与方向共同成分的时程。PMd两种汇总并不相同,提醒读者不同群体指标保留的活动成分也不同。
讨论与解读边界
本文提供的是有效的描述模型:方向和速度能同时进入单元活动,并可通过群体读出提取。高群体拟合部分得益于试次与细胞平均;它不证明神经系统只编码这些运动学变量,也不排除力、反馈、任务时序或群体动力学的贡献。
时滞是对整个平滑响应的最佳对齐量,不能单凭PMd、M1、EMG依次领先就确认一条串行因果通路。附录关于共同加性项抵消、乘性项改变向量长度的推导有调谐、归一化和群体分布前提,不适用于任意神经群体。与其他研究的负载加和效应比较,也不足以确定速度和负载经过两条独立神经通路。
思考问题
- 使用独立试次和未用于拟合的速度范围评估,并与力、肌电及任务时序模型比较,速度—方向模型还能保留多少预测能力?
- 同一批同步单元在负载、速度和反馈扰动下的最优时滞及调谐怎样变化?哪些结果能区分运动学描述与具体控制机制?