溫度對豬糞漿N2O轉(zhuǎn)化的影響

陸 揚,王 通,,成艷芬,夏 東?

(1上海市農(nóng)業(yè)科學(xué)院畜牧獸醫(yī)研究所,上海201106;2南京農(nóng)業(yè)大學(xué)消化道微生物實驗室,南京210095)

氧化亞氮(N2O)是農(nóng)業(yè)溫室氣體的重要來源,其溫室效應(yīng)是二氧化碳(CO2)的310倍,甲烷(CH4)的15倍[1]。畜牧源N2O排放量占農(nóng)業(yè)源N2O排放量的21.54%[2],主要來自畜禽糞便管理[3],其產(chǎn)生與消解涉及微生物的硝化和反硝化代謝反應(yīng)[4]。硝化反應(yīng)的第一限速酶氨單加氧酶(Ammonia monooxygenaes,AMO)的編碼基因amoA常被用于氨氧化菌(Ammonia oxidizing bacteria,AOB)的分子標(biāo)記;分別編碼不同類群反硝化菌的亞硝酸還原酶(Nitrite reductase,NIR)的nirK和nirS,以及編碼反硝化菌的氧化亞氮還原酶(Nitrous oxide reductase,NOS)的nosZ則常被用于反硝化菌的分子標(biāo)記。

溫度是影響N2O轉(zhuǎn)化的重要因素之一,但對此報道不一。有研究認(rèn)為,溫度升高會增加N2O排放[5-6]。但也有研究報道,N2O的排放與氣溫和糞污溫度均沒有關(guān)系[7]。其不一致的原因一方面可能是N2O濃度遠(yuǎn)低于其他兩種溫室氣體,檢測相對誤差高,且在現(xiàn)場測定時易受風(fēng)等自然因素的干擾;另一方面N2O是反硝化反應(yīng)的中間代謝產(chǎn)物,參與N2O轉(zhuǎn)化的各階段反應(yīng)相對速度的改變可能會直接影響N2O的排放。本試驗通過研究不同溫度條件下N2O的產(chǎn)生與消解規(guī)律及代謝底物和可能參與代謝的功能基因豐度變化,旨在探究溫度影響N2O產(chǎn)生與消解機制。

1 材料與方法

1.1 試驗材料、方法與分組

新鮮豬糞用去離子水以1∕2稀釋混勻制成糞漿,分裝于45個150 mL容積血清瓶。用橡膠塞和螺口蓋密封后,抽真空2 min,通高純空氣(78.09%N2,20.94%O2,0.93%Ar和0.033%CO2),重復(fù)抽、充氣3次后平衡氣壓。試驗按照培養(yǎng)溫度不同分為3組,每組15瓶,分別放入15℃、25℃和35℃的培養(yǎng)箱中培養(yǎng)。每組設(shè)5個采樣時間點(1 d、3 d、5 d、10 d和15 d),每個時間點取3個血清瓶進(jìn)行采樣,作為3個重復(fù)。取樣方法如下,首先用經(jīng)三通閥連接的氣密性注射器與壓力傳感器測定血清瓶中的氣體壓力和產(chǎn)氣量,然后將氣密式注射器取出的氣體保存在提前用高純氮氣反復(fù)抽洗3次以上的氣袋中保存,待測N2O濃度。取部分糞樣-20℃保存,用于測定三氮(氨態(tài)氮、亞硝態(tài)氮和硝態(tài)氮)濃度。

1.2 N2O濃度測定

氣相色譜法[俊齊儀器設(shè)備(上海)有限公司,GC-9890A]檢測N2O濃度[8]。檢測條件為Porapak Q填充柱(內(nèi)徑4 m×3 mm),電子捕獲檢測器(ECD),進(jìn)樣室100℃,柱溫50℃,檢測室320℃,載氣為99.999%氮氣,柱壓0.17 MPa,進(jìn)樣量5 mL。標(biāo)準(zhǔn)氣體購自上海偉創(chuàng)標(biāo)準(zhǔn)氣體有限公司。

1.3 三氮的提取與測定

