——————————————————————————————————————————————————————————————————————————

我國作為農業大國秸稈資源量豐富,但利用率較低造成了資源浪費和環境污染,利用秸稈生產沼

氣是提高秸稈利用率和緩解環境污染的有效途徑之  。新型沼氣池沼氣發酵是在多種微生物共同作用下完成的,

溫度,濃度,接種物,CN3,pH值和攪  拌“以及物料成分配比等對厭氧發酵有重要影  響。氮素是影響新型沼氣池沼氣發酵重要的營養元素之一,秸稈的碳氮比極高,這意味著以秸稈為原料的新型沼氣池沼氣發  酵體系中氮素對于發酵過程的影響可能更大,優化發酵條件有利于微生物迅速繁殖,從而提高反應器性能和新型沼氣池沼氣產量。

根據秸稈的分子式計算秸稈完全液化能夠提供厭氧發酵微生物所需的氮),但秸稈在實際發酵過程中很難全部被利用,且分解時間較長,微生物所需的氮含量不能在發酵各個階段滿足微生物的需求,尤其是在發酵初始階段,影響了新型沼氣池沼氣工程的發酵效率和最終產量。余少杰等通過添加不同濃度的氯化銨來研究外加氮源對秸稈厭氧發酵的影響,研究結果0顯示適量濃度的氮添加可以促進厭氧發酵體系中纖維素等物料的水解,而高濃度的氮素添加會抑制纖維素水解和甲烷生成。 Pia tek等通過模型分析發現氮含量是影響水解和甲烷生成過程動力學的關鍵因素。 Wagner等研究了添加不同種類的氮源在厭氧發酵過程中對甲烷產量的影響,發現  氮源中碳含量會影響甲烷的產率,在適宜的濃度范  圍內才能獲得較高的甲烷產量。雖然原料中碳素轉  化為CH4和CO2是碳素間的循環,但在厭氧發酵  中,微生物是將秸稈中的碳素轉化為CH4的功能執  行者。 Wagner等從多樣性指數上分析了不同氮  源對微生物多樣性的影響,但未對功能菌株的豐度

變化進行闡述。 Alsouleman'14-15等研究發現,高氨  氮脅迫下厭氧反應體系發生適應性變化,主導菌群  從擬桿菌轉變為梭菌。隨著宏基因組測序技術的快  速發展,應用宏基因組技術來研究樣本中直接參與  碳和氮循環的功能微生物類群的總量和多樣性是近年來環境微生物的研究熱點之一,但對添加氮素后微生物群落多樣性以及功能基因的變化以及如何影響發酵體系和甲烷產量的研究鮮有報道。因此,研究氮素添加對新型沼氣池沼氣發酵過程中微生物群落和關鍵功能基因的影響具有重要意義。

本研究解析玉米秸稈新型沼氣池沼氣發酵過程中添加氮素對料液中微生物多樣性的影響,探究產甲烷關鍵功能基因豐度的變化,通過計算纖維素等物料的降解情況并結合產氣量與產氣周期,從微生物和功能基因的角度揭示氮素添加與秸稈新型沼氣池沼氣產量之間的關系,解析發酵過程中不同階段細菌、古菌的多樣性與  功能特征,注釋功能菌株的豐度的變化,進一步探究

新型沼氣池沼氣發酵中關鍵因子的作用,為提高玉米秸稈新型沼氣池沼氣發酵產氣效率提供數據支撐。

1材料與方法

1、1試驗原料

試驗用玉米秸稈取自黑龍江八一農墾大學實驗實習基地,在10月份收獲籽粒后,秸稈風干粉碎備用,接種污泥源于本實驗室運行良好的厭氧發酵反應器,主要性質見表1,其中玉米風干秸稈的纖維素、半纖維素和木質素含量分別為30.52%,27.77%和10.99%。

1.2間歇式秸稈新型沼氣池沼氣發酵體系建立

試驗裝置選用25L自制厭氧發酵反應器,上設有出氣口和進料口,下設有出料口。出氣口與15L集氣袋相連,用于收集氣體并測定體積(見圖1)。甲烷含量通過GA2000便攜式新型沼氣池沼氣分析儀( GeotechBiogas Check)測定。發酵體積為20L,接種污泥12L,總TS設定為8%,啟動時秸稈一次性進料,不足蒸餾水補齊。

根據前期單因素試驗結果,選取添加有機氮尿素1gL作為處理,不加氮為對照,設置3個重復,控制發酵溫度為35℃,反應時間60d。隨著發酵的  進行,每10d進行1次投料,物料添加量為總秸稈添加的1/6

1.3試驗方法

1.3.116 SrrNA測序

