青少年大脑结构与微结构图谱绘制：来自ABCD研究的洞见

《Developmental Cognitive Neuroscience》：Mapping the structure and microstructure of the adolescent brain: Insights from the ABCD study

　　

摘要专业翻译 青少年脑认知发育（Adolescent Brain Cognitive Development, ABCD）研究为刻画青少年大脑成熟及其相关因素提供了前所未有的机遇。本综述综合了ABCD研究中结构磁共振成像（structural magnetic resonance imaging, sMRI）和弥散磁共振成像（diffusion magnetic resonance imaging, dMRI）的发现，重点阐述了方法学进展及其对心理健康、肥胖和环境风险的启示。研究人员首先概述了ABCD研究的多中心协调采集方案、高完成率及纵向覆盖范围，强调了质量控制与扫描仪协调（scanner harmonization）的影响。结构MRI研究结果趋于一致：从童年晚期到青少年早期，皮层厚度（cortical thickness, CT）呈大体线性变薄；皮层表面积（surface area, SA）在青少年早期达到峰值后出现适度收缩；与认知、青春期、性别和精神病理学存在区域性模式化但较小的关联。遗传与双生子分析表明，CT和SA具有高度遗传性但部分可分离，且稳定的个体特异性方差并未完全由遗传因素捕获。随后，研究人员综述了弥散成像，重点聚焦于限制谱成像（Restriction Spectrum Imaging, RSI），该方法将弥散信号分解为受限（restricted）、受阻（hindered）和自由水（free-water）三个隔室，并在白质和灰质中提供对神经突和胞体敏感的指标。ABCD的RSI研究揭示了皮层-纹状体-丘脑环路中受限弥散随年龄广泛增加，并证明这些微结构信号具有高度遗传性。RSI还增强了对环境和代谢暴露的推断：空气污染、中心性肥胖和体重增加与细胞内弥散改变表现出稳健关联，尤其是在奖赏环路中。最后，研究人员描述了快速混合效应算法（Fast and Efficient Mixed Effects Algorithm, FEMA），该框架支持在ABCD规模下进行顶点级和体素级建模，并概述了将结构、弥散和基因组数据整合到发展知情的大脑-行为模型中的机遇。

论文解读文章

一、研究背景与问题

青春期是神经发育的关键时期，大脑结构、微结构和连接性在此阶段发生快速而动态的变化。这些变化支撑着认知能力、社会情绪加工和自我调节功能的精细化发展，这些能力对成年期功能至关重要，同时也与心理健康障碍、物质使用及其他不良结局的易感性升高相伴随。然而，在ABCD研究之前，受限于样本规模、跨中心协调不足和纵向覆盖有限，学界对青少年大脑宏观结构与微观结构发育轨迹的精确刻画仍存在诸多争议。传统弥散张量成像（Diffusion Tensor Imaging, DTI）指标如各向异性分数（fractional anisotropy, FA）和平均弥散率（mean diffusivity, MD）虽为白质成熟提供了重要见解，但其单张量高斯假设限制了生物学特异性，并混淆了交叉纤维和多组织隔室的效应。此外，结构MRI指标（如皮层厚度、表面积和体积）的发育轨迹受头部尺寸校正方式、扫描仪效应、运动伪影和社会经济地位（socioeconomic status, SES）等混杂因素的显著影响，导致不同研究间的结论存在分歧。因此，亟需在大规模、协调化、纵向设计的队列中，结合先进的多隔室弥散模型和严格的统计方法学，系统解析青少年大脑发育的宏观与微观结构特征及其与认知、心理健康、肥胖和环境暴露的关联。

二、研究内容与主要结论

ABCD研究于2016年启动，是美国迄今最大的青少年纵向神经影像研究，招募了来自21个研究地点的近12000名9至10岁儿童，结合协调化影像采集、深度表型分析和纵向随访。该研究通过统一的多中心采集方案（涵盖Siemens、GE和Philips扫描仪）、严格的预处理流程和持续的体模质量保证，实现了跨站点、跨波次的标准化数据采集。结构MRI采用1.0mm各向同性高分辨率T1和T2加权扫描，支持基于表面的形态测量和皮层下分割；弥散MRI采用多壳层采集方案，包含7个b=0体积和96个弥散加权方向，分布在b=500、1000、2000和3000 s/mm2四个壳层，共103个体积，为RSI等多隔室模型提供了优化的高b值角度覆盖。

