沙棘黄酮对泌乳中期荷斯坦奶牛瘤胃发酵参数、细菌菌群结构及代谢产物的影响

赵慧颖 余诗强 贺李莹 熊本海 王 慧* 蒋林树*

(1.北京农学院动物科学技术学院，奶牛营养学北京市重点实验室，北京102206；2.中国农业科学院北京畜牧兽医研究所，北京100193)

沙棘是胡颓子科、沙棘属落叶性灌木，又称醋柳果、酸刺等，具有耐旱、抗风沙和在盐碱土地上生存的特性，因此在中国西北部常被用于沙漠绿化。同时沙棘也是一种药食同源的植物，含有黄酮类、皂苷类、维生素、糖类和有机酸等多种功能性成分，具有抗氧化、抗应激、促生长、抗肿瘤、降血压和降血脂等生理作用，其发挥作用的主要成分为黄酮类化合物[1-2]。沙棘黄酮还具有很强的体外抑菌作用，对大肠杆菌、枯草芽孢杆菌和金黄色葡萄球菌等常见的病原菌都有抑制作用[3]。沙棘黄酮在人类生产生活中应用广泛，并且其在绵羊、鸡和猪上的应用已有报道[4-5]。有研究表明，沙棘黄酮可提高肉鸡平均日增重，降低腹脂率[6]；提高断奶仔猪血清中瘦素含量，有助于改善猪的胴体品质[7]。在反刍动物方面，刘艳丰[8]在阿勒泰羊的饲粮中添加沙棘叶黄酮，发现其能显著提高瘤胃挥发性脂肪酸(VFA)的浓度，降低氨态氮(NH3-N)浓度，促进瘤胃代谢和微生物蛋白的合成，并提高其免疫性能。白齐昌等[9]通过体外产气法研究表明，沙棘黄酮可以改善绵羊的体外发酵，提高干物质降解率、有机物降解率和代谢能，添加水平过高会对纤维降解菌具有抑制作用，适宜添加水平为0.3%。目前，沙棘黄酮在反刍动物上多以羊为研究对象，并且大多进行体外试验，无法真正反映出沙棘黄酮在动物体内的作用效果。同时，在以往的研究中可知，黄酮对奶牛瘤胃发酵的改善及瘤胃微生物菌群的影响具有显著效果[10-12]，而沙棘黄酮在奶牛上的应用却鲜有文献报道。本实验室前期已探究了沙棘黄酮对奶牛生长性能、乳中生物活性成分、血清免疫指标的影响[13]，发现沙棘黄酮能够提高奶牛产奶量和乳成分含量以及抗氧化能力。因此，本研究在前期研究的基础上，进一步探讨沙棘黄酮对泌乳中期荷斯坦奶牛瘤胃发酵参数、细菌菌群结构及代谢产物的影响，以期为今后沙棘黄酮在反刍动物上的相关研究提供理论基础。

1 材料与方法

1.1 试验材料

本试验所用沙棘黄酮购自西安某科技有限公司，纯度为40%，其他成分为沙棘多糖(20%)、灰分(20%)、粗蛋白质(10%)、水分(3%)、其他(7%)。

1.2 试验设计与饲养管理

本试验在北京奶牛中心延庆基地良种场进行，选择产奶量[(35.5±3.6) kg/d]、泌乳日龄[(91.0±22.1) d]、胎次(1.3±0.5)相近的健康中国荷斯坦奶牛30头，随机分为3组，每组10头，在饲喂全混合日粮(TMR)的基础上，分别以口腔投喂方式饲喂0、60、100 g/(d·头)的沙棘黄酮，依次标记为BC组、BL组和BM组(沙棘黄酮的添加量根据本实验室前期对竹叶黄酮的研究经验及查阅沙棘黄酮的相关文献[9,14-15]确定)。每组的试验牛均单独采用散栏带卧床式饲养，每天07:00、13:00、18:00对奶牛进行投料，09:30、15:30、21:30进行挤奶。试验期间奶牛自由采食，自由饮水，牛舍每天清扫、清粪。试验期共计75 d，其中预试期15 d，正试期60 d。基础饲粮参照NRC(2001)奶牛饲养标准配制，以TMR形式饲喂，TMR组成及营养水平见表1。