在前期連續進料的新型沼氣池沼氣發酵試驗中,對照和處

理分別在發酵的第1.10,40天和60天進行取樣,每

個樣品平行取3次,混合均勻后加入緩沖液( Extrac

tion Buffer)-20℃保藏。利用改進的氯化芐法提取  樣品總DNA。采用1%瓊脂糖凝膠電泳檢測,使

用 Nanodrop200e(Them)檢測DNA的純度和濃

度。將合格的樣品總DNA送至四川博貝特生物科

技有限公司進行高通量測序。PCR擴增引物對應

區域:細菌16V4區引物515F(5′ GTGCCAGO

MGCCGCGGTAA3)和909R(5′ CCCCGYCAATTC  MTTTRAGT3)。古菌引物344F(5′ ACGGGGYG

CAGCAGGCGCGA3)和806R(5′ GGACTACH

VGGGTWTCTAAT-3)

1.3.2統計分析

研究采用 oniginpro2017軟件對試驗數據進行處理和分析。

1.3.3宏基因組學分析

在第一次產氣高峰第10天進行取樣,處理與對照樣品平行取3次,混合均勻后加入緩沖液( Extraction Buffer),將樣品送至上海派森諾公司進行DNA提取和宏基因組測序。

1.3.3.1物種注釋

將每個樣本的 Scaffolds/ Scanties序列與NCBlNT數據庫中的細菌、古菌、真菌和病毒序列進行  BLASTN相似性比對(E值設定為<0.001)。由于每一條目標序列可能匹配多條參考序列,而這些匹配的參考序列又分屬不同的分類單元,為使分析嚴謹可靠,同時又不丟失生物學意義,我們在 MEGAN軟件中采取“最低共同祖先( Lowest Common Ancesor,LCA)”算法),將參考序列分化為不同物種分枝前的最后一級共同分類,作為目標序列的物種分類注釋信息。結合 Scaffolds/ Scaftigs序列在各樣本中的豐度數據,獲得各樣本在所需分類等級上的相對豐度分布表

1.3.3.2功能注釋

基因的功能注釋是將非冗余蛋白序列集同常用的蛋白數據庫進行比對,從而對各樣本中的基因功能進行注釋分析。KEGG注釋是指通過將蛋白序列和KEGG代謝通路數據庫進行比對,對宏基因組預測得到的基因根據代謝通路進行注釋和分類。我們將上述非冗余蛋白序列集上傳至KAS( KEGG AutomaticAnnotation Server)進行功能注釋,并對返回的  g  注釋結果進行匯總統計,獲取各等級的注釋結果及

對應的豐度信息。

2結果與分析

2.1添加氮素對間歇式秸桿新型沼氣池沼氣發酵產氣量的影響

在發酵的前28d,處理的日產氣量穩定且高于對照(見圖2),在第10天達到第1個產氣高峰,為17750mLd-。處理和對照均在第25天達到產氣高峰,分別為15077mLd-和11707mL·d-。從第28天到第41天,對照恢復其產氣潛力,日產氣量高于處理,發揮產氣優勢。發酵前30d進料對產氣量的變化并不明顯,從第40天開始,處理和對照在產氣量上無顯著性差異,產氣趨勢相似。產氣量在每次進料后均上升,之后開始下降,到下一次進料產氣量又出現了顯著提升。測定發酵結束后纖維素和半纖維素降解率,處理達到33%和44%,顯著高于對照的6%和28%。

日產甲烷量變化趨勢與日產氣量的變化趨勢相似,筆者對發酵進行不同天數后的累積甲烷產量進行統計分析(見圖3),發酵前10d,20d,30d處理的甲烷累積產量分別是64943mL,137048mL,199459mL,顯著高于對照的14110mL,50343mL,123720mL。發酵前40d,處理的甲烷累積產量仍高于對照。發酵結束后,處理的甲烷總產量為308200mL,仍然高于對照的284728mL。處理在第10天達到第1個產甲烷高峰,為10880mL·d-1。對照在第20天達到一個峰值,為6103mL·d-1。處理和對照均在第23天達到日產甲烷產量的最高峰,分別為1388mL·d-和8164mL·d。從第26天到第41天,對照日產甲烷量高于處理。從第41天到發酵結束,處理和對照的甲烷產量無顯著差異。

