导语： 神经血管耦合将大脑活动与局部血流变化联系起来，构成了非侵入性脑功能成像的基础。然而，血管反应的空间特异性如何与神经元活动模式相对应，仍不完全清楚。2024年9月4日，《自然·神经科学》在线发表了一项来自加拿大蒙特利尔大学等机构的重要研究（2024年第27卷2101–2114页）。利用多尺度成像技术，研究者探究了由单根胡须刺激在小鼠皮层不同柱和层中诱发的血管活动如何与神经元活动在空间上相关联。他们发现，介观血流动力学信号在空间上定量反映神经元活动，但由单个血管节段上高度异质的反应模式组成，且这种异质性无法被局部神经元活动所预测。相反，这种异质性依赖于血管方向性，特别是在丘脑皮质输入的第4层，毛细血管沿其下游灌注域优先响应神经元活动模式。因此，毛细血管基于远端活动精细调节血流，并编码层特异性活动模式。这些发现意味着血管解剖结构为功能成像信号设定了分辨率极限：单个血管不精确报告其紧邻区域的神经元活动，而是整合沿血管树的活动模式。

研究背景：神经血管耦合的空间特异性之谜

功能成像的基础与困惑

• 神经血管耦合是fMRI、光学成像等技术的生理基础：神经元活动增加 → 局部血流增加（血流动力学反应）

• 核心假设：血流变化精确反映局部神经元活动

• 挑战：血流反应的空间特异性与神经元活动模式之间的关系尚未在单血管水平系统量化

已有矛盾的提示

• 远隔血管反应：刺激可诱发远离激活神经元区域的血管扩张

• 静息态相关性弱：血流与神经元活动的相关性在静息态较弱

• 血管网络结构：皮层血管网络是高度互连的网状结构，而非独立的柱状单元

本研究的问题

• 单个血管节段的扩张是否精确反映其邻近神经元的活动选择性？

• 如果否，什么因素决定了血管反应的空间模式？

核心发现之一：血管反应选择性差于局部神经元活动

实验范式（图1）

• 模型：小鼠桶状皮层（S1），单根胡须（C2或D2）刺激

• 技术：同时记录神经元钙信号（jRGECO1a或GCaMP6f）和血管直径（荧光染料标记血浆）

• 分析：计算选择性指数（SI）：(R_associated - R_neighbor) / (R_associated + R_neighbor)

关键结果

• L2/3（图1b-c）：

  • 局部神经元：对关联胡须高度选择性

  • 血管节段：扩张反应选择性差，许多血管对非关联胡须同样反应

• L4（丘脑输入层，图1d-f）：

  • 局部神经元：高度选择性（比L2/3更强，Extended Data Fig. 2）

  • 血管节段：仍缺乏选择性

结论：单个血管的扩张反应不能精确反映其邻近神经元的活动选择性。

核心发现之二：血管反应选择性依赖于血管方向性，而非分支等级

分析框架（图2-3）

• 将血管按分支等级（从穿通动脉起1-8级）和方向（相对于桶状柱的位置）分类

• 使用广义线性模型评估影响选择性的因素

关键结果

• 分支等级无关（Extended Data Fig. 4-5）：

  • 在L2/3和L4中，不同分支等级的毛细血管选择性无显著差异

  • 动脉和小动脉选择性略高，但仍远低于神经元

• 方向性相关（图3e-f）：

  • 沿C2-D2轴的血管位置与选择性显著相关

  • 朝向关联桶状柱中心的血管选择性更高

• 层特异性（图4）：

  • L4中，方向性对选择性的影响更强

  • L4血管更严格地沿灌注域（下游区域）反应

核心发现之三：毛细血管沿下游灌注域整合远端活动（图5-6）

血管网络追踪

• 标记单个穿通动脉及其整个毛细血管下游树

• 沿血管树追踪从动脉到静脉的直径变化

关键结果

• L4特异性模式（图5c-e）：

  • 毛细血管沿其下游灌注方向逐渐改变选择性

  • 例如：靠近C2桶的毛细血管段对C2选择性高，远离C2桶的同一毛细血管段对D2选择性升高

• 不是分支等级，而是“位置”：

  • 同一毛细血管的不同节段（属同一分支等级）选择性不同

  • 取决于该节段在网络中的位置（距离关联桶状柱的远近）

• L2/3中模式较弱（图5f）：

  • L2/3血管网络更各向同性，方向性影响较小

机制解释：血流动力学整合

• 毛细血管直径变化不仅影响局部血流，还通过影响下游阻力改变远端血流

• 毛细血管对远端活动的反应可能是为了优化整体灌注，而非仅服务局部神经元

核心发现之四：血流速度/流量的选择性也受方向性调控（图3g-h, Extended Data Fig. 6）

• 测量红细胞速度和流量（线扫描）

• 速度/流量选择性：同样受血管方向性影响，而非分支等级

