脑桥苔藓纤维输入至小脑皮层的高维双向群体编码

《Cell Reports》：High-dimensional bidirectional population coding in pontine mossy fiber inputs to the cerebellar cortex

　　

摘要：小脑通过脑桥核收集来自新皮层的感知运动信息，对其进行混合和扩展，以形成用于预测运动及其感觉后果的关联。然而，这一主要脑桥输入至小脑皮层的功能特性在很大程度上仍未被探索。研究人员利用三维双光子钙成像技术，发现脑桥苔藓纤维采用了一种高维双向群体编码。然而，仅有一小部分神经活动维度需要用于编码特定的运动和小须运动相关变量。值得注意的是，具有相反响应特性的轴突在颗粒细胞层内相互交织。受实验约束的回路模型表明，颗粒细胞接受不同组合的相反符号输入，能够有效地将高维双向群体编码传播至下游回路。研究结果确立了脑桥输入至小脑皮层的群体水平特性，并表明其神经表征为小脑扩增重编码提供了出乎意料的高维底物。

论文解读

研究背景与科学问题：小脑在运动控制与感觉预测中发挥关键作用，其通过接收来自大脑皮层运动指令的副本（内反馈/伴随放电）形成内部预测模型，从而协调动作并区分自生与外部感觉输入。感觉运动皮层-小脑回路对于适应环境变化至关重要，但小脑内进行的感运动计算本质仍不明确，主要原因在于对输入至小脑皮层的突触信号特性及其时间动态缺乏了解。苔藓纤维轴突（MFA）是向小脑皮层输入层传递感觉和运动信息的主要途径，其中最大的来源为基底脑桥核（BPN）。尽管已有单纤维记录显示个别MFA可编码肢体运动或感觉刺激，但缺乏群体水平的记录，使得无法理解复杂感觉运动信息在输入层如何被表征。小脑理论中的“扩增重编码”假说认为，MFA群体活动被投射到更高维的颗粒细胞（GrC）活动空间，以增加线性维度、去相关并分离重叠的活动模式。然而，近期GrC群体记录发现，简单学习行为中表征维度较低，而复杂自发行为中维度较高。若无MFA输入的群体功能特性，则难以判断GrC的高维表征是源于经典扩增重编码，还是仅仅继承了MFA群体中高度结构化的感觉运动活动。为填补这一空白，研究人员利用声光透镜（AOL）三维双光子显微镜，在清醒行为小鼠的颗粒细胞层记录了鉴定来源的脑桥MFA群体活动。

研究内容与结论：研究人员在Slc17a7-IRES-Cre转基因小鼠的BPN内注射AAV9-CAG-Flex-GCaMP6f，特异标记脑桥MFA，并在Crus I/II区域植入慢性颅窗，结合高速视频追踪与深度学习算法量化自发行为（运动、前后肢轨迹、 whisking参数）。通过改进的轴突分组方法，平均每记录会话获得191±4个MFA。结果显示，MFA群体活动在自发行为中高度异质：52.3%的轴突在活跃状态（A-state）正调节（MFA+），27.0%负调节（MFA?），20.7%无显著调节（MFA?）。主成分分析（PCA）表明，A-state与安静清醒状态（Q-state）的活动流形在嵌入子空间内近似正交，且这种正交性不因排除强调节轴突而消失。归一化特征谱显示群体活动方差分散在多个模式中，共享维度分析估计每维度仅需2.4–2.9个MFA，即嵌入子空间利用了37%的神经活动状态空间，呈现高维编码。然而，交叉验证偏最小二乘（PLS）回归显示，仅需3–5个潜在变量即可解释所测运动与小须变量的主要方差，表明行为编码子空间维度远低于嵌入子空间。单个最信息性MFA可解释特定行为变量交叉验证解释方差（CVEV）的8%–37%，但信息分布在一小组轴突中，且MFA+比MFA?携带更多信息。空间分析发现，MFA+、MFA?与MFA?的膨体（boutons）在颗粒细胞层内相互交织，在数十至数百微米尺度上无空间聚类或相关性结构。基于实验约束的三维尖峰模型（混合模型）预测，GrC群体平均发放率低于MFA，但群体稀疏度更高，且GrC+、GrC?、GrC?比例得以保留；当考虑MFA至GrC的局部数值扩增（3.7倍）时，GrC嵌入子空间维度增加。相比之下，将相反符号输入空间隔离的集群模型虽改变了GrC比例平衡，但低维特征谱与有效维度仍与混合模型类似，提示双向MFA编码的传播对空间混合程度具有稳健性。

