李 清，林光輝

(1.清華大學(xué) 地球系統(tǒng)科學(xué)中心 地球系統(tǒng)數(shù)值模擬教育部重點(diǎn)實(shí)驗(yàn)室，北京 100084;2.清華大學(xué) 深圳研究生院，海洋材料與技術(shù)學(xué)部，深圳 518055)

在經(jīng)歷了約十年的平穩(wěn)期后，2007年開(kāi)始全球大氣甲烷濃度又出現(xiàn)增長(zhǎng)，而且增長(zhǎng)主要來(lái)自北半球[1]。Kai[2]通過(guò)大氣甲烷穩(wěn)定同位素監(jiān)測(cè)、大氣反演、地球系統(tǒng)模型預(yù)測(cè)等多種研究手段綜合分析認(rèn)為，全球甲烷過(guò)去十年排放量增速出現(xiàn)下降的主要原因是北半球微生物資源減少，其中，亞洲特別是中國(guó)水稻田甲烷排放量的下降可能是主導(dǎo)因素[2]。Bousquet[3]研究發(fā)現(xiàn)，濕地甲烷排放量的變化是過(guò)去20年全球甲烷排放出現(xiàn)年際間波動(dòng)最重要的原因。1999年以來(lái)，全球濕地甲烷排放量的下降抵消了人類(lèi)活動(dòng)引起的大氣甲烷濃度上升，如果濕地甲烷的排放量恢復(fù)到1990年的高排放量水平，全球大氣甲烷的濃度將在經(jīng)歷這短暫的穩(wěn)定之后重新開(kāi)始增長(zhǎng)[3]。古氣候?qū)W家的研究[4]也發(fā)現(xiàn)，最后一次冰期的消失更可能是由于北半球濕地的甲烷排放引起而非甲烷冰的融化。這些研究均說(shuō)明了全球濕地特別是北半球濕地甲烷排放對(duì)于調(diào)節(jié)大氣甲烷濃度變化具有十分重要的作用。

濕地僅占陸地面積6%～7%，但其碳含量卻占陸地土壤碳儲(chǔ)量的68%[5]，凈初級(jí)生產(chǎn)力占全球的10%。濕地由于較低的有機(jī)質(zhì)分解率和較高的生產(chǎn)力而成為重要的碳匯，但同時(shí)也是重要的甲烷排放源。根據(jù)政府間氣候變化專(zhuān)門(mén)委員會(huì)(Intergovernmental Panel on Climate Change,IPCC)2007年的評(píng)估結(jié)果，濕地甲烷排放量占全球甲烷排放量的20%～39%。對(duì)于全球濕地甲烷的排放，目前不同的估算結(jié)果相差較大[6]。例如，Cao[7]利用過(guò)程模型研究發(fā)現(xiàn)，全球濕地的甲烷排放量為145 Tg·a-1，其中，天然濕地的排放量為92 Tg·a-1，水稻田的排放為53 Tg·a-1[7]。然而，Walter[8]對(duì)天然濕地的甲烷排放計(jì)算的結(jié)果則達(dá)到260 Tg·a-1，其中約60%釋放到大氣中，其他部分在排放到大氣之前就已經(jīng)在土壤中被氧化了[8]。而B(niǎo)astviken[9]研究發(fā)現(xiàn)，盡管未被計(jì)算在全球溫室氣體通量中，陸地湖泊、水庫(kù)和河流的總甲烷排放量事實(shí)上達(dá)到了103 Tg·a-1，大約抵消了目前陸地碳匯總量的25%[9]。

這種差異產(chǎn)生的重要原因是不同類(lèi)型和不同區(qū)域濕地甲烷排放的空間格局以及隨時(shí)間的動(dòng)態(tài)變化信息比較缺乏。而要得到這些結(jié)論，除了需要建立更多的甲烷通量數(shù)據(jù)監(jiān)測(cè)網(wǎng)絡(luò)外，還需要對(duì)濕地甲烷產(chǎn)生的微生物過(guò)程進(jìn)行深入的研究，從而為甲烷過(guò)程建模和通量數(shù)據(jù)分析提供支撐。