用1 mol∕L KCl溶液浸提,固液比為1∶10,浸提糞漿中的銨態(tài)氮、亞硝態(tài)氮和硝態(tài)氮,加1%硫酸鋅溶液,0.5%氫氧化鈉溶液和10%氫氧化鋁懸浮液,混勻后,12 000 r∕min離心5 min,上清液用于3種不同形態(tài)氮的測定。納氏試劑法用于測定銨態(tài)氮;格利斯試劑用于亞硝態(tài)氮測定;雙波長分光光度法用于硝態(tài)氮測定。

1.4 DNA提取與qPCR

取0.3 g糞漿,經(jīng)PBS(120 mmol∕L,pH 8.0)洗3次后,加800 μL CTAB溶液(十六烷基三甲基溴化銨),用珠磨法機械破碎樣品(上海凈信科技有限公司),而后用酚和氯仿∕異戊醇提取其總DNA[9]。糞漿中氨氧化菌代謝酶(amoA)、反硝化菌代謝酶(nirK、nirS和nosZ)及總菌(tot)用qPCR方法定量(引物及反應(yīng)條件見表 1)。 反應(yīng)體系 20 μL:10 μL SYBR Green PCR Master Mix(Qiagen),1 μmol∕L 上、下游引物和80 ng DNA模板。所有qPCR反應(yīng)在ABI 7500 Real-time PCR儀上進(jìn)行。相應(yīng)目的片段PCR產(chǎn)物經(jīng)克隆測序驗證后,其質(zhì)粒經(jīng)梯度稀釋后用作標(biāo)準(zhǔn)曲線模板,無菌水做陰性對照模板。

1.5 數(shù)據(jù)分析

產(chǎn)氣量用氣壓進(jìn)行校正。所有試驗數(shù)據(jù)經(jīng)Excel進(jìn)行初步計算后,用SPSS 13.0軟件中的單多因子方差分析法(One Way Anova)和一般線性模型(GLM Multivariate)進(jìn)行方差統(tǒng)計,用最小顯著差異法(LSD)進(jìn)行多重比較,P＜0.05為差異顯著,P＜0.01為差異極顯著。所得數(shù)據(jù)均以“平均值±標(biāo)準(zhǔn)差(Mean±SD)”表示。用Pearson法進(jìn)行相關(guān)性分析,P＜0.05相關(guān)性顯著,P＜0.01為相關(guān)性極顯著。

表1 定量PCR各引物序列及溫度Table 1 Primer sequences and thermal cycling conditions used for qPCR

2 結(jié)果與分析

2.1 溫度對糞漿的產(chǎn)氣量及N2O轉(zhuǎn)化的影響

如圖1所示,35℃、25℃和15℃三個溫度處理組,分別在第1天、第3天和第10天開始檢測到有氣體產(chǎn)生,其后產(chǎn)氣量隨時間呈線性增加;同一時間點,以35℃組的產(chǎn)氣量最高,25℃次之,15℃最低,差異極顯著(P＜0.01)。

圖1 儲存糞漿的產(chǎn)氣量以及N2O濃度、含量和速率的變化Fig.1 Variation of emission of N2O from pig maure slurry with time

試驗開始時,糞漿上方為高純空氣,其中N2O濃度均在檢測下限,N2O含量可忽略不計。試驗開始后,3個組的N2O的濃度和含量都經(jīng)歷了先升高后降低的過程(圖1)。3組的N2O累積量都在第3天達(dá)到峰值,按照溫度由低到高,N2O累積量依次為1.70 mg∕kg、2.57 mg∕kg和4.65 mg∕kg,3組之間差異極顯著(P＜0.01);從10 d開始,15℃處理組的N2O濃度就維持在0.6 mg∕L左右,該值接近大氣中N2O的背景值,而25℃和35℃處理組N2O濃度則低于大氣N2O濃度,并能繼續(xù)下降。

在試驗前3 d速率為正值(圖1),說明N2O的產(chǎn)生速度快于還原速度,N2O的累積主要發(fā)生在第1天,以35℃組最快,25℃次之,15℃最慢,3組之間差異極顯著(P＜0.01);在試驗后期速率為負(fù)值,說明N2O的還原速度快于產(chǎn)生速度,N2O的消解主要發(fā)生在5—10 d,仍以35℃組最快,25℃次之,15℃最慢,3組之間差異極顯著(P＜0.01)。

2.2 糞漿的三氮變化

