張建龍，王亞妮，董樂輝，陳國忠*，姜琳琳，于馨，朱洪偉，鄭維波，劉燕

試驗(yàn)研究

短時(shí)高溫好氧固態(tài)發(fā)酵對(duì)豬糞細(xì)菌群落結(jié)構(gòu)的影響

張建龍1,2,3，王亞妮1，董樂輝4，陳國忠1,2,3*，姜琳琳1,2,3，于馨1,2,3，朱洪偉1,2,3，鄭維波1,2,3，劉燕1,2,3

（1.魯東大學(xué)生命科學(xué)學(xué)院，山東 煙臺(tái) 264025；2.山東省養(yǎng)殖環(huán)境控制工程實(shí)驗(yàn)室，山東 煙臺(tái)；3.煙臺(tái)市動(dòng)物病原微生物與免疫學(xué)重點(diǎn)實(shí)驗(yàn)室，山東 煙臺(tái)；4.山東益生種畜禽股份有限公司，山東 煙臺(tái)）

為了研究畜禽糞便固態(tài)發(fā)酵前后微生物的變化，采用高通量測序技術(shù)，通過可操作分類單元分析、α多樣性分析和物種組成分析，比較了豬糞短時(shí)高溫好氧固態(tài)發(fā)酵前后的細(xì)菌群落組成，結(jié)果表明，發(fā)酵前后可操作分類單元總數(shù)分別為151個(gè)和126個(gè)，發(fā)酵后細(xì)菌群落多樣性指數(shù)降低；發(fā)酵前相對(duì)豐度最高的分別是變形菌門59.48 %、擬桿菌門23.97 % 和厚壁菌門13.79 %，發(fā)酵后變形菌門和擬桿菌門分別減少至0.66 % 和1.00 %，而厚壁菌門增加至99.01 %，是發(fā)酵豬糞的優(yōu)勢菌群。本研究表明發(fā)酵后豬糞細(xì)菌群落多樣性下降，抗生素抗性基因潛在攜帶者相對(duì)豐度降低，耐高溫的有機(jī)質(zhì)降解菌成為優(yōu)勢菌屬，該研究對(duì)豬糞的資源化處理具有理論指導(dǎo)意義。

豬糞；高通量測序；細(xì)菌群落結(jié)構(gòu)

集約化養(yǎng)殖導(dǎo)致畜禽糞便大量堆積，嚴(yán)重污染環(huán)境[1]。豬糞中含有大量的病原微生物，處理不當(dāng)會(huì)導(dǎo)致對(duì)環(huán)境的二次污染[2]，暴發(fā)大規(guī)模的傳染性疾病[3]。高溫好氧堆肥可防止抗生素、抗生素抗性基因和抗性細(xì)菌從糞肥釋放到農(nóng)田土壤中[4]，改善土壤，提高農(nóng)作物產(chǎn)量[5]。在高溫好氧堆肥過程中微生物群落結(jié)構(gòu)會(huì)發(fā)生顯著的變化[6]。據(jù)報(bào)道Kato等[7]研究發(fā)現(xiàn)堆肥加速了牛糞中微生物群落的演替，使微生物豐富度下降的同時(shí)保持了微生物的多樣性。Mao等[8]在研究堆肥前后豬糞中微生物的群落演替時(shí)發(fā)現(xiàn)厚壁菌門是嗜熱階段的絕對(duì)優(yōu)勢菌門，保持良好的微生物群落結(jié)構(gòu)對(duì)豬糞堆肥起到了重要的指導(dǎo)性作用。

高通量測序技術(shù)能夠通過從環(huán)境樣本中直接獲取總的DNA進(jìn)行文庫構(gòu)建，從而提高測序效率[9]。本研究采用高通量測序的手段，對(duì)短時(shí)高溫好氧固態(tài)發(fā)酵前后的豬糞進(jìn)行微生物多樣性分析，闡明短時(shí)高溫好氧固態(tài)發(fā)酵前后豬糞細(xì)菌群落及物種的組成和相對(duì)豐度，進(jìn)而可以確定豬糞發(fā)酵過程中的優(yōu)勢菌群，并通過對(duì)比總結(jié)出豬糞發(fā)酵前后細(xì)菌群落的演替規(guī)律，對(duì)豬糞發(fā)酵制備有機(jī)肥的資源化處理技術(shù)具有理論指導(dǎo)意義。

1 材料與方法

1.1 材料