• L4中更明显：毛细血管更严格地沿灌注域调节流量

核心发现之五：在清醒动物中验证（图7, Extended Data Fig. 9-10）

清醒 vs. 麻醉

• 麻醉（右美托咪定）：神经元活动空间特异性更高（可能因抑制性输入减少？）

• 清醒：神经元活动和血管反应空间特异性降低（Extended Data Fig. 10）

• 关键发现（图7d-f, Extended Data Fig. 9）：

  • 在清醒动物中，血管反应的选择性随时间变化

  • 早期相（刺激后0.5-1.5秒）：选择性差（类似麻醉下的整体反应）

  • 晚期相（刺激后3-4秒）：选择性提高（但仍低于神经元）

  • 提示：早期反应可能由快速传导的血管信号（通过内皮/周细胞网络）介导，晚期反应反映局部代谢需求

综合模型：血管解剖结构设定的成像分辨率极限

核心概念

1. 血管是“整合器”，而非“局部传感器”：

   • 单个血管节段的扩张不反映其紧邻区域的神经元活动

   • 而是反映沿其下游灌注域的整合活动模式

2. 层特异性：

   • L4（输入层）：血管网络更各向异性，方向性影响强 → 毛细血管编码层特异性活动模式

   • L2/3（层内处理）：血管网络更各向同性，整合范围更广

3. 对功能成像的意义：

   • 分辨率极限由血管网络结构设定，而非仅由成像硬件决定

   • 单个血管（或fMRI体素）的信号反映的是血管树下游区域的整合活动，而非其正下方的神经元活动

   • 解释为何fMRI可检测到远离激活脑区的血流变化（如默认模式网络负激活）

对功能成像和神经血管耦合研究的启示

对fMRI数据分析的实践建议

• 体素信号解释：某一体素的BOLD信号可能反映下游静脉引流区域的整合活动，而非体素内的神经元活动

• 层特异性fMRI：L4的血管方向性更强，可能解释为何层特异性fMRI能检测到丘脑输入（即使L4体素也包含来自其他层的血管）

• 静息态网络：血管的远端整合特性可能贡献于静息态功能连接（即使无直接神经元连接，共享血管下游区域也会产生相关血流）

资源可及性

• 数据集：https://doi.org/10.6084/m9.figshare.26121076

• 标准化ROI掩膜（S1亚区和单个桶状柱）包含在数据集中

• 分析代码：可向通讯作者索取

局限性与未来方向

1. 胡须刺激范式：仅测试了单根胡须刺激，未测试更复杂的自然感觉输入（如多根胡须、纹理探索）

2. 麻醉影响：主要实验在右美托咪定镇静下进行，清醒实验中神经元选择性和血管反应选择性均降低，提示麻醉可能放大了层特异性差异

3. 血管直径测量：仅测量了外径（而非内径），但外径变化与内径变化高度相关（因血管壁体积相对恒定）

4. 缺乏分子机制：未探究介导远端整合的细胞机制（内皮传导？周细胞网络？）

5. L4特异性机制未解：为何L4血管网络更具方向性？是否因L4毛细血管更密集或周细胞覆盖不同？

未来方向：

• 机制探究：使用内皮特异性敲除（如Cx40/Cx37）或周细胞特异性操控，测试内皮传导或周细胞在远端整合中的作用

• 自然刺激：使用多根胡须刺激或主动触觉任务，测试血管反应是否编码复杂感觉特征

• fMRI验证：在超高场fMRI（11.7T或更高）中，同时测量层特异性BOLD和血管直径，验证血管方向性是否预测BOLD信号的“错位”

• 人类验证：使用7T fMRI和血管成像，测试人类皮层中是否存在类似的方向性依赖血管反应

• 计算模型：构建血管网络-神经元活动耦合模型，预测给定血管树的血流反应，为fMRI信号解释提供生物物理基础

结语

这项研究挑战了神经血管耦合的“局部”假说，揭示了血管反应的空间特异性由血管解剖结构（特别是方向性）决定，而非简单地反映局部神经元活动。毛细血管沿其下游灌注域整合远端活动模式，在L4（丘脑输入层）中尤其显著。这意味着功能成像信号（fMRI、光学成像）的分辨率受限于血管网络结构，单个血管或体素反映的是沿血管树整合的活动模式，而非其紧邻区域的神经元活动。该发现为解释fMRI信号的空间模式（如层特异性、远隔反应、静息态网络）提供了新的生物物理框架，并提示在分析高分辨率功能成像数据时需考虑血管方向性和网络位置。

原始论文：Martineau, É., Malescot, A., Elmkinssi, N. et al. Distal activity patterns shape the spatial specificity of neurovascular coupling. Nat Neurosci 27, 2101–2114 (2024). https://doi.org/10.1038/s41593-024-01756-7

数据：https://doi.org/10.6084/m9.figshare.26121076