基于ABCD数据的研究揭示了以下主要发现：皮层厚度从童年晚期到青少年早期呈广泛且近单调的线性变薄，是四个宏观结构指标中年龄效应最显著者；表面积呈倒U形轨迹，在青少年早期达到峰值后出现适度收缩，且变化方向存在区域异质性；皮层体积作为CT和SA的复合指标，从童年中期起持续下降。纵向稳定性分析显示，CT的组内相关系数（intraclass correlation, ICC）均值仅为0.76，低于SA（0.93）和CV（0.93），表明CT是最具年龄敏感性但也最不稳定的宏观结构信号。在认知关联方面，结构指标与认知表现的相关性通常较小（r≈0.01，上限极少超过0.16），且社会经济地位调整可使皮层区域化所解释的认知方差降低约四至八倍。整合灰质和白质多模态数据的正则化结构方程模型可解释19.0%的认知表现方差，高于单独使用灰质（15.4%）或白质（12.4%）的模型，但该估计可能包含环境混杂因素的贡献。在青春期与性别差异方面，青春期状态和速率可解释超出实际年龄的全局CT和SA方差，女性中基线青春期状态可预测随后的皮层变薄（β=-0.08）。遗传分析表明，CT和SA高度遗传但部分可分离，其遗传架构在不同皮层区域存在差异；ENIGMA全基因组关联研究（genome-wide association study, GWAS）证实，SA相关变异聚集在胎儿中期皮层发生过程中活跃的调控区域，而CT信号则与突触组织和重塑等出生后过程更相关。

在弥散成像方面，ABCD的RSI研究揭示了皮层-纹状体-丘脑环路中受限弥散随年龄广泛增加，这些微结构信号具有高度遗传性。RSI指标在环境与代谢暴露研究中展现出更强的推断能力：空气污染、中心性肥胖和体重增加与细胞内弥散改变表现出稳健关联，尤其是在奖赏环路中。在方法学方面，ABCD的规模使得极小的关联（|r|≈0.04）即可在99%功效下被检测，这要求研究者将解释权重放在效应量和协变量选择上，而非统计显著性。扫描仪效应评估显示，制造商和个体扫描仪通常是结构及弥散特征中扫描仪诱导方差的主要来源，ComBat经验贝叶斯校正可显著提升统计功效，尤其是在较小亚样本中。FEMA框架支持在ABCD规模下进行顶点级和体素级的快速混合效应建模，为大规模纵向分析提供了高效工具。

三、主要关键技术方法

研究人员使用了以下关键技术方法：多中心协调化MRI采集方案（涵盖Siemens、GE和Philips三种扫描仪平台）；基于FreeSurfer的表面形态测量管线（用于皮层厚度、表面积和体积的提取）；多壳层弥散加权采集方案（四个壳层，b值分别为500、1000、2000和3000 s/mm2）；限制谱成像（RSI）多隔室弥散模型（将弥散信号分解为受限、受阻和自由水隔室）；ComBat、CovBat和RELIEF等回顾性协调化方法（用于消除站点和扫描仪效应）；双生子模型和基于全基因组单核苷酸多态性（single nucleotide polymorphism, SNP）的遗传分析方法；以及快速混合效应算法（FEMA）框架（用于大规模顶点级和体素级纵向建模）。样本队列来源于ABCD研究，包括基线近12000名9至10岁儿童，纵向随访至2年、4年和6年。

四、研究结果详述

典型大脑发育轨迹。 基于6323名9至13岁儿童（基线及2年随访两次扫描）的纵向分析显示，皮层厚度广泛变薄和皮层体积广泛下降是该年龄段两个最大的年龄效应；表面积变化幅度较小且方向区域混合，额颞叶和内侧枕叶增加伴随顶颞叶减少。值得注意的是，个体层面的模态模式是变薄伴随扩张而非收缩（男性39.3%、女性37.6%的区域-参与者观察）。

