鉴定的使用神经影像学数据的与疾病-相关的的空间协方差模式

Summary

多元技术，包括主成分分析（PCA）已经脑功能图像用于识别签名区域变化模式。我们已经开发了一种算法来识别重复的网络神经退行性疾病的诊断，评估病情恶化，患者人群的治疗效果和客观评价的生物标志物。

Abstract

的缩放的子配置文件模型（SSM） ^显示1-4，是一个多变量的基于PCA-的算法，该算法识别及确定其主要中的变异来源患者组和对照组群脑的图像数据在拒绝较小的的组件（ 图1）的同时。直接向应用于的稳态-是一张专属图片状态多模态图片包含的体素的-由-的体素的的协方差数据的，整个一个的的组的影像设定可以降低到一个几个的显着的线性独立的协方差图案和相应的的主题的分数。每个模式，都被称为，的一组不变量的子配置文件（GIS）的，是一个正交的的的本金代表一个空间上分布的的的网络的功能上相互关联的的的的大脑区域的组件，它。大型的全球性的均值的标量的的影响，可以掩盖较小的网络-特定的的贡献都由所固有的的对数转换中移除，，，并的意思的数据的定心^2,5,6。的的主题表达每个的这些模式给一个变量的表示的程度通过一个简单的的可以与独立的CLIN相关联的的的标量得分，的iCal描述符或心理测量学的描述符^7,8。通过使用的主题的分数（ 即图案表达式的值）的因素logistic回归分析，线性相关系数可以推导出，以结合起来多个的主要组件集成到与改进的的从健康的的控制科目^5,6的患者的歧视的单一的疾病-相关的的空间的协方差的的图案， 即的复合网络的。使用自举的重采样技术^9时，可以执行在的推导集的之内的交叉验证的。正向的验证是很容易证实由直接得分的评价的的派生的的图案在的前瞻性的数据集^10中 。可以用来一旦通过验证，疾病关连的的的型态，，以得分从个人的的例相对于到一个固定的的参考样品进行患者，往往是一套的"这是使用的的名健康受试者"（与这个疾病组）在原来的柄的推导^第11的 。这些标准化处理的值可以在反过来又被用来，以协助在的差分诊断^{12,13中，，，}并，以评估疾病的的进展和治疗效应在网络级别^7,14-16。我们提出了一种的应用程序此的方法来的帕金森氏病患者和正常对照组，使用我们的的的在-的房子的软件的情况，，以推导出一个的特征的协方差图案疾病的的生物标志物的的FDG PET显像的数据的的例子。

Introduction

神经退行性紊乱的已被广泛研究，本地化和量化异常的的大脑代谢作为以及作为非-推理的的那项研究区域相互作用^】17的方法的的的使用技术。 （，如本金成分分析^{（PCA）1,2,4,18}和独立分量分析^{（ICA）19,20，}如以及作为监督的的技术（，如偏至少的广场^（PLS）21和序数的发展趋势规范的））的以数据-驱动的的的的多变量的分析战略变元分析（ORT / ^CVA）22不只可以揭示的特征的的的图案或的相互关联的活动的"网络"。的多变量的的的程序，尤其是缩放后的的子配置文件模型（SSM）的的基础知识^1,2,4-6,18先前已经描述在朱庇特的^3中 。这家基于PCA-的做法被最初开发它来审查异常的功能性的协方差大脑区域之间的的关系，在的脑的血流量和新陈代谢ACQ的的稳定-状态的单一的的磁碟区影像中uired在静息状态的（如PET和SPECT），表现出高的的信号-，以-噪声的特性的网站的详细的方式，。疾病-特定的的的的SSM图案是反映了整体的的差异在的区域的地形中在予正常的科目^7,16相比，，并可以反映一个单一的网络过程或的同化的几个复杂的的异常的功能^23的患者中的成像生物标志物，。代谢的协方差图案大脑网络被相关联的与，可以区分正常的组和疾病的组之间的的，，，并提供（与临床的疾病的严重程度的评级吻合）的基于网络-的措施，，的的表达式的值（学科的分数）。通常情况下，这样的模式的科目比分为会增加与疾病进展，，，，并可能甚至可以被表达的在症状，发病时间之前^14,24。事实上，疾病关连的的的网络的生物标志物已被被其特征在于为神经退行性疾病的（，如帕金森氏病的^10条 （PD）的），亨廷顿舞蹈症（Huntington的病^25（"房署"），并阿耳茨海默氏的疾病的^8）的 "（AD）的。更重要的是，与疾病-相关的的的代谢拓朴图已经也被确定为非典型的帕金森病的的运动障碍性疾病，，如多个系统萎缩（，MSA）和渐进式的性核上性上性麻痹（PSP）的。这些模式已被用来在演唱会为个人的鉴别诊断中具有临床上相似的"，看看-的一致好评"的证候^12,13,26。

