塔克拉玛干沙漠防护林对土壤古菌群落与甲烷氧化潜力的影响

《Microorganisms》：Effects of Shelterbelts on Soil Archaea Community and Methane Oxidation Potential in the Taklimakan Desert

　　

摘要：塔克拉玛干沙漠极端干旱，土壤碳和水资源匮乏。尽管防护林建设能够改变土壤环境，但在植被恢复过程中，古菌群落与甲烷氧化指标之间的关系尚不清楚。为探究这一关系，研究人员从裸沙区和防护林地采集了土壤样品，包括0–10 cm表土层和10–40 cm亚表层。测定了土壤理化性质、古菌群落组成、pmoA丰度以及甲烷单加氧酶（MMO）相关指标。研究了防护林引起的土壤环境变化、古菌功能潜力与甲烷氧化之间的关联。防护林地土壤环境发生了显著变化；表土pH值从8.21降至7.73，而土壤有机碳（SOC）增加了146.88%。亚表层土壤含水量（SWC）增加了846.15%。防护林地检测到更高的古菌多样性。表土的香农指数上升了26.21%，与碳氮循环相关的微生物多样性也同时增加。防护林土壤中的pmoA基因拷贝数和MMO相关指标显著高于裸沙区，表明其甲烷氧化的功能潜力升高。随机森林分析显示，与目标功能指标相关的环境因子在不同土层间存在差异，且未发现单一强相关因子。这些发现增进了对沙漠人工植被建立后土壤古菌及甲烷循环相关功能的理解，并为沙漠公路防护林的生态恢复及干旱区土壤碳循环研究提供了参考。

塔克拉玛干沙漠防护林对土壤古菌群落与甲烷氧化潜力的影响——论文解读

一、研究背景与科学问题

沙漠约占全球陆地面积的41%，是陆地生态系统的重要组成部分。中国是世界上沙漠面积最大的国家之一，沙漠化土地超过国土总面积的34%。在干旱地区，土地沙漠化是一种普遍的环境现象，严重威胁区域乃至全球的生态安全与可持续发展。沙漠土壤结构松散，保水保肥能力差，生态功能脆弱，自我再生能力极为有限。作为干旱半干旱沙漠地区生态恢复和荒漠化治理的关键措施，防护林走廊的建设能够人工创建稳定的植被群落，改善土壤条件、水文过程和生物多样性，并在遏制沙漠化蔓延、维持区域生态平衡方面发挥至关重要的作用。

甲烷（CH₄）和二氧化碳（CO₂）都是重要的温室气体。在百年时间尺度上，甲烷的全球增温潜能约为二氧化碳的28–34倍，占全球总增温潜能的约20%。土壤甲烷循环由产甲烷和甲烷氧化两个耦合过程构成。产甲烷作用由古菌域中的产甲烷古菌完成，而甲烷氧化主要由细菌域中的甲烷氧化菌介导。沙漠土壤通常是好氧环境，不利于产甲烷古菌生存，但这些微生物仍能在局部厌氧微生境中定殖并参与CH₄的生物地球化学循环。甲基辅酶M还原酶（MCR）是产甲烷古菌特有的关键酶，催化产甲烷的最后一步反应，因此编码该酶α亚基的mcrA基因常被用作阐明产甲烷古菌群落丰度和活性特征的核心功能标记。相反，甲烷氧化主要由甲烷氧化菌介导，在有氧条件下将CH₄转化为CO₂、H₂O和生物量。甲烷单加氧酶（MMO）是催化该过程的关键酶，编码该酶的pmoA基因已成为表征甲烷氧化菌群落的重要分子标记。然而，pmoA基因丰度和MMO相关指标仅代表甲烷氧化的潜力，并不等同于原位甲烷氧化通量。

现有研究主要集中在植被恢复引起的土壤理化性质变化和细菌群落响应上，但很少有系统性分析关注超干旱沙漠，特别是塔克拉玛干沙漠的防护林建设如何影响土壤古菌群落和土壤微生物的甲烷循环功能潜力。这一知识空白制约了对沙漠植被恢复下与土壤甲烷循环相关的微生物响应模式的理解。本研究选择塔克拉玛干沙漠作为研究区域，从防护林地和未植被的裸沙区采集土壤样品，通过系统比较两种样地间的土壤理化性质、古菌群落组成以及与甲烷循环相关的关键指标，旨在实现以下目标：（i）刻画防护林建设下沙漠土壤古菌群落的结构和多样性特征；（ii）通过定量关键功能基因和相关功能蛋白的丰度，分析与CH₄循环相关的微生物关联特征；（iii）识别与古菌群落结构和CH₄循环相关指标显著关联的主要环境因子。本研究系统表征了沙漠防护林土壤中甲烷循环相关指标与微生物群落之间的关联，为评估防护林的生态效益以及揭示极端沙漠土壤中甲烷循环的响应方式提供了基础数据和参考。

