IF31.7｜HUNT大队列鉴定可重复的宿主遗传-肠道菌群关联

12,652人发现+21,976人复制，锁定6个基因座12个菌种关联，并揭示BMI对菌群的因果效应。

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

研究背景

肠道菌群参与消化、代谢和免疫调节，与多种疾病相关，但因果关系大多未明。菌群组成受饮食、药物、年龄和宿主遗传等多因素影响。双胞胎研究提示宿主遗传有贡献，但此前GWAS仅重复出LCT和ABO两个基因座与菌群组成相关，且样本量不足、缺乏复制或依赖16S rRNA测序，分类分辨率低。此外，尚无GWAS重复出宿主遗传对菌群功能潜力的影响。观察性研究报道菌群与BMI等宿主性状相关，但因果方向不清。孟德尔随机化（MR）可推断因果，但需要强大的菌群遗传工具和足够样本量的菌群GWAS汇总数据，而此前基于宏基因组的宿主GWAS样本量不超过7,738，难以满足MR需求。

科学问题

卡在哪里： 此前肠道菌群GWAS发现的宿主遗传位点极少且缺乏复制，尤其缺乏对菌群功能潜力的可重复遗传信号；同时，BMI与菌群组成的因果方向尚未明确，主要受限于菌群遗传工具不足和样本量小。

本文要回答： 在大型队列中系统鉴定影响肠道菌群组成和功能的宿主遗传变异，并通过复制验证；进一步利用遗传数据探究BMI与菌群的因果关联，以及新发现位点与疾病的潜在联系。

技术路线

作者主要用了：宏基因组测序、GWAS、KEGG功能模块、PheWAS、共定位分析、孟德尔随机化、RNAscope。

作者首先在HUNT队列12,652名欧洲血统个体中，对546个肠道菌种（相对丰度）和461个KEGG功能模块进行GWAS，筛选出达到研究显著性阈值的SNP-菌种/功能模块关联。随后在四个瑞典队列（n=16,017）和FINRISK队列（部分关联）中复制。对复制成功的位点，开展疾病PheWAS、共定位分析、代谢物关联和MR分析，以探索潜在机制和因果方向。整体设计从发现到复制再到功能解读，层层递进，确保遗传关联的可靠性。

图 1　Fig. 1。图1展示研究设计、遗传力和曼哈顿图。a显示HUNT发现队列和瑞典、芬兰复制队列的流程；b显示546个菌种的SNP遗传力估计，均值6.8%，部分菌种遗传力高达25%；c曼哈顿图展示全基因组关联信号，突出7个位点13个SNP-菌种关联，其中LCT和ABO位点信号最强。

核心亮点

亮点 1｜12个SNP-菌种关联在6个基因座成功复制

在HUNT发现队列中，13个SNP-菌种关联达到P<1.3×10^-10，其中12个在瑞典和芬兰复制队列中方向一致且P<3.8×10^-3，涉及6个基因座：LCT、ABO、HLA-DQB1、MUC12、SLC37A2和FUT2。其中LCT和ABO为已知位点，FUT2此前仅达到全基因组显著性但未复制，而HLA-DQB1、MUC12和SLC37A2为全新位点。敏感性分析（排除抗生素使用者、调整肠道蠕动、仅用无关个体、CLR转换、排除同居）显示效应估计基本不变，表明关联稳健。

图 2　Fig. 2。图2总结各基因座内独立遗传信号的数量。圆圈内数字表示研究显著性（SWS）信号数，灰色表示全基因组显著性（GWS）信号数。7个位点共13个SWS信号被选入复制，LCT和ABO位点信号最多，显示宿主遗传对菌群影响的复杂性。

亮点 2｜LCT乳糖酶位点关联Bifidobacterium adolescentis及功能模块

