1. 为什么我最终把手头的fMRI和PET数据交给JuSpace去做跨模态关联先说一个我实际工作中反复遇到的窘境。手里有一张统计 parametric map——比如精神分裂症患者相对健康对照的灰质体积差异图或者某个认知任务下全脑激活图。传统上我们的报告方式是激活簇位于前扣带、脑岛……然后讨论部分引几篇文献说前人报道前扣带与多巴胺系统高度相关。这种讨论其实相当模糊因为它本质上是在做视觉层面的主观对齐根本经不起定量推敲。另一个方向是拿AAL模板或Hammers模板把大脑切成若干ROI在每个ROI里同时提取fMRI效应值和某些PET示踪剂的结合电位然后做几十个点的相关分析。这种方法有一个非常现实的问题ROI的数量直接决定了相关分析的自由度而且ROI平均本身会抹掉脑回水平上的精细空间信息。JuSpace全称是Spatial Correlation between Neurotransmitter Maps and Neuroimaging Data解决的就是这个问题。它在这个领域里做了一件称得上范式转换的事把神经递质系统富集区域明确定义为一张张存储在MNI空间里的三维分布图谱然后将用户提供的目标脑图激活图、差异图、连接图、形态学图都可以与这些递质图谱在体素水平做空间相关再用置换检验估算显著性。这样一来这个激活模式与多巴胺系统是否在空间上共定位就从一句印象流表述变成了一个带p值的定量结论。需要说明的是JuSpace不是用来做PET定量分析的它本身并不重构PET数据而是作为一个脑图与已知递质分布先验图谱之间建立统计关联的桥梁。这个定位使得它在没有任何新PET采集的情况下就能让常规的fMRI或sMRI研究获得分子水平上的解释维度。对于没有核医学设施、拿不到PET数据的课题组这套方案几乎是当前唯一低门槛的路径。我觉得每一个准备挑战脑影像分子解释这个问题的人都应该先把JuSpace的底层设计逻辑和统计框架看清楚。它不是一个把数据丢进去就能自动出结论的黑盒子恰恰相反它对输入图的空间基准、平滑程度、掩模范围都有要求。如果这些细节处理不好算出来的相关性很可能只是平滑伪影造成的假象。2. JuSpace背后那套递质分布先验图谱到底靠不靠谱以及它包含哪些系统2.1 图谱的出处与预处理过程JuSpace目前封装的神经递质图谱主要来源于一些公开的PET示踪剂成像研究其中相当一部分数据是通过JuSpace作者所在团队与多个欧洲核医学中心合作获取的。博客里常有人问这些图谱具体是哪来的其实在JuSpace的README和原始论文里写得很清楚所有PET图谱都做了MNI空间标准化部分图谱经过部分容积效应校正partial volume correction并且以统一的体素分辨率存储。这一点非常关键。PET原始图像的空间分辨率通常只有3-6mm而且受头动、示踪剂非特异性结合、部分容积效应影响很大。JuSpace作者团队把不同来源、不同示踪剂、不同扫描仪的数据统一重新处理到MNI空间用一套相对固定的流程空间归一化、平滑、mask约束来降低数据之间的异质性。据我所知不同递质图谱之间的平滑核并不完全一致有些图谱源自原始分辨率较高的研究有些则经过了额外平滑以匹配群体模板。这种不一致在实际使用中直接表现为个别递质图谱参与相关分析时更容易得到显著结果这不是因为生物学效应更强而是因为空间平滑度不同导致统计功效不同。后面我在踩坑章节会专门展开这个问题。2.2 目前内置的递质系统一览JuSpace v1.x版本内置了多个神经递质系统图谱覆盖了精神科和神经科研究里最常被追问的那些系统。我建议初次使用者先在自己的电脑上打开JuSpace_Distribution_Map.m或者在MATLAB路径里翻一下/Maps目录把图谱文件名完整看一遍因为不同版本里图谱数量略有增减。截至我目前常用的版本主要包含血清素系统5-HT1A受体、5-HT1B受体、5-HT2A受体、5-HT4受体、血清素转运体SERT多巴胺系统D1受体、D2受体、多巴胺转运体DAT、多巴胺合成相关图谱如F-DOPA乙酰胆碱系统毒蕈碱型受体M1、烟碱型受体α4β2*通常用F-A-85380或2-FA示踪、乙酰胆碱转运体VAChT去甲肾上腺素系统去甲肾上腺素转运体NET氨基酸递质系统GABA-A受体/苯二氮卓位点如在癫痫研究中常用的FMZ代谢型谷氨酸受体mGluR5其他大麻素CB1受体等有趣的是你仔细观察这套清单就会发现它恰恰对应了精神药理学中最主要的药物靶点家族。精神分裂症的D2/5-HT2A假说、抑郁症的5-HT转运体/SERT假说、焦虑障碍的GABA-A/5-HT1A假说、阿尔茨海默病的胆碱能假说、成瘾的多巴胺/大麻素系统假说——JuSpace可以说是一站式把这几个最受瞩目的靶点全部摆在用户面前。表格大概长这样系统代表性图谱常见示踪剂相关疾病/药物领域多巴胺D1、D2、DAT、F-DOPA[11C]SCH23390、[11C]raclopride、[11C]PE2I精神分裂症、帕金森、成瘾血清素5-HT1A、5-HT1B、5-HT2A、5-HT4、SERT[11C]WAY100635、[18F]altanserin、[11C]DASB抑郁症、焦虑障碍、自杀行为乙酰胆碱VAChT、M1、α4β2*[18F]FEOBV、[11C]P943阿尔茨海默病、路易体痴呆去甲肾上腺素NET[11C]MRB抑郁症、ADHDGABA/谷氨酸GABA-A/BZFMZ、mGluR5[11C]flumazenil、[11C]ABP688癫痫、焦虑、双相、精神分裂症大麻素CB1[11C]OMAR、[18F]FMPEP-d2奖赏、食欲、应激相关疾病2.3 分辨率和坐标空间所有图谱都是3D NIfTI文件一般位于3mm或2mm的MNI152空间。在做相关分析前JuSpace内部会把你传入的目标图也重采样到与图谱一致的空间网格上。这个机制看起来是自动处理但我的经验是——最好不要完全依赖它的自动重采样。如果输入图的分辨率是1mm的T1像衍生图重采样到3mm网格后精细边缘信息会被抹掉反过来如果输入图是8mm平滑后的SPM激活图问题不大因为平滑后低频空间信息占主导。比较稳妥的做法是在SPM12的Coregister里面先把目标图与某一标准模板重采样到2mm或3mm然后检查一下重采样结果与图谱的边缘对齐情况再交给JuSpace做正式分析。这个细节在英文论坛里经常被提及国内分享相对少。很多初学用户报错image dimension mismatch之后简单地把所有图像reslice成同样分辨率就重新跑其实并没有意识到reslice过程中插值方法的选择三线性还是最近邻以及是否先做mask约束会直接影响最终空间相关值。对脑内比较光滑的图三线性插值问题不大对非常稀疏的图最近邻插值能避免边缘扩散带来的假相关。3. 从NIfTI文件到第一张相关表JuSpace的运行流程与参数取舍3.1 输入File Folder的组织方式JuSpace在Github上以MATLAB工具包的形式发布目前没有打包成独立的APP也没有Python版本社区有一些非官方的封装但我不建议初学者直接使用验证程度参差不齐。把工具包解压到磁盘上然后在MATLAB里addpath(genpath(/yourpath/JuSpace/))就可以开始用了。官方Demo运行的是一个交互式GUI在MATLAB命令行输入JuSpace会弹出主界面左侧选择图谱类型Neuromaps / JuSpace Legacy、目标统计图、协变量图covariates可以是年龄、性别编码图等右侧勾选需要的递质图谱、设置相关方法和置换次数。如果你要批量处理大量受试者或多种疾病组对比图界面操作会非常低效更推荐直接用命令行函数调用% 其中spm_pet_map是要做关联的统计图NIfTI路径 % condition是可选协变量图没有就填空数组 % outputdir是结果输出目录 [results] JuSpace(map, spm_pet_map.nii, ... maps, {5HT1A, D2, VAChT}, ... covariates, [], ... outputdir, /results/, ... n_perm, 1000);这段代码是我自己根据常用数据结构记忆写的不同版本参数名可能微调但基本思路一致。强烈建议先跑一遍官方自带的sample数据确认环境、路径、输出目录权限都没问题再开始处理自己的真实数据。3.2 相关方法选Pearson还是SpearmanJuSpace提供了Pearson相关和Spearman相关两种空间相关算法。它们之间的差异其实不算大但有一点必须提前想清楚如果你预期自己的效应图与递质图谱之间存在单调但非线性的关系Spearman更稳健如果效应图数值分布本身就接近正态且光滑Pearson统计功效更高。我的个人习惯是只要不是极端离群值很明显主分析用Pearson敏感性分析里把Spearman也跑一遍。如果两种方法结论一致报告里写起来就很硬气如果不一致那通常说明相关关系主要由少数体素驱动需要具体看散点图。3.3 吻侧-尾侧分割选项的意义JuSpace界面上有一个容易被忽略但很重要的设置是否将纹状体分割为吻侧rostral/尾侧caudal亚区。为什么要单独提这个因为在多巴胺系统研究里纹状体并不是一个功能均质的整体。腹侧纹状体包括伏隔核与奖赏、动机关系更密切尾状核参与认知控制壳核更多涉及运动功能。尾侧纹状体的D2受体密度和背侧感觉运动环路的联系紧密吻侧则更多与边缘-前额叶环路交互。如果不做分割直接将全纹状体的D2图谱与某个认知任务激活图做相关很可能会因为内部亚区方向的信号互相抵消而得到不显著结果。这个选项本质上是在区域特异性和统计功效之间做一个折中。分割后每个亚区的体素数量下降相关分析的稳定性会略降但如果生物学信号确实集中在某个亚区收益远大于代价。我见过一篇关于强迫症的研究把全纹状体D2相关做成阴性结果后来用吻侧/尾侧分割重新分析发现是吻侧亚区显著这算是个值得注意的教训。3.4 掩模和协变量处理再强调一次掩模问题。许多用户从SPM输出中直接拿到的统计图往往带有全脑灰质mask或白质限制。理论上传给JuSpace的图应该是全脑范围内每个体素都有明确数值的图而不是大部分体素被NaN或0占满的稀疏图。原因在于空间相关计算把所有体素都当作观测单位如果目标图里非脑区域的数值是0而递质图谱里对应区域是非0少数图谱的边缘mask不够精细就会人为制造出一批包含大量零点的点簇把相关值推向虚假的显著。协变量图Covariates则对应偏相关分析。比如你想排除灰质体积对相关结果的混淆可以在生成协变量图时使用VBM得到的调制后灰质密度图。JuSpace会在偏相关框架下把灰质密度的空间变异作为控制变量然后计算目标图与递质图谱的偏相关系数。这个功能在比较两个年龄跨度大的组别时几乎必备因为灰质萎缩本身与递质密度下降高度相关不控制的话说不清楚关联是来自递质系统还是脑结构。4. 空间相关的显著性问题为什么不能把体素当独立样本以及Spin Test在这里的作用4.1 空间自相关带来的自由度膨胀如果在做空间相关时直接把所有体素作为独立样本代入Pearson相关公式自由度会是几千甚至上万那么几乎任何微弱的空间共定位都会因为样本量巨大而变得高度显著。这明显是错的因为相邻体素的值绝对不是独立的——空间平滑、局部解剖一致性、数据插值都让相邻体素共享了大量信息。JuSpace提供的解决办法是基于置换检验的Null分布。具体来说它会在保持目标图空间结构不变的情况下对递质图谱或对目标图进行某种形式的空间重排重复数千次后得到相关系数的零分布然后考察真实相关系数落在零分布的哪个位置从而得到经验p值。这个过程绕开了参数检验对独立性的要求在思想上与permutation test一致。4.2 Spin Test在大球面上旋转数据的思路JuSpace的置换策略里一个重要的实现是Spin Test。这个思路最早来自脑连接领域Alexander-Bloch等人提出用球谐函数旋转皮层表面的脑图保持图谱自身的空间平滑结构只打破图谱与空间位置的对应关系。你可以想象一个地球仪你在球面上画了一套大陆海岸线的图案现在把这个图案整体绕球心旋转一个随机角度如果不看经纬度标签图案内部的相对结构一点都没变但它相对于球面上某个固定坐标比如某个脑区的对应关系已经被完全打乱了。把这个操作应用到脑图谱上就相当于生成一个空间结构相同、但与解剖位置不对应的零分布样本。每一次旋转都会重新计算一次目标图与该旋转后递质图谱的空间相关几千次旋转得到的就是一个纯空间结构驱动的零分布。JuSpace对皮层表面数据以及涉及皮层区域的图建议使用这种策略因为它能最大程度地保留空间自相关的程度。4.3 多重比较校正当你同时检验20个递质图谱时另一个容易被忽视的点是多重比较。JuSpace默认按照用户选择的递质图谱数量计算相关每个图谱都会得到一个相关值和p值。如果一次性选了20个图谱按0.05的阈值会期望有1个假阳性。这时候建议对p值做多重比较校正。常见的做法是FDRBenjamini-Hochberg校正——不需要编程把p值列表粘贴到R或者Python的statsmodels里就能做from statsmodels.stats.multitest import multipletests pvals [0.001, 0.02, 0.03, 0.11, 0.20, 0.35] fdr_p multipletests(pvals, methodfdr_bh)[1] print(fdr_p)更严格的做法是对所有图谱同时做置换检验然后用置换的极大值分布做校正即max-stat permutation。它可以控制family-wise error rate代价是计算量大得多。在论文里我一般建议至少汇报FDR校正后的结果因为这已经是当前领域的惯例。4.4 效应量比p值更值得关注我在实际审稿过程中发现很多人过分看重p值反而忽略了一个更重要的指标相关效应量。JuSpace输出结果里的r值通常在0到0.5这个区间浮动0.2左右就属于比较稳定的空间共定位了。看单张图和单张图谱的相关时r0.3在几千个置换后p0.001这是很亮眼的结果但在同一批分析里同时出现10个图谱、只有1个r0.3其余都是0.05这时候就算FDR后仍然显著也要警惕是不是20种图谱中恰好有一个因为平滑参数或来源差异碰上了。稳健的做法是把所有结果整理成表格按效应量排序看显著结果是否在生物学上构成一组有意义的系统而不是孤立的一条。5. 用JuSpace能讲出什么新故事三类最常见的研究设计5.1 第一类疾病组间差异图的分子解读这是JuSpace用得最多的一类设计。研究者在独立样本中计算患者组减去对照组的脑影像差异图如灰质体积、功能连接强度、低频振幅等然后把差异图喂给JuSpace看它是否在空间上与某种神经递质系统的分布重叠。举例来说一项关于精神分裂症的研究可能先得到患者相比对照在额叶-颞叶网络内部功能连接降低的差异图JuSpace分析显示该差异图与D2受体分布的空间相关显著r0.28, p_spin0.004。这就为精神分裂症的功能连接改变与多巴胺系统拓扑结构相关提供了一个定量证据链。要注意的是这个分析里的相关指空间共定位不代表个体水平的剂量-效应关系更不代表因果关系。论文里必须把这个边界讲清楚。这种研究设计我个人的体会是它非常适合用来做机制假说生成。它无法替代真正的PET上市后研究但能在大规模公开数据集如UK Biobank、ADNI、HCP里先打出候选信号为后续申请PET课题提供前期数据。5.2 第二类个体差异图与递质系统的关联如果你想研究某认知维度的个体差异与多巴胺系统的关系可以先把每个受试者的影像数据回归到行为指标上得到一张与行为相关的统计图比如工作记忆容量的全脑相关图再做JuSpace分析。HCP数据集的很多工作就是这个套路。这类设计里我特别强调要注意行为变量的分布如果存在离群受试者回归系数会被少数人主导全脑相关图会变得不稳定进而直接影响JuSpace的空间相关结果。建议先画散点图确定没有极端值再进入正式分析如果行为指标不是正态最好做一次秩变换。另外这类设计天然地存在一个天坑最终得到的个体差异图如果全脑只有两个很小的簇显著JuSpace跑出来的空间相关其实对这个稀疏度非常敏感。稀疏示意图与平滑递质图谱之间的空间相关值会比较低但置换检验仍然可能显著——这本身没错只是报告时建议同时说明目标图中显著体素的比例避免读者认为整个差异图都高度吻合多巴胺拓扑。5.3 第三类作为刺激靶点或手术靶点选择的辅助工具JuSpace还有一类比较出圈的应用用于神经调控靶点选择。rTMS、tDCS和DBS的靶点选择传统上主要基于解剖坐标和经验证据。如果你有一个初步的目标坐标区域可以把该区域作为种子点的静息态功能连接图或该区域在任务态下的激活图作为输入然后交给JuSpace评估它与哪种递质系统在空间上重合度最高。这能为为什么选择这个靶点提供更多分子层面的解释尤其适合需要在论文引言或讨论中说清楚靶点特异性的场景。当然我不会建议任何临床医生只凭JuSpace的结果来决定手术靶点但作为一个辅助验证工具它的价值是实打实的。尤其当你想在多项候选靶点里做选择时JuSpace能给出一个可量化的比较维度。6. 我在用JuSpace时实实在在踩过的坑以及绕坑建议6.1 平滑核导致伪显著怎么判断相关结果到底是真实信号还是平滑假象我先讲一个我自己的失败经历。某次我把一张6mm平滑的SPM脑激活图与5-HT2A图谱做空间相关结果r高达0.6p0.001。当时我一度以为发现了一个非常强的血清素-任务激活共定位。后来我拿同一张图做了多次下采样和平滑处理后发现只要目标图的空间平滑程度与递质图谱接近相关值就会显著上升把目标图平滑扩大到12mm后r几乎不变。这让我意识到空间相关的显著性本质上一部分由两张图的低频空间结构是否相似决定而平滑核正是塑造低频结构的关键因素。为了避免这个问题我现在的标准流程是先对目标图做与递质图谱一致的平滑处理然后同时跑原始分辨率目标图与平滑后目标图两组结果。如果两组的显著模式一致那比较可信如果只有平滑后显著论文里必须如实报告并谨慎解释。还有一个更隐秘的问题如果你输入的递质图谱和目标图分别来自不同的模板空间或mask边缘错位也会产生虚假相关。JuSpace官方提供的图谱mask是统一的但用户自己的目标图mask往往来自SPM的灰质先验边缘与图谱不一定吻合。运行前用spm_check_reg把所有图谱和目标图叠加起来肉眼看一遍边缘对应情况这30秒的花费远远少于出问题后排查的时间。6.2 图谱之间的间接依赖JuSpace里的递质图谱并非彼此独立5-HT1A和5-HT1B虽然受体亚型不同但表达区域有大面积重叠GABA-A/BZ图谱与全脑灰质分布的相关系数可能非常高。因此当你发现目标图与5-HT1A显著相关不要立刻说目标图特异地与5-HT1A系统相关因为没有控制其他图谱的共变。一个稳妥的做法是做偏相关把灰质密度图或者全脑平均摄取图作为协变量纳入模型看相关性是否稳定存在。如果加入灰质密度协变量后相关性跌到不显著那说明驱动信号的其实是目标图与灰质分布的关系不是递质特异性。6.3 统计结论的外推边界JuSpace结果能告诉我们的是在MNI空间里某张效应图与某种递质受体分布的空间共定位程度。它不告诉我们某个受试者个体身上这个效应是否与受体密度变化相关也不告诉我们药物干预这块脑区就一定会改变这条通路。生态学谬误是这类空间相关分析最容易犯的解读错误——群体水平的空间相关不能直接映射到个体水平的机制推断。写论文时一定要用spatially associated而不是causally linked这类表述。6.4 版本与运行时环境的建议JuSpace是在MATLAB中运行的对版本有一系列要求官方README有写明需要SPM12在路径上。我建议在Linux服务器上跑因为置换检验对CPU核数有较高需求JuSpace使用的是MATLAB Parallel Computing Toolbox的parfor。记得在调用前parpool开好工作池不然默认单线程跑1000次置换可能要几小时。另一点是内存如果你的目标图是1mm分辨率的大矩阵图谱重采样到相同网格后会消耗大量RAM。我一般统一重采样到2mm网格这既能跑得动又不会明显丢失空间信息。6.5 输出结果的归档JuSpace会生成一系列输出文件包括每个图谱的相关系数、p值、置换分布图等。我建议每次运行后把如下信息记录在分析日志里JuSpace版本号、所有输入图的具体路径和分辨率、重采样插值方法、平滑核大小、置换次数和相关方法、是否使用了协变量图、FDR校正版本。这样做的原因很实际——你半年后回来复现分析时如果没有这批元信息几乎不可能精确复现当时的参数组合。学术写作中把附录里的分析参数写清楚对审稿人来说也是一个极大的加分项。写在最后的一点体会做完几个项目之后我对JuSpace的定位有了更清晰的认识它不会直接证明某种递质系统导致了某群患者的影像异常但它能把这个影像异常在空间上更像哪种递质分布这个问题变成可检验的定量结果而且这个结果可以放进一篇论文的机制讨论里作为从宏观影像向分子层面延伸的重要佐证。对我个人来说最有价值的是它逼着我把空间共定位和分子机制两件事在写作上严格分界防止过度推断。如果你也正在处理一批积累了许久的脑图数据建议花一个下午把Demo跑通然后拿你自己的一张效应图做一次探索性分析——我猜你也会对这个工具能提供的解释力度感到意外。