低剂量Asparagopsis taxiformis联合3-硝基氧基丙醇通过微生物与酶学调控实现瘤胃甲烷生成近完全抑制

《Journal of Animal Science and Biotechnology》：Near-complete inhibition of rumen methanogenesis via microbial and enzymatic modulation using a low dose of Asparagopsis taxiformis combined with 3-nitrooxypropanol

　　

反刍动物肠道甲烷（CH₄）排放是农业温室气体排放的主要来源之一。红海藻Asparagopsis taxiformis（A. taxiformis）是一种高效的CH₄排放抑制剂，但其大规模应用受限于生物量不足。本研究旨在评估降低A. taxiformis添加水平（0.32%干物质，DM）并联合3-硝基氧基丙醇（3-nitrooxypropanol，3-NOP；0.05% DM）是否仍能维持高抑制效率，并通过体外发酵与宏基因组学分析阐明其潜在的微生物学机制。联合处理使CH₄产量降低98.21%（P<0.01），且未损害干物质降解率，同时显著促使瘤胃发酵向丙酸方向转变，降低乙酸/丙酸比值（1.59 vs. 2.65；P<0.01）。宏基因组学分析显示，Methanobrevibacter和Ruminococcus的丰度显著降低，而与丙酸相关的细菌如Prevotella、Treponema、Eubacterium和Selenomonas的丰度显著增加（P<0.01）。在功能层面，联合处理下调了氢营养型和甲基营养型产甲烷途径中的关键酶，包括甲基辅酶M还原酶（methyl-coenzyme M reductase，EC:2.8.4.1）和四氢甲烷蝶呤S-甲基转移酶（tetrahydromethanopterin S-methyltransferase，EC:2.1.1.86），从而阻断了产甲烷的终末步骤。综上所述，A. taxiformis与3-NOP联合添加可在大幅降低海藻用量的同时实现产甲烷的近完全抑制，该效应通过发酵模式、微生物群落结构和产甲烷酶学途径的协同改变而实现。这一策略为克服A. taxiformis的生物量限制提供了实用方案，并需在长期体内试验中进一步验证。

大气温室气体（greenhouse gas，GHG）浓度持续升高是全球变暖的主要驱动因素，其中甲烷（CH₄）是仅次于二氧化碳（CO₂）的第二大影响力温室气体。反刍动物瘤胃发酵产生的CH₄约占全球人为温室气体排放的14.5%以及畜牧业排放的40%，因此开发有效的减排策略对缓解全球变暖至关重要。在众多干预措施中，红海藻Asparagopsis taxiformis（A. taxiformis）作为一种高效饲料添加剂备受关注，研究表明添加≤5% DM的A. taxiformis可抑制高达99%的CH₄排放。其主要机制在于生物活性化合物（如溴仿）作为甲基辅酶M类似物，竞争性抑制产甲烷终末步骤的关键酶——甲基辅酶M还原酶（methyl-coenzyme M reductase，MCR），同时改变瘤胃代谢通路，使氢气（H₂）从产甲烷转向丙酸发酵等替代途径。然而，A. taxiformis的自然分布局限于热带和暖温带海域，大规模养殖困难，生物量有限制约其广泛应用。另一广泛研究的合成甲烷抑制剂3-硝基氧基丙醇（3-nitrooxypropanol，3-NOP）同样靶向抑制MCR，但通常仅能使肠道甲烷排放降低26%-54%，效果远低于A. taxiformis。如何在保持高抑制效率的同时降低A. taxiformis用量，成为该领域亟待解决的问题。基于此背景，研究人员在《Journal of Animal Science and Biotechnology》发表了这项研究，旨在评估低剂量A. taxiformis联合3-NOP对CH₄排放和瘤胃发酵参数的影响，并通过宏基因组学方法阐明其微生物学机制。