盡管?chē)?guó)際上關(guān)于產(chǎn)甲烷相關(guān)微生物的研究非常多，以產(chǎn)甲烷菌(Methanogens)或甲烷氧化菌(Methanotrophs) 作為主題詞搜索web of knowledge，查找到的文獻(xiàn)超過(guò)10 000篇，但大多集中在微觀(guān)生物學(xué)領(lǐng)域。與全球變化以及溫室氣體產(chǎn)生機(jī)制等主題相關(guān)的研究主要包括甲烷相關(guān)微生物的種群多樣性和群落結(jié)構(gòu)隨環(huán)境因子變化響應(yīng)情況[10]，包括對(duì)水位變化的響應(yīng)[11]，對(duì)溫度變化的響應(yīng)[12]，對(duì)營(yíng)養(yǎng)元素濃度變化的響應(yīng)[13]，以及甲烷相關(guān)微生物在土壤中的垂直分布情況[14]，產(chǎn)甲烷細(xì)菌的數(shù)量隨著季節(jié)變化以及土壤深度而變化[15]等。對(duì)于甲烷通量的空間格局與甲烷相關(guān)微生物的群落結(jié)構(gòu)之間的直接關(guān)系研究還比較缺乏。而事實(shí)上，如果將產(chǎn)甲烷菌的數(shù)量、種群變化規(guī)律與甲烷排放量的變化直接相關(guān)聯(lián)將更具有科學(xué)與現(xiàn)實(shí)意義[16]。

本文在簡(jiǎn)要介紹濕地甲烷產(chǎn)生的基本微生物過(guò)程后，重點(diǎn)闡明如何利用穩(wěn)定同位素技術(shù)研究產(chǎn)甲烷的不同途徑和不同微生物種類(lèi)在甲烷產(chǎn)生過(guò)程的作用，并分析該領(lǐng)域的最新研究進(jìn)展和未來(lái)研究重點(diǎn)，以期為我國(guó)開(kāi)展這方面的研究提供一些信息。

1 甲烷產(chǎn)生的基本微生物過(guò)程及其對(duì)甲烷排放計(jì)算的意義

目前的研究表明，濕地排放的甲烷主要是由產(chǎn)甲烷微生物在厭氧條件(Eh值小于-200 mV)下把乙酸、CO2或甲醇等作為電子受體，乙酸或H2作為電子供體，發(fā)生氧化還原反應(yīng)而產(chǎn)生的，有機(jī)物分解產(chǎn)生甲烷的過(guò)程示于圖1。

其中，乙酸和H2是產(chǎn)甲烷微生物所利用的主要基質(zhì)，其他的基質(zhì)對(duì)于甲烷產(chǎn)生的貢獻(xiàn)一般不超過(guò)5%[18]。根據(jù)使用基質(zhì)的不同，產(chǎn)甲烷微生物可以分為CO2型、甲基型和乙酸型三大類(lèi)，產(chǎn)甲烷微生物的三種類(lèi)型列于表1。

圖1 有機(jī)物分解產(chǎn)生甲烷的過(guò)程[17] Fig.1 Methane production pathways during organic matter decomposition[17]

類(lèi)型所利用的基質(zhì)代表性微生物CO2型CO2、CO和甲酸鹽甲酸甲烷桿菌(Methanobacterium)產(chǎn)甲烷球菌(Methanococcus)甲基型甲醇、甲胺、二甲胺、三甲胺甲硫醇和二硫醇甲烷八疊球菌(Methanosarcins)甲烷八疊球菌(Methanosarcina)乙酸型醋酸鹽甲烷絲菌(Methanothrix)鬃毛甲烷菌(Methanosaeta)

