生物医学网络科学Biomedical Network Science
经典论文难度:进阶

人类连接组的富人俱乐部组织

Rich-club organization of the human connectome

van den Heuvel MP, Sporns O · Journal of Neuroscience 2011 · DOI: 10.1523/JNEUROSCI.3539-11.2011

富人俱乐部枢纽连接组弥散张量成像零模型网络韧性

一句话

在 21 人结构网络中发现 12 个双侧枢纽互连远超随机预期,构成承担全局通信的富人俱乐部。

相关专题:脑连接组:结构与功能连接、网络可控性与稳定性

一句话看懂

大脑里的枢纽区域不仅各自连接多,它们彼此之间也连得特别紧密,紧密程度超出“仅凭度就能预期”的水平。这群互连的枢纽,包括楔前叶、额上回、顶上皮层以及海马、壳核、丘脑,构成人脑的“富人俱乐部”。大多数远距离区域之间的最短路径都要经过它们。

研究背景与问题

此前的研究(如 Hagmann 等, 2008)已经找到若干皮层枢纽,并发现它们在多种脑疾病中容易受累。但“有枢纽”只是节点层面的描述,还没有回答枢纽之间如何组织:它们是分散在各模块里各管一片,还是彼此相连、组成一个更高层次的骨干?

网络科学中的富人俱乐部系数(Zhou & Mondragón 2004;Colizza 等 2006 提出归一化版本)正好刻画这种更高阶的拓扑。本文把它系统地用于人脑结构网络,并同时纳入皮层下结构。

核心思想与方法

  1. 数据:21 名健康被试的弥散张量成像(DTI),确定性纤维追踪。
  2. 节点:FreeSurfer 分区得到 82 个区域(68 个皮层 + 14 个皮层下);另用 1170 个皮层小区域的高分辨率版本做验证。
  3. 边:二值网络,以及按纤维数、各向异性分数(FA)等加权的网络。
  4. 富人俱乐部检验:对每个度阈值 kk 计算富人俱乐部系数,与保持度序列的随机网络比较。
  5. 功能意义分析:统计最短路径经过富人俱乐部的比例;模拟定向破坏富人俱乐部与随机破坏对全局效率的影响。

富人俱乐部系数:设度大于 kk 的节点有 N>kN_{>k} 个,它们之间共有 E>kE_{>k} 条边,则 ϕ(k)=2E>kN>k(N>k−1)\phi(k)=\frac{2E_{>k}}{N_{>k}(N_{>k}-1)} 即“富人”子网络的密度。由于度高的节点本来就更可能相连,必须做归一化:ϕnorm(k)=ϕ(k)/ϕrand(k)\phi_{\text{norm}}(k)=\phi(k)/\phi_{\text{rand}}(k),ϕrand\phi_{\text{rand}} 取自保持度序列的随机网络。ϕnorm(k)>1\phi_{\text{norm}}(k)>1(且统计显著)才说明存在富人俱乐部。

主要发现

  1. 人脑结构网络存在富人俱乐部:在二值网络中,kk 约为 11 到 17 的范围内 ϕnorm(k)\phi_{\text{norm}}(k) 明显大于 1;以纤维数或 FA 加权、以及 1170 区高分辨率网络中也观察到显著的富人俱乐部组织。
  2. 成员是 12 个双侧区域:楔前叶、额上回、顶上皮层,以及皮层下的海马、壳核和丘脑(左右各一)。
  3. 富人俱乐部是全局通信的中转站:非富人俱乐部区域之间的最短路径中,约 89% 至少经过一个富人俱乐部节点;随机选取同样数量的节点时,这一比例约为 32%。
  4. 定向破坏代价更大:全部切断富人俱乐部内部连接时,全局效率最多下降约 18%;随机破坏同等数量连接时最多下降约 5%,前者约为后者的三倍。
  5. 皮层的富人俱乐部节点大多是“连接枢纽”,即其连接跨越多个模块,而不是局限在单个模块内部。

怎样读这篇论文

先读引言最后一段和方法中富人俱乐部系数的定义,确认你理解为什么要除以随机网络的系数。然后看展示 ϕ(k)\phi(k)、ϕrand(k)\phi_{\text{rand}}(k) 和 ϕnorm(k)\phi_{\text{norm}}(k) 随 kk 变化的曲线图,这是全文核心证据;注意曲线在哪个 kk 区间显著。接着看富人俱乐部成员的脑图和最短路径统计。加权版本与高分辨率版本是稳健性检验,读结论即可。讨论部分关于“高成本、高容量骨干”的推测可以和后续论文对照阅读。

局限与后续工作

  • 零模型与阈值选择:ϕnorm\phi_{\text{norm}} 依赖随机化方式,结论只在一段 kk 区间内成立,“俱乐部成员”的确定也带有选择性。不同分区、不同密度会改变成员名单。
  • DTI 与确定性追踪:交叉纤维和长程连接的追踪误差会直接影响高度节点的识别;样本只有 21 人。
  • 攻击模拟是理论推演:全局效率下降说明的是拓扑上的脆弱性,不等于真实脑损伤后的功能后果。
  • 后续发展:van den Heuvel 等(2012, PNAS)进一步把富人俱乐部描述为“高成本、高容量”的通信骨干;van den Heuvel 等(2013, JAMA Psychiatry)报告精神分裂症患者富人俱乐部节点之间的连接密度降低;Crossley 等(2014, Brain)发现多种脑疾病的病灶倾向集中在枢纽区域。这些临床研究多为横断面相关证据:富人俱乐部改变可能是病因、结果或伴随现象,不能直接推出“刺激或保护枢纽就能治疗”。

关联知识

延伸思考

  1. 为什么未经归一化的 ϕ(k)\phi(k) 随 kk 增大几乎总会上升?请用一个度异质的随机网络解释。
  2. 如果患者组的富人俱乐部连接减少,你能想到哪些非神经因素(如头动、纤维追踪参数、脑萎缩)可能制造这种差异?
  3. 定向破坏富人俱乐部使全局效率下降更多。这能否说明富人俱乐部区域是最好的神经调控靶点?还需要哪些证据?