《mSystems》：Host genotype-microbiota-metabolome interactions modulate semen quality in Holstein bulls
编辑推荐：
精液品质是奶牛繁殖性能的关键决定因素，然而其微生物、代谢与遗传基础仍不清楚。研究人员整合了55头按精液品质记录划分为相对高质（RSQH）与相对低质（RSQL）组的荷斯坦公牛的瘤胃微生物组、代谢组与全基因组变异数据，以探究与繁殖力相关的宿主-微生物组-代谢组互作
精液品质是奶牛繁殖性能的关键决定因素，然而其微生物、代谢与遗传基础仍不清楚。研究人员整合了55头按精液品质记录划分为相对高质（RSQH）与相对低质（RSQL）组的荷斯坦公牛的瘤胃微生物组、代谢组与全基因组变异数据，以探究与繁殖力相关的宿主-微生物组-代谢组互作。瘤胃微生物主要由拟杆菌门（Bacteroidota）、厚壁菌门（Firmicutes）与变形菌门（Proteobacteria）主导，其中纤维杆菌属（Fibrobacter，纤维杆菌门）在RSQH公牛中富集，并与精子活力及密度呈正相关。代谢组可区分两组，共鉴定64个判别代谢物，RSQH组中N-乙酰-L-苯丙氨酸与前列腺素H2升高，RSQL组中癸二酸（sebacic acid）升高。研究人员利用机器学习模型优先排序微生物与代谢特征，识别出纤维杆菌属、血栓烷B2（thromboxane B2）与N1-[4-(氰甲基)苯基]-4-氯苯甲酰胺可作为精液品质候选指示物，上述指标得到机器学习、高遗传力与全基因组关联分析（GWAS）信号支持。总体鉴定到511个单核苷酸多态性（SNP）共同关联纤维杆菌属与上述代谢物，其中6个SNP经GWAS分析与精子活力显著相关。这些SNP分为两类相反趋势：C1型变异在精液品质相关性状中呈下降趋势，C2型变异呈上升趋势。结构方程建模（SEM）显示，C2型变异通过微生物与代谢物富集正向关联精液性状，而C1型变异则呈负向关联。该结果凸显了宿主遗传、微生物与代谢对荷斯坦公牛精液品质的协同贡献，并为繁殖力相关研究及育种提供了候选生物标志物。
研究背景与立项依据
精液品质直接决定奶牛育种产业的受精成功率、群体繁殖效率与经济回报，传统研究多聚焦遗传因子，对瘤胃微生物组（rumen microbiota，栖息于瘤胃内的细菌、古菌、真菌等微生态集合）、代谢组与宿主基因组如何协同调控精子发生与精液性状仍缺乏系统认知。已有线索提示肠道-睾丸轴可能通过能量代谢、免疫与激素途径影响生殖健康，但公牛中“宿主遗传-微生物-代谢物”三级联动的机制尚未阐明，难以支撑可落地的育种生物标志物开发。
研究人员以55头4至6岁、标准化全混合日粮（TMR）饲喂的山东某育种场荷斯坦公牛为队列，按精子活力与密度为主划分为相对高质（RSQH）与相对低质（RSQL）两组，整合瘤胃16S rRNA扩增子测序、液相色谱串联质谱（LC-MS/MS）代谢组、Illumina BovineSNP150芯片基因分型，开展多组学关联与因果推断，旨在识别可遗传的微生物与代谢候选标志物，并解析其遗传调控路径。
主要技术方法概览
研究人员构建单中心公牛样本队列，依次开展瘤胃内容物的16S rRNA（V3-V4区）扩增子测序与ASV（amplified sequence variant，扩增子序列变体）分类学注释；采用LC-MS/MS非靶向代谢组检测并做分位数归一化；基于BovineSNP150芯片（UMD3.1参考基因组）进行质控后SNP分析；以随机森林（Random Forest）与极端梯度提升（XGBoost）做分组分类与特征优先排序；用GCTA（genome-wide complex trait analysis，全基因组复杂性状分析）估计窄义遗传力（h2）；通过mGWAS（microbiome/genome-wide association study，微生物或代谢物全基因组关联分析）定位共享SNP；最后用结构方程模型（SEM）评估基因型-微生物/代谢物-精液表型的直接与间接效应。
研究结果浓缩
瘤胃微生物组成表征
通过对55头公牛瘤胃16S测序，研究人员发现群落以Bacteroidota、Firmicutes、Proteobacteria为优势核心门，属水平以Prevotella与Rikenellaceae_RC9_gut_group为主。α多样性两组无差异，β多样性明显分离，说明精液质量差异伴随微生物结构重组而非单纯丰富度变化。
微生物多样性与差异类群关联精液性状
差异丰度分析鉴定34个显著类群，Fibrobacter在RSQH组富集，且与精子活力（r_s=0.307）和密度（r_s=0.276）正相关；RSQL组富集未培养拟杆菌与环细菌科相关类群。共现网络显示RSQH中Fibrobacter与Ruminobacter连接更紧密，提示高质组具更协同的纤维降解与能量代谢网络。
代谢组差异与功能指向
非靶向代谢组鉴定64个差异代谢物：RSQH组N-乙酰-L-苯丙氨酸、前列腺素H2升高，RSQL组癸二酸、sulcatol升高。相关性显示癸二酸负向关联精子密度，N-乙酰-L-苯丙氨酸正向关联活力与密度，提示氨基酸衍生物与脂质介质参与维持精子膜稳定与抗氧化防御。
机器学习分类与标志物优先排序
以差异微生物或代谢物训练Random Forest与XGBoost，训练集AUC近1.0，测试集微生物模型AUC达0.953，代谢模型达0.951。跨模型交集锁定11个关键菌属与11个关键代谢物，Fibrobacter、血栓烷B2、N1-[4-(氰甲基)苯基]-4-氯苯甲酰胺被优先判为精液质量判别标志。
宿主遗传控制与共享SNP轴
遗传力估计显示Fibrobacter（h2=0.63）、unclassified_Alphaproteobacteria（h2=0.66）及多个代谢物具中高遗传力。mGWAS发现Fibrobacter与N1-[4-(氰甲基)苯基]-4-氯苯甲酰胺共享468个SNP，与血栓烷B2共享43个SNP，合计511个。注释基因涉及PHGDH、PTPN13、PRKG1、MAPK10等，富集于能量代谢、氨基酸合成、线粒体蛋白加工与精子信号通路；基因型频率分析显示0/0纯合型偏向RSQH，1/1或0/1偏向RSQL。
SNP-微生物-代谢轴与因果建模
SEM将SNP分为C1与C2两类：C2型通过提升Fibrobacter及代谢物丰度正向间接作用于精子活力，C1型呈反向直接/间接效应。总效应为直接路径β?与间接路径α?×β?之和，证明宿主遗传可经瘤胃微生态中位变量塑形精液表型。
讨论与结论翻译
讨论指出，该研究首次在公牛中以多组学联动方式证明精液品质并非单一遗传性状，而是宿主基因型、瘤胃可遗传微生物（如Fibrobacter）与代谢介质（如血栓烷B2、特定苯甲酰胺类分子）共同调度的结果。共享SNP大多位于内含子与基因间区，说明调控多通过表达定量或远端增强机制发挥作用，需在独立大群体中验证。
研究结论如下：荷斯坦公牛精液品质受宿主遗传-瘤胃微生物组-代谢组协同调控；Fibrobacter、血栓烷B2与N1-[4-(氰甲基)苯基]-4-氯苯甲酰胺构成高遗传力候选生物标志物组合；C2型SNP经微生物与代谢物富集促进精液性状，C1型SNP起抑制作用；该框架为基于遗传选育与微生物信息辅助的公牛繁殖力评估及育种策略提供了理论与标记储备。