关键技术方法：1. 在Slc17a7-IRES-Cre转基因小鼠的脑桥核内立体定位注射Cre依赖性AAV表达GCaMP6f，实现对脑桥MFA的特异性标记（样本来自4只成年小鼠，雌雄兼有）。2. 使用具有实时三维运动校正的声光透镜（AOL）三维双光子显微镜，以8–20 Hz采样率记录250×250×100 μm体积内数百个MFA膨体的钙活动。3. 采用改进的相关性分析与线性偏差比（LDR）对同一轴突上的膨体进行分组，并推断MFA活动。4. 结合高速视频追踪（Lightning Pose）与深度学习量化前肢运动、小须角度、幅度和设定点等行为参数。5. 应用主成分分析（PCA）、交叉验证PCA、偏最小二乘回归（PLS）和岭回归分析群体活动的几何结构、维度和行为编码。6. 使用受实验约束的三维尖峰模型（NeuroML2/jLEMS模拟）模拟颗粒细胞层回路，包含MFA-膨体-颗粒细胞连接、短期突触可塑性及AMPA/NMDA/GABA受体电导。

研究结果分述：
（1）自发行为期间脑桥MFA群体活动异质性：通过三维钙成像记录鉴定出719条轴突（其中561个单膨体轴突），发现MFA活动 profiles 在运动和 whisking 期间高度多样，且正负调节轴突并存。
（2）相反的MFA响应产生异质性群体活动：通过二元状态分类（A-state vs Q-state），明确MFA+、MFA?和MFA?三个亚群；成对相关分析显示A-state期间相关性分布变宽，出现正负相关共存；空气 puff 刺激也诱发类似双向响应。
（3）群体活动的几何结构：PCA显示Q-state与A-state对应嵌入子空间中分离且近似正交的流形（流形间距显著大于内部距离，最大主角度约1.52弧度）；逐步排除强调节轴突后正交性仍保持，表明并非仅由MFA+和MFA?的简单分离造成。
（4）群体活动维度与观察行为的关系：特征谱显示方差分散，有效维度为12±2，显著高于7个行为参数的有效维度；共享维度外推估计100个轴突需要37个维度（约3轴突/维），说明MFA群体使用高维编码。
（5）嵌入子空间与行为变量的关系：PLS交叉验证显示，最佳分量数3–5个即可预测运动和小须变量，单个PLS分量对运动相关变量解释较多，而小须幅度和设定点需要更多维度；全部行为变量CVEV为11%–48%，远低于嵌入子空间维度。
（6）部分轴突对特定行为高度信息性：利用岭回归解码，每个会话中单一最信息性MFA可解释相应行为CVEV的8%–37%；信息在排名靠前的10个轴突中逐渐下降而非断崖式，说明信息分布在不同轴突亚群中。
（7）功能不同MFA在颗粒细胞层的空间分布：X和Y维度上MFB+与MFB?分布无显著差异，Ripley's L函数显示空间随机性；膨体间相关系数与距离无关，表明相反符号输入在颗粒细胞树突整合尺度（~40 μm）上充分混合。
（8）生物约束模型模拟颗粒细胞群体活动：混合模型在统计特性上与实验数据匹配，预测GrC平均发放率低于MFA但群体稀疏度更高；GrC+和GrC?比例受MFA+和MFA?局部占比影响；共享维度在MFA与GrC间类似，但计入MFA-膨体-颗粒细胞数值扩增后GrC嵌入子空间维度显著增加，支持扩增重编码。集群模型虽改变局部输入比例，但仍保留高维双向编码的传播能力，证明该机制对空间组织程度不敏感。

讨论与结论总结：研究首次在群体水平刻画了脑桥MFA输入至小脑皮层的编码特性。讨论部分指出：双向MFA群体编码（MFA+/MFA?）可能与新皮层ON/OFF编码策略类似，具有代谢优势和高效解码性能；高维嵌入子空间与低维行为编码子空间的分离表明存在大量未观测的感觉运动变量或内部状态；脑桥MFA的响应异质性可能继承自新皮层，并受到脑桥核本身三维结构、亚皮层输入及小脑核反馈的调制；相反符号轴突的空间交织为GrC提供了多样的输入组合，有助于产生双向GrC群体代码；模型预测与经典Marr-Albus理论一致，但MFA本身已具备高维嵌入子空间，实际扩增重编码的作用体现在MFA到GrC的数值扩增上（局部扩展约2–3倍维度）。局限性包括钙成像时间分辨率有限、仅记录脑桥来源MFA、模型缺乏动态反馈抑制等。

结论：研究人员结合群体钙成像、行为量化和生物约束模型，揭示了脑桥MFA输入采用高维双向群体编码，行为变量仅占用其中一小部分维度；相反符号轴突在颗粒细胞层内空间交织，并通过MFA-GrC数值扩增实现扩增重编码，产生GrC中观察到的双向高维活动。该研究确立了脑桥-小脑输入路径的群体编码原则，为理解小脑如何利用高维感觉运动信息进行预测性控制提供了基础。论文发表在《Cell Reports》。