微生物組學(xué)及其在厭氧消化中的研究進展

李葉青 景張牧 江皓 徐泉 周紅軍 馮璐

（1. 中國石油大學(xué)（北京）新能源與材料學(xué)院 生物燃氣高值利用北京市重點實驗室，北京 102249；2. 丹麥奧胡斯大學(xué)工程系，切勒 8830）

目前，全國每年產(chǎn)生數(shù)以億計的有機廢棄物。現(xiàn)階段常用的有機廢棄物處理方法包括：熱處理（如氣化、熱解等）以及生物處理（如厭氧消化、堆肥等）［1］。我國城市有機廢棄物的含水率高、熱值低（熱值為3 MJ/kg，含水率約80%［2］），導(dǎo)致熱處理的效率不高。而堆肥由于產(chǎn)生大量溫室氣體（CO2、CH4和N2O）、占地面積大以及周期較長等缺點，限制了其直接應(yīng)用［3］。不過，好氧堆肥可以作為厭氧消化的后續(xù)處理使沼渣深度腐熟，轉(zhuǎn)化為腐殖質(zhì)改良土壤［4］。將厭氧消化和好氧堆肥聯(lián)合使用，既解決了厭氧消化后產(chǎn)生的沼渣不宜直接還田的問題，還可縮短好氧堆肥周期、提高有機肥肥力。此外，從能源結(jié)構(gòu)、碳減排、城鎮(zhèn)化發(fā)展等多個角度考慮，沼氣工程在中國的發(fā)展已經(jīng)成為了必然趨勢［5］。到2020年，中國沼氣年利用量已達到440×1010m2［5］。盡管厭氧消化技術(shù)在我國發(fā)展迅速，但其系統(tǒng)中種群結(jié)構(gòu)等信息尚未被完全了解［6］，這在一定程度上阻礙了它們在生物能源生產(chǎn)方面的廣泛應(yīng)用［7］。如果僅依賴于微生物培養(yǎng)的技術(shù)，微生物活性和功能的生物因素（如種間關(guān)系和生物多樣性）等對消化過程產(chǎn)生的影響會被忽視。到目前為止，我們?nèi)栽谂捬跸到y(tǒng)中的群落結(jié)構(gòu)和功能與環(huán)境等因素聯(lián)系起來，并期望通過改變微生物群落來優(yōu)化過程效率和提高穩(wěn)定性［8］。微生物組學(xué)分析能夠極大的提高我們對這些復(fù)雜系統(tǒng)的理解，被認為是了解厭氧消化系統(tǒng)中的群落結(jié)構(gòu)和構(gòu)建其與外部環(huán)境變化關(guān)聯(lián)性的有力工具。本文系統(tǒng)回顧了厭氧消化過程的微生物學(xué)研究，并強調(diào)各種生物因素對厭氧消化表現(xiàn)的重要影響。

1 厭氧消化的過程

目前，對于厭氧消化過程的認識普遍基于“三階段、四種群”的基本原理［10］?！叭A段”依次為水解階段、產(chǎn)酸階段和產(chǎn)甲烷階段?！八姆N群”分別是水解酸化菌、產(chǎn)氫產(chǎn)乙酸菌、同型產(chǎn)乙酸菌和產(chǎn)甲烷菌。水解階段一般認為是復(fù)雜的有機化合物轉(zhuǎn)化為簡單的有機化合物的過程，簡單的有機化合物分為兩類：產(chǎn)甲烷物質(zhì)和其他消化產(chǎn)物。產(chǎn)甲烷物質(zhì)是可以直接被產(chǎn)甲烷菌利用并產(chǎn)生甲烷的底物，主要為乙酸（70%）［11］。其他消化產(chǎn)物則能夠被產(chǎn)氫產(chǎn)乙酸菌代謝，轉(zhuǎn)化為乙酸、氫氣和二氧化碳。與此同時，同型產(chǎn)乙酸菌將氫氣和二氧化碳轉(zhuǎn)化為乙酸（此階段可逆）。最后，在產(chǎn)甲烷階段，不同營養(yǎng)型的產(chǎn)甲烷菌將產(chǎn)甲烷物質(zhì)和H2轉(zhuǎn)化為氣態(tài)甲烷，這是整個厭氧消化的限速階段［12］。多種微生物發(fā)揮不同的功能，共同組成了厭氧消化系統(tǒng)［13］。厭氧消化的基本原理如圖1所示。

