脑连接组:结构与功能连接
从 dMRI 与 fMRI 构建结构、功能与有效连接网络,解读模块、枢纽、梯度与布线代价,并检验组间与个体差异。
在网络结构指标与随机图、小世界与无标度网络模型、社团结构与社区发现和网络可控性与稳定性中,度、聚集系数、路径长度、模块度和鲁棒性等工具都有一个前提:网络是“给定的”。在脑科学里,网络不是给定的,而是从影像数据里构建出来的。节点取多少个、边用什么统计量、阈值取多高,每一步都是研究者的选择,每一步都会改变最后的指标。
Bassett & Sporns(2017, Nature Neuroscience)把这一领域概括为“网络神经科学”:用网络描述从分子、细胞到脑区、行为的多个尺度,并把网络当作可以建模、预测和干预的对象,而不只是描述性的图。本专题聚焦宏观尺度的人脑连接组,沿三条线展开:
构建——怎样从弥散磁共振(dMRI)和功能磁共振(fMRI)得到结构网络、功能网络和有效连接模型,以及每一步引入了哪些偏倚:
- 连接组的概念:结构连接、功能连接与有效连接的区别;
- 脑图谱与节点定义;
- 结构连接:弥散 MRI、纤维追踪及其偏倚;
- 功能连接估计:相关与偏相关;
- 预处理对功能连接的影响;
- 有效连接:动态因果模型与格兰杰因果;
- 阈值与加权;
- 加权脑网络指标与归一化。
解读——人脑网络有哪些反复出现的组织特征,它们各意味着什么:
推断——怎样从网络得到关于人群、个体和疾病的结论:
- 组间比较与统计验证;
- 个体指纹与个体差异;
- 病变网络映射;
- 从 fMRI 到结论:分析流程与选择,用一张表汇总每一步的选择及其对结果的影响。
在此基础上加入时间维度、讨论动态功能连接、多层网络与结构–功能耦合的内容见动态功能连接与多层脑网络;不依赖特定模态的网络构建与统计验证原则见网络构建与统计验证。系统性的教科书参照是 Fornito、Zalesky & Bullmore 的《Fundamentals of Brain Network Analysis》(2016),其各章(节点与边、连接矩阵、度与强度、中心性与枢纽、成分/核与俱乐部、路径与导航、基序/小世界与网络经济、模块性、零模型、统计连接组学)与本专题的页面大体一一对应。
Contents知识页面
VIDEO讲解视频
文字稿
大家好,从这一期开始,我们进入应用专题的第一站,大脑。在基础理论部分,网络都是现成给定的,我们只管计算指标。但在脑科学里,网络是研究者从影像数据里一步步构建出来的。这一期我们就沿着构建和解读两条线,把脑网络的来龙去脉讲清楚。
第一个问题,结构连接和功能连接到底有什么不同。第二个问题,脑网络的节点和边,分别是从哪里来的。第三个问题,同一份数据,换一个阈值,结论会不会跟着变。
结构连接是解剖上的道路网,用弥散磁共振和纤维追踪来重建。功能连接是两个脑区信号一起波动的程度,最常用的是皮尔逊相关。有效连接更进一步,问的是谁影响谁,需要借助因果模型来估计。请记住,两个脑区功能连接强,未必有直接的纤维相连。它们可能通过第三个脑区间接相连,也可能只是接收了共同的输入。
第一步,用脑图谱把大脑划分成几十到几百个脑区,作为节点。第二步,估计每对脑区之间的连接强度,作为边的权重。结构网络常用流线数,功能网络常用相关或偏相关。第三步,决定是保留全部权重,还是按阈值变成二值网络。第四步,才轮到我们熟悉的聚集系数、路径长度和模块度。这里每一步都是研究者的选择,而每一个选择都会影响最后的数字。
我们来看一个模拟的脑区相关矩阵,现在逐步提高阈值。阈值越高,保留的边越少,网络密度一路下降。注意平均聚集系数也在跟着变,几乎所有指标都依赖密度。这就带来一个陷阱,如果患者组整体连接偏弱,固定阈值会让网络更稀疏。于是所谓效率降低,可能只是密度降低的副产品。反过来,如果固定密度,连接偏弱的一组又会被迫纳入更多接近噪声的弱边。比较稳妥的做法,是扫描一段密度范围,或者直接分析加权网络。
大脑既要局部分工,又要全局整合,人脑网络正好呈现小世界结构。高聚集意味着局部分工,短路径意味着信息能快速汇总到全脑。全局效率是所有脑区对之间距离倒数的平均,断开的脑区对记为零。识别枢纽时,度并不总是可靠,大模块里的脑区度天然就高。参与系数衡量一个脑区的连接是否均匀分布在多个模块之间。接近1的脑区是跨系统的连接枢纽,接近0的只是模块内部的枢纽。
2011年范登赫费尔和斯波恩斯发现,人脑结构枢纽之间彼此紧密相连。判断富人俱乐部时,一定要与保持度序列的随机网络比较。这些枢纽间的连接数量不多,却承载了大量最短路径,像城市间的高速公路。这些高代价的枢纽,在阿尔茨海默病和精神分裂症研究中常被报告受累。但要注意,组水平的差异显著,并不等于能用来诊断某一个病人。
第一,功能连接是统计相关,不是解剖连接,更不是因果。第二,图谱和阈值都会改变指标,方案要在分析之前确定。第三,判断枢纽要结合多种指标,并检查它在不同阈值下是否稳定。第四,头动等伪迹会制造虚假的连接差异,必须认真控制。建议打开网站上的阈值演示和代码示例,亲手体会阈值的影响。
INTERACTIVE交互可视化
CODE计算实现
10 个完整可运行的 Python 示例,每个都给出代码、输出图和解读;可在浏览器中直接运行,或下载 Jupyter Notebook · Python 脚本。
Papers论文精读
复杂网络中的同步
Synchronization in complex networks
Arenas A, Díaz-Guilera A, Kurths J, Moreno Y, Zhou C
系统综述网络拓扑如何影响同步:Kuramoto 模型的同步起点与路径、主稳定函数与可同步性及应用。
绘制人类大脑皮层的结构核心
Mapping the structural core of human cerebral cortex
Hagmann P, Cammoun L, Gigandet X, Meuli R, Honey CJ, Wedeen VJ, et al.
用弥散谱成像构建全皮层高分辨率结构网络,发现位于后内侧与顶叶皮层的结构核心,并显示结构连接强度能预测静息态功能连接。
复杂脑网络:结构与功能系统的图论分析
Complex brain networks: graph theoretical analysis of structural and functional systems
Bullmore E, Sporns O
系统梳理用图论分析脑结构与功能网络的流程与发现,总结小世界、枢纽、模块化等跨尺度、跨模态反复出现的特征。
脑连接的复杂网络指标:用法与解读
Complex network measures of brain connectivity: uses and interpretations
Rubinov M, Sporns O
按神经科学含义整理脑网络常用图论指标及其解读,并发布 Brain Connectivity Toolbox。
人类连接组的富人俱乐部组织
Rich-club organization of the human connectome
van den Heuvel MP, Sporns O
在 21 人结构网络中发现 12 个双侧枢纽互连远超随机预期,构成承担全局通信的富人俱乐部。
动态功能连接:希望、问题与解释
Dynamic functional connectivity: promise, issues, and interpretations
Hutchison RM, Womelsdorf T, Allen EA, Bandettini PA, Calhoun VD, Corbetta M, et al.
综述动态功能连接的证据与方法,系统指出窗口、噪声、统计检验和神经解释上的问题。
结构脑网络的可控性
Controllability of structural brain networks
Gu S, Pasqualetti F, Cieslak M, Telesford QK, Yu AB, Kahn AE, et al.
把线性控制理论用于人脑结构网络,提出平均、模态、边界三种区域可控性及其系统分布。
网络控制原理预测秀丽隐杆线虫连接组中的神经元功能
Network control principles predict neuron function in the Caenorhabditis elegans connectome
Yan G, Vértes PE, Towlson EK, Chew YL, Walker DS, Schafer WR, et al.
用目标可控性分析线虫连接组,预测运动控制所需的神经元,并以单神经元消融实验验证 PDB 与 DD04/DD05。
共同接收者中的突触组织预测输入神经元之间的连接
Synaptic organization in shared receivers predicts connectivity between input neurons
Zhang X
在两个成年果蝇连接组中,相连输入源的突触在共同接收神经元上靠得更近,这种空间信息可帮助预测未见过的输入神经元之间的连接。
References核心文献
Bullmore & Sporns 2009:Complex brain networks: graph theoretical analysis of structural and functional systems
- 背景:2000 年代中期,图论方法开始大量用于 MRI、EEG/MEG 数据,但流程和指标各不相同。
- 核心发现:这篇综述(Nature Reviews Neuroscience,DOI: 10.1038/nrn2575)系统梳理了从数据到网络的流程,总结了结构与功能脑网络中反复出现的小世界、模块、枢纽和经济布线等特征,并讨论了在疾病中的应用。
- 阅读提示:重点看构建流程示意图和“结构–功能关系”部分;文中关于无标度度分布的讨论要结合后来的研究理解——宏观脑网络的度分布多被描述为指数截断的重尾分布,而非严格幂律。
- 思考题:文中提出大脑在布线成本与拓扑效率之间取舍。如果用流线长度作为成本,你会如何设计一个指标同时衡量这两者?
- 精读:Bullmore 2009
Hagmann 等 2008:Mapping the structural core of human cerebral cortex
- 背景:此前的人脑结构网络研究多用粗分区和弥散张量成像,对交叉纤维和全皮层组织的刻画有限。
- 核心发现:该文(PLoS Biology,DOI: 10.1371/journal.pbio.0060159)用弥散谱成像在 998 个皮层 ROI 上构建个体结构网络,综合多种中心性指标识别出位于后内侧与顶叶皮层的结构核心,并显示结构连接强度与静息态 FC 相关。
- 阅读提示:注意边权对 ROI 面积和纤维长度的校正,以及 k-core/s-core 分解如何补充单一的度指标;样本只有 5 人,结论的普适性需要后续大样本检验。
- 思考题:如果只用流线数而不做长度校正,结构核心的位置可能怎样变化?为什么?
- 精读:Hagmann 2008
Rubinov & Sporns 2010:Complex network measures of brain connectivity: uses and interpretations
- 背景:随着指标数量激增,研究者需要一份“指标字典”说明每个指标的定义、加权和有向版本,以及在神经科学中的合理解释。
- 核心发现:该文(NeuroImage,DOI: 10.1016/j.neuroimage.2009.10.003)把指标分为分离、整合、中心性、韧性等类别,给出统一公式,并配套发布了 Brain Connectivity Toolbox(BCT)。
- 阅读提示:特别注意作者对每个指标“能解释什么、不能解释什么”的告诫,例如功能网络中的路径并不对应真实的信息传递路线。
- 思考题:在相关矩阵构建的功能网络上计算“最短路径”有什么概念上的问题?在什么样的网络上这个指标更有意义?
- 精读:Rubinov 2010
van den Heuvel & Sporns 2011:Rich-club organization of the human connectome
- 背景:已知结构网络中存在枢纽,但枢纽之间是否彼此紧密相连、形成一个核心,尚未被系统检验。
- 核心发现:该文(J Neurosci,DOI: 10.1523/JNEUROSCI.3539-11.2011)在弥散成像结构网络中用归一化富人俱乐部系数证实了人脑富人俱乐部的存在,并展示其在全脑通信中的中心地位。
- 阅读提示:注意归一化步骤——不经随机网络归一化的 随 上升几乎是必然的。
- 思考题:富人俱乐部依赖于度的阈值 。如果你的研究在患者组发现 下降,你会做哪些额外分析来排除“患者组整体连接密度更低”的解释?
- 精读:van den Heuvel 2011
Bassett & Sporns 2017:Network neuroscience
- 背景:到 2010 年代中期,连接组数据已覆盖从分子、细胞到全脑与行为的多个尺度,图论方法也从描述走向建模与控制。
- 核心发现:这篇综述(Nature Neuroscience,DOI: 10.1038/nn.4502)提出“网络神经科学”的整体框架,强调跨尺度、跨模态整合网络数据,并把网络模型用于解释和预测脑功能及其疾病改变。
- 阅读提示:可作为本专题与动态功能连接与多层脑网络、网络可控性之间的“地图”来读。
- 思考题:文中讨论的哪些方向(如多层网络、网络控制、生成模型)最依赖于结构连接的准确性?纤维追踪的假阳性会如何影响这些方向的结论?
补充方法文献(本专题新增页面的主要依据)
| 主题 | 文献 | DOI |
|---|---|---|
| 连接组概念 | Sporns, Tononi & Kötter 2005, PLoS Comput Biol | 10.1371/journal.pcbi.0010042 |
| 结构预测功能 | Honey 等 2009, PNAS | 10.1073/pnas.0811168106 |
| 纤维追踪比赛 | Maier-Hein 等 2017, Nat Commun | 10.1038/s41467-017-01285-x |
| 追踪的解剖准确性 | Thomas 等 2014, PNAS | 10.1073/pnas.1405672111 |
| 假阳性与假阴性 | de Reus & van den Heuvel 2013, NeuroImage | 10.1016/j.neuroimage.2012.12.066 |
| 一致性阈值 | Roberts 等 2017, NeuroImage | 10.1016/j.neuroimage.2016.09.053 |
| SIFT2 | Smith 等 2015, NeuroImage | 10.1016/j.neuroimage.2015.06.092 |
| 节点选择的影响 | Zalesky 等 2010, NeuroImage | 10.1016/j.neuroimage.2009.12.027 |
| 网络估计方法比较 | Smith 等 2011, NeuroImage | 10.1016/j.neuroimage.2010.08.063 |
| 头动伪相关 | Power 等 2012, NeuroImage | 10.1016/j.neuroimage.2011.10.018 |
| 去噪策略比较 | Ciric 等 2017, NeuroImage | 10.1016/j.neuroimage.2017.03.020 |
| 有效连接综述 | Friston 2011, Brain Connectivity | 10.1089/brain.2011.0008 |
| 动态因果模型 | Friston, Harrison & Penny 2003, NeuroImage | 10.1016/S1053-8119(03)00202-7 |
| DCM 使用原则 | Stephan 等 2010, NeuroImage | 10.1016/j.neuroimage.2009.11.015 |
| 格兰杰因果 | Granger 1969, Econometrica | 10.2307/1912791 |
| GC 在神经影像中的应用 | Seth, Barrett & Barnett 2015, J Neurosci | 10.1523/JNEUROSCI.4399-14.2015 |
| fMRI 有向连接的电生理验证 | David 等 2008, PLoS Biol | 10.1371/journal.pbio.0060315 |
| BOLD GC 与血管 | Webb 等 2013, PLoS ONE | 10.1371/journal.pone.0084279 |
| 比例阈值的偏倚 | van den Heuvel 等 2017, NeuroImage | 10.1016/j.neuroimage.2017.02.005 |
| 加权聚集系数 | Onnela 等 2005, Phys Rev E | 10.1103/PhysRevE.71.065103 |
| 带符号模块度 | Rubinov & Sporns 2011, NeuroImage | 10.1016/j.neuroimage.2011.03.069 |
| 功能制图 | Guimerà & Amaral 2005, Nature | 10.1038/nature03288 |
| 7/17 网络 | Yeo 等 2011, J Neurophysiol | 10.1152/jn.00338.2011 |
| 264 功能区 | Power 等 2011, Neuron | 10.1016/j.neuron.2011.09.006 |
| 模块化脑网络 | Sporns & Betzel 2016, Annu Rev Psychol | 10.1146/annurev-psych-122414-033634 |
| 连接枢纽损伤 | Warren 等 2014, PNAS | 10.1073/pnas.1322173111 |
| 网络经济性 | Bullmore & Sporns 2012, Nat Rev Neurosci | 10.1038/nrn3214 |
| 布线优化(线虫) | Chen, Hall & Chklovskii 2006, PNAS | 10.1073/pnas.0506806103 |
| 非最优布局与长程投射 | Kaiser & Hilgetag 2006, PLoS Comput Biol | 10.1371/journal.pcbi.0020095 |
| 指数距离规则 | Ercsey-Ravasz 等 2013, Neuron | 10.1016/j.neuron.2013.07.036 |
| 几何对连接组的贡献 | Roberts 等 2016, NeuroImage | 10.1016/j.neuroimage.2015.09.009 |
| 生成模型(功能) | Vértes 等 2012, PNAS | 10.1073/pnas.1111738109 |
| 生成模型(结构) | Betzel 等 2016, NeuroImage | 10.1016/j.neuroimage.2015.09.041 |
| 空间置换检验 | Alexander-Bloch 等 2018, NeuroImage | 10.1016/j.neuroimage.2018.05.070 |
| 主梯度 | Margulies 等 2016, PNAS | 10.1073/pnas.1608282113 |
| 扩散映射 | Coifman 等 2005, PNAS | 10.1073/pnas.0500334102 |
| 梯度工具箱 | Vos de Wael 等 2020, Commun Biol | 10.1038/s42003-020-0794-7 |
| 皮层下梯度与图谱 | Tian 等 2020, Nat Neurosci | 10.1038/s41593-020-00711-6 |
| 连接组指纹 | Finn 等 2015, Nat Neurosci | 10.1038/nn.4135 |
| 差异可识别性 | Amico & Goñi 2018, Sci Rep | 10.1038/s41598-018-25089-1 |
| 稳定的个体因素 | Gratton 等 2018, Neuron | 10.1016/j.neuron.2018.03.035 |
| 基于连接组的预测建模 | Shen 等 2017, Nat Protoc | 10.1038/nprot.2016.178 |
| FC 重测信度 | Noble, Scheinost & Constable 2019, NeuroImage | 10.1016/j.neuroimage.2019.116157 |
| 病变网络映射 | Boes 等 2015, Brain | 10.1093/brain/awv228 |
| 连接组与症状定位综述 | Fox 2018, NEJM | 10.1056/NEJMra1706158 |
| LNM 的方法学质疑 | van den Heuvel 等 2026, Nat Neurosci | 10.1038/s41593-025-02196-7 |
| 对质疑的回应 | Pini, Salvalaggio & Corbetta 2026, Nat Neurosci | 10.1038/s41593-026-02319-8 |
| LNM 的样本量与阈值 | Sperber & Dadashi 2020, Brain | 10.1093/brain/awaa094 |
| 卒中失连接预测 | Salvalaggio 等 2020, Brain | 10.1093/brain/awaa156 |
| 分析灵活性(NARPS) | Botvinik-Nezer 等 2020, Nature | 10.1038/s41586-020-2314-9 |
| HCP 概述 | Van Essen 等 2013, NeuroImage | 10.1016/j.neuroimage.2013.05.041 |
| HCP 静息态 | Smith 等 2013, NeuroImage | 10.1016/j.neuroimage.2013.05.039 |
| HCP 最小预处理 | Glasser 等 2013, NeuroImage | 10.1016/j.neuroimage.2013.04.127 |
| 疾病连接组学综述 | Fornito, Zalesky & Breakspear 2015, Nat Rev Neurosci | 10.1038/nrn3901 |
Data & tools相关数据与软件
Further延伸资源
| 类型 | 资源 | 说明 |
|---|---|---|
| 视频 | 功能连接是什么?(B 站) | 中文,7 分钟,零基础入门 |
| 视频 | Introduction to Brain Network Analysis – Part 1/2:B 站转载;YouTube 原始链接(需境外网络) | 英文,约 31 分钟,从相关矩阵、阈值到图论指标;B 站版为转载,授权情况不明 |
| 教材 | Fornito, Zalesky & Bullmore, Fundamentals of Brain Network Analysis, Academic Press, 2016,doi.org/10.1016/C2012-0-06036-X | 系统教材:第 2 章节点与边、第 3 章连接矩阵、第 4 章度与强度、第 5 章中心性与枢纽、第 6 章核与俱乐部、第 7 章路径与导航、第 8 章基序/小世界与网络经济、第 9 章模块性、第 10 章零模型、第 11 章统计连接组学 |
| 教材 | Sporns, Networks of the Brain, MIT Press, 2011 | 从复杂网络角度讲解连接组、结构–功能关系与脑网络演化的专著 |
| 综述 | Bullmore & Sporns 2009,doi.org/10.1038/nrn2575 | 脑网络图论的奠基综述,建议首先阅读 |
| 综述 | Bassett & Sporns 2017,doi.org/10.1038/nn.4502 | 网络神经科学的整体框架 |
| 方法 | Rubinov & Sporns 2010,doi.org/10.1016/j.neuroimage.2009.10.003 | 指标字典,撰写方法部分时查阅 |
| 综述 | Fornito, Zalesky & Breakspear 2015,doi.org/10.1038/nrn3901 | 图论指标在脑疾病中的发现与解读陷阱,配合“生物医学应用”一节阅读 |
| 工具 | Brain Connectivity Toolbox(MATLAB;Python 移植版 bctpy) | 加权/有向/带符号版本的聚集、效率、模块度、参与系数、富人俱乐部与零模型 |
| 教程 | Nilearn:从图谱提取时间序列并构建功能连接组 | Python 标准流程,ConnectivityMeasure
支持相关、偏相关、切空间等;可在 Colab 中运行 |
| 工具 | BrainSpace | 连接组梯度的计算、对齐与可视化(Python/MATLAB) |
| 工具 | GRETNA / BrainNet Viewer / DPABI | 国内团队开发的 MATLAB 图形界面工具,不熟悉编程也能完成建网、图论分析和出图 |
| 数据 | Human Connectome Project | 多模态 MRI 与行为数据;需注册并同意数据使用条款 |
| 数据 | ENIGMA Toolbox 内置 HCP 组平均连接组 | 一行代码读取组平均结构/功能连接矩阵(多种图谱),适合作为分析起点 |
| 交互网站 | FlyScope · 果蝇连接组浏览器 | 基于 MaleCNS 与 FlyWire 两个成年果蝇全脑连接组,可在单个神经元的分支上查看输入落在哪里;配合论文精读理解“节点内部结构”也携带网络信息。中英文界面 |
Open questions开放问题
- 假设你要比较抑郁症患者与健康对照的静息态功能网络全局效率。请写出一份分析前就确定下来的方案:节点图谱、是否做全局信号回归、阈值策略、头动控制、多重比较校正。哪一步你最没有把握?为什么?
- 结构连接相对稳定,功能连接却随状态变化。结合你自己的研究方向(例如影像、电生理、康复工程),讨论在什么问题上应优先使用结构网络、什么问题上应优先使用功能网络,以及什么情况下二者都需要。
- 一项研究用 fMRI 格兰杰因果报告“患者组前额叶对杏仁核的自上而下控制减弱”。请列出至少三种不涉及神经因果方向改变、却能产生这一结果的解释,并设计一个能区分它们的补充分析。
- 在结构网络中,长程连接既是追踪最不可靠的连接,又是承载富人俱乐部与网络效率的关键连接。这对“疾病选择性损害长程/枢纽连接”的结论意味着什么?你会用什么零模型或对照分析来检验?
- 连接组梯度与模块划分对同一份 FC 数据给出了“连续”与“离散”两种描述。设想一个研究问题,说明哪一种描述更合适,以及为什么。