發酵料液中提取的7個樣品總DNA條帶清晰,雜質較少,無嚴重拖尾現象。 Nanodrop檢測出的DNA濃度在197.8~380.9ng·L之間,OD260/280在1.90~2.07之間,均符合送樣要求(濃度>10ng·L-;OD260/280在1.8-2.2之間)。對7個樣品中的細菌和古菌的16 SrRNA V4區序列進行高通量測序。對比 Shannon指數可以看出(見表2),古菌的 Shannon指數均低于細菌,說明細菌的微生物多樣性要高于古菌的多樣性。比較處理和對照中的細菌 Alpha多樣性指數,發酵第10天處理高于對照,隨著發酵的進行,對照 Shannon指數增加,說明其多樣性有所增加,而處理的 Shannon指數和simpson指數均呈下降趨勢,微生物多樣性下降。從古菌的 Alpha多樣性指數變化規律可以看出,處理中古菌多樣性的變化規律同處理中細菌變化一致。

2.2.2發酵料液中細菌和古菌群落組成分析  將細菌的各菌群豐度占比1%以上的25類菌  屬作圖(見圖4),從屬的分類水平進行分析,發酵原  料中 Clostridium占比最高。對照中,普氏菌屬(Pre  wdla)、擬桿菌屬( Bacteroides)、瘤胃球菌屬(Rumi  m)和 Sphaerochaeta占比較高。處理中Ramnofilibacter,螺旋體屬( Treponema)和 Sphaerocho占比較高。比較發酵各時期對照和處理均表現出不同的細菌菌群組成,二者優勢菌群各不相同,菌群多樣性存在差異。比較對照和處理發現, Sphaerochaea菌屬在發酵的各時期均存在。

古菌的各菌群豐度在屬的分類水平上分析(見圖5),在整個發酵過程中,甲烷鬃毛菌屬( Methano-ea)(36.7%~68.7%)一直占據最高比例。對照中占比較高的菌群包括 Methanosaeta(45.8%5.1%), Methanoculleus(5.6%~16.2%)。處理中占比較高的菌群包括 Methanosaeta(36.7%~68.7%),

Methanoculleus(7.2%-14.2%)和甲烷螺旋菌屬Methanospirillum)(6.6%-85%)

在發酵第10天時,對照中 Methanosaeta的豐度為57.1%,處理為36.7%,優勢菌群 Methanosaeta的數量隨氮濃度的增加而減少。隨著發酵的進行,處理中的 Methanosaeta豐度逐漸增加,但對照中卻

在中期開始減少。說明 Methanosaeta在氮素充足的環境中缺少競爭力,豐度下降。只有在環境中營養成分缺少的情況下 Methanosaeta的競爭力才開始顯現,成為環境中的優勢菌屬。 Methanoculleus的豐度也受氮素的影響,其豐度隨著氮濃度的升高而增加,Methanoculleus能在氮素充足的環境中快速成為優勢菌屬,隨著發酵的持續和氮素的消耗,對照中Methanoculleus的豐度要高于處理。說明 Methaoculeus可以充分利用秸稈分解后的營養物質成為環境中優勢菌屬,在環境中有較強的競爭力和適應能力。在發酵初期,氮素處理產氣顯著高于對照,可以看出在玉米秸稈厭氧發酵過程中, Methanoculleus對產氣的貢獻更大。隨著發酵的進行,對照產氣量持續提升,較處理產氣率高,其中的甲烷桿菌屬( Methanobacterium)豐度逐漸提升,對產氣提升做出了較大貢獻。

2.3玉米秸稈新型沼氣池沼氣發酵體系宏基因組學研究  2.3.1物種豐度差異分析

前期16 SrRNA gene測序分析已經證明在發酵  10d處理和對照在豐度上存在差異,利用宏基因組  數據可以進一步分析每個樣本在各分類學水平的組  成豐度分布,逐一比較每個分類單元在兩個樣本之間的豐度差異,并通過統計檢驗評價差異是否顯著。  處理與對照在門水平上有27個門具有顯著差異,屬  水平上有475個屬具有顯著差異,在種水平上有777個種具有顯著差異(見表3)。物種豐度的變化直接影響新型沼氣池沼氣發酵體系群落的功能性,尤其是功能

菌種豐度的顯著變化,通過以上數據可知,氮素添加  能夠顯著影響發酵體系中微生物的豐度。

2.3.2功能基因差異分析

在功能基因水平上,我們分析了能量代謝中各功能基因的豐度(見圖6),發現氮素處理中甲烷代謝( methane metabolism)的豐度高于處理,而氧化磷酸化等有關能量代謝的基因的豐度變化不明顯。甲烷代謝基因是產甲烷的關鍵基因,氮素添加提升了甲烷代謝基因的豐度。

討論

