神经元活动增加时，局部代谢需求通常上升，血流、血容量和血氧状态随之改变。fMRI正是通过这些血流动力学变化，间接推断脑活动。这个过程被称为神经血管耦联（neurovascular coupling）。
间接测量带来一个长期存在的解释难题：BOLD升高，究竟对应多少神经元在兴奋？是否也可能同时存在活动下降的细胞？不同细胞类型、血管结构和局部回路，会不会让相同强度的神经活动产生不同的BOLD信号？
传统钙成像能以细胞分辨率追踪神经元活动，却通常只能覆盖有限的皮层区域；fMRI能覆盖全脑，却难以识别信号来自哪些细胞。过去也有研究把光纤记录与fMRI结合，但光纤收集的是一团混合信号，难以区分单个神经元、神经网状结构（neuropil）和非神经细胞。
要回答上述问题，关键不是把两类数据分别采集后再比较，而是在同一只动物、同一次刺激、同一段时间里同步记录。
强磁场、射频脉冲和快速切换的梯度，对普通显微镜并不友好。金属、电子元件和常见光学材料可能引起磁场畸变；MRI脉冲甚至会让相机产生微米级位移。
研究团队因此设计了一套MRI兼容单光子显微镜（MRI-compatible single-photon microscope），将不兼容的光源和电子部件移到磁场影响较小的位置，并通过光纤传递激发光。
镜头材料也经过重新选择。常见BK7玻璃和梯度折射率镜头会造成明显磁敏感伪影，研究人员改用丙烯酸非球面镜片制作小型物镜，并用重水（D₂O）作为浸润介质，减少空气与组织界面造成的信号空洞。测试显示，相机距离磁场中心超过6厘米后，射频诱发的位移基本消失；最终光路将这一距离延长到7.1厘米。
这套系统在小鼠桶状皮层（barrel cortex）获得了1.45平方毫米视野，像素分辨率1.1微米，实测光学分辨率为2.5±0.24微米，同时没有观察到明显的BOLD伪影或信噪比下降。换句话说，它把“局部单细胞视角”和“全脑血流动力学视角”稳定地放到了同一时间轴上。
实验使用表达钙指示蛋白GCaMP6f的小鼠，记录桶状皮层第2/3层神经元。为减少麻醉对神经血管耦联的干扰，小鼠在清醒状态下接受成像，并在实验前经历至少10天的适应训练。
每只小鼠接受12轮刺激，每轮用3赫兹气流刺激对侧胡须10秒。全脑fMRI显示，丘脑、纹状体、初级和次级躯体感觉区等多处出现显著BOLD反应。
值得注意的是，听觉皮层也出现强烈激活。论文认为，这可能来自气流装置产生的可听声音，而不只是胡须感觉本身。这个细节很有启发性：全脑成像捕捉到的信号，不仅反映研究人员想研究的过程，也会记录实验装置、运动和应激带来的额外输入。对任何fMRI研究而言，控制条件并不是技术附属品，而是解释结论的一部分。
更值得关注的是四类信号的时间顺序。
神经网状结构的钙信号变化约为2.7%，达到80%峰值只需0.2秒；单细胞群体平均钙信号变化约为1.7%，对应时间为1.1秒；血管直径增加约3.8%，时间为1.4秒；BOLD信号最终变化约4.9%，时间为1.7秒。
这个顺序符合神经血管耦联的基本逻辑：神经活动先发生，血管随后扩张，BOLD再以更慢的速度出现。
但这些数字也提醒我们，BOLD不是神经活动的即时复制。它经过了血管反应和血氧变化的“加工”，时间上更慢，空间上也更容易混合多个来源。
神经网状结构和单细胞群体的峰值时间也并不相同。前者包含大量轴突、树突和突触活动，后者是被识别细胞的平均结果。两条曲线之间的差异说明，即使在光学成像内部，“总体活动”也取决于测量对象和汇总方式。
研究人员最终从5只小鼠中提取了236个第2/3层神经元，平均每只约46个。随后，他们用支持向量机（support vector machine, SVM）根据单个神经元的钙活动预测局部BOLD变化。
模型的归一化均方根误差显著低于时间打乱后的基线，说明局部神经元群体确实包含可用于解释BOLD的信息。
但最有意思的结果出现在“空间位置”上。
研究团队标注显微视野中的血管，计算每个神经元到最近血管的欧氏距离，再用递归特征消除（recursive feature elimination, RFE）评估哪些神经元对BOLD预测最重要。
对胡须刺激反应增强、在模型中获得正权重的神经元，与血管的距离分布并没有明显特殊性，和总体分布相比P值为0.8919。
相反，那些在刺激期间活动下降、对BOLD预测呈负权重的神经元，更靠近血管，P值为0.0440。
这意味着，在一个整体BOLD升高的区域内，并非所有细胞都同步兴奋。某些靠近血管的神经元活动下降，局部血管和BOLD信号却仍在上升。
BOLD更像是神经活动、细胞类型、局部回路和血管结构共同加权后的输出，而不是一张简单的“神经元活跃度地图”。
论文提出了一种可能解释：这些神经元或许处在与血管活性中间神经元相关的局部回路中。某些抑制性神经元可通过释放血管活性物质调节血流，例如血管活性肠肽阳性（VIP⁺）细胞倾向于促进血管扩张，而生长抑素阳性（SOM⁺）细胞可能参与血管收缩。
因此，神经元活动下降与血管扩张同时出现，并不一定矛盾。局部回路中可能存在“神经活动方向”和“血管反应方向”不完全一致的组合。
不过，这仍是机制假设，而不是该研究已经直接证明的结论。当前成像主要针对第2/3层锥体神经元，并未在同一实验中同时分辨多种抑制性细胞亚型。
后续需要多色钙成像（multicolor calcium imaging）、细胞类型特异性标记或光遗传学（optogenetics）干预，才能确认哪些细胞在驱动这种反向关系。
研究团队又把问题反过来：既然一个局部回路受到全脑网络输入，那么远处脑区的BOLD信号，能否预测桶状皮层的神经元群体活动？
结果显示，不仅局部桶状皮层，其他参与胡须感觉处理的皮层和皮层下区域，其BOLD信号也能显著优于时间打乱基线地预测局部钙活动。分析覆盖躯体感觉区、运动区、听觉区、丘脑和纹状体等区域；各区域模型与打乱基线相比的差异达到P≤0.01。刺激同侧与对侧脑区的预测能力则没有显著差异，P值为0.18。
这并不表示远处脑区“决定”了局部神经活动。预测关系不能自动等同于因果关系。
更合理的理解是：局部桶状皮层并非孤立工作，它嵌入一个分布式感觉—运动网络中。多个脑区的BOLD变化共享了与该局部回路相关的信息。fMRI看到的全脑图谱，可能包含局部细胞活动在网络中的投影。
它没有证明人类fMRI过去的解释都是错误的，也没有建立从BOLD精确还原单个神经元活动的通用公式。
实验规模仍然有限：共使用12只小鼠，其中3只因神经活动不足、MRI信噪比较低或植入区信号缺失而被排除；BOLD分析主要基于9只小鼠，单细胞分析只有5只。靠近血管的负响应神经元差异P值为0.0440，而且该比较没有进行多重检验校正，因此更适合被视为需要独立重复的线索。
此外，单光子显微镜难以解析深度，神经元—血管距离是在二维视野中估计的；血管标注主要依据可见的大型软脑膜血管；显微镜视野与MRI体素还没有实现精确共配准；实验采用的也是相对简单的胡须刺激。小鼠桶状皮层的结果不能直接外推到人类复杂认知任务或疾病状态。
但这项研究提供了一个重要的观念修正：解释BOLD时，不能只问“这个脑区是否激活”，还应进一步追问：信号由哪些细胞群贡献？这些细胞与血管如何排列？局部回路中是否同时存在兴奋和抑制？远处网络又提供了多少共享信息？
未来，当多色成像、双光子显微镜与fMRI进一步整合，研究人员可能同时区分兴奋性神经元、抑制性神经元和胶质细胞，并在更深层、更复杂行为中建立细胞活动与全脑信号的映射。
到那时，fMRI图上的一个亮点，或许不再只是“哪里更活跃”，而会逐渐变成一个更具体的问题：究竟是哪一类细胞、通过怎样的血管与网络机制，让这里变亮？
参考文献
