减氮配施有机肥对燕麦土壤细菌群落特征影响与生态功能预测.pdf

山西农业科学 2023，51（10）：1196-1202Journal of Shanxi Agricultural Sciences减氮配施有机肥对燕麦土壤细菌群落特征影响与生态功能预测郑敏娜，邹璨阳，康佳惠，杨富（山西农业大学 高寒区作物研究所，山西 大同 037008）摘要：为了探讨减氮配施有机肥处理下燕麦根际土壤细菌群落多样性的差异状况，以晋燕 17号为试验材料，采用高通量测序技术，分析不施肥（CK）、仅施用尿素（T1，180 kg/hm2）和尿素与有机肥分别按照 9 1（T2）、8 2（T3）、7 3（T4）、6 4（T5）和 5 5（T6）的比例配施等 7个处理对土壤养分状况和细菌群落多样性的影响，并利用PICRUSt法预测其生态功能。结果表明，减氮配施有机肥处理下，随着有机肥添加比例的增加，土壤中有机质含量逐渐增加，且在 T6 处理时达到最大（14.86 g/kg），籽粒产量在 T4 处理下达到最高（4 357.61 kg/hm2）；与不施有机肥土壤相比，增施有机肥后显著提升了变形菌门、厚壁菌门的相对丰度，丰度前 5的优势细菌门依次为变形菌门、厚壁菌门、拟杆菌门、放线菌门、绿弯菌门；影响燕麦田土壤细菌群落结构的主要环境因子是速效氮（AN）、有机质（OM）、速效钾（AK）、全氮（TN）。PICRUSt预测结果表明，减氮配施有机肥后明显降低了土壤细菌细胞转化功能基因的丰度，明显增加了生物体系统功能基因的丰度。综上，T4处理（126 kg/hm2尿素+54 kg/hm2有机肥）下生物产量最高，且其优势菌群更有利于土壤的可持续生产。关键词：减施氮肥；有机肥；高通量测序；PICRUSt法；土壤细菌群落特征；生态功能预测中图分类号：S512.6 文献标识码：A 文章编号：10022481（2023）10119607Effects of Nitrogen Reduction Combined with Organic Fertilizer and Prediction of Ecological Functions of Bacterial Community Characteristics in Oat SoilZHENG Minna，ZOU Canyang，KANG Jianhui，YANG Fu（Institute of High Latitude Crops，Shanxi Agricultural University，Datong 037008，China）Abstract：In order to investigate the difference of bacterial community diversity in oat rhizosphere soil treated with nitrogen reduction combined with organic fertilizer，in this study,Jinyan 17 was taken as the experimental material,high-throughput sequencing technology was adopted,the effects of the seven treatments including no fertilization(CK),carbamide only(T1,180 kg/ha）and carbamide combined with organic fertilizer in the ratios of 9 1(T2),8 2(T3),7 3(T4),6 4(T5)and 5 5(T6)on soil nutrient status and bacterial community diversity were analyzed,and the ecological functions were predicted by PICRUSt method.The results showed that the content of organic matter in soil gradually increased with the increase of the proportion of organic fertilizer under the treatments of nitrogen reduction combined with organic fertilizer,reaching the maximum under T6 treatment(14.86 g/kg),and the grain yield reached the maximum under T4 treatment(4 357.61 kg/ha).Compared with the soil without organic fertilizer,the relative abundance of Proteobacteria and firmicutes was significantly increased after increased application of organic fertilizer.The top 5 dominant bacterial phyla in abundance were Proteobacteria,Firmicutes,Bacteroidetes,Actinobacteria,and Green Curved Bacteria in order.The main environmental factors affecting the soil bacterial community structure in oat fields were available nitrogen(AN),organic matter(OM),available potassium(AK),and total nitrogen(TN).The PICRUSt prediction results indicated that the abundance of transforming functional genes of bacterial cells was significantly decreased after nitrogen reduction combined with organic fertilizer,and the abundance of functional genes of biological systems was significantly increased.In conclusion,under T4 treatment(126 kg/ha carbamide+54 kg/ha organic fertilizer),the biological yield was the highest,and the dominant bacterial community was more conducive to sustainable production of soil.Key words：nitrogen reduction;organic fertilizer;high-throughput sequencing technology;PICRUSt method;soil bacterial community characteristics;ecological function prediction燕麦（Avena sativa L.）为 1 年生禾本科粮饲兼用型作物，在全球 42个国家和地区广泛栽培，居世doidoi:10.3969/j.issn.1002-2481.2023.10.11收稿日期：2023-07-31基金项目：山西农业大学育种工程项目（YZGC085）；财政部和农业农村部：国家现代农业产业技术体系（CARS-07-A-2）作者简介：郑敏娜（1983-），女，甘肃通渭人，副研究员，硕士，主要从事燕麦栽培等研究工作。1196郑敏娜等：减氮配施有机肥对燕麦土壤细菌群落特征影响与生态功能预测界谷物生产第 6 位，主要分布于欧洲、亚洲和非洲的温带地区1。我国是燕麦的起源地之一，主要分布在华北、西北和西南等地区。燕麦具有耐贫瘠、耐盐碱、耐寒、抗旱、抗风沙等特点2-3，适应性广泛，且栽培措施简易，在营养价值、生产性能、利用方式、开发潜力等方面优势突出，是山西省传统优势杂粮作物之一，尤其是大同地区的特色优势杂粮作物，很受市场的青睐，具有很大的发展空间和利用潜力，在山西区域经济建设中发挥着重要作用。施肥是燕麦田间管理中非常重要的措施之一，对燕麦田生产力的提升和土壤生态具有重要作用。肥料（包括无机肥和有机肥）输入土壤后，会对土壤生态环境产生巨大影响，合理施用不但可以直接提升作物产量，而且可以改变土壤微生态，有利于土壤微生物的正向演替；但在实际生产中，为了追求更高的产量，化肥施用量偏高，有机肥投入不足，造成土壤质量下降等问题，不但造成了土壤资源严重污染，而且会威胁食品安全4。因此，在燕麦生产过程中，探索合理的减氮配施有机肥比例，将有利于促进土壤和生产的可持续发展，对晋北地区燕麦田土壤生态系统的可持续发展具有重要意义。在土壤生态系统中，微生物是重要的参与者，其对土壤养分循环、能量流动和物质转化等起到至关重要的作用5-9，是生态系统健康可持续发展的重要保障。在土壤生态系统中，相同的微生物种群可以行使不同的土壤功能，相同的土壤过程也可由不同微生物种群完成，其种群会因环境变更或管理分异产生微生物多样性及其生态功能的变化9。目前，众多学者在施肥与土壤微生物及其生态系统方面进行了大量研究，大多试验结果都表明，配施一定比例的有机肥料后土壤微生物量碳、氮、磷和生物活性均有不同程度的提高，具体效果会因配施比例或施肥种类而有所差异10-14。然而，也有研究显示，长期氮磷肥配施对土壤微生物生物量的影响不大15。本研究通过田间小区试验，研究在不同剂量无机有肥料配施下，燕麦田土壤环境及其微生物的生态差异，揭示不同剂量肥料的输入对燕麦田土壤生态功能的影响，优选出最佳施肥模式，为晋北地区燕麦的种植和提质增效提供参考意见。1 材料和方法1.1试验地概况试验于 2022 年在山西农业大学高寒区作物研究所毛皂试验基地（N11234 11434，E3903 4044）进行。该地海拔 1 050 m，属中温带干旱区，土壤有机质含量为 12.61 g/kg，土壤全氮含量为0.81 g/kg，土壤有效磷含量为 5.03 mg/kg，土壤 pH值为 8.37。1.2试验材料选用燕麦品种晋燕 17 号为试验材料，由山西农业大学高寒区作物研究所提供。试验所用肥料为尿素和自制有机肥，其中，尿素 N 含量46%；自制 有 机 肥 以 牛 粪 为 主 要 原 料，其 有 机 质 含 量 为15.8%，N含量为 0.046 g/kg。1.3试验设计试验采用随机区组排列。播前进行精细耙耱镇 压，平 整 后 人 工 开 沟 条 播，小 区 面 积 为 10 m2（2 m5 m），3 次重复，种植密度为 400 万株/hm2，播种深度 35 cm，行间距约 25 cm。试 验 共 设 7 个 处 理，以 不 施 肥 处 理 为 对 照（CK），分别设置仅施用氮肥（T1）和氮肥与有机肥分别按照 9 1（T2）、8 2（T3）、7 3（T4）、6 4（T5）和5 5（T6）的比例配施，共计 21个小区，其中，T1处理N 的施用量为 180 kg/hm2，T2T6等 5个处理 N 的总用量与T1处理一致，对应有机肥施用量按照配施比例计算，P、K施用量各小区均相同，试验期间各处理田间管理保持一致，具体施用量如表1所示。1.4测定项目及方法1.4.1土壤样品采集与测定于燕麦成熟期，各小区按对角线 5点采集燕麦根际土壤，采用四分法取0.2 kg土壤样品，3次重复。取样后，一部分土样，风干后测定常规土壤化学指标16，包括全氮（TN）、速效氮（AN）、速效磷（AP）、速效钾（AK）、有机质（OM）；另一部分土样-80 保存进行土壤微生物分析。1.4.2籽粒产量测定于燕麦成熟期，收获测定各小区的种子产量。1.4.3土壤微生物总 DNA 提取及生物信息学分表 1 各处理的肥料施用时期和施用量Tab.1 The fertilizer application period and dosage of each treatment处理TreatmentCKT1T2T3T4T5T6CO（NH2）2/（kg/hm2）播种前Pre-planting010897.286.475.664.854拔节期Jointing stage07264.857.650.443.236有机肥/（m3/hm2）Organic fertilizer播种前Pre-planting0018365472901197山西农业科学 2023 年第 51 卷第 10 期析使用 MN NucleoSpin 96 Soil（Germany）试剂盒 进 行 DNA 的 提 取，每 个 样 品 重 复 3 次，采 用PacBio 测 序 平 台 的 通 用 引 物（正 向 27F 和 反 向1492R），以细菌 16SrRNA 基因全长序列（V3-V4）为目的片段进行 PCR 扩增，扩增体系参照艾锋等17，扩增合格的样品在北京百迈克生物有限公司进行测序分析，对得到的高质量 CCS 序列进行生物信息学分析，具体方法参照艾锋等17文献。最后，使用 PICRUSt软件进行细菌功能类群预测18。此外，利用 STAMP和 R软件进行方差分析等。1.5数据分析利用 Microsoft Excel 2009进行数据整理，采用SPSS 22.0软件用单因素方差分析不同处理间的差异显著性（P0.05），采用Pearson法进行相关性分析。2 结果与分析2.1不同处理下土壤养分及籽粒产量的变化不同处理下土壤特性和籽粒产量如表2所示。由表 2 可知，与 CK 相比，各处理的土壤全氮（TN）和速效氮（AN）含量均有所增加，但各处理间差异不显著；对土壤速效钾（AK）、速效磷（AP）和有机质（OM）含量则有显著影响，其中，速效钾（AK）和速效磷（AP）含量均在 T3 处理时最高，分别为251.14、10.99 mg/kg，与其他处理间差异显著（除T4处理外，P0.05）；随着施入有机肥肥量的增加，土壤中有机质的含量逐渐增加，在 T6处理时达到最大，为 14.86 g/kg，与其他各处理间差异显著（除 T5、T4处理外，P0.05）。此外，籽粒产量 T4处理最大，为4 357.61 kg/hm2，与其他处理间差异显著（P0.05）。2.2土壤微生物特征的变化2.2.1土壤细菌群落多样性分析本研究计算了减氮配施有机肥处理下燕麦根际土壤细菌群落的Alpha多样性。丰富度指数评估（图 1）发现，Chao1和 ACE 指数具有相似的变化规律，表现为 T6处理高于其他处理；多样性评估表明，T3 处理 Shannon多样性指数最高（8.78），与最低的 T2 处理间差异显著（P0.05）；T1处理的 Simpson指数最高，仍以T2处理的最低，二者间差异显著（P0.05）。表 2 不同处理下土壤特性和籽粒产量比较Tab.2 Comparison of soil characteristics and seed yield under different treatments处理TreatmentCKT1T2T3T4T5T6TN/（g/kg）1.010.22a1.250.21a1.510.09a1.520.28a1.470.16a1.560.26a1.400.17aAN/（mg/kg）162.469.24a165.6210.76a157.672.13a174.3214.32a167.263.69a159.326.90a150.365.98aAK/（mg/kg）178.3621.02c208.3445.36bc199.3425.31bc251.1456.34a248.1216.98a209.4615.24b228.9930.67bAP/（mg/kg）8.590.69b9.980.25b9.361.36b10.990.96a9.230.90ab8.211.24b8.760.06bOM/（g/kg）10.180.19e11.251.26d12.300.64c12.851.39bc13.901.28ab14.851.22a14.861.35a籽粒产量/（kg/hm2）Seed yield3 001.0093.12c3 957.00101.23b3 818.00110.26b4 026.78102.13b4 357.61102.13a3 966.0093.04b3 899.42102.13b注：不同小写字母表示处理间差异显著（P0.05）。表中数值为平均值标准差。Note:Different lowercase letters indicated significant difference between treatments(P0.05).The values in the table were the mean standard deviation.图 1不同处理下土壤细菌群落 Alpha多样性Fig.1Alpha diversity of soil bacterial communities under different treatments 1198郑敏娜等：减氮配施有机肥对燕麦土壤细菌群落特征影响与生态功能预测采用主坐标分析（Principal co-ordinates analysis，PCoA）结果表明（图 2），T5、T6 处理落在第 1 象限，T4 处理落在第 1 象限和第 4 象限，T3 处理落在第 1 象限和第 2 象限，T1 和 CK 落在第 2 象限，而 T2 处理则落在第 3 象限，与其他处理间无交集产生，且 T2 处理与其他 6 个处理间存在显著差异（P0.05）。2.2.2细菌群落组成分析结果显示（图 3-A），相对丰度占比较高的是变形菌门（Proteobacteria）、厚壁菌门（Firmicutes）、拟杆菌门（Bacteroidetes）、放线菌门（Actinobacteria）、绿弯菌门（Chloroflexi），它们占细菌群落相对丰度的 76.6%以上。在这些占主要优势的门中，CK、T1、T2 这 3个处理中的优势菌为变形菌门，它的相对丰度占比最高，分别为 44.21%、45.64%和 43.93%，且 CK、T1、T2 与 T3、T4、T5、T6 处理间差异显著（P0.05）；随着有机肥施入量的增加，变形菌门的相对丰度占比出现下降，厚壁菌门成为优势菌，其中，以 T4 处理中相对丰度占比最高，达到 36.12%，与 CK、T1和 T2处理间差异显著（P0.05）；拟杆菌门作为第三丰富的优势门，在 T2 和 T4 处理中相对丰度较高，分别为 23.16%和 20.64%，与其他 5 个处理间差异显著（P0.05）；放线菌门在 T3处理中最为丰富（17.53%），但与 CK（16.21%）间差异不显著；绿弯菌门则在 T3 处理中相对占比最高（8.26%），与CK（7.88%）间差异不显著，与其他 5个处理间差异显著（P1%）有不可培养c-Subgroup-6（uncultured-bacterium-c-Subgroup-6）、鞘 脂 单 胞 菌 属（Sphingomonas）、不可培养bacterium（uncultured-bacterium）、乳酸菌属（Lactobacillus）、拟杆菌属（Bacteroides）等。其中，鞘脂单胞菌属（Sphingomonas）为第 1 优势菌属，在 CK 中相对丰度占比最丰富（5.42%）；乳酸菌属（Lactobacillus）和拟杆菌属（Bacteroides）两类菌属则主要分布在T4、T5、T6处理中，相对丰度显著优于其他4个处理（P0.05）；不可培养c-Subgroup-6（uncultured-bacterium-c-Subgroup-6）则 在 CK、T1处理中相对丰度较高，占比均超过 5.20%。2.2.3不同处理的细菌群落结构与其环境因子的关系在属水平上，冗余分析（RDA）结果表明（图图 2不同处理下土壤细菌群落的 PCA分析Fig.2The PCA analysis of soil bacterial communities under different treatments图 3不同处理下细菌门（A）和属（B）的组成结构Fig.3Composition and structure of bacterial phyla（A）and genus（B）under different treatments1199山西农业科学 2023 年第 51 卷第 10 期4），第 1种和第 2种 RDA 成分能够累积解释细菌菌群 30.11%的变化，且不同环境因子对属水平上细菌群落的影响有所差异，其中，AN 与 Lamia、鞘脂单胞菌属（Sphingomonas）相对丰度呈显著正相关（P=0.004），而与 WCHB1-32、假单胞菌属（Pseudomonas）呈显著负相关（P=0.003）；OM、AK、TN与红球菌属（Rhodococcus）、芽孢杆菌属（Bacillus）、巨单胞菌属（Megamonas）、乳酸菌属（Lactobacillus）相 对 丰 度 呈 显 著 正 相 关（P=0.003），亦 与WCHB1-32呈显著负相关（P=0.004）。2.3土壤细菌 PICRUSt 基因功能预测为了解减氮配施有机肥处理下沙壤土细菌代谢功能，使用 PICRUSt 软件比对测序数据获得的物种组成信息，进行功能基因预测（图 5-A），经与KEGG 数据库比对发现，不同肥料处理下的土壤菌群基因序列，在一级功能代谢通路下注释到的功能主要分为新陈代谢（Metabolism）、环境信息处理（Environment information processing）、遗传信息处理（Genetic information processing）、细 胞 转 化（Cellular processes）、人类疾病（Human diseases）和生物体系统（Organismal systems）等 6 类，其余为PICRUSt目前无法鉴定的细菌类群。与不施肥处理（CK）相比，细胞转化功能基因的丰度表现为T3、T4、T5、T6 处理有显著降低的趋势，而遗传信息处理和生物体系统功能基因的丰度表现为 T4、T5、T6 处理有显著增加的趋势，但新陈代谢、环境信息处理和人类疾病功能基因的丰度在各处理间无显著差异。在二级功能预测中（图 5-B），不同肥料处理下的土壤菌群差异基因主要集中在 10 个通路中，其中，全 局 和 概 览 通 路（Global and overiew maps，41.17%43.51%）、碳水化合物代谢（Carbohydrate metabolism，8.26%9.71%）、氨基酸代谢（Amino acid metabolism，6.27%8.01%）、能 量 代 谢（Energy metabolism，3.95%5.88%）、辅助因子和维生素代谢（Metabolism of cofactors and vitamins，3.37%5.09%）为 主 要 子 功 能。与 不 施 肥 处 理（CK）相比，各处理的氨基酸代谢和辅因子和维生素代谢功能基因显著增加（除 TI、T2外，P0.05），氨基酸代谢功能显著降低（除 T2 外，P0.05），而其他功能基因未有显著差异。图 4不同处理下环境因子和土壤细菌的相对丰度间的 RDA分析Fig.4Redundancy analysis（RDA）of environmental factors and relative abundance of soil bacteria under different treatments 图 5不同施肥处理下土壤细菌功能基因 KEGG丰度（一级功能层（A）和二级功能层（B）Fig.5KEGG abundance of soil bacterial functional genes under different fertilization treatments（primary functional layer（A）and secondary functional layer（B）1200郑敏娜等：减氮配施有机肥对燕麦土壤细菌群落特征影响与生态功能预测3 结论与讨论在燕麦生产中，施肥是向农田系统输入养分元素的最主要途径，也是人为调控、管理养分不可替代的技术手段。但传统的燕麦田施肥是以作物高产为唯一目标，年复一年投入大量的化学肥料，不但造成营养元素利用率较低，而且还引起土壤有机质含量减少、土壤污染等问题。因此，科学、合理的肥料配施（有机肥和无机肥）有助于燕麦田养分高效利用。本研究发现，随着施入有机肥肥量的增加，有机质含量明显增加，与刘佳遥等18、商丽荣等19、XU 等20的研究结果基本一致。此外，T4 处理下籽粒产量达到最大，这也说明在该处理下不但可以保证燕麦正常生长发育所需的氮素，而且对土壤生态最有利。土壤生态系统中所有的微生物种类以及它们与环境之间相互作用的多样化程度被称为土壤微生物多样性9，包括了物种多样性、功能多样性、结构多样性及遗传多样性等不同层面8。其中，土壤生态系统中微生物的物种丰富度和均一度，是微生物 多 样 性 的 最 直 接 表 现 形 式21。本 研 究 表 明，Chao1 和 ACE 指数表现为 T6 处理最高，与 T4、T5处理间差异不显著。此外，肥料措施作为对土壤扰动的一种重要方式，也将会给微生物群落结构和成分带来显著的影响4,17。在本研究中，变形菌门（Proteobacteria）、厚壁菌门（Firmicutes）、拟杆菌门（Bacteroidota）、放线菌门（Actinobacteriota）、绿弯菌门（Chloroflexi）等 5 个菌门最具优势，它们占细菌群落相对丰度的 76.6%以上，这与前人的研究基本一致22-24。而且不同减氮配施比例引起的土壤环境因子变化对不同微生物群体的影响各不相同17。本试验 RDA 分析结果表明，AN、OM、AK、TN 对处理中的细菌群落有着更为显著的影响作用，例如，OM、AK、TN 与红球菌属（Rhodococcus）、芽孢杆菌属（Bacillus）、巨单胞菌属（Megamonas）、乳酸菌属（Lactobacillus）相对丰度呈显著正相关。PICRUSt 功能预测表明，燕麦田减氮配施有机肥后土壤细菌群落的主要功能是新陈代谢，表明代谢是细菌群落主要的核心功能，同时具有环境信息处理、遗传信息处理、细胞转化、人类疾病和生物体系统等 5 个功能，其中，减氮增施有机肥提高了遗传信息处理和生物体系统功能基因丰度，这可能是由于有机肥添加显著富集了与碳循环相关功能基因，有助于提升有机质的转化及碳源的利用率，后续将结合宏基因组学测序进一步分析燕麦田减氮配施有机肥后土壤细菌的群落功能。综上所述，减氮配施有机肥处理下，随着有机肥比例的增加，土壤中有机质的含量逐渐增加，在T6处理时达到最大，籽粒产量在 T4处理下达到最大。与不施有机肥土壤相比，增施有机肥后显著提升了变形菌门、厚壁菌门的相对丰度。并且减氮配施有机肥后明显降低了土壤细菌细胞转化功能基因 的 丰 度，明 显 增 加 了 生 物 体 系 统 功 能 基 因 的丰度。参考文献：1 赵秀芳，戎郁萍，赵来喜.我国燕麦种质资源的收集和评价J.草业科学，2007，24（3）：36-40.ZHAO X F，RONG Y P，ZHAO L X.The collection and evaluation of oat（Avena sative）in ChinaJ.Pratacultural Science，2007，24（3）：36-40.2 陈宝书.牧草饲料作物栽培学M.北京：中国农业出版社，2001：10-23.CHEN B S.Forage feed crop cultivationM.Beijing：China Agriculture Press，2001：10-23.3 李刚，郑敏娜，李荫藩.饲用燕麦品种在晋北农牧交错区的生产性能和营养价值研究J.中国农业科技导报，2021，23（12）：42-53.LI G，ZHENG M N，LI Y F.Comprehensive evaluation of production performance and nutritional value of forage oat varieties in northern Shanxi provinceJ.Journal of Agricultural Science and Technology，2021，23（12）：42-53.4 靳建刚，田再芳，郑敏娜，等.不同施肥措施对饲用燕麦土壤细菌群落多样性的影响J.中国农业科技导报，2023，25（3）：152-160.JIN J G，TIAN Z F，ZHENG M N，et al.Effect of different fertilization measures on the diversity of soil bacteria communities in fed oats（Avena sativa L.）J.Journal of Agricultural Science and Technology，2023，25（3）：152-160.5 ZHANG Y，JI G D，WU T，et al.Urbanization significantly impacts the connectivity of soil microbes involved in nitrogen dynamics at a watershed scaleJ.Environmental Pollution，2020，258：113708.6 WANG H T，MARSHALL C W，CHENG M Y，et al.Changes in land use driven by urbanization impact nitrogen cycling and the microbial community composition in soilsJ.Scientific Reports，2017，7：44049.7 CREAMER R E，HANNULA S E，VAN LEEUWEN J P，et al.Ecological network analysis reveals the inter-connection between soil biodiversity and ecosystem function as affected by land use across EuropeJ.Applied Soil Ecology，2016，97：112-124.8 黄昌勇，徐建明.土壤学M.北京：中国农业出版社，1999.HUANG C Y，XU J M.Pedology M.Beijing：China Agriculture Press，1999.9 林先贵.土壤微生物研究原理与方法M.北京：高等教育出版社，2010.1201山西农业科学 2023 年第 51 卷第 10 期LIN X G.Principles and methods of soil microbiology researchM.Beijing：Higher Education Press，2010.10 侯龙鱼，朱泽义，杨杰，等.我国饲草用燕麦现状、问题和潜力J.西南民族大学学报（自然科学版），2019，45（3）：248-253.HOU L Y，ZHU Z Y，YANG J，et al.Current status，problems and potentials of forage oat in ChinaJ.Journal of Southwest Minzu University（Natural Science Edition），2019，45（3）：248-253.11 VOURLITIS G L，ZORBA G，PASQUINI S C，et al.Chronic nitrogen deposition enhances nitrogen mineralization potential of semiarid shrubland soilsJ.Soil Science Society of America Journal，2007，71（3）：836-842.12 ZEGLIN L H，STURSOVA M，SINSABAUGH R L，et al.Microbial responses to nitrogen addition in three contrasting grassland ecosystemsJ.Oecologia，2007，154（2）：349-359.13 于铁峰.西北干旱灌区紫花苜蓿施肥模型构建及养分高效生理机制研究D.兰州：甘肃农业大学，2018：13-36.YU T F.Study on the construction of alfalfa fertilizer model and physiological mechanism of nutrient efficiency in the arid irrigation area of northwest ChinaD.Lanzhou：Gansu Agricultural University，2018：13-36.14 LEE Y M，JIN Y K，SHIN S C.Complete genome sequence of Antarcticibacterium flavum JB01H24T from an Antarctic marine sedimentJ.Marine Genomics，2020，50：100695.15 SUN J，ZHANG Q，ZHOU J，et al.Pyrosequencing technology reveals the impact of different manure doses on the bacterial community in apple rhizosphere soilJ.Applied Soil Ecology，2014，78：28-36.16 鲍士旦.土壤农化分析M.3 版.北京：中国农业出版社，2000.BAO S D.Soil and agricultural chemistry analysisM.3rd Edition.Beijing：China Agriculture Press，2000.17 艾锋，张生，李强，等.蚯蚓肥复配土壤调理剂对盐碱地土壤性质及中科羊草生长的影响J.山西农业科学，2023，51（5）：531-538.AI F，ZHANG S，LI Q，et al.Effects of earthworm fertilizer mixed with soil conditioner on soil properties in saline alkali land and the growth of Leymus chinensisJ.Journal of Shanxi Agricultural Sciences，2023，51（5）：531-538.18 刘佳遥，魏尊苗，程艳，等.轮作豆角对油莎豆生长及土壤微生物多样性的影响J.