撂荒地亚硝酸还原酶基因 nirK nirS丰度动态
呼和1, 陈先江2, 程云湘2
1.甘肃省农业科学院土壤肥料与节水农业研究所,甘肃 兰州 730070
2.草地农业生态系统国家重点实验室 兰州大学草地农业科技学院,甘肃 兰州 730020
通讯作者:程云湘(1977-),女,山东菏泽人,副教授,博士,主要从事草地生态、草地温室气体排放研究。E-mail:chengyx@lzu.edu.cn

第一作者:呼和(1976-),男(蒙古族),内蒙古通辽人,副研究员,博士,主要从事环境污染的微生物修复、农牧业固体废弃物循环利用和微生物生态研究。E-mail:huhecyxw@yahoo.co.jp

摘要

使用荧光定量PCR法测定亚硝酸还原酶基因( nirK nirS)遗传特异性片段,研究了草原撂荒地土壤反硝化微生物丰度随着撂荒时间的变化。结果表明,轻度放牧草原和3种不同撂荒时间地的土壤中 nirK基因型反硝化微生物丰度都显著高于 nirS基因型反硝化微生物丰度( P<0.05)。土壤 nirK nirS基因型反硝化微生物丰度撂荒地和轻度放牧地之间均有显著差异( P<0.05),但3种撂荒地均未发生显著变化( P>0.05)。此外,这两种基因型反硝化微生物丰度之间有极显著线性负相关关系( P<0.001)。以上结果说明, nirK nirS基因型反硝化微生物对环境因子的响应有差异。

关键词: 反硝化微生物; nirK基因; nirS基因; 撂荒地; 轻度放牧草原
中图分类号:S812.29 文献标志码:A 文章编号:1001-0629(2016)7-1253-07 doi: 10.11829/j.issn.1001-0629.2015-0751
The dynamics for abundance of nitrite reductase genes nirK and nirS in abandoned land
Huhe1, Chen Xian-Jiang2, Cheng Yun-Xiang2
1.Institute of Soil and Fertilizer and Save Water Agricultural, Gansu Academy of Agricultural Sciences, Lanzhou 730070, China
2.State Key Laboratory of Grassland Agro-ecosystems, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou 730020, China
Corresponding author: Cheng Yun-xiang E-mail:chengyx@lzu.edu.cn
Abstract

Determination of nitrite reductase gene ( nirK and nirS) specific fragment by real-time quantitative PCR method, we studied the dynamics of denitrification microbial abundance in abandoned land soil with different abandoned years. The results showed that the nirK genotype denitrifier abundance were significantly higher than nirS genotype denitrifier abundance in 3 kinds of abandoned land for 1, 5, and 25 years and light grazing grassland. Both of genotypes denitrifier abundance were significant difference between the light grazing grassland and the three abandoned land soils ( P<0.05), but there was no significant change among the 3 kinds of abandoned land soil. Moreover, there were extremely significantly negative linear relationship between the two genes abundance in the investigation soils ( P<0.001), indicating that the nirK and nirS show different response to the environmental factors on their growth.

Keyword: denitrifier; nirK gene; nirS gene; abandoned land; light-grazing grassland

