基于皮层编码的白质BOLD多变量预测：岭回归β指纹与边际功能连接G1分歧轴

《NeuroImage》：Multivariate Prediction of White Matter BOLD Activity and Its Relationship to Marginal Connectivity

　　

皮层活动可适度但可靠地预测留出的白质（WM）血氧水平依赖（BOLD）活动。
岭回归（ridge regression）系数指纹可重现边际灰质（GM）-白质（WM）功能连接（FC）。
岭回归梯度与功能连接（FC）梯度在静息态采集期次间呈现可重复差异。
梯度差异可组织预测残差与白质纤维束解剖。
干扰项重投影降低预测幅度，但保留空间模式。

该研究关注血氧水平依赖（BOLD）功能磁共振成像（fMRI）中长期以灰质（GM）为中心、白质（WM）信号常被视为伪迹或干扰的问题。尽管已有证据提示WM低频BOLD波动具有可重复功能信息，既往WM功能组织研究多依赖WM-GM边际功能连接（FC）。边际FC受广泛皮层同步、空间自相关和大尺度网络共线性限制，难以判断某一WM体素与皮层的关联究竟反映特定皮层贡献，还是与其他同步皮层系统共享方差。因此，研究人员以分布式皮层活动作为条件编码预测WM BOLD，并与边际FC比较，旨在检验皮层活动能否预测留出WM BOLD、岭回归β指纹与边际FC的关系，以及二者梯度差异是否形成可重复的分歧轴并组织预测残差、皮层网络贡献和WM纤维束解剖。研究基于人类连接组计划青年成人队列（HCP-Y）中经典unrelated-100子集的81名参与者，利用严格WM掩膜和Schaefer-400皮层分区开展体素级岭回归编码分析，并将β指纹与边际FC指纹、梯度组织、空间与信号质量对照及跨采集期次验证相结合。研究结论是：分布式皮层BOLD活动可适度但可靠地预测严格WM掩膜内BOLD波动；岭回归β指纹与边际FC共享强骨架，但二者低维梯度差异形成稳定的β-minus-FC G1分歧轴；该轴与经FC均方根校正的z标准化R²残差（FC-RMS-adjusted zR² residual）、皮层预测贡献谱及WM纤维束解剖相关。其意义在于，将WM fMRI从单纯边际FC描述推进到“边际FC骨架＋协方差校正条件编码”的互补框架，但并未确立WM信号的特异性神经来源或机制性主动作用。

关键技术方法：样本来自HCP-Y经典unrelated-100子集，要求四段静息态run均具完整生理记录、运动回归文件和预处理信息，最终81名。采用HCP最小预处理，以TAPAS PhysIO生成生理回归量，并与HCP运动回归量、线性趋势和截距共同回归，行0.01–0.1 Hz带通滤波。以Schaefer-400皮层分区时间序列为预测子、严格WM95体素为目标，行体素级岭回归，采用嵌套留一run交叉验证，正则化参数λ仅在外层训练数据内选择。比较岭回归β指纹与边际FC指纹；以BrainSpace扩散嵌入构建G1；采用环形移位零模型、跨采集期次、留一被试（LOSO）模板、共同坐标重建、空间块置换及距离与信号质量协变量控制。

