侯路钊 王增林 孙雪丽 刘桃桃 李秋凤 曹玉凤 李建国 高艳霞 陈攀亮

(河北农业大学动物科技学院，保定 071001)

反刍动物甲烷(CH4)排放与全球气候变暖和饲料能损失密切相关。甲烷能损失占饲料总能的2%～12%[1]。研究表明，反刍动物甲烷排放量占全球甲烷排放量的15%[2]。除此之外，畜牧业产生的甲烷是人为排放甲烷的重要来源之一，占全球人类活动甲烷排放量的28%[3]。因此，减少反刍动物的甲烷排放量对提高饲粮能量利用率和保护环境有重要作用。

反刍动物甲烷主要通过瘤胃产甲烷菌经二氧化碳(CO2)/氢气(H2)还原途径生成[4]。产甲烷菌属于瘤胃古菌，可利用多种降解底物生成甲烷。瘤胃产甲烷菌受多种因素的影响，其中饲粮组成对产甲烷菌影响较大[5-6]。研究表明，高纤维饲粮可促进产甲烷菌的增殖，并提高反刍动物甲烷日排放量[7]。与高纤维饲粮比，高淀粉饲粮可降低产甲烷菌多样性和甲烷菌代谢活性[8]。此外，研究表明饲粮非纤维性碳水化合物(NFC)/中性洗涤纤维(NDF)比值对甲烷排放量有显著影响[9]，随着饲粮NDF/NFC比值的降低，甲烷日排放量逐渐降低。这可能是低NDF/NFC比值的饲粮改变了瘤胃发酵模式，提高了丙酸利用菌与产甲烷菌竞争氢的能力，减少了甲烷菌合成甲烷的底物含量，并降低了产甲烷菌的活性，从而抑制了甲烷的生成[10]。目前对反刍动物甲烷减排的研究主要集中在调整饲粮组成；但是，通过调控饲粮NDF/淀粉比值来减少荷斯坦公牛甲烷排放量的相关研究较少。因此，本试验通过研究饲粮NDF/淀粉比值对荷斯坦公牛瘤胃古菌和甲烷排放量的影响，确定饲粮NDF/淀粉比值与甲烷排放量的相关关系，旨在为反刍动物甲烷减排提供理论依据和生产指导。

1 材料与方法

1.1 试验时间与地点

饲养试验于2018年2月初至2018年5月中旬在保定市宏达奶牛场完成。

1.2 试验设计

试验采用单因素完全随机试验设计，选择36头健康、膘情正常、体重[(282±37) kg]相近的荷斯坦公牛作为试验牛，随机分为3组，每组12头，散栏饲养。Ⅰ、Ⅱ、Ⅲ组试验牛分别饲喂NDF/淀粉比值为1.13、1.38、1.63的全混合日粮(TMR)。各组TMR能量水平相同，粗蛋白质水平接近。预试期7 d，正试期90 d。试验饲粮组成及营养水平见表1。

表1 试验饲粮组成及营养水平(干物质基础)Table 1 Composition and nutrient levels of experimental diets (DM basis) %

续表1项目Items组别 GroupsⅠⅡⅢ磷 P0.310.380.38淀粉 Starch 29.1926.9925.02中性洗涤纤维/淀粉 NDF/starch1.131.381.63

1.3 饲养管理

试验牛每天饲喂2次TMR(07:00和17：00)，自由采食与饮水。每天打扫料槽，每周清洗水槽2～3次，保持水槽和料槽干净、卫生。

1.4 瘤胃液的采集与处理

于正试期第88天晨饲前用口腔导管通过真空泵抽取瘤胃液，每组采集8头试验牛的瘤胃液，将采集的瘤胃液用4层纱布进行过滤后，装于已灭菌的离心管，立即投入液氮中，送回实验室后于-80 ℃冰箱保存。

1.5 绝对定量PCR

1.5.1 总DNA的提取

瘤胃液总DNA的提取主要参考郭嫣秋等[15]、Qiao等[16]、李长青等[17]的方法，具体理论依据参考李长青等[17]、李萌等[18]。将样品涡旋1 min，取200 μL瘤胃液样品，按Solarbio通用基因组DNA试剂盒(北京索莱宝科技有限公司)说明书提取瘤胃液总DNA。

1.5.2 总DNA的质量检测

使用Thermo NANo DROP 2000分光光度计检测瘤胃液总DNA的吸光度(OD)260 nm/280 nm在1.77～1.95，DNA浓度在25.0～97.7 ng/μL，并于-20 ℃保存。

1.5.3 引物序列

