任文義 李婭楠 楊雙鳴 張力莉 徐曉鋒

（寧夏大學(xué)農(nóng)學(xué)院，寧夏 銀川 750021）

飼料成本是制約現(xiàn)代養(yǎng)殖業(yè)發(fā)展的重要因素之一[1]。我國的稻草資源豐富，但稻草的品質(zhì)差，反芻動物利用效率低。稻草主要由纖維素和半纖維素組成，還含有木質(zhì)素、含氮化合物和粗灰分等[2]。若能找到合適的處理方法，能夠提高稻草的利用率。

微貯技術(shù)是利用一些酶制劑和益生菌，在合適條件下對秸稈進行厭氧發(fā)酵處理，此過程產(chǎn)生熱量、酸和酶；產(chǎn)生熱量會提高發(fā)酵的溫度，酸可以使稻草柔軟蓬松，降低發(fā)酵pH 值，較低的pH 值和較高的溫度可為纖維素分解酶、半纖維素分解酶和木質(zhì)素分解酶提供適宜的酶解環(huán)境[3]。趙林果等[4]發(fā)現(xiàn)，利用白腐菌及黑曲霉所產(chǎn)的纖維素復(fù)合酶處理稻草48 h，能夠使稻草中的纖維素、半纖維素和木質(zhì)素分別降解53.8%、57.8%、44.5%。微貯飼料化加工可以最大程度保存稻草的營養(yǎng)價值，在一定程度上提高稻草適口性。植物乳桿菌是目前微貯飼料生產(chǎn)中最常用的微生物制劑之一。乳酸菌可以抑制酵母菌和霉菌的生長，避免微貯飼料好氧降解造成的營養(yǎng)損失[5-7]。地衣芽孢桿菌是一種常見于土壤的革蘭氏陽性嗜熱細菌[8]，具有獨特的生物氧捕獲機制，能夠快速創(chuàng)造厭氧環(huán)境，產(chǎn)生抗菌活性物質(zhì)，抑制病菌的生長和繁殖[9]。但地衣芽孢桿菌與植物乳桿菌組合利用的研究相對較少，對微貯的作用有待進一步研究。本試驗旨在研究地衣芽孢桿菌與植物乳桿菌復(fù)合微貯稻草，基于16S rRNA技術(shù)找到發(fā)酵過程中細菌群落動態(tài)和預(yù)測的代謝途徑，分析微貯過程中菌群變化，探討其分子機制，為實際生產(chǎn)中稻草資源的利用提供參考。

1 材料與方法

1.1 試驗材料

稻草采自寧夏，干物質(zhì)（DM）含量為94.49%，粗蛋白質(zhì)（CP）含量為3.45%，中性洗滌纖維（NDF）含量為59.51%，酸性洗滌纖維（ADF）含量為33.59%。

植物乳桿菌為凍干型粉劑，活菌數(shù)為1.39×1010CFU/g，地衣芽孢桿菌為凍干型粉劑，活菌數(shù)為1.30×1011CFU/g，均購自山東中科嘉億生物工程有限公司。

1.2 試驗設(shè)計

試驗采用完全隨機設(shè)計，微貯稻草隨機分為5 組，每組6個重復(fù)。試驗設(shè)計見表1。試驗期30 d。

表1 微貯稻草試驗設(shè)計

1.3 微貯稻草的制作

以粉碎2～3 cm 左右的稻草為原料，利用小規(guī)模發(fā)酵法進行微貯，蔗糖、菌劑、蒸餾水均按照表1進行配制，使用噴壺均勻噴灑至稻草樣品，在乙醇消毒的塑料容器中充分混勻。裝入35 cm × 45 cm的聚乙烯發(fā)酵袋中，每袋500 g，每組6個重復(fù)，使用真空包裝機抽真空，密封。室溫條件（28±3）℃下避光貯藏。于微貯第30 d 后開袋，測定營養(yǎng)指標(biāo)。

