材料与方法
本研究于2023年11月在中国新疆喀什地区塔什库尔干塔吉克自治县（海拔约3800米）的自然牧场进行。研究对象为六头健康的雄性帕米尔牦牛，平均体重260±13.29公斤，年龄约3岁。屠宰后立即解剖，采集八个胃肠道区段的内容物：四个胃室（瘤胃Rum、网胃Ret、瓣胃Oma、皱胃Abo）和四个肠段（十二指肠Duo、空肠Jej、盲肠Cec、直肠Rec）。样品经过滤、分装后液氮速冻，用于后续分析。
挥发性脂肪酸（VFA）测定采用液相色谱-串联质谱法。样品经过预处理和衍生化后，使用超高效液相色谱-电喷雾电离-串联质谱进行分析，定量总VFA及乙酸、丙酸、丁酸、异丁酸、异戊酸、戊酸、己酸等单体浓度。
微生物群落分析通过16S rRNA基因高通量测序进行。使用通用引物对V3-V4可变区进行PCR扩增，在Illumina NovaSeq 6000平台上进行测序。原始数据经过质控、去噪、合并和去除嵌合体后，生成扩增子序列变体（ASV）并进行物种分类（基于SILVA数据库）。使用QIIME 2软件评估Alpha多样性（Shannon、Simpson、Ace、Chao1指数）和Beta多样性（基于Bray-Curtis距离的主坐标分析PCoA）。使用LEfSe分析识别不同区段的微生物标志物。使用PICRUSt2对16S rRNA基因序列进行功能预测，并通过STAMP软件比较不同组间预测功能通路的丰度差异。
结果
胃肠道挥发性脂肪酸分析
研究发现，帕米尔牦牛胃肠道不同区段的VFA浓度存在显著差异（P< 0.05）。总VFA浓度在瘤胃和网胃最高，皱胃和小肠最低，后肠（盲肠和直肠）呈现中等水平。其中，乙酸和丙酸是主要成分，其分布趋势与总VFA基本一致。丁酸和异丁酸在瘤胃中显著富集，其次是网胃。异戊酸和戊酸也在瘤胃和网胃最高，后肠次之。己酸则主要集中于瘤胃，在后肠含量显著降低。
测序结果与微生物多样性
经过质控，共获得2,582,129条有效序列。ASV韦恩图分析显示，随着食糜在消化道中前移，微生物群落经历了明显的空间分化和生态筛选。例如，瘤胃与网胃共享748个ASV，而空肠与盲肠仅共享42个ASV。Alpha多样性分析表明，网胃的Shannon和Simpson指数显著高于瘤胃（P< 0.05），皱胃的各项指数显著低于瓣胃（P< 0.01）。空肠的Alpha多样性指数在所有区段中最低（P< 0.01），而盲肠和直肠的多样性与前胃部分区段无显著差异。这证实了牦牛胃肠道微生物群沿着消化道纵轴呈现多样性梯度分布的趋势。
胃肠道微生物群落的分类与组成
主坐标分析（PCoA）显示，瘤胃-皱胃组与十二指肠-直肠组明显分离，瓣胃和皱胃的微生物组成相似。PC1和PC2分别解释了22.48%和11.37%的变异，表明不同胃肠道区域的微生物群落组成存在差异。
在门水平上，厚壁菌门（Firmicutes）、拟杆菌门（Bacteroidota）和变形菌门（Proteobacteria）是八个胃肠道区段的优势菌门。拟杆菌门在瘤胃、直肠和瓣胃中更丰富，而变形菌门在皱胃和十二指肠中含量最高。
在属水平上，共鉴定出347个属。前10位的优势属包括：普雷沃菌属（Prevotella）、Rikenellaceae_RC9_gut_group、UCG-005、毛螺杆菌科（Muribaculaceae）、克里斯滕森菌科_R-7群（Christensenellaceae_R-7_group）等。研究发现，普雷沃菌属（Prevotella）在瘤胃中丰度最高，在空肠和后肠中显著降低。毛螺杆菌科（Muribaculaceae）和F082在胃中的丰度显著高于肠道。UCG-005、UCG-010和拟杆菌属（Bacteroides）在盲肠和直肠中的丰度显著高于其他部分。而Romboutsia、狭义梭菌属1（Clostridium_sensu_stricto_1）和Turicibacter则仅在空肠中显著富集。
聚类热图显示，不同胃肠道部位的微生物群落呈现清晰的聚类模式。UCG-010、Monoglobus和Prevotellaceae_UCG-004主要富集于盲肠和直肠。狭义梭菌属1（Clostridium_sensu_stricto_1）、Romboutsia等则特异性地在空肠中高丰度存在。共现网络图分析揭示了优势菌属之间存在三个明显的模块化结构，表明微生物群落内部可能存在稳定的互作关系。
胃肠道各部位的生物标志物
LEfSe分析用于识别不同胃肠道区域的潜在微生物标志物。结果在瘤胃、网胃、瓣胃、皱胃、十二指肠、空肠、盲肠和直肠分别鉴定出2、1、5、2、8、4、3和4个属水平的生物标志物。例如，瘤胃的标志物是NK4A214_group和普雷沃菌属（Prevotella）；网胃是F082；瓣胃的优势标志物包括克里斯滕森菌科_R-7群（Christensenellaceae_R-7_group）、Rikenellaceae_RC9_gut_group等；皱胃是Ruminobacter和毛螺杆菌科（Muribaculaceae）；十二指肠富集了Capnocytophaga、瘤胃球菌属（Ruminococcus）等8个属；空肠以Romboutsia、Turicibacter等4个属为标志；盲肠富集UCG_010等3个属；直肠富集Monoglobus、拟杆菌属（Bacteroides）等4个属。
基于PICRUSt的功能基因预测分析
基于PICRUSt2的功能预测，比较了牦牛不同胃肠道区域微生物的功能通路。在瘤胃到网胃的转变中，戊糖和葡萄糖醛酸相互转化、各种植物次级代谢物的生物合成等21条通路在网胃显著富集。从网胃到瓣胃，鞭毛组装、细菌趋化性和硫代谢通路显著增强。从瓣胃到皱胃，与生物膜形成相关的通路以及阳离子抗菌肽（CAMP）抗性显著减少。从皱胃到十二指肠，核糖体、核苷酸代谢和嘧啶代谢等通路显著富集。从十二指肠到空肠，双组分系统、丙氨酸、天冬氨酸和谷氨酸代谢等通路显著减少。从空肠到盲肠，甘油磷脂代谢、硫代谢、D-氨基酸代谢等通路显著增加。从盲肠到直肠，磷酸转移酶系统（PTS）通路显著减少，而脂多糖生物合成、泛醌和其他萜醌生物合成通路显著增加。
挥发性脂肪酸与微生物的相关性分析
曼特尔检验评估了不同胃肠道区段微生物群落结构与VFA谱之间的关联。结果显示，瘤胃、瓣胃、皱胃和直肠的微生物群落与多种VFA呈显著正相关（P< 0.01）。相关性热图分析表明，大多数微生物属与VFA显著相关。其中，克里斯滕森菌科_R-7群（Christensenellaceae_R-7_group）、Rikenellaceae_RC9_gut_group等与VFA呈显著负相关（P< 0.05）。而Lautropia、Capnocytophaga等属与乙酸、丁酸、异戊酸等呈极显著负相关（P< 0.001）。此外，VFA之间也存在显著相关性，如丙酸与异戊酸、丙酸与异丁酸之间均存在强相关（r= 0.91）。
讨论
本研究系统分析了帕米尔牦牛胃肠道不同区域的微生物群落结构、VFA分布和功能通路特征，旨在阐明高原反刍动物消化道不同区域发酵代谢功能的空间异质性。VFA在牦牛胃肠道中呈现明显的空间梯度，前胃（瘤胃、网胃、瓣胃）浓度最高，皱胃和小肠显著下降，后肠出现“第二次峰值”。这种分布模式与优势菌属和关键代谢通路密切相关，表明肠道不同区域在能量获取和代谢调控中扮演不同角色。
前胃，尤其是瘤胃和网胃，作为主要的发酵场所，富含普雷沃菌属（Prevotella）、NK4A214_group等关键产酸菌，并富集了碳水化合物代谢相关通路，是宿主主要的能量来源。网胃中特有的F082可能与丙酸代谢密切相关。瓣胃的微生物群落和代谢活性显示出过渡性特征。
皱胃由于pH值较高、食糜停留时间短等特点，发酵活性有限，导致VFA浓度从前胃到此显著下降。小肠环境发生剧烈变化，微生物多样性降低，但十二指肠仍保持较高的代谢活性（如核苷酸代谢、蛋白质合成）。空肠的微生物群落最为简化，代谢通路也显著减少。
后肠（盲肠和直肠）接收前段未完全消化的物质进行“二次发酵”，其VFA浓度回升，并富集了UCG-005、拟杆菌属（Bacteroides）等特有菌属，以及丁酸代谢、免疫相关等通路，显示出潜在的代谢补偿功能。
相关性分析进一步证实了微生物群落与VFA代谢之间的紧密联系，不同区段的微生物可能通过协同或拮抗作用影响VFA的积累与转化。
结论
本研究揭示了帕米尔牦牛胃肠道微生物群落和VFA代谢存在显著的空间异质性。前胃是VFA生产的主要部位，富含产酸菌，是宿主的能量核心。小肠微生物多样性低但代谢功能活跃。后肠具有独特的细菌定植模式和潜在的代谢补偿特征。功能预测显示代谢通路从前胃到结肠存在分区化趋势。微生物与VFA的相关性分析表明，不同区段的微生物群落可能通过相互作用影响VFA代谢。总之，该研究揭示了帕米尔牦牛胃肠道微生物群落和代谢功能的空间分异，提示这种区域特化可能有助于其在高原条件下的能量获取和营养适应，但这一假说仍需进一步的功能验证。
