《Nature Genetics》：Genome-wide association analyses highlight the role of the intestinal molecular environment in human gut microbiota variation
编辑推荐：
本研究针对宿主遗传因素如何塑造肠道菌群组成这一长期悬而未决的问题，通过对大规模、高分辨率宏基因组数据的全基因组关联分析（GWAS），首次系统地揭示了多个与菌群组成显著相关的宿主遗传位点，特别是首次在研究中确立了肠内分泌细胞表达的短链脂肪酸（SCFA）化学传感器（OR51E1–OR51E2位点）与肠道菌群丰富度（richness）的关联。该研究不仅成功复现了先前已知的LCT、ABO等位点，还新发现了MUC12、CORO7–HMOX2、SLC5A11、FOXP1和FUT3–FUT6等多个全新位点，将宿主遗传与肠道生理功能（如肠内分泌脂肪酸感知、胆汁酸代谢和黏膜层构成）紧密联系起来，为理解宿主-微生物组互作机制开辟了新方向，并指出了潜在的精准医疗干预靶点。
肠道中居住着数万亿的微生物，它们与我们的健康息息相关，影响着从新陈代谢到免疫功能的方方面面。尽管环境因素，尤其是饮食，已被公认是塑造个体肠道菌群（gut microbiome）的关键力量，但一个更为根本的问题始终萦绕在科学家心头：我们的基因，即遗传自父母的“蓝图”，在其中扮演了怎样的角色？过往的双胞胎研究表明宿主遗传存在一定贡献，然而，通过全基因组关联研究（GWAS）明确识别出具有统计学显著性的具体基因位点却困难重重，仅有的少数重复发现也多局限于如LCT（乳糖酶）和ABO（血型）等个别基因。这仿佛在暗示，宿主的遗传密码对肠道这个“微观生态系统”的塑造力可能相当微弱或极为复杂。
那么，真相是否果真如此？为了更清晰地描绘宿主遗传对肠道菌群的贡献图谱，一个由多国科学家组成的团队进行了一项迄今为止规模最大、分辨率最高的探索。他们整合了来自瑞典四个基于人群队列研究的16,017名成年人的宏基因组数据，并在一项独立的挪威队列（HUNT研究，12,652人）中进行了严格的复现验证。这项雄心勃勃的研究成果于顶级遗传学期刊《自然·遗传学》（Nature Genetics）上发表，为我们揭开了宿主遗传调控肠道菌群的新篇章。
该研究主要采用了以下几种关键技术方法：首先，研究人员对来自瑞典SCAPIS、SIMPLER（包含V?stmanland和Uppsala子队列）和MOS四个队列的样本进行了高质量的全基因组基因分型和填充，并利用标准化的生物信息学流程（CHAMP分析器，基于HMR05参考目录）处理了所有队列的粪便宏基因组测序数据，以获得可比的物种和功能信息。其次，他们对微生物α多样性（丰富度、香农指数、逆辛普森指数）、921个物种、652个更高层级的分类群以及117个功能模块进行了全基因组关联分析（使用REGENIE软件），并进行了固定效应模型荟萃分析和独立队列复现。此外，为了解析关联位点的潜在生物学机制，他们还结合了单细胞RNA测序（scRNA-seq）数据来分析关键候选基因在肠道不同细胞类型中的表达，并利用血浆和粪便代谢组学数据、孟德尔随机化（Mendelian randomization, MR）分析、共定位（colocalization）分析等方法来探索遗传变异、菌群组成与宿主表型（如血脂、胆汁酸水平、腰臀比）之间的潜在因果关联。
研究结果
GWAS of deep shotgun metagenomic data from four Swedish studies profiled with a standardized pipeline (利用标准化流程分析的四个瑞典队列深度鸟枪法宏基因组数据的全基因组关联研究)
研究人员首先确认了参与者的基线特征和微生物组构成在不同队列中具有可比性。他们发现肠道微生物α多样性的遗传力估计值（heritability）为9%（香农指数）至20%（丰富度），低于双胞胎研究中的报告值。通过大规模GWAS，他们识别出数百个与微生物物种在基因组水平显著相关的位点。
A locus including genes encoding EEC receptors is implicated in gut microbial richness (一个包含编码肠内分泌细胞受体基因的位点与肠道微生物丰富度相关)
这是本研究最引人注目的发现之一。研究人员首次在研究中发现了位于染色体11上的 OR51E1–OR51E2基因座（lead variant rs10836441-T）与肠道菌群丰富度显著相关（每个T等位基因减少约5.7个物种）。该位点在挪威HUNT队列中成功复现。这个基因座编码的蛋白质是小鼠嗅觉受体Olfr558和Olfr78的人类同源物，近期被鉴定为肠内分泌细胞（enteroendocrine cells, EECs）感知肠道菌群衍生的短、中、支链脂肪酸的传感器。EECs是胃肠道上皮中分泌激素的细胞，在进食后的生理反应（如肠道蠕动和饱腹感）中发挥重要作用。单细胞RNA测序分析证实，OR51E1在整个肠道道的EECs中均有表达，而 OR51E2主要在结肠的EECs中表达。此外，在编码其他脂肪酸受体（如FFAR2和FFAR3，同样参与感知微生物源性SCFA）的 FFAR1–FFAR2–FFAR3基因座也观察到了与特定物种的近乎研究水平显著的关联。这些发现共同指向了肠内分泌细胞通过化学感应（chemosensing）机制在调控微生物组构成中的潜在关键作用。
Meta-analysis identified eight genetic loci associated with 14 microbial species at study-wide significance (荟萃分析确定了八个与14个微生物物种在研究水平显著相关的遗传位点)
通过设定严格的、经多重检验校正的研究水平显著阈值（P< 5.4 × 10?11），研究团队共识别出15个单核苷酸多态性-物种关联，涉及8个基因座和14个常见细菌物种。其中11个关联在HUNT队列中成功复现。这些基因座包括先前被多次报道的 LCT（乳糖酶）和 ABO（血型）位点，以及在本研究中首次达到该严格标准的 FUT2（分泌状态）位点，并新发现了五个位点：MUC12、CORO7–HMOX2、SLC5A11、FOXP1和 FUT3–FUT6。例如，LCT位点（与乳糖耐受相关）的变异与多种双歧杆菌（Bifidobacteriumspp.）和 Phocea massiliensis等物种的丰度变化相关。
Several associations support an important role of fucosylated glycans in microbiome regulation (多项关联支持岩藻糖基化聚糖在微生物组调控中的重要作用)
研究证实了与组织血型抗原（histo-blood group antigens）相关的基因在塑造菌群中的核心作用。除了已知的 ABO基因，研究发现 FUT2（决定分泌者状态，secretor status）的变异与 Blautia A obeum、Clostridiumsp900540255等物种在研究水平显著相关。这些抗原是存在于细胞表面和分泌物（包括胃肠道黏液层）中的岩藻糖基化聚糖，可作为许多肠道细菌的碳源和结合位点。研究还发现，ABO A抗原（A或AB血型）对 Mediterraneibacter torques丰度的遗传预测效应依赖于个体的分泌者状态，凸显了宿主糖生物学与特定细菌（如能降解A抗原的细菌）之间复杂的相互作用。此外，FUT3–FUT6基因座（与Lewis血型相关）也显示出与特定菌种的关联，进一步强调了岩藻糖基化聚糖在调控肠道微生物景观中的重要性。
Genes involved in the mucosal layer implicated in gut microbiome composition (涉及黏膜层的基因与肠道微生物组构成相关)
研究发现 MUC12基因（位于一个富含黏蛋白基因的区域，邻近 MUC3A和 MUC17）的遗传变异与 Coprobacillus cateniformis的丰度显著相关。黏蛋白（Mucins）是肠上皮糖萼的主要成分，是高度O-糖基化的糖蛋白。宿主聚糖在决定哪些细菌能够定植宿主方面起着至关重要的作用，同时也是肠道微生物的重要营养来源。该位点的变异此前已被报道与排便频率相关，共定位分析也支持 C. cateniformis与排便频率之间存在共享的遗传信号，这为遗传因素如何通过影响肠道黏膜环境进而调节菌群和肠道功能提供了新的线索。
Shared genetic background of Turicibactersp., Clostridium saudiense, Intestinibactersp900540355, adiposity traits and bile acids (Turicibactersp., Clostridium saudiense, Intestinibactersp900540355、肥胖特征和胆汁酸之间的共享遗传背景)
研究新发现了两个位点（CORO7–HMOX2和 FOXP1）与一组共享的细菌（包括 Turicibacter sanguinis、Clostridium saudiense和 Intestinibactersp900540355）显著相关。值得注意的是，CORO7–HMOX2位点先前被报道与调整体质指数后的腰臀比（waist-to-hip ratio adjusted for BMI, WHRadjBMI）相关。共定位分析显示，T. sanguinis和 C. saudiense的丰度与WHRadjBMI以及多个次级胆汁酸（secondary bile acids）的血浆水平存在共享的遗传信号。某些 Turicibacter菌株能够编码并产生胆汁盐水解酶，参与次级胆汁酸的生成。孟德尔随机化分析提示 Intestinibactersp9005540355 的丰度可能对低密度脂蛋白胆固醇（LDL cholesterol）水平有正向影响。这些发现将特定的细菌类群、宿主脂肪分布、胆汁酸代谢和血脂水平通过共同的遗传背景联系了起来，揭示了潜在的、由宿主遗传介导的微生物-宿主代谢轴。
Variants in the SLC5A11locus associated with a butyrate-producing bacterium (SLC5A11位点的变异与一种产丁酸盐的细菌相关)
研究发现 SLC5A11基因座的变异与产丁酸盐的细菌 Agathobaculum butyriciproducens的丰度相关。该基因编码钠/肌醇协同转运蛋白2（sodium/myo-inositol cotransporter 2），介导肠道顶端肌醇的吸收。研究中确认了该位点的lead variant rs55808472-A等位基因与较低的血清肌醇水平相关。该位点内的另一个基因 ARHGAP17编码的蛋白参与维持紧密连接和囊泡运输，其缺陷会增加肠道通透性。这一发现将宿主的营养吸收、肠道屏障功能与一种具有潜在神经保护作用的产丁酸菌联系了起来。
Loci associated with microbial functions suggest genetic links to microbial carbohydrate and amino acid catabolism (与微生物功能相关的基因座提示遗传与微生物碳水化合物和氨基酸分解代谢的联系)
在微生物功能层面，研究虽然没有发现达到研究水平显著性的关联，但在基因组水平上识别出11个候选基因座与11个微生物功能模块相关，这些功能大多涉及碳水化合物和氨基酸的分解代谢，表明宿主遗传也可能影响微生物群落的代谢潜能。
研究结论与讨论
本研究通过大规模、高分辨率的宏基因组GWAS，极大地扩展了我们对宿主遗传如何调控肠道菌群组成的理解。研究不仅首次确立了肠内分泌脂肪酸化学感应（通过OR51E1–OR51E2位点）与微生物组丰富度的遗传关联，还将稳健复现的微生物组相关基因座从两个（ABO和 LCT）增加到了六个（新增 FUT2、MUC12、CORO7–HMOX2和 SLC5A11），并为另外两个位点（FUT3–FUT6和 FOXP1）提供了强有力的支持证据。这些发现共同指向了胃肠道关键生理机制在微生物组调控中的核心作用：肠内分泌细胞的脂肪酸化学感应、胆汁组成、黏膜层构成以及细胞表面聚糖的呈现和分泌。
研究的优势在于跨队列的数据处理标准化、严格的统计校正以及独立队列的稳健复现。局限性则包括参与者主要来自北欧的欧洲血统人群，限制了结果的普适性，以及对稀有变异或低丰度物种的检测效力有限。此外，如同大多数GWAS研究一样，确定因果基因和阐明具体的分子机制仍然是未来的挑战。
总之，这项工作标志着在解析宿主遗传对肠道微生态影响的道路上迈出了重要一步。它强调，未来的研究应继续关注与胃肠道生理相关的基因，并考虑已知影响微生物组的因素（如抗生素、药物和饮食）。这些发现为未来开发同时考虑宿主遗传特征和微生物组谱的个体化治疗或干预策略奠定了重要的科学基础。