由于人口迅速增长, 我国对粮食及经济作物的需求量日益增大, 因此, 需要多种途径扩大耕地面积, 如开垦荒地、围海造田、围湖造田等, 这就造成了草地退化等现象, 如开垦已成为呼伦贝尔草原沙化、退化的主要原因之一[1]。近年来, 国家开始实施退耕还林、还草等工程, 以改善生态环境[2]。目前, 国内外关于退耕地的研究主要集中于植被群落的演替[3]、土壤理化性质的变化[4]以及种植植物的影响[5, 6]等方面, 但随着退耕时间的延长, 土壤微生物群落如何变化, 尤其是关于微生物氮循环重要环节之一— — 反硝化基因丰度变化及其相互关系的研究鲜见报道。生物反硝化作用是氮循环的主要环节之一, 反硝化作用在土壤中普遍存在, 其主要是通过微生物逐步地还原氮氧化合物, 释放出一氧化氮(NO)、氧化亚氮(N2O)和氮(N2)等, 从而导致土壤中氮元素的减少, 而且, N2O和NO的释放可能导致全球变暖[7]。反硝化微生物广泛分布于细菌和古菌中, 在几种真菌的线粒体中也发现了反硝化作用[8]。反硝化作用主要指在硝酸还原酶(Nar)、亚硝酸还原酶(Nir)、氧化亚氮还原酶(Nor)和一氧化二氮还原酶(Nos)等的作用下将硝酸盐(N O3-)中的氮通过一系列中间产物(N O2-、NO、N2O)还原为氮气(N2)的生物化学过程。其中, 亚硝酸盐向NO的还原酶和其它反硝化酶不同。首先, 与其它反硝化反应相比, 亚硝酸还原酶把亚硝酸盐还原为NO的反应是将亚硝酸盐还原成气体的最初阶段。其次, 这个反应由两种不同类型的亚硝酸盐还原酶(Nir)催化, 这两种酶存在于不同的细胞内且结构不同[8]。一种是nirS编码的含细胞色素cd1的酶, 一种是由nirK编码的含铜的酶。据报道[9, 10], 两种亚硝酸盐还原酶功能上没有显著差异, 但是, 如果土壤理化特性不同(如, 氧、碳含量等), 有可能观察到不同的反硝化微生物[11, 12]或不同反硝化微生物群落[13, 14]。然而, 目前针对nirSnirK基因丰度在不同特性土壤中的变化及其之间相互关系的研究较少。因此, 本研究使用实时荧光定量PCR法初步探讨在撂荒时间不同的耕地土壤中的nirSnirK基因丰度的变化规律, 以利于更好地理解撂荒地土壤反硝化作用, 从而为管理撂荒地提供理论依据。

1 材料和方法
1.1 研究区概况

研究区域位于内蒙古自治区东北部的呼伦贝尔草原新巴尔虎右旗境内(116° 40' 18''-117° 01' 56'' E, 48° 32' 00''-48° 38' 50'' N), 为半干旱草原。试验共设4个处理, 3个撂荒农田和1个轻度放牧草地。3个撂荒农田的撂荒时间分别为1、5和25年(Y1、Y5和Y25), 轻度放牧草地(LG)设为对照, 放牧率为1.4只羊· hm-2。在Y1、Y5和Y25样点, 玉米(Zea mays)、向日葵(Helianthus annuus)、圆柱披碱草(Elymus cylindricus)无规律轮作约40年, 各样地之间距离< 10 km, 且各样点地势相似, 均为平地。各样地2000-2009年平均温度和降水量分别为1.6 ℃和213 mm, 海拔在545~568 m, 均为沙质土壤。

1.2 调查取样及分析

在每个样点随机设5个1 m× 1 m的样方, 于2010年8月进行植被调查和土壤取样。记录每个样方内的植物物种组成, 测定样方的盖度、高度和密度, 并计算Shannon-Wiener多样性指数。剪取植物地上部分, 按物种分类, 60 ℃下干燥48 h, 称重, 测定地上植物干重。土壤湿度用TRIME-FM (Ettlingen, Germany)测定。另外, 在每个样方随机选择5个点收集0-10 cm土样, 混合为一个样品, 并移去表面有机质和细根。把每个混合样品分为两部分, 一部分过2 mm筛后密封于15 mL离心管中, 放在液氮罐里运回实验室, 储存于-80 ℃冰箱中待用于微生物试验, 另一部储存于4 ℃下, 用于测定土壤理化特性。土壤KCl浸提物中的NO3-N和NH4-N分别用锌还原-萘乙二胺法[15]和靛酚蓝比色法[16]测定, 土壤磷含量用Truog法测定[17], 有机碳和总氮使用NC分析仪(Sumigraph NC-900; Sumika Chemical Analysis Service, Tokyo, Japan)采用干燃烧法测定[18]。土壤pH采用电位法测定。

1.3 土壤DNA提取

每个土样分别称取3份, 每份0.5 g, 采用FastDNA Spin Kit(MP Biomedicals, Illkirch, France)提取土壤DNA, 实验步骤按照说明书进行。采用SmartSpec Plus spectrophotometer (Bio-Rad Laboratories, United States)测定DNA的纯度及浓度。样品保存于-20 ℃冰箱中备用。

1.4 荧光实时定量PCR(qPCR)

