清醒态α爆发：外侧膝状体丘脑皮层细胞模型中的动态工作状态

《PLOS Computational Biology》：Awake alpha bursting emerges as the dynamic working state in a lateral geniculate thalamocortical cell model

　　

已知清醒态丘脑能够在有无外部调制的情况下过滤初级感觉输入。研究人员通过构建并分析一个新型外侧膝状体丘脑皮层神经元计算模型，证明了感觉视网膜输入的处理如何受到潜在丘脑动态状态的影响。该模型仅使用经公开数据集表达数据验证的电流，是首个能够产生五种实验确立的不同动态放电模式（dynamic firing regimes）的模型。研究人员证明，丘脑皮层细胞响应来自皮层的谷氨酸能信号或来自脑干的胆碱能 arousal 信号而在这些动态状态之间转换。研究人员重点关注与清醒态丘脑α节律相关的模型动态状态中的信号处理，发现视网膜输入产生丘脑 spike 的能力是视网膜 spike 时序、M电流（M-current）施加的兴奋性间断以及L型钙电流（L-type calcium current）衰减时间之间的平衡。最后，研究人员探索了这两种电流如何帮助丘脑处理视网膜外节律性输入，表明该模型能够夹带（entrainment）至较慢的抑制性和兴奋性节律，并详细说明了将较快频率节律嵌套于慢周期内对于成功的丘脑传递的重要性。研究结果表明，清醒态α节律通过充当这两种电流相互作用的标志而间接具有因果性。这一生物物理约束的外侧膝状体丘脑皮层细胞模型产生了关于不同 arousal 状态下节律动力学、丘脑对视网膜-膝状体传递（retinogeniculate transmission）的控制，以及精神分裂症等神经疾病可能对丘脑处理产生影响的预测。该模型的变体可用于探索高阶丘脑核团的功能，从而将其应用扩展到研究更复杂的认知过程。

清醒态丘脑α节律：外侧膝状体丘脑皮层细胞模型中的动态工作状态与信号处理机制解读

一、研究背景与意义

丘脑长期以来被视为感觉信息向大脑皮层传递的被动中继站。然而，近年来的实验研究彻底改变了这一认知，表明丘脑实际上是一个具有复杂动态特性的主动处理枢纽。背外侧膝状体核（dorsal lateral geniculate nucleus, LGN）是这一认知转变中的关键丘脑核团。LGN中的丘脑皮层（thalamocortical, TC）细胞接受来自多个来源的调制性输入，包括丘脑网状核、视觉皮层和膝状体内中间神经元，但其放电的主要驱动来自少数（1-6个）突触连接的视网膜神经节细胞（retinal ganglion cells, RGC）。视网膜-膝状体传递（retinogeniculate transmission）的效率已被大量实验和计算研究探讨，以帮助理解丘脑在处理感觉输入中的作用。

早期实验发现，arousal 水平的增加会增强视觉诱发的LGN活动。更详细的实验研究通过使用s-电位（s-potentials，代表视网膜兴奋性突触后电位的细胞外事件）或双单元记录，证明较短的视网膜峰间间隔（interspike intervals, ISIs）更成功地驱动LGN放电。然而，大多数具有这种精细细节的数据来自麻痹、麻醉动物，因此关于非麻醉状态下视网膜-膝状体传递仍存在许多未解问题，包括：不同arousal状态下控制传递的生理特征是什么，为什么在清醒与麻醉动物中较长视网膜ISI的传递概率不同，以及如何解释清醒动物中不同LGN细胞传递效率的异质性。

此外，在非麻醉动物中，其他丘脑功能与视网膜传递同时发生并可能相互影响，其核心是节律性振荡。在自然清醒状态（relaxed but alert）下，丘脑产生α节律。然而，该节律如何产生、如何与感觉信号及皮层节律性反馈相互作用，以及它是促进还是阻碍视网膜传递，仍不清楚。同时，丘脑可因感觉信号、神经调质信号或神经精神疾病导致的节律失常而偏离α节律，这种偏离对视网膜传递特性的影响也是开放问题。

二、研究内容与结论

研究人员构建了一个新型的基于电导的（Hodgkin-Huxley公式）LGN丘脑皮层细胞模型，该模型能够产生与实验观察到的arousal状态相关的5种不同体外放电模式。通过该模型，研究人员证明了丘脑如何响应神经调质信号在状态间转换，并研究了不同arousal状态下视觉信息如何通过丘脑被门控或传递。研究最关注的5种状态中与自然清醒行为状态相关的动态区域，即与清醒α节律重合的"高阈值爆发"（high-threshold bursting）放电类型。在这些动态区域内，研究表明丘脑具有过滤视网膜神经节细胞放电序列的灵活性，表现出优先传递较短视网膜ISI的倾向，与实验结果一致。模型表明，M电流和L型钙电流的加入对于产生这些观察到的体外动力学以及清醒状态下的视觉信号过滤至关重要。研究还表明，该模型能够预测丘脑解释可能来自V1视觉皮层或丘脑网状核的视网膜外节律性输入的能力。最终，研究提出了对清醒丘脑α节律的新理解：该节律作为丘脑处于工作状态的标志，相关的过滤和夹带能力由M电流和L型钙电流的相互作用决定。