由于產(chǎn)甲烷微生物不能直接利用復(fù)雜化合物，因而其生存必須依賴(lài)其他至少三類(lèi)微生物的存在，包括水解、發(fā)酵以及產(chǎn)氫產(chǎn)乙酸細(xì)菌。

以上分析可見(jiàn)，甲烷的產(chǎn)生主要由以下三個(gè)反應(yīng)決定的：

CH3COOH→CO2+CH4

(1)

CO2+4H2→2H2O+ CH4

(2)

4 CH3OH→H2O+ CO2+3CH4

(3)

因而甲烷產(chǎn)生總量m可以表示為：

m=ma+mc+mi

(4)

或者

m=(fma+fmc+fmi)m，同時(shí)fma+fmc+fmi=1

(5)

其中，m代表甲烷產(chǎn)生的總量或者速率，ma代表來(lái)自于醋酸的甲烷，mc代表來(lái)自于CO2的甲烷，mi代表來(lái)自于其他化合物比如甲醇的甲烷；fma則表示來(lái)自于醋酸的甲烷的占總甲烷的比例，fmc表示來(lái)自于CO2的甲烷占總甲烷的比例，fmi則表示來(lái)自于甲基化合物等其他途徑的甲烷占總甲烷的比例。

由于反應(yīng)(3)的發(fā)生主要是在高鹽的環(huán)境中，且甲烷的產(chǎn)生量很低[19]，因而在大多數(shù)情況下甲烷的產(chǎn)生來(lái)自反應(yīng)(1)和反應(yīng)(2) ，因而上述方程可以進(jìn)一步簡(jiǎn)化為：

ma=(1-fmc)m

(6)

已有研究發(fā)現(xiàn)，fmc在不同的環(huán)境條件下不一樣，而且受到溫度、植被種類(lèi)和富營(yíng)養(yǎng)化程度等多個(gè)因素的影響[20]，因而fmc是計(jì)算甲烷產(chǎn)生的各種過(guò)程模型中十分重要的參數(shù)。

2 穩(wěn)定同位素技術(shù)在產(chǎn)甲烷途徑貢獻(xiàn)分析中的應(yīng)用

目前，對(duì)于fmc的計(jì)算主要有3種方法[19]：(1) 放射性同位素示蹤，例如用14C標(biāo)記區(qū)分碳酸氫根和醋酸中的碳然后測(cè)定產(chǎn)生的甲烷中的14C豐度進(jìn)而計(jì)算碳酸氫根和醋酸各自的轉(zhuǎn)化速率；(2) 利用特異性的抑制劑比如CH3F抑制利用醋酸產(chǎn)生甲烷的反應(yīng)，此時(shí)產(chǎn)生的甲烷則認(rèn)為是來(lái)自于CO2；(3)利用穩(wěn)定同位素的分餾效應(yīng)，通過(guò)測(cè)定產(chǎn)物和底物的同位素比值來(lái)計(jì)算不同產(chǎn)甲烷途徑對(duì)于總甲烷的貢獻(xiàn)率。其中，方法(3)由于操作簡(jiǎn)單，受實(shí)驗(yàn)條件以及環(huán)境的影響小且更加精確而成為目前計(jì)算fmc最重要的方法之一[19]。

不同的產(chǎn)甲烷途徑中的同位素分餾過(guò)程如圖2所示，因而在不考慮反應(yīng)(3)的發(fā)生時(shí)，可以得到以下質(zhì)量平衡方程：

δCH4=fmcδmc+ (1-fmc)δma

(7)

解上述方程可得到:

fmc= (δCH4-δma) / (δmc-δma)

(8)

圖2 甲烷產(chǎn)生過(guò)程中的同位素分餾[19] Fig.2 Isotopic fractionation in methane producing processes[19]

如果甲烷產(chǎn)生反應(yīng)的同位素效應(yīng)是確定的，則：

δma=δac+Δma

(9)

δmc=δCO2+ Δmc

(10)

而對(duì)于反應(yīng)A→B，通常需要定義一個(gè)同位素分餾因子αAB:

αAB= (δA + 103) / (δB + 103)

