Summary
本文提出了一种研究人脑空间导航功能网络的综合方法。这种方法结合了大规模神经影像元分析数据库、静息态功能磁共振成像以及网络建模和图论技术。
Abstract
空间导航是一项复杂的功能，涉及多感官信息的整合和操作。使用不同的导航任务，在各种大脑区域（例如海马体、内嗅皮层和海马旁位置区域）的特定功能上取得了许多有希望的结果。最近，有人提出，涉及多个相互作用的大脑区域的非聚集网络过程可以更好地表征这种复杂功能的神经基础。本文提出了一种构建和分析人脑空间导航功能特异性网络的综合方法。简而言之，这种综合方法包括三个主要步骤：1）识别对空间导航很重要的大脑区域（节点定义）;2）估计这些区域中每一对之间的功能连通性，并构建连通性矩阵（网络构造）;3）研究所得网络的拓扑特性（例如，模块化和小世界性）（网络分析）。从网络的角度来看，所提出的方法可以帮助我们更好地理解我们的大脑如何在复杂和动态的环境中支持灵活的导航，并且揭示的网络拓扑特性还可以为指导临床实践中阿尔茨海默病的早期识别和诊断提供重要的生物标志物。
Introduction
功能特异性是人脑的基本组织原则，在塑造认知功能方面起着至关重要的作用1.功能特异性组织的异常可以反映标志性的认知障碍和主要脑部疾病（如自闭症和阿尔茨海默病）的相关病理基础 2,3。虽然传统理论和研究倾向于关注单个大脑区域，例如用于面部识别的梭形面部区域 （FFA）4 和用于场景处理的海马旁位置区域 （PPA）5，但越来越多的证据表明，复杂的认知功能，包括空间导航和语言，需要跨多个大脑区域的协调活动6.研究支持复杂认知功能的相互作用背后的机制是一个关键的科学问题，将有助于阐明大脑的功能结构和操作。本文以空间导航为例，提出了一种对人脑空间导航功能网络进行建模的综合方法。
空间导航是一种复杂的认知功能，它涉及多个认知组件的整合和操作，如视觉空间编码、记忆和决策7.通过功能性磁共振成像 （fMRI），许多研究在理解潜在的认知处理和神经机制方面取得了重大进展。例如，使用各种导航任务将特定功能与不同的大脑区域联系起来：场景处理与 PPA 特别相关，导航策略的转换与脾后皮层 （RSC） 相关8,9。这些研究为空间导航的神经基础提供了重要的见解。然而，导航是一种内部动态和多模态功能，单个区域的功能不足以解释通常观察到的空间导航10中的巨大个体差异。
随着基于fMRI的连接组学的出现，研究人员开始探索一些关键的大脑区域如何相互作用以支持空间导航。例如，已发现内嗅和后扣带回皮层之间的功能连接是高危阿尔茨海默病11 导航差异的基础。在另一项研究中，我们首次提出了一种网络方法，将连接组方法和几乎所有功能相关的区域（节点）整合在一起进行空间导航，结果表明该网络的拓扑特性与导航行为具有特定的关联12。这项研究为多个大脑区域如何相互作用以支持灵活导航行为的理论提供了新的见解10,13。
本工作展示了用于功能网络建模的集成方法的更新版本。简而言之，包括两项更新：1）虽然原始研究中定义的节点是基于更早和较小的数据库（55项研究，2,765次激活，2014年访问），但目前的定义是基于最新的数据库（77项研究，3,908次激活，2022年访问）;2）为了增加每个节点的功能同质性，除了原始的解剖学AAL（解剖学自动标记）图谱14外，我们还应用了一种新的脑包裹，它具有更精细的分辨率和更高的功能同质性（见下文）。我们预计这两项更新都将改进功能网络的建模。这个更新的协议提供了一个详细的过程，从网络的角度研究空间导航的神经基础，并有助于理解健康和疾病中导航行为的个体差异。类似的程序也可用于其他认知结构（例如，语言和记忆）的网络建模。
Subscription Required. Please recommend JoVE to your librarian.
Protocol
注意：此处使用的所有软件都显示在 材料表中。本研究中用于演示目的的数据来自人类连接组项目 （HCP： http://www. humanconnectome.org）15。所有实验程序均已获得华盛顿大学机构审查委员会（IRB）的批准。HCP 数据集中的成像数据是使用带有 32 通道头线圈的改进型 3T 西门子 Skyra 扫描仪获取的。其他图像采集参数详见早先的论文16。为演示下载了最少的预处理数据，该演示已完成以下预处理步骤：梯度失真校正、运动校正、场图预处理、空间失真校正、蒙特利尔神经学研究所 （MNI） 空间空间归一化、强度归一化和偏置场去除。也可以使用来自研究人员项目的静息态fMRI数据。
1. 数据预处理
- 检查数据质量并排除缺少重测数据和头部运动过度（平移 3 毫米和旋转 3°）的参与者。
注：5名受试者被移除，38名年轻人（22-35岁）被纳入主要分析。 - 在 MATLAB 中打开图理论网络分析 （GRETNA） 工具箱17 ，执行进一步的预处理步骤。单击“ FC Matrix Construction”批次。选择功能数据集的路径以加载 NIFTI 文档并执行以下步骤，如 图 1 中的管道选项所示：
- 通过双击“删除第一个图像”中的“要删除的时间点编号”并输入 10，删除前 10 个图像。
- 空间平滑（半高全宽 [FWHM] = [4 4 4]，在空间平滑中双击 FWHM （mm） 并输入 [4， 4， 4]）。
- 回归协变量。选择“白质信号”、“脑脊液信号”和“头部运动”为“TRUE”。根据实际体素大小选择合适的蒙版，例如，此处的 2 mm 蒙版，头部运动选择 Friston-24 参数。
- 时间过滤器。根据 MRI 扫描的重复时间输入 TR 的值（例如，此处为 720 ms），并通过双击 频段 （Hz） 并输入 [0.01 0.1] 来消除高频和低频噪声。
注意：有和没有全脑信号回归的结果如下所示。当使用未预处理的数据时，还建议使用完善的管道，例如 fMRI-prep18 和静息态 fMRI 数据处理助手 （DPARSF）19 。
图 1：Rs-fMRI 预处理和功能网络连接估计。 预处理的设置（去除前 10 张图像、使用 4 mm 的 FWHM 进行空间平滑、线性时间去趋势、回归白质信号、脑脊液 （CSF） 信号和头部运动的 24 个参数，过滤 0.01-0.1 HZ 的波段）以及与 fisher' Z 的静态相关性变换。缩写：Rs-fMRI = 静息态功能磁共振成像;FWHM = 半高全宽;CSF = 脑脊液。 请点击这里查看此图的较大版本.
2. 网络构建与分析
注：导航网络构建和分析的一般工作流程总结为三个主要步骤（图2）。
图 2：导航网络构建和分析的一般工作流程。 （A） 选择导航作为要在 Neurosynth 数据库中搜索的术语。（B） 可以生成激活坐标列表。（C） 使用 Neurosynth 中的函数运行荟萃分析，以获得多个大脑图谱。（D，E）通过结合元分析图谱和全脑包裹图谱 （AICHA），可以生成节点 （ROI）。（F） 利用由此产生的导航节点及其功能连通性构建导航网络（连通性估计和网络分析）。缩写：ROI = 感兴趣区域;AICHA = 同位素区域的内在连通性图集。 请点击这里查看此图的较大版本.
- 网络节点定义
- 通过在 Python 中键入命令下载最新的 Neurosynth 数据库 （neurosynth.org）20：
>import neurosynth as ns
>ns.dataset.download （path='./'， unpack = True）
注意：数据集存档 （'current_data.tar.gz'） 包含两个文件：“database.txt”和“features.txt”。这些包含来自神经影像学文章和荟萃分析标签的所有激活坐标，这些标签分别在该文章中高频率出现。 - 从database.txt生成新的 Dataset 实例，并通过键入以下命令向这些数据添加要素：
从 neurosynth.base.dataset 导入数据集>
> 数据集 = Dataset（'data/database.txt'）
> dataset.add_features（'数据/features.txt'）
- 通过键入以下命令，使用感兴趣的术语（即“导航”）运行荟萃分析：
> ids = dataset.get_ids_by_features （'navigation'， threshold=0.01）
> 马 = 元。MetaAnalysis（数据集，ids）
> 马.save_results（'.'， '导航'）
注意：荟萃分析得出了几张 NIFTI 格式的大脑图谱。误发现率 （FDR） 阈值为 0.01，用于控制误报率。在此步骤中需要归档知识，以确保将常见报告的区域包含在荟萃分析图中。类似的步骤可以应用于对其他认知功能（如语言和记忆）进行荟萃分析。 - 通过键入 FSL 中的命令，将元分析图谱和全脑包裹图谱合并，从而定义感兴趣的聚类：
>fslmaths navigation_0.01.nii.gz -bin navi_bin.nii.gz
>fslmaths navi_bin.nii.gz -mul AICHA/AAL.nii.gz navi_label_aicha/aal.nii.gz
>fslmaths navi_label_aicha/aal.nii.gz -thr n -uthr n 标签_n.nii.gz
>集群 -i 标签 _n.nii.gz -t 0.2 -o cluster_n.nii.gz
>fslmaths cluster_n.nii.gz -thr m -uthr m cluster_n_m.nii.gz
>fslmaths cluster_n_m.nii.gz -bin -mul x node_x.nii.gz
>fslmaths node_1.nii.gz -添加 ... -添加 node_x.nii.gz navi_AICHA/AAL_mask.nii.gz
注：这里使用了两个图集：AAL 和 AICHA。AAL 是原始研究中用于节点定义12 的图集。该图集是根据解剖剖图创建的。同位素区域内在连通性图谱 （AICHA）21 具有更精细的分辨率和更高的功能同质性。我们使用每个图集定义了感兴趣的区域。 - 在 Python 中键入脚本以检查地图中每个区域的大小：
>for i in np.arange（n）+1：
>____region_list.附录（i）
>____size1_list.append（np.sum（img_dat==i））
>____size2_list.append（np.sum（aicha_img_dat==i））
>____pct_list.append（np.sum（img_dat==i）/np.sum（aicha_img_dat==i））
注意：脚本中的整数 n 表示 AICHA 和 AAL 宗地中的区域总数（分别为 384 和 128）。为了避免虚假聚类的影响，建议可以删除尺寸相对较小的聚类（例如，100 个体素）。这里使用的AICHA图集是使用功能连通性数据生成的，每个区域都显示出自身功能时间活动的同质性。
- 通过在 Python 中键入命令下载最新的 Neurosynth 数据库 （neurosynth.org）20：
- 网络连接估计
注： GRETNA 工具箱用于连通性估计和网络分析。- 单击“ FC Matrix Construction”批次。通过选择功能数据集的路径来加载预处理的 rs-fMRI 数据。单击 静态关联 选项。将上一步获得的节点上传为图谱，计算每对区域的rs-fMRI信号的静态相关性，并将其转换为Fisher's z评分以提高正态性。
注意： 详细操作如 图 1 所示。每个参与者的 N × N （N 表示节点数）的导航网络矩阵将以.txt格式获得。 - 通过以下步骤获得正加权网络，如 图 3 所示。
- 单击 网络分析批次。将网络矩阵添加到 “Brain Connectivity Matrix”（大脑连接矩阵 ）窗口中，然后选择一个输出目录进行准备。
- 对于网络配置的流水线选项，在矩阵的符号中选择正，这会将函数连接矩阵中的负连接设置为 0，并消除歧义连接22。选择网络类型作为加权，以获取无向加权网络。
注意：除了加权网络之外，还可以对网络进行二值化，以创建用于后续分析的二进制网络（使用不同的方法），但加权网络通常被认为具有更高的可靠性23,24。
- 单击“ FC Matrix Construction”批次。通过选择功能数据集的路径来加载预处理的 rs-fMRI 数据。单击 静态关联 选项。将上一步获得的节点上传为图谱，计算每对区域的rs-fMRI信号的静态相关性，并将其转换为Fisher's z评分以提高正态性。
- 网络分析
- 将小世界、全局效率、聚类系数、最短路径长度、度中心性和局部效率添加到 GRETNA 网络度量分析管道中，如图 3 所示。
注意：小世界和全球效率是两个全球网络指标。具体而言，具有小世界的网络可以以相对较低的布线成本最大限度地提高信息传输效率。全球效率反映了运输网络中并行信息的传输效率。对于节点网络指标，度中心度衡量连接到节点的链路数。顾名思义，最短路径长度是测量集成度的基础。聚类系数表示节点的相邻节点相互关联的程度。局部效率是与节点及其邻居的通信效率（详细公式和用法在这些论文中显示）17,25。Brain Connectivity Toolbox （BCT）25 和其他工具箱也可用于计算网络指标。 - 在阈值方法中选择 网络稀疏 性，排除杂散连接的混杂效应，并输入一组 阈值序列 （即此处使用 0.05、0.1、0.15、0.2、0.25、0.3、0.35、0.4、0.45、0.5 ），根据统计结果进一步确定合适的阈值。
注： 边数与网络中最大边数与节点数样本数的比率称为稀疏性阈值。稀疏性阈值保证不同的个体具有相同数量的边。我们选择探索不同的阈值进行验证，这可以为未来研究选择最佳阈值提供有用的数据。 - 将随机网络编号 设置为 1,000，以使用马尔可夫布线算法26 生成随机网络。单击“运行”，在设置完所有步骤后，在 GRETNA 中运行管道。
注意：与真正的大脑网络类似，随机网络保持相同数量的节点、边和度分布。为了确定它们是否是相当非随机的拓扑结构，将它们与大脑网络进行比较。运行管道后，将获得每个阈值的一组网络指标分数，以便进行进一步的统计分析。 - 通过四个步骤确定网络中的最佳模块数量。
- 计算平均导航网络。单击 “指标比较 ”的浴池，然后选择“ 连接”。加载上面获得的网络矩阵，然后选择 “平均（功能）” 操作。选择输出方向以保留平均网络矩阵;有关详细信息，请参阅 图 4 。
- 使用 MATLAB 中的函数 spectralcluster 将从上述步骤获得的平均网络分为 2、3、4 和 5 个模块。
- 使用脚本 procrustes_alignment.m 对齐模块划分后，计算 REST 1 和 REST 2 中划分为同一模块的节点比例。使用比例作为模块分区的重复性指标。
- 选择具有最高重复性的模块数量。
- 将小世界、全局效率、聚类系数、最短路径长度、度中心性和局部效率添加到 GRETNA 网络度量分析管道中，如图 3 所示。
- 统计分析
注意：以下分析主要用于验证，在将本协议应用于个体变异研究时不需要。- 检查具有不同类型节点定义策略的两个网络之间这些网络指标的相似性（即，本研究中产生的新策略称为 NaviNet_AICHA 和 Kong 等人的早期策略称为 NaviNet_AAL）12。使用 MATLAB 中的函数 corrcoef 计算 Pearson 相关性，并对每个稀疏性阈值重复分析。
注意：提取网络指标后，可以进行他们感兴趣的任何统计分析。 - 使用 MATLAB27,28 中的函数 ICC 检查这些网络指标的重测可靠性，该函数实现了类内相关系数的计算。
注：原始未校正的 p 值在代表性结果部分报告。0.2 < ICC < 0.4 被解释为表示公平的重测信度，而 ICC > 0.4 被解释为中等到良好的重测信度29,30。阴性 ICC 分数设置为零，因为阴性 ICC 的存在毫无意义且难以解释31.
- 检查具有不同类型节点定义策略的两个网络之间这些网络指标的相似性（即，本研究中产生的新策略称为 NaviNet_AICHA 和 Kong 等人的早期策略称为 NaviNet_AAL）12。使用 MATLAB 中的函数 corrcoef 计算 Pearson 相关性，并对每个稀疏性阈值重复分析。
图 3：网络指标分析。 该分析定义了具有 10 个阈值的加权正网络。计算小字和效率两个全局网络度量，聚类系数、最短路径长度、效率和度中心度四个节点网络度量。 请点击这里查看此图的较大版本.
图4：平均导航网络的计算。 平均（函数）操作有助于计算所有参与者的平均网络。 请点击这里查看此图的较大版本.
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
导航网络
本研究通过结合最新的荟萃分析神经影像数据库和 AICHA 图谱，确定了 27 个与空间导航相关的大脑区域。这些区域包括内侧颞区和顶叶区，这些区域在导航神经影像学研究中很常见。这些区域的空间分布如图5A和图5C所示。作为比较，我们还在图5B和图5D中可视化了空间导航区域的早期定义。AAL地图集中的20个地区被纳入比较范围。这两组区域显示出很大的重叠。
图 5：模块化平均导航网络。 （A） REST1 中NaviNet_AICHA的模块化。（b） REST1 中NaviNet_AAL的模块化。（C，D）分别表示REST1中NaviNet_AICHA和NaviNet_AAL的模块化，分别回归了整个大脑信号。不同的节点颜色表示每个网络中标识的不同模块。NaviNet_AICHA和NaviNet_AAL显示了两个模块，其中包括内侧颞模块和顶叶模块。缩写：AICHA = 同位素区域的内在连通性图谱;AAL = 解剖学自动标记。 请点击这里查看此图的较大版本.
有趣的是，这两个网络显示出相似的群落分布（图5）。具体而言，模块化和可重复性分析显示NaviNet_AICHA和NaviNet_AAL都有两个模块（一个腹侧模块包括内侧颞区，一个背侧模块包括顶叶区域）（表1）。NaviNet_AICHA 和 NaviNet_AAL 的腹侧和背侧模块相似，尽管前者的节点数量更大，因为 AICHAI 图谱的大脑包裹更精细。这些结果与预处理过程中用于处理全局信号的策略无关（图5）。此外，在REST2数据集中也观察到类似的群落分布（补充图S1）。
两个导航网络拓扑特性的相似性
接下来，我们检查了两个网络之间每个网络度量的相似性。相似性分析的目的有两个：（1）评估使用不同定义策略时结果的泛化性，以及（2）确定网络分析的最佳网络阈值。
总的来说，除聚类系数外，六个指标中有五个指标与网络分析中使用的大多数网络稀疏性阈值之间存在显著相关性（图 6）。除了平均节点度外，所有指标的相似度值都随着稀疏性阈值的增加而快速增加，该阈值与任一阈值都显示出出色的相似度值。小世界指标在 0.30 和 0.40 之间的阈值处显示最高的相似度，而其他指标也显示最高的相似度。这些结果表明，网络级分析可以反映稳定的个体差异，而与节点定义选择无关，并且 0.30-0.40 的稀疏阈值将导致导航网络分析具有更好的泛化性。请参阅 补充图 S2 ，了解与 REST2 的更多相似之处。
图 6：两个网络拓扑特性的相似性。 显示了 REST 1 （A） 没有和 （B） 回归全脑信号的结果。y 轴上的 Pearson 相关系数表示两个网络拓扑特性的相似性。稀疏性阈值范围为 0.05 到 0.5。图中的星号表示显著性级别，其中 *p < 0.05， **p < 0.01， ***p < 0.001。 请点击这里查看此图的较大版本.
重测信度
我们还评估了导航网络拓扑测量的重测信度。在计算这些网络度量值时，使用了介于 0.05 和 0.50 之间的各种稀疏阈值，以消除网络中潜在的虚假连接（详见 补充表 S1 ）。在这里，鉴于上述相似性结果，我们主要报告了阈值为0.40的结果。大多数网络指标在网络NaviNet_AICHA和NaviNet_AAL方面都显示出一般到良好的可靠性（ICC > 0.2），而NaviNet_AAL的可靠性相对高于NaviNet_AICHA。此外，我们发现在fMRI数据预处理中包括全局信号回归可以带来更高的可靠性（图7）。NaviNet_AAL网络中的聚类系数、最短路径长度和小世界表现出最高的重测信度，而聚类系数和小世界NaviNet_AICHA的重测信度也高于其他度量。这些结果表明，聚类系数和小世界是这些指标中最可靠的。
图 7：导航网络拓扑特性的重测信度。 （A） 在不回归整个大脑信号的情况下对数据进行重测信度。（B）全脑信号回归数据的重测信度。缩写：ICC = 类内相关系数;Cc = 聚类系数;Lp = 最短路径长度;Sw = 小世界;Nd = 平均淋巴结度;Eg = 全局效率;Eloc = 局部效率;AICHA = 同位素区域的内在连通性图谱;AAL = 解剖学自动标记;cNGS = 没有回归整个大脑信号的数据;cWGS = 回归整个大脑信号。 请点击这里查看此图的较大版本.
| 模块编号 | |||||
| 2 | 3 | 4 | 5 | ||
| NaviNet_AICHA | cNGS系列 | 1 | 1 | 0.96 | 0.67 |
| cWGS系列 | 1 | 0.96 | 0.78 | 0.89 | |
| NaviNet_AAL | cNGS系列 | 1 | 0.95 | 0.95 | 0.65 |
| cWGS系列 | 1 | 0.95 | 0.95 | 0.95 | |
表 1：REST 1 和 REST 2 之间模块分区的可重复性。 第一行表示模块数。cNGS 表示 rs-fMRI 而不回归全局信号，cWGS 表示 rs-fMRI 回归全局信号。数字越大表示重复性越高，本文选择了两个模块进行NaviNet_AICHA和NaviNet_AAL。缩写：AICHA = 同位素区域的内在连通性图谱;AAL = 解剖学自动标记。
补充图S1：REST 2中的模块化平均导航网络。 （甲，乙）REST 2 数据集中NaviNet_AICHA和NaviNet_AAL的模块化，而不会回归整个大脑信号。（三、四）REST 2 数据集中 NaviNet_AICHA 和 NaviNet_AAL 的模块化与全脑信号回归。不同的节点颜色表示每个网络中标识的不同模块。这两个模块都以 NaviNet_AICHA 和 NaviNet_AAL 显示。 请点击这里下载此文件。
补充图S2：REST 2中两个网络拓扑特性的相似性。 显示了没有/有回归全脑信号的 REST 2 的结果（分别为 A 和 B）。y 轴上的 Pearson 相关系数表示两个网络拓扑特性的相似性。稀疏性阈值范围为 0.05 到 0.5。图中的星号表示显著性水平，其中 *p < 0.05， **p < 0.01， ***p < 0.001。 请点击这里下载此文件。
附表S1：不同稀疏阈值导航网络拓扑特性的重测信度。 这些值表示具有不同稀疏性阈值的NaviNet_AICHA和NaviNet_AAL的类内相关系数。缩写：Cc = 聚类系数;Lp = 最短路径长度;Sw = 小 woldness;Nd = 平均淋巴结度;Eg = 全局效率;Eloc = 局部效率;cNGS = 没有回归整个大脑信号的数据;cWGS = 回归整个大脑信号。 请点击这里下载此文件。
Subscription Required. Please recommend JoVE to your librarian.
Discussion
网络神经科学有望帮助理解大脑网络如何支持人类认知功能32.该协议展示了一种研究人脑空间导航功能网络的综合方法，这也可以激发其他认知结构（例如语言）的网络建模。
该方法包括三个主要步骤：节点定义、网络构建和网络分析。虽然网络构建和网络分析与一般的全脑网络研究相同，但节点定义是该协议最关键的一步。此步骤利用与空间导航相关的功能激活的大规模荟萃分析来定位导航行为最重要的大脑区域。因此，我们可以对功能有意义的网络进行建模，这有助于从网络的角度理解复杂处理的神经基础。请注意，前额叶区域在节点定义结果中缺失，而越来越多的导航研究表明这些区域具有关键作用33.这可能是由于数据库内导航相关研究中缺乏这些区域的激活，导致meta分析的数据有限。当有更多的数据可用于定位这些与导航相关的前额叶区域时，在未来的研究中研究它们在导航网络中的作用将是一个有趣的问题。当可以定位单个大脑区域时，研究人员还可以应用该协议来研究其他认知功能。需要现场知识来识别感兴趣的区域，以确保与特定功能的密切关联。
在该协议中，我们专注于空间导航网络，并展示了空间导航研究中报告的各种大脑区域的高覆盖率。鉴于缺乏普遍同意的支持空间导航的大脑区域定义，演示使用了两组区域，一组通过结合最大的荟萃分析和 AICHA 图谱生成，另一组通过 AAL 图谱生成。基于两个定义的网络拓扑属性总体上表现出较高的相似性，在一定程度上支持了功能特定网络建模的有效性。
我们注意到，相似性强度随着网络分析中使用的稀疏性阈值的增加而增加，结果表明，0.30-0.40 的稀疏性阈值将是一个合适的选择，因为所有网络指标都显示出与这些阈值的最高相似性。有了这样的阈值，网络指标也显示出良好到良好的重测可靠性，特别是在数据预处理中包含全局信号回归的情况下，对于最短路径长度和小世界性。这些结果在很大程度上支持在个体差异和相关脑部疾病研究中使用这些指标。
由于缺乏适当的行为数据，我们无法在该协议中呈现网络指标与空间导航的行为相关性。基于之前对导航相关区域功能连接指标的大脑行为关联的一些研究12,34，我们预计使用该协议的网络建模将显示出与空间导航的特定关联。需要具有神经影像学和行为数据的大规模样本来进一步研究这些关联。此外，虽然重测信度结果不是很高，但强度与先前的 fMRI 研究报告的强度相当35.
未来的研究可以应用该协议，从网络角度更好地理解空间导航的神经基础，并探索其在人类中的变化。例如，研究人员可以利用该协议来研究导航网络的发展和老化轨迹，在临床实践中，网络特性为指导阿尔茨海默病等脑部疾病的早期识别和诊断提供了重要的生物标志物。此外，未来的研究还可以应用类似的协议来构建其他认知结构的网络模型。
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
作者声明不存在利益冲突。
Acknowledgments
孔祥珍获得国家自然科学基金（32171031）、科技创新2030重大专项（2021ZD0200409）、中央高校基本科研专项（2021XZZX006）、浙江大学信息技术中心等项目资助。
Materials
| Name | Company | Catalog Number | Comments |
| Brain connectivity toolbox (BCT) | Mikail Rubinov & Olaf Sporns | 2019 | The Brain Connectivity Toolbox (brain-connectivity-toolbox.net) is a MATLAB toolbox for complex-network (graph) analysis of structural and functional brain-connectivity data sets. |
| GRETNA | Jinhui Wang et al. | 2 | GRETNA is a graph theoretical network analysis toolbox which allows researchers to perform comprehensive analysis on the topology of brain connectome by integrating the most of network measures studied in current neuroscience field. |
| MATLAB | MathWorks | 2021a | MATLAB is a programming and numeric computing platform used by millions of engineers and scientists to analyze data, develop algorithms, and create models. |
| Python | Guido van Rossum et al. | 3.8.6 | Python is a programming language that lets you work more quickly and integrate your systems more effectively. |
| Statistical Parametric Mapping (SPM) | Karl Friston et.al | 12 | Statistical Parametric Mapping refers to the construction and assessment of spatially extended statistical processes used to test hypotheses about functional imaging data. |
References
- Kanwisher, N. Functional specificity in the human brain: a window into the functional architecture of the mind. Proceedings of the National Academy of Sciences of the United States of America. 107 (25), 11163-11170 (2010).
- Coughlan, G., Laczo, J., Hort, J., Minihane, A. M., Hornberger, M. Spatial navigation deficits - overlooked cognitive marker for preclinical Alzheimer disease. Nature Reviews Neurology. 14 (8), 496-506 (2018).
- Gilbert, S. J., Meuwese, J. D., Towgood, K. J., Frith, C. D., Burgess, P. W. Abnormal functional specialization within medial prefrontal cortex in high-functioning autism: a multi-voxel similarity analysis. Brain. 132 (4), 869-878 (2009).
- Kanwisher, N., McDermott, J., Chun, M. M. The fusiform face area: a module in human extrastriate cortex specialized for face perception. Journal of Neuroscience. 17 (11), 4302-4311 (1997).
- Epstein, R., Harris, A., Stanley, D., Kanwisher, N. The parahippocampal place area: recognition, navigation, or encoding. Neuron. 23 (1), 115-125 (1999).
- Epstein, R. A., Patai, E. Z., Julian, J. B., Spiers, H. J. The cognitive map in humans: spatial navigation and beyond. Nature Neuroscience. 20 (11), 1504-1513 (2017).
- Baumann, O., Mattingley, J. B. Extrahippocampal contributions to spatial navigation in humans: A review of the neuroimaging evidence. Hippocampus. 31 (7), 640-657 (2021).
- Auger, S. D., Mullally, S. L., Maguire, E. A. Retrosplenial cortex codes for permanent landmarks. PloS One. 7 (8), e43620 (2012).
- Zhen, Z., et al. Quantifying the variability of scene-selective regions: Interindividual, interhemispheric, and sex differences. Human Brain Mapping. 38 (4), 2260-2275 (2017).
- Ekstrom, A. D., Huffman, D. J., Starrett, M. Interacting networks of brain regions underlie human spatial navigation: a review and novel synthesis of the literature. Journal of Neurophysiology. 118 (6), 3328-3344 (2017).
- Coughlan, G., et al. Functional connectivity between the entorhinal and posterior cingulate cortices underpins navigation discrepancies in at-risk Alzheimer's disease. Neurobiology of Aging. 90, 110-118 (2020).
- Kong, X. Z., et al. Human navigation network: the intrinsic functional organization and behavioral relevance. Brain Structure & Function. 222 (2), 749-764 (2017).
- Weisberg, S. M., Ekstrom, A. D. Hippocampal volume and navigational ability: The map (ping) is not to scale. Neuroscience & Biobehavioral Reviews. 126, 102-112 (2021).
- Tzourio-Mazoyer, N., et al. Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage. 15 (1), 273-289 (2002).
- Van Essen, D. C., et al. The WU-Minn human connectome project: an overview. Neuroimage. 80, 62-79 (2013).
- Smith, S. M., et al. Resting-state fMRI in the human connectome project. Neuroimage. 80, 144-168 (2013).
- Wang, L., et al. GRETNA: a graph theoretical network analysis toolbox for imaging connectomics. Frontiers in Human Neuroscience. 9, 386 (2015).
- Esteban, O., et al. fMRIPrep: a robust preprocessing pipeline for functional MRI. Nature Methods. 16 (1), 111-116 (2019).
- Yan, C., Zang, Y. DPARSF: a MATLAB toolbox for" pipeline" data analysis of resting-state fMRI. Frontiers in Systems Neuroscience. 4, 13 (2010).
- Yarkoni, T., Poldrack, R. A., Nichols, T. E., Van Essen, D. C., Wager, T. D. Large-scale automated synthesis of human functional neuroimaging data. Nature Methods. 8 (8), 665-670 (2011).
- Joliot, M., et al. AICHA: An atlas of intrinsic connectivity of homotopic areas. Journal of Neuroscience Methods. 254, 46-59 (2015).
- Murphy, K., Birn, R. M., Handwerker, D. A., Jones, T. B., Bandettini, P. A. The impact of global signal regression on resting state correlations: are anti-correlated networks introduced. Neuroimage. 44 (3), 893-905 (2009).
- Fox, M. D., Zhang, D., Snyder, A. Z., Raichle, M. E. The global signal and observed anticorrelated resting state brain networks. Journal of Neurophysiology. 101 (6), 3270-3283 (2009).
- Xiang, J., et al. Graph-based network analysis of resting-state fMRI: test-retest reliability of binarized and weighted networks. Brain Imaging and Behavior. 14, 1361-1372 (2020).
- Rubinov, M., Sporns, O. Complex network measures of brain connectivity: uses and interpretations. Neuroimage. 52 (3), 1059-1069 (2010).
- Maslov, S., Sneppen, K. Specificity and stability in topology of protein networks. Science. 296 (5569), 910-913 (2002).
- Shrout, P. E., Fleiss, J. L. Intraclass correlations: uses in assessing rater reliability. Psychological Bulletin. 86 (2), 420 (1979).
- McGraw, K. O., Wong, S. P. Forming inferences about some intraclass correlation coefficients. Psychological Methods. 1 (1), 30 (1996).
- Andellini, M., Cannatà, V., Gazzellini, S., Bernardi, B., Napolitano, A. Test-retest reliability of graph metrics of resting state MRI functional brain networks: A review. Journal of Neuroscience Methods. 253, 183-192 (2015).
- Cao, H., et al. Test-retest reliability of fMRI-based graph theoretical properties during working memory, emotion processing, and resting state. Neuroimage. 84, 888-900 (2014).
- Rousson, V., Gasser, T., Seifert, B. Assessing intrarater, interrater and test-retest reliability of continuous measurements. Statistics in medicine. 21 (22), 3431-3446 (2002).
- Bullmore, E., Sporns, O. Complex brain networks: graph theoretical analysis of structural and functional systems. Nature Reviews Neuroscience. 10 (3), 186-198 (2009).
- Patai, E. Z., Spiers, H. J. The versatile wayfinder: prefrontal contributions to spatial navigation. Trends in Cognitive Sciences. 25 (6), 520-533 (2021).
- Wegman, J., Janzen, G. Neural encoding of objects relevant for navigation and resting state correlations with navigational ability. Journal of Cognitive Neuroscience. 23 (12), 3841-3854 (2011).
- Braun, U., et al. Test-retest reliability of resting-state connectivity network characteristics using fMRI and graph theoretical measures. Neuroimage. 59 (2), 1404-1412 (2012).