采用实时定量PCR(Applied Biosystems StepOne Plus 96-well, 美国)与预混液(Applied Biosystems Power SYBR Green, 美国)相对定量nirKnirS基因片段, 各反应的引物、qPCR反应条件、扩增片段大小、标准微生物及混合体系组成如表1表2所示。qPCR反应开始前先在95 ℃下进行10 min的酶活化。扩增产物的特异性通过观察熔解曲线上的单一溶解峰和溴化乙锭染色的2%琼脂糖凝胶上特异条带的出现来检测。连续稀释(10~107拷贝数· μ L-1)标准细菌对应基因PCR扩增产物制作各基因的标准曲线, 使用Eff=[10(-1/slope)-1]公式计算nirKnirS基因扩增效率(Eff), 分别为92.0%和98.2%。每个反应3个重复, 平均拷贝数转换为每克干土样中的基因拷贝数。

表1 nirKnirS基因的qPCR条件及引物 Table 1 qPCR Condition of nirK and nirS gene and primer
表2 nirKnirS基因的qPCR反应体系 Table 2 qPCR reaction system of nirK and nirS gene
1.5 数据分析

qPCR数据和环境因子特征采用SPSS 19.0(SPSS, Inc., Chicago, IL, 美国)进行方差分析(ANOVA)。ANOVA分析采用了Tukey’ s test检验。如果每组数据不是正态分布, 在分析前将数据进行对数变换。采用SIMCA-P 11.5(Umetrics AB, Umea, 瑞典)分析软件, 偏最小二乘法(PLS)分析nirSnirK基因拷贝数与环境因子之间的回归系数和变量投影重要性(VIP)。

2 结果与分析
2.1 nirKnirS基因丰度

开垦耕种活动显著影响了nirKnirS基因丰度比例(图1)。撂荒地nirK基因的丰度是nirS基因的1 000倍左右, 然而, 在LG土壤中, nirK基因丰度是nirS基因10倍左右。随着撂荒时间的延长, 这两种基因丰度都没有显著变化(P> 0.05), 但是, 3个撂荒地的nirK基因的丰度均显著高于LG的(P< 0.05)。另外, 4个样地亚硝酸还原酶基因nirK的拷贝数和nirS基因的拷贝数存在极显著线性关系(图2), 回归方程为nirS=-0.696 (nirK)+9.650 7(R2=0.597 7, P< 0.001)。

图1 撂荒地和轻度放牧草原土壤nirKnirS基因丰度
注:Y1撂荒1年; Y5撂荒5年; Y25撂荒年25年; LG轻度放牧地。不同小写字母表示nirSnirK基因丰度在不同样地间差异显著(P< 0.05)。
Fig.1 Abundance of nirK and nirS gene of soil of abandoned land and light-grazing steppe grassland
Note: Y1, abandoned land with 1 year; Y5, abandoned land with 5 years; Y25, abandoned land with 25 years; LG, litght grazing grassland. Different lower case letters indicate significant difference of the same gene abundance among different plots at 0.05 level.

图2 nirKnirS基因丰度之间关系Fig.2 Relation between abundances of nirK and nirS gene

2.2 环境因子与亚硝酸还原酶基因的相互作用

两种亚硝酸还原酶基因(nirKnirS)都与NH4-N、有机碳、总氮、植物多样性和地上生物量显著相关, VIP值均大于1(图3)。由回归系数值可知, nirK基因拷贝数与NH4-N、有机碳、总氮和地上生物量呈负相关, 与植物多样性呈正相关, 而nirS基因的拷贝数与nirK基因的相关性分析结果正好相反。此外, pH、NO3-N、有效磷和H2O等环境因子对nirKnirS基因的拷贝数无显著影响, VIP值都小于1。

图3 各环境因子和nirK, nirS基因丰度之间的回归系数及VIP值Fig.3 Regression coefficient between each environmental factors and abundance of nirK and nirS, and VIP value

3 讨论与结论

土壤氮循环主要依赖于微生物的氮循环作用, 即固氮、亚硝化、硝化及反硝化等作用把不同类型的氮素化合物改变为亚硝酸根、硝酸根、NO及N2O的形式。其中, 具有亚硝酸还原酶基因的微生物能把亚硝酸盐还原为NO, 是生物氮循环这一过程的重要影响环节之一, 与其它参与生物氮循环基因相同, 亚硝酸还原酶基因也是研究者关注的重要反硝化功能基因。因此, 本研究以nirKnirS基因作为研究对象, 采用实时荧光PCR法探讨了不同类型亚硝酸还原酶基因丰度在呼伦贝尔草原不同撂荒时间耕地土壤中的动态变化以及其与环境因子之间的关系。