豬糞由養(yǎng)殖戶提供，發(fā)酵菌劑AT酵素由百威國際生物股份有限公司提供，花生殼粉采購自煙臺(tái)市農(nóng)貿(mào)市場。

1.2 方法

1.2.1 短時(shí)高溫好氧固態(tài)發(fā)酵 將豬糞和花生殼粉按比例混拌均勻，碳氮比調(diào)節(jié)為25:1，控制含水量60 %，從不同位點(diǎn)隨機(jī)采樣5份，合并為1份，記為BEFE，作為發(fā)酵前樣本。將物料加入百威高速發(fā)酵設(shè)備（百威國際生物股份有限公司）內(nèi)，每1 t原料添加1 kg發(fā)酵菌劑，啟動(dòng)電源加熱至50 ℃，啟動(dòng)攪拌每1 h轉(zhuǎn)1圈，曝氣速率為15 L/min。發(fā)酵10 h后加熱至80 ℃ 維持1 h，加熱至100 ℃ 再維持1 h，結(jié)束發(fā)酵，出料，在室外建成一個(gè)長2 m、寬1 m、高0.8 m的堆體，每天進(jìn)行1～2次翻堆，腐熟6 d。當(dāng)發(fā)酵結(jié)束后，在堆肥不同位點(diǎn)隨機(jī)采集樣品5份，合并為1份，記為AEFE，作為發(fā)酵后樣本。每個(gè)樣本采集3份，置于-80℃冰箱保存?zhèn)溆谩?/p>

1.2.2 高通量測序 將樣品送北京諾禾致源生物信息科技有限公司進(jìn)行16S rRNA高通量測序，測序引物515F/907R，序列為515F：5'-GTGCCA GCMGCCGCGG-3'，907R：5'-CCGTCAATTCM TTTRAGTTT-3'，測序的區(qū)域是V4-V5區(qū)，測序平臺(tái)是Illumina Miseq。

1.2.3 數(shù)據(jù)處理 對(duì)得到的測序原始數(shù)據(jù)進(jìn)行篩選處獲得優(yōu)化序列，然后在97 % 的相似水平上進(jìn)行可操作分類單元聚類分析，用Mothur軟件計(jì)算得到豐富度指數(shù)、均勻度指數(shù)、香農(nóng)指數(shù)、辛普森指數(shù)和覆蓋率，利用R語言工具分析樣本在不同分類水平上的群落結(jié)構(gòu)。用SPSS統(tǒng)計(jì)軟件進(jìn)行單因素方差分析，顯著性差異結(jié)果表示為<0.05。

2 結(jié)果

2.1 豬糞樣品的測序質(zhì)量

通過引物對(duì)樣品的16S進(jìn)行測定，豬糞發(fā)酵前后高通量測序總的堿基數(shù)分別為70 750和80 333，其中合并的堿基數(shù)分別為66 493和75 328，各占所測堿基總數(shù)的93.98 % 和93.77 %。從序列長度來看，所測序列長度分別在264 bp ～ 402 bp和265 bp ～ 394 bp之間，平均測序長度為372 bp和371 bp。

2.2 豬糞樣品菌群的α多樣性

2.2.1 稀釋曲線 為了判斷本樣品測序數(shù)據(jù)量是否合理，筆者對(duì)發(fā)酵前后豬糞的測序序列進(jìn)行隨機(jī)抽樣，然后利用稀釋曲線進(jìn)行評(píng)估。發(fā)酵后和發(fā)酵前的曲線隨著測序深度的增加最終接近平臺(tái)期趨于平穩(wěn)，說明本次取樣的合理性，測序深度已經(jīng)基本覆蓋樣本中的所有物種。

2.2.2 群落多樣性指數(shù) 香農(nóng)指數(shù)和辛普森指數(shù)代表了微生物的多樣性。由表1可以看出發(fā)酵后細(xì)菌群落的香農(nóng)指數(shù)與辛普森指數(shù)相比于發(fā)酵前降低了，表明發(fā)酵后豬糞的物種多樣性下降但是群落分布多樣性卻有微量的增長。發(fā)酵后的豐富度指數(shù)和均勻度指數(shù)明顯低于發(fā)酵前，這說明發(fā)酵后的豬糞中物種豐富度顯著降低。發(fā)酵前細(xì)菌可操作分類單元為151個(gè)，發(fā)酵后細(xì)菌可操作分類單元為126種，兩種狀態(tài)下都存在的為85種，表明發(fā)酵前后的細(xì)菌類別存在較大差異。

表1 群落多樣性指數(shù)

注：*表示差異顯著（<0.05）

