集约经营对毛竹林土壤反硝化细菌丰度的影响*

曹林桦，刘彩霞，刘 茗，方 伟，梁辰飞，秦 华，陈俊辉，徐秋芳†

（1. 浙江省森林生态系统碳循环与固碳减排重点实验室，浙江农林大学，浙江临安 311300；2. 浙江农林大学环境与资源学院，浙江临安311300）

反硝化作用是土壤氮素循环过程中的重要一环，是指反硝化微生物在厌氧的条件下，通过一系列的酶催化还原反应，先将转化为，然后再进一步转化为气态氮（NO，N2O，N2）的过程[1]。反硝化作用不仅会导致氮肥损失，造成严重的经济危害。同时反硝化作用过程释放的N2O 是三大温室气体之一，其产生的温室效应约是CO2的300 倍[2]。此外，N2O 还会对臭氧层造成极大的破坏，增加皮肤癌的患病概率[3]。虽然硝化和反硝化过程均能产生N2O。但与硝化过程不同的是，N2O 是反硝化过程中的必要产物，甚至是终产物[4]，而且有大量研究表明亚热带酸性森林土壤中由反硝化作用产生的N2O 大于硝化作用[5]。虽然反硝化微生物主要分布在水体、沉积物和土壤等环境[6]，但由于土壤的环境更适合微生物生存，微生物多样性更为丰富，而且土壤中氮含量普遍较水体、沉积物高，因此，研究土壤微生物反硝化作用对于减少氮肥损失，降低N2O 排放均有着重要的意义。

反硝化作用是一个由四步反应所构成的生物地球化学过程，包括还原、还原、NO 还原和N2O 还原。其每个步骤均由相应细菌的代谢合成酶催化完成。

毛竹林土壤属于旱地土壤，虽然总体为好氧环境，但土壤中存在大小不等的孔隙，微小孔隙的局部环境也可能是厌氧的。故在此种土壤环境中也会有相应的反硝化过程存在。然而已有的大多数研究均围绕着不同树种森林[8]、沙漠[9]、草原[10]等其他生态系统展开。有关毛竹林反硝化过程微生物的变化情况少见报道。

毛竹是我国重要的经济和生态竹种，具有生长周期短、产量高、经济价值高等特点。为追求更高的经济效益，人们对原来以粗放经营为主的毛竹林进行了集约化经营。集约经营的主要措施有施肥、去除林下植被、翻耕等。这种集约经营模式虽然提高了毛竹的产量，但人为耕作导致水土流失加剧，使土壤出现不同程度的退化[11]。与此同时，长期集约经营降低了土壤微生物的活性[12]，这也促使人们加大了对集约经营过程中土壤微生物相关生态变化的关注[13-15]。然而，与碳氮循环密切相关的反硝化微生物的研究少见报道。许多研究证明氮肥施用后因硝化作用产生大量N2O[16]，集约经营毛竹林土壤是否存在反硝化微生物？是否也会导致排放N2O 等问题值得探索。本文以nirK、nirS和nosZ3 个基因为研究对象，运用定量PCR 的方法初步探究集约经营过程中毛竹林土壤3 种功能基因丰度的变化，研究结果对毛竹林经营管理具有重要指导意义。

1 材料与方法

1.1 研究区概况

研究区遂昌县位于浙江省西南部（28°13′—18°49′N，118°41′—119°30′E）。属中亚热带季风类型，全年平均气温16.8 ℃，年降水量1 510 mm，降水天数172 d，年无霜期251 d。土壤母岩为由花岗岩变质而成的片麻岩石，土层厚度大于 1 m，坡度为20°～25°。毛竹密度约为 1 800 株·hm-2，1989 年之前为竹林开发阶段，仅开垦、劈山，不施肥。该阶段毛竹林下主要有 9 种植被，其中狗脊蕨（Woodwardia japonica） 和 兔 儿 伞 （Syneilesis aconitifolia）盖度为 90%。1989—2000 年施肥方法为：大年（留养新竹年）6 月在每株毛竹的上方开环形沟，施用尿素（或碳铵）450 kg·hm-2；2000—2010 年施用配方肥，每年用量为尿素450 kg·hm-2、过磷酸钙 380 kg·hm-2、氯化钾 75 kg·hm-2；2010 年之后每年施用毛竹专用肥（福建中化生产）750 kg·hm-2和尿素 450 kg·hm-2。集约经营过程毛竹林下基本无灌木杂草，毛竹平均密度约3 250 株·hm-2。