(11)

所以αma和αmc可以分別表示為：αma= (δac+ 103) / (δma+ 103)

(12)

αmc= (δCO2+ 103) / (δmc+ 103)

(13)

即δma和δmc可以重新表示為：

δma=(1/αma)(δac+ 103-αma103)

(14)

δmc=(1/αmc)(δCO2+ 103-αmc103)

(15)

由以上分析可見(jiàn)，不考慮反應(yīng)(3)的情況下，根據(jù)方程(8)、 (14)和(15)，如果已知αma和αmc的值，則fmc可以很容易通過(guò)測(cè)定CH4、CO2和乙酸甲基中的δ13C來(lái)確定。利用穩(wěn)定同位素信號(hào)進(jìn)行產(chǎn)甲烷途徑貢獻(xiàn)分析的研究總結(jié)列于表2，由表2可見(jiàn)，目前關(guān)于δac的信息十分缺乏，僅在Cape Lookout Bight和日本的水稻田土壤中有所研究，還有一些研究測(cè)定了總乙酸的δ13C，由于乙酸羧基和乙酸甲基的δ13C相差超過(guò)14‰[19]，總乙酸的δ13C對(duì)于fmc的準(zhǔn)確計(jì)算意義不大。

表2 利用穩(wěn)定同位素信號(hào)進(jìn)行產(chǎn)甲烷途徑貢獻(xiàn)分析[19]

另一方面，αma和αmc在研究中更難被確定。由于很多中間值不容易得到，通常用簡(jiǎn)單的參數(shù)αc來(lái)對(duì)不同產(chǎn)甲烷途徑的貢獻(xiàn)進(jìn)行區(qū)分。定義αc= (δCO2+ 103) / ( δCH4+ 103)，通過(guò)對(duì)αc進(jìn)行計(jì)算得出：當(dāng)αc>1.065時(shí)，CO2途徑占主導(dǎo)；而當(dāng)αc<1.055時(shí)，乙酸途徑占主導(dǎo)[21]。然而，這種方法只能進(jìn)行粗略的定性區(qū)分，無(wú)法得到fmc的精確值。

為了精確計(jì)算fmc，可以利用特定培養(yǎng)基分別培養(yǎng)利用不同途徑產(chǎn)甲烷的微生物，希望對(duì)αma和αmc進(jìn)行精確計(jì)算[19]。然而，αma和αmc受微生物種類(lèi)、溫度以及可利用的H2濃度等因素影響變化很大，其中αmc在1.031～1.077，αma在1.007～1.027[19]。盡管范圍變化很大，但是能夠很清楚的看出αmc要顯著的大于αma。這些結(jié)果部分驗(yàn)證了前文中對(duì)αc的計(jì)算得出的結(jié)論，同時(shí)也指出αma和αmc在各種特定類(lèi)型生態(tài)系統(tǒng)中的確定還需要更多深入研究。

3 同位素標(biāo)記技術(shù)在產(chǎn)甲烷微生物功能研究中的應(yīng)用

由于αma和αmc受微生物種類(lèi)的影響很大，因而對(duì)不同類(lèi)型濕地中的微生物群落結(jié)構(gòu)進(jìn)行分析，特別是對(duì)產(chǎn)甲烷微生物種類(lèi)、種群大小以及其功能情況進(jìn)行研究將對(duì)αma和αmc的計(jì)算提供重要信息。利用分子生物學(xué)技術(shù)能夠?qū)ν寥牢⑸锏姆N群數(shù)量、種群大小甚至功能基因的表達(dá)豐度等進(jìn)行定性或半定量分析，比如Godin[22]利用末端限制性片段長(zhǎng)度多樣性(Terminal Restriction Fragment Length Polymorphism ，T-RFLP)分析方法發(fā)現(xiàn),泥炭地甲烷產(chǎn)生速率與產(chǎn)甲烷菌群落豐富度之間在春季沒(méi)有相關(guān)性，但是生長(zhǎng)季節(jié)存在明顯相關(guān)關(guān)系，但這種分析更多的是起到指示和推測(cè)微生物在環(huán)境中的功能[22]，難以將這些參數(shù)與微生物具體代謝反應(yīng)過(guò)程間建立直接聯(lián)系。同位素標(biāo)記技術(shù)與分子生物學(xué)技術(shù)的結(jié)合為建立微生物的群落結(jié)構(gòu)、功能結(jié)構(gòu)與代謝功能之間的關(guān)系研究提供了可能[22]。