认知、青春期与性别差异。 横断面分析显示CT、SA与认知存在可靠但较小的关联，SA和CV比CT更能预测认知表现，白质体积优于DTI模型指标。青春期状态和速率可解释超出实际年龄的全局CT和SA方差，且该调节作用在正确进行脑尺寸校正后呈现区域模式化特征。

遗传与双生子分析。 CT和SA高度遗传但部分可分离，稳定的个体特异性方差并未完全由遗传因素捕获。CT遗传度从童年期到青春期增加，CT和SA的遗传架构在很大程度上相互独立，最强的遗传影响定位于不同皮层区域。FEMA衍生的A/C/S/E方差分解可分离加性遗传（A）、共享环境（C）、稳定个体特异性（S）和残差（E）对方差的贡献。

RSI微结构发现。 ABCD的RSI研究显示，皮层-纹状体-丘脑环路中受限弥散随年龄广泛增加，这些微结构信号具有高度遗传性。RSI指标在环境与代谢暴露研究中展现出更强的推断能力：空气污染、中心性肥胖和体重增加与细胞内弥散改变表现出稳健关联，尤其是在奖赏环路中。

扫描仪效应与协调化。 扫描仪诱导方差存在于每个结构及弥散MRI特征中，制造商和个体扫描仪通常是主要来源。ComBat在较小亚样本中可显著提升统计功效，但在大样本设置中仅将扫描仪作为协变量即可产生可比结果。RELIEF在平衡样本中常达到最低（最佳）的站点识别曲线下面积（area-under-the-curve, AUC），但在样本不平衡时性能下降。

完成率与缺失数据。 基线总体完成率为T1 99.2%、T2 95.7%、dMRI 95.2%；至2年、4年和6年随访时，T1完成率分别降至68.4%、53.6%和34.5%。白人和高收入参与者完成率略高于黑人和低收入群体，COVID-19相关中断显著加剧了流失。

五、讨论与结论

讨论部分强调，ABCD研究的前瞻性标准化减少了但并未消除站点和厂商效应。对于大样本单变量特征分析，将扫描仪/软件作为协变量建模可能已足够；但对于较小或不平衡样本，经验贝叶斯协调化（如ComBat）可恢复功效并减少偏倚。依赖跨特征结构的多元设置（如连接性、纤维追踪）可能受益于协方差感知方法（CovBat、RELIEF）或应用于原始弥散数据的信号级协调化。研究建议，最佳实践是应用协调化或显式建模扫描仪/软件协变量，并可能报告按制造商分层的敏感性分析以验证程序有效性。关于协变量选择，研究强调应基于原理和假设驱动而非反射性选择，建议明确每个协变量的假设作用，并报告含与不含该协变量的模型，比较效应量及其置信区间而非p值。全局脑尺寸建模是结构表型最关键的决策，因为SA和CV随头尺寸缩放而CT基本尺寸不变；建议分别报告CT和SA，将CV视为复合指标而非单一测量，并记录对替代缩放规格的敏感性。社会经济地位调整可显著减少脑结构对认知的表观预测方差，未调整的关联应被解释为脑特异性效应的上限。扫描仪内运动和噪声控制同样关键，因为运动与年龄或临床状态非随机相关，可诱导虚假的发育模式。研究结论指出，结构脑-行为效应量在ABCD中较小且理应如此，典型脑-行为相关接近r≈0.01，上限极少超过r≈0.16；小效应不一定不重要，因其可能仍具有临床或公共卫生相关性，可能被测量误差衰减，或可在众多因素中累积。ABCD的结构和弥散数据集能够以规模化的方式整合宏观和微观结构绘图，具备充分功效来检验竞争性发育模型、量化个体内变化，并解析与认知、心理健康症状、肥胖和环境暴露的微小但可靠的关联。未来研究可整合结构、弥散、功能和基因组数据，以刻画个体神经发育特征，并在这些儿童进入成年期后重新审视本文综述的模式、预测因子和机制。