洪义国, 吴佳鹏
广州大学大湾区环境研究院, 珠三角水质安全与保护教育部重点实验室, 广东 广州 510007
随着世界石油、天然气资源的日渐耗尽, 各国科学家正在努力寻找新的替代能源。目前的研究已经发现, 海底深部富含甲烷天然气水合物, 其含量大约是目前地球化石能源的两倍(Etiope et al,2014; Suess, 2014; Joye, 2020)。天然气水合物被称为是21 世纪最具有商业开发前景的战略资源, 或将改变全球能源结构, 正受到各国科学家和各国政府的重视。此外, 甲烷又是具有强烈温室效应的气体, 其温室效应大约是CO2的25 倍。在全球变暖和能源短缺的背景下, 富含水合物的海相沉积物中古菌的厌氧甲烷氧化(anaerobic oxidation of methane, AOM)作用日益受到科学家关注(Chen et al, 2019)。近年来, AOM 过程的研究取得了很多进展, 本文将从AOM 过程的研究历史、甲烷氧化古菌(methane-oxidizing archaea, MOA)的系统进化、AOM 的代谢过程机制、AOM 在深海冷泉碳酸盐形成以及控制海底甲烷渗漏的生物地球化学功能方面进行评述。
长期以来, 甲烷氧化被认为是在有氧条件下由嗜 甲 烷 细 菌( 主 要 为 α-Proteobacteria 和 γ-Proteobacteria)介导完成的, 细菌从这一代谢过程中获得能量(Hanson et al, 1996)。然而, 在厌氧条件下能够氧化甲烷的微生物却一直没有被发现, 因此在20 世纪的大部分时间里, 很多微生物学家坚持认为甲烷在厌氧条件下是惰性的。直到20 世纪70 年代, 一些开创性的研究使得这一研究领域有了新突破。在1976 年, Reeburgh (1976)发现, 在卡里亚克海槽的厌氧环境中, 甲烷代谢与硫酸盐还原菌(sulfate-reducing bacteria, SRB)密切相关; 在同一年,Barnes 等(1976)的研究发现, 在海洋沉积物中, 甲烷浓度从沉积物到水体方向渐逐降低, 而硫酸盐浓度从海体到沉积物方向也逐渐降低, 这表明硫酸盐可能是AOM 过程的电子受体。在20 世纪80 年代,一系列研究表明这种甲烷和硫酸盐在海洋沉积物中的分布具有普遍类似的特征, 形成了这样的假说:AOM 是被甲烷氧化古菌和硫酸盐还原菌共同作用的结果(Zehnder et al, 1980), 但是对于其作用机制还不明确。到了20 世纪初, 由于Delong 研究小组的出色工作, 使得这一领域的研究有了重大突破。他们通过稳定同位素分析、基因文库构建以及荧光探针杂交(fluorescence in situ hybridization, FISH),为这个假说的成立提供了有力的证据, 明确了AOM 是由微生物作用的生物地球化学过程(Boetius et al, 2000)。Orphan 等使用FISH 与二级离子质谱技术(fluorescence in situ hybridization-secondary ion mass spectroscopy, FISH-SIMS)以及FISH 偶联免疫磁细胞捕获技术, 使得海洋冷泉古菌厌氧甲烷氧化的研究获得了重要突破, 成为与能源和环境问题密切相关的研究热点(Orphan et al, 2001a, b, 2002;Pernthaler et al, 2008)。
到目前为止, 还没有获得厌氧甲烷氧化菌的纯培养菌株, 所以对于海底冷泉厌氧甲烷氧化的研究都是基于分子生态学的方法进行的。综合目前的16S rRNA 基因和mcr-A 功能基因研究结果, 参与AOM 过程的微生物包括以硫酸盐为电子受体的甲烷氧化古菌、以亚硝酸盐为电子受体的甲烷氧化古菌, 以及以金属离子为电子受体的甲烷氧化古菌。
以硫酸盐为电子受体的甲烷氧化古菌包括三个类群: ANME-1、ANME-2 和ANME-3。其中ANME-1是最早鉴定的厌氧甲烷氧化古菌。研究者对加利福尼亚北海岸Eel 河盆地的研究发现, 冷泉区沉积物中存在古菌特异性的古脂(Archaeol), 而且其 δ13C值强烈亏蚀, 用 16S rRNA 基因分析鉴定为ANME-1。 ANME-1 类群在系统进化上分为ANME-1a 和ANME-1b 两个亚群, 与甲烷微菌属(Methanomicrobiales) 和 甲 烷 八 叠 球 菌 属(Methanosarcinales)关系较近, 但在进化树上却有独立的分支(Boetius et al, 2000)。随后的研究进一步发现, ANME-2 和ANME-3 两个甲烷古菌类群也参与了 AOM 的代谢过程, 而且 ANME-2 又分为ANME-2a、ANME-2b 和ANME-2c 三个亚群, 在进化上与甲烷八叠球菌属较近; ANME-3 与甲烷球菌属 (Methanococcoides)较近(Orphan et al, 2001a,2002; Pernthaler et al, 2008)。这三个甲烷氧化古菌类群的 16S rRNA 基因序列差异较大, 其同源性在80%左右, 表明具有逆向甲烷生成的甲烷古菌向着多个方向进化。最近的研究进展也表明多种古菌参与了厌氧甲烷氧化的过程, 例如嗜热古菌Candidatus Ethanoperedens (Hahn et al, 2020)和Asgard 古菌(Seitz et al, 2019)。甲基辅酶M 还原酶(Methyl coenzyme M reductase, MCR)是甲烷生成过程中的末端酶, 催化辅酶 B(H-S-CoB)和甲基辅酶M(CoM-S-CH3)还原, 产生甲烷和CoM-S-S-CoB。2002 年, Luton 等比较了甲烷古菌的mcrA 基因和16S rDNA 基因的系统进化关系, 发现甲烷古菌mcrA 基因与16S rDNA 基因在系统进化上的关系非常吻合, 发展了mcrA 基因作为检测甲烷古细菌多样性和分布的功能基因标记(Luton et al, 2002)。2003年以后的大量研究表明, mcrA 基因包含a、b、c、d、e 五个亚群。由于mcrA 基因与MOA 密切相关, 因此成为研究厌氧甲烷氧化古菌的一个有力分子工具。通过综合的系统进化对比分析, 发现c、d 的mcrA 亚群与ANMEⅡ古菌相对应, a、b 的mcrA 亚群与ANMEⅠ古菌相对应, e 的mcrA 亚群可能与ANME Ⅲ或其它甲烷古菌相对应(Hallam et al,2003; Dhillon et al, 2005; Chistoserdova et al, 2005;Nunoura et al, 2006)。
以金属离子(Fe3+、Mn4+)作为电子受体的甲烷氧化古菌目前还不是很明确, 初步认为ANME-2 及其ANME-2a、ANME-2c、ANME-2d 均可参与此类反应(Ettwig et al, 2016)。随着研究的继续深入, 应用16S rRNA 基因结合mcrA 功能基因标记, 可能会陆续鉴定出更多类群的厌氧甲烷氧化古菌。
根据地球化学和群体细胞外学特征, 科学家建立了“逆甲烷产生”过程的假说, 认为相关的产甲烷古菌在特殊环境条件下进化出反向甲烷产生的能力,甲烷氧化过程为菌体细胞的生长和繁殖提供必要的碳源和能源。
Delong 的研究小组利用环境基因组的方法对这一假说的验证提供了证据(Hallam et al, 2004。他们利用密度梯度离心的方法, 从加利福尼亚海岸Eel河盆地的冷泉样品中富集菌体细胞, 然后建立了全基因组鸟枪DNA 文库(whole-genome shotgun, WGS)和fosmid DNA 文库。对这两种文库分析后发现, 在样品微生物群体中, ANME-1、ANME-2 和SRB 菌体细胞占绝对优势。而且除mer 基因(Methenyltetrah ydromethanopterin reductase)外, 在文库中包含相对丰富的mcr、mtr、mtd、mch、ftr 和fmd 编码与甲烷产生相关酶的所有基因。另外在文库中发现了有编码F420依赖的呼吸链相关蛋白基因的存在。这些分析结果为“逆甲烷产生”过程的假说的成立提供了有力证据, 这种在热力学上很难进行的反应可能通过偶联F420依赖的呼吸链电子传递产生能量而得以克服。基因组学的研究角度为这种特殊海洋环境中的AOM 生物地球化学过程提供了更深入的理解(Chistoserdova et al, 2005; Evans et al, 2015)。
生物化学的研究也为这一假说提供了进一步的证据。2003 年, Krüger 等报道了一个镍蛋白参与了AOM 的过程, 他们分析了一个黑海西北陆架的冷泉样品, 发现样品中有约70%的ANME-1 菌群, 25%的SRB, 其余5%为其他微生物。在12℃条件下进行培养, 该菌群依赖硫酸盐的甲烷氧化率为0.1nmol·min-1·mg protein-1。色谱分析显示, 在样品菌体抽提物中包含分子量为905kDa 和951kDa 的两种镍复合物, 组成比为1: 4。较轻的镍复合物经鉴定与F430相同, 但较重的镍复合物是以前没有发现的(Krüger et al, 2003)。而且在北太平洋Hydrate Ridge和波罗的海的Eckernfiöde Bay 的冷泉样品中也同样发现了951kDa 的镍复合物, 然而在一个产甲烷的污泥中却只有905kDa 的镍复合物。经过硫酸铵沉淀和离子交换色谱, 一种包含951kDa 镍复合物的蛋白被分离纯化(Ni-Ⅰ蛋白), 这一蛋白占到样品中抽提蛋白的7%。进一步分析发现, 与MCR 蛋白(Ni-Ⅱ蛋白)一样, 这一蛋白由三个亚基组成, 其氨基酸序列与从该样品中鉴定的一个34.8kb 的DNA 片断上的三个ORF 的编码高度一致, 而且与产甲烷古菌的MCR 酶高度同源。其中β 和γ 亚基非常保守, 其活性位点与MCR 基本一致。但α 却有显著差异, 在MCR α 亚基上的一个活性中心位点(VX2CCX4CX5C)与Ni-Ⅰ蛋白α 亚基相应的序列又很大不同。Ni-Ⅰ蛋白α 亚基与其它MCR 的α 亚基的系统进化分析发现, 这一蛋白在进化上与MCR α 有了明显较大的距离。这也表明, 在ANME 古菌进化过程中, Ni-Ⅰ蛋白的 α 亚基也随之产生了平行进化(parallel evolution)。综合以上分析得出结论, ANME 古菌的Ni-Ⅰ蛋白参与了AOM 的代谢过程, 从分子水平为AOM 的反甲烷产生的假说提供了有力的证据。
一般认为甲烷厌氧氧化是由甲烷氧化古菌和硫酸盐还原细菌共同完成的一个两步反应过程, 其反应方程可以表示为: CH4+ SO2- 4→ 5CO2+ HS+H2O, 这种微生物细胞之间的相互作用, 称之为互养作用(syntrophy)。通过互养, 几个代谢过程被同时驱动, 因此有人认为互养作用在生物细器的进化中具有中心作用(Martin et al, 1998)。2015 年, 美国Orphans 实验室报道了有关AOM 过程中甲烷氧化古菌和硫酸盐还原细菌互养过程机制的研究, 从多方面提供了重要证据(Green-Saxena et al, 2014), 而且在冷泉微生物群体中还发现了真菌的存在(Sapir et al, 2014)。在从美国俄勒冈州海岸卡斯卡底古陆海底甲烷冷泉区的样品具有高的依赖甲烷的硫酸盐还原率。16S rRNA 分析表明, 样品中主要为甲烷古菌和硫酸盐还原菌。通过16S rRNA 特异性的寡核苷酸荧光探针定点杂交分析, 发现菌群中有明显的结构分布: 大约100 个古菌细胞聚居在菌团的中央, 周围被SRB 细菌包围, 形成近似球形结构(Boetius et al, 2000)。Orphan 等使用FISH、FISH-SIMS 以及FISH-IMCC 的方法检测到其中的古细菌群体δ13C值高度亏损(-96‰)表明古菌群体是甲烷营养型的,从而说明甲烷氧化是由菌群中的古菌负责完成的;而且发现嗜甲烷菌与硫酸盐还原菌在代谢上相互依赖(Dekas et al, 2016)。
在MOA 与SRB 的互养过程中, 甲烷氧化与硫酸盐还原是两个相对独立的过程, 那么把这两个过程联系到一起形成共养作用(syntrophy)机制是什么呢?有人推测, MOA 氧化甲烷产生一种未知的中间电子载体, 然后SRB 以硫酸盐为氧化剂还原这种电子载体, 这种中间电子载体可能为H2、甲酸盐、乙酸盐或者其他化合物(图1)。对这一问题目前还没有统一的结论, 需要更多生物化学以及热力学的理论探索。除了硫酸根作为电子受体外, 在理论上硝酸根、三价铁、锰等也都可以作为末端电子受体氧化甲烷。在2006 年, Marc Strous 团队报道一个微生物菌群能够厌氧甲烷氧化偶联反硝化。这一微生物菌群是从厌氧环境中分离到的, 包含两种菌: 一种为到目前为止还没有纯培养的细菌, 另外一种为海洋甲烷古菌。海洋甲烷古细菌进行AOM 代谢, 而细菌进行反硝化反应, 这个偶联反应可以表示为:5CH4+ 8NO3+ H+→ 5CO2+ 4N2+14H2O (每升甲烷氧化释放的自由能ΔG0= -765kJ·mol-1CH4), 这个反应直接把碳和氮的物质循环紧密联系起来, 推测在整个自然界的生物地球化学物质循环中发挥重要作用(Raghoebarsing et al, 2006)。最新的单细胞实验证明了MOA 与SRB 的互养过程中发生了直接的电子转运(McGlynn et al, 2015), 而且发现铁氧化能显著促进MOA 与SRB 的互养代谢活性(Sivan et al, 2014;Glass et al, 2014)。
图1 厌氧甲烷氧化过程Fig. 1 Schematic diagram of reverse methanogenesis,X-electron transporter
对于金属离子驱动的AOM 的反应机理目前还不是很明确, 根据研究推测可能的机制有三种: 第一种认为微生物还原金属离子直接耦合甲烷厌氧氧化, 与硫酸盐还原耦合甲烷厌氧氧化机制相类似(Norði et al, 2013); 第二种是在硫化物存在条件下,形成零价硫, 零价硫通过歧化反应生成硫酸盐和负二价硫, 其中硫酸盐可以驱动甲烷氧化的发生(Sivan et al, 2011; Riedinger et al, 2014); 第三种是反应体系中的H2被Fe3+还原菌利用至低浓度后, 通过反向产甲烷作用氧化甲烷(Sivan et al, 2011)。
在ANME 和SRB 互养的过程中, 其ΔG 值随着甲烷浓度的变化呈现很大的差异。当甲烷压力在0.1MPa 时, 释放的自由能为ΔG 值为-30kJ·mol-1;而当甲烷的浓度为 8MPa 时, ΔG 值达-40kJ·mol-1(Boetius et al, 2000)。根据吉布斯(Gibbs)自由能计算, 由ADP 和PO34-合成ATP 需要的能量为 50~70kJ·mol-1, 生命最小的吉布斯自由 能为20kJ·mol-1(Schink, 1997)。这样, 发生在浅水中的AOM 反应过程就很难达到生命所需要的最小量值。但是在深海高压条件下, AOM 过程的ΔG 值会升高,而且如果甲烷氧化的产物迅速得到转移和消耗, 就会使得反应顺利进行。尽管这样, 偶联SO2- 4 还原的 AOM 过程产生的能量仍然很低, 而且又要被ANME 和SRB 共享, 所以生长都会非常缓慢, 是已知的维持生命活动的最小代谢方式。总之, 对于ANME 与SRB 互养的AOM 过程机制仍然不清楚,有待进一步研究。
在富含天然气水合物的冷泉区域, 微生物的生物地球化学代谢过程相当活跃。在许多天然气渗漏发育区域, AOM 作用引起环境碱度增加, 进而沉淀出冷泉碳酸盐岩, 在海底表层沉积物中形成块状碳酸盐岩结壳(Michaelis et al, 2002; Marlow et al, 2014;Meister et al, 2018)。因此, 海底冷泉碳酸盐岩是天然气渗漏经微生物作用的结果, 是海底天然气渗漏作用的标志。
一般而言, 间隙水甲烷和硫酸盐浓度、CO2及甲烷碳同位素、甲烷厌氧氧化速率和硫酸盐还原速率等参数在沉积剖面上的变化可以指示海洋沉积物中存在AOM 反应带, 该反应带亦称作为硫酸盐-甲烷过渡带(sulfate-methane transition, SMT)。SMT 通常有以下特征: (1)甲烷浓度剖面在硫酸盐还原带的上部都非常低, 直到在界面附近由于AOM 作用使硫酸盐被消耗殆尽后才快速变高; (2)由于AOM 作用, 短时间内甲烷被快速消耗, 使甲烷浓度在剖面上呈现下凹状; (3)SMT 界面附近, 由于AOM 过程不断消耗12CH4, 致使13CH4不断被富集, 使得CH4的δ13C 值不断升高, 从约-90‰升高到-20‰, 所以MOA 菌体细胞的δ13C 值一般较低, 而SRB 的δ13C值较高, 这些已经被大量实验所证实(Nauhaus et al,2007; Wegener et al, 2008)。2013 年, 在东沙西南深水海域发现3 处冷泉碳酸盐岩区, 并采集到了大块体冷泉碳酸盐岩, 发现了活体双壳类冷泉生物及冷水珊瑚, 初步证实在该海域存在现代冷泉生态系统。中国科学院深海科学研究与工程研究所的彭晓彤研究员团队应用纳米离子探针(Nano-SIMS)以及同步辐射X 射线近边吸收谱等技术手段, 对冲绳海槽海底冷泉区富铁碳酸盐烟囱体开展了系统研究,在微米尺度上获取了富铁碳酸盐烟囱体中的硫同位素信息并判别了硫的来源, 同时结合铁同位素和碳同位素的数据, 从冷泉沉积环境中识别了由铁氧化驱动的厌氧甲烷氧化过程(Fe-driven AOM)过程形成的碳酸盐烟囱体, 分析了Fe-driven AOM 过程对碳酸盐烟囱体形成的影响, 提出了海底冷泉区富铁碳酸盐烟囱体的形成模式, 揭示了Fe-driven AOM 过程对于海洋冷泉环境中甲烷氧化的重要性, 表明在高铁输入的海域, Fe-driven AOM 是阻碍沉积物中甲烷向上层水体扩散的另一道重要屏障(Peng et al,2017)。此外, 最新的研究表明锰还原与甲烷的氧化也密切关联(Leu et al, 2020)。
基于以上的结论, 在深海地区寻找甲烷冷泉的探索中, 碳酸盐是其存在的重要标志性物质。而由于 AOM 过程与碳酸岩的形成密切相关, 因此ANME 的存在也是甲烷冷泉存在的重要间接证据之一。ANME 的特异性古脂、16S rDNA 以及mcr-A功能基因标记都可以作为分子工具辅助寻找深海天然气水合物。
在过去的200 多年中, 地球大气中的甲烷含量在逐渐上升, 其浓度是200 年前的两倍, 已经占到大气温室气体效应的20%(Frankenberg et al, 2005)。然而, 据1992 年IPCC 报告, 海洋对大气圈中甲烷的贡献仅为2%(Houghton et al, 1992), 但海洋沉积物实际产生的甲烷量远高于这个比例。Reeburgh(1996)的研究表明, 黑海沉积物中由于生物过程产生的甲烷每年约为2.9×1011mol, 而释放进入大气的仅为4.1×109mol, 99%的甲烷在沉积物中被消耗; Blair 等(1995)对亚马逊大陆架沉积物的研究也发现, 几乎全部的甲烷在沉积物中被消耗。那么“失踪”的甲烷去哪里了?诸多研究证实AOM 过程是海洋沉积物中 CH4消耗的主要途径(Burns,1998; Moore et al, 2004; Treude et al, 2005)。AOM 过程有效阻止海洋沉积物产生的甲烷进入沉积物-水体界面。据Reeburgh 和Alperin 估算, 海洋AOM 过程年消耗甲烷约为25~94Tg, 占大气圈各类甲烷源年排放总量的5%~20%, 故其在甲烷的全球收支平衡及大气温室效应过程中扮演极为关键的角色(Reeburgh, 2007)。因此AOM 在全球大洋海底缺氧带(oxygen minimum zone, OMZ)广泛分布, 在全球甲烷循环和全球碳循环中发挥重要作用。另外, 由于AOM 过程同时消耗沉积物间隙水中甲烷和硫酸盐, 产生无机碳和挥发性硫, 继而影响海洋沉积物碳和硫的生物地球化学循环(Hensen et al, 2003)。
AOM 是涉及能源和环境的重要古细菌代谢过程, 对这一现象已经进行了30 多年的研究, 特别是近10 年的研究, 对AOM 的认识才有所突破。尽管AOM 的意义重大, 但对于其生物学作用机制知道的还很少, 仍然存在诸多需要解决和完善的问题,未来的研究可以从以下几方面来进行。
1) AOM 菌的纯培养研究。到目前为止, 还没有获得ANME 的纯培养菌株, Girguis 等(2003)在持续流生物反应器中得到了ANME 的富集培养, 为进一步研究AOM 过程的生物化学特征提供了基础。只有得到了ANME 的纯培养菌, 才能从细胞学水平深入研究AOM 的作用机制。
2) AOM 的“逆向甲烷生成”和与SRB 的互养机制理论的进一步完善。虽然从多方面为AOM 的作用机制提供了证据, 但要更清楚的弄清这一过程还需要更多细胞水平上的直接证据。
3) AOM 过程的生态作用估价。AOM 在冷泉区碳酸盐生成和控制甲烷渗漏的自然生态过程的作用,是AOM 过程研究的一个核心问题。尽管对于黑海、波罗的海和墨西哥的加利福尼亚湾等几个海区的AOM 过程进行了研究, 但还需要其他相关海域的研究来提供更多数据。
4) 进一步研究AOM 菌的生态分布与天然气水合物分布的地化关系, 建立有效的检测ANME 菌的分子生物学方法。
国内有关AOM 的研究还不多。西沙海槽位于南海北部陆坡区的新生代被动大陆边缘型沉积盆地,新生代最大沉积厚度超过7000m。综合资料表明,南海陆坡和陆隆区应有丰富的天然气水合物矿藏,估算其总资源量达643.5~772.2 亿吨油当量, 大约相当于我国陆上和近海石油天然气总资源量的1/2(Chen et al, 2005)。然而, 天然气水合物在南海区的确切分布及特征目前仍不清楚。由于AOM 过程与冷泉碳酸盐的形成密切相关, 因此研究ANME 菌在我国沿海海区的分布可以为寻找天然气水合物储藏区提供科学参考。