LCT位点rs182549-T（与乳糖酶持续型rs4988235-A强连锁，r2=0.996）与Bifidobacterium adolescentis相对丰度降低显著相关（P=2.1×10^-45），并降低四个KEGG功能模块的潜力，包括SenX3-RegX3双组分系统（P=5.5×10^-37）和谷氨酸转运系统（P=2.7×10^-37）。B. adolescentis丰度与这些模块高度相关（解释60-61%方差），提示乳糖不耐受宿主中该菌的代谢变化驱动功能信号。

图 3　Fig. 3。图3展示共定位和MUC12表达定位。a,b显示C. cateniformis与痔病共享因果变异（PP4=99.8%）；c,d显示C. cateniformis与直肠MUC12表达共享因果变异（PP4=99.96%）；e,f RNAscope显示MUC12 mRNA在结肠细胞和杯状细胞中表达，蛋白位于腔面，提示MUC12在肠道屏障中的作用。

亮点 3｜HLA-DQB1位点关联Agathobacter sp000434275与乳糜泻风险

rs28407950-T（位于HLA-DQB1附近）与Agathobacter sp000434275丰度增加相关，PheWAS显示该等位基因降低1型糖尿病（P=1.2×10^-183）、哮喘（P=4.3×10^-118）和乳糜泻（P=2.1×10^-84）风险。在HUNT中，Agathobacter sp000434275丰度与乳糜泻患病率呈负相关（OR=0.53，P=3.7×10^-12）。MR提示乳糜泻可能因果性降低该菌丰度（beta=-0.036，P=9.9×10^-3），但反向因果证据不足。

图 4　Fig. 4。图4比较BMI与α多样性的观察性和MR关联。a显示Shannon指数和丰富度与BMI的观察性负相关，MR因果效应方向一致且效应量相似；b显示546个菌种的观察性关联beta与MR因果beta高度相关（r=0.76），支持BMI对菌群组成的总体因果效应。

亮点 4｜MUC12位点共定位Coprobacillus cateniformis与痔病

rs4556017-T与Coprobacillus cateniformis丰度增加相关，且是直肠MUC12的eQTL。共定位分析显示MUC12位点共享因果变异影响C. cateniformis丰度和直肠MUC12表达（PP4=99.96%），并与痔病风险共享因果变异（PP4=99.8%）。RNAscope显示MUC12在结肠细胞和杯状细胞中高表达，蛋白定位于结肠细胞腔面，提示MUC12可能通过糖萼调节菌群并影响痔病风险。

亮点 5｜FUT2与ABO位点交互影响Mediterraneibacter torques

FUT2位点rs2287921与Clostridium sp900540255关联（P=3.4×10^-21），该SNP与功能变异rs601338（无分泌者）连锁（r2=0.65）。ABO位点rs550057与FUT2位点rs2287921对Mediterraneibacter torques丰度存在显著交互（P=5.2×10^-7），支持FUT2依赖的ABO抗原表达影响特定菌种。FUT2无分泌者状态还与心血管风险标志物相关，共定位分析提示Clostridium sp900540255与高胆固醇共享因果变异（PP4=80%）。

亮点 6｜MR支持BMI对菌群组成的因果效应

观察性分析显示Shannon指数和丰富度与BMI负相关。两样本MR表明遗传预测的BMI升高降低Shannon指数和丰富度，效应量与观察性关联相似。对546个菌种的MR分析发现104个菌种与BMI有名义显著因果关联（39个增加，65个减少），显著多于随机预期（P=3.3×10^-49）。观察性关联beta与MR因果beta高度正相关（Pearson r=0.76，P=7.0×10^-104），表明BMI对菌群组成有总体因果效应。

生信可带走

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

• 方法栈： 宏基因组测序、GWAS、KEGG功能模块、PheWAS、共定位分析、孟德尔随机化、RNAscope
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事肠道菌群GWAS、宿主-微生物互作遗传学、孟德尔随机化或疾病微生物组关联研究的研究者。

带走一句

大规模宏基因组GWAS+严格复制是发现可靠宿主-菌群遗传关联的关键；MR分析提示BMI是菌群组成的因果决定因素，为后续机制研究提供方向。