研究中使用的A. taxiformis采集于中国广东省湛江市硇洲岛沿海水域，经清洗、冻干、研磨后保存；3-NOP购自上海阿拉丁生化科技股份有限公司。体外发酵实验采用两处理设计：对照组（CON，无添加剂）和联合处理组（AN，添加0.32% A. taxiformis和0.05% 3-NOP），每组6个重复，发酵48小时。瘤胃液来源于3头荷斯坦奶牛（3岁泌乳牛，平均体重550±30 kg），采集后过滤并与缓冲液按1:2混合。发酵底物为与奶牛饲粮相同的总混合日粮。气体产量及成分通过气相色谱测定；干物质降解率（dry matter degradation，DMD）通过尼龙袋法测定；挥发性脂肪酸（volatile fatty acid，VFA）浓度采用气相色谱分析。宏基因组学分析中，总DNA toolbar DNA提取后构建约350 bp片段的测序文库，在Illumina NovaSeq 6000平台行双端测序。原始数据经fastp质控、BWA去宿主后，使用MEGAHIT进行组装，Prodigal预测开放阅读框（open reading frame，ORF），CD-HIT构建非冗余基因集，SOAPaligner进行reads映射。 taxonomic和功能注释分别采用DIAMOND比对NCBI NR数据库和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes，KEGG）数据库。统计分析采用SPSS 27.0进行独立样本t检验。

气体产生与发酵参数方面，联合处理组（AN）的总产气量显著低于对照组（240.02±6.40 vs. 293.03±12.02 mL/g DM，P<0.01）；CH₄产量近乎完全抑制，降至0.48±0.24 mL/g DM，较对照组（26.89±1.56 mL/g DM）降低98.21%（P<0.01）；CO₂产量亦降低（161.31±12.05 vs. 199.02±10.99 mL/g DM，P<0.01）；同时H₂产量显著升高（8.58±0.64 vs. 0.09±0.03 mL/g DM，P<0.01）。DMD未受影响（85.08%±2.31% vs. 85.87%±2.13%，P=0.55）。总挥发性脂肪酸（total volatile fatty acid，TVFA）浓度降低（57.73±2.59 vs. 72.57±7.82 mmol/L，P<0.01），其中乙酸、丁酸和异丁酸浓度降低，而丙酸和戊酸浓度升高（P<0.01），导致乙酸/丙酸比值（acetate-to-propionate ratio，A:P ratio）从2.65±0.06降至1.59±0.03（P<0.01）。

微生物群落组成方面，细菌群落α多样性显示，AN组的ACE、Chao1和Simpson指数显著降低，Shannon指数显著升高（P<0.01），表明物种丰富度降低但均匀度增加。主坐标分析（principal component analysis，PCA）显示两组明显分离（P<0.01，分析相似性分析analysis of similarity，ANOSIM检验）。门水平上，Bacteroidota、Bacillota和Pseudomonadota为优势门；属水平上，Prevotella、Treponema和Ruminococcus为优势属。差异丰度分析显示，AN组Prevotella、Treponema、Eubacterium、Hallella、Oscillibacter、Selenomonas、Bacteroides和Clostridium丰度显著增加，而Ruminococcus和Ruminobacter丰度显著降低（P<0.01）。线性判别分析效应大小（linear discriminant analysis effect size，LEfSe）结果（LDA>2）进一步确认Prevotella、Treponema、Hallella和Succinatimonas在AN组显著富集。

古菌群落呈现相似趋势，ACE、Chao1和Simpson指数降低，Shannon指数升高（P<0.01）。PCA显示两组显著分离（P<0.01）。Euryarchaeota为绝对优势门，Methanobrevibacter、Methanocorpusculum和Methanosphaera为优势属。AN组Methanocorpusculum、Methanosphaera和Methanomethylophilus丰度显著增加，Methanobrevibacter丰度显著降低（P<0.01）；LEfSe分析显示Methanocorpusculum和Methanosphaera在AN组显著富集。

功能注释与通路分析方面，KEGG level 3通路比较显示，AN组中微生物代谢于多样环境、代谢通路、甲烷代谢、碳代谢、次级代谢物生物合成、糖酵解/糖异生的相对丰度显著降低（P<0.01），而原核生物碳固定通路、丙酮酸代谢、丙酸代谢、牛磺酸和次牛磺酸代谢的相对丰度显著升高（P<0.01）。在甲烷代谢通路中，氢营养型产甲烷途径的关键酶基因丰度显著降低，包括EC:1.8.98.6、EC:1.8.7.3、EC:1.2.7.12、EC:1.8.98.1、EC:2.1.1.86、EC:2.8.4.1、EC:2.3.1.101、EC:1.5.98.2（P<0.01）和EC:1.8.98.5（P<0.05）。乙酸裂解型产甲烷途径的酶基因（EC:2.7.2.1和EC:2.3.1.8）丰度显著增加（P<0.01）。甲基营养型产甲烷途径中，利用甲醇和甲胺的酶基因（EC:2.1.1.90、EC:2.1.1.246和EC:2.1.1.248）丰度显著降低（P<0.01），而EC:2.1.1.249和EC:2.1.1.250丰度显著增加（P<0.01）。

