大脑皮层结构变异与一般认知能力：来自孟德尔随机化的证据

《Human Brain Mapping》：Cerebral Cortical Structural Variation and General Cognitive Ability: Evidence From Mendelian Randomization

　　

理解一般认知能力（general cognitive ability，GCA）个体差异背后的皮层架构仍是认知神经科学的核心问题。既往工作已确立全脑体积与GCA的关联，但区域效应及关系的方向性仍存争议。研究人员利用11,289名英国生物银行（UK Biobank，UKB）参与者的遗传信息驱动的皮层分区，考察了皮层表面积（surface area，SA）、皮层厚度（cortical thickness，CT）与通过言语–数字推理测量的GCA之间的关联。总SA与GCA呈稳健正相关。在区域水平上，背外侧前额叶和颞上回SA表现出最强的正向关联，且在调整全脑SA后依然显著。相比之下，CT的关联相对较弱。利用全基因组显著遗传工具变量进行孟德尔随机化（Mendelian randomization，MR）分析，研究人员观察到总SA与GCA之间存在符合双向关系的证据。在区域水平，背外侧前额叶和颞叶SA显示出MR推断的对GCA的方向性效应，而GCA显示出对总SA和perisylvian厚度MR推断的方向性效应。这些发现支持GCA之下存在多区域SA架构，其中前额叶和颞叶联合皮层的贡献尤为突出。结果细化了全脑–GCA模型，并凸显了遗传信息分区在识别区域皮层贡献中的价值。

研究背景与立项依据

一般认知能力（GCA）个体差异的遗传率可高达80%，但其神经基础仍不清楚。过去三十年的磁共振成像（MRI）研究确立了脑体积与GCA约0.24的正相关，近十年研究焦点由脑体积转向皮层表面积（SA）与皮层厚度（CT）——二者遗传率高、发育轨迹不同且解剖正交。证据表明脑体积–GCA联系主要由总SA而非平均CT驱动。大样本队列（如UKB、ABCD）报道了与GCA的区域关联，但基于汇总统计的GWAS研究也提示全局脑测量贡献认知。然而传统解剖图谱的分区方案、影像指标选择会影响检出力与可重复性；且全局测量虽最强，区域贡献的一致证据仍缺乏。考虑到皮层结构与GCA均受强遗传影响，采用遗传定义的皮层分区有望更敏感地捕捉区域关联。为此，研究人员在UKB欧洲裔样本（N=11,289）中，以言语–数字推理（field 20016）作为GCA代理，使用遗传信息驱动的SA与CT分区（各12区加全局测量），结合关联分析与双向孟德尔随机化（MR），细化GCA的皮层架构。论文发表于《Human Brain Mapping》。

主要关键技术方法

研究人员提取UKB中34,842名有神经影像的欧洲裔个体，经亲缘去除（GRM>0.1移除859人）与认知得分过滤得11,289人。基因型用UKB v3填补数据，按缺失率、HWE、MAF质控。T1像由三站点同型号Siemens Skyra 3T采集，FreeSurfer v5.3 recon-all产出顶点级SA（白质面）与CT（白质–软膜最短距），再以Chen等建立的遗传模糊聚类分区（12 SA区+12 CT区+全局）合并双半球。GCA取首次访视fluid intelligence原始分。形态测量先残差化去除站点、Euler数、年龄、性别、诊断，区域分析再加残差化去除对应全局测量（得相对测量），秩逆正态变换。多重比较用matSpD估有效独立表型数t_e=22，阈值p<2.3e?3。MR采用GSMR，暴露为Makowski等皮层GWAS、结局为Savage等智力GWAS（非重叠样本），工具变量为r²<0.05/250kb且p<5e?8的SNP，HEIDI-outlier滤水平多效，F>10。

研究结果

3.1 样本描述性人口学

男女在年龄、教育年数无临床差异；男总SA显著大于女（d=-1.27），女CT略厚（d=0.27）；男GCA略高（d=-0.12）。符合SA与身高相关、CT与SA弱负相关的已知模式。

3.2 与GCA相关的绝对与相对皮层形态

绝对模型中总SA与GCA正关联（β=0.36，p<1e?10，R²=2.4%）；所有区域SA均正关联，背外侧前额叶最强（β=0.35）。平均CT名义正关联弱（β=0.06）。调整总SA后，背外侧前额叶（β=0.07）与颞上（β=0.09）仍正关联，枕叶负关联；调整平均CT后，腹侧额叶、中颞、内侧前额叶CT呈负关联。说明SA贡献超越全局大小，CT呈异向梯度。

3.3 皮层形态与GCA的双向孟德尔随机化

正向MR：总SA→GCA方向性效应显著（β_xy=0.13，p=2.08e?10）；区域层背外侧前额叶（β_xy=0.09）、后外侧颞（β_xy=0.06）、眶额（β_xy=0.06）SA→GCA正向，顶上小叶负向；CT无全局/区域方向效应。反向MR：GCA→总SA（β_xy=0.31，p=1.02e?3）与GCA→perisylvian CT（β_xy=0.33，p=5.19e?4）显著，提示双向关系。多方法敏感性分析方向一致。

讨论浓缩

全局总SA而非平均CT复制为GCA最强相关（解释方差仅2.4%，属部分因子）。遗传分区下背外侧前额叶与颞上SA在控全局后留存，支持区域配置贡献超脑大小。SA–GCA呈前额叶-外侧颞联合皮层梯度，CT呈背–腹梯度且效应弱，契合扩展梯度与遗传地形。背外侧前额叶凸显与工作时间记忆/执行控制一致；遗传分区比Desikan更准确定位。MR支持总SA与GCA双向：总SA→GCA与GCA→总SA并存，区域层背外侧前额叶/后外侧颞SA→GCA，GCA→perisylvian厚度（经验依赖可塑性解释需谨慎）。结果细化P-FIT等分布式框架而非矛盾。SA与CT机制分离：SA反映早期皮层扩展，CT反映后期成熟/修剪。局限含双半球合并遮掩偏侧化、未做跨站统计调和、相对模型用表型匹配全局而非ICV、方差层扫描质量未建模、群体级小方差解释不支个体预测与生物决定论。

结论翻译

综上，结果支持GCA的多区域模型，特征为高级联合皮层内分布式SA扩展。总SA仍是最强全局相关，但区域配置——尤其背外侧前额叶与颞叶皮层——提供额外解释力。通过整合遗传信息分区与MR建模，本研究推进了关于全局与区域皮层组织如何共同贡献认知能力的理解。