2.3 細(xì)菌群落結(jié)構(gòu)變化

2.3.1 在門水平的細(xì)菌群落結(jié)構(gòu) 高通量測序結(jié)果表明（表2），發(fā)酵前豬糞細(xì)菌群落在門水平上的主要類群為變形菌門、擬桿菌門、厚壁菌門、放線菌門，相對(duì)豐度分別為59.48 %、23.97 %、13.79 %、2.55 %。變形菌門是其優(yōu)勢菌門。發(fā)酵后厚壁菌門是優(yōu)勢菌門，并且占有絕對(duì)優(yōu)勢，相對(duì)豐度達(dá)到99.01 %。

表2 發(fā)酵前后細(xì)菌在門水平的相對(duì)豐度

表3 發(fā)酵前后細(xì)菌在屬水平的相對(duì)豐度

2.3.2 在屬水平的細(xì)菌群落結(jié)構(gòu) 由表3可見，發(fā)酵前在屬水平上相對(duì)豐度排名前7的菌屬為不動(dòng)桿菌屬40.85 %、香味菌屬17.58 %、毛球菌屬7.88 %、橄欖桿菌屬4.19 %、弧菌屬3.65 %、假單胞菌屬3.63 %、腸球菌屬2.63 %。發(fā)酵后這7個(gè)屬的相對(duì)豐度值都顯著降低，尤其是優(yōu)勢菌屬不動(dòng)桿菌屬在發(fā)酵后相對(duì)豐度僅為0.66 %，遠(yuǎn)低于發(fā)酵前。嗜熱糞桿菌屬、中度嗜熱微菌屬和中度嗜熱厭氧菌屬成為發(fā)酵后的優(yōu)勢菌屬，占比分別達(dá)到66.95 %、7.41 % 和4.35 %。

2.3.3 細(xì)菌群落各分類單元占比差異 對(duì)發(fā)酵前后豬糞細(xì)菌群落各分類單元的分布差異進(jìn)行了分析，發(fā)酵前細(xì)菌在綱、目、科、屬中的分布廣泛，發(fā)酵后則主要集中于厚壁菌門梭菌綱。發(fā)酵前豬糞中的細(xì)菌主要分布在變形菌門、擬桿菌門和厚壁菌門，發(fā)酵后則是厚壁菌門占據(jù)絕對(duì)優(yōu)勢。對(duì)于厚壁菌門而言，發(fā)酵前僅有芽孢桿菌綱乳桿菌目的肉桿菌科、氣球菌科和腸球菌科。發(fā)酵后的豬糞細(xì)菌群落中有芽孢桿菌綱乳桿菌目的肉桿菌科、腸球菌科、乳酸細(xì)菌科和芽孢桿菌目，以及梭菌綱的梭菌目和鹽厭氧菌目。

3 討論

