黃貯對玉米秸稈品質及微生物多樣性的影響

劉鑫陽 田瑞華

摘要 [目的]研究黃貯對玉米秸稈的品質及微生物多樣性的影響。[方法]采用現行有效的測定方法及高通量測序，對黃貯玉米秸稈飼料的營養品質及微生物進行了研究。[結果]經過98 d的黃貯，玉米秸稈的總酸含量提高了31.62%，粗纖維含量降低了39.92%，粗蛋白含量提高了8.53%，細菌群落豐富度和多樣性明顯降低，細菌群落組成由初始的藍藻科某些屬為主導變為以乳酸菌為主導。[結論]黃貯有效改善了秸稈飼料的營養品質及微生物群落組成。

關鍵詞 黃貯;高通量測序;微生物多樣性;營養品質

中圖分類號 S816.6文獻標識碼 A

文章編號 0517-6611（2019）19-0197-04

Abstract [Objective] To study the effects of yellowstorage on the quality and microbial diversity of corn straw. [Method] The nutritional quality and microbial diversity of yellowstorage corn straw feed were analyzed by using current effective methods and high throughput sequencing.[Result] The results showed that the total acid content of maize straw increased by 31.62%， cellulose content decreased by 39.92%， crude protein content increased by 8.53% after 98 days of yellowstorage. The richness and diversity of bacterial community decreased obviously， and the composition of bacterial community changed from initial miscellaneous bacteria to lactobacillusdominated structure distribution. [Conclusion] The nutritional quality and microbial community composition of straw feed were effectively improved by yellowstorage.

Key words Yellowstorage;High throughput sequencing;Microbial diversity;Nutritional quality

隨著人們生活水平的不斷提高，對肉制品的需求也不斷增加，使得畜牧業的發展越來越快，但同樣也面臨著自然環境變化、天然草場產量下降等問題，開發新的能源已成為學者們熱切關注的問題。我國秸稈資源非常豐富，年產量高達7億多t[1]，產量占世界總產量的20%～30%，但用于制作飼料的量不到20%。若將這么多的秸稈焚燒會造成嚴重的環境污染，堆積會造成不必要的浪費，秸稈還田也會增加農作物的倒伏以及病蟲害的危害。這些秸稈如果能被合理利用，制成飼料，不僅解決了資源浪費和秸稈焚燒帶來的環境污染問題，而且解決了畜牧發展中存在的一些問題。玉米秸稈主要是由木質素、纖維素和半纖維素組成，其質地粗硬，營養價值低，影響家畜的消化率和采食量。微生物發酵可以顯著降低玉米秸稈中半纖維素和纖維素的含量[2]，提高飼料作物的消化率[3]。通過微生物發酵既改善了玉米秸稈的適口性，也提高了玉米秸稈的營養品質。

高通量測序技術是近年來新興發展起來的免培養分子生物學技術，又稱新一代測序技術[4]。因為高通量測序具有讀數長、效率高、靈敏度高等特點，可最大限度地保留菌群原有的群落組成和分布特征，是一個研究微生物多樣性的重要手段[5-7]。目前，高通量測序技術在人類基因組、青貯微生物、腸道微生物、土壤微生物、植物根根際及內生菌等多種環境微生物生態的研究中都有應用[8-13]。高通量測序技術雖然已被廣泛應用，但該技術在黃貯中應用得較少，且目前對黃貯飼料中菌群結構的研究基本處于空白。筆者利用高通量測序技術分析了黃貯對玉米秸稈中微生物多樣性的影響，結合現行有效的方法分析了黃貯對玉米秸稈中總酸含量、粗纖維含量及粗蛋白含量的影響。

1 材料與方法

1.1 儀器與設備

超低溫冰箱 DW-86L388A（青島海爾特種電器有限公司）、移液器 Eppendorf N13462C（Eppendorf）、小型離心機 Eppendorf 5430 R（Eppendorf）、高速臺式冷凍離心機 Eppendorf 5424R（Eppendorf）、超微量分光光度計 NanoDrop2000（Thermo FisherScientific）、電泳儀 DYY-6C（北京市六一儀器廠）、MISEQ 測序儀 Illumina Miseq（Illumina）、酶標儀 BioTek ELx800（Biotek）。

1.2 主要試劑

DNA 抽提試劑盒、2% agarose gels、FastPfu Polymerase、AxyPrep DNA GelExtraction Kit、MiSeqplatform、FastDNASPINKitforSoil、酚酞、氫氧化鈉、纖維濾袋、濃硫酸、無水葡萄糖、纖維濾袋、混合催化劑、混合指示劑、硼酸、鹽酸、DNS試劑。

1.3 試驗方法

1.3.1 采樣方法。