试验引物选用β-肌动蛋白(β-actin)[19]为内参基因，目的基因为反刍兽甲烷短杆菌(Methanobrevibacterruminantium)[15]，引物合成由生工生物工程(上海)股份有限公司完成，基因引物序列见表2。

表2 基因引物序列Table 2 Primer sequences of genes

1.5.4 普通PCR

打开Eppendorp PCR仪预热20 min，反应体系和条件见表3。

表3 反应体系和条件Table 3 Reaction system and conditions

1.5.5 琼脂凝胶电泳

配制1%琼脂糖凝胶[称取0.2 g琼脂糖加20 mL 0.5×硼酸(TBE)]放入带盖的瓶子中，置于微波炉中(调到中低档，1～2 min)，待瓶身微热加2 μL无毒核酸染料Super GelRed(产品货号：S2001，Everbright Inc，美国)，倒胶，20～30 min后拔梳子，倒0.5×TBE，没过胶面，加样2 μL，电泳条件设为：120 V，100 mA。琼脂凝胶电泳图谱结果见图1。

M：标记物 marker；M.r：为反刍兽甲烷短杆菌 Methanobrevibacter ruminantium；β：β-肌动蛋白 β-actin。图2同 the same as Fig.2。图1 琼脂凝胶电泳图谱Fig.1 Agar gel electrophoretogram

1.5.6 胶回收

称取0.25 g琼脂糖加25 mL 0.5×TBE(1%琼脂糖凝胶)放入带盖的瓶子中，置于微波炉中(调中低档，1～2 min)，待瓶身微热加2.5 μL无毒核酸染料Super GelRed，倒胶，20～30 min后拔梳子，倒0.5×TBE，没过胶面，加样20 μL，电泳条件设为：120 V，100 mA，切胶。使用SanPrep柱式DNA胶回收试剂盒[生工生物工程(上海)股份有限公司]回收做标准品，测其浓度在10～30 ng/μL，于-20 ℃保存。胶回收产物电泳图谱结果见图2。

图2 胶回收产物电泳图谱Fig.2 Gel recovery product electrophoretogram

1.5.7 标准曲线

10倍稀释标准品(5 μL标准品加45 μL Easy dilution)后做8个梯度(目的基因：1.13×1013～1.13×107；内参基因：5.74×1014～1.13×108)，每次加Easy dilution需涡旋3 min，离心30 s(保证混匀)，全部溶解后，做好标记，于-20 ℃保存。使用Bio-Rad荧光定量PCR仪检测样品，反应体系及条件见表4，选用6个点做标准曲线。熔解曲线有明显高峰，出峰温度在81 ℃，说明扩增产物特异性高。目的基因标准曲线：y=-3.543x+50.063(y是Ct值，x是log拷贝数)，扩增效率为91.5%，r2=99.6%；内参基因标准曲线：y=-3.029x+49.217(y是Ct值，x是log拷贝数)，扩增效率为113.0%，r2=98.4%。扩增曲线见图3，熔解曲线见图4，目的基因标准曲线见图5，内参基因标准曲线见图6。

图3 扩增曲线Fig.3 Amplification curve

图4 熔解曲线Fig.4 Melting curve

图5 目的基因标准曲线Fig.5 Standard curve line of target gene

图6 内参基因标准曲线Fig.6 Standard curve line of reference gene

1.5.8 样品测定

同一批次96孔板目的基因、内参基因各设标准曲线7个点、24个待测混池样品、空白样品各2个重复，待测样品目的基因、内参基因各3个重复，反应体系及条件同表4。

表4 反应体系及条件Table 4 Reaction system and conditions

1.6 指标测定

1.6.1 瘤胃古菌多样性测定

对古菌16S V9区进行Illumina HiSeq测序，由北京诺禾致源科技有限公司代测。数据库测定方法：参照Qiime V1.9.1软件的Tags质量控制流程对下机样本数据进行拼接、过滤后为有效数据(clean data)。利用Uparse v7.0.1001软件[20]以97%的一致性对有效数据进行操作分类单元(OTUs)聚类，对具有OTUs代表序列使用Mothur方法与SILVA 132数据库[21]进行物种注释。

1.6.2 样品六氟化硫(SF6)和甲烷浓度的测定以及甲烷排放量的计算

