《PLOS Computational Biology》：Whisker stimulation reinforces a resting-state network in the barrel cortex: Nested oscillations and avalanches
编辑推荐：
大脑皮层即使在缺乏外部刺激的情况下，也处于一种持续不息的活动状态。诸如神经雪崩（neural avalanches）和同步振荡等群体性神经元活动，同样可在静息条件下观察到，并且这些特征被认为有助于感觉信息处理、脑对快速反应的准备状态，以及计算效率。大鼠桶状皮层
大脑皮层即使在缺乏外部刺激的情况下，也处于一种持续不息的活动状态。诸如神经雪崩（neural avalanches）和同步振荡等群体性神经元活动，同样可在静息条件下观察到，并且这些特征被认为有助于感觉信息处理、脑对快速反应的准备状态，以及计算效率。大鼠桶状皮层与丘脑回路因其对来自触须感觉输入的躯体拓扑组织（somatotopic organization），为探索这种相互作用提供了强有力的系统。为表征这些静息态回路，研究人员在大鼠桶状皮层和丘脑中进行了同步多电极记录。在自发活动期间，检测到中心频率约为11 Hz的振荡，与低于4 Hz的慢振荡以及服从幂律分布（power-law distributed）的雪崩同时出现。较低频振荡的相位似乎调制较高频振荡的振幅，并在控制雪崩发生中发挥作用。随后，研究人员在受控触须运动期间记录神经活动，并观察到静息时处于活跃状态的11 Hz桶状回路，确实也是参与触须刺激反应的回路。最后，研究显示，一个由丘脑驱动的放电率模型（firing-rate model）能够描述静息态下观察到的全部现象学特征，并预测桶状皮层对受控触须运动的反应，这提示支撑静息态活动的同一本征动力学，也塑造了感觉反应。
该论文发表于《PLOS Computational Biology》，聚焦于静息态脑活动与具体感觉加工之间的关系，研究对象为大鼠丘脑-桶状皮层（barrel cortex）系统。长期以来，大脑在无外界刺激时仍保持持续自发活动，这类静息态活动不仅消耗大量代谢能量，而且呈现出丰富的群体动力学特征，如神经振荡与神经雪崩。然而，静息态活动究竟是否仅是背景噪声，还是具有为感觉输入预备皮层网络的功能意义，仍缺乏直接证据。既往研究提示，某些静息态活跃回路会在任务或刺激到来时被重新招募并增强，但在具体感觉通路中，尤其是在桶状皮层与触须感觉之间，这一命题仍主要停留在理论层面。基于此，研究人员希望阐明：静息态下表现出的复杂内在动力学，是否与须刺激诱发的感觉反应共享相同的网络基础与动力学机制。
为解决这一问题，研究人员在4只乌拉坦麻醉大鼠中，对桶状皮层及丘脑腹后内侧核区域进行了同步多通道电生理记录，分别分析自发活动与机械控制触须刺激后的诱发活动。研究重点放在两类典型集体动力学现象：一类是局部场电位（local field potential, LFP）中的瞬态振荡，另一类是多单位活动（multi-unit activity, MUA）中的神经雪崩。结果表明，在静息态下，桶状皮层存在以约11 Hz为中心的较快振荡，以及低于4 Hz的慢振荡，两者之间形成明确的相位-振幅耦合（phase-amplitude coupling, PAC），即慢振荡的相位调节快振荡的振幅；与此同时，桶状皮层中的神经雪崩呈幂律分布，且其发生时序也受慢振荡相位调控。更为重要的是，在给予触须刺激后，静息态下已存在的11 Hz振荡显著增强，提示静息时活跃的回路正是感觉刺激到来时被招募的回路。最后，研究人员构建并验证了一个丘脑输入驱动的简化兴奋-抑制群体模型，该模型不仅重现了静息态中的嵌套振荡和雪崩统计特征，也成功预测了刺激后11 Hz振荡增强的现象。由此，论文提出：静息态活动并非无意义的背景噪声，而是一种结构化、富动力学特征的内在网络状态，可为快速而高效的感觉处理提供准备。
本研究主要采用以下关键技术方法：其一，基于4只P25–P35龄Wistar大鼠建立丘脑-桶状皮层同步电生理记录体系，使用32通道线性探针采集LFP与MUA信号，并设置机械化、可控的触须刺激范式；其二，应用经验模态分解（empirical mode decomposition, EMD）与Hilbert-Huang变换分析非平稳LFP信号，提取本征模态函数（intrinsic mode functions, IMF）并量化边际Hilbert谱及PAC；其三，基于MUA离散事件，采用平均事件间隔（inter-event interval, ISI）分箱方法识别神经雪崩，并对其大小与持续时间分布进行幂律拟合；其四，构建丘脑驱动的Wilson-Cowan类兴奋-抑制放电率模型，用实验实测丘脑放电率作为输入，对静息态和刺激后皮层动力学进行统一模拟。
在结果部分，论文首先报告了“Resting state neural activity in the barrel cortex and thalamus”。在这一部分中，研究人员围绕静息态下丘脑与桶状皮层的振荡和雪崩活动展开系统分析。
在“EMD analysis of spontaneous activity oscillations”中，研究人员对LFP进行了EMD分解，重点关注IMF 2与IMF 4两个成分。IMF 2对应频率位于须刺激反应相关范围，中心约11 Hz；IMF 4则为低于4 Hz的慢振荡。结果显示，在静息态下，桶状皮层与丘脑都可见这两类振荡，但IMF 2功率峰在桶状皮层显著高于丘脑，提示11 Hz节律在皮层中更为突出。进一步PAC分析表明，IMF 4的相位与IMF 2的振幅之间存在稳定的非均匀关系，呈近似正弦型依赖，说明快慢振荡之间形成嵌套振荡结构。与打乱试次对应关系构建的保守零模型相比，该耦合明显偏离偶然水平，支持该现象具有真实的神经生理学意义。
在“Avalanche dynamics during spontaneous activity”中，研究人员转而使用MUA分析神经雪崩，以减少阈值化连续信号带来的偏差。结果发现，桶状皮层中雪崩大小与持续时间分布均符合幂律，并且平均雪崩大小与持续时间之间满足标度关系，反映出复杂的跨时间与跨空间组织特征。相比之下，丘脑中的对应分布更符合指数分布，平均雪崩大小与持续时间关系更接近平凡过程，说明具有临界样特征的雪崩动力学主要在桶状皮层中形成，而非直接由丘脑输入预成。进一步地，研究人员计算了雪崩密度随时间的变化，发现桶状皮层雪崩发生概率也受到IMF 4慢振荡相位调制，并且这种调制模式与LFP中IMF 2振幅的相位偏好相一致。该结果支持这样一种解释：快振荡与神经雪崩并非彼此对立，而是同一底层动力学的互补表现，前者体现局部节律性爆发，后者体现跨通道活动级联，而慢振荡则为两者提供共同的时间门控框架。
随后，论文进入“Neural activity after whisker stimulation”，检验静息态动力学是否会在感觉加工时被增强和招募。
在“EMD Analysis of post-stimulation oscillations”中，研究人员对刺激后2 s内的活动进行同样的EMD分析。结果表明，须刺激后桶状皮层中IMF 2功率显著高于静息态，说明约11 Hz振荡被明显强化；而在丘脑中，这一频段的变化并不显著。与此同时，IMF 4低频振荡在刺激前后及丘脑、皮层之间都较为稳定，未显示类似IMF 2的显著增强。该结果直接说明，静息态已经存在的11 Hz皮层节律并非偶然背景活动，而是在感觉输入到来后被进一步放大的内在动力学模式。
在“Avalanche statistics after stimulation”中，研究人员分析刺激后的雪崩统计特征。结果显示，桶状皮层中的雪崩大小分布总体仍保持幂律特征，但在大雪崩区域出现轻度隆起，对应于与刺激相关的增强活动；这些增强雪崩与IMF 2振荡共现。丘脑中则仍主要表现为与静息态类似的指数分布。该结果提示，感觉刺激并未破坏桶状皮层原有的临界样动力学框架，而是在其基础上叠加并放大了特定规模的活动级联。
在“Modelling the neural dynamics in the barrel cortex”部分，论文进一步通过计算模型统一解释实验现象。
在“Oscillatory activity in the model in resting state”中，研究人员基于Pinto和Ermentrout提出、并结合Kyriazi和Simons实验参数校准的简化兴奋-抑制群体模型，对桶状皮层进行建模。模型以实验记录的丘脑MUA放电率作为输入，仅对丘脑到抑制性群体的连接权重进行适度调整，使系统动力学保有以抑制为主导的特征，同时允许出现10–15 Hz范围的极限环（limit cycle）行为。模拟结果显示，在静息态丘脑输入驱动下，模型能够同时生成与实验一致的IMF 2和IMF 4振荡，并重现实验中观察到的PAC关系，说明桶状皮层内在非线性动力学与外部丘脑驱动相结合，即可产生嵌套振荡。
在“Avalanche dynamics in the model”中，研究人员通过将模型输出的皮层放电率转化为人工神经元事件序列，采用与实验完全一致的流程分析雪崩。结果表明，模型能够再现实验中桶状皮层雪崩大小和持续时间的幂律分布，以及平均雪崩大小与持续时间的标度关系；同时，也能重现慢振荡相位对雪崩密度的调制。由此说明，实验中观察到的神经雪崩及其与嵌套振荡的耦合，并不需要额外复杂机制，单凭丘脑时变输入与皮层兴奋-抑制网络的相互作用即可产生。
在“Model response to whisker movement”中，研究人员将模型输入替换为刺激后实测的丘脑放电率，而不修改任何静息态拟合参数。模型成功预测到IMF 2振荡在刺激后显著增强，且统计检验与实验结果一致。该结果具有重要意义，因为它表明静息态与刺激态并不依赖两套独立机制，而是同一动力学系统在不同输入驱动下的不同表现形式。
讨论部分围绕静息态活动的功能意义、振荡与雪崩的统一框架，以及模型的生理学解释展开。研究人员指出，桶状皮层在静息态下已存在11 Hz附近的瞬态振荡与低于4 Hz的慢振荡，二者之间存在明确的PAC，同时伴随幂律分布的神经雪崩。这些复杂而富结构的动力学特征，在须刺激后并未消失，而是表现为11 Hz振荡增强及相关大雪崩的增加，因此支持“静息态回路为感觉处理提供预备状态”的观点。研究提出，皮层并非持续驻留于固定临界点，而是在时变丘脑输入作用下动态徘徊于分岔附近，并短暂越过Hopf分岔阈值，从而产生间歇性振荡与临界样雪崩行为。该框架有助于化解“振荡是有特征尺度现象，而雪崩是无特征尺度现象”的表面矛盾：由于生理振荡本身是瞬态和间歇性的，其爆发持续时间具有随机性，因此完全可以与幂律分布的雪崩共存。
论文也讨论了11–15 Hz频段与纺锤波样活动（spindle-like activity）的关系，指出该频带既往曾在睡眠、发育及麻醉状态下的桶状皮层中被观察到。虽然本研究未进行正式纺锤波检测，但结果提示这一频段的集体活动可能与须相关感觉处理密切相关。研究同时承认麻醉状态可能影响低于4 Hz的慢振荡，并指出未来需要在清醒动物中验证这些发现。不过，11 Hz振荡在刺激后显著增强这一现象，支持其并非单纯麻醉伪影，而更可能是与感觉输入相关的潜在动力学模式。
研究结论部分可概括为：桶状皮层静息态活动中存在由慢振荡调制的11 Hz瞬态振荡及幂律型神经雪崩，这些动力学特征在须刺激后被进一步增强；一个由丘脑输入驱动的简化皮层模型能够统一再现静息态与刺激态中的主要现象，表明塑造静息态活动的同一本征动力学也参与并组织感觉反应。整体而言，该研究支持这样一个结论：自发神经活动并非闲置背景，而是具有组织性的动力学网络状态，可使桶状皮层为触须感觉输入的高效加工做好准备。