1.2 样品采集与处理

于2013 年9 月按照生态控制原则，选择生长历史一致但集约经营时间不同的毛竹林，采用随机取样法采集不同经营年限毛竹林表层（0～20 cm）和亚表层（20～40 cm）的土壤样品。该毛竹林地均位于同一山头，并由同一农户经营，采取的施肥管理措施也均一致。集约经营的时间分别为：0 a（对照组）、10 a、15 a、20 a 和25 a。每个年份分别选择3个样地，即3 个重复，每个样地采用S 型取样法选取5 个样品点就地混合为一个样品，分别采集表层（0～20 cm）、亚表层（20～40 cm）的土样，过筛（2 mm）后装入密封袋，将样品放入冰盒带回实验室。样品共分为2 份，1 份放入-70℃冰箱，经冷冻干燥后，提取土壤细菌总DNA，用于分子生物学分析；另一份于室内自然风干，充分研磨过筛后用于土壤基本理化性质的测定。

1.3 分析方法

土壤理化性质分析：土壤理化性质分析方法参考《土壤农化分析》[17]。土壤pH 采用酸度计测定，土水比为 2.5︰1；有机碳含量采用重铬酸钾-浓硫酸油浴外加热法测定；土壤全氮含量采用半微量凯式滴定法测定；碱解氮含量采用碱解扩散法测定；有效磷含量采用盐酸-氟化铵溶液浸提，钼锑抗比色法测定；速效钾含量采用醋酸铵提取，火焰光度计测定；硝态氮含量采用氯化钾浸提，紫外分光光度计测定；铵态氮含量采用氯化钾浸提，靛酚蓝比色法测定。

土壤总 DNA 提取：采用 Power Soil DNA Isolation Kit 试剂盒提取土壤总DNA。称取0.50 g 保存于-70 ℃冰箱，经冷冻干燥后，按试剂盒说明书进行提取。DNA 提取成功后经1%（m/v）的琼脂糖凝胶电泳检测DNA 片段大小，并用微量分光光度计检测其浓度和纯度。提取后的DNA 样品保存于-40℃。

土壤反硝化细菌质粒制备：以土壤样品 DNA为模板分别对nirK、nirS、nosZ三种功能基因进行PCR 扩增，扩增条件如表1 所示。PCR 产物用纯化试剂盒（Takara）纯化后与PEASY-T3 Vector（全式金生物技术公司，北京）连接，然后将连接产物转化到Trans1-T1 感受态细胞中，并涂在含有X-gal、IPTG 和 Amp 的 LB 琼脂平板上进行蓝白斑筛选，用菌落PCR 方法检测阳性克隆子。对阳性克隆子用LB 液体培养基继续培养，最后用质粒抽提试剂盒提取质粒。

表1 聚合酶链式反应中的引物及反应条件Table 1 Primers and PCR condition used for the PCR amplification