表1 TMR组成及营养水平(干物质基础)

1.3 瘤胃液的采集

在正试期最后1天于晨饲前1 h使用口腔采液器采集瘤胃液。每头奶牛采集200 mL瘤胃液，采集的前50 mL瘤胃液舍弃(避免奶牛口腔黏液和采液器中的残留液体产生影响)，采集后立即用4层纱布过滤，测定pH，然后将滤液分装到2 mL冻存管和10 mL离心管(按照瘤胃液体积与25%偏磷酸体积为10∶1的比例加入25%偏磷酸)。2 mL冻存管立即置于液氮中，于-80 ℃保存，用于测定瘤胃细菌菌群结构及代谢物含量；10 mL离心管置于-20 ℃冰箱中保存，用于测定NH3-N及VFA浓度。

1.4 瘤胃发酵参数的测定

pH使用PHS-10型便携式pH计测量。NH3-N浓度参照冯宗慈等[17]的比色法测定。使用安捷伦7890B气相色谱仪检测瘤胃中VFA的浓度，气相色谱条件如下：硝基对苯二甲酸改性的聚乙二醇(DB-FFAP)填充柱，进样量为2 μL，氮气(N2)流速为30 mL/min，氢气(H2)流速为40 mL/min，空气流速为400 mL/min，总流量为54 mL/min，吹扫流量为3 mL/min，平均线速度为27.8 cm/s，分流比为50∶1，柱温为180 ℃，检测温度为280 ℃，进样口温度为250 ℃。

1.5 瘤胃细菌多样性的测定与分析

本实验室前期研究结果表明，在饲粮中添加60 g/d的沙棘黄酮对生产性能的提高效果更佳[13]。因此，本研究只探究了0与60 g/d的沙棘黄酮对奶牛瘤胃细菌菌群及代谢物的影响。

瘤胃液在冰上解冻后使用PowerSoil DNA Isolation Kit试剂盒进行DNA提取，使用PCR进行扩增(引物：V3～V4区域，338F：5’-ACTCCTACGGGAGGCAGCAG-3’，806R：5’-GGACTACHVGGGTWTCTAAT-3’)，将特定的接头序列拼接在目的DNA两端，构建瘤胃微生物16S rRNA基因文库。基于Illumina HiSeq双端测序得到原始数据，经碱基识别后转为sequencedreads，结果以FASTQ(fq)格式保存。

将HiSeq测序得到的数据拼接(merge)成1条序列Tags，并对reads质量和merge效果进行质控过滤。具体步骤参照余诗强等[11]的方法进行分析。

1.6 瘤胃代谢产物的测定与分析

1.6.1 代谢物萃取及衍生化

瘤胃液室温解冻后涡旋10 s混匀，取50 μL瘤胃液加入300 μL分析甲醇，涡旋3 min，12 000 r/min、4 ℃下离心10 min；取上清液200 μL，20 ℃下静置30 min，12 000 r/min、4 ℃下离心3 min；取150 μL上清液用于上机分析。

1.6.2 代谢组学数据分析

数据采集仪器系统主要包括超高效液相色谱(UPLC)和四级杆飞行时间质谱(AB SCIEX Triple TOF6600)。代谢物的定性定量分析使用三重四级杆质谱的多反应监测模式(MRM)完成，得到不同样品的液相色谱数据，对所有代谢物的提取离子色谱峰的峰面积进行积分，并校正不同样品中同一代谢物的色谱峰。

1.7 数据统计与分析

采用Excel 2021进行数据整理，用SPSS 25.0中one-way ANOVA程序进行单因素方差分析，使用Duncan氏多重比较对组间差异进行分析，P≤0.05为差异显著。16S rRNA高通量测序结果采用SPSS 25.0统计软件中的独立样品t检验进行分析，P<0.05为差异显著。