（1）采用高通量測序技術(shù)研究了短時(shí)高溫好氧固態(tài)發(fā)酵對(duì)豬糞細(xì)菌多樣性的影響。樣品的稀釋曲線末端趨于平坦說明樣本的測序數(shù)據(jù)量是足夠的。發(fā)酵后豬糞中細(xì)菌的多樣性下降，其原因在于短時(shí)高溫好氧發(fā)酵的溫度高于50 ℃ 甚至達(dá)到80～100 ℃，在此階段大多數(shù)微生物因?yàn)槟蜔嵝暂^差而失活[10]，嗜熱細(xì)菌生存下來成為優(yōu)勢菌群[8]。（2）通過種群結(jié)構(gòu)分析，豬糞發(fā)酵前后細(xì)菌群落變化顯著，發(fā)酵前厚壁菌門、變形菌門、放線菌門和擬桿菌門是堆肥微生物區(qū)系中占主導(dǎo)地位的門類，發(fā)酵后厚壁菌門成為了豐度最高的優(yōu)勢菌門且占比超過99 %，該結(jié)果與王旭輝等[11]報(bào)道的牛糞堆肥過程中微生物的變化是一致的。動(dòng)物糞便是人畜共患病原體的主要來源，其中的變形菌門是細(xì)菌中最大的門，含有大量的病原性菌種。發(fā)酵前變形菌門是豬糞細(xì)菌群落的優(yōu)勢菌門，符合畜禽腸道微生物菌落的一般規(guī)律。大多數(shù)病原菌在50 ℃ 以上24 h內(nèi)就會(huì)被滅活，因此豬糞中的病原菌經(jīng)過高溫發(fā)酵后絕大部分被殺死，原來占優(yōu)勢的變形菌門在發(fā)酵后僅占細(xì)菌總數(shù)0.66 %。變形菌門中很多菌種是抗生素抗性基因的攜帶者，變形菌門豐度降低和抗生素抗性基因的消除具有一定的相關(guān)性[12-13]。發(fā)酵后厚壁菌門成為豬糞中的絕對(duì)優(yōu)勢菌群，其原因是多方面的。首先，厚壁菌門細(xì)胞壁的肽聚糖含量約在50 % ～ 80 %之間，細(xì)胞壁較厚具有較強(qiáng)的耐熱性；其次，發(fā)酵菌劑大多為厚壁菌門菌種，其本身就是優(yōu)勢菌群；最后，厚壁菌門能利用現(xiàn)成可發(fā)酵的碳水化合物并參與纖維素外圍側(cè)鏈的初始定植[14]，同時(shí)分泌多種果膠酶來完成降解纖維素等難消化的碳水化合物的目的，發(fā)酵過程中加入的輔料是含有大量木質(zhì)素和纖維素為主要成分的花生殼，對(duì)厚壁菌門具有促生作用。此外，豬糞中氮源含量高，而厚壁菌門具有能夠利用氮源和促進(jìn)氨基酸運(yùn)輸代謝的能力[15]，這也是厚壁菌門成為優(yōu)勢菌群的主要原因之一。（3）從屬水平分析可知，發(fā)酵前不動(dòng)桿菌屬、香味菌屬、毛球菌屬、橄欖桿菌屬、弧菌屬、假單胞菌屬、腸球菌屬7個(gè)相對(duì)豐度較高的菌屬在發(fā)酵后都顯著降低，嗜熱糞桿菌屬、中度嗜熱微菌屬和中度嗜熱厭氧菌屬成為發(fā)酵后的優(yōu)勢菌屬。已有報(bào)道表明畜禽糞便中磺酰胺類抗生素抗性基因的相對(duì)豐度較高[16]，而不動(dòng)桿菌屬和假單胞菌屬被認(rèn)為是主要載體[17]，這兩個(gè)菌屬在發(fā)酵后豐度降低有利于抗生素抗性基因的消減，即有利于畜禽糞便的無害化。嗜熱糞桿菌屬、中度嗜熱微菌屬和中度嗜熱厭氧菌屬在纖維素和有機(jī)物的降解中起主導(dǎo)作用[18-21]，都屬于嗜熱菌，這可能是這三類菌高溫發(fā)酵后在豬糞的微生物菌落中成為優(yōu)勢菌種的主要原因。（4）本研究對(duì)采用短時(shí)高溫好氧固態(tài)發(fā)酵前后的豬糞樣本進(jìn)行了高通量測序，通過可操作分類單元分析、α多樣性分析和物種組成分析發(fā)現(xiàn)：發(fā)酵后豬糞細(xì)菌群落物種多樣性下降，變形菌門、擬桿菌門相對(duì)豐度降低，厚壁菌門相對(duì)豐度增加成為優(yōu)勢菌門，抗生素抗性基因攜帶者不動(dòng)桿菌屬和假單胞菌屬相對(duì)豐度降低，耐高溫的有機(jī)質(zhì)降解菌嗜熱糞桿菌屬、中度嗜熱微菌屬和中度嗜熱厭氧菌屬成為優(yōu)勢菌屬。發(fā)酵前后豬糞細(xì)菌群落結(jié)構(gòu)的變化主要與發(fā)酵溫度、輔料和發(fā)酵菌劑特性有關(guān)，短時(shí)高溫好氧固態(tài)發(fā)酵能夠有利于豬糞無害化處理。

[1]金香琴, 段麗杰, 馬繼力, 等. 不同微生物菌劑對(duì)畜禽糞便資源化堆肥效果的影響[J]. 科學(xué)技術(shù)與工程, 2015, 15（7）: 280-283.

[2] Niu Q, Hojo T, Qiao W, et al. Characterization of methanogenesis, acidogenesis and hydrolysis in thermophilic methane fermentation of chicken manure[J]. Chemical Engineering Journal, 2014, 244: 587-596.

[3] Erickson M C, Liao J, Ma L, et al. Inactivation of Salmonella spp. in cow manure composts formulated to different initial C:N ratios[J]. Bioresource Technology, 2009, 100（23）: 5898-5903.

[4] Lin H, Zhang J, Chen H, et al. Effect of temperature on sulfonamide antibiotics degradation, and on antibiotic resistance determinants and hosts in animal manures[J]. Science of The Total Environment, 2017, 607-608: 725-732.