1.4 測定指標(biāo)及方法

1.4.1 營養(yǎng)成分

微貯30 d后的稻草開袋取樣，65 ℃烘干48 h，粉碎，過40 目標(biāo)準(zhǔn)篩。采用烘干法（GB/T 6435—2014）測定干物質(zhì)（DM），采用凱氏定氮法（沛歐SKD-100）（GB/T 6432—2018）測定粗蛋白（CP）含量，采用Van Soest 洗滌法測定中性洗滌纖維（NDF）、酸性洗滌纖維（ADF）含量。樣品開袋后，混勻，從每袋中取微貯稻草20 g，加入180 mL蒸餾水，使用多功能攪拌機均質(zhì)30 s，使用4 層紗布過濾至三角瓶中，使用雷磁臺式酸度計測定濾液pH值。

1.4.2 隸屬函數(shù)分析

將DM、CP、NDF、ADF、pH 值采用模糊數(shù)學(xué)的隸屬函數(shù)理論，綜合評價其微貯品質(zhì)[10-11]。其中DM、CP為正隸屬函數(shù)，NDF、ADF、pH值為反隸屬函數(shù)。隸屬函數(shù)均值越大，說明微貯稻草品質(zhì)越好。

1.4.3 微生物測序及生物信息學(xué)分析

根據(jù)模糊隸屬函數(shù)分析將最優(yōu)組合即復(fù)合菌組Ⅲ的1、3、7、15、30 d微貯樣品進行細菌變化規(guī)律研究及功能預(yù)測分析。樣品的DNA提取、PCR擴增、構(gòu)建文庫及測序由廣州基迪奧生物科技有限公司完成，測序平臺為Illumina Hiseq 2500測序平臺。測序得到的原始數(shù)據(jù)經(jīng)過質(zhì)控得到有效數(shù)據(jù)，進行后續(xù)分析。

1.5 數(shù)據(jù)統(tǒng)計與分析

試驗數(shù)據(jù)采用Excel 2019 處理，SAS 9.4 統(tǒng)計軟件進行單因素方差分析，Duncan's 法進行多重比較，結(jié)果以平均值和標(biāo)準(zhǔn)誤表示，P＜0.05表示差異顯著。

根據(jù)OTU 的物種注釋和豐度信息，使用PICRUSt 2軟件對細菌的KEGG 進行功能注釋，并統(tǒng)計每個通路豐度信息。使用Origin 2022軟件作圖。

2 結(jié)果與分析

2.1 地衣芽孢桿菌與植物乳桿菌復(fù)合微貯對稻草營養(yǎng)成分的影響（見表2）

表2 地衣芽孢桿菌與植物乳桿菌復(fù)合微貯對稻草營養(yǎng)成分的影響

由表2 可知，微貯30 d 后，植物乳桿菌組、復(fù)合菌組Ⅰ、復(fù)合菌組Ⅱ、復(fù)合菌組Ⅲ微貯稻草的NDF 含量分別比對照組降低了8.6%、4.6%、6.3%、5.6%（P＞0.05），ADF含量分別比對照組降低了5.9%、3.9%、5.8%、6.8%（P＞0.05）。與對照組相比，復(fù)合菌組Ⅰ、復(fù)合菌組Ⅱ、復(fù)合菌組Ⅲ微貯稻草的pH 值分別降低了6.2%、11.9%、10.6%（P＜0.05）。

2.2 地衣芽孢桿菌與植物乳桿菌復(fù)合微貯對稻草品質(zhì)的隸屬函數(shù)值的影響（見表3）

表3 地衣芽孢桿菌與植物乳桿菌復(fù)合微貯對稻草品質(zhì)的隸屬函數(shù)值的影響

由表3 可知，復(fù)合菌組Ⅲ稻草微貯的隸屬函數(shù)平均值綜合表現(xiàn)最優(yōu)。

