科研 | NPJ BIOFILMS MICROBI：牦牛适应高海拔生活过程中，食物-肠道微生物相互作用的季节变化（国人佳作）( 三 )

R=0.47 ， p<0.0001；adonis分析：R2=0.16 ， p<0.0001），图基于Bray-Curtis不相似性绘制在非计量多维标度（NMDS）上。相似性分析（ANOSIM）、adonis分析和PERMANOVA分析用于相似性处理的统计检验。虚线椭圆边框表示95%的置信区间。

本文插图

图2 季节性迁移放牧和自由放牧条件下牦牛食物组成的季节变化。 Stream图显示了春季、夏季、秋季和冬季在季节性迁移放牧（a）和自由放牧（c）条件下植物科水平类群的相对丰度。低丰度类群（<5%）被归为“其他类群” 。在季节性迁移放牧（b）和自由放牧（d）区域，使用Sankey图跟踪与每个季节相关的指示科。线条表示指示植物科和季节之间的关联，它们由植物科着色。线宽按比例缩放以反映指标值（族的指标值越高，与季节的关联越强）。指示值如图S7所示。统计p值表示与季节相关的家系, *p<0.05,**p<0.01, ***p<0.001 。在OCG ，温暖季节（6~9月）地上生物量（AGB）最高， 9月份为2891±148 kg DM/ha ，在寒冷季节（10月~5月）最低， 5月份为99±10.8 kg DM/ha（图3a）。中性洗涤纤维（NDF）和酸性洗涤纤维（ADF）在寒冷季节分别比温暖季节高出10%和16% ，而粗蛋白（CP）和乙醚提取物（EE）分别比温暖季节低53%和25%（图3b）。具体来说，最高CP（166± 10.2 g/kg DM）出现在6月份，最低（39± 5.5 g/kg DM）出现在4月份。在TH区也观察到类似的趋势（图S3a, b）。因此，牦牛在寒冷季节不得不面临低氮的食物。

本文插图

图3 自由放牧草地全年地上生物量（AGB）和食物化学组成（基于干物质）的季节变化。折线图a表示AGB（kg DM/ha）。虚线是全年平均AGB 。折线图b表示全年食物的化学成分（粗蛋白（CP）、乙醚提取物（EE）、酸性洗涤纤维（NDF）和中性洗涤纤维（ADF））。虚线是全年各化学成分的平均值，并由各化学成分着色。数值和误差线显示为平均值±标准差。 CP平均含量与NDF（c）（R22=0.78 ， p<0.001）呈负相关。统计检验采用FDR（错误发现率）校正p值的t-test ， *p<0.05 ， **p<0.01 ， ***p<0.001 。3 不同季节肠道菌群的多样性及组成肠道微生物群落包括25375771个高质量序列、1125235个去除单子后的独特序列和14239个OTUs（表S3）。在TH和OCG ，肠道微生物群落的季节变化都很明显（图1c ， d）。共识别出18个肠道菌门，在各季节和各放牧方式下， Firmicutes和Bacteroidetes都最为丰富，（图S6、图S7c、d）。在TH和OCG ，肠道菌群属水平上表现出明显的季节动态，其中

科研 | NPJ BIOFILMS MICROBI： 牦牛适应高海拔生活过程中，食物-肠道微生物相互作用的季节变化（国人佳作）( 三 )

科研 | NPJ BIOFILMS MICROBI：牦牛适应高海拔生活过程中，食物-肠道微生物相互作用的季节变化（国人佳作）( 三 )