Nature Methods | 单细胞钙成像与全脑fMRI，能否在清醒小鼠上同步完成？|fmri|体素|信号|小家鼠|显微镜|细胞|血管

认知神经科学前沿文献分享

基本信息

Title:Simultaneous single-cell calcium imaging of neuronal population activity and brain-wide BOLD fMRI

发表时间:2026-07-15

发表期刊:Nature Methods

获取原文:

1. 添加小助手:PSY-Brain-Frontier即可获取PDF版本

引言

BOLD fMRI与单细胞分辨率光学成像之间存在长期难以跨越的方法学鸿沟。BOLD fMRI能无创记录全脑血流动力学响应，但只能提供局部神经活动的慢速、间接均值估计；单光子或双光子成像可在单细胞层面记录Ca²⁺或电压信号，却需要侵入式手术，且只能覆盖少量皮层区域（Fig. 1a–j）。两种模态各有取舍，都只能回答脑功能的一个尺度。

过去十余年里，研究者尝试用fiber photometry把Ca²⁺成像与BOLD fMRI缝合在同一只麻醉或清醒动物上。Schulz等（2012）在麻醉大鼠上做了开创性联合记录，指出BOLD信号不仅反映神经元源，还与星形胶质细胞的慢振荡Ca²⁺波相关；Lake等（2020）用多束光纤把全皮层Ca²⁺信号与BOLD并行采集，证明BOLD可在高空间精度下由局部神经活动部分重建。然而，fiber-based方法无法分辨单细胞，把神经元、神经胶质与neuropil信号混在一起，丢失了细胞类型特异性信息（Fig. 1a–j）。电生理联合fMRI虽然能拿到单细胞分辨率，但电磁屏蔽工程复杂且代价高昂。

本文的目标正是填补这一空白：组装一台MRI兼容的单光子显微镜，能在清醒行为小鼠的7T扫描仪内同步采集单细胞Ca²⁺活动与全脑BOLD信号，并把遗传靶向、低成本、易构建三个性质放在同一套设备上。作者随后做了两个概念验证——一个聚焦神经元与血管的空间关系如何决定局部BOLD预测的符号方向，另一个聚焦脑内多个whisker处理相关的远端区域如何协同预测SSp-bfd的群体活动。

实验设计与方法逻辑

整体实验沿两条线推进：先做显微镜硬件与成像范式的开发，再做两个分析层面的应用（Fig. 1a–j）。硬件部分需要解决三件事——消除MRI磁场对显微镜的干扰、保持BOLD信号的高SNR、实现单细胞分辨率的荧光成像。作者用acrylic aspheric lenses（Edmund Optics 15-271/36-629）替代传统BK7或GRIN物镜（Fig. 1c），用HDT180树脂（Loctite 3D 3860）做3D打印外壳（Extended Data Fig. 1），并通过两片消色差透镜把相机与等中心距离拉大到71.1 mm以避开RF激发脉冲引发的微米级位移（Fig. 1b,d）。CCD相机选用MRC HighResolution（MRC Systems GmbH），激发光由488 nm激光器通过400 μm光纤从扫描仪外耦合进显微镜。成像窗口为3 mm直径的cranial window，置于SSp-bfd（AP=2.0 mm, LR=3.0 mm）。样品与物镜之间用D2O（99.9% D-content）作浸没介质，以压制气-组织界面造成的磁敏感伪影（Fig. 1f inset）。最终成像FOV为1.45 mm²，像素分辨率1.1 μm，NA 0.363，光学分辨率2.5 μm FWHM（Fig. 1h–j）。

动物与手术：12只4–6月龄Rasgrf2-2A-dCre × Ai148D小鼠在L2/3锥体神经元表达GCaMP6f（Extended Data Fig. 2），其中2只因无神经活动或MRI SNR过低被排除，1只因植入体周围信号严重丢失被排除，最终保留9只用于BOLD分析、5只用于单细胞分析。术后进行≥10天的MRI环境习惯化，以framewise displacement < 0.05 mm为达标阈值（Extended Data Fig. 3）。fMRI扫描同时以10 Hz帧率采集单细胞Ca²⁺活动。

数据采集与分析：MRI使用7T BioSpec 70/20，GE-EPI序列（TE/TR=15/1500 ms，7片冠状面，体素0.2×0.2×0.8 mm）。Ca²⁺成像经TurboReg配准后用CNMF-E做细胞提取，得到236个L2/3 SSp-bfd神经元（47±14 per animal, N=5）。BOLD信号用FSL做GLM分析，显著激活阈值为z-score ≥ 3.1, FWE P ≤ 0.00001。SVM解码用linear epsilon-insensitive regression，NRMSE作为性能指标；通过SVM-RFE按三种权重准则（最低绝对权重、最低权重、最高权重）对神经元排序，量化单细胞对局部BOLD的预测贡献。血管距离定义为神经元中心到手动标注血管结构的最小Euclidean距离（2D，受单光子成像无深度分辨能力所限）。

核心发现

发现一：MRI兼容单光子显微镜在保持BOLD SNR的前提下实现单细胞分辨率