土壤反硝化细菌相关功能基因荧光定量PCR：采用引物 FlaCu：R3Cu、Cd3aF：R3cd、Lb：Rb 分别对nirK、nirS、nosZ三种功能基因进行基因片段扩增，定量 PCR 在荧光实时定量 PCR 仪（Bio-Red,USA）上进行。反应体系为 20 μL，反应体系如下：2×Premix Ex Taq 10 μL，50 μmoL·L-1上游和下游引物各0.2 μL，模板DNA1.0 μL，无菌双蒸水8.6 μL。三种功能基因荧光定量PCR 条件相同，均为:95℃预变性3 min，95℃变性10 s，56 ℃退火30 s，72℃延伸 20 s，35 个循环，溶解曲线过程（65℃至95℃，0.5℃·5 s-1）。标准曲线制作：三种功能基因质粒 DNA 经 Nanodrop®ND-1000 测定，nirK、nirS和nosZ基因重组质粒浓度分别为105.6、102.5、135.2 ng·μL-1，计算其拷贝数分别为 2.74×1010、2.70×1010、3.74×1010copies·μL-1，以 10 倍梯度对重组质粒进行梯度稀释（ 10-3～10-8），每个梯度3 次重复，nirK、nirS和nosZ决定系数（R2）分别为 0.995、0.998 和 0.999。

1.4 数据处理

采用 Microsoft Excel 2016 软件对数据进行处理，采用Origin 7.5 软件作图，采用SPSS 21.0 软件对数据进行分析，采用单因素方差分析（One-way ANOVA）和最小显著差数（LSD）法进行显著性和多重比较（α=0.05），采用皮尔逊（Pearson）法进行相关性分析。通过Canoco 4.5 结合环境因子分析微生物功能基因丰度的变化。

2 结 果

2.1 集约经营毛竹林土壤理化性质

有关毛竹集约经营土壤理化性质变化已有报道[15]。简要概括：毛竹集约经营过程中土壤理化性质变化，表层土壤 pH 呈下降趋势；有机碳、全氮含量于经营10 a 达到最大（P＜0.05），其余年份无显著差异；碱解氮呈现先上升再下降后又上升的波动变化规律；碳氮比经营初期缓慢上升，15 a 时达到最大，20 a 显著下降（P＜0.05）；有效磷和速效钾总体呈上升趋势；硝态氮出现先减少后增加的变化规律，而铵态氮在整个经营过程中变化不大。

与表层土壤不同的是，集约经营前期亚表层毛竹土壤pH 没有发生明显变化，而在第25 年时显著下降 0.42 个 pH 单位（P＜0.05）；全氮和碱解氮均在经营25 a 之后显著增加（P＜0.05）；有机碳含量在经营10 a 和25 a 显著增加（P＜0.05），而在经营20 a 显著下降（P＜0.05）。其余理化因子变化与表层相同。

2.2 集约经营对毛竹林土壤反硝化细菌功能基因丰度的影响

采用荧光实时定量PCR分析集约经营毛竹林土壤反硝化细菌丰度结果表明，nirK基因数量在2.62×107～1.28×108copies·g-1干土之间，随着集约经营时间的持续，表层土壤nirK基因数量呈现周期性波动态势（图1A），上升过程发生在集约经营10 a和20 a 的2 个时间，其基因丰度显著大于其他年份（P＜0.05），集约经营15 年时丰度达到最低值，显著小于其他年份（P＜0.05），第 25 年时的基因丰度与对照组无显著性差异。而亚表层土壤的nirK丰度变化呈V 字型，集约经营15 a 时才开始下降、20 a 时达到低谷，土壤nirK基因数量在15 a 和20 a 分别较对照下降43%和63%（P＜0.05），25 a 时又恢复至对照水平，集约经营10 a 和25 a 与对照无显著性差异。nirS基因数量在 1.45×106～4.96×107copies·g-1干土之间，总体低于nirK基因。表层土壤基因数量在集约经营 10 a 显著上升 41%（P＜0.05），之后一直保持较高的水平（图 1B），说明集约经营促进表层土壤中携带nirS基因微生物的繁殖。土壤亚表层nirS基因数量则在前20 a 的集约经营过程中呈下降趋势，第20 年最大下降幅度为63%，而第25 年则恢复到对照水平。nosZ基因数量在 3.81×107～3.03×108copies·g-1干土之间。表层土壤基因丰度在集约经营10 a 明显增加，第15 年时开始下降，集约经营 20 a 较对照低 70%，之后则变化不大（P＞0.05）。亚表层土壤nosZ基因丰度前 10 年变化不大（图1C），第15 年和20 年分别降低34%和48%，显著低于对照和 10 a 土壤（P＜0.05），第 25 年时nosZ基因丰度恢复、且数量显著高于对照（增加35%，P＜0.05）。集约经营过程中表层土壤的（nirK+nirS）/nosZ比率为 0.34～1.08，在亚表层为1.63～2.27。表层（nirK+nirS）/nosZ的比率在集约经营20 a 显著增加，而亚表层中显著降低。