2 结果与分析

2.1 沙棘黄酮对泌乳中期荷斯坦奶牛瘤胃发酵参数的影响

由表2可知，饲粮中添加沙棘黄酮对泌乳中期荷斯坦奶牛瘤胃pH及NH3-N、丙酸、异丁酸、丁酸、异戊酸和总挥发性脂肪酸浓度均无显著影响(P>0.05)。与BC组相比，BL组的乙酸浓度及乙酸/丙酸显著升高(P<0.05)，但BL组及BM组的戊酸浓度显著降低(P<0.05)。

2.2 沙棘黄酮对泌乳中期荷斯坦奶牛瘤胃细菌菌群的影响

2.2.1 稀释曲线及分类操作单元(OTU)的Venn图分析

由图1可以看出，随着测量的样本量的增加，稀释曲线趋于平缓。由图2可知，BC组和BL组泌乳中期荷斯坦奶牛瘤胃液中共含有8 399个OTU，2组共有5 933个OTU，BC组特有1 155个OTU，占总OTU的13.75%，BL组特有1 311个OTU，占总OTU的22.10%，2组之间OTU数目存在一定的差异，且BL组含有的OTU数目较多。

2.2.2 瘤胃细菌菌群的α多样性指数分析

由表3可知，BL组的Chao1指数和Shannon指数与BC组相比均无显著差异(P>0.05)，与BC组相比，BL组的Simpson指数有下降的趋势(P=0.06)，Coverage指数有上升的趋势(P=0.09)。

2.2.3 瘤胃菌群结构分析

2.2.3.1 门水平分析

由图3可知，BC组、BL组瘤胃液中共检测出24个细菌菌门，其中拟杆菌门(Bacteroidota)、厚壁菌门(Firmicutes)、变形菌门(Proteobacteria)、螺旋体门(Spirochaetota)及Patescibacteria相对丰度>1%；拟杆菌门、厚壁菌门的相对丰度占比前2，为门水平上的优势菌门。与BC组相比，BL组的厚壁菌门、Patescibacteria相对丰度较高，其他菌门相对丰度较低，但差异均不显著(P>0.05)。

表2 沙棘黄酮对泌乳中期荷斯坦奶牛瘤胃发酵参数的影响

横坐标代表随机抽取的测序数据量；纵坐标代表观测到的物种数量。

2.2.3.2 属水平分析

图4可知，BC组、BL组瘤胃液中共检测出334个细菌菌属，其中18个菌属的相对丰度>1%；普雷沃氏菌属(Prevotella)、解琥珀酸菌属(Succiniclasticum)和Muribaculaceae_norank的相对丰度占比前3，为属水平上的优势菌属。与BC组相比，BL组泌乳中期荷斯坦奶牛瘤胃液中拟杆目菌RF16菌群(Bacteroidales_RF16_group)、理研菌科_RC9_肠道群(Rikenellaceae_RC9_gut_group)的相对丰度显著增加(P<0.05)，F082_norank的相对丰度有增加的趋势(0.050.05)。

BC代表BC组；BL代表BL组。下图同。

2.2.3.3 线性判别分析(LEFSE)

