发育期大鼠新皮层上行神经调制系统的定量解剖学与生物物理建模

《PLOS Computational Biology》：Quantitative anatomy and biophysical modeling of ascending neuromodulatory systems in the developing rat neocortex

　　

摘要：研究人员利用免疫细胞化学染色与体视学(stereology)技术，定量测定了两周龄Wistar大鼠体感觉皮层中胆碱能(cholinergic)、儿茶酚胺能(catecholaminergic)及5-羟色胺能(serotonergic)上行神经调制系统在新皮层神经毡(neuropil)内的每单位皮层体积纤维长度及曲张体(varicosity)密度，并提出一套基于上述数据生成生物物理详细(biophysically detailed)计算模型的方法学框架。所获数据被整合入体感觉皮层生物物理详细计算模型中，以探究神经调制支配的解剖学组织方式及其功能意义。研究人员发现，尽管神经调制支配稀疏，仍对网络活动产生显著影响。网络模拟支持乙酰胆碱(ACh)抑制慢振荡并促进皮层网络去同步化的假说，与已有文献广泛相符。此外，ACh调制的时序特征符合突触释放(synaptic release)而非容积传递(volume transmission, VT)。研究人员还发现，感觉皮层内多巴胺(DA)与5-羟色胺(5-HT/serotonin)的释放通过抑制δ振荡诱导网络去同步化，且5-HT还可诱发θ振荡的出现，揭示了二者调控皮层网络活动中此前未被充分探讨的方面。实验数据与生物物理计算模型均以开放获取形式提供作为社区资源。

论文解读：

本文发表于《PLOS Computational Biology》，研究背景如下：上行神经调制(neuromodulatory, NM)系统——主要包括乙酰胆碱(acetylcholine, ACh)、多巴胺(dopamine, DA)和5-羟色胺(5-hydroxytryptamine, 5-HT/serotonin)——通过基底前脑(basal forebrain, BF)、腹侧被盖区(ventral tegmental area, VTA)及中缝核(dorsal raphe, DR)等核团投射至新皮层，调控神经元通讯及脑状态（注意、睡眠、觉醒等）。ACh的去同步化作用已有较多记载，但其机制尚未明确；DA与5-HT在感觉皮层（而非前额叶）的微环路效应及多尺度影响数据匮乏。此外，上行调制系统究竟以经典点对点突触传递(synaptic transmission, ST)还是以非突触弥散性容积传递(volume transmission, VT)发挥作用长期存在争议，传统二维超微结构分析难以解决此问题。因此，研究人员通过在P14 Wistar大鼠初级体感觉皮层(primary somatosensory cortex, S1，后肢代表区hindlimb representation, HLS1)定量测定三种NM系统的纤维长度密度与曲张体(varicosity)密度，并将之整合入已有的新皮层微环路生物物理模型中，分别构建ST与VT模式的NM投射，模拟激活后观察其对网络慢振荡及去同步化的影响，以阐明解剖—功能关系并探讨ST与VT的作用差异。

主要关键技术方法：取出生后14天(P14) Wistar大鼠(n=11)，取S1 HLS1区进行胆碱乙酰转移酶(choline acetyltransferase, ChAT)、酪氨酸羟化酶(tyrosine hydroxylase, TH)及5-HT免疫染色，采用体视学Space Ball探针估算每单位体积纤维长度，光学分馏器(optical fractionator)估算曲张体密度，并进行组织收缩校正。将实验测得的层状(laminar)曲张体密度分布导入基于Markram等建立的大鼠S1新皮层微环路模型（含163,528个细胞），使用Projectionizer工具按层状密度与树突几何分布放置突触释放位点(synaptic release site, sRS)，虚拟纤维垂直贯穿各层；VT模式以释放位点为球心设影响半径R_max（ACh: 5 μm；DA: 2 μm；5-HT: 3 μm），在球内树突段随机添加随距离衰减的调制电导。NM效应以受体介导的附加电导(phenomenological conductance-based implementation)实现，依据文献报道的细胞类型特异性去极化/超极化方向设置。网络模拟采用NEURON环境下的Neurodamus，施予~20 Hz训练刺激，分析功率谱密度(power spectral density, PSD)及δ带(1.5–3 Hz)功率变化，重复30次取均值。

