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

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

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

微貯技術(shù)是利用一些酶制劑和益生菌，在合適條件下對(duì)秸稈進(jìn)行厭氧發(fā)酵處理，此過程產(chǎn)生熱量、酸和酶；產(chǎn)生熱量會(huì)提高發(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)養(yǎng)價(jià)值，在一定程度上提高稻草適口性。植物乳桿菌是目前微貯飼料生產(chǎn)中最常用的微生物制劑之一。乳酸菌可以抑制酵母菌和霉菌的生長(zhǎng)，避免微貯飼料好氧降解造成的營(yíng)養(yǎng)損失[5-7]。地衣芽孢桿菌是一種常見于土壤的革蘭氏陽(yáng)性嗜熱細(xì)菌[8]，具有獨(dú)特的生物氧捕獲機(jī)制，能夠快速創(chuàng)造厭氧環(huán)境，產(chǎn)生抗菌活性物質(zhì)，抑制病菌的生長(zhǎng)和繁殖[9]。但地衣芽孢桿菌與植物乳桿菌組合利用的研究相對(duì)較少，對(duì)微貯的作用有待進(jìn)一步研究。本試驗(yàn)旨在研究地衣芽孢桿菌與植物乳桿菌復(fù)合微貯稻草，基于16S rRNA技術(shù)找到發(fā)酵過程中細(xì)菌群落動(dòng)態(tài)和預(yù)測(cè)的代謝途徑，分析微貯過程中菌群變化，探討其分子機(jī)制，為實(shí)際生產(chǎn)中稻草資源的利用提供參考。

1 材料與方法

1.1 試驗(yàn)材料

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

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

1.2 試驗(yàn)設(shè)計(jì)

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

表1 微貯稻草試驗(yàn)設(shè)計(jì)

1.3 微貯稻草的制作

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

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

1.4.1 營(yíng)養(yǎng)成分

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

1.4.2 隸屬函數(shù)分析

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

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

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

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

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

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

2 結(jié)果與分析

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

3 討論

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

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

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

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

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

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

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

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

4 結(jié)論

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