IF31.7｜休闲体力活动的遗传基础及其健康保护作用

休闲时间体力活动与工作/家务体力活动遗传上截然不同，前者对多种疾病具有保护作用。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

体力活动（PA）对健康有广泛益处，包括降低癌症、高血压和2型糖尿病等慢性病风险。个体PA水平受遗传因素影响，遗传度估计在48%至71%之间。既往全基因组关联研究（GWAS）已发现一些与PA相关的遗传位点，但多数研究集中于休闲时间PA，且样本量有限，发现的位点较少。此外，不同情境下的PA（如工作、家务与休闲）往往被等同对待，但流行病学证据提示其健康效益可能不同。利用大规模生物样本库（如百万退伍军人计划MVP和英国生物银行UKB）进行遗传分析，有助于揭示PA的遗传结构及其与健康结局的因果关系。

科学问题

卡在哪里： 此前GWAS发现的PA相关位点数量有限，且主要基于欧洲血统人群，缺乏跨血统分析。更重要的是，不同情境（休闲、工作、家务）的PA遗传结构差异及其对健康结局的差异化影响尚未被系统研究。

本文要回答： 利用MVP和UKB数据，鉴定休闲时间PA的遗传位点，并比较不同情境PA的遗传差异及其与多种健康结局的因果关系。

技术路线

作者主要用了：GWAS、跨血统Meta分析、LDSC、MAGMA、TWAS、MR、MVMR、PheWAS、genomic-SEM。

作者首先在MVP中按血统（EUR、AFR、AMR）对休闲时间体力活动（All-PA–leisure）进行GWAS，并与UKB的剧烈运动（SSOE）表型进行Meta分析以增加功效。随后通过LDSC计算遗传相关，利用MR探究因果关系，并用MVMR校正BMI的混杂。此外，作者对剧烈体力活动按情境（休闲、工作、家务）分别进行GWAS，通过遗传相关和genomic-SEM揭示其遗传异质性。最后，通过PheWAS和生存分析验证PA与健康结局的关联。

图 1　Fig. 1。图1：跨血统All-PA–leisure+SSOE Meta分析的曼哈顿图。横轴为染色体位置，纵轴为-log10(P)。红线为全基因组显著性阈值（P=5×10−8）。最显著位点位于3号染色体CADM2基因附近（rs62253088），其次为MST1R。绿色高亮表示与EUR All-PA–leisure Meta分析相比新增的位点，显示跨血统分析增加了发现。

核心亮点

亮点 1｜休闲PA的跨血统GWAS发现70个新位点

对MVP中189,812名欧洲血统、27,044名非洲血统和10,263名拉丁美洲血统个体的All-PA–leisure表型进行GWAS，并与UKB的SSOE表型进行跨血统Meta分析（总样本量548,617），鉴定出70个独立显著位点（56个基因座）。最显著位点为CADM2rs62253088（P=6.1×10−20），其次为MST1Rrs3733134（P=6.9×10−18）。SNP遗传度在EUR中为0.083±0.005，AFR中为0.034±0.017，AMR中为0.083±0.036。基因集分析富集于突触前活性区等神经元相关通路，组织特异性分析显示小脑显著富集。

图 2　Fig. 2。图2：a) EUR All-PA–leisure+SSOE与多种性状的遗传相关条形图，误差线为标准误。所有显示的性状均通过Benjamini-Hochberg FDR校正。可见与卒中、胃食管反流病等呈强负相关，与HDL-C、教育程度等呈正相关。b) 基于多基因评分的PheWAS结果，横轴为phecode类别，纵轴为-log10(P)。红线上方为Bonferroni显著阈值（P=3.99×10−5）。三角形方向表示关联方向，显示与糖尿病、心血管疾病等负相关。

亮点 2｜休闲PA与多种疾病呈负遗传相关

LDSC分析显示，All-PA–leisure+SSOE与多种疾病呈显著负遗传相关，最强为卒中（rg=−0.62±0.11，P=3.54×10−9）和胃食管反流病（rg=−0.482±0.023，P=9.9×10−97）。此外，与2型糖尿病、冠心病、肺癌、心力衰竭、腹主动脉瘤、骨关节炎等均呈负相关。局部遗传相关分析在61个位点发现91个显著局部相关，其中教育程度相关位点最多（31个），BMI次之（15个）。

图 3　Fig. 3。图3：a) 三种剧烈PA情境（休闲、工作、家务）与多种性状的遗传相关热图或条形图，星号、楔形和加号分别表示休闲与工作、休闲与家务、家务与工作之间的显著差异。可见休闲PA与收入、教育程度等呈正相关，而工作PA呈负相关。b) genomic-SEM双因子模型路径图，显示因子1（休闲/喜好PA）和因子2（职业/家务PA）的载荷及因子间相关（0.50）。

亮点 3｜MR证实休闲PA对多种疾病的保护作用

两样本MR分析显示，All-PA–leisure对BMI、2型糖尿病、休闲屏幕时间、每日吸烟量和胃食管反流病具有双向负效应；对腹主动脉瘤、心力衰竭、COVID-19住院、甘油三酯、表型年龄加速和骨关节炎具有单向保护效应。MVMR校正BMI后，对COVID-19住院的保护效应依然显著（β=−0.067±0.016，P=2.8×10−5），对HDL-C、甘油三酯、骨关节炎和胃食管反流病的保护效应也仍存在。

亮点 4｜不同情境剧烈PA遗传结构差异显著

对剧烈PA按情境（休闲、工作、家务）分别进行GWAS，发现Vig-PA–leisure与Vig-PA–work的遗传相关仅为0.46±0.07，而Vig-PA–home与Vig-PA–work高达0.97±0.05。与健康相关性状的遗传相关分析显示，Vig-PA–leisure与收入、教育程度、休闲屏幕时间等的相关模式与工作/家务PA显著不同。genomic-SEM分析支持双因子模型：因子1（休闲/喜好PA）载荷Vig-PA–leisure、All-PA–leisure+SSOE和喜好PA；因子2（职业/家务PA）载荷Vig-PA–home和Vig-PA–work，两因子相关为0.50。

亮点 5｜生存分析显示休闲剧烈PA降低全因死亡风险

Cox比例风险模型显示，Vig-PA–leisure显著降低全因死亡风险（β=−0.1173±0.0037，P≤2×10−16），Vig-PA–work效应微弱（β=−0.0061±0.0029，P=0.036），而Vig-PA–home无显著效应（β=−0.0052±0.0030，P=0.082）。即使校正收入后，休闲PA的保护效应依然显著。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： GWAS、跨血统Meta分析、LDSC、MAGMA、TWAS、MR、MVMR、PheWAS、genomic-SEM
• 公开数据： GSE76381
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事复杂性状遗传学、流行病学或运动科学的研究者，尤其是关注基因-环境交互和因果推断的人。

带走一句

做PA遗传研究时，必须区分情境；休闲PA的遗传工具变量可用于因果推断，而工作/家务PA可能混杂社会经济因素。