利用信息化连通性来衡量的功能磁共振成像多体素穿越时空信息的同步崛起

摘要

Informational connectivity measures the correspondence between time courses of multi-voxel information across different brain regions. Multi-voxel pattern discriminability time series are extracted from regions and compared, revealing networks that are not identified in a typical functional connectivity approach.

摘要

人们现在认识到，条件相关的信息可以是分布式功能性磁共振成像（fMRI）大脑活动模式中存在，即使是单变量的激活水平相近的条件。多体素模式（MVP）的分析已被用于这些信息与巨大的成功解码。 FMRI研究人员还经常设法了解大脑区域在互联网络的交互方式，并使用功能连接（FC），以确定有相关的反应随着时间的推移地区。正如单变量分析中可以不敏感的信息的MVP，FC可能无法完全表征大脑网络，处理与特征的MVP签名的条件。这里介绍的方法，信息连通性（IC），可识别区域与MVP-辨别力与时间相关的变化，揭示连接是无法访问的FC。该方法可以是探索性的，用探照灯来确定种子CONNected区，或计划，预先选择的区域 - 的利益之间。该结果可以阐明处理的MVP相关条件的地区网络，可以击穿MVPA探照灯地图为独立的网络，或者可以在任务和患者群体进行比较。

引言

这里所描述的分析方法的目的是测量的基础上在其多体素信息的波动的大脑区域之间的连接。在功能磁共振成像（fMRI）分析技术的进展已经表明，大量的信息可以被包含在被分布在多个个体素^1-3的血液氧合水平依赖（BOLD）的活动模式。一组技术，是多元信息敏感-被称为多体素模式分析（MVPA） -已被用来证明条件可以有区别的最有价值球员，尽管有区分单变量的反应^1,2,4。标准的分析，该比较单变量的响应，可以是不敏感的这种多体素信息。

多个脑区当人类处理刺激和进行认知操作订婚了。功能连接（FC）是常用来调查个方法门这样的功能网络^5,6。在其最基本的形式，FC量化共激活，或同步性，不同的体素或区域之间。 FC已经被用于鉴定功能性连接大脑网络取得很大的成功。对于许多地区和条件，但是，单变量的反应并不反映BOLD活动中的所有可用信息。跟踪动态变化的单变量响应级别的FC技术可能缺乏敏感性，在多体素信息，共同波动。这里的分析方法介绍，信息连通性（IC，在最近的一篇论文⁷首次描述），桥梁MVPA和FC之间的差距，通过测量连接性是跨时多体素敏感信息的度量。虽然FC跟踪动态变化的单因素激活，IC音轨动态变化的MVP鉴别力 - 如何以及一个MVP的真实状况可以从（不正确）的替代品加以区别的措施。重要的是，在T他同样的方式，不同的区域可以显示单变量的反应相似水平的条件，尽管执行不同的计算（ 例如 ，视觉处理或行动计划，当一个人观看人造物体），不同的区域也可以有类似的（和同步）的MVP级别而他们加工条件不同鉴别力。最近的一项调查表明，IC可以揭示区域间的连接是无法检测用标准的FC方法^7。因此，研究人员可以使用的IC来探测大脑区域之间的互动为参与者有特征的分布式模式的条件或刺激作出反应。 IC是从最近的几个连接的应用程序，在有关分类结果^8,9检验单因素激活的波动显着。不同于这些方法，IC检测区域之间同步多体素模式鉴别力。

研究方案

1，准备功能磁共振成像数据

注：进行的功能磁共振成像扫描后，前处理使用功能磁共振成像大多数软件包开始本协议之前，可用的工具（但应避免或减少维护多体素模式空间平滑）收集的数据。一个合适的数据集的例子中的方法⁷的先前申请中描述。

1. 除去运动和平均脑白质信号从时间序列预处理的fMRI数据的创建与预测回归模型运动参数（横滚，俯仰，偏航，X，Y，Z）和平均脑白质信号。进行下面对所得的残差（ 即剩余的方差）的分析。
2. 导入生成的残差成分析软件包（ 如 MATLAB，Python的）。开源信息化连通性工具箱（http://www.informationalconnectivity.org）可以导入fMRI数据到MATLAB中。
3. Z评分​​EACħ体素的时间序列。
4. 独立数据集的时间点到独立集（'折叠'），如不同的扫描仪运行。注意：使用扫描仪的运行，确保褶皱，可以以其他方式难以保证（例如，依赖关系可以运行的时间点之间的前处理过程中产生）之间的独立性。运行可被组合在一起，以减少褶皱的数量（ 例如 ，偶数和奇数奔跑^2），虽然采用单次运行将给予更多的训练数据。
5. 创建与时间点通过生成的条件标签的向量是N久时间点相关联的条件标签的记录。
6. 向前移动的状态标签中的每个通过数次到重复（TRS）相当于5秒运行，为了解释事件和记录的fMRI信号之间的血流动力学滞后。

