耦合神经元振荡系统中频率失谐控制振幅的理论与证据

《PLOS Computational Biology》：Theory and evidence of amplitude control by frequency detuning in a coupled neuronal oscillator system

　　

摘要
产生节律性运动的神经元振荡回路必须灵活且可靠地运作，以产生动物自然表现出的多样化运动模式。啮齿动物节律性地“拂动”其触须（vibrissae）以进行触觉感知，并动态调整拂动范围以服务于不同的感知目标。拂动由与呼吸耦合的脑干振荡回路控制，然而拂动振幅如何被调制仍属未知。本研究提出并评估了一种基于耦合振荡器同步原理的振幅控制机制。具体而言，对大鼠行为数据的重新分析表明，拂动表现出通过高频呼吸模式“嗅探”（sniffing）的夹带（entrainment）实现放大的运动学特征和相位动力学。拂动振荡回路的神经元网络模型表明，拂动振幅可通过将振荡器固有频率相对于嗅探频率进行偏移来调制，类似于工程中的“失谐”（detuning）技术。基于这些结果，研究人员提出耦合神经元振荡器之间的失谐可能代表了神经系统中增益控制的一种通用计算策略。

论文解读

**研究背景与问题**

动物产生多种节律性运动行为，如呼吸、运动、进食和主动感知，这些行为需要神经系统产生时间结构化的多肌群激活以形成协调、目标导向的运动指令。这些基本运动序列被认为由中枢模式发生器（central pattern generators, CPGs）产生，即在没有感觉反馈的情况下也能以正确时序产生节律性运动的神经元网络。尽管这些运动被视为高度结构化且自动化的，但其运动学特征可大幅变化以满足动物的目标。一个核心问题是：CPG的神经元架构如何使其在遵守运动系统生物力学约束的同时，具备适应性地产生动物自然表现出的多样化运动模式的能力？

啮齿动物的触须（vibrissae）拂动（whisking）作为一种探索性行为，与快速呼吸和鼻孔运动（即“嗅探”，sniffing）共同构成探索行为的主要模式。拂动和呼吸在全部呼吸频率范围内高度同步，且拂动CPG已被定位在延髓的触须中间网状结构（vibrissa intermediate reticular formation, vIRt），该区域接收来自附近呼吸振荡回路的输入。尽管拂动的肌肉激活序列高度刻板化，其运动学特征却是动态的：啮齿动物在不同行为情境下改变拂动的频率、振幅和设定点。然而，拂动与呼吸振荡器之间的相互作用如何支持如此宽的动态范围仍属未知。本研究聚焦于啮齿动物神经系统如何动态调制拂动振幅这一具体问题。

**研究内容与结论**

研究人员提出假说：拂动振幅的调制可通过耦合振荡器系统中的“失谐”（detuning）机制实现，即通过调整振荡器固有频率相对于驱动频率的偏移来控制振幅。为验证这一假说，研究人员重新分析了公开的大鼠拂动和呼吸同步记录数据库，识别出与夹带（entrainment）假说一致的运动学特征。由于呼吸始终在进行，vIRt网络在离体状态下的固有动力学难以在实验动物模型中直接观察，研究人员构建了一个基于生理约束的vIRt拂动CPG神经元网络模型，并模拟了该模型对不同频率呼吸控制信号的响应。结果显示，模拟拂动振幅在vIRt固有频率与呼吸频率匹配时达到最大，且产生的拂动信号具有与行为观察相似的关键运动学特征和相位动力学。这些发现表明频率失谐是实现耦合神经元振荡系统振幅调制的生物学合理算法。该研究发表在《PLOS Computational Biology》。

**关键技术方法**

研究人员使用自适应指数整合发放（Adaptive Exponential Integrate-and-Fire, AdEx）模型神经元构建vIRt拂动CPG的神经元网络模型，模型包含preB?tzinger复合体（preB?tzinger complex, PB）抑制性吸气神经元群、vIRt回缩（vIRt-retraction, ）和伸出（vIRt-protraction, ）抑制性神经元群以及面神经核振动毛运动神经元（facial motor neurons, ）四组神经元群。神经元数量基于先前实验观察估算。模型采用基于电流的突触连接，突触电流呈指数衰减。效应器模型通过将运动神经元群体放电活动与alpha函数进行卷积产生振动毛位置信号。输入电流采用Ornstein-Uhlenbeck过程模拟突触噪声。行为数据来自公开发表的大鼠（Long-Evans品系）同步拂动与呼吸记录数据集，共7只动物25个session，467个epoch。使用多窗谱方法进行频谱和相干性分析，采用Wilcoxon符号秩检验和Student's t检验进行统计推断。