相关性网络分析揭示了关键瘤胃发酵参数、微生物物种与KEGG通路之间的关联模式。CH₄产量与Methanobrevibacter、Mogibacterium和Ruminococcus丰度呈显著正相关，与Treponema、Eubacterium、Prevotella、Hallella、Oscillibacter和Selenomonas丰度呈显著负相关（P<0.05）。乙酸浓度与Methanobrevibacter、Mogibacterium和Ruminococcus呈正相关，与Eubacterium、Prevotella和Treponema呈负相关（P<0.01）；丙酸浓度则呈现相反模式。Methanobrevibacter、Ruminococcus和Mogibacterium的功能贡献在AN组降低，而Prevotella、Treponema和Eubacterium的功能贡献在AN组升高，表明Methanobrevibacter丰度的降低是CH₄减少的主要驱动因素，Prevotella和Treponema丰度的增加也发挥重要作用。

讨论部分，研究人员首先指出联合处理实现了与5.00% A. taxiformis相当的甲烷抑制效果（98.21% vs. 98.53%），但海藻用量大幅降低，解决了生物量限制的瓶颈问题。DMD未受影响但VFA产生减少，表明发酵模式发生改变——乙酸和丁酸产生减少降低了H₂形成，从而限制产甲烷底物；丙酸适度增加但主要由乙酸降低驱动A:P比值变化，提示H₂向丙酸途径的 redirect 有限，不足以完全解释CH₄的大幅下降。值得注意的是，AN组H₂显著积累，但其最终命运尚不明确，这可能部分反映了体外系统的固有限制。

微生物群落重塑方面，Prevotella、Treponema、Eubacterium和Selenomonas等丙酸相关细菌的富集，以及Ruminococcus等纤维降解菌的减少，共同促成了发酵模式的转变。尽管Ruminococcus减少，DMD仍维持稳定，归因于Prevotella和Treponema等功能冗余微生物的补偿作用。Methanobrevibacter的抑制是CH₄减少的最直接机制；Methanosphaera和Methanocorpusculum的相对富集可能反映了古菌群落的适应性重组——前者通过甲基营养型途径以较低CH₄产率产甲烷，后者亦为"低CH₄排放者"。

酶学机制层面，联合处理不仅降低了Methanobrevibacter丰度以破坏氢营养型产甲烷途径，还下调了多种关键酶基因，包括上游的甲酰甲烷呋喃脱氢酶（EC:1.2.7.12）和终末步骤的四氢甲烷蝶呤S-甲基转移酶（EC:2.1.1.86）及MCR（EC:2.8.4.1），实现了从上游到终末步骤的多层次抑制。这种"双管齐下"的酶学抑制可能部分解释了低剂量联合处理为何能达到与高剂量A. taxiformis相当的抑制效果。A. taxiformis的卤代代谢物（如溴乙酸、溴酚衍生物）与3-NOP均靶向MCR，协同作用增强了抑制效率。

研究人员也指出研究的局限性：体外发酵系统的H₂积累可能高于体内实际情况；瘤胃微生物生态系统具有高恢复力，长期体内效果有待验证；未涉及原生动物和厌氧真菌等其他重要微生物类群，限制了对瘤胃微生物生态系统的全面解读。

研究结论表明，A. taxiformis（0.32% DM）与3-NOP（0.05% DM）联合添加可实现体外CH₄排放的近完全降低（98.21%）。其抑制机制涉及：（1）重塑瘤胃细菌群落以利于促丙酸发酵模式（如Prevotella增加、Ruminococcus减少），降低A:P比值，使H₂远离产甲烷途径；（2）直接抑制产甲烷作用，减少Methanobrevibacter丰度及氢营养型和甲基营养型途径关键酶（如EC:2.1.1.86、EC:2.8.4.1和EC:2.1.1.246）水平，最终限制MCR活性并阻碍CH₄产生的终末步骤。这些发现证实3-NOP可与A. taxiformis协同作用，在大幅降低A. taxiformis有效用量的同时维持最大功效。该策略为反刍动物CH₄减排提供了更为可行和经济的方法，但仍需在长期体内研究中进一步验证其对动物健康、生产性能和产品安全性的影响。