图1 长期集约经营毛竹林土壤反硝化细菌nirK、nirS、nosZ 功能基因丰度Fig. 1 Abundance of nirK，nirS，nosZ genes in the soils of the moso bamboo plantations under long-term intensive management

2.3 土壤理化性质对反硝化细菌功能基因丰度的影响

集约经营主要通过土壤性质的改变来影响土壤微生物。为揭示土壤性质对反硝化细菌丰度的影响，将表层与亚表层土壤的三个功能基因的丰度分别与对应的土壤理化性质进行相关性分析，结果如表 2所示。在表层土壤中，nirK基因数量与全氮呈显著正相关（P＜0.05），nirS基因与有效磷呈显著正相关（P＜0.05），而nosZ基因则与有机碳和碳氮比呈显著正相关（P＜0.05）、与速效钾和硝态氮含量呈显著负相关（P＜0.05），说明影响表层土壤三个功能基因丰度的主要因素并不一致。作为 N2O 排放量指标[21]的（nirK+nirS）/nosZ比值与碳氮比呈显著负相关（P＜0.05）、与速效钾和硝态氮呈显著正相关（P＜0.05）。在亚表层土壤中，三种功能基因均与有机碳呈显著正相关（P＜0.05），nosZ基因还与全氮和碱解氮含量呈显著正相关（P＜0.05）。（nirK+nirS）/nosZ的比值与硝态氮呈显著负相关（P＜0.05）。

表2 毛竹林土壤反硝化功能基因与环境参数之间的皮尔逊相关系数Table 2 Pearson’s correlations（r）between denitrifying functional gene abundance and environmental parameters in the moso bamboo plantation

为了进一步揭示集约经营过程中土壤性质对反硝化功能基因丰度变化的影响，将表层以及亚表层土壤的反硝化功能基因丰度分别与对应土层的理化性质作为两组变量进行冗余分析（RDA，Redundancy analysis）。根据蒙特卡洛检验，选取影响较大的 4个环境因子分析结果表明，表层土壤中第一排序轴和第二排序轴分别解释了 25.8%和 23.6%的变异（图 2A），第一排序轴将 CK 和10 年样品与20 年和25 年样品很好地区分开来，而 15 年样品则在第二坐标轴附近，说明集约经营过程中土壤反硝化功能基因丰度变化分为三个阶段。第一排序轴与碳氮比呈正相关、与硝态氮和速效钾呈负相关，第二排序轴与全氮呈负相关。土壤硝态氮（F=3.855，P=0.024）、碳氮比（F=3.832，P=0.022）、速效钾（F=2.844，P=0.050）对土壤反硝化细菌功能基因丰度的影响均达显著水平，而全氮（F=2.646，P=0.054）则接近显著水平。进一步分析可得，硝态氮和速效钾对 20 年和 25 年样地反硝化细菌影响最大。nirK与全氮呈显著正相关，nosZ与硝态氮和速效钾呈显著负相关，这与相关性分析结果相同。

图2 集约经营毛竹林土壤反硝化功能基因丰度与环境参数的冗余分析Fig. 2 Redundancy analysis of abundance of denitrifying functional genes and environmental parameters in soils of the moso bamboo plantation under intensive management