**研究结果**

**拂动与呼吸相关的运动学特征**
对大鼠同步拂动和呼吸数据的定量重新分析表明，拂动和呼吸信号具有重叠的功率谱。在嗅探期间，拂动和呼吸频谱对齐，主导拂动频率为6.1 Hz；而在基础呼吸期间，主导拂动频率更高（7.8 Hz）且频谱更宽，提示基础呼吸期间vIRt固有主导频率与典型嗅探频率存在偏移。对基础呼吸期间连续拂动振幅的分析显示，吸气伴随的拂动振幅大于后续“间隔拂动”（intervening whisks），且后续拂动振幅在前4个周期内单调衰减。前4次拂动的第N次振幅与第N+1次振幅显著相关（Pearson相关系数r = 0.23–0.33, p < 0.001），间隔拂动振幅在5–6个周期后达到稳态。从基础呼吸向嗅探转变时，拂动振幅在第3次拂动后达到稳态。96%的嗅探段在p < 0.05水平表现出显著的嗅探-拂动相干性。嗅探频率与拂动振幅无显著相关（r = 0.05, p = 0.27），但拂动-嗅探相干相位（_φ）与拂动振幅显著相关（circular-linear correlation r = 0.41, p ≈ 0），在振幅范围10–80°内相位变化-0.57弧度。

**拂动CPG的自适应指数整合发放模型**
研究人员构建了包含200个PB抑制性吸气神经元、800个vIRt回缩抑制性神经元、400个vIRt伸出抑制性神经元和100个振动毛运动神经元的网络模型。模型采用稀疏随机连接，抑制性突触电流时间常数约10 ms。vIRt神经元表现出短时程适应（short-term adaptation），这是模型半中心振荡器网络中相变的基础机制。在无外部周期驱动时，模型产生单一频谱峰位于6.4 Hz的节律性效应器振荡。模型神经元的CV2值聚集在1附近，与先前清醒小鼠和红藻氨酸诱导拂动大鼠的vIRt单单元记录一致，表明波动驱动的发放机制。

**模型对周期性PB刺激在不同频率下的响应**
在1.25 Hz低频刺激下，效应器信号主导频率（5.8 Hz）保持在固有频率（6.4 Hz）附近。在6 Hz刺激下（接近模型vIRt固有频率），效应器输出振幅和频谱纯度相对于基础呼吸频率刺激增加，表现出频率选择性放大。在10 Hz刺激下，输出信号振幅降低且功率谱呈双峰，主导频率接近驱动频率（9.3 Hz）。驱动拂动振幅在vIRt固有频率附近出现峰值（对应1:1模式），固有频率整倍频比处出现较低振幅峰（2:1、3:1等）。对1–10 Hz PB脉冲频率的相干相位分析显示，在1:1放大频率范围（约4–9 Hz）内，拂动-到-PB相位随频率近似线性偏移。增加PB到vIRt神经元的突触驱动（_PB→P）可加宽1:1放大频率范围并降低相位偏移斜率。Monte Carlo模拟（30次随机连接矩阵实例化）表明，vIRt固有频率附近的振幅放大峰对不同随机连接拓扑高度稳健。

**vIRt固有频率对输入电流统计量的依赖**
为探究固有频率能否通过下行神经元输入改变，研究人员对vIRt神经元施加具有不同均值（μ）、标准差（σ）和时间常数（5 ms）的Ornstein-Uhlenbeck过程噪声输入电流，同时对运动神经元施加恒定电流F。系统变化μ、σ、F值并计算稳态固有频率、振幅和设定点。升高中μ可提高固有频率但也影响振幅和设定点。主成分分析（principal component analysis, PCA）表明，在保持振幅（15–30°）和设定点（25–75°）在目标范围内的参数组合中，单一主成分（PC1）解释了参数空间中超过80%的方差，μ、σ、F与PC1的相关系数分别为0.97、0.83和0.96，且PC1与固有频率正相关（r = 0.75, p < 0.001）。这表明通过同时增加输入电流的均值和宽带波动幅度以及运动神经元的去极化水平，可以在保持振幅和设定点不变的情况下控制vIRt固有频率。