甲烷浓度测定：每次测定样品前应先用2.74 mg/m3的甲烷标准气体进行校正。色谱条件：载气为99.999%高纯氮气，燃气为99.999%高纯H2，空气为助燃气，柱温为70 ℃，前进样口温度为50 ℃，后进样口温度为375 ℃，火焰离子化检测器(FID)温度为200 ℃，进样量为1 mL，空气流速为400 mL/min，H2流速为30 mL/min，氮气流速为25 mL/min，分流比为20∶1。SF6送北京中科院大气物理研究所中国生态系统研究网络大气分中心代为检测，SF6标准气体浓度为1.99×10-4mg/m3。SF6测定试验条件：载气为氩气(Ar)-甲烷(10%)，载气流速为30 mL/min，电子捕获检测器(ECD)温度为350 ℃，分析柱温度为55 ℃，进样量为5 mL。反刍动物甲烷排放量计算公式如下：

RCH4=(RSF6/6.518)×([CH4]/[SF6])×103。

式中：RCH4是甲烷排放速率(L/d)；RSF6是SF6渗透速率(mg/d)；6.518是SF6的密度(kg/m3)；[CH4]是样品气体中甲烷浓度(×10-6摩尔分数)；[SF6]是样品气体中SF6浓度(×10-12摩尔分数)。

1.7 未知样品拷贝数

未知样品拷贝数计算公式[15]如下：

未知样品拷贝数=[10(Ct1-50.063)/-3.543/10(Ct2-49.217)/-3.029]×1010。

式中：Ct1为目的基因Ct值；Ct2为内参基因Ct值；指数(10)为胶回收产物中目的基因标准曲线拷贝数与内参基因标准曲线拷贝数的比值。

1.8 统计数据

应用SPSS 19.0统计软件的单因素方差分析对试验数据进行统计分析，试验结果用“平均值±标准差”表示。

2 结果与分析

2.1 饲粮NDF/淀粉比值对荷斯坦公牛甲烷排放量的影响

由表5可知，Ⅰ组甲烷排放量较Ⅱ、Ⅲ组极显著降低了33.13%、67.89%(P<0.01)。

表5 荷斯坦公牛甲烷排放量Table 5 Methane production of Holstein bulls L/d

2.2 饲粮NDF/淀粉比值对荷斯坦公牛瘤胃古菌的影响

由图7可知，大部分样本量的稀释曲线在序列数为60 000～80 000时趋于平坦，说明测序数据量渐进合理，更多测序数据对于产生新的OTUs边际效应较小。这说明该样本量可反映瘤胃古菌的真实情况。

由表6可知，各组的覆盖率(goods-coverage)均为1.000，说明测序深度高，足以反映古菌的菌群多样性和丰富度。Ⅲ组Chao1指数显著高于Ⅰ组(P<0.05)，系统发育树(PD_whole_tree)极显著高于Ⅰ组(P<0.01)。

表6 Alpha多样性指数Table 6 Alpha diversity index

由表7可知，各组甲烷短杆菌属(Methanobrevibacter)在瘤胃古菌中占绝对优势，其相对丰度≥97%。Ⅲ组甲烷短杆菌属相对丰度高于Ⅰ、Ⅱ组(P<0.05)。

表7 荷斯坦公牛属水平瘤胃古菌组成 Table 7 Composition of rumen archaea of Holstein bulls at genus level %

Beta多样性是对不同样品的微生物菌落进行比较分析，非加权组平均法(unweighted pair-group method with arithmetic mean，UPGMA)属于Beta多样性的一种分析方法，可以用于研究不同样品的相似性，构建样品的聚类树，对样品进行聚类分析。由图8可知，Ⅰ、Ⅱ组样品瘤胃古菌菌群相似度较高，故聚于一类。

Lg1.1～Lg1.8、Lg2.1～Lg2.8、Lg3.1～Lg3.8分别为Ⅰ、Ⅱ、Ⅲ组瘤胃液样品稀释曲线。Lg1.1 to Lg1.8, Lg2.1 to Lg2.8 and Lg3.1 to Lg3.8 were the dilution curves of rumen fluid samples in groups Ⅰ, Ⅱ and Ⅲ, respectively.图7 稀释曲线Fig.7 Dilution curve

2.3 饲粮NDF/淀粉比值对荷斯坦公牛瘤胃甲烷菌的功能预测

由图9可知，Ⅰ组荷斯坦公牛瘤胃甲烷菌产甲烷途径有4种，包括CO2/H2还原途径、甲烷生成途径、暗氢氧化途径、氢营养型途径；Ⅱ组荷斯坦公牛瘤胃甲烷菌产甲烷途径有5种，包括CO2/H2还原途径、甲烷生成途径、暗氢氧化途径、氢营养型途径、乙酸发酵途径；Ⅲ组荷斯坦公牛瘤胃甲烷菌产甲烷途径有7种，包括乙酸发酵途径、甲基化合物/H2还原途径、甲酸生成甲烷途径、固氮途径、化能异养、甲基原子团失衡生成甲烷途径、其他途径。这说明各组甲烷产生途径存在差异。

