编辑推荐：
人类肠道微生物组成的改变已被证实与多种慢性疾病相关。然而，仅依靠物种水平的组成分析会掩盖物种内广泛的遗传多样性与功能变异。这推动了学界向菌株解析、遗传层面的微生物组分析转变，即利用微生物遗传变异来识别与疾病相关的基因与通路。研究人员在6个独立队列共10,781
人类肠道微生物组成的改变已被证实与多种慢性疾病相关。然而，仅依靠物种水平的组成分析会掩盖物种内广泛的遗传多样性与功能变异。这推动了学界向菌株解析、遗传层面的微生物组分析转变，即利用微生物遗传变异来识别与疾病相关的基因与通路。研究人员在6个独立队列共10,781名个体中，对1,280万个位点水平的肠道微生物单核苷酸变异（Single-Nucleotide Variant，SNV）与244项宿主健康表型及环境暴露因子开展了全表型关联研究（Phenome-Wide Association Study，PheWAS）。研究共鉴定出940个独立关联信号，覆盖69项宿主因素，包括年龄、性别以及心血管与代谢参数。结构分析显示，与宿主表型相关的SNV位于蛋白质功能重要区域，为微生物遗传改变如何参与环境与宿主生理因素提供了机制假说。例如，研究人员发现普拉梭菌（Faecalibacterium prausnitzii）中存在一个与年龄相关的SNV，影响结合谷胱甘肽（glutathione）的蛋白质口袋内的氨基酸。为支持后续研究，所有汇总统计量均提供于在线目录。总之，这些结果在碱基对分辨率上识别了宿主与微生物遗传变异之间新的生物学联系，并为假设生成提供资源，助力理解肠道微生物组的菌株水平功能。
研究背景与立项依据
人类肠道微生物组的组成与遗传景观具有高度变异性，物种内遗传变异决定了菌株在适应度、代谢能力、致病性等层面的差异。既往微生物组研究多停留在物种丰度或结构变异层面，忽视了核心基因组中大量单核苷酸变异（SNV）对微生物功能与宿主互作的调控作用。已有个例研究表明，大肠杆菌blc基因、幽门螺杆菌丝氨酸蛋白酶及拟杆菌糖苷水解酶相关SNV可分别影响屏障功能、胃癌风险与2型糖尿病，但缺乏在全表型尺度上系统解析微生物SNV与宿主表型关系的研究。同时，微生物群体遗传结构复杂，若不充分控制谱系效应，关联分析极易产生假阳性。为此，研究人员在《Nature》报道了这项跨队列、碱基对分辨率的SNV水平全表型关联研究，旨在建立宿主表型—微生物遗传变异的系统关联图谱，并为机制假说生成提供资源。
主要技术方法
研究人员整合荷兰5个队列与坦桑尼亚1个队列共11,831份粪便宏基因组测序样本，覆盖10,781名个体。利用GT-Pro进行基于参考面板的代谢基因型（metagenotyping）判定，保留核心基因组中1,280万个常见双等位SNV。针对每个物种、表型与队列，采用线性混合模型（Linear Mixed Model，LMM）将物种内样本间遗传连锁矩阵作为随机效应以控制种群结构，并将替代等位基因频率编码为0至1的剂量值。各队列独立分析后，按样本量加权进行荟萃分析，再通过连锁不平衡分组提取独立信号，并以置换检验评估跨队列稳健性。功能解读上结合SNV注释、AlphaFold与ESM结构预测及分子对接。
研究结果
肠道菌群SNV景观
研究人员在881个物种核心基因组中检出超5,090万SNV，过滤后保留433个物种的约1,270万常见SNV，其中93.02%位于蛋白编码区。不同队列SNV密度与物种覆盖度存在差异，坦桑尼亚队列物种与SNV丰富度最高。该部分为后续关联分析建立了高分辨率的微生物基因型基线。
关联信号图谱
通过LMM荟萃分析并经过置换检验，研究人员获得49,346条显著且稳健的关联，涉及101个物种的39,663个SNV与69项表型；经连锁不平衡分组去冗余后定义940个独立信号。其中65.74%的信号在物种丰度层面无显著关联，说明SNV分析可提供超出丰度与结构变异的互补信息。年龄信号最多，达99个；部分信号在荷兰与坦桑尼亚队列间可重复，提示跨人群收敛性。
年龄相关SNV提示谷胱甘肽转运机制
年龄是最强关联表型。研究人员发现普拉梭菌_I中sarA_2基因（oppA同源基因）内一组错义SNV与年龄显著相关，携带替代等位基因个体平均年龄高约4.9年。结构建模显示该蛋白属ABC型寡肽通透酶（OPP）复合体，分子对接表明其口袋可结合谷胱甘肽。年龄相关突变可能通过改变转运效率，帮助菌株抵御老化肠道的氧化应激。
错义SNV的结构属性
针对6,799个关联错义SNV，研究人员系统注释其相对溶剂可及性、至配体距离（DTL）与至蛋白互作界面距离（DTPPI）。752个SNV位于结构重要区，涉及40种表型。除OPP外，年龄相关SNV还富集于serpin结构域蛋白及marR同源转录调控二聚体界面，计算机丙氨酸扫描提示其可改变二聚化自由能，进而影响氧化应激响应。
讨论与结论翻译
过去十年，全基因组关联研究（GWAS）与基于宏基因组学的微生物组研究增进了人们对宿主遗传与共生微生物组的认识，但将GWAS应用于宏基因组数据仍有限。本研究利用大规模深度表型队列，通过肠道微生物SNV水平全表型关联研究，系统刻画了244项表型因子，识别出101个肠道微生物物种SNV与69项人类表型及暴露间的940个独立荟萃分析关联信号，全部结果收录于在线目录。为克服低丰度物种分型缺失与菌株混合导致的群体结构偏差，研究人员过滤低信息物种，保留准可定相（quasi-phaseable）样本以构建基于亲缘关系的种群结构。在940个信号中，51.3%含至少一处错义改变；大量编码区同义SNV可能具有调控或作为因果变异代理的作用。本研究提供的微生物SNV目录是验证与假说生成的重要资源。但样本以荷兰人群为主，存在胜者诅咒风险；未来需在多祖先人群中扩展，并结合宏转录组、宏蛋白组、代谢组与培养验证，从关联信号推进至机制结论。
研究结论如下：本研究在碱基对分辨率上识别了宿主与微生物遗传变异之间新的生物学联系，并提供可用于假说生成的资源，将助于未来理解肠道微生物组的菌株水平功能。