2.3 地衣芽孢桿菌與植物乳桿菌復(fù)合微貯對稻草發(fā)酵過程中細菌Alpha多樣性的影響（見表4）

表4 地衣芽孢桿菌與植物乳桿菌復(fù)合微貯對稻草發(fā)酵過程中細菌Alpha多樣性的影響

由表4 可知，微貯第7 d 時，微貯稻草的Shannon 指數(shù)最高，微貯第30 d 的Shannon 指數(shù)顯著低于第1、3、7、15 d（P＜0.05）。微貯第7 d 時，微貯稻草的Simpson指數(shù)最低，微貯第7、15 d的Simpson指數(shù)顯著低于第1 d（P＜0.05）。研究表明，稻草微貯過程中細菌群落多樣性在第7 d時最高，第30 d時最低。

微貯稻草的Chao 1和Ace指數(shù)呈先升高后降低趨勢，第7 d最高，微貯第7、15 d的Chao 1、Ace指數(shù)顯著高于第1 d（P＜0.05）。研究表明，稻草微貯過程中細菌群落豐富度呈先升高后降低趨勢，第7 d時豐富度最高。

2.4 地衣芽孢桿菌與植物乳桿菌復(fù)合微貯對稻草發(fā)酵過程中細菌組成的影響（見表5、圖1）

圖1 稻草發(fā)酵過程中細菌門水平和屬水平和弦圖

表5 稻草發(fā)酵過程中細菌群落的相對豐度 單位：%

由表5、圖1可知，在門水平上，稻草發(fā)酵過程中的優(yōu)勢菌門為厚壁菌門、變形菌門、放線菌門；厚壁菌門隨發(fā)酵時間增加而增加；變形菌門隨發(fā)酵時間增加而減少；放線菌門呈先增加后減少趨勢，在第7 d時相對豐度最高。在屬水平上，稻草發(fā)酵過程中的優(yōu)勢菌屬為乳桿菌屬、泛菌屬、芽孢桿菌屬、魏斯氏菌屬；發(fā)酵第1～3 d，乳桿菌屬和泛菌屬增加，芽孢桿菌屬和魏斯氏菌屬減少；第3～7 d，乳桿菌屬繼續(xù)增加達到最高，而泛菌屬、芽孢桿菌屬開始減少，魏斯氏菌屬開始增加；第7～30 d，乳桿菌屬減少，泛菌屬、芽孢桿菌屬和魏斯氏菌屬增加。

2.5 稻草發(fā)酵過程中微生物功能預(yù)測（見圖2）

圖2 稻草發(fā)酵過程中微生物功能預(yù)測

用PIRUSt的方法與KEGG數(shù)據(jù)庫比對預(yù)測細菌功能選取6條Level 2層級代謝通路和14條Level 3層級代謝通路繪制豐度圖。由圖2 可知，稻草發(fā)酵過程主要涉及的Level 2 層級功能通路包括碳水化合物、氨基酸、輔酶因子和維生素、脂質(zhì)和能量代謝等；Level 3 層級功能通路包括脂肪酸生物合成、D-丙氨酸代謝和磷酸戊糖途徑等。稻草發(fā)酵過程中，溶酶體、鞘脂類代謝、其他多糖降解、糖胺聚糖降解、磷酸鹽和磷酸鹽代謝等途徑存在明顯差異。

3 討論

3.1 地衣芽孢桿菌與植物乳桿菌復(fù)合微貯對稻草營養(yǎng)成分的影響

