Urochloa brizantha提取物对绵羊瘤胃发酵、产甲烷古菌多样性及代谢组学的影响

《Metabolites》：Ruminal Fermentation, Methanogenic Archaeal Diversity, and Metabolomic Responses to Urochloa brizantha Extracts in Sheep

　　

背景：植物源性饲料添加剂正被探索作为传统瘤胃调节剂的替代品，然而仅通过常规发酵参数评估时，其对瘤胃生态系统的影响仍未得到充分表征。本研究评价了Urochloa brizantha乙醇和水-乙醇提取物对绵羊瘤胃发酵、产甲烷古菌多样性及瘤胃代谢组的影响。方法：八只装有瘤胃瘘管的去势公羊被分配到重复的4 × 4拉丁方设计中，包含四个28天周期。处理组包括无添加剂组（CTL）、U. brizantha乙醇提取物组（EE；50 mL/天）、水-乙醇提取物组（HE；80 mL/天）和莫能菌素组（MON；25 mg/kg DMI）。对瘤胃液进行了发酵性状、古菌16S rRNA基因多样性和非靶向代谢组学分析。结果：干物质采食量和发酵参数未受处理影响。古菌α多样性表现出依赖于对比组的响应，接受植物提取物的绵羊其Shannon多样性低于接受莫能菌素的绵羊（p = 0.034），且Methanimicrococcus在EE组中被选择性富集（p = 0.012）。非靶向代谢组学呈现最清晰的响应：尽管整体PCA图谱存在大量重叠，PLS-DA仅在EE处理的绵羊中显示出第0天与第21天之间处理特异性的时间分离，分类错误率的置换检验支持该结果（p = 0.027），但预测性能应谨慎解读。该EE特异性响应以羧酸亚类发生显著变化为特征（FDR = 0.022）。VIP排名靠前的判别代谢物（天冬氨酸、琥珀酸、谷氨酸、苯丙氨酸、亮氨酸和胆碱）呈名义上的下降，但未通过多重检验校正（FDR = 0.0602–0.2580）。CTL、HE或MON均未显示出时间分离。结论：Urochloa brizantha提取物对常规瘤胃发酵参数产生有限影响，但与选择性古菌响应相关。特别是乙醇提取物与处理特异性的时间代谢组模式相关，包括羧酸亚类的显著变化，而单个代谢物的变化仍属探索性结果，需进一步研究。

论文解读：Urochloa brizantha提取物对绵羊瘤胃发酵、产甲烷古菌多样性与代谢组学响应的影响

瘤胃发酵虽可将纤维性生物质转化为高价值畜产品，但伴随以甲烷形式损失日粮能量的固有低效性，而甲烷同时是重要温室气体。莫能菌素等离子载体作为饲料添加剂被广泛用于改善瘤胃能量效率，但消费者对天然畜产品的需求及对抗菌添加剂的监管趋严推动了植物源性添加剂的研发。植物次生代谢物如皂苷等具有调节瘤胃微生物群和发酵谱的潜力，然而其作用机制尚不充分明确。既往研究多仅依赖pH、氨态氮和总短链脂肪酸等常规发酵终点评估此类添加剂，可能无法反映微生物生态和代谢的功能性变化。由于瘤胃微生物组具有显著功能冗余，常规参数变动有限时仍可能发生生物学重要改变。Urochloa brizantha作为巴西广泛种植的饲草，含有甾体皂苷（如原薯蓣皂苷）等活性化合物，体外试验显示其提取物可调节瘤胃发酵，但其体内效应未明。因此本研究旨在综合评估其乙醇和水-乙醇提取物对绵羊瘤胃发酵常规参数、产甲烷古菌多样性和瘤胃代谢组的影响，并与莫能菌素和对照比较。研究论文发表于《Metabolites》。

研究人员采用8只装有瘤胃瘘管的Santa Inês × Dorper去势公羊（体重41.1 ± 8.1 kg），按重复4 × 4拉丁方设计开展4个28天周期的实验。处理包括对照（CTL，无添加剂）、U. brizantha乙醇提取物（EE，50 mL/天）、水-乙醇提取物（HE，80 mL/天）和莫能菌素（MON，25 mg/kg DMI）。提取物通过加压液体萃取获得，EE和HE原薯蓣皂苷浓度分别为7.50和7.48 mg/mL，每日供给量分别为375.0和598.4 mg原薯蓣皂苷。日粮精粗比为30:70，试验周期包括第0天基线采样、第1至21天添加剂补充和第22至28天洗脱期。第21天采集瘤胃液用于发酵参数分析、古菌16S rRNA V8区扩增子测序和1H NMR非靶向代谢组学分析；代谢组学另在第0天采集周期特异性基线样本。发酵参数包括pH、氨态氮和短链脂肪酸；采用DADA2流程处理序列并对SILVA数据库进行分类；使用MetaboAnalyst和R进行PCA、PLS-DA、VIP分析及化学分类学亚类比较；相关性分析采用Spearman秩相关。统计分析使用SAS GLIMMIX，设有先验正交对比。