由图5可知，丰佑菌目(Opitutales)、希瓦氏菌科(Shewanellaceae)、紫红球菌科(Puniceicoccaceae)、普雷沃氏菌属_9(Prevotella_9)及菖苣属(Psychrobium)这5种细菌在BC组中发挥重要作用；无胆甾原体目(Acholeplasmatales)、Ethanoligenenaceae、瘤胃梭菌属(Ruminiclostridium)、梭状芽孢杆菌属_vadinBB60_菌群(Clostridia_vadinBB60_group)、鞘氨醇杆菌科(Sphingobacteriaceae)、厌氧原体属(Anaeroplasma)、理研菌科_RC9_肠道菌群(Rikenellaceae_RC9_gut_group)、理研菌科(Rikenellaceae)、胆质菌科(Acholeplasmataceae)、拟杆菌目_RF16_菌群(Bacteroidales_RF16_group)、拟杆菌门_BD2_2(Bacteroidetes_BD2_2)、鞘氨杆菌属(Sphingobacterium)、毛螺菌科_FCS020_菌群(Lachnospiraceae_FCS020_group)、鞘氨醇杆菌目(Sphingobacteriales)、噬糖厌氧产氢杆菌属(Hydrogenoanaerobacterium)、毛螺菌科_FE2018_菌群(Lachnospiraceae_FE2018_group)、UCG-004、Allorhizobium_Neorhizobium_Pararhizobium_Rhizobium这18个细菌在BL组中发挥重要作用。

表3 沙棘黄酮对泌乳中期荷斯坦奶牛瘤胃细菌多样性指数的影响

Bacteroidota:拟杆菌门;Cyanobacteria:蓝菌门;Desulfobacterota:脱硫杆菌门; Fibrobacterota:纤维杆菌门; Firmicutes:厚壁菌门; Others:其他; Proteobacteria:变形菌门; Spirochaetota:螺旋菌门; Unclassified:未分类序列。

2.3 沙棘黄酮对泌乳中期荷斯坦奶牛瘤胃代谢物的影响

2.3.1 正交偏最小二乘法判别分析(OPLS-DA)

从图6可以看出，BC组、BL组之间的分离是较好的，代表组间聚类差异显著。

2.3.2 差异代谢产物分析

由表4可知，2组共筛选出25种差异代谢物。与BC组相比，BL组6种差异代谢物含量显著下调(P<0.05)，包括(2S)-2-氨基-5-戊酮酸[(2S)-2-amino-5-oxopentanoic acid]、3-脱氢奎尼酸(3-dehydroquinic acid)、D-来苏糖(D-lyxose)、磷酸高丝氨酸(o-phosphohomoserine)、5,6-二羟基-1,3-二烯-1-羧酸(5,6-dihydroxycyclohexa-1,3-diene-1-carboxylate)以及N-acetylpuromycin；19种差异代谢物含量显著上调(P<0.05)，包括苯丙酮酸(phenylpyruvic acid)、邻氨基苯甲酸(anthranilic acid)、岩白菜素(bergenin)、烟胺比林(nifenazone)、天竺葵色素(pelargonidin)、芦丁(rutin)、7,8,4′-三羟基异黄酮(7,8,4′-trihydroxyisoflavone)、西阿尼醇(cianidanol)、顺-3,4-白天竺葵甙元(cis-3,4-leucopelargonidin)、2-羟基-2,3-染料木素(2-hydroxy-2,3-dihydrogenistein)、(2S)-2-azaniumyl-3-[4-(3-methylbut-2-enyl)-1H-indol-3-yl]propanoate2-phenyl-1H-benzimidazole-5-sulfonic acid、大黄素(emodin)、2,3,7,9-tetrahydroxy-1-methyl-6H-benzo[c]chromen-6-one、2,6,7,4′-tetrahydroxyisoflavanone、2-[2-(1H-indol-3-yl)acetyl]amino]propanoate、香橙素(aromadendrin)、4-[(E)-2-carboxyethenyl]-2,6-dimethoxyphenolate及黄颜木素(fustin)。