玉米秸稈黃貯池共3個，均為長、寬、高各2 m的池子，玉米秸稈黃貯時間為98 d，分別在玉米秸稈黃貯封窖前及98 d黃貯結束過后進行取樣。在黃貯池的1 m深處分別取5個樣，將取出的玉米秸稈樣品進行充分混合，然后平均分成3份放入自封袋中，置于-80 ℃冰箱中保存。

1.3.2 黃貯秸稈樣品的品質指標測定。

玉米秸稈黃貯過程中品質指標的測定參照現行有效的分析方法[14-16]。

1.3.3 Illumina Miseq 高通量測序。

1.3.3.1 DNA抽提與PCR擴增。

參照FastDNA SPIN Kit for soil試劑盒說明書對玉米秸稈樣品進行總DNA提取，并利用NanoDrop2000對DNA的純度和濃度進行檢測，利用806R（5-GGACTACHVGGGTWTCTAAT-3）和338F（5-ACTCCTACGGGAGGCAGCAG-3）引物對V3～V4可變區進行PCR擴增。

1.3.3.2 Illumina Miseq測序。

利用Illumina公司的Miseq PE300平臺進行測序（上海美吉生物醫藥科技有限公司），將原始數據上傳至NCBI網站進行比對。

1.3.3.3 數據處理。

原始測序的序列使用Trimmomatic 軟件進行質控，并使用FLASH軟件進行拼接。

2 結果與分析

2.1 黃貯玉米秸稈的品質指標測定

為了解黃貯玉米秸稈的表型特征，對黃貯玉米秸稈的總酸含量、粗纖維含量、粗蛋白含量這3個重要品質指標進行測定，測定結果見表1。

由表1可以看出，通過黃貯發酵的玉米秸稈總酸含量由初始的8.57 g/kg升高到11.28 g/kg。與未進行黃貯發酵的玉米秸稈相比，經過黃貯發酵的玉米秸稈總酸含量提高了31.62%，這說明黃貯顯著提高了秸稈中的總酸含量;粗纖維含量由初始的34.47%降低到20.71%，與未進行黃貯發酵的玉米秸稈相比，經過黃貯發酵后的玉米秸稈粗纖維含量降低了39.92%，說明黃貯顯著降低了玉米秸稈中粗纖維的含量;粗蛋白含量由初始的3.87%提高到4.20%，與未進行黃貯發酵的玉米秸稈相比，經過黃貯發酵的玉米秸稈粗蛋白含量提高了8.53%，說明黃貯有效提高了玉米秸稈中粗纖維的含量。

2.2 玉米秸稈黃貯的 Illumina Miseq 高通量測序

2.2.1 黃貯玉米秸稈樣品中細菌16S rRNA的Illumina測序的稀釋曲線。

稀釋曲線（rarefaction curve）主要利用各樣本在不同測序深度時的微生物α多樣性指數構建曲線，以此反映各樣本在不同測序數據量時的微生物多樣性，也可以用來說明樣本測序數據量的合理性。稀釋曲線采用對序列進行隨機抽樣的方法，以抽到的序列數與它們對應的物種（OTU）數目或多樣性指數，構建稀釋曲線。該研究采用Shannon-Wiener指數，選擇97%相似度的OTU水平，利用mothur計算不同隨機抽樣下的Shannon-Wiener多樣性指數，利用R語言工具繪制稀釋曲線。當曲線趨向平坦時，說明測序數據量合理，更多的數據量只會產生少量新的物種，反之則表明繼續測序還可能產生較多的新物種。

圖1為黃貯玉米秸稈樣品的細菌稀釋曲線，隨著樣本序列數的增加，檢測到OTU數量的增加量逐漸減小，曲線趨于平坦，這說明該試驗測序的數據量比較合理，通過檢測獲得了樣品細菌群落中絕大多數OTU水平上的物種。

2.2.2 黃貯玉米秸稈樣品中細菌的α-多樣性。

通過chao指數、Shannon-Wiener指數、物種覆蓋度指數對黃貯玉米秸稈樣品中細菌的α多樣性進行分析。chao指數是描述物種豐富度chao1指數，是用chao1算法估計樣本中所含OTU數目的指數，用來估計物種總數;Shannon-Wiener指數是反映多樣性的指數，其值越高說明群落多樣性越高;物種覆蓋度指數（coverage）是反映群落覆蓋度的指數，其數值越高則樣本中序列被測出的概率越高，而沒有被測出的概率越低，該指數反映此次測序結果是否代表樣本中微生物的真實情況。

從表2可以看出，黃貯玉米秸稈樣品中細菌的物種覆蓋度均在0.99以上，說明樣本中序列被測出的概率很高，且此次試驗數據代表了樣本中微生物的真實情況。經過98 d黃貯后，3個黃貯池中玉米秸稈細菌的Shannon-Wiener指數和chao指數均顯著降低，說明黃貯降低了玉米秸稈中細菌的豐富度和多樣性。