圖1 厭氧消化的基本原理（三階段四種群）

2 微生物組學(xué)類型

微生物組學(xué)包括基因組學(xué)、蛋白組學(xué)以及代謝組學(xué)，研究微生物從遺傳基因到表達代謝的全過程（表1）。基因組學(xué)包括16S rRNA基因組學(xué)［14］、宏基因組學(xué)［15］、宏轉(zhuǎn)錄組學(xué)［16］。其檢測手段為高通量測序技術(shù)，又稱為下一代測序（Next-Generation Sequencing，NGS）。高通量測序技術(shù)至今為止已經(jīng)經(jīng)歷了三代的發(fā)展，第1代高通量測序技術(shù)以Sanger測序為代表，第2代以Illumina測序平臺為代表，目前正在開發(fā)第3代——單分子測序技術(shù)［17］。近年來隨著測序平臺吞吐量的提升、參考數(shù)據(jù)庫的完善以及成本不斷下降［18］，高通量測序技術(shù)已經(jīng)被廣泛用于調(diào)查和研究各種厭氧消化系統(tǒng)的生物多樣性［19］，成為研究各種環(huán)境中微生物群落的有力工具。目前最廣泛使用的高通量測序系統(tǒng)主要基于二代 Illumina 平臺所開發(fā)的MiSeq、HiSeq和NovaSeq等［17］。基因組學(xué)的研究為微生物群落提供了詳盡的物種信息［20］。

總體而言，各種基因測序方法都有其特定的局限性，而利用多種測序方法共同分析可以彌補其不足。以16S rRNA基因測序為例，由于通用引物無法匹配所有16S rRNA靶點，研究的高變區(qū)數(shù)量有限（大多為V4）以及擴增過程中產(chǎn)生嵌合體等原因，其精度會受到制約，還會使得特定物種的豐度發(fā)生改變，例如Euryarchaeota和Spirochaetes［21］。而宏基因組測序則不受系統(tǒng)發(fā)育標(biāo)記基因的限制，能夠檢測出多種物種（細菌、真菌、古菌、病毒等），并可能獲得新基因乃至新物種信息。因此被認為是評估分類學(xué)組成最準確的方法［22］。而宏基因組的局限性在于其不能區(qū)分表達和非表達的基因，因此不能反映實際的代謝活動［23］。而宏轉(zhuǎn)錄組可被用于識別復(fù)雜生態(tài)系統(tǒng)中mRNA表達的生物學(xué)特征［16］。但宏轉(zhuǎn)錄組的研究也存在包括從環(huán)境樣本中難以回收高質(zhì)量mRNA，物種的mRNA半衰期短，以及需要從其他RNA物種中分離mRNA等問題［22］。因此，從宏基因組分析推斷出的基因豐度高并不一定表明具有高的轉(zhuǎn)錄或代謝活性。Maus等［24］通過對嗜熱沼氣工程的微生物組進行研究，發(fā)現(xiàn)其中包含大量迄今為止未知或未充分了解的物種，只有18%-25%的16S rRNA基因轉(zhuǎn)錄序列能被分類到屬級，這表明宏轉(zhuǎn)錄組的測序廣度還有待提高。截止到目前，被測序和研究的16S rRNA序列，5S和23S rRNA序列較為豐富，mRNA序列則相對較少。盡管宏轉(zhuǎn)錄分析的應(yīng)用中存在需要進一步細化的技術(shù)問題，但其已經(jīng)提供了一些新的生物學(xué)發(fā)現(xiàn)。Wang等［25］通過宏轉(zhuǎn)錄分析和關(guān)鍵酶的測定發(fā)現(xiàn)磁鐵礦誘導(dǎo)的直接種間電子傳遞（Direct interspecific electron transfer，DIET）可以部分取代種間氫轉(zhuǎn)移（Interspecies Hydrogen Transfer，IHT），清除電子轉(zhuǎn)移阻斷，改善乙酸營養(yǎng)型甲烷生成效率。