研究发现, 撂荒地和轻度放牧草原的土壤理化特征及植被状况差异较大[19]。3个撂荒地的地上生物量有显著差异(P< 0.05), 但铵态氮、有效磷、有机碳和植物多样性无显著差异(P> 0.05)[19]。此外, Y5样地的土壤pH和硝酸盐氮(NO3-N)显著高于Y1和Y25, Y1和Y25之间无显著差异。Y5的总氮显著高于Y1的(P< 0.05), 但Y25的与Y1和Y5的均无显著差异。Y25的土壤水分含量(H2O)显著高于Y1的, Y5的与Y1、Y25的均无显著差异[19]

在3个不同撂荒时间的撂荒地土壤中, nirK基因丰度显著高于nirS基因, 前者均约为后者的103倍, 而在轻度放牧草原土壤中, 前者约为后者的10倍, 这说明该研究地土壤中具有nirK基因的微生物是主要的亚硝酸还原微生物。在森林土壤(壤土和粉壤土)[20]和稻田土壤[21]中也得到了相同的结果, 即nirK基因丰度高于nirS基因丰度。此外, 在撂荒地土壤中nirKnirS基因丰度分别显著高于和低于轻度放牧草原。根据PLS分析的结果, 植物多样性和nirK基因丰度之间有显著正相关(VIP=1.05), 表明nirK类型反硝化微生物有可能很容易受不同组成和数量的根系分泌物的影响, 在沙壤土中, 随着植被多样性的变化其土壤微生物群落结构也变化[22], Bremer等[23]使用狐尾草(Alopecurus pratensis)、黄花茅(Anthoxanthum odoratum)、燕麦草(Arrhenatherum elatius)、绒毛草(Holcus lanatus)、老鹳草(Geranium pratense)、车前(Plantago lanceolata)、毛茛(Ranunculus acris)和蒲公英(Taraxacum officinale)的微宇宙试验中发现, 不同植物微宇宙土壤有不同的nirK群落结构。Baneras等[24]也指出植被类型的改变对nirK类型反硝化微生物群落具有显著影响。除此之外, NH4-N、有机碳、总氮和地上生物量对于nirKnirS基因丰度分别呈显著负相关和正相关关系(VIP > 1)(图3), 说明亚硝酸还原细菌的丰度变化受到土壤不同理化性质的共同影响, 导致两种类型亚硝酸还原细菌丰度呈此消彼长的趋势, nirKnirS基因丰度之间呈极显著线性关系(P< 0.001)(图2)。许多研究[25, 26, 27]表明, 土壤微生物的分布、活动、多样性以及丰度与土壤理化性质密切相关。Maag和Vinther[28]在森林土壤中发现, 反硝化速率受多个环境因子联合调控。 尽管nirKnirS基因与NO3-N的含量密切相关, 而本研究未发现nirKnirS基因丰度与NO3-N含量具有显著相关性, 这有可能是在自然环境土壤中不同环境因子共同影响亚硝酸还原微生物的结果。许多研究者也在不同环境的土壤中发现了反硝化细菌群落与NO3-N含量具有很弱的相关性[29, 30, 31]。此外, Kong等[32]认为, 在亚硝酸还原微生物生存的微生境中, 其变化有可能与NO3-N含量变化无关。

总之, 在本研究区的撂荒地, nirKnirS基因型反硝化微生物丰度随着撂荒时间无有显著变化, 然而它们之间有显著线性关系, 两种基因型反硝化微生物随着环境因子的变化, 有可能互补调节亚硝酸盐还原为NO的作用。但本研究只分析了DNA水平上的反硝化微生物丰度, 其活性及表达的规律还有待RNA样品来分析说明。另外, 尚需开展nirKnirS基因型反硝化微生物群落结构和分布, 以便了解随着撂荒时间的延长反硝化微生物群落的演替规律, 更科学地通过调节反硝化作用管理撂荒地。

The authors have declared that no competing interests exist.