（2）选择和分析一个地区种子

1. 通过隔离解剖选择种子区域人区，在功能上局部区域或表现最出色的“信息脑地形图'射灯^10。
   注：步骤2.2至4.2以下可以由开源信息化连通MATLAB工具箱进行（ http://www.informationalconnectivity.org ）。
2. 比较每个时间点的最有价值球员为原型的MVP每个条件（注：这是因为在流行的相关性为基础的最近邻分类²用同样的方法）。 图1（上）给出了一个例子从收集真实的数据与会者认为，四种类型的人造物体块。
   1. 通过平均所有，但一折的每个条件的时间点计算出一个原型（平均值）最有价值球员为每个条件。这是'训练'数据为每个倍（ 例如 ，5倍2，均值的MVP是由褶皱1，3，4个时间点计算，和5）。
   2. 每一个相关时间点“;最有价值球员从训练数据中每个条件的均MVP。这将使每一个时间点的每个条件一个相关值（注：这里有最高相关性的条件将是受欢迎的相关基础MVPA分类²的预测）
   3. 费雪变换的r值z分数。
3. 量化'MVP的辨别力'为每个时间点：首先确定从2.2.3，表示该时间点的MVP和该时间点的状态的均值-MVP之间的关系的相关性，然后减去最高的剩余相关性（ 即 '相关的正确条件'减'有一个不正确的条件“最大相关）。其结果是，时间点的MVP鉴别力。一种替代（和有效的）的方法是减去的不正确的条件下的平均相关性。
   注：建议的方法具有直观的优势的GE与正面的鉴别力值时间点的相关​​性为基础的分类是正确的分类，而时间点与负的值不正确预测。所得到的值的一个例子示于图1（下部）。步骤至此为止被捕获在下面的公式。
   
   x是归一化的1×M个行向量，在时间点n×m个像素的激活值，Y是正确的（c）或不正确的（ⅰ）的条件下的平均训练模式的归一化的1×M个行向量时间点 n。该artanh功能应用费希尔Z变换。

3，计算时间序列的MVP区分度为每个探照灯

1. 进行分析探照灯^10：将三维C周围依次在每个体素（A'探照灯'）光泽。
2. 重复步骤2.2和2.3为每个探照灯，让每一个探照灯具有时间序列的MVP鉴别力值（每个时间点一个）。

4，计算信息化连接的种子和探照灯之间

1. 相关种子的MVP鉴别力的时间序列（2.3）使用Spearman秩相关各探照灯的鉴别力的时间序列（3.2）。由此产生的_R S值是种子和探照灯之间的IC。
2. 每个探照灯的IC值分配给探照灯的中心像素，并写出所产生的个体的大脑地图。

5，计算集团统计地图

1. 如果数据是不是已经在标准化的空间（ 例如 ，塔莱拉什或MNI），变换参与者的集成电路映射到相同的空间中。
2. 可选的平滑和个人＃39;探照灯地图。
3. 创建使用单程t-检验对每个探照灯的集成电路值是否显著大于零的一组统计图。

6，测试的意义

注：许多方法来确定的功能磁共振成像组地图统计学意义存在。作为置换检验可以用最少的假设确定的意义，而占平滑（因为每个置换组映射经过相同的处理）的数据集的水平，这个选项概述如下。

1. 对于每1,000排列，随机整个时间序列洗牌种子的MVP，鉴别力值。保持相邻的时间点与时间自相关性（如相同的块内的时间点）在一起（例如，通过混洗块，而不是相邻的TR）。
2. 计算个体的集成电路地图的每个置换（上面的步骤4）。
3. 产生1,000置换群图：RandomlŸ选择从每个参与者一个置换IC图，并在此随机设置（上面的步骤5）进行一组测试。
4. 阈值（ 例如 ，P <0.001），每置换组映射在所需的阈值，并从地图中提取的最大簇大小。
5. 排序产生的1,000最大类卷和识别该集群的大小在95 ^百分位（ 如 50 ^个最大的1000排列）。
6. 为了使群集校正为p <0.05适用于6.4中使用的阈值（ 例如 ，P <0.001）和最小群集大小从6.5到真正的（非重排的）集成电路组的映射。因为每一个置换图绘制在同一MVP的鉴别力值（以不同的顺序），这个地图的意义与突出偶然比预期的更同步鉴别力值区域。

结果

该IC结果现在可以使用研究者的首选功能磁共振成像分析软件包装上显示。 图2显示芯片的结果，从视觉上呈现男人的块计算的人造物体（在相关的出版物⁷全部细节）。