宏蛋白質(zhì)組學(xué)被定義為在特定位點的環(huán)境微生物群落的所有蛋白質(zhì)組成進行大規(guī)模表征［26］，為蛋白質(zhì)結(jié)構(gòu)、定位和蛋白質(zhì)-蛋白質(zhì)相互作用的研究提供最直接的證據(jù)。由于不同微生物群落的復(fù)雜性、自然環(huán)境的異質(zhì)性和干擾化合物的存在，很難提取出合適的蛋白質(zhì)部分進行分析［27］。盡管存在這些困難，蛋白質(zhì)組學(xué)在將微生物群落的遺傳多樣性和活動與其對生態(tài)系統(tǒng)功能的影響相關(guān)聯(lián)方面發(fā)揮著巨大的潛力［22］。例如，Jing等［28］首次用于混合培養(yǎng)的iTRAQ定量蛋白質(zhì)組學(xué)分析表明，在丙酸鹽甲烷化過程中，磁鐵礦誘導(dǎo)了涉及各種途徑的蛋白表達水平的變化，并發(fā)現(xiàn)丙酸代謝相關(guān)蛋白上調(diào)?，F(xiàn)階段宏蛋白質(zhì)組學(xué)仍是一個新興領(lǐng)域，其在未來對復(fù)雜微生物群落的研究可能得到進一步發(fā)展和廣泛應(yīng)用。

近些年代謝組學(xué)雖然發(fā)展迅速，但是仍然遠遠落后于基因組學(xué)和蛋白質(zhì)組學(xué)，尚未應(yīng)用于厭氧消化，本文不做討論。

表1 應(yīng)用于厭氧消化的微生物組學(xué)類型

3 常用生物信息學(xué)分析手段

常用的厭氧微生物群落生物信息學(xué)分析如圖2所示。從圖中看出，生物信息學(xué)分析測序結(jié)果包括幾個方面，物種組成分析，α多樣性分析，OTU聚類分析，多元統(tǒng)計學(xué)分析（PCA/ PCoA/ NMD，RDA/CCA），共現(xiàn)關(guān)系研究和代謝功能與通路研究。

圖2 厭氧微生物群落生物組學(xué)分析（部分圖片經(jīng)參考文獻［31-32］許可轉(zhuǎn)載，版權(quán)歸愛斯維爾、美國化學(xué)學(xué)會所有）

3.1 物種組成分析

物種組成分析是厭氧微生物群落最常用的生物信息學(xué)分析方法之一?？梢允褂肞hylophlan［37］、MetaPhlAn［39］、PyNAST［7］或者 RDP Classifier［46］等軟件對核心OTU序列和參考序列進行比對，并將比對后的序列進行系統(tǒng)發(fā)育組成推斷。常用的基因組 參 考 序 列 數(shù) 據(jù) 庫GenBank［46］、EBI、DDBJ［39］、IMG/M［47］。常用的物種注釋數(shù)據(jù)庫包括Greengene數(shù)據(jù)庫［31］、Silvia數(shù)據(jù)庫［41，47］。以圖2-a為例，微生物群落的系統(tǒng)發(fā)育分析可以從界、門、綱、目、科、屬、種等不同層次進行，運用Krona軟件［40］能夠展現(xiàn)出整個微生物群落的多級物種分類圖。

目前，學(xué)者們對不同操作條件下厭氧消化池微生物組成進行了大規(guī)模的分析。其中最常見的研究內(nèi)容是細菌和古菌的相對豐度。不同底物、接種物以及反應(yīng)參數(shù)下，反應(yīng)器的微生物組成各不相同［31-32］。但可以確定的是，隨著時間的推移，相同反應(yīng)器微生物豐度在一定時間內(nèi)的總體概況是相似的，群落演替遵循確定的趨勢［31］。