2.2.3 黃貯玉米秸稈樣品中細菌群落分類學分析。

2.2.3.1 黃貯玉米秸稈中細菌群落在門水平上的結構分布。

通過與Silva數據庫進行對比發現，玉米秸稈初始的細菌群落由19個門組成，經過98 d黃貯后玉米秸稈中細菌群落變為18個門。

從圖2可以看出，玉米秸稈初始的細菌群落中藍藻門（Cyanobacteria）和變形菌門（Proteobacteria）占據明顯的優勢地位，此外還有擬桿菌門（Bacteroidetes）和放線菌門（Actinobacteria）等一些相對豐度較低的細菌門類。黃貯98 d后的細菌群落中厚壁菌門（Firmicutes）和變形菌門占據明顯的優勢地位，此外還有擬桿菌門和放線菌門等相對豐度較低的門類。這說明黃貯明顯改變了細菌在門水平上的結構分布。

2.2.3.2 黃貯玉米秸稈中細菌群落在屬水平上的結構分布。

通過與Silva數據庫進行對比發現，玉米秸稈初始的細菌群落由349個屬組成，經過98 d黃貯后玉米秸稈中細菌群落變為323個屬。

從圖3可以看出，玉米秸稈初始的細菌群落中藍藻科的某些屬占據明顯的優勢地位，此外還有泛菌屬、鞘脂單胞菌屬、根瘤菌屬等一些相對豐度較低的屬。黃貯98 d后的細菌群落中，乳桿菌屬代替了藍藻科某些屬的優勢地位，成為優勢菌屬，此外還有片球菌、乳球菌等一些相對豐度較低的菌屬。由黃貯玉米秸稈中細菌群落在屬水平上的結構分布可以看出，經過98 d的黃貯，秸稈飼料中細菌群落的組成由一些雜菌變為乳酸菌主導。乳酸菌分為同型發酵乳酸菌、異型發酵乳酸菌和兼性異型發酵乳酸菌。其中，同型發酵乳酸菌可以快速降低pH、增加乳酸菌數量和乳酸含量;專性異型發酵乳酸菌具有增加有氧穩定性的作用;兼性異型發酵乳酸菌具有二者聯合的作用[17]。

3 結論

總酸含量是衡量玉米秸稈黃貯發酵秸稈品質的一個重要指標，相對較高的總酸含量不僅可以抑制黃貯發酵秸稈中雜菌的生長，而且可以增加黃貯發酵秸稈的有氧穩定性。玉米秸稈飼料中纖維素的存在會影響飼料的適口性，含有很多不能被動物消化和利用的物質，因此降低玉米秸稈中的纖維含量是必要的;黃貯飼料中蛋白的存在可以有效改善飼料的營養價值，從而提高飼料的飼用價值，因此飼料中蛋白含量是衡量黃貯飼料品質的一個重要指標。Huang等[18]研究發現秸稈經青貯后總酸的含量顯著提高，有效抑制了秸稈中有害菌的生長，提高飼料的品質。Hu等[19]、Nsereko等[20]、Zayed等[21]研究發現青貯可以顯著提高秸稈飼料的干物質、粗蛋白等含量以及纖維素的降解率。該研究發現玉米秸稈經黃貯后有效提高了玉米秸稈的總酸及粗蛋白的含量，顯著降低了秸稈的粗纖維含量，提高了玉米秸稈的營養品質，這與上述研究結果相似。

經過黃貯發酵，玉米秸稈中的細菌豐富度和多樣性都有所降低。Li等[22]對丁香和黑麥草進行青貯，結果表明青貯降低了丁香和黑麥草飼料的細菌豐富度，與該研究結果相似。Ogunade等[9]研究發現苜蓿青貯飼料細菌群落中74.1%和20.4%的細菌分別屬于厚壁菌門和變形菌門。Ogunade等[23]對玉米秸稈進行青貯，通過對青貯飼料的細菌多樣性分析表明青貯可以改變飼料中的細菌群落組成。結果表明，玉米青貯飼料中涉及乳酸發酵的細菌群落中，98.3%屬于厚壁菌門，96.5%屬于乳酸桿菌屬，與該研究結果相似。

通過對黃貯玉米秸稈飼料的品質指標進行分析，發現黃貯有效改善了秸稈飼料中總酸、粗纖維和粗蛋白的含量。通過對黃貯玉米秸稈飼料的微生物進行分析，發現黃貯顯著降低了秸稈飼料中細菌群落的豐富度和多樣性，并明顯改變了玉米秸稈飼料中細菌群落的結構組成。

參考文獻

[1] 霍麗麗，吳娟娟，趙立欣，等.華北平原地區玉米秸稈連續供應模型的建立及應用[J].農業工程學報，2016，32（19）：203-210.