根据蒙特卡洛检验，选取亚表层土壤中影响较大的4 个环境因子分析结果表明，第一排序轴和第二排序轴分别解释了66%和1.7%的变异，第一排序轴上三个功能基因与有机碳相关性很强（图 2B）。不同集约经营时间的土壤均聚焦在一起（一个对照重复外），但分组情况与表层土壤不同，15 年和20年聚在一组、位于第一排序轴正值区，其他3 个年份组位于第一排序轴的负值区，且25 年又与CK 和10 年在第二排序轴分开。说明集约经营时间对3 个反硝化功能基因影响显著，其中影响最大是有机碳和铵态氮，与第一排序轴呈负相关，有机碳（F=14.717，P=0.002）达到显著水平，而铵态氮（F=2.795，P=0.092）则没有；其次为碱解氮和碳氮比，分别与第二排序轴呈正相关和负相关，但均未达到显著水平（碱解氮 F=2.315，P=0.114、碳氮比F=1.811，P=0.181）。进一步分析发现，碱解氮对集约经营25 年样地有积极影响，有机质与3 种功能基因均呈显著正相关，全氮和碱解氮与nosZ功能基因呈正相关。

3 讨 论

前人的研究结果证明，影响土壤反硝化作用的因素比较多，而且各因素之间互为影响，主要影响因子包括速效氮、有机碳、土壤水分、温度、pH 以及土壤质地等[22]。理论推测以及研究证明，单因素的速效氮[23]、有机碳[24]、土壤水分[23]、温度[25]、pH[26]与反硝化作用均有正相关关系，而自然或人为耕作土壤的土壤性质变化不一致（不同步），因此，实际测量得到的反硝化基因丰度变化是多个因素综合作用的结果。由于不同因素的作用强度不同，正反作用可能会抵消、或者作用强的因子将弱的因子掩盖掉。

3.1 集约经营对毛竹林表层和亚表层反硝化基因的总体影响

毛竹林土壤的反硝化基因nirK、nirS和nosZ丰度总体较报道的结果高（图1），除旱地农田nirS外，水稻田、杉木林和旱地农田中的三种功能基因丰度均显著低于毛竹林或与毛竹林持平[27-29]。这表明毛竹林土壤含有更高的反硝化基因丰度，因此相比其他土地利用方式可能N2O 排放量较高。集约经营措施主要措施有施肥、去除林下植被、翻耕等。由此，集约经营过程中土壤速效养分水平提高，土壤通气性增加，pH 下降，有机碳含量不断波动。其中 pH下降的原因在于集约经营过程中大量氮肥的施用，而有机碳含量不断变化的原因在于一方面集约经营初期施肥结合翻耕促进土壤中地下根鞭的分解，另一方面施肥促进毛竹生长，增加毛竹根系分泌物量，从而提高了土壤的有机质含量。而10 年后不断翻耕施肥导致有机碳含量下降。25 年有机质又有小幅回升，这可能与毛竹生产量提高后，枯落物和根系分泌物的增加有关。随着集约经营时间的增加，三种反硝化基因的变化表现出一定的规律性，无论是表层还是亚表层土壤，三种功能基因丰度在集约经营10 年时增加或保持不变，说明在前 10 年期间毛竹林集约经营措施对反硝化细菌没有产生不良影响，有的还有促进作用。尽管有研究表明酸性范围内土壤 pH 越低，反硝化细菌丰度也会随之降低[26]。而在本研究中有机碳、全氮和碱解氮明显增加对反硝化作用的正作用抵消了 pH 下降的负面影响，从而掩盖了 pH 下降带来的不利影响。土壤硝态氮下降是反硝化作用增强导致底物减少的结果（表1）。这与斐自伟等[30]的研究结果相似。除表层土壤nirS基因外，三种功能基因在集约经营过程的某一阶段出现不同程度的下降，但集约经营持续25 年时，除了nosZ表层功能基因外，其他反硝化功能基因丰度均恢复到对照甚至超过对照的水平（图1）。10 年之后土壤有机碳及pH 明显低于对照和10 年土壤，这可能是导致反硝化基因丰度下降的主要原因。而在25年反硝化基因丰度回升可能是由于有机碳含量的上升。由于土壤因子变化较复杂，将土壤性质与基因丰度进行相关性和冗余分析以揭示内在规律，相关分析表明，nirK和nirS基因丰度与不同土层土壤性质的相关性没有一致的规律，而nosZ基因丰度在表层和亚表层均与有机碳显著相关。表层和亚表层土壤pH 和铵态氮与所有基因均没有显著相关性（表2），其他指标则分别或同时与某一土层的某个基因表现出显著相关性，其中单因子有机碳、全氮以及它们的比值（C/N）与反硝化基因相关的频率最高。冗余分析表明，表层土壤的硝态氮、C/N、速效钾以及全氮含量等4 个因子对3 个反硝化基因影响的显著水平依次递减，而对亚表层影响较大的4 个因子依次为：有机碳、铵态氮、碱解氮以及 C/N，不同土层驱动基因丰度变化的主要土壤性质不同，而C/N 同时影响两个土层，且两个土层分别同时出现2 个含氮指标，说明集约经营措施主要通过土壤氮和有机碳综合影响反硝化微生物的活动和功能。由于所有土壤的质地相同，集约经营措施是每年去除林下植被、翻耕，故不同集约经营时间的土壤疏松度差别不大，因此，土壤化学性质变化是影响基因丰度的主要因素。