核心挑战不是单纯把显微镜放进MRI腔，而是要在不破坏BOLD信号质量的前提下实现细胞分辨率的荧光成像。作者用acrylic aspheric lenses替代常规BK7/GRIN物镜后，ex vivo脑组织BOLD EPI几乎看不到信号丢失（Fig. 1c, 右）；而BK7和GRIN物镜引入显著的磁敏感伪影。这一替换是整套装置得以成立的基础。

Figure 1. MRI-compatible microscope design for combined BOLD fMRI and calcium imaging

第二项关键设计是相机-磁场距离。MRI兼容CCD相机对RF脉冲仍存在微米级磁响应，作者定量测试后确定相机边缘与等中心距离d > 6 cm可避免位移（Fig. 1d），并通过两片消色差透镜把成像路径延长到71.1 mm，使最终成像稳定（Fig. 1b）。第三项是材料筛选——HDT180树脂在打印分辨率、耐久性与热稳定性间提供了最优平衡，SNR与含黑色颜料的聚合物相当（Extended Data Fig. 1, Extended Data Table 1）。三管齐下，最终FOV 1.45 mm²、像素1.1 μm、光学分辨率2.5 ± 0.24 μm FWHM（Fig. 1h–j），既不损失BOLD SNR，也不损失细胞分辨率。

这一步的工程意义在于，过去fMRI-光学联合记录只能在fiber photometry下做，把神经-血管问题压缩为群体均值；本文的装置首次把研究单元下沉到单细胞水平，并且保留遗传靶向能力（Ca²⁺指示剂由Cre-line控制）。Extended Data Fig. 2的组织学验证确认GCaMP6f确实表达在L2/3锥体神经元（neurogranin阳性）中，而非GABA能中间神经元（Extended Data Fig. 2c,d）。

Figure 2. Simultaneous BOLD fMRI and calcium imaging of neural population activity at single-cell resolution

发现二：SSp-bfd局部BOLD信号可由单细胞群体Ca²⁺活动解码，且神经元-血管距离决定解码权重的符号

概念验证一聚焦局部尺度。清醒小鼠在3 Hz contralateral whisker刺激下（12 trials，10 s pneumonic stimulation, 40 ms duty cycle），9只动物（6雄3雌）的BOLD fMRI在SSp-bfd、SSs、SSp-n、TH、STR、SSp-ul、MOp、MOs、AUDp、AUDd等区域出现显著激活（Fig. 2a；z-score ≥ 3.1, FWE P ≤ 0.00001），与已有whisker fMRI图谱一致。BOLD信号在SSp-bfd的trial-averaged响应达到~4.9%信号变化，TTP80% 1.7 s（Fig. 2c）；同步Ca²⁺成像上neuropil变化~2.7%（TTP80% 0.2 s），L2/3群体活动~1.7%（TTP80% 1.1 s），量级与时程与既有报道吻合（Fig. 2b,c,e）。

用linear SVM把236个L2/3神经元的ΔF/F Ca²⁺ trace解码为同区SSp-bfd平均ΔBOLD信号，NRMSE显著低于时序打乱基线（two-sided Wilcoxon signed rank test, P ≤ 0.005；Fig. 3a,b）。这一步本身只证明群体活动含有局部BOLD信息；进一步的发现来自SVM-RFE权重排序。

Figure 3. Spatial location of neurons with respect to the vasculature has an effect on the predictive relationship between BOLD and single-cell calcium activity

当用"最低绝对权重"准则排序时，神经元的血管距离分布与整体没有差异（Fig. 3d）；但改用"最高权重（保留负权重）"准则后，排名靠前的神经元到血管的距离显著小于总体（two-sided Wilcoxon rank-sum test, P = 0.0440；Fig. 3d,e）。也就是说，离血管近的神经元在whisker刺激下表现出负向ΔF/F Ca²⁺响应（Fig. 3e,f），而远离血管的神经元保留正向响应。Extended Data Fig. 7使用deconvolved ΔF/F（constrained FOOPSI, decay_time=0.4）后结论保持稳定——血管距离效应不依赖于ΔF/F的具体处理方式。

Figure 7. Extended Data Fig

发现三：SSp-bfd局部群体Ca²⁺活动可由全脑BOLD fMRI在多个相关区域解码

概念验证二聚焦全脑尺度。作者反向使用SVM：把每个区域的BOLD体素经过PCA（20成分解释≥75%方差）后解码为SSp-bfd的低通（0.01 Hz）ΔF/F群体活动。每个区域随机抽取显著激活体素（样本量因动物而异，如Subject 3为325 voxels），10× cross-validation取均值，再叠加双向时间延迟的BOLD特征。