干物质采食量、瘤胃pH、氨态氮、总短链脂肪酸浓度及乙酸盐、异丁酸盐和丁酸盐摩尔比例在各处理间无显著差异。正交对比显示，MON组氨态氮趋于低于植物提取物组（p = 0.093），EE组氨态氮趋于高于HE组（p = 0.083）；MON组丙酸盐比例趋于高于植物提取物组（p = 0.067）；添加剂组异戊酸盐比例（p = 0.062）和乙酸/丙酸比（p = 0.094）趋于高于对照组；EE组戊酸盐比例趋于高于HE组（p = 0.070）。

古菌α多样性无整体处理效应，但对比显示植物提取物组Shannon多样性低于MON组（p = 0.034），Simpson多样性呈相同趋势（p = 0.065）。属水平上，Methanimicrococcus是唯一受处理显著影响的属（p = 0.012，Bonferroni校正后p = 0.048），其在EE组相对丰度（1.63%）显著高于CTL（0.46%）、HE（0.60%）和MON（0.07%）。主要氢营养型属Methanobrevibacter以及Methanosphaera和Methanoplanus无处理相关差异。

代谢组学PCA显示所有样本在整体谱图上大量重叠，前两主成分解释82.32%方差。针对单个处理的PLS-DA仅显示EE组第0天与第21天样本间存在部分分离（四组分模型R2X = 0.713，R2Y = 0.921，Q2 = 0.504）；分类错误率置换检验显著（p = 0.027），但Q2置换检验不显著（p = 0.058），故该时间分离被解释为探索性结果。CTL、HE和MON均无显著时间区分。化学分类学分析显示，EE组羧酸亚类占瘤胃代谢组的比例由基线80.57%显著升至第21天的85.57%（p = 0.002，FDR = 0.022）；碳水化合物、二羧酸及其衍生物、氨基酸肽及类似物等亚类呈数值下降，但未达FDR显著性。VIP排名前列的判别代谢物为天冬氨酸（VIP = 1.91）、琥珀酸（VIP = 1.82）、谷氨酸（VIP = 1.79）、苯丙氨酸（VIP = 1.73）、苏氨酸（VIP = 1.53）和亮氨酸（VIP = 1.45）。EE补充后第21天相对基线，天冬氨酸（-49.4%）、琥珀酸（-54.0%）、谷氨酸（-41.0%）和苯丙氨酸（-37.9%）呈名义显著下降（原始p < 0.01），但FDR校正后未达显著（FDR = 0.0602）；亮氨酸和胆碱的数值下降亦未通过FDR校正。第21天高丰度代谢物比较显示，MON组谷氨酸、麦芽糖、葡萄糖和黄嘌呤浓度倾向于高于植物提取物组，CTL组亮氨酸趋于高于添加剂组。

EE组内的相关性分析显示，氨态氮与Methanobrevibacter负相关（r = -0.74，p = 0.046）、与Methanoplanus正相关（r = 0.77，p = 0.027）；戊酸盐比例与Methanimicrococcus强负相关（r = -0.81，p = 0.027）。代谢物-发酵性状关联中，氨态氮与2-苯基丙酸盐正相关（r = 0.89，p = 0.003），丁酸盐与丙酮（r = 0.87，p = 0.011）和琥珀酸（r = 0.86，p = 0.013）正相关，异丁酸盐与苯丙氨酸、甘油酸、尿苷、黄嘌呤和尿嘧啶正相关，异戊酸盐与2-羟基异丁酸盐负相关。所有相关性均未通过FDR校正，仅作为假设生成依据。

讨论部分指出，研究核心发现是常规发酵测量与多组学揭示的功能响应明显分离。原薯蓣皂苷日摄入量（0.375–0.598 g/kg DM）远低于通常用于体内甲烷减排研究的皂苷浓度（可达5.0 g/kg DM），故而常规参数稳定。Methanimicrococcus的选择性富集与其甲基营养型代谢特征相关，该菌缺乏完整Wood–Ljungdahl途径，依赖甲胺类和甲醇等甲基化合物。EE与HE的差异受提取溶剂极性和皂苷剂量双重混淆，无法单独归因于溶剂效应。代谢组学中天冬氨酸、谷氨酸和琥珀酸的下降提示微生物中间代谢和交叉互养的变化，而非发酵活性受损。莫能菌素显示广谱抗菌机制，而EE提取物呈现选择性代谢调节特征。甾体皂苷的膜相互作用可能是潜在机制，但直接证据有限。局限性包括单时间点采样、提取物比较存在剂量混杂、短扩增子分辨率有限、未直接定量甲胺及代谢通量等。未来需结合稳定同位素示踪、靶向代谢组学和宏转录组学验证。

结论部分指出，乙醇提取物与瘤胃中间碳氮代谢差异及甲基营养型古菌Methanimicrococcus相对富集相关，其模式更符合靶向生物响应而非广泛微生物干扰。由于提取溶剂、投喂体积和原薯蓣皂苷剂量同时存在差异，提取极性对植化生物活性的具体贡献需在未来剂量匹配研究中确定。本研究凸显了将代谢组学与微生物群落分析相结合以检测常规发酵测量难以发现的响应模式的价值，并为后续剂量-效应和体内验证研究提供了基础。