生物行移动版

大脑的功能是由从分子到宏观网络的多个尺度共同决定的。然而，如何将特定神经回路的功能动态（通过全脑功能磁共振成像测量）与其潜在的分子和细胞特征（通过空间转录组学测量）直接联系起来，一直是一个巨大的挑战。2026年4月18日发表于 Nature Communications 的一项研究中，中国科学院生物物理研究所的王晓群团队（Zhang C., Zou Y., Zhuo L.等人共同第一作者）开发了一种创新框架：在同一只清醒小鼠中，结合光遗传fMRI（opto-fMRI） 和空间转录组学（opto-ST），对内侧前额叶皮层（mPFC） 或海马下托（SUB） 进行刺激，并分析它们诱发的不同丘脑功能磁共振成像（BOLD）反应模式。通过多变量空间建模和主成分分析，他们发现能量代谢和核心生化过程是这些基于一般线性模型（GLM）激活模式的分子基础。更重要的是，这些结果提示常规线性模型可能主要捕获线性的基因-过程关系，因此需要结合非线性计算方法来进行功能磁共振成像分析。该研究提供了一个将局部转录组学与大规模功能成像特征相连接的通用框架。

研究亮点速览

1. 方法论创新：首次在同一批小鼠中完成清醒opto-fMRI（全脑BOLD响应）和空间转录组学（10x Visium，从刺激靶点及其下游丘脑核团取材），实现了真正的“跨尺度同体整合”。

2. 双回路对比：选择两个已知的皮质-丘脑回路：mPFC → 丘脑（参与认知控制，与抑郁症/焦虑相关）和SUB（海马下托）→ 丘脑（参与记忆和空间处理）。光遗传刺激这两个脑区，在全脑范围内（尤其是丘脑核团）产生了不同的BOLD响应模式。

3. 多变量建模：使用偏最小二乘（PLS）回归，将丘脑中每个体素的BOLD激活强度与空间转录组数据（基因表达矩阵）相关联，识别出与mPFC-和SUB-诱发的BOLD信号显著相关的基因集合。

4. 富集的生物过程：与GLM衍生的BOLD激活相关的基因，显著富集于能量代谢通路（如氧化磷酸化、三羧酸循环、ATP合成）和核心细胞过程（如RNA剪接、蛋白质翻译）。这些结果在数据驱动的主成分分析（PCA）中得到交叉验证。

5. 对“神经血管耦合”的启示：BOLD信号依赖于神经血管耦合。本研究的分子富集表明，代谢途径的基因表达水平是决定特定脑区BOLD响应幅度（和效率）的基础因素之一。

6. 方法学警示：GLM假设线性，而神经血管耦合可能是非线性的。作者主张将数据驱动的非线性模型（如PLS）与GLM结合使用，以获得更完整的生物学图景。

背景：脑连接组与转录组的“脱节”

已知：全脑功能连接（例如来自静息态或任务态fMRI）与基因表达模式存在相关性（Allen人脑图谱）。某些精神病/神经疾病的易感基因在功能连接异常的区域“富集”。

局限：通常fMRI和转录组数据来自不同个体，甚至不同物种（如小鼠连接数据 vs. 人尸检转录组），混杂了群体差异。无法建立因果性的回路-基因关系（相关性不代表回路激活直接驱动的转录特征）。宏观BOLD信号如何与局部的转录谱对接，机制上不透明。

本研究的解决方案：光遗传：因果激活特定回路（mPFC-丘脑或SUB-丘脑）。清醒动物功能磁共振成像：在不使用麻醉剂的情况下测量全脑BOLD效应（避免麻醉干扰神经血管耦合）。空间转录组：从相同动物的相同脑区（丘脑）提取组织进行测序，保持空间完整性。

方法学核心：从回路激活到分子注释

1. 动物准备：ChR2病毒注入mPFC或SUB，光纤植入刺激位点。

2. 清醒功能磁共振成像：在9.4T MRI中，对清醒小鼠进行光遗传刺激（蓝光，10秒开/关），全脑BOLD采样。

3. 成像后处理：计算mPFC或SUB刺激的GLM激活图（BOLD信号 % 变化，以体素为单位），重点关注丘脑核团（前核、腹核、髓板内核等）。

4. 空间转录组：对同一批功能磁共振成像后的新鲜脑组织进行速冻、切片，进行10x Visum空间转录组（探针覆盖全转录组）。从与功能磁共振成像相同的丘脑亚区捕获基因表达数据。

5. 多变量关联（PLS）：输入矩阵X：体素级基因表达矩阵（每体素数千个基因）。输入矩阵Y：体素级GLM BOLD激活值（来自mPFC或SUB刺激）。PLS提取潜在变量，最大化X与Y之间的协方差。通过特征选择，识别出与BOLD激活最相关的基因模块。

6. 功能富集：Metascape和ClueGO进行基因本体（GO）和通路富集。

关键结果

1. mPFC和SUB刺激产生不同的丘脑BOLD图谱

• mPFC刺激：在丘脑束旁核（PF）、中央内侧核（CM） 等中线/髓板内核团激活强烈，这与mPFC调控觉醒和认知的已知功能一致。

• SUB刺激：强激活丘脑前核（AD、AV） 、外侧背核（LD） 等边缘系统相关核团，与海马-丘脑记忆回路一致。

• 共激活区域：某些丘脑核团（如腹侧基底复合体）对两种刺激均有响应，但模式不同。

2. 与GLM激活模式相关的基因富集于能量代谢

工作流程：PLS识别出与BOLD激活体素模式显著相关的“特征基因”列表。

富集结果（mPFC和SUB回路均一致）：氧化磷酸化 (OXPHOS)：复合物I-V亚基基因（如Ndufv1、Uqcrb、Cox6a1）。三羧酸循环（如Idh3a、Sdha）。ATP合成（如Atp5a1）。RNA剪接和蛋白质翻译相关基因。

3. 数据驱动的主成分分析（PCA）验证

对丘脑空间转录组数据进行PCA，发现第一主成分（PC1） 的基因载荷与PLS识别的特征基因高度重叠，且PC1得分与BOLD激活强度显著相关。这表明能量代谢基因的表达梯度是丘脑BOLD响应异质性的主要内在决定因素。

4. 回路特异性分子特征

尽管能量代谢是共同基础，但PLS也识别出回路特异性基因：mPFC-丘脑回路：富集于突触传递和神经递质受体基因（如Gria1、Grin1）。SUB-丘脑回路：富集于细胞骨架和轴突导向基因（如Map2、Dcx）。这提示不同回路的功能特异性也体现在分子水平。

讨论与展望

本研究通过同体整合光遗传fMRI与空间转录组学，首次在因果水平上揭示了特定神经回路激活与下游脑区分子特征的直接关联。主要发现：能量代谢基因表达是BOLD信号强度的主要分子基础，而回路特异性基因则赋予功能特异性。这一框架为理解脑疾病（如抑郁症、阿尔茨海默病）中回路功能障碍与分子病理的联系提供了新工具。未来可扩展至疾病模型，探索病理状态下回路-基因关系的改变。