在微貯飼料中加入添加劑有助于乳酸菌的增殖，阻礙需氧菌和厭氧菌的增加，從而降低飼料pH 值，減少DM損失[12]。本研究中，與對照組相比，稻草微貯中單獨添加植物乳桿菌并未顯著降低發(fā)酵pH值；與對照組和植物乳桿菌組相比，地衣芽孢桿菌與植物乳桿菌復(fù)合后顯著降低了稻草發(fā)酵pH值，可能與地衣芽孢桿菌獨特的生物氧捕獲機制快速創(chuàng)造了厭氧環(huán)境有關(guān)[8]。pH 值是評估發(fā)酵飼料發(fā)酵程度的重要參數(shù)，較低的pH值可確保更好的有氧穩(wěn)定性，并抑制進一步發(fā)酵，pH值降至低于4.5，有助于抑制病原微生物。范子豪[13]研究發(fā)現(xiàn)，地衣芽孢桿菌在酸性環(huán)境中可以正常生長，在有氧狀態(tài)下的生長速度遠大于無氧狀態(tài)，在無氧狀態(tài)下可以利用營養(yǎng)物質(zhì)產(chǎn)生乳酸。因此，在微貯飼料中加入地衣芽孢桿菌可以快速消耗氧氣，有助于乳酸菌生長，產(chǎn)生的乳酸可以快速降低發(fā)酵pH值。

Mara 等[14]研究發(fā)現(xiàn)，地衣芽孢桿菌具有很強的降解纖維的能力，具有促進ADF 分解的作用。本試驗中，添加地衣芽孢桿菌與植物乳桿菌組微貯稻草的NDF和ADF含量均低于對照組，原因可能是地衣芽孢桿菌與植物乳桿菌在降解植物纖維成分方面具有協(xié)同作用。NDF 和ADF 有助于生成可溶性碳水化合物，是一種良好的乳酸菌生長促進劑底物，可促進更多乳酸產(chǎn)生，促進發(fā)酵過程。此外，結(jié)構(gòu)碳水化合物組分的分解可以提高微貯飼料的消化率。

3.2 地衣芽孢桿菌與植物乳桿菌復(fù)合微貯對稻草發(fā)酵過程中菌群變化規(guī)律的影響

微貯飼料微生物是影響其發(fā)酵品質(zhì)的關(guān)鍵因素，原料表面附著的微生物還是外源添加的微生物在微貯發(fā)酵過程中均對微生物群落的演替起重要作用[15]。本試驗中，稻草微貯細菌群落的Shannon、Chao1 和Ace 指數(shù)在發(fā)酵7 d時最高，表明此時稻草發(fā)酵中細菌群落豐富度和多樣性最高。在本研究中，厚壁菌和變形菌是稻草微貯飼料中豐度最高的前兩個門，與之前的研究結(jié)果一致[16-19]。厚壁菌是革蘭氏陽性菌，主要包括產(chǎn)孢子菌、不產(chǎn)孢子菌和支原體菌，能夠降解大分子物質(zhì)，如纖維素、淀粉、蛋白質(zhì)等，增加厚壁菌的豐度會提高產(chǎn)酸能力，從而降低pH值，進一步分解ADF[20]。Li等[21]研究發(fā)現(xiàn)，使用乳酸菌和糖蜜對本地牧草進行微貯，在微貯的早期階段（前7 d），乳酸菌的相對豐度增加，隨著微貯時間的延長，細菌群落的豐富度下降。Ni 等[22]研究發(fā)現(xiàn)，大豆在微貯前7 d乳酸菌數(shù)量顯著增加，在14 d后減少。本研究中，稻草在微貯前7 d，乳桿菌屬豐度隨著發(fā)酵時間的增加而增加，可能是因為地衣芽孢桿菌與植物乳桿菌復(fù)合微貯降低了pH值，抑制了有害微生物的生長，促進了乳桿菌的生長。添加植物乳桿菌可以直接增加乳桿菌屬的豐度，而地衣芽孢桿菌可以通過其獨特的生物氧捕獲機制，快速創(chuàng)造厭氧環(huán)境，促進乳桿菌屬的快速生長繁殖[23]。魏斯氏菌是發(fā)酵飼料中最常見的異型發(fā)酵細菌之一，能夠產(chǎn)生與乳桿菌屬等量的有機酸。研究表明，泛菌屬、芽孢桿菌屬和魏斯氏菌屬等是發(fā)酵飼料中的主要菌屬，但其豐度與環(huán)境和地理位置密切相關(guān)[24-26]。