[5] Huang G F, Wong J W C, Wu Q T, et al. Effect of C/N on composting of pig manure with sawdust[J]. Waste Management, 2004, 24（8）: 805-813.

[6] Huang G, Wang X, Han L. Rapid estimation of nutrients in chicken manure during plant-field composting using physicochemical properties[J]. Bioresource Technology, 2011, 102（2）: 1455-1461.

[7] Kato K, Miura N. Effect of matured compost as a bulking and inoculating agent on the microbial community and maturity of cattle manure compost[J]. Bioresource Technology, 2008, 99（9）: 3372-3380.

[8] Mao H, Lv Z, Sun H, et al. Improvement of biochar and bacterial powder addition on gaseous emission and bacterial community in pig manure compost[J]. Bioresource Technology, 2018, 258: 195-202.

[9] 秦楠, 栗東芳, 楊瑞馥. 高通量測序技術(shù)及其在微生物學(xué)研究中的應(yīng)用[J]. 微生物學(xué)報(bào), 2011, 51（4）: 445-457.

[10] Reuveni R, Raviv M, Krasnovsky A, et al. Compost induces protection against Fusariumoxysporum in sweet basil[J]. Crop Protection, 2002, 21（7）: 583-587.

[11] 王旭輝, 徐鑫, 寶哲, 等. 高通量測序分析玉米秸稈與牛糞聯(lián)合發(fā)酵階段微生物多樣性變化[J]. 食品與發(fā)酵工業(yè), 2019, 45（3）: 47-55.

[12] Cui E, Wu Y, Zuo Y, et al. Effect of different biochars on antibiotic resistance genes and bacterial community during chicken manure composting[J]. Bioresource Technology, 2016, 203: 11-17.

[13] Zhang J, Lin H, Ma J, et al. Compost-bulking agents reduce the reservoir of antibiotics and antibiotic resistance genes in manures by modifying bacterial microbiota[J]. Science of The Total Environment, 2019, 649: 396-404.

[14] Pitta D W, Parmar N, Patel A K, et al. Bacterial diversity dynamics associated with different diets and different primer pairs in the rumen of kankrej cattle[J]. PLOS ONE, 2014, 9（11）: e111710.

[15] Zhang L, Li L, Pan X, et al. Enhanced growth and activities of the dominant functional microbiota of chicken manure composts in the presence of maize straw[J]. Frontiers in Microbiology, 2018, 9:1131.

[16] Urra J, Alkorta I, Lanzén A, et al. The application of fresh and composted horse and chicken manure affects soil quality, microbial composition and antibiotic resistance[J]. Applied Soil Ecology, 2019, 135: 73-84.

[17] Yang C W, Hsiao W C, Chang B V. Biodegradation of sulfonamide antibiotics in sludge[J]. Chemosphere, 2016, 150: 559-565.

[18] Awasthi M K, Chen H, Wang Q, et al. Succession of bacteria diversity in the poultry manure composted mixed with clay: Studies upon its dynamics and associations with physicochemical and gaseous parameters[J]. Bioresource Technology, 2018, 267: 618-625.

[19] Sekiguchi Y, Imachi H, Susilorukmi A, et al. Tepidanaerobactersyntrophicus gen. nov., sp. nov., an anaerobic, moderately thermophilic, syntrophic alcohol- and lactate-degrading bacterium isolated from thermophilic digested sludges[J]. 2006, 56（7）: 1621-1629.

[20] Slobodkin A I, Tourova T P, Kostrikina N A, et al. Tepidimicrobium ferriphilum gen. nov., sp. nov., a novel moderately thermophilic, Fe（III）-reducing bacterium of the order Clostridiales[J]. 2006, 56（2）: 369-372.

[21] Yokoyama H, Wagner I D, Wiegel J.Caldicoprobacter oshimai gen. nov., sp. nov., an anaerobic, xylanolytic, extremely thermophilic bacterium isolated from sheep faeces, and proposal of Caldicoprobacteraceae fam. nov[J]. 2010, 60（1）: 67-71.

山東省現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)體系（SDAIT-13-011-10）；山東省重大科技創(chuàng)新工程（2017GGH5129）；山東省農(nóng)業(yè)重大應(yīng)用技術(shù)創(chuàng)新項(xiàng)目；煙臺(tái)市重點(diǎn)研發(fā)計(jì)劃（2021YT06000636）

（2021–02–08）

S852.6

A

1007-1733（2021）05-0001-04