3.2 α多樣性分析

通常使用Mothur［20，34］、R語言中的Vegan包［31］或RDP［48］來研究產(chǎn)甲烷系統(tǒng)的微生物群落多樣性。常用的多樣性分析的指數(shù)（圖2-b）包括：Chao1豐富估計量，基于覆蓋度的物種豐富度（ACE），香農(nóng)多樣性指數(shù)（Shannon- Wiener index）和辛普森多樣性指數(shù)（Simpson diversity index）。其中Chao1、ACE和Shannon- Wiener index與微生物群落生物多樣性呈正相關(guān)，Simpson diversity index 與微生物群落生物多樣性呈負相關(guān)。Chao1和ACE指數(shù)可直接用于計算微生物群落中物種總數(shù)，Shannon-Wiener和Simpson diversity index還與群落均勻度有關(guān)，可用于估算生物多樣性。群落的均勻度越高，其產(chǎn)甲烷活性就越高。Wittebolle等［49］觀察到群落的均勻度與功能穩(wěn)定性相對應(yīng)。微生物多樣性結(jié)構(gòu)的均勻性確保了群落的多樣化代謝的能力。此外，物種的覆蓋率（Good’s coverage）越高，測序結(jié)果越真實。α多樣性分析各指標(biāo)的計算公式如表2所示。α多樣性的下降可作為反應(yīng)器運行狀況惡化的早期預(yù)警［47］。

表2 α多樣性分析參數(shù)的公式和符號意義［50］

3.3 OTU相似性分析

維恩圖（VennDiagram）可以用來表示OTU的共同和獨特的組成部分［32］。在通常情況下進行分析時，常選用相似水平為97%的OTU 樣品?？梢允褂肦 語言VennDiagram包［32］、Venny工具［47］或在線工具（http：//www.biovenn.nl/）繪制出圖。如圖2-c可以直觀地反映出厭氧反應(yīng)器中OTU的相似部分，3個反應(yīng)器的相同OTU比例為60.9%。結(jié)合其他討論，維恩圖的公共部分可能是該系統(tǒng)中發(fā)揮特定功能的物種［32］。

3.4 多元統(tǒng)計學(xué)分析

最常見的多元分析方法包括主成分分析（Principal Component Analysis，PCA）［43，47］，主坐標(biāo)分 析（Principal Coordinate Analysis，PCoA）［34，51］，非度量多維尺度分析（Non-metric Multidimensional Scaling，NMDs）［35，42］，冗 余 分 析（Redundancy Analysis，RDA）［7，21，37，48］和典范對應(yīng)分析（Canonical Correspondence Analysis，CCA）［51］。這些多元技術(shù)的共同之處在于，每一種技術(shù)本質(zhì)上都被看作是一種排序分析，其目的是描述相鄰對象的相似性。常用的計算工具包括R語言的Vegan包［7］、PC-ORD軟件［20］和Matlab［45］等。

PCA與PCoA以及NMDs是3種典型的非約束排序分析。其中PCA與PCoA的一個顯著區(qū)別在于，PCA基于相似系數(shù)矩陣（如歐式距離）確定主成分，而PCoA基于距離矩陣（如Bray-Curtis，UniFrac 和Jaccard）確定主坐標(biāo)（圖2-d）。另一個非約束分類分析NMDS算法能夠?qū)ο蟮木嚯x進行排序，并使用這些排序?qū)ο蠓蔷€性地映射到一個簡化的二維排序空間上，從而保持它們排序的差異而非原始距離［52］。與PCA和PCoA相比，NMDS算法相似度高的數(shù)據(jù)集更密集，相似度低的數(shù)據(jù)集更疏遠，這使得被測細菌和古菌群落剖面之間關(guān)系視圖更加清晰［52］。無約束排序分析在分析復(fù)雜宏基因組數(shù)據(jù)，或者微生物群落結(jié)構(gòu)差異方面的應(yīng)用越來越多［7，31，34］。Feng等［31］運用三維- PCoA分析發(fā)現(xiàn)不同污泥停留時間（SRT）的反應(yīng)器內(nèi)的微生物群落迅速從接種菌群中發(fā)生分化，且微生物群落動態(tài)在不同反應(yīng)器之間同步。Ning等［35］將NMDS結(jié)合參數(shù)聯(lián)合解釋分析發(fā)現(xiàn)，酸和的協(xié)同脅迫是抑制機制的主要驅(qū)動因素，并進一步發(fā)現(xiàn)細菌群落對相關(guān)負載蛋白更為敏感，有可能引起產(chǎn)甲烷群落從乙酸營養(yǎng)型產(chǎn)甲烷向氫營養(yǎng)型產(chǎn)甲烷的轉(zhuǎn)變。