该IC分析是已知有关联的MVP情况尤为可贵：条件与特点的MVP，但没有在单因素反应的差异，更可能有集成电路和FC之间的区别（说明与被记录为参与者的数据观察不同类型的男人在图3中制造的物品）， 图4示出了与显著多体素信息可具有高的集成电路，但不能很好地在FC结果表示射灯。

。图1图案辨别力随时间的实例上： 从在着Haxby受试者等计算MVP识别性的基片^{（2001）2，}如在Coutanche与汤普森-希尔（2013）分析了^7。蓝线表示的时间点“最有价值球员，平均（'培训'）正确的类的模式之间的z得分的相关性。绿线代表了最有价值球员“的相关性有三个不正确的类。底部：模式鉴别力是正确的类和最高不正确的类的相关性之间的差异。时间点与正面图案辨别力的值将通过基于相关的分类器被正确地分类。图最初发表于Coutanche＆汤普森-希尔（2013年^）7。 请点击此处查看该图的放大版本。

图2实施例的连接图。 每行显示显著连接到一个种子（显示为蓝色）的区域。意义是由一组的t检验（P <0.001）与由置换检测最小的簇大小来确定。该IC结果显示使用AFNI ^{1 1}与FreeSurfer ^{1 2}生产表面地图。图是从Coutanche＆汤普森-希尔（2013年）修改^7。 请点击此处查看该图的放大版本。

数字3，同步MVP鉴别力相比，意味着激活MVP的鉴别力在两个地区与同步MVP的鉴别力（ 即信息连接）不同步的平均激活（ 即功能连接）的例子;数据来自从哈克斯比等人的主题^{。（2001）2，}如Coutanche＆汤普森希尔（2013年）分析^7。这些数据点采集，而受观察的人造物体，它是可区分由多体素模式，但并不意味着响应的可视化演示。 请点击此处查看该图的放大版本。

图4例的IC和LEF种子之间的FC值吨梭状回和整个大脑探照灯，信息性 ​​和功能性的连接强度（Z轴）显示种子和探照灯，对于每个探照灯的平均响应（x轴）和MVPA分类精度（y轴）到四之间类型的人造物体（几率= 25％）。探照灯与种子区域共享的体素被拆除。该IC图包括具有较高的分类性能，但较低的平均响应级别，这是不拿起一个典型的FC方式（通过在右图的左上角八分的差距可见）具有较强的连接探照灯的例子。图最初发表于Coutanche＆汤普森-希尔（2013年^）7。 请点击此处查看该图的放大版本。

讨论

信息连接有MVPA的灵敏度分布模式信息，并给出了通过连接的方式来研究区域之间的相互作用的能力。 MVPA和标准的单因素分析可以揭示每个不同区域的参与，有时他们的结果^{1 3}之间几乎没有重叠。正如预期的那样进行，吸引了这些分析方法的方法，集成电路和FC也给相辅相成的结果^7。是否采用集成电路的决定最终将取决于调查的条件和所提出的理论问题。设计考虑这种影响是否MVPA是对数据集进行的，也会影响芯片是否被使用。研究在设计时考虑到IC明确将要遵循建议的MVPA ^{1 4，}同时确保审判级别的数据可以从整个扫描的时间过程来提取。

当检查和报告集成电路的效果，重要的是，探照灯与种子重叠被去除，以避免圆形度。此外，如果直接比较集成电路和FC的效果，建议也比较基于均值活化探照灯，而不是仅仅一个体素的FC分析。这个附加的分析，可以确保结果之间的差异是因为在信号与噪声在探照灯与体素的水平差异不。

这里所描述的步骤主要侧重于探索性分析采用探照灯。值得注意的是，通过与区域的感兴趣更换探照灯，IC还可以比较选择一个我之前的区域。目前的辨别力指标 - 一个比较MVP的相关性为“真”状态的相关性为最大的交替条件 - 也是可以修改的。许多机器学习分类有预测的权重差异erent类，它可以很容易地更换这里进行（ 例如 ，跟踪一个分类器在一段时间的'信心'）的相关性比较。集成电路具有多种潜在的用途。以及作为一个主要的分析，调查信息网络，集成电路可以进行二次随访分析到MVPA探照灯。 MVPA探照灯地图是为了解哪些区域可以区分不同条件有价值的，但通常不会分解成不同的网络。该IC的方法可以帮助在这里，通过揭示这台探照灯有同步的鉴别力。最后，从不同的任务IC图可以为了了解任务的网络进行比较，患者可与对照组相比，以更好地理解多体素差额^{1 5}是如何体现在网络层面。

材料

Name	Company	Catalog Number	Comments
MATLAB software package	MathWorks
AFNI software package	NIMH