3.2 集约经营过程中表层和亚表层反硝化基因的差异

为弄清土层深度对不同反硝化基因的影响，我们对不同土层之间各种基因丰度差异情况展开讨论。表层nirS和nosZ土壤基因丰度显著高于亚表层土壤，但nirK基因丰度出现相反现象。大部分研究报道表层土壤反硝化基因丰度高于亚表层土壤，nirK、nirS和nosZ丰度在不同类型土壤如耕地[31]、草地[32]和森林[33]中均随着深度的增加而降低。虽然亚表层土壤的氧化还原电位通常较表层低，有利于反硝化作用，但表层土壤有机碳高于亚表层，可能对反硝化作用的影响更大。为什么亚表层土壤的nirK基因丰度反而高于表层的异常现象？从前面分析得知，毛竹林集约经营土壤反硝化细菌丰度变化受到氮素水平和有机碳的影响最大，而亚表层土壤的氮素水平和有机碳均低于正常表层，从土壤性质解释不通。推测可能是分布在30 cm 亚表层的毛竹竹鞭分泌大量可溶性有机碳，对携带nirK基因的微生物促进效果较nirS和nosZ基因微生物更加显著。李振高等[34]也发现在小麦生长过程中随着根系分泌物的增加反硝化细菌数量不仅显著增加，同时根系分泌物还刺激了反硝化细菌的代谢活动。此外，研究证明pH 与nirK反硝化基因丰度呈显著正相关[26]，而在本研究中表层pH 显著低于亚表层，因此pH 也可能是导致nirK表层和亚表层丰度出现异常的原因之一。

3.3 集约经营过程中不同反硝化基因之间的差异