3.3 地衣芽孢桿菌與植物乳桿菌復(fù)合微貯對稻草發(fā)酵過程中菌群功能預(yù)測的影響

乳酸菌將稻草中可用的碳水化合物發(fā)酵為丙酮酸，在厭氧條件下丙酮酸在胞質(zhì)中還原生成乳酸，完成糖的無氧氧化。當(dāng)飼料中的葡萄糖和果糖水平隨著植物呼吸或微生物利用而下降時，只有多糖降解可以增加葡萄糖含量。許冬梅[27]研究發(fā)現(xiàn)，添加布氏乳桿菌后，玉米青貯發(fā)酵中期糖胺聚糖降解上調(diào)，且主要貢獻微生物為乳酸桿菌屬。本研究發(fā)現(xiàn)，在發(fā)酵第7 d時糖胺聚糖降解途徑顯著上調(diào)，此時乳酸桿菌屬達到最高。本研究中，發(fā)酵第7 d 時，D-精氨酸和D-鳥氨酸代謝途徑開始上調(diào)，表明此時氨基酸的脫氨作用加強，生成的酮酸開始作為菌群發(fā)酵的碳源。Keshri等[28]研究發(fā)現(xiàn)，氨基酸代謝水平可能反映微貯飼料中細菌種群重新合成氨基酸的能力。磷酸鹽和磷酸鹽代謝途徑的變化可能是由乳酸菌引起，乳酸菌發(fā)酵產(chǎn)生的有機酸可溶解磷復(fù)合物而具有解磷作用[29]。

乳酸菌和低氧壓力產(chǎn)生的有機酸會阻止許多霉菌和有害細菌的生長。霉菌和有害細菌在發(fā)酵飼料生長，對線粒體功能產(chǎn)生毒性。Rezabakhsh 等[30]研究發(fā)現(xiàn)，線粒體和溶酶體功能障礙，細胞內(nèi)腺嘌呤核苷三磷酸（ATP）耗竭和過多的活性氧（ROS）生成是霉菌產(chǎn)生毒素的機制。ROS 生成的增加可以使呼吸鏈和三羧酸循環(huán)受到抑制，破壞細胞結(jié)構(gòu)[31]。硫辛酸是一種二硫醇，可通過其兩個巰基部分保護機體免受氧化應(yīng)激。研究表明，硫辛酸及其還原形式二氫硫辛酸可直接清除ROS 和活性氮（RNS），保護細胞免受氧化應(yīng)激[32-34]。本研究中，硫辛酸的代謝在發(fā)酵7 d上調(diào)，15 d時開始下調(diào)。鞘脂可調(diào)節(jié)生物膜的生物物理特性，鞘脂類代謝的大部分路徑都發(fā)生在溶酶體中。溶酶體中鞘脂及其代謝產(chǎn)物會抑制脂質(zhì)循環(huán)，這直接影響細胞膜（包括線粒體內(nèi)膜）的脂質(zhì)組成[35]。本研究中，溶酶體新陳代謝和鞘脂類代謝途徑在發(fā)酵7 d上調(diào)，15 d時開始下調(diào)。本試驗條件下，稻草發(fā)酵的7～15 d，硫辛酸、溶酶體和鞘脂類的代謝與霉菌的變化有關(guān)。Rocchetti 等[36]研究基于霉菌毒素污染的玉米青貯飼料飼喂奶牛后牛乳中代謝組的變化，結(jié)果表明，鞘脂以及嘌呤和嘧啶衍生的代謝物是受影響最大的化學(xué)類別。許多研究表明，霉菌毒素污染的飼料會破壞鞘脂類代謝[37-38]。

4 結(jié)論

本研究條件下，植物乳桿菌復(fù)合地衣芽孢桿菌微貯稻草發(fā)酵優(yōu)于單獨添加植物乳桿菌的效果。