在約束排序分析中，RDA和CCA是兩種常用的微生物與環(huán)境因子關(guān)聯(lián)性分析方法。這兩種方法非常相似，都可以使用箭頭代表環(huán)境因素，解釋環(huán)境梯度影響了系統(tǒng)發(fā)育結(jié)構(gòu)的變化情況。不同的是，RDA基于線性模型，而CCA基于單峰模型。以用Anoca軟件進行去趨勢對應(yīng)分析（Detrend Correspondence Analysis，DCA）分析結(jié)果作為選擇RDA和CCA的標(biāo)準［20，47］。若DCA第一軸大于4.0，CCA擬合效果更好；而在3.0-4.0之間，RDA和CCA均可；小于3.0，RDA的結(jié)果可信度更高。Feng等［31］運用CCA分析表明，隨著OLR的增加，碳水化合物、脂類、蛋白質(zhì)等有機大分子向可溶性有機分子的水解和轉(zhuǎn)化增加。反之丙酸和乙酸等可溶性中間產(chǎn)物濃度的升高可能表示產(chǎn)甲烷過程被部分抑制。

3.5 共現(xiàn)關(guān)系研究

厭氧消化菌群并非以單個OTU來執(zhí)行功能，而是通過厭氧食物鏈形成一個復(fù)雜的群落。為了描述種間關(guān)系，共現(xiàn)關(guān)系網(wǎng)絡(luò)常用于揭示物種-物種的關(guān)聯(lián)，描述核心物種之間的相互作用［38，47］。在這些研究中，共現(xiàn)網(wǎng)絡(luò)的拓撲結(jié)構(gòu)主要在Gephi［31］、Cytoscape［47］軟件中使用皮爾遜或斯皮爾曼相關(guān)性分析［53］和存在缺失數(shù)據(jù)的超幾何分布［54］得到了可視化。密集連接的節(jié)點被劃分為簇，這些簇被解釋為具有重疊生態(tài)位［53-54］的類群。

以圖2-e為例，當(dāng)兩個節(jié)點（OTUs或任何分類相關(guān)的單位）在多個樣本中同時出現(xiàn)或顯示相似的豐度模式時，則為正相關(guān)；反之則為負相關(guān)。節(jié)點的大小代表與其相關(guān)物種的多少，節(jié)點之間各邊的寬度與相關(guān)強度成正比［34］。然而這些關(guān)系的生態(tài)相關(guān)性并不容易解釋。例如，由于交叉取食、生物膜內(nèi)的共聚集、共定殖、生態(tài)位重疊等原因，都可能形成正相關(guān)關(guān)系；而由于偏害共棲、捕食-被捕食關(guān)系、競爭等原因，則可能形成負相關(guān)關(guān)系［55］。此外，更復(fù)雜的生態(tài)相互作用形式容易被忽視，例如滯后關(guān)系和多物種的共生（競爭）關(guān)系。

3.6 代謝功能與通路研究