[2] GUO H W，CHANG J，YIN Q Q，et al.Effect of the combined physical and chemical treatments with microbial fermentation on corn straw degradation[J].Bioresource technology，2013，148：361-365.

[3] HUNT C W，KEZAR W，HINMAN D D，et al.Effects of hybrid and ensiling with and without a microbial inoculant on the nutritional characteristics of wholeplant corn[J].Journal of animal science，1993，71（1）：38-43.

[4] 令利軍，何楠，白雪，等.基于高通量測序的玉米秸稈自然發酵過程中細菌菌群結構特征[J].蘭州大學學報（自然科學版），2017，53（4）：526-533.

[5] 王玉榮，折米娜，劉康玲，等.基于MiSeq高通量測序技術內蒙古地區酸粥細菌多樣性研究[J].食品工業科技，2018，39（19）：124-129.

[6] 魏軍，趙志軍.下一代測序技術在分子診斷中的應用[J].分子診斷與治療雜志，2013，5（3）：145-151.

[7] HOLT K E，PARKHILL J，MAZZONI C J，et al.Highthroughput sequencing provides insights into genome variation and evolution in Salmonella Typhi[J].Nature genetics，2008，40（8）：987-993.

[8] REUTER J A，SPACEK D V，SNYDER M P.Highthroughput sequencing technologies[J].Molecular cell，2015，58（4）：586-597.

[9] OGUNADE I M，JIANG Y，PECH CERVANTES A A，et al.Bacterial diversity and composition of alfalfa silage as analyzed by Illumina MiSeq sequencing：Effects of，Escherichia coli，O157：H7 and silage additives[J].Journal of dairy science，2018，101（3）：2048-2059.

[10] LONG L L，GUO J J，LI P，et al.Bacterial diversity in Bercaea Cruentata gut described using highthroughput sequencing[J].Forensic science international genetics supplement series，2015，5：e479-e481.

[11] DUAN D Y，LIU G H，CHENG T Y，et al.Microbial population analysis of the midgut of Melophagus ovinus via highthroughput sequencing[J].Parasites & vectors，2017，10（1）：1-7.

[12] FAN T L，SUN Y X，PENG J J，et al.Combination of amplified rDNA restriction analysis and highthroughput sequencing revealed the negative effect of colistin sulfate on the diversity of soil microorganisms[J].Microbiological research，2018，206：9-15.

[13] HUANG Y H.Comparison of rhizosphere and endophytic microbial communities of Chinese leek through highthroughput 16S rRNA gene Illumina sequencing[J].Journal of integrative agriculture，2018，17（2）：359-367.

[14] 中國食品發酵工業研究院.食品中總酸的測定：GB/T 12456—2008[S].北京：中國標準出版社，2008.

[15] 中國農業科學院農業質量標準與檢測技術研究所（國家飼料質量監督檢驗中心）.飼料中粗蛋白的測定 凱氏定氮法：GB/T 6432—2018[S].北京：中國標準出版社，2018.

[16] 農業部飼料質量監督檢驗測試中心.飼料中粗纖維的含量測定 過濾法：GB/T 6434—2006[S].北京：中國標準出版社，2006.

[17] ?；溃＝埽诟?，等.青貯微生物的發掘及其應用的研究進展[J].動物營養學報，2018，30（11）：4279-4285.

[18] HUANG J L，WANG L K，DAI S F.Effects of previously fermented juice on nutritive value and fermentative quality of rice straw silage[J].Journal of northeast agricultural university，2013，20（2）：48-52.

[19] HU W，SCHMIDT R J，MCDONELL E E，et al.The effect of Lactobacillus buchneri 40788 or Lactobacillus plantarum MTD1 on the fermentation and aerobic stability of corn silages ensiled at two dry matter contents[J].Journal of dairy science，2009，92（8）：3907-3914.

[20] NSEREKO V L，SMILEY B K，RUTHERFORD W M，et al.Influence of inoculating forage with lactic acid bacterial strains that produce ferulate esterase on ensilage and ruminal degradation of fiber[J].Animal feed science and technology，2008，145（1/2/3/4）：122-135.

[21] ZAYED M S.Enhancement the feeding value of rice straw as animal fodder through microbial inoculants and physical treatments[J].International journal of recycling of organic waste in agriculture，2018，7（2）：117-124.

[22] LI P，ZHANG Y，GOU W L，et al.Silage fermentation and bacterial community of bur clover，annual ryegrass and their mixtures prepared with microbial inoculant and chemical additive[J].Animal feed science and technology，2019，247：285-293.

[23] OGUNADE I M，JIANG Y，KIM D H，et al.Fate of E.coli O157：H7 and bacterial diversity in corn silage contaminated with the pathogen and treated with chemical or microbial additives[J].Journal of dairy science，2017，100（3）：1780-1794.