《Human Brain Mapping》：A Wavelet-Based Framework for Mapping Long Memory in Resting-State fMRI: Age-Related Changes in the Hippocampus From the ADHD-200 Dataset
编辑推荐：
功能磁共振成像（functional magnetic resonance imaging, fMRI）时间序列已知表现出超越传统建模方法能力的长期时间依赖性。本研究提出一种新颖的计算管线，利用长记忆（long-memory, LM）框架刻画和解释这些依赖性，
功能磁共振成像（functional magnetic resonance imaging, fMRI）时间序列已知表现出超越传统建模方法能力的长期时间依赖性。本研究提出一种新颖的计算管线，利用长记忆（long-memory, LM）框架刻画和解释这些依赖性，该框架捕获静息态fMRI（resting-state fMRI, rs-fMRI）信号中自相关的慢幂律衰减。该管线通过基于小波的贝叶斯方法进行体素级LM参数估计，生成反映全脑时间依赖性的空间图。随后，这些图通过复合基投影到低维空间，并通过组级回归与个体水平协变量关联。研究人员将该框架应用于ADHD-200数据集。分析在调整ADHD症状严重程度和用药状态后，识别出海马体中年龄与长记忆参数之间存在显著正相关。在纳入采集站点后，估计效应量减弱且不再达到统计学显著性，尽管空间模式基本保留，最强估计效应仍定位在相同解剖区域。这些发现提示长程时间依赖性与记忆相关脑区的发育变化之间存在关联，同时强调采集站点变异性对效应量和统计显著性的影响。总体而言，所提出的方法能够详细绘制rs-fMRI中的内在时间动态，并为功能信号记忆与脑发育之间的关系提供新见解。
静息态功能磁共振成像（resting-state functional magnetic resonance imaging, rs-fMRI）通过测量血氧水平依赖（blood oxygen level dependent, BOLD）信号变化间接反映神经活动，已广泛用于刻画大脑内在功能网络。BOLD时间序列并非独立同分布的白噪声，而是表现出持续数分钟的长程时间依赖。传统自回归（autoregressive, AR）模型假定短程相关结构，难以充分描述这种跨越多时间尺度的信号组织。虽然已有研究通过相关噪声建模或熵指标来刻画时间复杂性，但尚未有在全脑体素水平系统估计长记忆（long-memory, LM）参数并将其与个体临床协变量关联的分析框架。理解LM特性具有重要意义：其对应的谱密度在双对数坐标中的线性衰减（幂律行为）反映了神经活动内在的时间组织，并且与赫斯特指数（Hurst exponent, H）线性相关。为此，研究人员提出了一种基于离散小波变换（discrete wavelet transform, DWT）和贝叶斯推断的体素级LM建模完整管线，并利用ADHD-200公开数据集验证其应用价值。
研究人员提出两阶段分析流程。第一阶段（个体水平分析）：将每个体素的BOLD时间序列建模为零均值高斯长记忆过程，应用DWT将原始密集协方差矩阵近似为近对角结构，利用方差递推公式将小波系数方差表达为LM参数与创新方差（innovation variance）的函数；采用贝叶斯马尔可夫链蒙特卡洛（Markov chain Monte Carlo, MCMC）算法（截断Metropolis–Hastings与Gibbs采样相结合）逐体素获得参数后验分布，生成个体特异性LM空间图。第二阶段（组级分析）：借助复合基（composite basis）降维策略——先按自动解剖标记（automated anatomical labeling, AAL）图谱将大脑划分为感兴趣区（region of interest, ROI），各ROI内用奇异值分解（singular value decomposition, SVD）提取局部基；再将投影后的局部特征进行第二级SVD以捕获全脑尺度的全局结构——将29,639个体素压缩至411个主成分，继而建立组级贝叶斯回归模型，将LM空间模式与年龄、ADHD严重程度指数、用药状态及其交互作用、海马灰质密度等协变量关联，并将估计系数反投影回体素空间。多重比较校正采用联合置信带方法控制实验错误率，并辅以至少50个连续体素的空间簇阈值增强结果可解释性。
将其应用于ADHD-200数据集（Athena预处理流程）的355名7–18岁参与者后，组级分析发现年龄与海马体及小脑前叶的LM参数呈显著正相关；贝叶斯与频率派回归结果在显著区域上存在重叠，且贝叶斯方法检测到更大的激活区域。模拟研究的五个场景表明，LM估计器在参数宽动态范围、创新方差变化、0.1强度差异、激活区域缩小及加性噪声条件下均保持良好的准确性与稳健性；较小的激活区域和高噪声水平会导致精度下降，但主要空间模式仍可恢复。
研究结果分别概括如下：4.1 个体水平分析：在个体水平估计每个体素的LM参数与创新方差，并展示不同年龄与用药状态亚组的平均LM空间图，显示出可辨识的空间差异；4.2 组级分析：通过复合基降维和贝叶斯回归，识别出年龄系数图在海马体及小脑前叶存在显著激活簇（簇大小>50体素），且基于空间正则化的贝叶斯方法相比频率派方法具有更大的检测功效；4.3 站点效应与敏感性分析：单因素方差分析显示年龄、海马灰质密度、ADHD严重程度指数及估计的LM参数在各采集站点间均存在显著差异；将采集站点作为固定效应纳入模型后，平均绝对年龄系数下降约60%，不再有体素通过多重比较校正，但空间模式基本保留、最强效应仍位于相同解剖区域；各站点内部单独分析均未发现显著年龄效应，这受到样本量和年龄范围缩小的制约，提示ADHD-200数据中年龄与采集站点的强混杂限制了区分生物学年龄效应与站点变异的能力。
讨论部分强调，海马是记忆巩固与空间导航的核心结构，其在发育期经历显著的体积和功能变化；研究人员将年龄增长伴随的LM参数升高（即时间持久性降低与Hurst指数下降）解读为朝向更复杂、更灵活神经信号动态的发育转变，该解释与近期关于谱斜率与Lempel–Ziv复杂度（Lempel-Ziv complexity, LZc）呈负相关的证据方向一致。研究的主要局限性包括两阶段推断未形式化传播第一阶段参数不确定性，但通过嵌套MCMC按总方差分解显示，缺失的方差在多数脑区相对较小，说明两阶段框架在计算可行性与统计严谨性之间取得了务实平衡。研究结论为：识别出海马体年龄与静息态fMRI信号长记忆属性之间的正相关，估计的线性回归系数表明，年龄增长与较高的长记忆参数值相关，提示该区域存在朝向较少持续的时域结构的发育趋势；未观察到长记忆参数与ADHD指数、用药状态或其交互作用之间的显著关联。该研究首次在全脑体素水平系统分析LM时间复杂性并关联个体发育协变量，为理解儿童青少年期神经活动的时间组织提供了新视角。