John[24]利用放射性同位素14C標(biāo)記法和T-RFLP法分別研究了鹽沼不同土壤深度的產(chǎn)甲烷微生物對(duì)不同類(lèi)型底物的利用速率和產(chǎn)甲烷微生物的種群變化，表層土壤甲胺的利用速率最快，與此同時(shí)，表層土壤微生物的分布則以利用甲胺的甲烷八疊球菌(Methanosarcinales)的種群數(shù)量占據(jù)主導(dǎo)。深層土壤則是乙酸的利用速率最快，而同時(shí)主要利用乙酸的鬃毛甲烷菌(Methanosaeta)的種群數(shù)量占據(jù)主導(dǎo)地位[23]。

目前同位素標(biāo)記中應(yīng)用最為廣泛的是穩(wěn)定同位素探測(cè)技術(shù)(Stable Isotope Probing，SIP)，是指一類(lèi)利用穩(wěn)定同位素(比如13C，15N)標(biāo)記微生物所利用的基質(zhì)，然后對(duì)微生物經(jīng)過(guò)代謝反應(yīng)后被標(biāo)記DNA或RNA等進(jìn)行分析的技術(shù)，根據(jù)檢測(cè)的分子不同進(jìn)而區(qū)分為DNA-SIP和RNA-SIP[24]。Lu等[26]在水稻田中加入13C標(biāo)記的CO2后研究發(fā)現(xiàn)，水稻根系土壤中一類(lèi)被稱(chēng)為Rice ClusterⅠ的產(chǎn)甲烷微生物RNA中13C的比例顯著升高，說(shuō)明其利用了水稻根呼吸產(chǎn)生的CO2，因而在水稻田甲烷產(chǎn)生中可能起著重要作用，盡管與其他幾種產(chǎn)甲烷微生物相比，其在水稻土壤中的相對(duì)數(shù)量并不是最高的[25]。這個(gè)結(jié)果顯示了同位素標(biāo)記技術(shù)相對(duì)分子生物學(xué)技術(shù)的優(yōu)勢(shì)，如果利用T-RFLP法對(duì)水稻根系周?chē)漠a(chǎn)甲烷菌的種群數(shù)量進(jìn)行研究，Rice ClusterⅠ會(huì)因?yàn)榉N群數(shù)量偏低而被認(rèn)為其對(duì)水稻田甲烷產(chǎn)生的貢獻(xiàn)也偏低，如果將不同產(chǎn)甲烷途徑的貢獻(xiàn)率和同位素標(biāo)記的結(jié)果結(jié)合起來(lái)，水稻田產(chǎn)甲烷微生物的種群數(shù)量和實(shí)際代謝功能與甲烷通量之間的關(guān)系將更加明確。Conrad 等[27]通過(guò)測(cè)定CO2和CH4的穩(wěn)定同位素比值來(lái)粗略估算CO2途徑和乙酸發(fā)酵途徑的相對(duì)產(chǎn)甲烷貢獻(xiàn)率并結(jié)合T-RFLP技術(shù)發(fā)現(xiàn)，溫度高于40 ℃時(shí)，土壤產(chǎn)甲烷以完全CO2還原途徑為主，而溫度低于30 ℃時(shí)，CO2還原途徑和乙酸還原途徑同時(shí)存在[26]。Conrad等的研究很好的結(jié)合穩(wěn)定同位素分餾效應(yīng)和分子生物學(xué)技術(shù)，但沒(méi)有測(cè)定土壤乙酸中甲基的碳同位素比值，因而只能對(duì)不同產(chǎn)甲烷途徑的貢獻(xiàn)率進(jìn)行粗略計(jì)算。