Wagner"等研究了多種氮源添加對產氣高峰期的影響,與本研究的氮素添加后產氣高峰提前且一段時間產氣量顯著提升的研究結果一致,但沒有對多周期發酵進行測定。我們在研究氮素添加的半連續發酵結果顯示,處理在中后期產氣能力低于對照。氮素添加雖然促進了纖維素的水解,為古菌的代謝提供了營養底物。但隨發酵進行,處理中其它微生物所需的元素和營養底物被大量消耗,后期雖然氮素含量和對照持平,但根據元素守恒定律,其它必需元素含量會低于對照。而對照初期由于氮素含量和其他營養物質含量低,細菌代謝和纖維素水解受到限制,隨著發酵的進行,纖維素水解產物等營養物質含量增加,細菌和古菌多樣性逐漸提高,產氣潛力開始顯現。但整個發酵過程中處理的纖維素等物料的水解率顯著高于對照,提供了更多的營養物質,處理的總甲烷產量仍高于對照。

在發酵起始階段 Clostridium豐度占比較高,這與 Weimer1)等發現 Clostridium能夠將纖維素和半纖維素等大分子水解成葡萄糖的結果一致。甲烷八

更球菌屬( Methanosarcina)等產甲烷古菌的豐度提  升可能和產酸菌群一樣受氮素含量和底物濃度影  響。產氣高峰期 Methanosarcina豐度增加,對產甲  具有較大貢獻,這與 Bharathi等的研究結論一  hmm在處理和對照中豐度一直占據優  致。  ,Cam2)等發現 Methanosaeta有較強的產甲烷能  力結論與本研究結果相似。發酵體系中還注釋到  Sphaerochaeta,AbtB(21等曾在白蟻的腸道中分離出  該均屬。 Ritalaht(2等曾在紅杉河底層沉積物中分  離出 Sphaerochaeta globosa,發現其能利用多種底物  進行發酵產甲酸鹽和乙酸鹽。但在新型沼氣池沼氣發酵體系  中,鮮有關于 zerochaeta的報道,其在玉米秸稈降  解產甲烷環境中可能具有重要作用。  張彩杰2等在研究了C/N對細菌群落的影響時發現碳氮比降低,細菌多樣性提升,優勢菌群豐度增高,這與本研究發現細菌的變化規律相同。本研究還分析了在產甲烷階段具有重要作用的古菌群落多樣性,發現古菌的群落多樣性變化和細菌相似對于食氫產甲烷菌和食乙酸產甲烷菌來說,豐富的

養底物可以被利用進行代謝,此時,產甲烷功能基因的豐度要比環境中微生物的多樣性更為重要Lin(21等發現功能基因的多樣性提升不利于代謝途徑的有序性,并且可能帶來更多的功能冗余,在探究有關能量代謝的基因豐度變化時發現,甲烷代謝相關基因豐度的增加和新型沼氣池沼氣產量呈正相關,與Qiang23等研究溫度對新型沼氣池沼氣發酵體系影響的結果致。雖然本研究沒有探究功能基因的多樣性與產氣量的關系,但處理中的甲烷代謝功能基因表達顯得更為集中。這說明氮素的添加和溫度一樣對甲烷代謝的功能基因具有調節作用。

本研究發現處理中古菌的多樣性雖然提升了,  但甲烷代謝的功能基因的豐度并沒有降低,說明微  生物多樣性提升和功能基因的集中表達不存在矛  盾。 Shade20等認為微生物的多樣性只是環境微生  物中一個單純的指標,我們同樣認為多樣性的高低  不能直接說明一定的問題,復雜的代謝過程需要多  種功能微生物共同參與,又或者環境中微生物多樣性低,而重要的功能微生物豐度又占據絕對優勢,或許這正是我們想要的,但自然環境往往不是這樣的,這同樣需要在具體的環境中分析。

4結論

(1)以玉米秸稈為原料的厭氧新型沼氣池沼氣發酵過程

中,利用16 SrRNA gene和宏基因組測序技術分析發現,氮素添加能夠提升厭氧發酵過程中細菌的多樣性和豐度,細菌多樣性的變化特別是具備纖維素降解功能菌株豐度的提升能夠加速纖維素、半纖維素等物質的水解和代謝過程。氮素添加和纖維素等物質的加速分解有效的提升了反應器內古菌代謝所需  的營養底物,提升了古菌的多樣性和豐度,促進了產甲烷古菌的代謝和甲烷代謝基因豐度的提升。  (2)氮素添加通過提升微生物的豐度以及產甲  烷功能基因的豐度,縮短了達到產氣高峰所需的周  期,并在一段時間內維持穩定的產氣量和甲烷含量,  最終在整個發酵周期內顯著提升新型沼氣池沼氣產氣量和TS甲烷產率。

摘自《中國沼氣》第5期 趙一全 馬茹霞 李家威 晏磊 羅濤 梅自力 王偉東