《Frontiers in Human Neuroscience》：Prefrontal cortical morphometric covariance and cognitive task performance in healthy young adults: an exploratory analysis of Human Connectome Project data
前额叶皮层（prefrontal cortex, PFC）参与多种认知任务，包括工作记忆、基于奖赏的决策、执行控制和元认知。本探索性分析考察了PFC皮层厚度与表面积协方差成分的个体表达是否与认知任务表现相关。样本包含来自人类连接组计划（Human Connectome Project, HCP）的743名健康成人（年龄22–36+岁）的MRI数据。研究对22个PFC区域的残差化厚度与表面积应用了Oblimin旋转主成分分析（principal component analysis, PCA），并控制了年龄、性别和颅内体积（intracranial volume, ICV）。成分得分作为预测变量分别纳入Flanker年龄校正得分（抑制控制）、赌博任务反应时差值（Gambling RT Difference，即奖赏与惩罚区块中相同“更大”猜测的反应时之差）以及2-Back中位反应时（工作记忆速度）的回归模型。平行分析支持两个厚度成分；研究按预设决定保留三成分解（全局/外侧PFC，RC1；眶额叶，RC2；前扣带回，RC3；解释53.7%方差），其中RC1与RC2在自举重采样中仅表现出中等稳定性。RC1与较低Flanker得分呈名义显著关联（β = ?0.103, p = 0.027），方向与预测相反；RC2与较小的Gambling RT差值相关（β = ?0.103, p = 0.025）；两者均未通过错误发现率（false discovery rate, FDR）校正（q = 0.054）。在评审阶段补充的未校正模型中，全皮层平均厚度单独稳健地预测了Flanker得分，但未预测Gambling RT差值；在PCA前回归掉平均厚度后，两个成分关联均不再显著。在表面积分析中，较大的眶额叶表面积（AC2）与更快的2-Back反应时相关（β = ?0.125, p = 0.003；在12项表面积检验中q = 0.032，在全部24项成分检验中q = 0.064），且不能用速度-准确性权衡解释。这一未预测的关联在多数敏感性分析中持续存在，但在PCA前回归掉总表面积后减弱且不再显著（β = ?0.065, p = 0.078）。总之，未发现稳健的前额叶特异性关联。全局皮层厚度更简约地解释了Flanker关联，Gambling RT差值关联在两种解释下均不稳健，而眶额叶表面积关联无法与全局表面积分离；所有关联均需重复验证。
一、研究背景与问题
前额叶皮层（prefrontal cortex, PFC）在认知神经科学中占据核心地位。PFC在灵长类谱系中的进化程度超过其他任何皮层区域，其功能支撑着适应性人类行为的核心调控过程，包括反应抑制、奖赏与代价评估、错误检测以及工作记忆。损伤研究、功能神经影像研究和发展研究均广泛支持这一观点。然而，PFC与健康成人认知功能表现相关的精确结构参数，即其形态学特征、皮层厚度、表面积及其相互关系，仍待阐明。PFC并非单一结构，而是由功能和细胞构筑异质性区域组成的联合体，各区域对认知具有不同贡献。外侧PFC大致可分为额上回、额中回和额下回，工作记忆和反应抑制是外侧PFC的一般主题，其中背外侧前额叶皮层深度参与需要注意力门控操作记忆信息的任务。眶额叶皮层（orbitofrontal cortex, OFC）及其腹内侧延伸参与编码未来结果的预期（奖赏价值）以及整合既往反馈以校准新选择。前扣带回皮层（anterior cingulate cortex, ACC）则倾向于在存在冲突或需要个体意识到当前策略失效时被激活。
PFC各亚区的功能已被广泛研究，但其结构组织（即形态特征协变模式）与功能的关系尚不清楚。结构协方差被认为具有多种生物学来源，包括发育时序的协调以及区域间使用相关活动的协调。皮层厚度与表面积具有 largely 独立的遗传影响：厚度更多与迁移后过程相关，如神经元分化、髓鞘化、树突树和突触；表面积更多与产前过程相关，主要由皮层祖细胞增殖产生的放射单位数量决定。因此，厚度和表面积结构协方差网络预期携带不同且可分离的脑功能信息。然而，目前厚度和表面积常被分别分析，整个PFC的形态组织及不同亚区如何协变很少被研究。
既往关于PFC结构与认知关系的研究结果 largely 不确定。虽有研究表明PFC结构改变与智力相关，但少有研究深入考察PFC。在成人中，主动修剪过程 largely 完成后，这种关系更不明确。尚无研究考察PFC的结构协方差（即群体水平的形态协变模式）及其与多个认知领域表现的关系，这正是本研究所探讨的问题。
二、研究内容与主要结论
本探索性分析考察了PFC内群体水平结构协方差成分的个体表达是否与认知表现变异相关。研究分析了人类连接组计划（Human Connectome Project, HCP）数据集参与者的结构MRI数据。通过对22个残差化皮层厚度变量和22个残差化表面积变量应用主成分分析（principal component analysis, PCA）来分析PFC。成分得分用于预测三项主要认知指标的表现：Flanker任务（抑制控制测量）、赌博任务（产生奖赏与惩罚区块间反应时测量）以及2-Back任务（工作记忆速度测量）。清洗后纳入743名健康成人。
研究旨在实现五个主要目标：（1）通过厚度和表面积测量考察PFC的结构协方差网络；（2）评估外侧PFC厚度网络是否预测抑制控制；（3）分析OFC结构对奖赏与惩罚区块间反应时的影响；（4）探索外侧PFC与ACC成分的联合表达是否预测反馈依从性选择行为；（5）确定PFC结构是否预测工作记忆功能速度。研究提出五项假设：H1，PFC皮层厚度将分解为与外侧执行、眶额评估和扣带监测系统相关的三个结构协方差网络；H2，Flanker任务上更好的抑制控制与外侧PFC厚度成分的更大表达相关；H3，更快的工作记忆反应时也与该外侧PFC成分相关；H4，增加的赢-留/输-换敏感性与外PFC和ACC成分表达的交互作用相关；H5，OFC厚度成分的更大表达预测更平衡的奖赏-惩罚决策行为。H4和H5因对赌博任务变量的解读有误而无法按原述检验。
研究主要结论为：前额叶厚度可分离为解剖布局类似既往文献描述的亚区划分的协方差成分；两个厚度成分与两个结局存在名义显著关联，但均未通过多重比较校正或全部敏感性分析；一个未预测的眶额叶表面积与更快2-Back反应时关联在其家族内通过校正并经受住多数敏感性分析，但在PCA前回归掉总表面积后减弱。总体而言，未发现稳健的前额叶特异性关联。
三、主要关键技术方法
研究使用HCP Young Adult开放获取表格数据（S1200样本，1,206名参与者），于2026年10月5日从ConnectomeDB导出。纳入标准为完成Full 3T Imaging Protocol且无编码QC问题，最终保留743名成人（年龄22–36+岁）。结构MRI数据在华盛顿大学定制3T Siemens Skyra扫描仪上采集，使用3D T1加权MPRAGE序列（TR = 2,400 ms，TE = 2.14 ms，TI = 1,000 ms，翻转角8°，体素大小0.7 mm各向同性）。皮层重建和体积分割使用FreeSurfer 5.3.0-HCP在HCP最小预处理流程中完成。数据按Desikan-Killiany图谱分割，获取每半球11个前额叶区域的厚度和表面积（共22个变量）。分析前，22个双侧前额叶皮层厚度变量和22个表面积变量经残差化去除年龄、性别和颅内体积（intracranial volume, ICV）作为干扰变量。对残差化变量分别进行独立PCA，使用psych包，采用Oblimin旋转以允许成分间相关。平行分析提示两个成分，但按先验假设（H1）保留三成分解。成分稳定性通过1,000次自举重采样的Tucker一致性系数评估。认知结局包括Flanker年龄校正得分、Gambling RT差值（奖赏区块减去惩罚区块中“更大”猜测的中位反应时）和2-Back中位反应时。采用多元线性回归，以成分得分为预测变量，年龄和性别为协变量。使用Benjamini-Hochberg错误发现率（false discovery rate, FDR）校正进行多重比较控制。敏感性分析包括加入平均皮层厚度或总表面积作为协变量、在PCA前回归掉全局测量、使用HC3稳健标准误、稳健M回归以及排除影响个案等。
四、研究结果
（一）描述性统计与缺失情况
完整案例样本（N = 743）的缺失率显示，赌博反应时变量缺失率最高：Gambling RT Reward Larger为2.23%，Gambling RT Punish Larger为1.44%，原因是参与者从未在某一区块类型内猜测“更大”则该单元格无定义中位反应时。Flanker得分无缺失，2-Back反应时缺失0.65%。FreeSurfer结构变量在筛选样本中全部764名参与者均完整。21名有缺失数据的参与者（占筛选样本2.75%）经列表删除移除。
（二）统计假设检验
八个回归模型（三项主要结局加Block Choice Score，各对应厚度和表面积）均经标准回归诊断。两个主要显著模型均未违反方差齐性（Breusch-Pagan检验：Flanker p = 0.568；Gambling RT p = 0.564）。Shapiro-Wilk检验显示所有模型残差非正态（所有p < 0.05），与Flanker得分的天花板倾向和反应时差值的重尾一致。N = 743时，普通最小二乘估计对中度非正态稳健。八个模型均未显示显著异方差性（所有Breusch-Pagan p ≥ 0.07）。多重共线性诊断显示所有模型VIF值均未超过1.7，远低于常规阈值5和10。Oblimin旋转产生的成分间相关（r = 0.22–0.60）表明成分间存在一定关联，但不足以引起有问题的多重共线性。
（三）主成分分析：结构协方差网络
皮层厚度成分：对22个残差化前额叶厚度变量的PCA揭示三个成分，共解释53.7%方差。RC1（全局/外侧PFC）SS载荷最高（7.11；32.3%方差），包括额上回、尾侧额中回、喙侧额中回和岛盖部厚度。RC2（眶额叶）SS载荷第二（2.87；13.1%方差），主要载荷于外侧OFC和内侧OFC以及眶部厚度。RC3（前扣带回）SS载荷第三（1.83；8.3%方差），载荷于双侧喙侧和尾侧ACC。Ward D2层次聚类 largely 支持PFC的三簇组织。自举重采样中，RC1中位Tucker一致性为0.99（95%区间0.70–1.00），RC2为0.97（0.49–0.99），RC3为0.98（0.87–0.99）；一致性低于0.85的比例在RC1为19.4%，RC2为28.0%，RC3仅0.7%。Oblimin旋转所得成分间相关：RC1与RC2中度相关（r = 0.60），RC1与RC3相关r = 0.29，RC2与RC3相关r = 0.22。
表面积成分：对22个残差化表面积变量的PCA产生三个成分，解释39.9%总方差，地形模式与厚度解部分不同。AC1（内侧/扣带表面积）载荷于内侧壁和扣带回，以及尾侧和喙侧额中回和额极表面积。AC2（眶额叶表面积）以OFC和眶部为中心，另有额上回表面积载荷。AC3（外侧额下表面积）以额下外侧PFC为中心，包括Broca区部分。厚度与表面积解仅部分对应（RC2与AC2的Tucker一致性最高为0.81，RC1与AC1为0.68），与其 largely 独立的遗传影响一致。
（四）回归结果：皮层厚度模型
抑制控制（Flanker任务）：Flanker年龄校正得分对RC1、RC2、RC3、年龄和性别的回归显著，F(5, 737) = 4.13, p = 0.001, R2 = 0.027。唯一名义显著的结构预测因子是RC1（全局/外侧PFC厚度）：B = ?1.02, 95% CI [?1.93, ?0.12], β = ?0.103, t(737) = ?2.22, p = 0.027, η2p = 0.007。该关联未通过FDR校正（四项假设检验中q = 0.054；全部12项厚度成分检验中q = 0.162）。负系数表明外侧PFC厚度成分的更大表达与较低抑制控制得分相关。RC1与全皮层平均厚度强相关（r = 0.85）。加入平均厚度作为协变量后，RC1系数下降约30%（B = ?0.72 vs. ?1.02），标准误从0.46升至0.72（VIF分别为4.00和4.90），关联不再显著（β = ?0.073, p = 0.32）；平均厚度在该模型中也不显著（p = 0.59），但单独可预测Flanker得分（B = ?11.08, p = 0.012）。在PCA前将平均厚度从区域厚度测量中回归掉后，RC1关联不显著（β = ?0.036, p = 0.33）。RC1关联在HC3标准误（p = 0.022）、稳健回归（p = 0.028）和排除30个影响个案后（p = 0.025）类似，但在未校正Flanker得分（β = ?0.085, p = 0.065）和残差化中省略ICV（p = 0.051）时较弱。RC1在Flanker模型中的效应量最大（η2p = 0.007），但低于常规小效应基准（η2p = 0.01）。RC2和RC3不显著（ps > 0.23）。性别是年龄校正抑制控制得分的进一步高度显著独立预测因子，β = +0.122, p < 0.001, η2p = 0.015，男性平均得分比女性高约2.4分。年龄不显著（p = 0.73）。
奖赏-惩罚区块情境（赌博任务）：Gambling RT差值（奖赏区块减惩罚区块“更大”猜测反应时）对厚度成分的回归显著，F(5, 737) = 2.58, p = 0.025, R2 = 0.017。唯一名义显著预测因子是RC2（OFC厚度）：B = ?7.35, 95% CI [?13.77, ?0.93], β = ?0.103, t(737) = ?2.25, p = 0.025, η2p = 0.007。该关联未通过FDR校正（q = 0.054）。OFC厚度成分表达更大的个体在奖赏和惩罚区块间反应时差异更小。“更小”猜测的对应关联方向相同但不显著（β = ?0.053, p = 0.25）。RC2关联在加入平均皮层厚度作为协变量时保持名义显著（β = ?0.098, p = 0.046；RC2与平均厚度r = 0.66），在HC3标准误（p = 0.027）和省略ICV残差化（p = 0.027）时亦保持，但在稳健回归（p = 0.13）或排除36个影响个案后（β = ?0.066, p = 0.17）不再显著。在PCA前将平均厚度回归掉后亦不显著（β = ?0.052, p = 0.16）。平均厚度单独也预测Gambling RT差值（B = ?69.35, p = 0.029），但在稳健回归下不显著（B = ?30.09, p = 0.31）。RC1和RC3不显著（ps > 0.69）。性别和年龄亦不显著。
Block Choice Score：因BCS无法检验H4，其模型在补充材料中报告。无厚度或表面积成分预测BCS（所有ps > 0.09），RC1 × RC3交互不显著（p = 0.615）。平均BCS为0.998（IQR = 0.92–1.08）。
工作记忆速度（2-Back任务）：前额叶厚度或厚度成分均未预测2-Back中位反应时（RC1: p = 0.893；所有ps > 0.43），模型不显著，F(5, 737) = 1.26, p = 0.281, R2 = 0.008。仅年龄预测2-Back中位反应时（β = +0.079, p = 0.034），表明年龄较大个体工作记忆反应时较慢。评审阶段补充的次要模型中，RC1与较低2-Back准确率相关（β = ?0.135, p = 0.003；六项准确率成分检验中q = 0.010）。该关联在控制平均厚度后保持（β = ?0.308, p < 0.001），该模型中平均厚度符号相反（β = +0.248, p = 0.002；VIF分别为4.00和4.90），模式与高度相关预测因子间的抑制一致；平均厚度单独不预测准确率（p = 0.37）。
（五）回归结果：表面积模型
表面积成分未预测Flanker得分、Gambling RT差值或BCS（所有ps > 0.06）。对于2-Back反应时，表面积模型显著，F(5, 737) = 2.86, p = 0.014, R2 = 0.019：较大的眶额叶表面积（AC2）与更快的2-Back反应时相关（B = ?17.80, 95% CI [?29.40, ?6.20], β = ?0.125, t(737) = ?3.01, p = 0.003, η2p = 0.012）。该关联在12项表面积成分检验中通过FDR校正（q = 0.032），但在全部24项成分检验中未通过（q = 0.064），且非预设。AC2与总皮层表面积相关r = 0.38，关联在加入总表面积作为协变量时保持（β = ?0.102, p = 0.019），在HC3标准误（p = 0.005）、稳健回归（p = 0.004）、排除46个影响个案后（β = ?0.126, p = 0.003）以及省略ICV残差化（β = ?0.151, p < 0.001）时均保持。AC2表达更大亦与较高2-Back准确率相关（β = +0.150, p < 0.001；六项准确率成分检验中q = 0.001），2-Back反应时关联在加入准确率作为协变量后减弱但保持（β = ?0.090, p = 0.028），用于排除速度-准确性权衡解释。然而，在PCA前将总表面积从区域表面积测量中回归掉后，重建的AC2与原始解高度一致（Tucker一致性0.91），但其与2-Back反应时的关联减弱约一半且不再显著（β = ?0.065, p = 0.078）。总表面积单独（与年龄和性别一同纳入，不含成分）与更快2-Back反应时相关（β = ?0.100, p = 0.032），在同时包含AC2的模型中系数为β = ?0.106（p = 0.070）；总表面积先经年龄、性别和ICV残差化后与2-Back反应时无关联（p = 0.37）。AC2关联因此无法与全局表面积干净分离，但全局表面积并未解释该关联。
（六）假设总结
H1仅部分支持：三个成分解剖上可解释，但平行分析支持两个成分，三成分解由决策而非平行分析指示保留，且RC1和RC2在相当比例的自举重采样中重组。H2未获支持：外侧PFC成分的更大表达与较低而非较高Flanker得分名义相关，该未预测发现未通过FDR校正。H3（外侧PFC预测2-Back反应时）未获支持。H4无法检验，因BCS为区块水平，不捕捉试次间对反馈的反应。H5无法按原述检验，因Gambling RT差值比较同一猜测跨区块情境；RC2与较小区块情境反应时差值名义相关，未通过FDR校正。表面积与2-Back反应时的关联不受任何假设覆盖。
五、讨论与结论
讨论部分提出三个要点。第一，前额叶厚度可分离为解剖布局类似既往文献描述的亚区划分的协方差成分，但未检验与功能网络的对应。第二，两个厚度成分与两个结局存在名义关联（外侧PFC厚度与Flanker得分，OFC厚度与Gambling RT差值），但均未通过多重比较校正或全部敏感性分析。第三，未预测的眶额叶表面积与更快2-Back反应时关联在其家族内通过校正并经受住多数敏感性分析，但在PCA前回归掉总表面积后减弱。两个名义关联各解释不到1%行为方差，未通过FDR校正，此类规模的脑-行为关联在该样本量下不稳定，且在选择统计显著性时倾向于膨胀，因此应视为待重复验证的候选效应，而非特定前额叶网络支持特定认知能力的证据。AC2关联略大（η2p = 0.012），但同样需谨慎。
关于结构组织，PCA识别出由大外侧PFC成分（RC1）、眶额成分（RC2）和前扣带成分（RC3）组成的三成分解，解释约53.7%方差。虽然成分数由决策而非平行分析设定，但皮层成分位置与已知PFC解剖和功能亚区拓扑一致。外侧与眶额成分间的中度相关（r = 0.60）在解剖上可预期，因外侧PFC和OFC通过钩束等回路 heavily 互联，且均与边缘区域共享连接。成分间相关未转化为统计冗余，回归诊断未显示高多重共线性（所有VIF ≤ 1.64）。
关于表面积与厚度，皮层表面积和厚度在本分析中仅部分重叠的协方差结构。表面积PCA仅解释39.9%方差，与厚度解仅部分对应（最大Tucker一致性0.81，位于RC2与AC2之间）。眶额叶表面积与2-Back反应时的关联在其家族内通过校正并经受住多数敏感性分析，且不能用速度-准确性权衡解释，但未预测、未在全部24项成分检验中通过校正、效应小（η2p = 0.012），在PCA前回归掉总表面积后减弱约一半且不再显著，在省略ICV残差化时增强（β = ?0.151），表明部分信号与整体皮层和头部大小共享方差。总表面积单独仅在ICV校正前与2-Back反应时相关，因此结果表明AC2关联无法与全局表面积干净分离，而非全局表面积解释了该关联。
关于外侧前额叶厚度与抑制控制，外侧PFC厚度成分（RC1）与Flanker任务表现呈名义负相关。但RC1与全皮层平均厚度强相关（r = 0.85）。加入平均厚度使RC1系数下降约30%但标准误增加超过一半，平均厚度单独也预测Flanker得分，因此数据无法区分全局厚度效应与外侧PFC特异性效应。在PCA前将平均厚度回归掉后RC1关联不显著。对于Flanker得分，平均厚度单独预测表现（包括稳健回归下），RC1关联未能在PCA前回归掉平均厚度后存活，因此更简约的解释是更大全局皮层厚度与本样本中稍差的Flanker表现相关。对于Gambling RT差值，平均厚度关联在稳健回归下消失（p = 0.31），RC2关联同样消失，因此全局或眶额解释均不受支持。两种情况下均无法证明前额叶特异性厚度效应。
关于眶额叶厚度与奖赏-惩罚情境效应，OFC厚度（RC2）与Gambling RT差值的负关联意味着RC2表达更大的个体对同一猜测在奖赏和惩罚区块间的反应速度差异更小。该关联在稳健回归（p = 0.13）或排除36个影响个案后（系数约减半，B = ?3.90）消失，对“更小”猜测不显著（p = 0.25），在PCA前回归掉平均厚度后亦不显著（p = 0.16）。一种与OFC维持结果预期并随反馈更新一致的可能性是，更大RC2表达与对近期结果历史较不敏感的反应速度相关；鉴于这些结果，该关联不应进一步解读。
关于无效发现，BCS无法检验H4：在区块未提示且结果预定的情况下，即使严格遵循赢-留/输-换规则的参与者预期得分也接近1，因此其无效结果无信息量。厚度与工作记忆反应时的无效结果亦不确定。处理速度是由分布广泛的脑组织网络（包括额顶网络和白质束）介导的网络功能。白质属性如纤维束完整性可能比局部皮层厚度更相关。
局限性方面：第一，横断面设计无法确定结构协方差对认知的因果效应，两者可能受共同遗传或环境因素影响，需纵向研究阐明时间关系。第二，HCP样本包含双胞胎和非双胞胎兄弟姐妹，因家族标识仅可在HCP受限数据发布中获取而未建模家族结构。将相关参与者视为独立可能低估标准误，因此报告的p和q值可能反保守，需使用混合效应模型或交换性块置换进行重复验证。第三，两个预设厚度关联均未通过FDR校正，各解释不到1%行为方差（各η2p = 0.007）。
结论：未发现稳健的前额叶特异性关联。全局皮层厚度更简约地解释了Flanker关联，Gambling RT差值关联在全局或眶额解释下均不稳健，眶额叶表面积关联无法与全局表面积分离；所有关联均需重复验证。