[Eubacterium]coprostanoligenes group_norank:产粪甾醇真杆菌属未分类群; [Eubacterium]ruminantium group:啮齿真杆菌属群; Anaeroplasma:厌氧枝原体属; Anaerovibrio:厌氧弧菌属; Bacteroidales RF16 group_norank:拟杆菌目RF16未分类菌群; Bacteroidales_norank:未分类拟杆菌目; Butyrivibrio:丁酸弧菌属; Christensenellaceae R-7 group:克里斯滕森菌科R-7菌群; Clostridia UCG-014_norank:梭状芽孢杆菌属UCG-014未分类;Fibrobacter:丝状杆菌属; Lachnospiraceae_NK3A20_group:毛螺菌科_NK3A20_菌群; Lachnospiraceae_Unclassified:毛螺菌科_未分类序列; Others:其他; Prevotella:普雷沃氏菌属; Prevotellaceae Ga6A1 group:普雷沃氏菌科Ga6A1菌群; Prevotellaceae UCG-001:普雷沃氏菌科UCG-001; Prevotellaceae UCG-003:普雷沃氏菌科UCG-003; Prevotellaceae UCG-004:普雷沃氏菌科UCG-004; Pseudobutyrivibrio:假丁酸弧菌属; Rikenellaceae_RC9_gut_group:理研菌科_RC9_菌群; Ruminococcus:瘤胃球菌属; Succiniclasticum:解琥珀酸菌属; Succinivibrionaceae UCG-001:琥珀酸弧菌科UCG-001; Succinivibrionaceae UCG-002:琥珀酸弧菌科UCG-002; Treponema:密螺旋体属; Unclassified:未分类序列; Veillonellaceae UCG-001:韦荣球菌科UCG-001。

2.3.3 差异代谢物KEGG功能注释及富集分析

从图7、图8可知，25种差异代谢物主要分布在与蛋白质代谢、糖类代谢相关的9条代谢途径上，其中代谢产物显著富集在苯丙氨酸、酪氨酸和色氨酸生物合成途径以及氨基酸生物合成途径中。相比于其他代谢途径，富集在上述2个代谢途径上的差异代谢物相对较多，分别为3和5个。

3 讨 论

3.1 沙棘黄酮对泌乳中期荷斯坦奶牛瘤胃发酵参数的影响

VFA包括乙酸、丙酸、丁酸、戊酸和异戊酸等，是淀粉等碳水化合物在瘤胃内发酵的终产物，也是反刍动物维持正常生命活动及生产所需的主要能量来源[18]。乙酸是动物合成脂肪的前体，也是合成乳脂的主要成分[19-20]。丙酸主要被转移到肝脏通过糖异生作用生成葡萄糖，为葡萄糖合成提供80%～90%的原料[21]。戊酸与异戊酸、异丁酸等含4～5个碳原子的短链VFA称为异位酸，异位酸可以提高瘤胃内微生物蛋白的合成，抑制脲酶活性，增加非蛋白氮的利用，同时，还可以促进瘤胃内乙酸的合成[22-24]。白齐昌等[9]通过体外产气法研究沙棘黄酮对绵羊体外产气、瘤胃发酵参数和微生物菌群的影响，结果发现添加0.3%的沙棘黄酮显著增加了乙酸和丙酸的浓度。此结果与本研究相似，添加60 g/d的沙棘黄酮使泌乳奶牛瘤胃内乙酸浓度及乙酸/丙酸显著升高，戊酸浓度显著下降。这说明添加60 g/d的沙棘黄酮促进了瘤胃微生物中有益菌的生长，促进乙酸生成菌的繁殖与代谢，有助于其对碳水化合物的发酵与分解，进而促进泌乳奶牛对营养物质的代谢与吸收；并且沙棘黄酮并没有改变奶牛的发酵类型，奶牛仍是乙酸型发酵，甚至使乙酸生成量升高，在后续的研究中可知沙棘黄酮可促进奶牛瘤胃内能生成乙酸的纤维素分解菌生长，并促进氨基酸代谢。而对于某些纤维素分解菌来讲氨是唯一合成自身蛋白质的氮源[25]。这说明沙棘黄酮有利于瘤胃内产生乙酸的纤维素分解菌的生长，进而使乙酸生成增加。同时作为异位酸的戊酸含量下降、乙酸含量升高，也可进一步推测沙棘黄酮可促进异位酸向乙酸转化，乙酸可进一步转化为乙酰辅酶A作为脂肪酸合成的前体参与合成脂肪酸，乙酸的增加可为乳脂的合成提供前体物质。