三、主要技术方法

研究人员构建了包含M电流（I_M）和高阈值L型钙电流（I_CaL）的新型LGN丘脑皮层细胞计算模型。模型中的所有电流均在丘脑内表达，参数值范围基于电流钳记录并参考先前建模研究。通过固定离子电导值并增加外加电流（I_app），模型产生五种动态放电类型。通过改变钾漏电导（g_KL）和外加电流模拟神经调质效应。视网膜输入通过兴奋性AMPA突触连接到模型LGN TC细胞，视网膜放电序列的ISI由伽马分布生成。研究还使用周期性（5-70Hz）和嵌套（25-70Hz嵌入4-12Hz）的兴奋性和抑制性节律性输入刺激模型，分析其夹带特性。

四、研究结果

2.1 新型丘脑皮层细胞模型对于捕捉实验观察到的清醒态动力学是必要的

通过固定离子电导值并增加外加电流，模型依次产生：低阈值爆发（I和II）、慢速紧张性放电（III）、高阈值爆发（IV）和快速紧张性放电（V）五种动态类型。这些状态可与arousal水平相关联：状态I和II对应慢波睡眠和δ睡眠节律，状态III和IV对应"自然清醒"状态，状态V对应注意力处理相关的高兴奋状态。

2.2 神经调质arousal信号将丘脑细胞从低阈值爆发切换到清醒态动态状态

通过绘制钾漏电导与外加电流的参数空间图，研究发现模型存在6个动态不同区域（0-V），其中区域II和III是双稳态的。谷氨酸能和胆碱能信号通过降低背景钾电导实现睡眠到清醒的转换。研究还发现，仅通过调制g_KL无法实现全部状态转换，达到慢波睡眠的<1Hz频率需要额外的外加电流变化。

2.3 LGN对视网膜信号的传递取决于模型arousal水平

在状态I和II中，视网膜放电不塑造输出模式，而是作为整体兴奋增加低阈值爆发频率。在自然清醒状态III和IV中，模型进入"过滤模式"，仅传递部分视网膜放电。在状态V中，视网膜放电几乎被忠实传递，与放电时序无关。

2.4 模型的新型清醒态动态状态重现实验观察到的平均视网膜传递变异性

传递频率和延迟时间统计表明，过滤主要限制在状态III和IV。模型产生的传递比分布与实验观察一致，大部分过滤参数空间的传递频率在0.25-0.75之间，高传递比（>0.75）仅在少数参数空间出现，表明高传递比应罕见。

2.5 高阈值爆发状态III和IV优先传递短视网膜ISI，与实验观察到的清醒动物视网膜传递相似

模型延迟时间分布显示，大多数（>90%）LGN放电发生在最近视网膜放电后0-5ms内，与实验一致。视网膜ISI分析表明，短ISI（0-10ms）几乎保证产生LGN输出，15ms后传递概率急剧下降，而长ISI（>30ms）在清醒状态下也有较高传递概率，这与麻醉动物中长ISI传递概率接近零不同。

2.6 不同视网膜ISI的过滤取决于I_M和I_CaL提供的时间尺度

M电流是产生状态III和IV的必要条件，其存在使过滤成为可能。I_CaL的失活时间常数和I_M的衰减时间是决定短和中等长度视网膜ISI传递概率的关键因素。增加g_M可降低中等ISI传递概率20-50%，而降低g_CaL可使短ISI传递概率降低约40%。

2.7 丘脑夹带至α和β周期性频率，但θ和γ作为嵌套振荡传递最佳

模型优先传递α和β兴奋性周期性频率（7-30Hz完美夹带），γ频率（30-70Hz）不被忠实传递，θ频率（4-7Hz）产生部分2:1夹带。兴奋性γ嵌套于θ频率振荡时传递最成功，在4-5.5Hz和7.5-9Hz有两个完美1:1夹带峰。抑制性周期性输入在5-28Hz产生三种夹带形式，但输出限制在8-16Hz。抑制性γ嵌套于慢振荡时，丘脑夹带至慢频率而非快频率。状态IV对于夹带较慢周期性兴奋性节律是最优的，从状态III到IV，夹带上限从约15Hz增加到25Hz以上。

五、讨论与结论

研究结果表明，清醒态α节律是丘脑处于工作状态的标志，而非直接因果机制。过滤和夹带能力由M电流和L型钙电流的相互作用决定，这两种电流恰好能在无输入时产生α节律。模型预测，丘脑可能传递较慢的周期性频率，而较快的γ节律只有在嵌入较慢频率振荡时才能成功传递。研究还提出，精神分裂症等神经疾病对丘脑α节律工作状态的扰动可能导致丘脑局部和全脑水平的显著处理错误。该模型为理解不同arousal状态下的节律动力学、丘脑对视网膜-膝状体传递的控制以及神经疾病对丘脑处理的影响提供了可测试的预测。