基于功能基因的主成分分析(PCA)表示，如果样品的功能组成越相似，则它们在PCA图中的距离越接近。由图10可知，Ⅰ、Ⅱ组功能组成较为接近。

methanogenesis by_CO2 reduction with H2：CO2/H2还原途径；methanogenesis：甲烷生成途径；dark_hydrogen_oxidation：暗氢氧化途径；hydrogenotrophic_methanogenesis：氢营养型途径；acetoclastic_methanogenesis：乙酸发酵途径；methanogenesis_by_reduction_of_methyl_compounds_with H2：甲基化合物/H2还原途径；methanogenesis_using_formate：甲酸生成甲烷途径；nitrogen_fixation：固氮途径；chemoheterotrophy：化能异养 methanogenesis_by_disproporionation of metyl groups：甲基原子团失衡生成甲烷途径；others：其他途径。图9 PICRUSt功能基因预测聚类热图Fig.9 PICRUSt functional genes predicted clustering heat map

图10 基于功能基因的PCAFig.10 PCA based on functional gene

2.4 饲粮NDF/淀粉比值对荷斯坦公牛反刍兽甲烷短杆菌数量的影响

由表8可知，各组之间反刍兽甲烷短杆菌数量存在极显著差异(P<0.01)。

表8 饲粮NDF/淀粉比值对荷斯坦公牛反刍兽甲烷短杆菌数量的影响Table 8 Effects of dietary NDF/starch ratio on Methanobrevibacter ruminantium number of Holstein bulls ×109个

3 讨 论

3.1 饲粮NDF/淀粉比值对荷斯坦公牛甲烷排放量的影响

反刍动物的甲烷排放量约占畜牧业总甲烷排放量的81%，其中90%是由瘤胃微生物甲烷菌产生[22]，甲烷作为饲料发酵的终产物也意味着机体能量损失，如果将这些能量在饲料转化过程中转化为瘤胃发酵产物的能量在能量守恒理论上可用于动物的生产。因此，降低甲烷排放量有助于提高动物的生产效率。

反刍动物甲烷生成量受采食量、饲料种类、饲养水平、饲喂次数等因素的影响[23]。合理调控饲粮精粗比是减少甲烷排放量的关键。研究表明，当饲粮精料比例从0提高到70%，瘤胃甲烷排放量可降低40%[24]。与高纤维饲粮相比，高淀粉饲粮能降低甲烷排放量[25]。丁静美等[10]研究表明，NDF/NFC比值为2.32饲粮组的肉羊甲烷排放量显著高于NDF/NFC比值为1.58和1.04饲粮组。淀粉是NFC主要成分，是机体重要的能量供应来源，淀粉较NFC测定准确，用饲粮NDF/淀粉比值来反映饲粮易发酵程度和瘤胃发酵特性更为直观。本试验研究结果表明，饲粮NDF/淀粉比值为1.63组的荷斯坦公牛甲烷排放量显著高于饲粮NDF/淀粉比值为1.13和1.38组，这与Popova等[8]研究结果一致，其原因是随着饲粮淀粉水平的升高，产甲烷菌的活性会降低，甲烷排放量随之减少。

ck1、ck2、ck3分别代表Ⅰ、Ⅱ、Ⅲ组。下图同。ck1, ck2 and ck3 represented for groups Ⅰ, Ⅱ and Ⅲ, respectively. The same as below.图8 UPGMA聚类树Fig.8 UPGMA cluster tree

3.2 饲粮NDF/淀粉比值对荷斯坦公牛瘤胃古菌多样性的影响

Alpha多样性是用于分析样品内物种的丰富度和多样性[26]。研究表明，纤维摄入能够提高产甲烷菌群多样性，表现为与高淀粉饲粮组比，高纤维饲粮组Shannon指数更高[27-28]，原因是饲粮组成会通过影响产甲烷菌产甲烷途径或基于基因水平的差异对产甲烷过程起重要作用[29]。研究表明，补饲精料组瘤胃甲烷菌群组成发生显著变化。本试验中，饲粮NDF/淀粉比值为1.63组Chao1指数显著高于饲粮NDF/淀粉比值为1.13组(P<0.05)，PD_whole_tree极显著高于饲粮NDF/淀粉比值为1.13组，这与杨硕[7]的研究结果一致。这说明随饲粮NDF/淀粉比值的增加瘤胃古菌物种多样性有所提高，原因可能是高饲粮NDF/淀粉比值促进了纤维分解菌、原虫、纤毛虫、真菌的增殖，促进了以不同底物为生的甲烷菌的生长繁殖和甲烷合成途径的多样性，菌群种类更加丰富；也可能是宿主动物个体差异对瘤胃产甲烷菌产生了较大影响[30]。