在SSp-bfd自身，NRMSE远低于打乱基线（two-sided Wilcoxon signed-rank test, P ≤ 0.005；Fig. 4a,b）。扩展到TH、SSp-bfd、SSp-n、SSp-ul、SSs、MOp、MOs、AUDp、AUDd、STR等多个whisker相关区域时，所有区域都达到显著水平（two-sided Wilcoxon rank-sum test, P ≤ 0.01；Fig. 4c-e, Extended Data Fig. 8）。预测能力可视化在冠状面magnitude map上（Fig. 4d, Extended Data Fig. 9），但与刺激对侧/同侧相关的laterality比较未达显著（two-sided Wilcoxon rank-sum test, P = 0.18）——这是与本文style-ref对象（Lake 2020类工作）的一个数值差异。

Figure 4. Predictive power of BOLD on the SSp-bfd population activity varies across brain regions

Figure 8. Extended Data Fig

Figure 9. Extended Data Fig

这一组结果验证了三件事：第一，局部Ca²⁺活动与脑内多个远端区域的BOLD之间存在可量化的耦合；第二，时序打乱基线排除了简单反射性血管动力学；第三，whisker系统的解剖连接（Aronoff 2010, Yamashita 2018）与功能连接（Esmaeili 2020）预测的区域，都显示出对局部Ca²⁺活动的高信息含量——单细胞层面的局部计算与全脑BOLD所刻画的宏观功能组织在同一套解剖骨架上耦合。

发现四：血管附近神经元的负向Ca²⁺响应与VIP/SOM中间神经元的血管活性环路假说一致

虽然原文未将这一观察单列为独立图，但Fig. 3d–f中的"距离-权重符号"耦合是全文最具机制含义的信号。Fig. 3e,f的trial-averaged ΔF/F曲线显示，最高权重（即负权重）排名靠前的25%神经元在whisker刺激窗口内呈现明显的负向Ca²⁺响应，与同侧SSp-bfd的正向ΔBOLD信号方向相反；这些神经元的欧氏距离到最近血管点显著小于群体中位数。

作者在Discussion中给出的解释是：血管活性interneuron亚型（vasodilatory VIP+，vasoconstricting SOM+）可能在血管附近释放vasoactive神经递质，并通过局部抑制性微环路使相邻excitatory神经元呈现负响应。这一假说与Uhlirova等（2016）、Anenberg等（2015）、Vazquez等（2018）的optogenetic抑制性神经元研究一致，也呼应了Schulz等（2012）发现的星形胶质细胞慢振荡Ca²⁺波对BOLD的贡献。换言之，BOLD信号的神经源不是单一兴奋性群体，而是excitatory、inhibitory与glial三个组分在血管附近的微环路整合结果（Cauli 2004, Lee 2020）。

这套装置的价值由此体现出来：它能以遗传靶向的方式分别记录excitatory与inhibitory亚型在血管附近的活动，并通过多色XCaMP成像（Inoue 2019, Han 2019）或光遗传学操纵（Lee 2010, Bernal-Casas 2017）直接测试VIP/SOM假说，这是fiber photometry与电生理联合fMRI都难以同时实现的组合。

归纳总结和点评

这项工作的主要贡献集中在方法学层面。作者把MRI兼容的物镜、CCD相机、3D打印外壳与acrylic浸没介质组装成一台能在清醒小鼠7T扫描仪内做单细胞Ca²⁺成像的显微镜，且不损失BOLD SNR。装置成本可控、构建简单、保留遗传靶向能力，相对电生理联合fMRI的电磁屏蔽工程更易普及。应用层面，作者用两个概念验证分别回答了"局部BOLD由哪些单细胞预测"和"局部Ca²⁺由哪些远端BOLD区域预测"，并发现神经元到血管的距离决定了Ca²⁺响应相对于BOLD的符号方向。

需要做的克制评价是，最终分析样本量相对较小（N=5只小鼠用于单细胞解码，N=9用于BOLD），且全部基于Rasgrf2×Ai148D一个transgenic line，结论在跨基因型与跨皮层区域的泛化性需要进一步验证。血管距离分析是2D的，受限于单光子显微无法解析深度——神经元的z轴位置同样可能影响与血管的关系。Fig. 3d–f中血管距离与权重符号的耦合虽具机制暗示，但作者自己也指出这是相关性观察，未做optogenetic或药理学因果操纵。SVM解码的因果方向（"哪些神经元预测BOLD" vs "BOLD被什么预测"）也不能等同于神经环路因果。最后，fMRI扫描期间的运动残余与stress残余虽经严格habituation控制（>10天，FWD < 0.05 mm），但仍可能在神经血管耦合数值上留下痕迹（Extended Data Fig. 3, 4）。

更广泛地说，这台显微镜为单细胞分辨率的neurovascular coupling研究提供了一个可扩展的硬件平台。接下来的工作可以沿三条线推进：把细胞类型特异性扩展到inhibitory亚型（VIP、SOM、PV），把成像模式扩展到多色XCaMP以同时记录两种细胞类型，把任务范式从passive whisker刺激扩展到Go/No-Go等行为任务，从而把"局部-全脑"耦合的研究从描述层面推到因果层面。作者也指出，未来可通过几何对齐策略把显微镜FOV与MRI体素精确共配准，并整合two-photon成像以获取layer-specific的神经血管结构信息（Cui 2017）。