4 展 望

全球甲烷排放在2007年前的十年中增長(zhǎng)速率下降，這引發(fā)了科學(xué)家們激烈的爭(zhēng)論，但參與討論的主要是大氣科學(xué)家，盡管濕地和微生物的作用可能是決定性因素[2]，而濕地生態(tài)學(xué)家和微生物學(xué)家參與程度卻并不高。微生物種群的復(fù)雜性、微生物種群大小與其實(shí)際代謝功能之間的非線(xiàn)性關(guān)系、研究手段的缺乏可能是主要原因。由于研究手段的限制，利用穩(wěn)定同位素分餾效應(yīng)對(duì)不同類(lèi)型的濕地中各種產(chǎn)甲烷途徑的比重進(jìn)行精確分析還需要付出巨大的努力，尤其是αma和αmc受到的影響因素十分復(fù)雜，因而對(duì)于不同類(lèi)型濕地中αma和αmc的確定將是甲烷產(chǎn)生相關(guān)微生物過(guò)程的研究重點(diǎn)。與此同時(shí)，盡管基因組學(xué)研究的發(fā)展和應(yīng)用使微生物種群多樣性有了更為深入的研究，但許多微生物的代謝功能仍然很難確定[22]，即便是已經(jīng)進(jìn)行了全基因組測(cè)序基因組的模式生物大腸桿菌(E. coli)基因組中仍然有大約50%的基因其功能是未知的[27]。而同位素標(biāo)記技術(shù)特別是最近發(fā)展的SIP技術(shù)能夠通過(guò)標(biāo)記微生物利用的底物來(lái)對(duì)微生物在環(huán)境中的功能進(jìn)行研究。因此，如果將分子生物學(xué)技術(shù)、SIP和穩(wěn)定同位素分餾效應(yīng)等結(jié)合起來(lái)，可以更為深入研究濕地甲烷產(chǎn)生的微生物過(guò)程以及其與環(huán)境因子之間的關(guān)系，為減小濕地甲烷排放量估算的不確定性提供機(jī)理上的信息。

參考文獻(xiàn):

[1] Rigby M, Prinn R, Frase P，et al. Renewed gro-wth of atmospheric methane[J]. Geophysical Research Letters, 2008， 35:L22805.

[2] Kai FM, Tyler SC, Randerson JT, et al. Reduced methane growth rate explained by decreased Northern Hemisphere microbial sources[J]. Nature, 2011， 476:194-197.

[3] Bousquet P, Ciais P, Miller JB, et al. Contribution of anthropogenic and natural sources to atmospheric methane variability[J]. Nature, 2006， 443:439-443.

[4] Fischer H, Behrens M, Bock M, et al. Changing boreal methane sources and constant biomass burning during the last termination[J]. Nature, 2008， 452:864-867.

[5] Schlesinger WH. Biogeochemistry: an analysis of global chang[M]. Academic press, 1997.

[6] Solomon S. Climate change 2007: the physical science basis: contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change[M]. Cambridge Univ Pr, 2007.

[7] Cao M, Gregson K, Marshall S. Global methane emission from wetlands and its sensitivity to climate change[J]. Atmospheric Environment, 1998， 32:3 293-3 299.

[8] Walter B, Heimann M, Matthews E. Modeling modern methane emissions from natural wetlands. I- Model description and results[J]. Journal of Geophysical Research D Atmos, 2001， 106:34.

[9] Bastviken D, Tranvik LJ, Downing JA, et al. Freshwater methane emissions offset the continental carbon sink[J]. Science, 2011， 331:50.

[10] Yavitt JB, Yashiro E, Cadillo-Quiroz H, et al. Methanogen diversity and community composition in peatlands of the central to northern Appalachian Mountain region[J].North America Biogeochemistry, 2012， 109:117-131.

