IF31.7｜瑞典万人宏基因组GWAS揭示肠道分子环境塑造菌群

16,017人宏基因组GWAS发现OR51E1–OR51E2等8个位点关联14个菌种，并在挪威HUNT中复制。

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

研究背景

人类肠道微生物组对健康至关重要，但宿主遗传变异如何塑造其组成仍不清楚。此前仅有LCT和ABO两个位点被稳健重复地与特定菌种关联，且多数研究基于16S rRNA测序，物种分辨率有限。本研究利用来自四个瑞典人群研究的16,017名欧洲血统成年人的深度宏基因组数据，采用标准化流程进行物种级分类，并在挪威HUNT研究的12,652人中复制，以系统鉴定影响肠道菌群组成和多样性的宿主遗传位点。

科学问题

卡在哪里： 此前肠道微生物组GWAS受限于样本量小、测序方法分辨率低或缺乏标准化处理流程，导致可重复的遗传关联极少。除LCT和ABO外，其他候选位点如MED13L和FUT2未能在严格阈值下重复。此外，宿主遗传变异如何通过肠道生理功能（如内分泌细胞脂肪酸感应、胆汁酸代谢、黏液层组成）影响菌群的具体机制尚不清楚。

本文要回答： 利用大规模、高分辨率、标准化处理的宏基因组数据，系统鉴定与肠道菌群多样性和物种丰度相关的宿主遗传位点，并在独立队列中验证，以揭示宿主遗传调控菌群的关键生理途径。

技术路线

作者主要用了：宏基因组测序、GWAS、scRNA-seq、代谢组、共定位分析、孟德尔随机化。

作者对四个瑞典队列的粪便宏基因组数据进行统一物种分类和功能注释，计算α多样性和921个物种的相对丰度。随后分队列进行GWAS，采用REGENIE回归模型并校正年龄、性别、遗传主成分等协变量，再通过逆方差加权固定效应荟萃分析合并结果。对达到全研究显著性（P<5.4×10⁻¹¹）的关联，在挪威HUNT队列中复制。进一步利用scRNA-seq验证候选基因在肠道细胞中的表达，结合血浆和粪便代谢组数据、共定位分析和孟德尔随机化，探索遗传变异影响菌群的潜在机制。

图 1　Fig. 1。图1展示了四个瑞典发现队列和挪威HUNT复制队列的年龄和BMI分布。年龄密度图显示各队列年龄范围覆盖18-96岁，中位年龄约60岁；BMI小提琴图显示中位BMI约26 kg/m²，分布相似。这为后续荟萃分析提供了人群可比性基础。

核心亮点

亮点 1｜OR51E1–OR51E2位点关联菌群丰富度

发现rs10836441-T等位基因与肠道菌群丰富度降低相关（每T等位基因减少5.7个物种，P=1.9×10⁻⁹），并在HUNT中复制（P=2.1×10⁻³）。该位点包含编码肠道内分泌细胞脂肪酸化学感受器的OR51E1和OR51E2基因。scRNA-seq显示OR51E1在肠内分泌细胞中表达，OR51E2在结肠内分泌细胞和免疫细胞中表达。此外，FFAR1–FFAR2–FFAR3位点也接近全研究显著关联Pullichristensenella excrementipullorum，支持脂肪酸感应通路在菌群调控中的作用。

图 2　Fig. 2。图2为曼哈顿图，展示全基因组范围内SNP与菌群物种关联的P值。虚线为全研究显著性阈值（P<5.4×10⁻¹¹），实线为常规全基因组显著性阈值（P<5×10⁻⁸）。三角形代表二元性状（存在/缺失），圆圈代表连续性状（相对丰度）。紫色标记为此前GWAS未报道的位点。图中可见LCT、ABO、FUT2等位点突出。

亮点 2｜8个位点关联14个菌种，11个关联被复制

在全研究显著性阈值下，共鉴定出15个SNP-物种关联，涉及8个位点和14个菌种，其中11个关联在HUNT中严格复制。新位点包括MUC12、CORO7–HMOX2、SLC5A11、FOXP1和FUT3–FUT6。这些基因均在人肠道中表达，涉及黏液层组成、胆汁酸代谢和岩藻糖基化聚糖等通路。

图 3　Fig. 3。图3a为系统发育树，展示921个物种在门、纲、目、科、属分类层级上的分布，有遗传关联的门被着色。图3b为每个位点与菌种关联的点图，每个点代表一个物种，展示8个位点与14个菌种的关联模式，部分位点关联多个菌种。

亮点 3｜岩藻糖基化聚糖位点ABO、FUT2、FUT3–FUT6

确认ABO与Faecalibacterium longum关联，并发现与Mediterraneibacter torques的新关联。FUT2位点与Blautia A obeum等三个物种达到全研究显著，且与炎症性肠病共享遗传信号。FUT3–FUT6位点与Clostridium sp900540255关联，并与LDL胆固醇共定位。ABO和FUT2的交互分析显示，分泌者状态下A抗原表达影响M. torques丰度。

亮点 4｜MUC12位点关联Coprobacillus cateniformis

发现MUC12内含子变异rs4556017与Coprobacillus cateniformis丰度关联（P=2.0×10⁻¹²），并在HUNT中复制。该位点与粪便频率共定位。MUC12在结肠肠上皮细胞和杯状细胞中高表达，提示黏液层糖蛋白在菌群定植中的作用。

亮点 5｜CORO7–HMOX2和FOXP1位点关联Turicibacter等菌

CORO7–HMOX2位点与Turicibacter sanguinis和Clostridium saudiense关联，且与腰臀比调整BMI共定位。FOXP1位点与Intestinibacter sp900540355关联，并与LDL胆固醇共定位。孟德尔随机化提示Intestinibacter sp900540355可能升高LDL胆固醇。这些菌与胆汁酸代谢相关，血浆代谢组显示该位点变异与次级胆汁酸水平相关。

亮点 6｜SLC5A11位点关联产丁酸菌

SLC5A11位点rs55808472-A等位基因与Agathobaculum butyriciproducens丰度增加相关，该等位基因降低SLC5A11表达并降低血浆肌醇水平。A. butyriciproducens是产丁酸菌，提示宿主肌醇转运可能影响该菌的定植。

生信可带走

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

• 方法栈： 宏基因组测序、GWAS、scRNA-seq、代谢组、共定位分析、孟德尔随机化
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事微生物组GWAS、宿主-菌群互作机制研究或肠道生理学的研究者。

带走一句

大规模标准化宏基因组GWAS能可靠鉴定宿主遗传对菌群的影响，且候选基因多涉及肠道生理功能，提示后续分析应整合单细胞表达和代谢组数据来解析机制。