Lindquist Martin A
Department of Biostatistics, Johns Hopkins University, Baltimore, MD.
bioRxiv. 2024 Apr 14:2023.10.06.561221. doi: 10.1101/2023.10.06.561221.
Resting-state fMRI (rs-fMRI) data is used to study the intrinsic functional connectivity (FC) in the human brain. In the past decade, interest has focused on studying the temporal dynamics of FC on short timescales, ranging from seconds to minutes. These studies of time-varying FC (TVFC) have enabled the classification of whole-brain dynamic FC profiles into distinct "brain states", defined as recurring whole-brain connectivity profiles reliably observed across subjects and sessions. The analysis of rs-fMRI data is complicated by the fact that the measured BOLD signal consists of changes induced by neuronal activation, as well as non-neuronal nuisance fluctuations that should be removed prior to further analysis. Thus, the data undergoes significant preprocesing prior to analysis. In previous work [24], we illustrated the potential pitfalls involved with using modular preprocessing pipelines, showing how later preprocessing steps can correlation with signal previously removed from the data. Here we show that the problem runs deeper, and that certain statistical analysis techniques can potentially interact with preprocessing and reintroduce correlations with previously removed signal. One such technique is the popular sliding window analysis, used to compute TVFC. In this paper, we discuss the problem both theoretically and empirically in application to test-retest rs-fMRI data. Importantly, we show that we are able to obtain essentially the same brain states and state transitions when analyzing motion induced signal as we do when analyzing the preprocessed but windowed data. Our results cast doubt on whether the estimated brain states obtained using sliding window analysis are neuronal in nature, or simply reflect non-neuronal nuisance signal variation (e.g., motion).
静息态功能磁共振成像(rs-fMRI)数据用于研究人类大脑的内在功能连接(FC)。在过去十年中,研究兴趣集中在短时间尺度(从几秒到几分钟)上FC的时间动态变化。这些时变FC(TVFC)研究已能够将全脑动态FC图谱分类为不同的“脑状态”,这些“脑状态”被定义为在不同受试者和实验环节中可靠观察到的反复出现的全脑连接图谱。rs-fMRI数据的分析较为复杂,因为测量到的血氧水平依赖(BOLD)信号既包含神经元激活引起的变化,也包含非神经元干扰波动,在进一步分析之前应将其去除。因此,数据在分析之前要经过大量预处理。在之前的工作[24]中,我们阐述了使用模块化预处理流程所涉及的潜在陷阱,展示了后续预处理步骤如何与之前从数据中去除的信号产生相关性。在此我们表明问题更为严重,某些统计分析技术可能会与预处理相互作用,并重新引入与之前去除信号的相关性。一种这样的技术就是用于计算TVFC的流行的滑动窗口分析。在本文中,我们从理论和实证两方面讨论了该问题在重测rs-fMRI数据中的应用。重要的是,我们表明在分析运动诱发信号时,我们能够获得与分析预处理后的加窗数据时基本相同的脑状态和状态转换。我们的结果让人怀疑使用滑动窗口分析获得的估计脑状态本质上是否为神经元的,还是仅仅反映了非神经元干扰信号变化(例如,运动)。