其他
导读人类的许多行为都是由在不同时间尺度上发生的共同过程所支配的。标准的事件相关电位分析假设有关实验事件的响应持续时间是固定的。然而,最近对动物的单个单元记录显示,在需要灵活计时的行为中,神经活动尺度跨越了不同的持续时间。本研究采用了一种通用的线性建模方法,使用固定持续时间和可变持续时间回归器的组合,以分离人类脑磁/脑电图(M/EEG)数据中的固定时间和缩放时间成分。研究者利用这一点来揭示人类头皮记录电位在四个独立EEG数据集中的一致时间尺度,包括时距知觉、产生、预测和基于价值的决策。时间尺度响应的试次间变化预测了被试反应时间的试次间差异,证明了这种时间尺度信号与行为时间变化的相关性。本研究结果为研究人脑中的灵活计时行为提供了一种通用的方法。前言现实世界中的行动和感知需要灵活的时间安排。我们可以快走或慢走,可以识别以不同节奏播放的同一首音乐,并在长短时间间隔内形成时间预期。在许多认知任务中,反应时间的变异性是根据内部证据积累来建模的,因此相同的动态过程在不同的试次中以不同的速度展开。灵活的时间安排在我们的生活中至关重要,然而,尽管经过几十年的研究,它的神经关联仍然受到广泛的争论。脑磁图和脑电图(M/EEG)由于具有高时间分辨率的特点,因而在理解计时的神经基础方面发挥了特别重要的作用,通常用于分析这类数据的方法是事件相关电位(ERP),它可以在多次重复刺激中进行平均事件锁定反应。例如,这种方法已被用于识别在一定时间间隔内的缓慢负向信号的存在。这个信号被称为关联性负变(CNV),被认为与时间有关,因为它的斜率与时间间隔的持续时间成反比。重要的是,ERP分析策略隐含地假设神经活动以相对于实验事件的固定时间延迟发生。然而,最近的研究表明,单个神经元水平上的大脑活动可以用时间尺度模型进行最好地解释,在该模型中,活动可以根据产生的间隔长度被拉伸或压缩的单个响应来解释。当猴子被提示产生不同长度的间隔时,时间尺度模型解释了内侧额叶皮层单个单元的神经反应的大部分差异。这表明,灵活的计时机制是通过调节一个共同神经过程的速度来实现的,这一观点可以从动力学系统理论的角度来看。与动力学系统在一系列神经计算中所起的广泛作用一致,最近对神经集群的研究表明,时间偏差是许多不同集群记录和行为任务的共同特性。例如,在感官和基于价值的决策过程中,时间尺度也隐含在证据整合的神经关联中(这在先前的研究中也被提出作为时间估计的机制)。成功地表征人类的时间尺度成分,可以为研究时间尺度在更复杂的层级任务中的作用打开大门,如音乐创作或言语感知,以及在时序受损的患者群体中。然而,目前尚不清楚神经反应的时间尺度是如何在人类头皮上(使用非侵入性记录)显现出来的。这是因为ERP分析策略的固定时间特性。而且ERP的一个成分,称为CNV,其在不同的时间间隔内以不同的速度变化,这表明存在时间尺度。但关键是,由于叠加问题,任何缩放的活动都将在头皮上与固定时间的成分混合。因此,本研究开发了一种方法来分离EEG中的缩放时间和固定时间成分(图1A)。本研究提出的方法建立在最近开发的基于最小二乘回归的方法之上,这些方法已被证明有助于分离相互重叠的固定时间成分,如刺激相关的活动和反应相关的活动。为了克服叠加问题,这些方法使用一般线性模型(GLM)以去卷积潜在重叠的神经反应。在此基础上,使用GLM来估计固定时间的ERP,其中设计矩阵中填充了时滞“stick函数”(一个回归量,在感兴趣的时间点周围的值为1,在其他时间点为0)。重要的是,stick函数可以捕获潜在神经响应中的重叠(图1B),并且可以通过向模型估计中添加正则化惩罚来提高对神经数据的拟合度。在没有任何重叠的情况下,GLM将完全返回传统的ERP。图1.基于回归的模拟数据分解成功地复原了缩放时间和固定时间成分。这里引入的关键创新是,除了固定持续时间回归器外,此模型中还允许使用可变持续时间回归器,以测试缩放时间响应的存在。特别是,该模型允许stick函数的持续时间根据刺激和反应之间的间隔而变化,这意味着相同的神经反应可以在不同的试次中跨越不同的时间。因此,所提出的方法不是对每种条件的平均间隔持续时间进行建模,而是捕捉试次到试次反应的变异性。返回的缩放时间电位是刺激和反应之间的时间流逝百分比的函数。研究者在单个EEG传感器上模拟了间隔计时任务的数据,该任务由两个固定时间成分(锁定线索和反应)和一个跨越线索和反应之间的缩放时间成分组成(图1A)。本研究提出的方法成功地复原了所有三个成分(图1C),而传统的ERP方法掩盖了缩放时间成分(图1A)。重要的是,在真实的EEG数据中,研究者在所有传感器上重复了这种方法,潜在地揭示了固定时间与缩放时间成分的不同头皮分布(以及不同的神经源)。通过分离固定和缩放成分,本研究方法超越了以往处理EEG实验中时序变化的方法。例如,可以通过平移整个波形或单个ERP成分来调整EEG试次的事件相关时序。另一方面,原始EEG可以通过升采样/降采样或动态时间变化将试次与一个共同时间框架对齐。然而,这些方法并不是设计来分解固定时间和缩放时间成分的。最后,可以使用另一种基于回归的方法——时间响应函数(TRF)来量化连续变量对EEG的影响。TRFs在捕捉不同类型的延迟活动方面特别灵活。本研究方法与TRFs有关,因为它涉及到待估计输入信号与连续回归量的卷积。然而,与TRFs不同的是,本研究方法包含一个额外的缩放步骤,其中输入信号被拉伸或压缩(图1)。时频分解也可以很容易地分离高频和低频响应。事实上,大量的神经振荡都与时间知觉有关。人们可能会合理地预期拉伸/压缩信号在时频域中的表现不同。然而,与本研究所提出的方法不同的是,时频分解并不容易用来寻找缩放时间响应的时间尺度。如果信号在同一频段,时频分解也无法将固定时间响应与缩放时间响应分开。本研究分析了四个独立的EEG数据集,包括三个间隔计时任务和一个决策任务。在第一个任务中,参与者根据提示生成目标间隔(短、中或长)(图2A)。参与者能够根据提供的反馈做出更接近目标时间间隔的反应。在第二个任务中,参与者评估了计算机生成的时间间隔(图2B)。生成的时间间隔越接近目标时间间隔,参与者越有可能判断反应是“准时”的。在第三个任务中,参与者根据节奏预测对即将发生的事件进行时间预测(图2C)。图2.三种时间估计和一种决策范式的任务示意图。在第四个任务中,参与者在成对零食中进行选择(图2D),在这个过程中,反应时间的变化可以被建模为随时间推移的内部证据积累过程,与证据积累相关的神经活动可在头皮上测量。快速、简单试次EEC的增长速度快于缓慢、困难试次EEG的增长速度,这表明内部证据积累的速率更高。因此,本研究预测EEG将包含一个与不同证据积累速率相关的潜在尺度成分。在所有四项任务中,可以观察到一个与前后固定时间成分不同的缩放时间成分(图3),这类似于传统的ERP。通常,ERP成分由其极性和头皮分布来定义。观察到的缩放时间成分具有共同的极性(负性)和头皮分布(中央)。在每项任务中,基于聚类的置换检验显示,缩放时间成分与零显著不同。这些差异是由产生任务中36-87%(P<0.001)、感知任务中42-100%(P<0.0001)、预测任务中18-27%(P=0.004)和决策任务中36-55%(P<0.001)的团簇驱动的。图3.在所有四种范式中都一致地观察到了缩放时间成分,其头皮地形图与固定时间成分不同。为了进一步验证本研究方法,研究者通过计算每个任务和参与者的缩放指数来量化时间尺度(图4)。为了计算这一点,拉伸/压缩每个epoch,以匹配每个任务中的最长间隔,按条件求平均,然后使用较短间隔的拉伸版本计算用于预测较长间隔的决定系数。首先对原始数据执行此操作,然后分别对仅包含固定时间成分和缩放时间成分的数据执行此操作。在所有四项任务中,缩放成分的缩放指数都超过了固定成分的缩放指数。图4.未混合信号在其缩放程度上存在定量差异。然后本研究考察了缩放时间成分与行为变异性的关系:缩放时间成分的潜伏期能否预测被试的反应时间?本研究关注的是时间产生和决策任务,其中间隔持续时间等于响应时间。由于不同试次的响应时间不同,建模的缩放成分也不同。为了测量成分潜伏期,本研究使用主成分分析(PCA)对时间产生任务中中央电极上的延迟活动进行建模。与简单的峰值检测不同,PCA可以解释一系列波形动态。首先回归出GLM识别的固定时间成分,得到的数据集只包含剩余的缩放时间活动。然后,计算了三种时距条件下的平均缩放时间活动(图5A-C)。将PCA分别应用于每种时距条件,并一致地揭示了与缩放时间成分形状匹配的第一个主成分和与它的时间导数匹配的第二个主成分。这一分析证实了本研究数据中缩放时间成分的存在,因为它是去除固定时间成分后的第一主成分。至关重要的是,添加或减去第二个主成分可以捕获缩放时间成分的潜伏期变化。在反应时间分位数中,可以发现PC2得分与反应时间显著相关(图5D)。这意味着,缩放时间成分达到峰值的时间越早,被试对该试次的反应就越快。该结果在决策任务中得到了验证(图5E)。图5.缩放时间成分的变化预测了时间估计中的行为变化。方法模拟。使用MATLAB