参考文献
[1] 张德平. 呼伦贝尔会消失吗——对草原区耕地保护的反思. 中国土地, 2009(4): 44-46. [本文引用:1]
[2] Li J H, Fang X W, Jia J J, Wang G. Effect of legume species introduction to early aband oned field on vegetation development. Plant Ecology, 2007, 191: 1-9. [本文引用:1]
[3] Cheng Y X, Nakamura T. Phytosociological study of steppe vegetation in east Kazakhstan. Grassland Science, 2007, 53: 172-180. [本文引用:1]
[4] 李强, 刘国彬, 许明祥, 张正, 孙会. 黄土丘陵区撂荒地土壤抗冲性及相关理化性质. 农业工程学报, 2013, 29(10): 153-159.
Li Q, Liu G B, Xu M X, Zhang Z, Sun H. Soil anti-scouribility and its related physical properties on aband oned land in the Hilly Loess Plateau. Transactions of the Chinese Society of Agricultural Engineering, 2013, 29(10): 153-159. (in Chinese) [本文引用:1]
[5] Van der Putten W H, Mortimer S R, Hedlund K, Van Dijk C, Brown V K, Lepš J, Rodriguez-Barrueco C, Roy J, Diaz Len T A, Gormsen D, Korthals G W, Lavorel S, Santa Regina I, Smilauer P. Plant species diversity as a driver of early succession in aband oned fields: A multi-site approach. Oecologia, 2000, 124: 91-99. [本文引用:1]
[6] 薛智德, 侯庆春, 韩蕊莲, 王胜琪. 黄土丘陵沟壑区白刺花促进生态恢复的研究. 西北林学院学报, 2002, 17(3): 26-29.
Xue Z D, Hou Q C, Han R L, Wang S Q. Trails and research on ecological restoration by Sophora viciifolia in Gullied Rolling Loess Region. Journal of Northwest Forestry University, 2002, 17(3): 26-29. (in Chinese) [本文引用:1]
[7] Houghton J T, Calland er B A, Varney S K. Climate Change 1992: The Supplementary Report to the IPCC Scientific Assessment. New York: Cambridge University Press, 1992. [本文引用:1]
[8] Zumft W G. Cell biology and molecular basis of denitrification. Microbiology and Molecular Biology Reviews, 1997, 61: 533-616. [本文引用:2]
[9] Priemé A, Braker G, Tiedje J M. Diversity of nitrite reductase (nirK and nirS) gene fragments in forested upland and wetland soils. Applied and Environmental Microbiology, 2002, 68: 1893-1900. [本文引用:1]
[10] Ruiz-Rueda O, Trias R, Garcia-Gil T J, Bañeras L. Diversity of the nitrite reductase gene nirS in the sediment of a free-water surface constructed wetland . International Microbiology, 2007, 10: 253-260. [本文引用:1]
[11] Abou-Seada M N I, Ottow J C G. Einfluß chemischer bodeneigenschaften auf ausmaß and zusammensetzung der denitrifikationsverluste drei verschiedener bakterien. Zeitschrift für Pflanzenernährung und Bodenkunde, 1988, 151: 109-115. [本文引用:1]
[12] Tiedje J M. Ecology of denitrification and dissimilatory nitrate reduction to ammonium. In: Zehnder A J B. (ed). Biology of Anaerobic Microorganisms. New York: John Wiley & Sons, 1988: 179-244. [本文引用:1]
[13] Cavigelli M, Robertson G P. The functional significance of denitrifier community composition in a terrestrial ecosystem. Ecology, 2000, 81: 1402-1414. [本文引用:1]
[14] Holtan-Hartwig L, Dörsch P, Bakken L R. . Comparison of denitrifying communities in organic soils: Kinetics of NO3- and N2O reduction. Soil Biology and Biochemistry, 2000, 32: 833-843. [本文引用:1]
[15] Leonardo M. Development and validation of a method for determination of residual nitrite/nitrate in foodstuffs and water after zinc reduction. Food Analytical Methods, 2009, 2: 212-220. [本文引用:1]
[16] Motsara M R, Roy R N. Guide to laboratory establishment for plant nutrient analysis. FAO Fertilizer and Plant Nutrition Bulletin No. 19. Rome: FAO, 2008: 17-76. [本文引用:1]
[17] Truog E, Meyer A H. Improvements in the deniges colorimetric method for phosphorus and arsenic. Inddustrial & Engineering Chemistry Analytical Edition, 1929, 1(3): 136-139. [本文引用:1]
[18] Ueno T, Osono T, Kand a H. Inter- and intraspecific variations of the chemical properties of high-Arctic mosses along water-regime gradients. Polar Science, 2009, 3(2): 134-138. [本文引用:1]
[19] Huhe, Borjigin S, Cheng Y X, Nomura N, Nakajima T, Nakamura T, Uchiyama H. Effect of aband onment on diversity and abundance of free-living nitrogen-fixing bacteria and total bacteria in the cropland soils of Hulun Buir, Inner Mongolia. PLoS One, 2014, 9(9): e106714. [本文引用:3]
[20] Levy-Booth D J, Winder R S. Quantification of nitrogen reductase and nitrite reductase genes in soil of thinned and clear-cut Douglas-fir stand s by using real-time PCR. Applied and Environmental Microbiology, 2010, 76: 7116-7125. [本文引用:1]
[21] Bannert A, Kleineidam K, Wissing L, Mueller-Niggemann C, Vogelsang V, Welzl G, Cao Z, Schloter M. Changes in diversity and functional gene abundances of microbial communities involved in nitrogen fixation, nitrification, and denitrification in a tidal wetland versus paddy soils cultivated for different time periods. Applied and Environmental Microbiology, 2011, 77: 6109-6116. [本文引用:1]
[22] Kowalchuk G A, Buma D S, De Boer W, Klinkhammer P G L, van Veen H. Effects of above-ground plant species composition and diversity on the diversity of soil-borne microorganisms. Antonie Van Leeuwenhoek, 2002, 81(1-4): 509-520. [本文引用:1]
[23] Bremer C, Braker G, Matthies D, Reuter A, Engels C, Conrad R. Impact of plant functional group, plant species, and sampling time on the composition of nirK-type denitrifier communities in soil. Applied and Environmental Microbiology, 2007, 73: 6876-6884. [本文引用:1]
[24] Bañeras L, Ruiz-Rueda O, López-Flores R, Quintana X D, Hallin S. The role of plant type and salinity in the selection for the denitrifying community structure in the rhizosphere of wetland vegetation. International Microbiology, 2012, 15: 89-99. [本文引用:1]
[25] Marschner P, Yang C H, Lieberei R, Crowley D E. Soil and plant specific effects on bacterial community composition in the rhizosphere. Soil Biology and Biochemistry, 2001, 33: 1437-1445. [本文引用:1]
[26] Pengthamkeerati P, Motavalli P, Kremer R. Soil microbial activity and functional diversity changed by compaction, poultry litter and cropping in a claypan soil. Applied Soil Ecology, 2011, 48: 71-80. [本文引用:1]
[27] Cruz-Martínez K, Rosling A, Zhang Y, Song M, Andersen G L, Banfield J F. Effect of rainfall-induced soil geochemistry dynamics on grassland soil microbial communities. Applied and Environmental Microbiology, 2012, 78: 7587-7595. [本文引用:1]
[28] Maag M, Vinther F P. Nitrous oxide emission by nitrification and denitrification in different soil types and at different soil moisture contents and temperatures. Applied Soil Ecology, 1996, 4(1): 5-14. [本文引用:1]
[29] Mergel A, Schmitz O, Mallman T, Bothe H. Relative abundance of denitrifying and dinitrogen-fixing bacteria in layers of a forest soil. FEMS Microbiology Ecology, 2001, 36(1): 33-42. [本文引用:1]
[30] Chen Y N, Zhou W, Li Y P, Zhang J C, Zeng G M, Huang A Z, Huang J X. Nitrite reductase genes as functional markers to investigate diversity of denitrifying bacteria during agricultural waste composting. Applied Microbiology and Biotechnology, 2014, 98: 4233-4243. [本文引用:1]
[31] Smith C J, Nedwell D B, Dong L F, Osborn A M. Diversity and abundance of nitrate reductase genes (narG and napA), nitrite reductase genes (nirS and nirK), and their transcripts in estuarine sediments. Applied and Environmental Microbiology, 2007, 73: 3612-3622. [本文引用:1]
[32] Kong A, Hristova K, Scow K, Six J. Impacts of different N management regimes on nitrifier and denitrifier communities and N cycling in soil microenvironments. Soil Biology and Biochemistry, 2010, 42: 1523-1533. [本文引用:1]