在此相反，典型的功能磁共振成像以体素的-基础的的的单变量的方法评估在孤立的的大脑集群在患者和对照组之间的的差异的重要意义。更多最近的是，的方法已被开发来衡量功能性的之间的连通性各种不同的定义的的大脑区域^{27日至29日} 。功能连通性了这个定义评论的被受限制，以〔人〕遭受和地区的，特定的的的相互作用和偏离从原来的的SSM /主成分分析的概念，是指，以的的内在的空间上分布的的的大脑网络区域^1,2,23,30的横截面的互联互通。要他们的的的优势，MRI检查平台一个重新很容易地安装，广泛可用的，非侵入性和通常会需要更短的的的传统的放射性示踪剂（如PET或SPECT）导致在一个在最近的文献的描述的潜在的的方法的热潮中的的成像方式的扫描的时间比。然而，由此产生的的依赖于时间-的的fMRI的信号提供了局部神经的活性^31,32的的间接措施的。的一般复杂的的分析受雇于的算法已被由的大型的大小（的的数据集），的生理的的噪声fMRI信号中固有的的，如以及作为的高的的的的可变性在的大脑活动中，之间存在着科目和地区的^19,23的限制。虽然有趣的的了解关于中的大脑组织可以的详细资讯，请参阅被方面的推断从功能磁共振成像的"网络"的属性，他们还没有一直足够的稳定，以被用来作为可靠的疾病的的生物标志物。此外，本因此而产生之的的网络拓朴图是不一定是相等于到使用已建立的的功能成像（例如作为SSM范围/ PCA）的方法的那些已识别的。对于THë键的大部分的的一部分，严谨的的跨-验证的的由此产​​生的的fMRI拓朴图的一直缺乏与在从单一的情况下，的的的前瞻性的扫描数据中的的派生图案的成功的向前的应用程序的的几个例子。

的PCA协方差分析的一个优点就在于在其能力，以找出的最显着的的的的数据变异来源在的第一个的的的少数几个主组件，，，但它是无效的，如果的突出的特征向量代表，而比实际的的的的内在的的网络响应的的的随机噪声因素的。通过选择只有成为第一个的几年的的特征向量，并。限制，以那些，显示在耐心相对于正常的的控制分数线的中的显着的的差异的的的中，我们大大降低其他的噪声元素的的影响力。然而，对于的基本的此处描述的方法中，这些措施的，可能会无法是足够的，以产生在一个典型的的功能性磁振造影中资料集与在下面的描述的的方式的的例外状况的稳健性估计量。

因此，，因为的的的区域的克的稳定的直接的的恋爱关系lucose的的代谢和突触神经的的活性^33，这种方法已已应用于的休息时间状态FDG PET显像的数据的分析主要是为了。然而，给定的，脑血流量（CBF）是紧密耦合，以代谢在的静息状态下^{10,11,34，SPECT [35,36]}和更最近的动脉自旋贴标（海拨）MRI灌注成像方法^37,38中的活性，已被用来，以评估在个别情况下的异常的的代谢活性。也就是说说，与静息状态下的fMRI（rsfMRI）的可靠的的空间的协方差图案的推导是如前面所指出的不直截了当的^31,32。即便如此的是，初步从PD患者和对照组受试者的的的rsfMRI数据的之基于SSM /对主成分分析（PCA）曾透露的一些之间的的地形的同源性与疾病-相关的的图案确定了，使用的两个方式，的BOLD功能磁共振成像^39,40的 PET和调幅的低-频率的波动的影响（ALFF）的振幅。最后，我们还注意到，这种做法已成功地应用于在基于体素的的形态中以及两者混合的起始点，（VBM）的结构性的MRI数据^41,42，揭示了与与年龄-相关的的的的的体积损失和在在的相同的科目^43中的 VBM的和ASL的图案的进一步的比较的中相关联的的鲜明的的空间的协方差图案。 SSM / PCA的空间使用不同的的分析的方法和成像平台的的的的协方差拓朴图和类似的确定的大脑网络的的之间的关系是正在进行的调查的一个话题。