二、关键技术方法

研究人员在新疆维吾尔自治区巴音郭楞蒙古自治州塔克拉玛干沙漠腹地的塔中地区设置了两个代表性样地：一个裸沙区和一个人工防护林地。两个样地均为砂质土壤。样品采集分为两个土层：表土层（0–10 cm）和亚表层（10–40 cm）。每个处理设置三个空间独立的1 m × 1 m样方作为生物学重复，采用五点法采集相应土层的土壤，混合后获得一个复合样品，共获得12个混合土壤样品。测定土壤理化性质的方法包括：电位法测定pH值，烘干称重法测定土壤含水量（SWC），干残渣法测定全盐量（TSS），外加热重铬酸钾氧化-硫酸亚铁铵滴定法测定土壤有机碳（SOC），紫外分光光度法和纳氏试剂比色法分别测定硝态氮（NO₃^?-N）和铵态氮（NH₄⁺-N）。微生物分析方面，使用高通量测序技术，针对古菌16S rRNA基因的V4–V5区进行扩增和测序，基于Silva138.2数据库进行分类。采用SYBR Green I实时荧光定量PCR（qPCR）对pmoA基因进行绝对定量。采用终点PCR扩增mcrA基因，但所有样品均未获得可测序的靶标条带。采用商业ELISA试剂盒定量MCR和MMO免疫反应性抗原蛋白的丰度。统计分析方法包括方差分析（ANOVA）、基于Bray–Curtis距离矩阵的ANOSIM检验、冗余分析（RDA）、皮尔逊相关性分析以及随机森林回归模型。

三、研究结果

3.1 土壤理化性质的差异

防护林增强了土壤水分、pH值和SOC的垂直分层。在裸沙区，表土层（0–10 cm）与亚表层（10–40 cm）之间无显著差异。而在防护林地块，不同土层表现出明显的特异性差异：表土pH值较低、SOC含量较高，而亚表层土壤含水量（SWC）更高。两类样地的土壤总体呈高碱性（pH > 8.0），S10的pH值最低（7.73）。S40的SWC最高（1.23%），较D40（0.13%）高出846.15%。S10的SOC含量最高（2.37 g/kg），较D10（0.96 g/kg）高出146.88%。各组间的TSS、NO₃^?-N和NH₄⁺-N无显著差异。

3.2 古菌群落组成与多样性

3.2.1 古菌群落组成

在门水平上，裸沙土和防护林土壤的古菌群落均以盐杆菌门（Halobacteriota）、热变形菌门（Thermoproteota）和热原体门（Thermoplasmatota）为主。裸沙土以嗜盐古菌为主，而防护林土壤中氨氧化和产甲烷谱系的相对丰度更高，尤其在亚表层。在属水平上，D10的优势类群为Candidatus Nitrosocosmicus（37.56%）、Marine Group II古菌（21.26%）等；D40主要为未分类的Haloferacaceae（28.69%）和Natronomonas（21.10%）。S10的特征是多样化的嗜盐古菌群落；S40的优势类群则转向氨氧化和产甲烷古菌，如Candidatus Nitrosocosmicus（16.04%）、Candidatus Nitrososphaera（8.20%）和Methanomassiliicoccales（8.13%）。

3.2.2 古菌群落多样性

Chao-1指数在各组间无显著差异，但防护林区域的数值略高于沙土区域。香农指数显示S10组最高（4.43），显著高于D10组（3.51）。辛普森指数显示S10组最低（0.02），显著低于S40组（0.08）。总体而言，防护林区域的表土（S10）拥有最高的多样性和均匀度。主坐标分析（PCoA）显示，裸沙区（D）与防护林（S）的古菌群落样品之间存在明显的空间分离，ANOSIM检验证实两组间群落组成差异显著（R = 0.6142，p = 0.003）。PC1和PC2共同解释了总变异的76.68%。

3.3 甲烷相关基因和功能蛋白丰度的变化