研究结果如下：

Cholinergic fibers（胆碱能纤维）：ChAT免疫染色显示各层均存在曲张纤维，纤维长度密度与曲张体密度在皮层厚度方向较均匀，仅第Ⅰ层与第Ⅳ层间有显著差异(p=0.048)，其余各层无显著差异。

Catecholaminergic fibers（儿茶酚胺能纤维/主要为DA）：TH免疫反应(TH-ir)纤维密度在第Ⅰ、Ⅳ层最高，第Ⅴ、Ⅵ层最低，但层间差异未达显著；曲张体呈相似层状分布趋势亦无统计差异。

Serotonergic fibers（5-羟色胺能纤维）：5-HT-ir纤维在上颗粒层(第Ⅰ、Ⅱ层)密度显著高于第Ⅴ层(p=0.007, 0.035)，第Ⅰ与第Ⅳ层有显著差异(p=0.048)；曲张体密度也以浅层较高。

Comparison of the three ascending NM systems（三种上行NM系统比较）：按单位体积纤维长度排序：胆碱能 > 儿茶酚胺能 > 5-羟色胺能（全层及各分层均显著）。胆碱能曲张体密度为儿茶酚胺能的1.4倍、为5-HT的2.3倍；儿茶酚胺能曲张体密度为5-HT的1.6倍。

In silico predictions about the organization of NM input（NM输入的计算机预测）：基于实测曲张体密度与重建神经元树突几何，预测单根NM纤维所支配的突触后神经元数目及接触数。所有系统中TPC（顶树突锥体细胞）和LBC（大篮状细胞）因树突总长最大而接受最多NM接触。胆碱能系统所及突触后靶细胞总数最多。DA系统影响最广，可支配各层兴奋性与抑制性神经元。

Simulating cholinergic network effects（ACh网络效应模拟）：ACh ST激活使δ带功率降低68%(p=6.07×10^?36)，实现持续去同步化；VT模式(R_max=5 μm)在仅10%纤维招募时δ功率降78%，90%纤维时降99%(p=1.85×10^?42)。与实验optogenetic数据更匹配的是ACh ST或弱VT，提示ACh在体作用更符合突触传递或受限容积传递特征。

Simulating dopaminergic network effects（DA网络效应模拟）：DA ST使δ功率降约49%(p=3.37×10^?32)；VT模式下10%纤维招募时δ降50%，90%时降68%(p=9.1×10^?38)。DA同样具去同步化效应，VT空间覆盖范围使其影响略强，但效果不及ACh持久。

Simulating serotonergic network effects（5-HT网络效应模拟）：5-HT ST使δ功率降约87%(p=7.5×10^?41)，并将功率谱向高频(5–8 Hz, θ频段)偏移；VT模式下10%纤维招募时δ降99%，同时诱发低θ频段振荡。表明5-HT除抑制慢振荡外可诱导θ节律出现。

讨论与结论翻译总结：研究人员定量确认新生大鼠S1中胆碱能纤维密度最高，5-HT纤维偏好浅层分布。计算模型预测胆碱能系统支配靶细胞最多，DA系统支配范围最广。网络模拟表明三种NM系统均可抑制慢(δ)振荡并促去同步化，其中ACh去同步化最强，5-HT独特地诱发θ振荡；比较模拟结果与实验数据提示皮层ACh作用更符合突触传递或弱VT而非广域VT。ST与VT在时空特征上存在差异——VT具更大空间足迹(spatial footprint)与较慢动力学，少量VT纤维即可产生近似全量ST纤维的网络效应。该框架将定量解剖数据与生物物理模型结合，为解析上行NM系统解剖—功能关系提供可检验假设与开放资源，未来可纳入突触前调制(presynaptic modulation)及受体亚型特异机制进一步完善。