[11] Yrj?l? K, Tuomivirta T, Juottonen H, et al. CH4production and oxidation processes in a boreal fen ecosystem after long‐term water table drawdown[J]. Global change biology, 2011， 17:1 311-1 320.

[12] Kim SY, Freeman C, Fenner N, et al. Functional and structural responses of bacterial and methanogen communities to 3-year warming incubation in different depths of peat mire[J]. Applied Soil Ecology, 2012, 57:23-30.

[13] Eriksson T, ?quist MG, Nilsson MB. Production and oxidation of methane in a boreal mire after a decade of increased temperature and nitrogen and sulfur deposition[J]. Global change biology, 2010, 16:2 130-2 144.

[14] Cadillo‐Quiroz H, Br?uer S, Yashiro E, et al. Vertical profiles of methanogenesis and methanogens in two contrasting acidic peatlands in central New York State, USA[J]. Environmental Microbiology, 2006, 8:1 428-1 440.

[15] 湯博, 唐杰, 吳俐莎, 等. 若爾蓋高原產(chǎn)甲烷菌數(shù)量的時(shí)空差異性[J]. 微生物學(xué)通報(bào), 2010, 37:1 706-1 711.

Tang Bo, Tang Jie, Wu Lisha, et al. Time-space change of methanogen amount from zoige plateau[J]. Microbiology China, 2010, 37: 1 706-1 711(in Chinese).

[16] 金城. 高原濕地產(chǎn)甲烷菌[J]. 微生物學(xué)通報(bào), 2010， 37:1 705.

Jin Cheng. Methanogensin zoige wetland[J]. Microbiology China, 2010, 37:1 705(in Chinese).

[17] Reddy KR， DeLaune RD. Biogeochemistry of wetlands: science and applications[M]. CRC, 2008.

[18] Segers R. Methane production and methane consumption: a review of processes underlying wetland methane fluxes[J]. Biogeochemistry, 1998， 41:23-51.

[19] Conrad R. Quantification of methanogenic pathways using stable carbon isotopic signatures: a review and a proposal[J]. Organic Geochemistry, 2005， 36:739-752.

[20] Conrad R. Contribution of hydrogen to methane production and control of hydrogen concentrations in methanogenic soils and sediments[J]. Fems Microbiology Ecology, 1999， 28:193-202.

[21] Hornibrook ERC, Longstaffe FJ, Fyfe WS. Evolution of stable carbon isotope compositions for methane and carbon dioxide in freshwater wetlands and other anaerobic environments[J]. Geochimica Et Cosmochimica Acta, 2000， 64:1 013-1 027.

[22] Godin A, McLaughlin JW, Webster K， et al. Methane and methanogen community dynamics across a boreal peatland nutrient gradient[J]. Soil Biology and Biochemistry, 2012, 48: 96-105.

[23] Chen Y, Murrell JC. When metagenomics meets stable-isotope probing: progress and perspectives[J]. Trends in microbiology, 2010, 18:157-163.

[24] John PR, Brock F, Banning N, et al. Changes in methanogenic substrate utilization and communities with depth in a salt-marsh, creek sediment in southern England. Estuarine[J]. Coastal and Shelf Science, 2012, 96: 170-178.

[25] Dumont MG, Murrell JC. Stable isotope probing - linking microbial identity to function[J]. Nature Reviews Microbiology, 2005， 3:499-504.

[26] Lu Y， Conrad R. In situ stable isotope probing of methanogenic archaea in the rice rhizosphere[J]. Science, 2005， 309:1 088-1 090.

[27] Conrad R, Klose M, Noll M. Functional and structural response of the methanogenic microbial community in rice field soil to temperature change[J]. Environmental Microbiology, 2009， 11:1 844-1 853.

[28] Serres MH, Gopal S, Nahum LA,et al. A functional update of the Escherichia coli K-12 genome[J]. Genome Biol, 2001， 2:1-35.