饲粮组成是影响瘤胃产甲烷菌种群组成的因素之一[31]。研究表明，高淀粉饲粮能够显著降低甲烷排放量，且瘤胃甲烷排放量与产甲烷菌密切相关[31]。饲粮精料比例增加能够提高优势菌群相对丰度[7]。本试验中，各组瘤胃甲烷短杆菌属在古菌中占绝对优势。饲粮NDF/淀粉比值为1.63组甲烷短杆菌属相对丰度高于其他2组，其他2组之间甲烷短杆菌属相对丰度无显著差异，而各组之间甲烷排放量均存在极显著差异，说明饲粮NDF/淀粉比值不仅对甲烷短杆菌属相对丰度产生了影响，也在种水平数量上产生了较大影响。

3.3 饲粮NDF/淀粉比值对荷斯坦公牛瘤胃甲烷菌的功能预测

PICRUSt能够根据标记物基因(如16S rRNA)对元基因组进行功能预测[32]。金磊等[33]基于高通量测序技术对山羊盲肠进行预测分析发现，盲肠菌群对山羊营养物质消化代谢影响较大。本试验中，各组之间甲烷产生途径存在明显差异，这与UPGMA聚类树、基于功能基因的PCA结果相似，高饲粮NDF/淀粉比值组的荷斯坦公牛瘤胃甲烷菌的甲烷生成途径多于中、低饲粮NDF/淀粉比值组；而从实测甲烷排放量可知，高饲粮NDF/淀粉比值组的甲烷排放量是中、低粮NDF/淀粉比值甲烷排放量之和的1.25倍。这说明瘤胃古菌多样性提高是甲烷排放量大幅提高的主要诱导因素。

3.4 饲粮NDF/淀粉比值对荷斯坦公牛反刍兽甲烷短杆菌数量的影响

高通量测序技术主要基于OTUs进行结果分析，但序列数并不能完全反映所代表细菌相对丰度[34]。研究表明，高通量测序技术可能会高估或低估一些微生物类群的丰度或多样性(如死的生物也被检测到)[35]。一项研究指出，实时定量PCR能够准确检测绵羊瘤胃未知细菌的种群大小及分布情况[36]。因此，实时定量PCR技术可能是检测古菌种数量的重要途径。

相对定量PCR是将待测靶基因片段和内参基因片段同时扩增，就两者Ct值之差进行测量[18]。陈小连[37]对3种瘤胃优势纤维菌绝对定量与相对定量的比较发现，2种方法都能准确测定瘤胃微生物数量。反刍兽甲烷短杆菌属于甲烷短杆菌属，能够从奶牛、公牛、鹿[38-39]瘤胃中分离得到。甲烷短杆菌属属于甲烷杆菌科，是氢营养型甲烷菌[40]，对甲烷产生贡献率较大[41]，且在荷斯坦公牛瘤胃古菌中为优势甲烷菌。因此，本试验采用绝对定量方法对反刍兽甲烷短杆菌数量进行检测。从绝对定量分析结果发现，各组之间反刍兽甲烷短杆菌数量有极显著差异。这说明饲粮NDF/淀粉比值对种水平甲烷短杆菌具有显著影响。杨硕[7]研究表明，采食高纤维含量牧草的绒山羊瘤胃以乙酸发酵为主，从而生成大量的H2，导致以此为生长底物的甲烷菌大量繁殖，且瘤胃产甲烷菌数量随纤维降解菌的升高而增多。本试验中，低饲粮NDF/淀粉比值组荷斯坦公牛瘤胃反刍兽甲烷短杆菌数量显著低于其他2组，说明低NDF/淀粉比值的饲粮改变了荷斯坦公牛瘤胃发酵模式，瘤胃丙酸利用菌与产甲烷菌竞争H2产生丙酸，减少了产甲烷菌可利用底物，导致反刍动物甲烷菌株受到抑制。

4 结 论

降低饲粮NDF/淀粉比值能够减少甲烷排放量，降低瘤胃古菌丰富度指数，降低甲烷短杆菌属相对丰度和反刍兽甲烷短杆菌数量，并减少了荷斯坦公牛瘤胃甲烷菌生成甲烷途径。