為了深入了解分類組成之外的信息，可以使用生物信息學(xué)平臺或工具對宏基因組數(shù)據(jù)進行功能注釋。常用讀取基因組的工具為MG-RAST服務(wù)器［39］。常 用 功 能 注 釋 數(shù) 據(jù) 庫 包 括GO、KEGG［36，40，42］、COG［36，41］、Pfam［41］、SWISS-PROT［45］、UniProt-KB［45］。趙一全［42］利用宏轉(zhuǎn)錄組學(xué)在研究預(yù)處理秸稈對微生物的影響時，通過基于KEGG數(shù)據(jù)庫比較甲烷代謝相關(guān)途徑豐度時發(fā)現(xiàn)未處理樣本微生物群落能夠進行多種途徑的營養(yǎng)代謝活動，而經(jīng)過預(yù)處理的樣本中的甲烷代謝途徑較為單一。之后將注釋到的數(shù)據(jù)映射到代謝通路圖（類似于圖2-f）中，發(fā)現(xiàn)預(yù)處理并未改變代謝途徑，而是能夠?qū)⒛承┩分械闹虚g代謝產(chǎn)物的合成途徑變得更為直接。Heyer等［45］首次對40份沼氣廠樣本進行了大規(guī)模宏蛋白組分析。該方法比起基因組學(xué)研究，可以過濾出更具體的核心物種和核心函數(shù)。在這些系統(tǒng)研究的基礎(chǔ)上，可以確定某種蛋白作為生物標(biāo)志物。如果這些生物標(biāo)記物的豐度迅速變化或表現(xiàn)出與選擇的工藝參數(shù)相矛盾，這可能是之后過程失敗的標(biāo)志。

4 厭氧消化微生物組學(xué)的研究進展

4.1 微生物群落與運行參數(shù)關(guān)聯(lián)分析

微生物組學(xué)分析將微生物群落組合模式與運行參數(shù)相聯(lián)系，可以指導(dǎo)厭氧消化參數(shù)的調(diào)節(jié)。但目前用來解釋微生物群落組合模式的理論仍存在一定的分歧。其中占主導(dǎo)地位的理論主要分兩種：生態(tài)位理論和中性理論。經(jīng)典的生態(tài)位理論強調(diào)了微生物的種間關(guān)系（如共生和競爭）生態(tài)位分化和底物營養(yǎng)等確定性過程的重要性［29］。中性理論則只考慮出生、死亡、漂移、移民等隨機過程對生態(tài)環(huán)境的影響［56］。大部分研究認為生態(tài)位理論和中性理論對構(gòu)建土壤和生物反應(yīng)器的微生物群落非常重要［57］。微生物群落的組成并不完全是由環(huán)境條件決定，而是在一定程度上依賴中性過程［56］。了解微生物群落組成對生態(tài)位理論和中性過程及對環(huán)境條件的依賴度，對設(shè)計和維護厭氧消化工程有著極大的意義。

由于群落組合模式的復(fù)雜性以及影響微生物群落結(jié)構(gòu)的參數(shù)的多樣性，理解環(huán)境參數(shù)與微生物的相互作用變得困難［58］。一般來說在所有厭氧消化系統(tǒng)中均存在一個核心菌群［37］。Tao等［59］從20個運行了長達兩年的實驗室規(guī)模的厭氧消化反應(yīng)器中采集了138個樣本，并得出了核心菌群的組成（Bacillus、Clostridium、Bacteroides、Eubacterium、Cytophaga、Anaerophaga和Syntrophomonas等），而產(chǎn)甲烷菌的組成因接種源、反應(yīng)器溫度或鹽度不同而不同。核心菌群與沼氣產(chǎn)量有顯著的相關(guān)性。產(chǎn)甲烷量高的反應(yīng)器，其細菌群落非常相似。相對應(yīng)的是，古菌群落由于功能的變異性較大而具有較強的多樣性［37］。此外，有些物種雖然稀有，但也是厭氧消化過程的重要驅(qū)動力，如Hyphommicium、Gordonia、Methanospirillum等［48，60］。

4.2 代謝途徑分析

最近的研究表明，厭氧消化所涉及的細菌群落具有高度的恢復(fù)力（Resilience）。水解、產(chǎn)酸的每一個主要步驟都有多個相關(guān)的系統(tǒng)發(fā)育群落，功能冗余在所難免［29］。較高水平的發(fā)酵菌群（如互營細菌）的動態(tài)變化與初級發(fā)酵菌群（如Clostridia和Bacteroidetes）豐度變化明顯不同，它們更依賴冗余來維持群落的整體功能［37］。均勻度較高的群落利用冗余功能通路的潛力更高，即具有更多并行通路的群落對底物代謝的反應(yīng)更有效［7］。這使得以特定代謝方式調(diào)控細菌菌群變得困難［61］。