Opitutales:丰佑菌目; Shewanellaceae:希瓦氏菌科; Puniceicoccaceae:紫红球菌科; Prevotella_9:普雷沃氏菌属_9; Psychrobium:菖苣属; Acholeplasmatales:无胆甾原体目; Ruminiclostridium:瘤胃梭菌属; Clostridia_vadinBB60_group:梭状芽孢杆菌_vadinBB60_菌群; Sphingobacteriaceae:鞘氨醇杆菌科; Anaeroplasma:厌氧原体属; Rikenellaceae_RC9_gut_group:理研菌科_RC9_肠道菌群; Rikenellaceae:理研菌科; Bacteroidales_RF16_group:拟杆菌目_RF16_菌群; Bacteroidetes_BD2_2:拟杆菌门_BD2_2; Sphingobacterium:鞘氨杆菌属; Lachnospiraceae_FCS020_group:毛螺菌科_FCS020_菌群; Sphingobacteriales:鞘氨醇杆菌目; Hydrogenoanaerobacterium:噬糖厌氧产氢杆菌属; Lachnospiraceae_FE2018_group:毛螺菌科_FE2018_菌群; Acholeplasmataceae:胆质菌科。

3.2 沙棘黄酮对泌乳中期荷斯坦奶牛瘤胃细菌菌群的影响

反刍动物瘤胃内的微生物复杂多样，包括细菌、古菌、真菌、原虫以及少量噬菌体，这些微生物可以帮助宿主消化和同化食物，为宿主提供能量和营养，与动物机体存在着共生关系[26-27]。本研究以泌乳中期荷斯坦奶牛为研究对象，利用16S rRNA高通量测序技术探讨饲粮中添加沙棘黄酮对奶牛瘤胃细菌多样性的影响。本研究中样本覆盖率高达98%，且稀释曲线趋于平缓，说明样本量充足能够真实地反映泌乳中期荷斯坦奶牛瘤胃细菌的种类和多样性。α多样性指数分析表明，BL组Simpson指数有降低的趋势，其他指数均无显著变化，说明饲粮添加60 g/d的沙棘黄酮可使瘤胃内细菌多样性增加，而对其丰富度无显著影响。研究表明，黄酮类化合物具有抑菌活性，能够抑制一些对动物有害的细菌、真菌等[10,28]。瘤胃细菌多样性增加可能是因为沙棘黄酮可以抑制一些病原菌生长，从而使一些与其竞争、拮抗的细菌数量增加。

前人研究表明，瘤胃细菌主要由拟杆菌门和厚壁菌门2个门组成，它们共同构成了瘤胃内容物中约90%的16S rRNA基因丰度，且不会随品种、饲粮和饲养环境的改变而改变[29]。本研究中泌乳中期荷斯坦奶牛的优势菌门也为拟杆菌门和厚壁菌门。拟杆菌主要通过分泌碳水化合物反应酶来分解饲粮中的碳水化合物，但对纤维素不起作用[30]。而厚壁菌门主要起到降解纤维素的作用[31]。本研究中，与BC组相比，BL组厚壁菌门的相对丰度较高。Zhan等[32]的研究表明，在初产荷斯坦母牛的饲粮中添加苜蓿黄酮提取物可以改变瘤胃微生物组成，使拟杆菌门相对丰度降低，厚壁菌门相对丰度升高，但差异不显著，此结果与本研究相似。De Nardi等[33]饲喂100 g/d的多酚精油(主要成分为黄酮)给荷斯坦奶牛，增加了厚壁菌门和拟杆菌门的丰富度，说明黄酮可在一定程度上促进厚壁菌的生长与增殖，进一步促进瘤胃内纤维素的降解，提升荷斯坦奶牛对营养物质的消化与吸收功能。本试验中添加沙棘黄酮后，瘤胃中厚壁菌门和拟杆菌门等优势菌门的丰富度并没有显著变化，说明沙棘黄酮对瘤胃菌群无不良影响。