糞漿中的銨態(tài)氮和硝態(tài)氮濃度在試驗期內(nèi)持續(xù)增加,35℃處理組的銨態(tài)氮濃度極顯著高于其他兩組(P＜0.01),15℃處理組的硝態(tài)氮濃度在5 d以后顯著低于其他兩組(P＜0.05,圖2)。糞漿中亞硝態(tài)氮濃度僅為1.5—4.5 mg∕kg,各組亞硝態(tài)氮在5 d達(dá)到峰值,隨后逐漸降低,不同溫度處理組亞硝態(tài)氮沒有顯著差異(P＞0.05,圖2)。

圖2 糞漿的三氮變化情況Fig.2 Variation of concentrations of NH4+-N,N2O and NO3--N of pig manure slurry with time

2.3 糞漿中總菌及硝化、反硝化反應(yīng)功能基因的變化

糞漿中總菌及硝化、反硝化反應(yīng)功能基因的變化見表2。試驗前3 d的總菌、amoA、nirK和nirS豐度均顯著高于5 d及以后的相應(yīng)基因豐度(P＜0.01),僅nosZ豐度差異不顯著(P＞0.05)。35℃組的總菌和nirK豐度顯著低于15℃和25℃組(P＜0.05),后兩組之間差異不顯著(P＞0.05)。溫度對nosZ豐度的影響有顯著的時間效應(yīng)(P＜0.05),表現(xiàn)為3 d時,25℃組nosZ豐度顯著低于其他兩組(P＜0.05),15℃和35℃組之間差異不顯著(P＞0.05);15 d時15℃組nosZ豐度顯著高于其他兩組(P＜0.01),25℃和35℃組之間差異不顯著(P＞0.05),其他時間3組之間無顯著差異(P＞0.05)。溫度對amoA和nirS豐度的影響則不顯著(P＞0.05)。相關(guān)性分析結(jié)果表明,N2O的排放速率與總菌、amoA、nirK、nirS和nosZ豐度都顯著相關(guān),相關(guān)系數(shù)分別為0.64、0.62、0.66、0.58和0.49(P＜0.01,表2)。由于nosZ是唯一的N2O還原酶,用其他各基因豐度與nosZ豐度的差值進(jìn)行相關(guān)性分析,結(jié)果表明僅nirK∕nosZ與N2O的排放速率顯著相關(guān),相關(guān)系數(shù)為0.57(P＜0.01,未列出)。

表2 每克糞漿中總菌及硝化、反硝化反應(yīng)功能基因拷貝數(shù)的變化Table 2 Variations of total bacteria and amoA,nirK,nirS and nosZ gene copy numbers per gram manure slurry

3 討論

畜牧源N2O排放主要來自糞污后期管理,包括糞污的儲存與處理。Wang等[5]報道,沼液的存儲溫度從5℃上升到25℃,其總氮排放從8.9%上升到41.9%,25℃是N2O排放顯著增加的臨界溫度。陸日東等[6]報道,5—35℃范圍內(nèi),溫度越高牛糞N2O排放越多。本試驗開始時,糞漿上方是高純空氣,因此,前期結(jié)果與上述在開放條件下測定的結(jié)果基本一致。由于是密閉的試驗體系,隨著氧的耗竭,血清瓶內(nèi)形成嚴(yán)格厭氧環(huán)境,試驗前期累積的N2O被進(jìn)一步還原成N2,N2O濃度和含量開始下降。溫度升高不但會加速N2O累積,還會加速N2O還原,3個處理組的N2O濃度和含量都同時在第3天達(dá)到峰值,并在10 d后,N2O濃度都下降到甚至低于大氣的本底值。以上結(jié)果說明,過短的厭氧處理時間會增加N2O的排放,當(dāng)厭氧處理時間超過10 d時,N2O排放可忽略不計。提示,為了減少N2O排放,養(yǎng)殖場糞污厭氧處理時間不可少于10 d。