对所有土壤样品进行了mcrA基因的PCR扩增，琼脂糖凝胶电泳显示所有样品均未获得可用于进一步测序的mcrA靶标扩增子，表明研究区土壤中mcrA基因丰度低于该PCR系统的检测限。pmoA基因在所有样品中均有检出，防护林表土（S10）的pmoA基因拷贝数最高。S10和S40的pmoA拷贝数显著高于裸沙表土（D10）和裸沙亚表层（D40），相对增幅分别为2735.56%和220.58%，但其绝对丰度仅为2.2 × 10⁵ copies g^?1 dry soil，整体丰度仍然较低。MCR蛋白丰度在四组样品间无显著差异。对于MMO蛋白丰度，S10与其他三组差异显著；与D10和D40相比，S10和S40分别增加了1150.54%和220.57%。

3.4 环境因子与古菌群落及甲烷循环指标的关联

冗余分析（RDA）显示，RDA1和RDA2分别解释了古菌群落总变异的33.52%和16.57%，累计解释量为50.09%，表明环境因子与古菌群落组成之间存在显著关联（置换检验，p < 0.05）。裸沙土与较高的pH值、铵态氮和硝态氮相关，而与较低的SOC和SWC相关；防护林土壤则与较高的SOC和SWC相关。相关性分析显示，SOC与Arch_Shan（古菌群落香农多样性指数）、pmoA基因丰度和MMO蛋白丰度呈显著正相关；pH值与SOC、pmoA基因丰度和MMO蛋白丰度呈显著负相关。pmoA基因丰度与SOC和Arch_Shan均呈显著正相关。MCR蛋白丰度与TSS、pH、NO₃^?-N和SWC呈正相关，与NH₄⁺-N和SOC呈负相关。随机森林模型显示，在表土层（0–10 cm），pH值对MMO蛋白丰度的相对重要性最高（OOB R² = 0.691）；在亚表层（10–40 cm），SWC的相对重要性最高，而pH值的关联很弱（OOB R² = 0.036）。

四、讨论与结论

讨论部分总结：防护林建设可以显著改善土壤理化性质，特别是提高土壤水分和有机碳含量。这可能源于植物生长对土壤物理结构的改良、增强土壤持水能力以及促进土壤有机碳积累。裸沙土的古菌群落以寡营养、耐盐和嗜热的极端微生物为主，其中盐杆菌门占绝对优势，但不参与甲烷循环。防护林土壤中与甲烷代谢相关的古菌类群包括norank_c__Thermoplasmata和norank_o__Methanomassiliicoccales，但仅有后者被认为具有产甲烷潜力。然而，ELISA结果显示MCR蛋白丰度普遍较低且组间无显著差异，终点PCR也未在任何土壤样品中检测到mcrA基因，表明研究区土壤可能缺乏功能活跃的产甲烷种群，或其丰度低于检测限。防护林表土中pmoA基因拷贝数和MMO蛋白丰度均显著升高，表明防护林建设改变了土壤理化性质，在表土形成了相对稳定的好氧-厌氧界面，有利于甲烷氧化微生物。相关性分析和RDA结果表明，SOC和pH是与目标群落和功能指标关联最强的环境因子。随机森林模型揭示了这些关联的土层特异性：表土中pH值与MMO蛋白丰度关联最强，而亚表层中SWC的关联最强。这些发现表明，环境条件、底物可用性和微生物群落结构与不同土层的功能指标均有关联。研究也存在局限性，包括恶劣野外条件限制了重复数、两种样地背景条件的潜在混杂效应、小样本量导致模型不确定性、pmoA qPCR未进行回收率校正、缺乏土壤甲烷通量实测数据等。

研究结论：研究人员从塔克拉玛干沙漠的裸沙和防护林地采集土壤样品，考察了防护林建设后古菌群落结构和土壤甲烷氧化相关功能指标的变化，并鉴定了与这些甲烷氧化相关指标强关联的环境因子。结果表明，防护林建设改善了沙漠土壤的理化性质，提高了表土有机碳浓度，降低了土壤pH值，并增加了土壤湿度。这些环境变化与古菌群落多样性的提升以及甲烷氧化相关功能指标值的升高相关联。环境变量与这些功能指标之间的关联呈现土层分层特征：表土中的指标与pH值关联更强，而亚表层指标则与土壤含水量和有机碳共同关联。尽管防护林的碳和水输入有利于甲烷氧化相关功能指标的升高，但塔克拉玛干沙漠总体上仍是一个寡营养系统，所测指标仍处于干旱生态系统的典型范围内。这些发现有助于阐明沙漠人工植被重建与土壤古菌甲烷循环相关功能之间的潜在关联，并为沙漠公路防护林的生态恢复和干旱区碳循环研究提供参考。