《Developmental Cognitive Neuroscience》：Mapping the developmental course of the multispectral dynamics serving executive functioning
尽管支持基本感觉加工的皮层区域在儿童早期经历显著的成熟过程，且该过程在青春期明显减缓，但支持高级执行功能的神经环路和神经动力学在整个青春期和成年早期持续发育。虽然许多研究描绘了伴随发生的皮层厚度和脑激活变化，但定量刻画这一发育窗口中可能支撑执行功能改善的皮层动力学变化的研究却少得多。在本研究中，研究人员采用一种基于高密度脑磁图（magnetoencephalography, MEG）的高级动态功能映射技术，在80名典型发育的11-16岁青少年完成一项经过验证的执行功能任务期间，量化其神经振荡反应。所有MEG数据均被转换至时频域，并采用波束形成（beamforming）方法对显著的振荡反应进行成像。随后，研究人员以实足年龄、性别及性别×年龄交互项为预测变量，对所得到的体素级全脑功能图像进行发育变化检验。主要结果表明，右侧背外侧前额叶皮层和前扣带回皮层的alpha振荡随年龄增长显著减弱（即相对于基线负向偏离减小）。同样地，左侧前额叶皮层、左侧顶上小叶皮层及右侧辅助运动区的gamma振荡强度也随年龄增长而下降。至关重要的是，在双侧顶叶皮层中识别出alpha和gamma频段的显著年龄×性别交互效应。这些发现揭示了支持高级认知功能的关键皮层区域中alpha和gamma振荡动力学的发育性变化，这些变化可能支撑了青春期和成年早期执行功能的精细化。
## 论文解读
### 一、研究背景与问题提出
童年期和青春期是感觉运动及高级认知过程获得与精细化的敏感时期。既有研究已较为充分地记录了支持这些认知功能的结构架构的发育性改变，包括以突触修剪程度（即皮层厚度）和髓鞘化水平作为皮层成熟度的指标（Giedd等，1999；Gilmore等，2018；Gogtay等，2004；Norbom等，2021；Walhovd等，2017）。这些结构变化通常呈序列性和层级性，遵循从单模态到跨模态的梯度，使得皮层层级中较早成熟的脑区为联合皮层提供支架，从而在整个童年期和青春期支持日益复杂的认知、情绪和社会功能。神经功能的出现和精细化遵循类似路径，其中执行功能是发育时程最为延长的认知领域。这种延长的神经成熟期为考察支撑人类行为诸多关键方面的认知功能的发育轨迹提供了独特窗口。
既往研究已表明，这些发育过程易受多种来源干扰的影响，包括创伤暴露、环境毒素、物质使用及其与神经内分泌和炎症系统的交互作用。此类损害可能导致支持执行功能的脑区发生发育性改变，并与长期不良后果相关联，如学业表现下降、情绪失调和精神病理风险增加。
前额叶皮层被认为是执行功能发育的核心区域，在童年期和青春期表现出区域特异的成熟模式。大量执行功能任务的表现从儿童早期到青少年晚期持续改善，提示潜在的神经机制在发育早期即已出现，但其效率和/或有效性低于成人系统。值得注意的是，扣带回皮层和顶叶皮层在此窗口期也经历显著的结构和功能性发育变化，并且是分布式网络（如突显网络和背侧注意网络）的重要组成部分，这些网络在完成高级功能时被持续募集。多项研究将上述区域激活的年龄相关变化与行为改善联系起来，且大量工作表明，执行功能的行为改善似乎早于神经变化趋于稳定。具体而言，行为表现在童年期和青春期显著改善，进入成年期后趋于平台期；而神经变化则往往持续至成年期。有趣的是，青少年在表现出成人水平行为表现的同时，对前额叶皮层的依赖程度更高，这可能反映了代偿性募集，与神经募集随任务难度增加的研究一致。此外，在相似行为表现水平下，这些区域神经募集的减少可能反映了支持执行功能的神经环路不断优化，使得成人在相同区域内完成任务所需的神经活动少于青少年。
虽然基于功能磁共振成像（functional magnetic resonance imaging, fMRI）的神经发育研究革新了领域内对青少年执行功能基础的理解，但对这些脑区内神经动力学的了解仍然有限。脑磁图和脑电图（magnetoencephalography/electroencephalography, M/EEG）等补充技术具有出色的时间分辨率，使得聚焦于支持这些高级认知功能的振荡动力学研究成为可能。事实上，多项使用M/EEG的发育研究已识别出支持工作记忆、注意和其他执行功能过程的振荡动力学中的年龄相关神经变化，其中部分变化表现出性别依赖效应。理解这些振荡动力学是连接神经发育研究的关键环节，因为神经振荡通常可与特定的神经生物学机制相关联，例如高频gamma振荡与局部GABA能神经环路及小清蛋白阳性篮状细胞直接相关。
然而，关于执行功能在儿童与成人之间的显著差异已有较多研究，但针对青少年早期这一特定阶段发生的更为细微的神经变化的研究却相对匮乏。因此，本研究旨在利用高密度MEG和先进的动态功能映射技术，系统考察11-16岁健康青少年在执行功能任务中神经振荡动力学的年龄相关变化及其性别差异。
### 二、研究目的
本研究的主要目的是：（1）量化11-16岁典型发育青少年在执行功能任务中的alpha（8-14 Hz）和gamma（>30 Hz）频段神经振荡反应；（2）考察这些振荡反应的年龄相关变化；（3）检验年龄与性别之间的交互效应。基于发育性fMRI文献和青少年高级认知功能的MEG研究，研究人员假设：广泛参与执行功能的关键脑区（包括前额叶和扣带回皮层）将在alpha和gamma频段表现出年龄相关的变化，具体表现为随着任务有效募集能力的提升，alpha振荡减弱（即相对于基线更正向）和gamma振荡减弱；同时，鉴于女性青春期启动普遍较早，年龄与神经振荡之间的关系将在高级脑区表现出性别差异。
### 三、研究方法概述
本研究从奥马哈都会区招募了80名健康参与者（11-16岁，平均12.80岁，标准差1.23岁，其中女性35名）。排除标准包括：任何神经或精神疾病的正式诊断、可能影响中枢神经系统功能的主要医学状况、显著头部外伤史、或可能构成MRI安全风险或对MEG数据质量产生负面影响的金属植入物。
参与者均在MEG室中完成了一项知觉特征匹配任务（perceptual feature matching task, PFMT）以探测执行功能（即认知灵活性）。该任务的 behavioral 表现与此前研究证实与标准化神经心理测量工具中的传统执行功能指标高度相关。每次试验中，参与者首先注视中央十字，经过可变的刺激间隔（1500-1900 ms）后，同时呈现三幅图像，持续2700 ms。其中一幅目标图像位于注视十字上方，周围环绕正方形或菱形，指示该次试验的规则集；另外两幅图像位于注视十字下方。参与者需根据规则集（基于形状或基于图案）选择与目标图像匹配的底部图像。每位参与者完成200次试验。反应时和准确率被记录用于行为分析。
MEG数据采集使用配备306个传感器的MEGIN TRIUX Neo系统（204个平面梯度计，102个磁强计），采样率为1 kHz，采集带宽为0.1-330 Hz。结构MRI使用Siemens Prisma 3-Tesla扫描仪配合32通道头线圈获取T1加权MP-RAGE序列。MEG数据预处理包括头动校正、信号空间分离降噪、心电和眼电伪迹的信号空间投影去除。无伪迹的epochs通过复解调转换为时频域，频率采样为1 Hz，时间采样为50 ms，带宽为3-100 Hz。传感器级别统计采用两阶段方法：首先对每个像素进行相对于基线的配对t检验，然后采用非参数聚类置换检验（10,000次置换）进行多重比较校正。随后使用动态成像相干源（dynamic imaging of coherent sources, DICS）波束形成器对统计上显著定义的时频窗口进行源成像。每个体素的噪声归一化源功率以4 mm各向同性分辨率计算，生成pseudo-t图。全脑体素级ANCOVA模型以年龄、性别和性别×年龄交互项为预测变量，显著性阈值为p < .005且簇大小≥8个连续体素（即>500 mm3）。
### 四、研究结果
#### 4.1 行为与传感器水平结果
在原始80名参与者中，有2名因伪迹校正和运动补偿后MEG数据噪声过大而被排除。其余78名参与者完成任务表现良好，平均准确率为81.26%（标准差11.52%），平均反应时为1716.17 ms（标准差246.51 ms）。反应时分析未显示年龄（F = 1.96, p = .165）、性别（F = .003, p = .955）或性别×年龄交互（F = .025, p = .876）的显著效应。同样，准确率与年龄（F = .370, p = .545）、性别（F = .053, p = .812）或性别×年龄交互（F = .092, p = .763）均无显著关联。
传感器水平时频谱分析检测到相对于基线显著的持续性alpha功率降低（12-18 Hz；150-1750 ms），该窗口被分为四个不重叠的400 ms时间箱。此外，还检测到早期gamma（72-90 Hz，100-250 ms）和晚期gamma（74-84 Hz，350-750 ms）功率增加的显著簇。纳入最终MEG分析的试验数量与实足年龄无显著相关（r = .03, p = .790），排除了信噪比差异对结果的潜在混淆。
#### 4.2 Alpha振荡的年龄主效应
在第四个（最后一个）alpha时间窗口中，右侧前扣带回皮层（anterior cingulate cortex, ACC）和右侧背外侧前额叶皮层（dorsolateral prefrontal cortex, dlPFC）的alpha振荡随年龄增长显著减弱（即相对于基线负向偏离减小）。其中，ACC区域效应量为F = 12.52, r = .39，dlPFC区域效应量为F = 10.97, r = .37，均达到校正后显著性水平（p < .005）。这一结果表明，随着年龄增长，青少年在执行功能任务中上述关键脑区的alpha去同步化程度降低，可能反映了神经募集效率的提升。值得注意的是，alpha的年龄相关变化仅出现在最晚的时间窗口，提示任务表现的早期感觉加工方面可能在童年期已经完成巩固，而接近决策点的高级执行功能成分仍在积极发育中。
#### 4.3 Gamma振荡的年龄主效应
在早期gamma窗口中，左侧前额叶皮层（F = 16.13, r = -.39）、左侧顶上小叶皮层（F = 16.30, r = -.41）和右侧辅助运动区（supplementary motor area, SMA）（F = 11.89, r = -.32）的gamma振荡强度随年龄增长显著下降（均p < .005，校正后）。gamma振荡的减弱可能反映了GABA能抑制性环路的逐步成熟。小清蛋白阳性篮状细胞在gamma振荡的产生中发挥关键作用，并在发育过程中经历显著的重塑和成熟。这种成熟可能改变了兴奋/抑制（E/I）平衡，抑制自发神经噪声背景，同时增强低阶区域向高阶皮层输入的相干信号，从而提高信噪比，使联合皮层能够更有效地编码和整合信息。
#### 4.4 年龄×性别交互效应
在第四个alpha窗口中，右侧 inferior parietal cortex（IPL）表现出显著的性别×年龄交互效应（F = 10.83, p < .005，校正后）。具体而言，女性参与者的alpha振荡随年龄增长而增强（r = -.37, p = .035），而男性参与者则呈现相反模式（r = .36, p = .016）。在早期gamma窗口中，左侧 superior parietal cortex（SPL）同样表现出显著的交互效应（F = 11.08, p < .005，校正后）：男性参与者的gamma反应与年龄呈负相关（r = -.36, p = .018），而女性参与者未观察到显著关系（r = .21, p = .241）。这些双侧顶叶皮层的性别分化轨迹，可能反映了在行为表现无显著性别差异的前提下，不同性别参与者采取了不同的神经加工策略来维持同等的执行功能水平。
#### 4.5 其他窗口的结果
在前三个alpha时间窗口及晚期gamma窗口中，未检测到显著的年龄主效应、性别主效应或交互效应。这一阴性结果进一步支持了上述发现的特异性，即年龄相关的振荡动力学变化主要集中于执行功能加工后期阶段和早期gamma反应成分。
### 五、讨论与结论
#### 5.1 行为表现缺乏年龄效应的解释
本研究中未观察到任务表现随年龄、性别或交互项的显著改善。既往类似任务的元分析显示儿童期执行功能准确性显著改善，但大多数研究对象为明显更年幼的样本（30-80个月），聚焦于执行功能的出现而非长期成熟。Ezekiel等人分析了与本研究年龄范围相似的样本，同样未观察到反应时的年龄相关改善。总体上，执行功能在童年期随年龄增长而改善，但此类变化在青春期逐渐变小。因此，包括PFMT在内的大多数执行功能任务可能不够敏感，难以检测青春期和成年早期常见的较小规模行为改善。本研究数据在数值上支持这一可能性——准确率和反应时均表现出随年龄增长的一致但微小的改善趋势。缺乏行为差异虽然可被视为局限性，但也简化了结果解释，因为可以更确信所观察到的年龄相关神经变化未受到表现差异的潜在偏倚影响。
#### 5.2 Alpha和Gamma振荡减弱的功能意义
Alpha振荡减弱：既往研究已将alpha振荡与冲突监控、注意的顶部-down调控和反应选择等功能相关联，涉及flanker、Simon、Posner、抽象推理和工作记忆等多种执行功能任务。fMRI文献表明，青少年表现出成人水平的行为表现，但对dlPFC的依赖程度高于成人，这可能反映了更努力性的注意和代偿性募集。因此，本研究中ACC和dlPFC中alpha振荡随年龄增长而减弱，可能表明年龄较大的青少年能够更有效地募集这些区域以完成PFMT。值得强调的是，alpha振荡的发育变化仅限于最晚的时间窗口，这提示较早的、更多基于感觉的任务加工方面可能在童年早期已完全巩固，而执行功能的高阶成分（即最接近决策点的部分）仍在积极发育中。
Gamma振荡减弱：gamma振荡随年龄增长而减弱可能反映了从童年期到青春期抑制性环路的逐渐成熟。GABA能环路和小清蛋白阳性篮状细胞在gamma振荡的产生中发挥关键作用，并在整个发育过程中经历显著的重塑和成熟。这被认为导致兴奋/抑制平衡的改变，可能在抑制自发神经噪声背景的同时增强从低阶区域输入到高阶皮层的一致信号。信噪比的提高可能使联合皮层能够更有效地编码和整合信息，可能通过经验驱动的突触修剪实现。与PFMT完成过程中alpha反应的下降类似，观察到的gamma功率降低可能反映了有助于成功执行功能的神经环路成熟和效率提高。
#### 5.3 年龄×性别交互的潜在机制
本研究中观察到的双侧顶叶皮层的性别×年龄交互效应与既往文献中关于灰质、灰白质比例及结构连接组的发育性别差异一致。已有研究在基本感觉加工以及支持抽象推理和工作记忆任务的高级脑区中识别出神经功能的发育性别差异。女性参与者在顶叶皮层中alpha和gamma振荡随年龄增长而增强，支持了既往关于执行功能任务中女性相对于男性参与者顶叶皮层BOLD反应增加的发现。已知青春期启动时间存在性别差异——女性青少年平均早于男性进入青春期。这些发现可能反映了某些大脑区域和过程中同一神经发育轨迹的左侧偏移（即更早启动），但整体模式可能在不同脑区以及不同认知和神经过程之间存在相当大的差异。文献中观察到的发育性别差异也存在显著异质性，这可能反映了性别内青春期转变初始时间的不一致性，以及部分可归因于青春期性激素变化的早期差异的显现。
#### 5.4 研究局限
本研究存在以下局限性：第一，横断面设计限制了因果推断，未来应采用包含个体内变异性测量的纵向分析。第二，该年龄段显著受到青春期和高度交互的神经内分泌系统影响，未来应整合性激素、皮质醇和多巴胺等元素以建立振荡动力学与其神经生物学基础之间的联系。第三，未观察到显著行为效应，未来研究应扩大纳入标准以包括更年幼参与者，或增加任务难度以提高对青少年早期执行功能细微变化的敏感性。第四，任务探测的是情感中性情境下的“冷”执行功能，未来研究应评估情绪突显性对执行功能的影响。第五，实足年龄作为成熟度量尚不完整，未来应纳入Tanner分期和性激素数据以更全面地考察青春期阶段及激素级联与脑发育的交互作用。
#### 5.5 研究结论
本研究通过对80名11-16岁典型发育青少年在执行知觉特征匹配任务过程中的高密度MEG数据进行时频分析和波束形成源成像，揭示了支持执行功能的高阶跨模态区域中频谱特异性的发育变化。alpha振荡在右侧前扣带回皮层和右侧背外侧前额叶皮层随年龄增长而减弱，gamma振荡在左侧前额叶皮层、左侧顶上小叶皮层和右侧辅助运动区随年龄增长而减弱，且双侧顶叶皮层表现出显著的年龄×性别交互效应。这些发现在行为表现无显著年龄或性别差异的背景下出现，表明青少年的神经动力学正在经历朝向更成熟、更高效表型的精细化重构。该研究为绘制执行功能神经动力学的规范发育轨迹提供了重要数据，有望为识别神经发育轨迹异常并在行为表现显现之前进行早期干预提供窗口。作者贡献声明显示，Jake J. Son、Caroline E. Howard、Elizabeth Santos、Abraham D. Killanin、Danielle L. Rice、Mikki Schantell、Thomas W. Ward、Grace C. Ende、Anna T. Coutant、Christine M. Embury、Giorgia Picci和Tony W. Wilson共同完成了本文的研究工作。本研究由美国国立卫生研究院（NIH）多项基金资助（R01-MH121101、P20-GM144641、S10-OD028751、F30-MH134713、F30-MH130150、F31-DA056296）。所有作者均声明不存在生物医学经济利益或潜在利益冲突。