横坐标表示预测成分得分值，横坐标方向可以看出组间的差距；纵坐标表示正交成分得分值，纵坐标方向可以看出组内的差距；百分比表示成分对数据集的解释度。

从属水平上分析，本次研究2组的优势菌属为普雷沃氏菌属，这与前人研究结果[34-36]一致。普雷沃氏菌属主要代谢瘤胃中的半纤维素、果胶、蛋白质等，主要产生乙酸、甲酸等[37]。本研究中，BL组较BC组的普雷沃氏菌属、解琥珀酸菌属的相对丰度有所下降，可能是由于厚壁菌门的相对丰度升高，厚壁菌门中纤维降解菌的相对丰度升高抑制了普雷沃氏菌的生长。有研究表明，解琥珀酸菌属中的产琥珀酸杆菌相对丰度与乙酸的生成呈负相关与戊酸的生成呈正相关[38]。与本试验的研究结果一致，产琥珀酸菌相对丰度的下降使乙酸的生成增加，戊酸的生成减少。毛螺菌科_FE2018_菌群、毛螺菌科_FCS020_菌群和Tuzzerella 3个菌属同属于毛螺菌科厚壁菌门，均为BL组中发挥重要作用的菌群，能够降解多糖、纤维生成乙酸[39]。并且与BC组相比，BL组的泌乳中期荷斯坦奶牛瘤胃细菌中氢氧杆菌属、瘤胃梭菌属及理研菌科_RC9_肠道菌群的相对丰度显著增加。研究表明，理研菌科主要在降解纤维多糖方面发挥作用，瘤胃梭菌属也具有降解纤维素的功能，二者均能生成乙酸，而氢氧杆菌属也具有降解葡萄糖生成乙酸的功能[40-43]。占今舜[10]在奶牛饲粮中添加苜蓿黄酮，发现普雷沃氏菌属、琥珀酸菌属的相对丰度随着苜蓿黄酮添加量的升高呈线性下降趋势，该结果与本研究相似。这说明沙棘黄酮可以提高瘤胃内细菌的丰富度和多样性，并且通过提高瘤胃内纤维降解菌的相对丰度，提高乙酸的浓度，进而提高乳脂率。

表4 差异代谢物

4ItemsVIP2) Fold change Type Nifenazone3.1311 067.32up2,3,7,9-tetrahydroxy-1-methyl-6H-benzo[c]chromen-6-one3.1349 084.33up7,8,4′-trihydroxyisoflavone1.902.36up Aromadendrin2.9716.99up Cianidanol2.7215.94up-3,4- cis-3,4-leucopelargonidin2.8012.63up Fustin1.0520.85up Pelargonidin2.9816.87up Rutin1.016.94up2,6,7,4′-tetrahydroxyisoflavanone2.9116.30up2--2,3- 2-hydroxy-2,3-dihydrogenistein1.0622.09upN-acetylpuromycin2.100.43down5,6--1,3--1-5,6-dihydroxycyclohexa-1,3-diene-1-carboxylate1.030.32down(2S)-2-azaniumyl-3-[4-(3-methylbut-2-enyl)-1H-indol-3-yl]propanoate3.0171.69up2-[2-(1H-indol-3-yl)acetyl]amino]propanoate2.8727.77up4-[(E)-2-carboxyethenyl]-2,6-dimethoxyphenolate3.139 216.32up

Tryptophan metabolism:色氨酸代谢; Phenylalanine, tyrosine and tryptophan biosynthesis:苯丙氨酸、酪氨酸和色氨酸的生物合成; Phenylalanine metabolism:苯丙氨酸代谢; Pentose and glucuronate interconversions:戊糖和葡糖醛酸的相互转化; Metabolic pathways:代谢途径; Glycine serine and threonine metabolism:甘氨酸、丝氨酸和苏氨酸代谢; Biosynthesis of amino acids:氨基酸的生物合成; Arginine and proline metabolism:精氨酸和脯氨酸代谢; 2-oxocarboxylic acid metabolism:2-氧代羧酸代谢。下图同 the same as below。