**拂动振幅与PB/vIRt频率失谐的关系**
选取8组单调递增的μ、σ、F参数集，产生固有频率范围为5–9 Hz、固有振幅和设定点均在目标范围内的振荡。对每组参数在不同PB脉冲频率（1–10 Hz）下进行频率扫描模拟，结果显示最大驱动拂动振幅在PB驱动频率接近vIRt固有频率时实现。放大因子（amplification factor, AF）定义为驱动拂动振幅与固有振幅之比，其最大值对应的PB驱动频率与vIRt固有频率紧密对齐。放大因子随频率失谐（_Δf = _fPB ? _fintrinsic）近似线性变化，在±3 Hz范围内AF约为1.5–3.5倍。在_Δf在-3至3 Hz范围内，相位-振幅曲线呈侧向“V”形：_Δf > 0时斜率为负，_Δf < 0时斜率为正。模型预测的相位-振幅依赖关系与行为观察（图2l）一致，表明大鼠在嗅探期间可能主动维持vIRt固有频率略低于嗅探频率以实现大幅度探索性拂动。

**模拟拂动与呼吸相关的运动学特征**
模拟30秒1.25 Hz基底呼吸后转变为6 Hz嗅探的协议。在30次Monte Carlo模拟中，模拟拂动表现出与自然拂动一致的运动学特征：基础呼吸期间连续拂动振幅间存在显著相关（前5个周期）；拂动节律对基础呼吸脉冲表现出强相位重置；间隔拂动振幅依次衰减，在8个拂动周期后收敛至约19°的稳态值（接近该网络的固有振幅18°）；从基础呼吸向嗅探转变后，前3个振荡周期的拂动振幅显著增加。这些模拟运动学特征与自然大鼠探索性拂动观察广泛一致，支持拂动与嗅探之间存在夹带和频率依赖性放大的假说。

**讨论与结论总结**

**讨论部分总结**：本研究以大鼠拂动为模型系统，研究了神经系统如何通过相互作用的神经元振荡回路控制节律性运动振幅。行为数据重分析表明基础呼吸期间拂动频率分布较宽且峰值略高于典型嗅探频率，提示vIRt固有频率具有可变性可被主动调制。模型预测神经元适应赋予网络节律发生能力和跨多个振荡周期的扰动效应持续性。行为观察中拂动-嗅探相位与嗅探频率无显著相关但与振幅相关，支持大鼠主动共调制PB和vIRt频率的假说。模拟表明仅通过改变输入电流即可在保持振幅和设定点不变的情况下改变vIRt固有频率。模型还预测拂动-嗅探相位随拂动振幅增加而减小，与行为观察一致。基于解剖学证据，研究人员提出延髓网状结构（medullary reticular formation, MY-Rt）可能是提供混合兴奋/抑制性下行驱动的来源，初级运动皮层（M1）和中脑网状结构（MB-Rt）是调制vIRt固有频率的候选脑区。

**研究结论**：基于上述发现和解剖学观察，研究人员提出大鼠能够相对独立控制探索性拂动的频率、振幅和设定点的机制：（1）频率调制通过直接调制呼吸振荡回路中的嗅探频率实现；（2）振幅调制通过下行兴奋性驱动至MY-Rt中混合兴奋/抑制性神经元群，使vIRt固有节律与呼吸失谐实现；（3）设定点调制通过对面部运动神经元的兴奋性输入实现。频率、振幅和设定点的分离控制能够支持在啮齿动物中观察到的全部灵活拂动行为。本研究发现支持频率差异不仅是噪声，而是控制神经元同步性方式的一般概念，将频率失谐和夹带的概念扩展到运动控制领域，证明两个网络级振荡回路间的频率失谐可实现节律性运动振幅的模拟控制。