反硝化作用是土壤氮素流失的关键途径之一,该过程由多种微生物协同完成,涉及一系列关键功能基因
[1],不仅加剧全球气候变暖,还对臭氧层具有显著的破坏作用
[2]。盐渍农田是我国重要的后备耕地资源,已有研究
[3-5]多集中于土壤改良,而对其中关键氮循环过程-反硝化作用的研究仍较为薄弱
[6],盐渍农田中氮肥利用率普遍低于常规农田
[7],高氮输入,一方面抑制固氮微生物功能基因
nifH表达
[8];另一方面,为反硝化过程提供过量底物(
),从而加速N
2O排放
[9],导致氮素气态损失大幅增加。此外,高施氮肥常因离子拮抗,降低植物体内GS/GOGAT途径对
的同化速率
[10],导致土壤中
累积;累积的
经硝化作用生成大量
,而
不易被土壤吸附,易随水淋失,进而导致地下水或地表水的面源污染。白芳芳等
[11]研究表明,减氮处理可显著降低农田土壤中反硝化微生物的功能基因丰度,并显著下调编码亚硝酸盐还原酶的关键功能基因(如
nirK和
nirS)及编码氧化亚氮还原酶的基因(
nosZ)的丰度
[12]。适度减氮不仅抑制微生物反硝化作用,还可激活植物的抗逆性,促进根系分泌柠檬酸、苹果酸等有机酸
[13]。该类有机酸既可螯合Na⁺离子缓解盐胁迫下的渗透胁迫,也能释放被碳酸钙固定的Ca
2⁺,提高土壤中钙的有效性
[14-15]。已有研究
[16-18]多聚焦于湿地、高原湖泊、河口沉积物等典型生态系统中的反硝化作用,然而在盐渍逆境作用下,根际土壤细菌反硝化途径与根系生长及产量形成对不同施氮量的响应尚不明确。为此,本研究以轻中度盐渍向日葵根际土壤为研究对象,基于田间小区定位试验与宏基因组测序技术,探究不同施氮量对盐渍农田根际土壤细菌反硝化途径、根系生长及产量的影响,以期为盐渍农田向日葵氮肥合理施用提供理论依据。
1 材料与方法
1.1 试验区概况
2023年和2024年5—9月于内蒙古农业大学海流图科技园区农业水土工程试验基地(40°41'30''N,111°02'30''E)开展田间小区试验,试验基地位于土默川平原灌区,属典型温带大陆性季风气候,多年平均降水量为421.6 mm,多年平均气温6.3 ℃,无霜期135 d。2023年向日葵生育期内(5—9月)平均气温21.2 ℃,年日照时间为2 861.3 h,生育期降水量314.6 mm,蒸发量为1 868.1 mm。2024年向日葵生育期平均气温22.9 ℃,年日照时间为2 746.2 h,生育期降水量589.9 mm,蒸发量为1 842.2 mm。试验地0~40 cm土层土壤质地为砂质壤土,体积质量(BD)为1.73 g/cm3,土壤平均全盐量为3.22 g/kg,属轻中度盐渍土。土壤表层基础养分有机质(SOM)为16.79 g/kg,全氮(TN)为990.3 mg/kg,水解性氮(HN)为63.51 mg/kg,有效磷(AP)为43.66 mg/kg,速效钾(AK)为216.38 mg/kg,向日葵整个生育期地下水埋深为3~5 m。
1.2 试验设计
供试作物为食用向日葵“HZ2 399”,播种前对盐渍农田进行春灌压盐(灌水定额900 m
3/hm
2),春灌采用畦灌,向日葵生育期采用膜下滴灌,底肥施用磷酸二胺(含N 18%),底肥折合纯氮50 kg/hm
2,现蕾期追施尿素(含N 46%)。以当地施肥量为对照CK(追施纯氮150 kg/hm
2),另设3个减氮梯度,减氮25%(N1∶追氮100 kg/hm
2)、减氮50%(N2∶追氮50 kg/hm
2)、减氮75%(N3∶0 kg/hm
2),该设计采用养分管理研究中的经典梯度设置方法,即通过设定基于常规施肥量的百分比减量梯度,旨在系统探究盐渍农田中作物-土壤系统对氮输入减少的非线性响应关系。试验采用单因素完全随机区组设计,小区规格为10 m×16 m,每个处理设置3次重复,共计15个小区,小区四周砌筑水泥挡板,以防止水、氮、盐的交互影响。每个小区铺设6膜,每膜种植2行,行距50 cm,株距50 cm,种植密度21 000株/hm
2,各小区种植密度一致,具体试验设计见
表1。
1.3 测定项目与方法
1.3.1 样品采集
试验采用五点取样法,灌水后分别于向日葵关键生育期使用根钻法取出根系并抖落根系附着土壤。5个点土样混和均匀并装入无菌自封袋并分为3份:1份置于-80 ℃下保存用于DNA提取和宏基因组测序;1份室内自然风干,研磨后过2 mm筛,用于土壤理化性质测试分析;1份用于根系分泌物的测定。使用pH计(梅特勒-托利多仪器,上海)测定土水质量比为1∶2.5的土壤酸碱度(pH),土壤铵态氮使用KCl(2 mol/L)浸提靛酚蓝比色法,土壤硝态氮按水土比10∶1加入50 mL KCl(1 mg/L)溶液,采用紫外分光光度计测定。土壤有机质(SOM)含量测定参照NY/T 1121.6—2006。根际有机酸获取采用乙醇浸提法,使用GC-MS联用仪(HP6890/5975C,美国)进行分析。
1.3.2 宏基因组测序及生物信息
称取500 mg glass beads和0.2~0.5 g的土壤,放入2 mL离心管并添加Buffer SLX Mlus进行振荡,加入80 μL Buffer DS孵育10 min,再离心5 min。转移600 μL上清液至2 mL离心管中,加入200 μL Buffer SP 2,振荡后加入100 μL HTR Reagent,冰浴5 min后离心5 min,移上清液400 μL至离心管中,加入40 μL的磁珠,混合均匀加入450 μL Binding Buffe,室温放置2 min。使用Qubit 4荧光计检测DNA的浓度,并用1%的琼脂糖凝胶电泳检测DNA的质量。调整DNA溶液浓度。多样品DNA文库均一化至10 nM后等体积混合,将混合好的文库(10 nM)稀释至合适的浓度后进行上机测序。
1.3.3 向日葵根系指标与产量测定
使用根钻采用五点取样法进行根系样本采集,将取得的根系样本用清水清洗干净后,采用根系分析系统(WinRHIZO),精准获取根系长度、根系表面积及根系体积等关键数据,进一步通过运算推导出根长密度(RLD)、根体积密度(RVD)及根表面积密度(RSAD)等参数。在向日葵成熟期时,在每个小区随机选取10株长势较为一致的向日葵,各小区分别测产。并对各小区定株的向日葵花盘上,选取最外侧3~4圈的位置摘取籽粒,测百粒重。
1.4 数据分析
采用Excel 2016软件进行数据的处理,使用SPSS.27软件进行单因素方差分析及差异显著性水平(p<0.05)分析,用最小显著差异法进行显著性检验。用Origin 2023软件绘制反硝化功能基因相对丰度图、门水平群落相对丰度图与根系指标图。文中图表数据均为2 a试验结果的平均值,用以体现其共性规律;2 a试验结果不一致之处则单独说明。
2 结果与分析
2.1 施氮对根际土壤反硝化细菌群落组成及功能基因的影响
施氮显著改变向日葵根际土壤反硝化细菌的群落组成(
图1)。与CK处理相比,N2与N3处理显著降低根际土壤中
norB、nirK型Proteobacteria的相对丰度,同时显著降低
narG、
nirS型Acidobacteria相对丰度。值得关注的是,N1与N2处理在调控微生物群落结构方面表现出积极效应,2023年N1处理显著提高
nosZ型Proteobacteria(30.03%),2024年N1与N2处理显著提高
nosZ型Proteobacteria(24.13%、47.20%)与
nosZ型Acidobacteria(40.32%、11.16%)相对丰度,该群落结构得到优化,尤其是完全反硝化细菌群落的富集,有助于推动根际土壤细菌介导的完全反硝化途径(N
2O→N
2),从而降低N
2O的排放潜力。
由
图2可知,施氮处理下向日葵根际土壤反硝化功能基因相对丰度。随着施氮量的减少,功能基因
nirK、
narG、
nirS等基因的相对丰度呈逐步下降趋势。然而施氮处理下编码氧化亚氮还原酶的
nosZ基因相对丰度由大到小排序为N1>CK>N2>N3,表明减氮25%(N1处理)在抑制根际土壤细菌不完全反硝化功能基因(
narG、
nirK、
nirS)表达的同时,促进完全反硝化功能基因(
nosZ)的富集,说明减氮25%处理不仅从源头阻控根际土壤氮素损失,还促进N
2O向N
2还原,降低盐渍农田氮素损失的同时减少温室气体N
2O的排放潜力。
2.2 施氮对向日葵根系生长与产量形成的影响
2 a试验表明,适量减氮对根系生长与产量形成具有积极促进作用(
图3、
表3)。与CK处理相比,N1、N2处理在开花期显著提高RSAD值,增幅分别为13.40%、10.70%,N1与N2处理还显著增加现蕾期与成熟期根系RLD值;从产量构成来看,N1处理增加向日葵产量(5.05%),并显著提升向日葵百粒重(9.07%)、单盘籽粒重(9.08%);相比之下,N2与N3处理则显著降低成熟期RSAD值(23.66%、29.46%),N2、N3处理还降低向日葵产量(6.14%、20.22%),并显著降低百粒重(21.67%、29.51%)、单盘籽粒重(21.66%、29.49%)。综上得出,减氮25%处理(N1)促进关键生育期根系发育,最终有利于产量形成。根系分泌物分析表明,各处理下根系分泌的优势有机酸为草酸、甲酸和苹果酸,也检测到乙酸和柠檬酸(
表2)。N1处理显著促进现蕾期根系分泌草酸与柠檬酸(59.07%、60.11%);N2处理显著促进现蕾期和成熟期根系分泌甲酸(14.17%)。然而,当氮素水平低至N3处理(50 kg/hm
2)时,成熟期甲酸分泌受到抑制,且现蕾期、开花期与成熟期根系分泌苹果酸与草酸分泌量均显著降低。总体而言,适度减氮(尤其是N1处理)能够激发根系有机酸的分泌,可能与其在促进根系生长和产量形成中发挥的积极作用有关。
2.3 施氮条件下根际土壤理化性质与细菌反硝化途径及根系产量间的综合分析
由
表4、
表5可知,施氮处理下向日葵生育期根际土壤的理化性质。根际土壤中
-N与
-N含量随生育进程推进呈逐渐下降趋势。与CK处理相比,N2和N3处理的现蕾期根际土壤
-N含量均显著降低(22.21%、26.89%)。N1处理显著提高现蕾期根际土壤
-N含量(13.57%),而N3处理显著降低现蕾期与成熟期
-N含量。N1与CK处理的现蕾期、开花期与成熟期根际土壤SOM含量无显著差异,然而N2、N3处理显著降低现蕾期与成熟期根际土壤SOM含量,说明过量减氮可能影响土壤养分的维持能力。
由
图4可知,施氮处理的根际土壤可溶性离子含量。阴离子方面,N1、N2、N3处理显著降低现蕾期、开花期和成熟期的根际土壤
和
含量;阳离子方面,N1和N2处理均显著提高现蕾期根际土壤Ca
2+(17.22%、10.56%)和Mg
2+(34.75%、25.40%)含量。N2与N3处理显著降低现蕾期、开花期与成熟期K
++Na
+含量;N1处理的现蕾期根际土壤K⁺+Na⁺的含量较CK处理无显著差异。综上,适量减氮(N1、N2)处理有助于降低现蕾期根际土壤中
、
含量,同时促进Ca
2+、Mg
2+等有益阳离子的累积,可能有助于改善根际离子环境,为根系生长发育提供更为均衡的矿质营养条件。
由
图5可知,根际土壤理化因子与向日葵根系生长指标皮尔逊相关性分析结果。现蕾期反硝化功能基因
narG与SOM、
-N和
-N均呈极显著正相关,与
、Ca
2+、TN及甲酸含量呈显著正相关。
nosZ与柠檬酸、甲酸及
-N呈极显著正相关,与Mg
2+、Ca
2+和苹果酸呈显著正相关,而与pH和
呈显著负相关,表明现蕾期反硝化功能基因与根际环境因子关系密切,可能与该阶段根系分泌物活跃、植株氮素需求较高有关。开花期功能基因
narG和
nirK与甲酸呈显著正相关,与
呈显著负相关,开花期反硝化功能基因与土壤环境因子的相关性较现蕾期有所减弱。成熟期基因
narG、
nirS和
nosZ均与Ca
2+和乙酸呈显著正相关。在根系形态方面,现蕾期与开花期的根系表面积密度(RSAD)和根系体积密度(RVD)均与K
++Na
+和
呈极显著正相关,至成熟期,RSAD、RVD及根长密度(RLD)与
、K
++Na
+、
呈显著正相关,表明根际土壤中的盐分离子(尤其是K
++Na
+和
)可能是调控向日葵根系生长的关键环境因子。
采用主成分分析对2 a施氮处理下成熟期根际土壤理化性质、细菌反硝化功能基因、作物根系与产量指标进行综合评价(
表6、
表7),前2个主成分的特征值>1,累计贡献率达到86.53%,数据具有可信度。因此,选择2个主成分作为综合评价指标并计算各处理的综合得分。结果显示,不同处理组综合得分排序为N1>CK>N2>N3,说明减氮25%(N1)在改善根际土壤环境、优化反硝化功能基因组成、促进根系生长及提高产量方面的综合效果最佳,而过度减氮(N3)则表现最差。
3 讨 论
3.1 施氮处理对根际土壤养分与向日葵产量及根系生长和产量的影响
适度减氮有利于提高根长密度(RLD)与根表面积密度(RSAD)
[19],本研究减氮25%与50%处理显著增加现蕾期RLD与开花期RSAD值,可能与现蕾期追施尿素水解后土壤中铵态氮含量增加有关
[20],过量施氮促使
大量同化,进而抑制根系发育
[21]。本试验减氮25%处理较常规施氮处理减弱根际土壤
的同化,显著提高现蕾期根际土壤
-N含量(13.57%),从而为根系生长提供更充足的氮源。此外,胡柯鑫等
[22]研究也表明,合理减氮能够改善土壤理化性质,增加速效养分含量,促进根系生长发育
[23]。本研究减氮25%处理显著提高向日葵百粒重和单盘籽粒重,并最终促进产量的形成,说明在适宜减氮条件下,能够实现籽粒物质的有效积累,从而获得更高的产量,与石丽等
[24]研究结果一致。
减氮还可能触发植物的胁迫抗逆响应,促进根系分泌有机酸(如柠檬酸、草酸等)
[25],根系分泌物是土壤微生物的重要碳源。本研究减氮25%处理显著促进现蕾期根系分泌草酸与柠檬酸(59.07%和60.11%),减氮50%处理显著促进现蕾期和成熟期根系分泌甲酸(14.17%)。草酸和柠檬酸不仅能增强作物的耐盐性,还可通过螯合作用促进Mg
2+、Ca
2+等阳离子的释放,提高养分有效性
[26],也为本研究减氮25%与50%处理显著提高现蕾期根际土壤Ca
2⁺(17.22%、10.56%)与Mg
2⁺(34.75%、25.40%)含量提供合理解释。值得注意的是,减氮50%处理在增加现蕾期根际土壤
与Mg
2⁺含量的同时,显著降低现蕾期K
++Na
+含量。可能是因为Ca
2+和Mg
2+可置换土壤胶体吸附的Na
+,促进胶体絮凝,改善土壤结构,增强土壤对有益阳离子的吸附能力,从而优化根际离子环境,缓解盐渍化风险。
3.2 施氮处理对根际土壤反硝化细菌群落及其介导的反硝化途径的影响
施氮显著影响土壤中反硝化微生物丰度与群落结构
[27-28],放线菌门与变形菌门在反硝化作用中发挥重要作用
[29],放线菌门作为反硝化途径中的优势菌门,自身就具有抵抗病害,保护植物健康的能力
[30],本研究现蕾期减氮25%处理显著提高放线菌门相对丰度促进向日葵根系生长并维持健康的根际土壤环境。变形菌门(Proteobacteria)中的反硝化细菌依赖高氮环境,可通过降低硝态氮等底物供应来抑制其活性,从而达到减少氮素淋失风险的目的,本研究减氮处理显著降低
narG、
nirK、
nirS型Proteobacteria相对丰度,可能直接使根际土壤反硝化潜势下降,有助于将更多的氮素以硝酸盐(
)的形式留存于土壤提高氮肥利用率。
Ouyang等
[31]通过Meta分析发现,农业上施用氮肥可显著提高反硝化基因(
nirK、
nirS)的丰度,本研究通过减氮处理降低
nirS、
nirK、
narG基因相对丰度,此结论与李志龙等
[32]与宋亚娜等
[33]研究相同,减量施氮减少反硝化细菌的底物抑制其生长繁殖。Wang等
[34]、Gao等
[35]研究表明,施氮显著影响土壤N
2O排放及
nosZ型反硝化细菌的数量和群落结构,黄莹等
[36]通过对强酸性茶园土壤的培养发现,施氮显著降低土壤细菌
nosZ基因的丰度,然而本研究减氮25%处理(150 kg/hm
2)却特异性地提升根际土壤细菌
nosZ基因的丰度,可能与减氮促进根系分泌有机酸密切相关;本研究分析得出,现蕾期根际细菌
nosZ基因与柠檬酸、甲酸、苹果酸等呈显著正相关,成熟期根际细菌
nosZ基因与乙酸呈显著正相关,进一步说明,减氮处理下细菌功能基因
nosZ的富集与根系有机酸的分泌密切相关。有机碳含量是影响
nosZ基因丰度的重要因子,且介导N
2O还原的反硝化细菌(
nosZ)对碳源尤为敏感,故本研究适度减施氮肥可能促使根系分泌有机酸,特异性地为
nosZ型微生物提供碳源,创造其竞争优势。
4 结 论
1)减氮处理(N1、N2、N3)显著抑制根际细菌反硝化功能基因narG、nirK、nirS的表达;减氮25%处理(N1)在抑制成熟期根际细菌narG、nirK、nirS表达的同时,促进nosZ的表达,不仅从源头上阻控盐渍农田根际土壤氮素损失,还能减小盐渍向日葵农田土壤中N2O的排放潜力。
2)N1处理显著促进现蕾期根系分泌草酸(59.07%)与柠檬酸(60.11%),并增加根际土壤Ca2⁺(17.22%)与Mg2⁺(34.75%)并降低、含量,有效改善盐渍农田根际土壤环境,并显著增加开花期RSAD值(13.40%)、百粒重(9.07%)与单盘籽粒重(9.08%),促进产量形成。
3)减氮50%(N2)处理显著提高现蕾期根际土壤Ca2⁺(10.56%)与Mg2⁺(25.40%),却显著降低成熟期根系RSAD值、百粒重与单盘籽粒重;N2与N3分别导致向日葵减产6.14%和20.22%。基于主成分分析综合评价成熟期根际土壤养分、反硝化功能基因、根系生长及产量指标,减氮25%(150 kg/hm2)处理表现最佳。
国家自然科学基金项目(U2443210)
国家自然科学基金项目(52469012)
国家自然科学基金项目(52069022)
内蒙古自治区自然科学基金项目(2024LHMS05022)
内蒙古自治区科技重大专项(2023JBGS0003)