牛与山羊瘤胃微生物对饲粮抑制甲烷生成的物种特异性响应

《Frontiers in Microbiology》：Species-specific rumen microbial responses to dietary inhibition of methanogenesis in cows and goats

　　

引言：甲烷（CH4）作为一种温室气体，由反刍动物瘤胃微生物组产生。多种营养策略可降低家畜肠道CH4排放，但饲粮是否以相似方式影响不同反刍动物物种的瘤胃微生物组尚不清楚。本研究旨在确定牛和山羊在对基于淀粉和/或脂质补充的饲粮策略的瘤胃微生物功能响应上是否存在差异，以及这些差异能否解释可比实验条件下与饲粮相关的CH4排放变异。方法：实验同时开展，两物种均接受以草地干草和精料为基础的控制饲粮（CTL），或补充玉米油与小麦淀粉（COS）、海藻粉（MAP）或氢化棕榈油（HPO）。为从瘤胃微生物宏转录组中识别具有生物学意义的特徵，研究人员遵循五步整合统计与网络分析流程，结合网络方法、聚类与监督模型拟合，将差异表达基因和分类单元与瘤胃CH4排放关联。结果：COS饲粮在牛和山羊中均降低CH4排放并改变瘤胃微生物组。然而COS对比CTL中鉴定出的差异表达KEGG同源物（KO）和差异丰度OTU数量，牛多于山羊。此外，聚类分析揭示两物种间网络拓扑存在差异。在山羊中，CH4减排与碳水化合物代谢及甲烷生成的甲基营养型和氢营养型途径相关基因强关联；而在牛中，氢营养型途径占主导。另外，稀疏偏最小二乘（sPLS）分析识别出物种特异性判别微生物特徵。讨论：总体而言，结果表明宿主物种调控瘤胃微生物对饲粮的响应，并提示牛与山羊间CH4减缓背后的微生物互作存在差异，尽管需在更大样本量的独立研究中验证。

研究背景与立项依据：反刍动物凭借瘤胃内专性微生物组将植物生物质转化为高营养食物，但肠道甲烷（CH₄）排放作为强效温室气体带来高碳足迹压力。奶牛是产奶系统主要CH₄源，而奶山羊、绵羊与水牛亦有贡献，且山羊因适应严酷气候其种群与产区可能扩张。瘤胃甲烷由产甲烷古菌经氢营养型（hydrogenotrophy）为主、甲基营养型（methylotrophy）为辅的途径生成，营养策略（如淀粉与脂质补充）可缓解排放，但既往研究多聚焦单一物种（尤指牛），关于跨反刍物种微生物信号是否保守或受物种生理/采食行为/食糜周转塑造知之甚少。Martin等（2021）前期发现玉米油+小麦淀粉（COS）饲粮在奶牛与山羊中均降CH₄，而海藻粉（MAP）与氢化棕榈油（HPO）无显著效应，为本研究提供同条件跨物种比较契机。研究人员假设宿主物种会因瘤胃环境差异导致微生物互作与产甲烷途径关联不同，进而采用整合宏转录组与表型性状分析验证。

主要技术方法概览：研究在法国Herbip?le设施同期开展牛（4头荷斯坦奶牛）与山羊（4头阿尔卑斯山羊）各一个4×4拉丁方试验，每物种每饲粮4个观测值共32份瘤胃液样本；饲粮为CTL、CTL+5%DM玉米油与小麦淀粉（COS）、CTL+1.5%DM裂壶藻粉（MAP）、CTL+3%DM氢化棕榈油（HPO）。表型端在呼吸测热舱测定CH₄产量/产率/强度及CH₄/CO₂，同步测挥发性脂肪酸（VFA）、氨态氮、原虫计数；核酸端取瘤胃液液氮速冻提总RNA行Illumina NovaSeq 4000宏转录组测序，序列经MetaTrans流程分rRNA/tRNA与mRNA部分，分类用OTU（99%相似聚团对SILVA v138.1注释），功能用非rRNA读段比对牛瘤胃基因目录并赋KEGG同源物（KO）。统计端构建五步管线：DESeq2做差异丰度（deKO/deOTU）、ccrepe算Pearson共现网络（Cytoscape可视化）、K-means聚类、KEGGREST通路富集、mixomics包sPLS回归关联微生物簇与CH₄表型，全部按物种分别分析。

研究结果

5.1 牛与山羊瘤胃微生物组宏转录谱：CTL下两物种共检出5235 KO与368 OTU；优势门均为厚壁菌门（Firmicutes，牛45%/山羊54%）与拟杆菌门（Bacteroidetes，牛24%/山羊22%）。PCoA显示OTU谱物种聚类明显（anosim R=0.7917, p=0.033），KO谱无分离（R=0.3021, p=0.057），说明分类结构有别而功能活性相似。

5.2 控制与补充饲粮下比较：CTL牛vs山羊仅16个deOTU，山羊富集Methanobrevibacter与VadinCA11（Methanomassiliicoccales），牛富集Fibrobacter succinogenes与Sharpea；400个deKO中多数低丰度。仅COS相关比较（COS vs CTL/MAP/HPO）出现显著deKO/deOTU，与CH₄降约28%一致。牛COS特异deKO 255个、deOTU 10个，山羊仅22个deKO、4个deOTU；牛网络（COS vs CTL）10902边，簇1与5富产甲烷途径（如甲基辅酶M还原酶mcr、四氢甲烷蝶呤S-甲基转移酶mtr、甲酰甲烷呋喃脱氢酶fwd），山羊网络仅518边、5簇，簇1/5关CH₄且日产量与产率同簇；通路富集示牛簇1甲烷代谢KO占16%，山羊簇5甲烷代谢KO占8%且多碳水代谢。

5.3 候选微生物标记：sPLS于牛簇1留25变量（含8个甲烷代谢KO如K00320/K00579/K03422），山羊簇5留30变量（含8个甲烷代谢KO如K00202/K04480/K14081），两物种仅共享3个判别deKO（K00320、K03422、K00402）。牛簇1解释CH₄日产量变异58.7%，山羊簇1解释84%。

讨论部分总结：研究人员指出虽COS在两物种表型减排幅度相近，但支撑表型的微生物功能响应不同——牛靠氢营养型通路广受扰动、KO级联响应强，山羊则甲基营养型（Mta甲醇辅酶M甲基转移酶）特异浮现且碳水代谢权重更大，这可能与山羊采食快、瘤胃周转高、晨采样前可发酵底物少有关，但采食行为未做因果检验。HPO/MAP低脂添加未改KO表达，说明仅高淀粉+脂质协同触发代谢重编程。Andersen等（2023）同批样本多组学佐证COS将氢导向琥珀酸/丙酸菌而非产甲烷弧菌。作者强调4×4拉丁方每物种n=4致高维分析属探索性，需大队列验证。

结论译文：本研究在相同饲粮与实验条件下整合CH₄排放性状与宏转录组，直接比较牛和山羊瘤胃微生物组；此前山羊瘤胃微生态刻画薄弱。结果显示两物种CH₄产率均与碳水化合物代谢基因紧密关联，凸显饲粮驱动底物可用性对产甲烷之重要；尽管COS下CH₄减排相似，牛与山羊产甲烷相关通路的候选微生物标记各异，提示两物种与饲粮诱导CH₄减缓关联的微生物特徵可能存在差异。但因高维分析属探索性且主要种内分析，结果应视为物种关联模式而非确定机制；需更大队列与专门检验种间差异设计以验证候选标记，并判定瘤胃生理、采食行为或食糜动力是否因果贡献于CH₄减缓微生物响应。论文发表于《Frontiers in Microbiology》。