作為厭氧消化過程中甲烷生產(chǎn)的最終步驟，產(chǎn)甲烷古菌代謝途徑較為單一。目前所有已知的產(chǎn)甲烷菌都表達甲基輔酶M還原酶（mcr）［62］。最常見的乙酸營養(yǎng)型產(chǎn)甲烷菌參與的乙酸代謝過程如圖3-A所示：首先，以乙酸受到乙酸激酶（ack）、磷酸乙酰轉(zhuǎn)移酶（pta）或乙酰輔酶A合成酶（acs）等酶催化形成乙酰輔酶A。然后乙酰CoA 再在乙酰輔酶A脫羰酶/合酶復(fù)合體（cdh）的作用下形成甲基四氫甲烷喋呤。此外，氫營養(yǎng)型產(chǎn)甲烷菌參與了完整的CO2還原途徑的所有通路（圖3-B）。涉及這些通路的重要酶包括：四氫甲烷蝶呤甲酰轉(zhuǎn)移酶（hfr）、次甲基甲酰四氫甲烷蝶呤環(huán)化水解酶（mch）、亞甲基四氫甲烷蝶呤脫氫酶（mtd）以及5，10-亞甲基四氫甲烷蝶呤還原酶（mer）等［63］。此后氫營養(yǎng)型產(chǎn)甲烷途徑與乙酸營養(yǎng)型產(chǎn)甲烷途徑相同，經(jīng)四氫甲烷蝶呤S-甲基轉(zhuǎn)移酶（mtr）和甲基輔酶M還原酶（mcr）的催化下形成甲烷［40］。另外，一些新的生化途徑被發(fā)現(xiàn)，包括DIET［25］，膜結(jié)合的反向電子轉(zhuǎn)移，以及以氫和/或甲酸鹽的形式處理電子的胞質(zhì)凝聚型酶［64］。厭氧消化的代謝過程的“面紗”正逐漸地被揭開。

圖3 產(chǎn)甲烷途徑

4.3 功能生物標(biāo)記物鑒定

隨著宏基因組測序技術(shù)的廣泛使用，越來越多的厭氧消化相關(guān)的功能生物標(biāo)記物被鑒定出來，從而簡化了基因組研究。mcrA基因的a亞基是研究最多的產(chǎn)甲烷菌的功能基因，其拷貝數(shù)［65］和轉(zhuǎn)錄本［66］與甲烷產(chǎn)量存在顯著正相關(guān)，被認為是厭氧消化中最重要的生物標(biāo)志物。輔酶F430作為mcrA基因的輔酶同樣可以量化產(chǎn)甲烷潛力［66］。Fe-Fe-氫化酶大亞基的基因（hydA）被用作產(chǎn)氫發(fā)酵細菌生物標(biāo)記物［67］。同型產(chǎn)乙酸細菌以甲酰四氫葉酸合成酶編碼基因（FTHFS）（Wood-Ljungdahl通路的關(guān)鍵酶）作為生物標(biāo)志物［68］。此外，基因組研究還可鑒定編碼關(guān)鍵微生物基因組的特定代謝途徑的酶，包括碳水化合物利用、脂肪酸降解、氨基酸發(fā)酵等［69］。通過16S rRNA基因測序與已培養(yǎng)的相似物種做比對，也可以推測出厭氧消化過程的功能［7］。但由于密切相關(guān)的生物可能在功能上有所差異，因此應(yīng)謹慎做出結(jié)論［70］。如果基于相關(guān)性分析來建立假設(shè)，并且使用補充技術(shù)進一步論證，得出的結(jié)果將更加可靠。

4.4 抗生素耐藥基因分析