研究结果如下。
3.1 基于皮层活动的严格WM预测：通过Schaefer-400皮层预测子对WM95目标的嵌套留一run岭回归，研究人员发现平均被试水平交叉验证R²为0.0199±0.0177，平均70.3%±7.1%的WM体素具有正留出R²；400分区模型优于17网络模型，滞后模型仅小幅增益；环形移位零模型显示预测超过时间自相关解释，二次干扰项重投影使平均R²降低31.1%但组图空间模式高度保留。结论为皮层活动对WM BOLD具有适度、可靠且空间组织的预测能力。
3.2 岭回归β指纹共享FC骨架但界定独特稳定梯度：通过比较β与边际FC指纹，研究发现二者高度相似，组水平符号轮廓相关为0.917、被试水平平均为0.948±0.018；原始高维FC指纹跨采集期次更稳定，但β第一梯度可重复性更高（r=0.965，FC为0.929），β与FC G1模板相关为0.701，其差异定义β-minus-FC G1轴；跨采集期次组水平β-FC对应仍达r=0.735，被试内为0.164；LOSO与G1映射显示共享组结构外仍存在可重复个体特异成分，K-means轮廓系数较弱，支持连续梯度而非强离散分区。结论为β保留边际FC骨架，但低维梯度差异稳定且可解释。
3.3 β权重广泛分布于皮层预测子：通过检查β分布及top-k修剪再拟合，研究发现β系数虽小但非均匀；保留200个预测子时留出R²几乎不变，保留100个时仅小降，进一步修剪至50、20、10个则R²下降、保留系数膨胀且全β轮廓相关性降低。结论为预测信息广泛分布于较宽皮层预测子集合，而非由少数主导分区驱动。
3.4 局部预测残差揭示超出边际FC强度的组织：通过将局部z标准化R²对局部FC均方根回归得到残差，研究发现残差图存在15个家庭wise错误率（FWER）显著簇，并与β-minus-FC G1相关（r=0.491）；岭回归相对仅使用完整400维边际FC指纹并统一缩放的基线具有平均ΔR²=0.01201且81名均正；ΔR²图保留残差组织并与G1相关；共同坐标重建高度复现原图，跨采集期次组图相关约0.32。结论为条件编码包含不能由局部边际FC强度充分概括的空间组织，并与G1分歧轴一致。
3.5 空间与信号质量控制支持WM特异效应：通过控制到GM距离、WM边界距离、WM出现率、时间标准差和空间坐标，研究发现残差与G1相关由0.491降至0.435但仍显著，控制后两图显著体素及同号联合体素超过12 mm空间块置换零分布；距离分层分析四层相关均正，后部丘脑辐射/视辐射区域去分层均值后仍正。结论为在测试控制下该效应不能仅由GM邻近、信号质量或平滑空间自相关解释，但作者未据此声称纯神经来源。
3.6 β-minus-FC G1分歧轴组织网络、纤维束与验证结果：通过沿G1十分位分组、皮层网络消融、JHU纤维束映射及跨采集期次验证，研究发现低G1端边际FC较高但残差为负，高G1端边际FC较低但残差为正；低G1端贡献较多来自显著性/腹侧注意、躯体运动与边缘网络，高G1端视觉中央与默认网络贡献较大；低G1端富集胼胝体及内/外囊，高G1端富集后部丘脑辐射、视辐射、后放射冠和矢状层，纤维束水平残差与G1相关为0.658；主要图谱呈双侧对称，跨采集期次组图可重复性较高，但探索性行为预测无校正显著关联。结论为该轴连续组织预测残差、皮层网络贡献和WM解剖，并具组水平可重复性。

讨论部分总结：作者认为边际FC与协方差校正编码提供互补描述；边际FC刻画主导同步骨架，条件岭回归刻画分布式皮层活动对局部WM信号的联合预测。β并非独立于边际FC的信息源，而是在给定皮层预测子相关矩阵和λ下对边际FC指纹的协方差与正则化重加权；二者高度相似说明编码结构锚定共同大尺度GM-WM耦合，而稳定低维G1差异比原始高维β更可解释。跨采集期次、共同坐标和空间对照支持空间组织的可重复性，但并非独立队列复制；干扰项重投影降低预测幅度却保留空间模式，提示幅度对预处理敏感而空间组织较稳定。局限包括交叉验证R²绝对値较小、β非因果或有效连接、fMRI为间接测量、样本为HCP-Y青年健康成人且最终样本量有限、行为关联未校正显著、未进行性别分层分析，且未完成ICA-FIX去噪复制或无平滑重分析。

研究结论翻译：总之，分布式皮层BOLD活动可适度但可靠地预测严格WM掩膜内测得的BOLD波动。岭回归β指纹与边际FC共享共同骨架，而β-minus-FC G1对比定义了一个可重复轴，该轴与FC-RMS校正的zR²残差、皮层预测贡献谱和WM纤维束解剖相关。这些发现表明，多变量皮层编码刻画了边际GM-WM连接与其协方差校正预测表征之间的可重复差异；但并未确立WM信号的特异性神经来源或机制性主动作用。