N2O產(chǎn)生與消解涉及微生物的硝化和反硝化代謝反應(yīng)[4]。定量PCR結(jié)果表明,在N2O減少的階段(5—15 d)糞漿中的總菌、amoA、nirK和nirS豐度均顯著低于N2O增加的階段(1—3 d),但兩個階段的nosZ豐度差異不顯著,而后者是將N2O還原的關(guān)鍵代謝酶,因此,在試驗后期N2O的含量下降可能是由于N2O的生成速度降低,還原速度超過了生成速度造成的。Levy-Booth等[12]報道,土壤中的amoA、nirK、nirS和nosZ豐度與N2O排放顯著相關(guān);Cui等[13]報道土壤中nirK和nirS與N2O排放都顯著相關(guān),但nirS相關(guān)性更強。本研究中amoA、nirK、nirS和nosZ豐度都與N2O排放顯著相關(guān),相關(guān)系數(shù)分別為0.62、0.66、0.58和0.49(P＜0.01),說明硝化和反硝化反應(yīng)的各階段都可影響糞漿N2O的排放。

Poh等[14]認(rèn)為,反硝化反應(yīng)各環(huán)節(jié)的還原速度受溫度影響不一致,是溫度升高引起N2O排放增加的原因之一,如從25℃升高到35℃時,反硝化反應(yīng)過程中的硝態(tài)氮、亞硝態(tài)氮和N2O還原速率分別上升了62%、61%和41%。本試驗并未發(fā)現(xiàn)各代謝酶基因隨溫度的升高而豐度增加,相反,35℃組的總菌和nirK豐度顯著低于15℃和25℃組(P＜0.05);35℃組的nosZ豐度也在3 d和15 d時顯著低于15℃組(P＜0.05),溫度對amoA和nirS豐度的影響則不顯著(P＞0.05)。出現(xiàn)以上結(jié)果的可能原因是在任何一個封閉的環(huán)境都存在群體效應(yīng),受營養(yǎng)物質(zhì)或代謝產(chǎn)物的限制,會出現(xiàn)衰退期,溫度的升高可能會加快菌的新陳代謝,在提高代謝效率的同時加速了菌的死亡。事實上,N2O的排放同時受其生成和還原的共同影響,不是由單一的菌或單一的酶決定的,本試驗將與N2O生成可能有關(guān)的代謝酶豐度與其唯一還原酶nosZ的豐度相減后進(jìn)行相關(guān)性分析,結(jié)果表明僅nirK∕nosZ與N2O的排放速率顯著相關(guān)(0.57,P＜0.01),說明N2O排放可能與反硝化菌中nirK與nosZ的相對數(shù)量關(guān)系更密切。

除溫度外,底物濃度也會影響N2O的排放。Huang等[15]報道,高氨氮會引起亞硝態(tài)氮的積累,增加N2O排放。Kim等[16]報道,硝酸鹽濃度為1 500 mg∕L時,污水中25%—40%的硝酸鹽轉(zhuǎn)化為N2O,當(dāng)硝酸鹽濃度降為750 mg∕L時,反硝化終產(chǎn)物中則檢測不到N2O。以往研究認(rèn)為,亞硝酸氮是引起N2O累積的重要原因,但近來研究表明,質(zhì)子化的亞硝酸氮,即游離亞硝酸(Free nitrous acid,FNA)才是污水處理過程中抑制細(xì)菌生長和增加N2O產(chǎn)生的真正原因[17]。當(dāng)FNA濃度達(dá)到0.42—1.72 mg N∕L時,AOB的活性降低50%;當(dāng) FNA濃度達(dá)到 0.011—0.070 mg N∕L時,硝化反應(yīng)開始受到抑制,當(dāng)達(dá)到 0.026—0.220 mg N∕L時,硝化反應(yīng)完全受到抑制;當(dāng)FNA濃度達(dá)到0.010—0.025 mg N∕L時,反硝化反應(yīng)下降40%,當(dāng)達(dá)到0.2 mg N∕L時,反硝化反應(yīng)完全停止。本試驗中,隨著有機氮被分解成無機氮,糞漿內(nèi)銨態(tài)氮和硝態(tài)氮都持續(xù)升高,在第15天時35℃的分別達(dá)到1 500 mg N∕L和100 mg N∕L,雖然亞硝酸氮在前5 d升高,隨后降低,似乎與N2O的累積的變化規(guī)律相似。FNA受硝酸氮濃度、溫度和pH的共同影響,通過公式[18]計算,本試驗的FNA濃度在0.08—0.69 μg N∕L,遠(yuǎn)低于以上報道閾值,說明厭氧處理糞漿時,不會因亞硝酸的累積而導(dǎo)致N2O的大量生成。

