《PLOS Computational Biology》：A generalized scaling law reveals synchrony-driven reorganization of brain dynamics across human lifespan
编辑推荐：
摘要专业翻译
生态系统遵循普遍标度律，如泰勒定律（Taylor's law），该定律将种群方差与样本均值联系起来。其他具有整合与分离网络拓扑结构的复杂系统，尤其是人脑，是否遵循类似的标度律，以及这些标度律在同步神经活动下如何表现，仍属未知。为解答这一问题，研究
摘要专业翻译
生态系统遵循普遍标度律，如泰勒定律（Taylor's law），该定律将种群方差与样本均值联系起来。其他具有整合与分离网络拓扑结构的复杂系统，尤其是人脑，是否遵循类似的标度律，以及这些标度律在同步神经活动下如何表现，仍属未知。为解答这一问题，研究人员最初尝试应用泰勒定律的努力受阻，因为所有测量尺度上的神经信号同时包含正值与负值，而泰勒定律仅限于非负生态数据[1]。为规避这一限制，研究人员用测量的均方根（root-mean-square, RMS）替代样本均值来扩展空间泰勒定律，并将这一广义定律应用于去趋势的多变量神经时间序列。去趋势操作是去除数据中趋势和时间均值的必要步骤。在此框架内，研究人员探讨了脑网络的同步或协调活动如何影响标度律这一后续问题。研究人员使用功能协调的专用度量来量化同步性，并通过解析推导及使用多变量泊松（Poisson）、负二项（negative binomial）、均匀（uniform）和伽马（gamma）分布的数值模拟表明，标度指数与同步性呈负相关。将该方法应用于三个大型队列（N=840，年龄18–88岁）的人脑功能磁共振成像（fMRI）数据，揭示了标度指数随年龄的显著变化轨迹，健康老化以该指数因同步性诱导的显著降低为特征。同步性–标度关系在静息态活动中保持稳定，但在自然主义任务中随生命周期逐渐变化，且在边缘系统、皮层下和小脑亚网络中尤为显著，表明亚网络水平上存在功能保留。最后，注意缺陷/多动障碍（ADHD）个体表现出改变的同步性–标度耦合，凸显了这一广义标度指标的潜在临床与转化价值。
论文解读
一、研究背景与问题
复杂系统常表现出尺度不变的波动特征，人脑亦不例外，但长期以来缺乏一个统一的理论框架，将这种变异性与大脑的协调动力学联系起来，尤其是跨不同心理状态、健康老化及神经疾病的背景下。生态学中的泰勒定律（Taylor's law, TL）作为普遍标度律，将非负生态计数数据的方差与均值联系起来，已在生态学、物理学、气候学、经济学、生命科学和神经科学等领域得到广泛验证。然而，经典泰勒定律仅适用于非负数据分布，而神经信号在去趋势后同时跨越正负值，这一根本性障碍使得泰勒定律无法直接应用于脑信号分析。此外，先前关于同步性与标度律关系的研究忽视了经验数据中常见的时间均值和趋势所带来的偏差。因此，本研究旨在通过将样本均值替换为空间均方根（RMS），将空间泰勒定律推广至含正负值的多元信号，并系统探究大脑广泛协调活动（以平均空间同步性为指标）如何影响标度指数，以及这一关系在健康老化、不同脑状态和神经精神疾病中如何变化。
二、研究内容与主要结论
研究人员提出广义空间标度律，以空间均方根（RMS）替代样本均值，将经典泰勒定律扩展至去趋势后含正负值的多元神经时间序列。通过解析推导、数值模拟及大规模fMRI数据验证，研究人员证明标度指数与平均空间同步性呈负相关：同步性越强，标度指数越低。在静息态下，同步性–标度关联保持稳定，但在自然主义任务（感觉运动任务和观影任务）中随年龄发生显著重组。健康老化伴随同步性诱导的标度指数显著降低，且该效应在边缘系统、皮层下和小脑亚网络中最为显著。ADHD患者表现出同步性–标度耦合的显著破坏。该研究引入了一个简洁、可解释且广泛适用的指标，将长期协调活动与神经活动的变异性联系起来，可作为健康老化和神经精神疾病的潜在生物标志物。
三、关键技术与方法
研究人员首先对多元BOLD时间序列（T=261个时间点，n=116个脑区，基于AAL图谱）进行去趋势操作，去除时间均值和低频趋势，使数据围绕零均值波动且近似正态分布。其次，在每个时间点计算空间RMS（m）和空间方差（v），通过对log m与log v进行线性回归拟合得到空间标度律（SL）指数b。平均空间同步性通过功能连接（FC）矩阵的平均值量化，即BOLD信号两两皮尔逊相关。研究人员通过解析推导（基于协方差矩阵理论）和数值模拟（使用多变量泊松、负二项、均匀和伽马分布，采用公共元素法生成相关数据）验证标度指数与同步性的理论关系。为分离边缘分布贡献与同步性贡献，研究人员采用相位随机化（phase randomization）生成替代时间序列，将标度指数分解为边缘分量和同步性驱动分量。数据分析基于三个大型队列：CamCAN队列（N=645，年龄18–88岁，含静息态、感觉运动任务和观影任务三种脑状态）、NKI/Rockland队列（N=89，用于交叉验证）和UCLA-CNP队列（N=106，含40名ADHD患者和66名健康对照，用于临床验证）。
四、研究结果
数值模拟：同步性对标度律的影响。 在近零同步性条件下，空间标度指数接近2，表明纯随机过程驱动。随着同步性增加（ρ=0.001、0.01、0.1），标度指数单调下降，该趋势与数据分布类型无关，经泊松、负二项、均匀和伽马分布验证均一致。
经验数据分布：去趋势的效应。 原始BOLD信号含有非零时间均值和低频趋势，不符合正态分布。去趋势操作后，所有脑区的信号围绕零均值波动，呈近似正态分布（95%置信区间内拟合良好），为后续标度分析奠定基础。
任务驱动的协调动力学与空间标度律。 在原始相关性分析中，三种脑状态下标度指数与同步性均呈显著负相关。经相位随机化去除边缘分布贡献后，静息态下同步性驱动的标度指数与同步性无显著相关（p>0.005），而感觉运动任务（ρ=0.3866，R2=0.0318，p<0.005）和观影任务（R2=0.0311，p<0.005）中呈显著正相关。半解析方法和线性拟合方法在100次独立实现中结果一致。
老化大脑中同步性–标度关联的减弱。 将CamCAN队列分为青年组（YA，18–36岁，n=138）、中年组（MA，46–60岁，n=155）和老年组（OA，72–88岁，n=146）。在三种脑状态下，同步性–标度关联斜率均随年龄增长呈顺时针旋转趋势，即关联强度进行性减弱。1000次相位随机化验证了结果的一致性。感觉运动任务中斜率从青年组的0.7921降至老年组的0.4417，观影任务中从0.5957降至0.2728。
功能亚网络中的同步性–标度关联。 将116个脑区划分为七个功能亚网络：默认模式网络（DMN）、额顶网络（FPN）、感觉运动网络（SMN）、视觉网络（VN）、边缘系统、皮层下区域和小脑。边缘系统、皮层下和小脑区域在三种脑状态下均表现出最强的同步性–标度耦合；视觉网络在观影任务中表现出显著关联；DMN和SMN在观影任务中的关联强于静息态和感觉运动任务；FPN在静息态下表现出相对较强的关联。
交叉验证、临床相关性与与其他指标的比较。 NKI/Rockland队列（N=89）重复验证了老化相关的同步性–标度负相关关系，青年组（18–30岁）与老年组（59–83岁）比较显示显著顺时针旋转。UCLA-CNP队列中，ADHD患者（N=40）相比健康对照（N=66）表现出显著降低的同步性–标度关联斜率，经1000次独立实现确认。与元稳定性（metastability）、功能连接动力学熵和Kuramoto序参量等传统指标的比较显示，这些指标在静息态下均表现出最小的年龄相关变化，而在观影和感觉运动任务中显著改变，与同步性–标度结果一致。
五、讨论与结论
理论扩展。 经典泰勒定律仅适用于非负数据，研究人员通过以空间RMS替代样本均值，将广义空间标度律扩展至含正负值的多元信号。在近零同步性下标度指数接近2，代表纯随机过程的零基线；同步性增强则指数降低，标志着从空间不相关状态向集体协调状态的转变。
随机过程与协调动力学的分离。 去趋势操作去除时间均值和趋势，相位随机化进一步将标度指数分解为边缘分量和同步性驱动分量，从而将真正的协调动力学与信号分布或外部驱动的伪迹分离。
静息态与任务诱发活动的特征化。 静息态下同步性–标度无显著关联，反映内在功能网络架构对协调强度相对不敏感；任务诱发状态下显著正相关，表明外部驱动或目标导向的脑状态招募灵活协调模式，从根本上改变波动标度特性。
生命周期敏感性与统一简洁指标。 同步性–标度关联斜率在三种脑状态下均随年龄进行性减弱，对任务条件和年龄均敏感，优于传统单区域时间标度或自相关方法。老化主要损害大脑在认知或感觉需求下灵活招募和维持大规模同步活动的能力，而静息态内在动力学相对保留。
亚网络特异性。 边缘系统、皮层下和小脑区域作为整合枢纽表现出最强的同步性–标度耦合；视觉网络在观影任务中特异性增强；FPN在静息态下相对较强。这种网络选择性表明框架捕捉的是生物物理上有意义的空间异质性。
临床意义。 ADHD患者同步性–标度耦合的显著破坏提示该指标对神经精神疾病敏感，可作为跨诊断标志物，适用于自闭症谱系障碍、癫痫、精神分裂症和神经退行性疾病等以异常大规模同步为特征的疾病。
结论。 本研究通过将协调动力学、波动标度与脑健康统一于单一框架，建立了同步性诱导的标度律作为复杂脑动力学的简洁而普适的标志物，为健康老化和神经精神疾病的评估提供了新的定量工具。