Protocol

1。数据的收集和预处理

1. 该SSM / PCA的方法可以被应用于来自的各种来源和方式的获得的，以单一的的磁碟区影像的。具体而言，对于上的-现场的PET显像对的新陈代谢中，准备一个合适的的的放射性核素的示踪剂例如为^[18 F]-氟脱氧葡萄糖脱氧葡萄糖（FDG），，和管理，以每一个病人。的患者通常是在休息扫描的与眼睛睁得大大，在至少12小时一个快速的的后，，关的药物。
2. 扫描个人或一组的评估的每个主题为。如需"格局的推导"中，扫描对性别问题-和年龄-匹配的的的患者和对照组的的的一个相等数量的。
3. 中传输一些数据，到一台工作站，，，，并将其转换在你的的平台上以供分析为适当的格式。我们的Windows PC为基础的:MATLAB的分析软件（scanvp，ssm_pca， www.feinsteinneuroscience.org ）的需要ANALYZE或nifti格式的的图像（梅奥诊所的，罗切斯特，MN的）。它提供了了一个转换例程的Ţ:Ø变换GE垫款（密尔沃基，WI的，USA）扫描仪和其他的格式的的图像，以ANALYZE的格式。
4. 每个主题的形象正常化到一个共同的的脑立体定位的空间（ 例如:蒙特利尔神经学研究所[MNI的]），使用的一个标准的的（，如的统计的参数化映射（SPM））的的神经影像学检查软件程序包（ http://www.fil.ion.ucl.ac.uk/spm中 ），所以，，有是一个one-to -酮的的对应的子民与子民之间的体素值（ 图2）。美纹纸，以限制分析，以的的灰质领域（ 图3）和日志转型顷中描述的接下来的步骤中。

2。执行多变量的SSM / PCA的

进行多变量的基于SSM / PCA（ 图4）的操作可以执行由外部软件。的步骤下文分项列示反映执行的的程序为在大多数的一部分下会自动由我们的在-的房子例程（ 图5）（ 扫描VP / ssmpca， www.feinsteinneuroscience.org也显示有空房的作为一个SPM工具箱ssm_pca）。

1. 使用一个可用的的0/1的图片此的的面膜数据以消除不必要的的领域的的体素的空间（，如脑白质及白质脑室的）。使用该ssmpca例程，用户将被给定的一个选项，以中选择的默认情况下或其他的外部掩码（ 图3），，或，以不得不的程式自动地建立一个面具，通过去除较低比一个选定的固定的的的每个主题的数据的最大中的百分之的值。然后，个别的主题口罩都相乘，以确定一个复合的面具。在的掩模之内的的领域定义了常见的一种非-零体元空间，用于分析。
2. 三维上的汇率换算每个检者的屏蔽的图像，体素数据，以一个连续的的的行向量，通过追加从连续的的平面的的表顺序扫描，线（ 图4a）。形成一组数据矩阵（D），以便，每个主题的数据对应到一个特定的R嗷嗷:的基体中。每个Column，然后选择代表一个特定的的的跨科目的体素的。理想的情况下，对于生物标志物的的推导中时，矩阵将被组成的一个的正常的控制的主题的的数据和疾病的的主题数据的的的行的相等数目的。
3. 变换每个的数据条目，以对数的形式。
4. 中心的数据通过中减去的每一行平均水平或主题的的矩阵的意思是从的行元素中的。的所有为中心的的行的的平均代表的特性基是什么意思日志映像的映像被称为一组均值设定档（ 的GMP，图4a）。中减去的的级的GMP的元素从相应的的矩阵column元素的顷的的列装置​​，，装置。的双重的为中心的矩阵的每一行都代表一个'残余的'的的形象被称为一个主题的残余设定档（SRP），，其代表从双方的主体s及体素的V族装置（ 图4b）的偏差的元素。
   SRP _SV = LOGD _SV -的意思_{S -} GMP _v