4 結(jié)論

溫度升高會加快糞漿N2O生成和還原,提高其排放峰值,但不影響其轉(zhuǎn)化的時間節(jié)律。硝化和反硝化菌可能都參與了N2O的轉(zhuǎn)化,但反硝化代謝酶nirK與nosZ可能在N2O排放中起關(guān)鍵作用。

[1]IPCC.Guidelines for national greenhouse gas inventories,Volume 4:Agriculture,forestry and other land use[M].Kanagawa,Japan:The IPCC National Greenhouse Gas Inventories Program,2006.

[2]譚秋成.中國農(nóng)業(yè)溫室氣體排放:現(xiàn)狀及挑戰(zhàn)[J].中國人口·資源與環(huán)境,2011,21(10):69-75.

[3]國家發(fā)展和改革委員會應(yīng)對氣候變化司.中華人民共和國氣候變化第二次國家信息通報[M].北京:中國經(jīng)濟出版社,2013.

[4]劉秀紅,楊慶,吳昌永,等.不同污水生物脫氮工藝中N2O釋放量及影響因素[J].環(huán)境科學(xué)學(xué)報,2006,26(12):1940-1947.

[5]WANG Y,DONG H,ZHU Z,et al.CH4,NH3,N2O and NO emissions from stored biogas digester effluent of pig manure at different temperatures[J].Agriculture,Ecosystems&Environment,2016,217:1-12.

[6]陸日東,李玉娥,萬運帆,等.溫度對奶牛糞便溫室氣體排放的影響[J].生態(tài)與農(nóng)村環(huán)境學(xué)報,2008,24(1):29-31.

[7]SOMMER S G,PETERSEN S O,SOGAARD H T.Greenhouse gas emission form stored livestock slurry[J].J Environ Qual,2000,29:744-751.

[8]王躍思,鄭循華,王明星,等.氣相色譜法檢測大氣中N2O濃度[J].分析測試技術(shù)與儀器,1994(2):19-24.

[9]ZOETENDAL E G,AKKERMANS A D,DE VOS W M.Temperature gradient gel electrophoresis analysis of 16S rRNA from human fecal samples reveals stable and host-specific communities of active bacteria[J].Appl Environ Microbiol,1998,64(10):3854-3859.

[10]CASTILLO M,MARTIN-ORUE S M,ROCA M,et al.The response of gastrointestinal microbiota to avilamycin,butyrate,and plant extracts in early-weaned pigs[J].J Anim Sci,2006,84(10):2725-2734.

[11]HALLIN S,JONES C M,SCHLOTER M,et al.Relationship between N-cycling communities and ecosystem functioning in a 50-year-old fertilization experiment[J].Isme Journal,2009,3(5):597-605.

[12]LEVY-BOOTH D J,PRESCOTT C E,GRAYSTON S J.Microbial functional genes involved in nitrogen fixation,nitrification and denitrification in forest ecosystems[J].Soil Biology and Biochemistry,2014,75:11-25.

[13]CUI P,FAN F,YIN C,et al.Long-term organic and inorganic fertilization alters temperature sensitivity of potential N2O emissions and associated microbes[J].Soil Biology and Biochemistry,2016,93:131-141.

[14]POH L S,JIANG X,ZHANG Z,et al.N2O accumulation from denitrification under different temperatures[J].Appl Microbiol Biotechnol,2015,99(21):9215-9226.

[15]HUANG T,GAO B,HU X K,et al.Ammonia-oxidation as an engine to generate nitrous oxide in an intensively managed calcareous fluvo-aquic soil[J].Sci Rep,2014,4:3950.

[16]KIM E W,BAE J H.Alkalinity requirements and the possibility of simultaneous hetemtmphic denitritifiation during sulfur utilizing autotrophic denitrification[J].Water Science and Technology,2000,42(34):233-238.

[17]ZHOU Y,OEHMEN A,LIM M,et al.The role of nitrite and free nitrous acid(FNA)in wastewater treatment plants[J].Water Res,2011,45(15):4672-4682.

[18]ANTHONISEN A C,LOEHR R C,PRAKASAM T B,et al.Inhibition of nitrification by ammonia and nitrous acid[J].J Water Pollut Control Fed,1976,48(5):835-852.