《Frontiers in Human Neuroscience》：Population-referenced multiregional brain-volume profiles associated with performance and concussion history in athletes: an exploratory study
编辑推荐：
引言：精英运动员可能表现出由长期训练、选拔进入精英运动、重复性头部冲击暴露或这些因素组合塑造的结构性脑特征。解释这一人群的脑结构具有挑战性，因为一般人群参考值可能无法区分与运动员相关的适应性变异与损伤相关或临床有意义的变化。在这项探索性横断面研究中，研究人员使
引言：精英运动员可能表现出由长期训练、选拔进入精英运动、重复性头部冲击暴露或这些因素组合塑造的结构性脑特征。解释这一人群的脑结构具有挑战性，因为一般人群参考值可能无法区分与运动员相关的适应性变异与损伤相关或临床有意义的变化。在这项探索性横断面研究中，研究人员使用年龄和性别参考的区域脑容量百分位数作为基于人群参考的框架，检验精英运动员的多区域结构特征是否与运动表现、感觉运动功能和报告的脑震荡史相关。方法：51名现役和退役精英运动员接受了3T结构磁共振成像（MRI）和机器人感觉运动测试（KINARM）。T1加权扫描使用BrainKey（brainkey.ai）进行分割，以获得18个脑区域的年龄和性别参考区域百分位数。主成分分析（PCA）直接应用于区域百分位数测量，贝叶斯模型评估与有序表现评级、KINARM成分和报告的脑震荡史的关联。结果：五个主成分解释了71.4%的方差。脑PC1与较高表现组隶属度呈正相关[β = 0.53, 95% CrI (0.04, 1.02), 方向后验概率 = 0.983, n = 51]。脑PC3与KINARM PC3呈正相关[β = 0.57, 95% CrI (0.01, 1.13), 方向后验概率 = 0.978, n = 44]。然而，这两个成分均表现出有限的自举稳定性，使得这一结果高度探索性。50名参与者提供了报告的脑震荡史，并显示出较低的脑PC2 [OR = 0.61, 95% CrI (0.36, 1.04), 方向后验概率 = 0.97]和较高的脑PC1 [OR = 1.57, 95% CrI (0.92, 2.71), 方向后验概率 = 0.95]的方向性证据。讨论：这些探索性发现表明，基于人群参考的多区域脑容量特征可能与运动表现、感觉运动功能和报告的脑震荡史相关。结果并未确立病理性或有益的结构性脑变化，也未构建或验证运动员特异性常模模型。相反，它们识别了候选的多区域特征，供在更大样本、纵向运动员队列中评估，并匹配非运动员对照组。
**论文解读：基于人群参考的多区域脑容量特征与运动员表现和脑震荡史的相关性**
**研究背景与问题**
精英运动员的大脑结构可能因长期训练、运动选拔、重复性头部冲击（如接触性运动中的脑震荡和次震级冲击）等因素而呈现独特特征。然而，使用一般人群的参考标准来评估运动员的脑健康存在局限性：一般人群参考值可能无法区分运动员因适应运动产生的正常变异与损伤相关或临床意义的变化。此外，重复性脑震荡和次震级冲击与脑结构或微结构变化的相关性研究结果不一，且精英表现、脑震荡史和累积头部冲击暴露的影响可能重叠甚至方向相反。现有常模建模方法（如基于年龄和性别的脑容量百分位数）虽能估算个体相对于预期分布的位置，但其有效性依赖于参考样本的组成、规模和可迁移性，而针对精英运动员是否系统偏离一般人群结构参考值的系统研究较少。因此，本研究旨在探索性地使用基于人群参考的区域脑容量百分位数，检验精英运动员的多区域结构特征是否与运动表现、感觉运动功能和报告的脑震荡史相关。该研究发表于《Frontiers in Human Neuroscience》。
**主要关键技术方法**
研究人员采用以下关键技术：1）**结构磁共振成像（MRI）**：使用3T Siemens Prisma Fit采集T1加权图像（MPRAGE序列，参数：TR/TE=2250/4ms，TI=925ms，翻转角9°，矩阵256×256，视野25.6cm，层厚1mm，分辨率1×1×1mm）。2）**自动脑分割与常模百分位数估计**：通过BrainKey（brainkey.ai）进行自动分割，将灰质、白质、脑干、小脑等区域体积归一化至总脑体积，并基于外部参考人群（3564名参与者，年龄6-93岁，约53%男性）生成年龄和性别参考的百分位数。3）**主成分分析（PCA）**：对18个脑区域的年龄-性别参考百分位数进行PCA降维，识别多区域结构变异的主成分模式（PC1-PC5），解释71.4%总方差。4）**贝叶斯建模**：以脑PC和性别为预测变量，对运动表现（有序累积-logit模型）、KINARM感觉运动成分（高斯回归模型）和脑震荡史（伯努利-logit模型）进行贝叶斯回归分析，使用弱先验（Normal(0,0.5)）并评估后验概率方向（Pr(direction)）。5）**自举稳定性评估**：对PCA解进行2000次重采样自举，评估成分负载的稳定性（基于Tucker同余系数）。样本来源：51名现役和退役精英运动员（2016年6月-2017年9月招募），其中43名男性，8名女性，来自多种运动项目（冰球最多，占45%）。
**研究结果**
**3.1 参与者特征**：表现组分为高（n=12）、中（n=15）、低（n=24）。性别比例组间差异不显著，年龄无显著差异（Kruskal-Wallis p=0.345）。冰球运动员在中等表现组中占比最高（80%）。
**3.2 脑区域成分**：PCA保留5个成分，PC1（22.3%）、PC2（21.6%）、PC3（10.6%）、PC4（9.6%）、PC5（7.3%），共解释71.4%方差。PC1主要负载脑干、双侧海马、苍白球、小脑和壳核的高百分位数，以及枕叶灰质和左侧侧脑室的低百分位数。
**3.3 脑成分的自举稳定性**：PC1和PC2的负载向量同余中位数分别为0.855和0.869（显示“尚可”相似性），PC3-PC5稳定性较弱（同余中位数0.664-0.749），提示PC1-PC2内部可重复性较好，PC3-PC5需谨慎解释。
**3.4 脑结构与运动表现**：贝叶斯累积-logit模型显示，脑PC1与较高表现组隶属度呈正相关（β=0.53，95% CrI [0.04, 1.02]，Pr(direction)=0.983，OR=1.70）。敏感性分析支持该关联在原始5级评分和二元高表现分类中持续存在，但在排除冰球运动员后略有减弱（β=0.31，95% CrI [-0.33, 0.96]，Pr(direction)=0.82），表明运动项目构成可能部分影响该关联。
**3.5 脑结构与KINARM感觉运动表现**：KINARM PCA保留5个成分（解释55%方差）。自举分析显示KINARM成分稳定性有限（中位数同余0.492-0.819）。贝叶斯回归中，脑PC3与KINARM PC3呈正相关（β=0.57，95% CrI [0.01, 1.13]，Pr(direction)=0.978，n=44），但两者均稳定性不足，该结果高度探索性。KINARM PC3主要负载“物体击打-回避”任务得分和手速，脑PC3特征为白质高百分位数及顶叶、枕叶灰质、小脑和丘脑低百分位数。
**3.6 脑结构与脑震荡史**：报告脑震荡史（23人有，27人无）与较低脑PC2（OR=0.61，95% CrI [0.36, 1.04]，Pr(direction)=0.97）和较高脑PC1（OR=1.57，95% CrI [0.92, 2.71]，Pr(direction)=0.95）呈方向性证据。脑PC2特征为双侧侧脑室、第三脑室和丘脑高百分位数，以及颞叶、顶叶、额叶灰质和杏仁核低百分位数，提示脑震荡史与相反模式（较低脑室/丘脑百分位数和较高灰质/杏仁核百分位数）相关。
**3.7 模型诊断**：所有贝叶斯模型收敛良好（R-hat最大值1.002，有效样本量最小值3326），后验预测检查未发现明显失配。
**讨论与结论**
**研究结论**：在这项探索性精英运动员队列中，基于人群参考的多区域脑容量特征与专家评定的运动表现、KINARM感觉运动表现和报告的脑震荡史显示了候选关联。这些横断面发现并未确立有益或病理性的结构脑变化，也无法区分先天性特质、训练、运动参与、头部冲击暴露、脑震荡史或抽样变异性的效应。该研究未构建或验证运动员特异性常模模型，但识别出的特征为在更大纵向队列中结合匹配非运动员对照、客观暴露测量和重复运动员内评估提供了外部评估框架。研究人员强调，结果应被视为假设生成，需进一步验证，尤其是脑PC1与表现和脑震荡史的双重关联可能反映运动专长、累积竞争暴露、运动类型和头部冲击机会的重叠。