研究沼渣中抗生素耐藥基因（Antibiotic Resistance Genes，ARGs），可以界定其給土壤等其他環(huán)境因素可能帶來的風(fēng)險［71］。沼渣中的ARGs可能會傳播給致病菌，從而降低致病菌在醫(yī)療過程中對抗生素的敏感性，對公眾健康構(gòu)成威脅。ARGs可通過水平基因轉(zhuǎn)移在不同微生物種群間傳播，形成具有抗生素耐受性的細菌。使用宏基因組測序不僅揭示了微生物水平基因轉(zhuǎn)移的證據(jù)，也證明了其耐藥性傳播機制。Lee等［72］對具有代表性的有機廢棄物（糞便、污泥和食品垃圾回收廢水）中ARGs的多樣性和數(shù)量進行研究，發(fā)現(xiàn)糞便中ARGs的總量最高，其次是污泥和食品垃圾回收廢水。不同基質(zhì)間ARGs的多樣性和作用機制存在較大差異。研究發(fā)現(xiàn)，厭氧消化能有效地降低ARGs拷貝數(shù)，但不同抗性機制下的相對豐度則存在明顯差異［72］。某些細菌（如Clostridium和Nitrosomonas）可能是多個ARGs的潛在宿主。Luo等［39］發(fā)現(xiàn)，在以糞肥和工業(yè)廢料為原料的反應(yīng)器中，ARGs的總豐度在7×10-3-1.1×10-1拷貝數(shù)/ 16S rRNA基因拷貝數(shù)之間變化，并且嗜熱型厭氧消化器在去除ARGs上有一定優(yōu)勢。除此之外，延長固體停留時間以及添加生物炭都可以降低ARGs的環(huán)境風(fēng)險［73］。

5 總結(jié)與展望

厭氧消化是一種生態(tài)友好的有機廢棄物處理及資源回收過程。在過去的十年中，隨著微生物組學(xué)的快速發(fā)展，我們對厭氧消化中微生物菌群的結(jié)構(gòu)和功能的認識也有了質(zhì)的提升。這不僅能夠更好地建立微生物群落模型，而且也將使人們更清楚地了解相關(guān)的環(huán)境條件、代謝途徑和演替規(guī)律，從而解釋這種格局。微生物組學(xué)研究為曾經(jīng)被認為是“黑盒”的復(fù)雜微生物群落提供了光明。本文綜述了目前微生物組學(xué)的類型和常用的生物信息學(xué)分析手段，以及最近的研究進展，有助于確定最佳的沼氣生產(chǎn)因素的組成，改進厭氧消化器的操作參數(shù)；也有助于充分利用厭氧微生物群落潛力，為設(shè)計最佳厭氧消化工藝提供理論依據(jù)。

目前為止，通過厭氧消化微生物組來改善厭氧消化產(chǎn)甲烷效率仍面臨許多挑戰(zhàn)。首先，微生物組的厭氧微生物群落研究仍處于起步階段，需要儲存已有數(shù)據(jù)作為參考。并且新的研究需要讀入大量的微生物組數(shù)據(jù)，為之前的研究結(jié)果提供統(tǒng)計支持，這導(dǎo)致需要為測序數(shù)據(jù)提供龐大的存儲空間；其次，大多數(shù)序列分析需要與參考數(shù)據(jù)庫中相似序列進行比對，但仍有部分基因在現(xiàn)有數(shù)據(jù)庫中無法匹配到相似序列。未來有必要擴大參考數(shù)據(jù)庫，建立能夠?qū)蚪M學(xué)數(shù)據(jù)進行全面和綜合分析的平臺；最后，未來的研究不僅需要考慮微生物的代謝，還需要考慮微生物代謝物如何改變自身或其他微生物的基因潛力與表達，整合多組學(xué)分析也是必不可少的?？偠灾胪ㄟ^厭氧微生物群落的微生物組學(xué)信息來精確調(diào)節(jié)厭氧消化反應(yīng)器的性能，仍有很長的路要走。而研究微生物群落如何響應(yīng)性能參數(shù)，已經(jīng)是向目標(biāo)邁出了第一步。