本研究中，反硝化基因nirK、nirS和nosZ等随着集约经营时间的变化规律不同。尽管表层土壤中nirK和nirS两种功能基因共同编码亚硝酸还原酶，且丰度上有着相同的数量级，但毛竹集约经营过程中两种nir基因变化趋势不同，且分别受不同的土壤理化因子的控制。关于这两种基因的同步[35]和反向[36]变化规律的结果均有报道。由于土壤中携带反硝化nirK和nirS的微生物属于不同的种[7]，在环境中有各自的生态位[27]，对同一环境因子的响应可能不同。有时可能存在竞争关系[36]、有时候存在正相关关系[35]，有时则互不影响[26]，三种关系均有可能发生。毛竹林表层土壤中nirK和nirS无显著相关性，属于互不影响的关系，这可能是因为nirK丰度在整个集约过程中的波动远大于nirS。研究表明nirK对环境的敏感程度要远大于nirS[37]。而集约经营过程中由于施肥、毛竹凋落物、翻耕等抗干扰活动导致微生物的生境不断变化，可以很好地解释nirK基因丰度波动大于nirS的结果。与nirK和nirS变化规律不同，集约经营过程中表层和亚表层土壤的nosZ基因丰度均表现为先增加后减少的一致规律，且与有机碳含量呈显著正相关；表层土壤nosZ丰度与速效钾和硝态氮含量呈显著负相关、而亚表层土壤nosZ丰度则与碳氮比呈显著正相关，说明有机碳是影响nosZ丰度的共性因子，研究表明[38]随着氮肥施用量增加，土壤碳氮比降低，土壤中反硝化微生物nosZ的丰度显著降低。毛竹林在集约经营初期由于施肥促进有机质分解以及促进毛竹生长使毛竹分泌根系分泌物增加，因而，土壤有机质含量明显增加，虽然施肥增加土壤氮水平，但增幅小于有机碳，因此，土壤碳氮比总体升高有利于携带nosZ基因的微生物活动。而长期集约经营碳氮比显著低于对照，nosZ基因丰度同碳氮比一样呈现同步变化的趋势。

3.4 集约经营过程中（nirK+nirS）/nosZ 比值变化

研究发现（nirK+nirS）/nosZ比值可预测 N2O排放的变化，相同基因丰度水平情况下比值越高产生的 N2O 越多[21]。表层土壤（nirK+nirS）/nosZ比值在集约经营初期变幅不大，而在集约经营后期显著上升（图 1D），说明随着集约经营时间可能会增加N2O 排放。表层后期nosZ基因降幅较大（图1C）而nirS基因则没有下降（图 1B），因此造成（nirK+nirS）/nosZ比值明显增加。集约经营后期硝态氮明显积累，为反硝化作用提供丰富的底物，反硝过程可以保持在较高水平，而nosZ基因数量显著减少使硝化过程不能彻底还原为N2，导致N2O 排放量增加。nosZ基因数量减少主要受到与硝态氮积累有关，研究表明[23]高浓度的硝态氮会抑制N2O 还原酶（nosZ）的活性，使 N2O 转化为 N2速率减慢。亚表层土壤由于nirK明显高于表层、而nosZ则明显低于表层，导致（nirK+nirS）/nosZ比值显著高于表层（图 1D），但随着集约经营时间呈缓慢下降趋势。从3 个基因丰度变化情况看，亚表层土壤对nosZ基因的影响大于（nirK+nirS），说明亚表层土壤反硝化步骤主要为还原形成NO，因而产生较多的N2O。然而，亚表层土壤（nirK+nirS）/nosZ的比值在集约经营后期显著降低，这正好与表层的变化规律相反、呈显著负相关。亚表层土壤（nirK+nirS）/nosZ与硝态氮浓度呈显著负相关，这可能是由于土壤亚表层相比表层硝态氮浓度更低，故在集约经营后期并未抑制N2O 还原酶活性。以上结果推测，整个集约经营过程亚表层土壤反硝化产生的N2O 均高于表层，表层土壤集约经营后期N2O 排放量增加。

4 结 论

反硝化作用是竹林土壤氮循环的重要一环。随着集约经营时间的增加，两层土壤的三种功能基因丰度在集约经营 10 年时表现增加或不变的一致规律。除表层土壤nirS基因外，三种功能基因在集约经营过程的某一阶段（15 或 20 年）出现不同程度的下降，集约经营持续25 年时，除nosZ表层土壤基因丰度仍低于对照，其他不同土层土壤基因丰度均恢复甚至超过对照水平，说明土壤反硝化细菌表现出对集约经营干扰的抵抗和恢复反应。集约经营措施主要通过土壤氮和有机碳变化综合影响反硝化细菌的活动和功能。毛竹林集约经营土壤的反硝化细菌积极参与氮循环过程，加剧N2O 温室气体的排放。建议在集约经营过程中，通过施用缓释肥料以减少反硝化细菌的底物硝酸根。