Nature Methods | 單細胞鈣成像與全腦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能無創記錄全腦血流動力學響應，但只能提供局部神經活動的慢速、間接均值估計；單光子或雙光子成像可在單細胞層面記錄Ca2?或電壓信號，卻需要侵入式手術，且只能覆蓋少量皮層區域（Fig. 1a–j）。兩種模態各有取舍，都只能回答腦功能的一個尺度。

過去十余年里，研究者嘗試用fiber photometry把Ca2?成像與BOLD fMRI縫合在同一只麻醉或清醒動物上。Schulz等（2012）在麻醉大鼠上做了開創性聯合記錄，指出BOLD信號不僅反映神經元源，還與星形膠質細胞的慢振蕩Ca2?波相關；Lake等（2020）用多束光纖把全皮層Ca2?信號與BOLD并行采集，證明BOLD可在高空間精度下由局部神經活動部分重建。然而，fiber-based方法無法分辨單細胞，把神經元、神經膠質與neuropil信號混在一起，丟失了細胞類型特異性信息（Fig. 1a–j）。電生理聯合fMRI雖然能拿到單細胞分辨率，但電磁屏蔽工程復雜且代價高昂。

本文的目標正是填補這一空白：組裝一臺MRI兼容的單光子顯微鏡，能在清醒行為小鼠的7T掃描儀內同步采集單細胞Ca2?活動與全腦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 mm2，像素分辨率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幀率采集單細胞Ca2?活動。

數據采集與分析：MRI使用7T BioSpec 70/20，GE-EPI序列（TE/TR=15/1500 ms，7片冠狀面，體素0.2×0.2×0.8 mm）。Ca2?成像經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 mm2、像素1.1 μm、光學分辨率2.5 ± 0.24 μm FWHM（Fig. 1h–j），既不損失BOLD SNR，也不損失細胞分辨率。

這一步的工程意義在于，過去fMRI-光學聯合記錄只能在fiber photometry下做，把神經-血管問題壓縮為群體均值；本文的裝置首次把研究單元下沉到單細胞水平，并且保留遺傳靶向能力（Ca2?指示劑由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信號可由單細胞群體Ca2?活動解碼，且神經元-血管距離決定解碼權重的符號

概念驗證一聚焦局部尺度。清醒小鼠在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）；同步Ca2?成像上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 Ca2? 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 Ca2?響應（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局部群體Ca2?活動可由全腦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

這一組結果驗證了三件事：第一，局部Ca2?活動與腦內多個遠端區域的BOLD之間存在可量化的耦合；第二，時序打亂基線排除了簡單反射性血管動力學；第三，whisker系統的解剖連接（Aronoff 2010, Yamashita 2018）與功能連接（Esmaeili 2020）預測的區域，都顯示出對局部Ca2?活動的高信息含量——單細胞層面的局部計算與全腦BOLD所刻畫的宏觀功能組織在同一套解剖骨架上耦合。

發現四：血管附近神經元的負向Ca2?響應與VIP/SOM中間神經元的血管活性環路假說一致

雖然原文未將這一觀察單列為獨立圖，但Fig. 3d–f中的"距離-權重符號"耦合是全文最具機制含義的信號。Fig. 3e,f的trial-averaged ΔF/F曲線顯示，最高權重（即負權重）排名靠前的25%神經元在whisker刺激窗口內呈現明顯的負向Ca2?響應，與同側SSp-bfd的正向ΔBOLD信號方向相反；這些神經元的歐氏距離到最近血管點顯著小于群體中位數。

作者在Discussion中給出的解釋是：血管活性interneuron亞型（vasodilatory VIP+，vasoconstricting SOM+）可能在血管附近釋放vasoactive神經遞質，并通過局部抑制性微環路使相鄰excitatory神經元呈現負響應。這一假說與Uhlirova等（2016）、Anenberg等（2015）、Vazquez等（2018）的optogenetic抑制性神經元研究一致，也呼應了Schulz等（2012）發現的星形膠質細胞慢振蕩Ca2?波對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掃描儀內做單細胞Ca2?成像的顯微鏡，且不損失BOLD SNR。裝置成本可控、構建簡單、保留遺傳靶向能力，相對電生理聯合fMRI的電磁屏蔽工程更易普及。應用層面，作者用兩個概念驗證分別回答了"局部BOLD由哪些單細胞預測"和"局部Ca2?由哪些遠端BOLD區域預測"，并發現神經元到血管的距離決定了Ca2?響應相對于BOLD的符號方向。

需要做的克制評價是，最終分析樣本量相對較?。∟=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）。

請打分

這篇剛剛登上 Nature Methods 的研究，是否實至名歸？我們邀請您作為“云審稿人”，一同品鑒。精讀全文后，歡迎在匿名投票中打分，并在評論區分享您的深度見解。

分享人：BQ

審核：PsyBrain 腦心前沿編輯部

你好，這里是「PsyBrain 腦心前沿」

專注追蹤全球認知神經科學的最尖端突破

視野直擊 Nature, Science, Cell 正刊 及核心子刊與頂級大刊

每日速遞「深度解讀」與「前沿快訊」

科研是一場探索未知的長跑，但你無需獨行。歡迎加入PsyBrain 學術社群，和一群懂你的同行，共同丈量腦與心智的無垠前沿。

點擊卡片進群，歡迎你的到來

一鍵關注，點亮星標 ? 前沿不走丟！

一鍵分享，讓更多人了解前沿