横坐标表示每个通路对应的Rich factor，纵坐标为通路名称，点的颜色为P值，越红表示富集越显著。点的大小代表富集到的差异代谢物的个数多少。

3.3 沙棘黄酮对泌乳中期荷斯坦奶牛瘤胃代谢物的影响

相对于一些传统的研究方法，代谢组学可以通过核磁共振、液相色谱/质谱(LC/MS)等高通量技术更加清晰直观地检测到机体代谢物和代谢通路的变化[44]。与BC组相比，BL组的泌乳中期荷斯坦奶牛瘤胃中6种差异代谢物含量显著下调，19种差异代谢物含量显著上调。对上述25种差异代谢物进行KEGG富集分析，结果表明，25种代谢物分别分布在9条代谢通路上，代谢物显著富集的通路有2条，分别是苯丙氨酸、酪氨酸和色氨酸生物合成通路、氨基酸的生物合成通路。

与氨基酸代谢相关的通路中共检测到5种代谢产物，其中2种代谢物在BL组中峰值更高，分别是苯丙酮酸以及邻氨基苯甲酸；剩余3种代谢物在BC组的峰值更高，分别是谷氨酰胺、3-脱氢奎尼酸以及高丝氨酸邻磷酸。芳香族氨基酸是高等动物所必需的，通常在细菌和植物中通过莽草酸途径由葡萄糖合成[45]。芳香族氨基酸包括苯丙氨酸、酪氨酸和色氨酸，在蛋白质生产和糖原合成等过程中都发挥着重要的作用[46]。苯丙氨酸为生糖生酮氨基酸，主要通过不可逆的羟基化反应生成酪氨酸，还可通过转氨基作用先生成对羟基苯乙酸之后再生成苯丙酮酸，苯丙酮酸可以转化为富马酸盐进入三羧酸循环供能，若生成过多则会损伤中枢神经系统或出现苯丙酮尿症[47]。

3-脱氢奎尼酸、谷氨酰胺和邻氨基苯甲酸均为莽草酸代谢途径的中间产物，在微生物合成芳香族氨基酸过程中发挥着重要作用[48]。首先，磷酸烯醇式丙酮酸(PEP)和赤藓糖-4-磷酸发生缩合作用生成3-脱氢奎尼酸，之后转化为莽草酸，再与PEP结合生成分支酸，一部分分支酸在羧甲基变位酶(CM)的催化下转化为预苯酸，另一部分则经脱羧、脱水、脱羟基等过程生成苯丙氨酸、酪氨酸等氨基酸。预苯酸可被邻氨基苯甲酸合成酶(AS)催化生成谷氨酰胺后转化为邻氨基苯甲酸，之后进一步生成色氨酸[49]。本研究中观察到3-脱氢奎尼酸和谷氨酰胺的生成下调，苯丙酮酸和邻氨基苯甲酸的生成上调，说明沙棘黄酮对芳香族氨基酸的代谢过程具有调节作用。苯丙酮酸含量适当上调有助于能量的生成，而邻氨基苯甲酸含量上调有助于色氨酸的生成，进而有助于动物性蛋白质的补充，有利于奶牛对蛋白质的吸收利用。

4 结 论

本研究表明，在泌乳中期荷斯坦奶牛饲粮中添加60 g/(d·头)沙棘黄酮提高了乙酸浓度和乙酸/丙酸，降低了戊酸浓度，提高了奶牛瘤胃中拟杆菌目RF16菌群、理研菌科_RC9_肠道菌群、瘤胃梭菌属的相对丰度，同时提高了泌乳中期荷斯坦奶牛瘤胃中与苯丙氨酸、酪氨酸和色氨酸生物合成通路相关的代谢物含量，为脂质的合成提供了更多的前体物。