5. 构建各方，不对译文的准确性，ECT-受试者由-协方差矩阵 C的的复合的双通过计算的非-每个主题的线对 i 的归一化之间的协方差时的数据中心化的矩阵中，J的SRP矩阵的行的作为一个产品的超过所有的体素v总结出的相应的的体素的元素进行评估（ 图4b） 。
   _{ÇIJ =ΣV（SRP IV}所述SRP _JV）
6. 套用PCA到的主题-由-须受到的的协方差矩阵 C。该的结果将是与相关联的的本征值的的主题的得分的特征向量了一套。赋权的每个向量，通过乘以由的平方的根目录中公司相应的的的本征值。上表示的得分特征向量的提取的一组被由所得到的矩阵 S的列[^]在图4c的中 。
7. 统为可以在矩阵的特征值的相同的一组的体素的的特征向量通过以下来确定上施加PCA +，以柱中用柱的协方差矩阵，或由这种替代计算上减要求苛刻的的的程序，URE描绘在（ 图4c）。左路乘以的先前的由来的的的由的转置建议零售价矩阵的的得分矢量矩阵的情况，，以派生在按降序请订购矩阵的特征值（ 图4c，图6，图7a）时的体素的的柄的特征向量一个数组P于。
   ，P = ^SRPŢX S
   每一个的列向量代表的一个本金组成部分（PC）的的图像图案的SSM /主成分分析（PCA）的归因于一个的总方差的百分之占为（VAF）对应到的相对的它的的本征值大小的。

3。图案生物标志物的推导

1. 检查的的查询结果。如前述的的分析，以确定具有高VAF的值的相关联的的的的个人电脑型态，顷。我们的的例程以内，体素的模式向量是的Z-转化的，以便，他们的的价值观代表正面和负面的的标准的从他们的平均值的偏差。然后，他们正在重塑了成图像格式的，在显示之前（ 图7a）。
在某些情况，与这个疾病相关的所在地区公知的的偏差（RSD）的型态与可以目视被，识别。的的数学制剂认为作为在这种情况下，的相等于的解决方案的两个的个人电脑"和他们的相关联的的破门得分由零下的一方，，以符合到一个物理上正确的的的解释的，可能被乘以的派生的个人电脑的正面和负面的的形式的。
2. 对应到每一台PC模式的的的分数被显示为棒图和散点图（ 图7）的中 。一个可选的的ROC情节可以被生成的（ 图7d）。若要识别的一种疾病-特定的的图案中，注意到反映的相应的的两样本t-测试和AUC的值的由的p-部分的的值的患者和对照组的科目之间的恩怨的分化程度。限制本分析，以那些与一个高的的VAF和为一些固定的的截止的值（ 例如:，对<0.05和VAF> 5％）的高的的分化相关联的的个人电脑。通常情况下，只有一个或两个台个人电脑满足该标准，的，，用于PET的的数据。
3. 有是VARI适用于PC的选择是可以考虑的^44的 OU的其他的的方式。在那里本征值开始退化时，的流石滩的情节的序贯的特征值（ 图6）可能会给予一个尖锐的的的表示的的截止值由一个弯位在的曲线中，并一个尖锐的的减少在的曲线的斜率中的。另一种的做法是，以包括所有，占为顶级的的方差的的50％的的的个人电脑。一条被广泛接受的的程序是要计算的阿凯克的信息准则^{（AIC）45，}以确定的个人电脑的哪个组合定义该机型配备的最低的AIC值，，可以区分患者和对照组之间的的。
4. 所选择的的PC（s）可以是的矢量归一化的，并呈线性相结合，以产生一个单一的的与疾病-相关的的图案。我们的软件可选地使用的MATLAB函数glmfit，，，以确定通过基于上的后勤事项或其他的应用到学科的分数的回归模型的出来的系数。虽然的病人组和对照组的的分化通常提高了与在去的考虑的额外的PC设备rivation组，将得到的的的型态顷一个复合的的的表示形式，，表示形式，可能会无法的通讯到一个单一的的的物理网络或者"可能会纳入离群值的偏差（ 图7a和图7c的的 ）。
5. 进一步的验证是必需的可靠性和的前瞻性重要的意义。可以执行Bootstrap的重采样的，如下面所讨论^8，以确定的最可靠的的的的体素在的原来的的的的推导与的至少为标准的在的重复的花纹的推导中的偏差相关联的的数据集，之内。正向，验证被执行，以测试为的敏感性和特异性的独立的群体的歧视，通过派生使用的单一的-学科的的得分的评价方法（ 图4d）在接下来的分部该协议的描述的的患者和对照组的的准群体的的比分为。

4。单-学科的的的前瞻性分数评价使用一个预定的的生物标志物的

1. 一旦一个显着性的的SSM-GIS生物标志物的格局已被确定，一个的得分为及其在一项前瞻性主题的中的的表达可以进行评估从那个从个人的扫描，并使用一个简单的的计算的内部中断向量的商品被检者的建议零售价，矢量和所述地理信息系统图案矢量（ 图4d次，图7d）时中。
   SCORE =（乐谱）SRPS•病毒码
2. 先前操作的是自动化的由我们的的的TPR例行程序中。然而，为了独立地派生的主题的的SRPS矢量使用相关联的的内在的的GIS掩码上的日志变换后的数据，，并减去双方的的的主题的均值和在步骤2.4中，的相应的体素的值的的prederived的参考组的GMP作为。这确保了的得分情况可使用相比，到所预定的的的的参考值范围。一个新的组的遇分数核给可以被同样地进行评价作为一组的准的单一的主题破门得分。对于一个不同的的参考群组与参数设定的一一起使用，GMP，而模式是不变，可以重新校准。

Results

一个简单的的应用程序的多变量的SSM /主成分分析（PCA），，以推导出一个为PD人仕的神经影像学检查的生物标志物图案的就是如下图所示。的十的PET的FDG的图像临床确诊的帕金森氏症患者（6M/4F，59Y±部7y的sd）的持续时间可变的害病的的（9y的±5Y的sd）和十个年龄和性别相匹配的了正常对照组（6M/4F，58Y±部7y的sd）进行了分析用我们的的ssmpca例程。该二十个相应的空间上预先-归一的的图像被最初选定的的类别下，的疾病的科目或控件沿与的预设的遮罩为（ 图3）。该方案的输出的第一个的16个本金分量图像的文件和的他们的相关联的的分数沿与一个的序贯的特征值（ 图6）的流石滩的情节一个文件中房源用。是一张专属图片序贯的PC图片包含和相关联的的的的的酒吧和散点图被显示在屏幕上对之进行初期审查，如在图7a和7b中描绘的。需要注意的是的第一台PC成分占为的最大的VAF，为（24％）和它的相关联的的主题的分数显著的歧视患者从控件（p值= 0.0002）的基础上两样本t-测试对比和...

Discussion

SSM / PCA模型最初提出默勒等 ⁴已经演变成一个简单和强大的神经影像学数据分析技术^1-3。然而，已经有含糊之处，在应用这种方法，我们试图在这里澄清，并在以前的出版物^5-7,10。其中一些问题已得到解决在文本中，但在这里再次强调，因为他们的重要性。 SSM / PCA的详细介绍，主要是在休息状态FDG PET数据有效，但也已经成功地使用其他单量成像技术和扫描平台，包括H ₂ ^15O PET和鎝-99m-ECD SPECT，以及灌注翔升的体积VBM MRI。它的成功应用到时间序列rsfMRI必须大大增强了技术的未来发展^48。

SSM的主要力量在于，它能够识别源中的数据的变化，将它们分开成在空间上的正交分量。在大多数研究中，只有少数几个重要的组成部分是由于疾病的影响，从而显着降低了相关的数据集。要找出疾病的特定组件，我们最好能包括疾病的受试者性别和年龄相匹配的健康对照组，在联合用药组分析^2,5的数目相等。在这...

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
Image Acquisition
PET扫描仪	GE医疗系统	GE前进	、西门子和飞利浦PC工作站的任何PET、PET/CT和PET/MRI扫描仪
	联想	Any	http://www.lenovo.com/us/en/
[18F]-氟脱氧葡萄糖	范斯坦医学研究所	常规生产	也由 Cardinal Health 分销 http://www.cardinal.com/
Software
ScanVP	范斯坦医学研究所	版本 5.9.1，版本 6.2，即将发布	www.feinsteinneuroscience.org
SPM	伦敦大学学院神经病学研究所	spm99-spm8	http://www.fil.ion.ucl.ac.uk/spm
Windows	Microsoft	任何
Matlab	Mathworks	Matlab 版本 7.0、7.3	http://www.mathworks.com/
JMP	SAS	版本 5	http://www.jmp.com/