大豆连作根际土壤微生物群落构建特征

2024-03-20何德鑫张雪静李佳宜张鹏宇姚兴东谢甫绨

沈阳农业大学学报 2024年1期

何德鑫，赵翔，张雪静，李佳宜，张鹏宇，姚兴东，谢甫绨

（1.沈阳农业大学大豆研究所，沈阳 110161；2.内蒙古自治区农牧业技术推广中心，呼和浩特 010051）

大豆作为一种重要的油料作物，对全世界粮食安全都具有举足轻重的影响[1-2]。近年来受耕地面积缩减、种植结构调整、过度依赖进口等多方面影响，在有些地区大豆连作已经成为一种普遍的种植模式[3-4]。大豆忌连作，连作会导致大豆土壤理化性质失调、土壤微生物群落发生异常改变、加剧土传病害蔓延，继而导致大豆产量和品质下降[5]。

土壤微生物群落对陆地生态系统循环至关重要，如土壤碳循环、氮循环等[6]。近年来研究表明，土壤微生物群落组成及多样性的变化可引起土壤理化性质和土壤酶活性的变化[7-8]。季尚宁等[9]对连作大豆土壤进行杀菌消毒后种植，结果表明杀菌消毒能从根本上消除大豆连作对土壤造成的不利影响，从而保证大豆良好的生长发育和增产效果。这表明土壤微生物群落的改变是导致大豆连作障碍的主要因素之一[4]。大豆连作会使土壤中致病细菌和病原真菌的富集丰度增加，土壤理化性质发生异常变化，土壤微生物多样性减少继而引起大豆的产量和品质的下降[10]。

群落构建过程对于微生物生态学研究起着至关重要的作用，对于解释微生物组成及多样性形成过程有着重要的指导作用[11]。微生物群落构建过程分为确定性和随机性两个过程，二者主要在研究尺度、群落类型等方面存在差异[12]。目前主要是通过βNTI（β-Nearest Taxon Index,最近种间指数）模型来构建土壤微生物群落，并通过βNTI模式的构建发现确定性过程对土壤生态环境起主导作用[13]。郭晗玥等[14]对西瓜连作土壤群落构建过程的研究表明，西瓜连作细菌群落构建过程以随机性为主导，真菌构建过程由确定性过程主导。然而，对于连作大豆根际土壤生态环境中土壤微生物群落的构建过程机制，至今仍缺乏深入细致的探究。

本研究选择连作大豆根际土壤微生物群落为对象，采用16S rRNA 测序、ITS 测序和βNTI 模型构建，分析了连作大豆、玉米-大豆轮作根际土壤微生物群落组成、多样性和群落构建过程的差异，试图阐释连作大豆土壤微生物组装机制，为缓解大豆连作障碍研究提供新的见解。

1 材料与方法

1.1 试验地点和取样

2015年在沈阳农业大学辽中定位试验基地（41°52’N，122°72’W，海拔5.5~23.5 m）建立定位试验，采用大豆连作和玉米-大豆轮作的种植制度进行种植。本试验于2021-2022年进行取样测定。采用大豆连作（LCR）和玉米-大豆轮作（LRR）两种处理方法，玉米-大豆轮作为对照处理，玉米：东单6531，大豆：辽豆14。小区行数6 行，垄长15 m，垄距0.6 m，小区面积54 m2，密度15 万株·hm-2，3 个重复，每穴保苗2株，常规田间管理。

在大豆鼓粒期（R6），按照5点法取样，即从每个处理中随机选择5个取样点，用无菌刷子将大豆根系表面5 mm以内的土壤轻轻取下，作为根际土壤，混合充分，在-80 ℃保存，用于土壤微生物的测序。

每个小区按照5点法取5个点（直径5 cm×深度20 cm），然后充分混合成1个复合样品，并用2 mm 的筛子筛分土壤样品以除去叶子和根等作物残茬，然后每个样分成两部分：1 份保存在-80 ℃用于土壤酶活性的测定，1份自然条件下风干用于土壤化学指标分析。

1.2 土壤化学指标测定

使用pH 计测定土壤pH 值;土壤含水量和容重的测定分别使用铝盒烘干法和环刀法;总碳、总氮使用元素分析仪测定;铵态氮(NH4+-N)和硝态氮(NO3--N)使用间断化学分析仪测定。

1.3 土壤细菌16S和真菌ITS高通量测序

将待测样品送到武汉迈特维尔生物科技有限公司，提取土壤总DNA 和后续分析。采用CATB 方法进行土壤样品DNA 提取；DNA 提取完成后，采用琼脂糖凝胶电泳检测提取的DNA 浓度以及纯度。用无菌水将所取的样品稀释到1 ng·μL-1；稀释后得到DNA 模板，进行PCR 扩增；选用琼脂糖凝胶对扩增后的引物进行电泳检测；对16S r RNA 的V4 区进行PCR 扩增，使用特异性引物515F（GTGCCAGCMGCCGCGGTAA）和806R（GGACTACHVGGGTWTCTAAT）[15]；ITS1 区引物ITS5-1737F(5′-GGAAGTAAAAGTCGTAACAAGG-3′)和ITS2-2043R(5′-GCTGCGTTCTTCATCGATGC-3′)[16]。

1.4 数据分析

2 结果与分析

2.1 土壤环境因子对大豆连作的响应

由表1 可知，两种种植模式下土壤的化学性质呈现不同变化。与玉米-大豆轮作相比，大豆连作土壤全氮、全磷、全钾、全碳、有效磷、硝态氮、铵态氮、有机质和速效磷含量均不同程度下降，其中土壤全氮、全磷、硝态氮含量下降较大，下降幅度分别为40.35%、56.45%、33.13%。玉米-大豆轮作降低了土壤速效钾含量和pH值。

表1 大豆轮作和连作处理对土壤理化性质参数的影响Table 1 Effects of soybean rotation and continuous cropping on soil physicochemical properties parameters

2.2 土壤酶活性对大豆连作的响应

由表2可知，对比玉米-大豆轮作处理，连作显著提高土壤蔗糖酶和过氧化氢酶活性，提高幅度分别为30.15%和54.30%；显著降低土壤脲酶、亮氨酸氨肽酶活性和N-乙酰氨基葡萄糖苷酶活性，下降幅度分别为12.32%、13.60%和20.28。两处理间土壤蔗糖酶、过氧化氢酶、亮氨酸氨肽酶活性及N-乙酰氨基葡萄糖苷酶活性差异显著。

表2 大豆轮作和连作处理对土壤酶活性参数的影响Table 2 Effects of soybean rotation and continuous cropping on soil enzyme activity parameters

2.3 根际土壤微生物多样性和群落结构对大豆连作的响应

由图1a可知，在属水平下（相对丰度>0.01%），玉米-大豆轮作和大豆连作及根际土壤细菌群落的优势属均为鞘氨醇单胞菌属（Sphingomonas，LRR：6%；LCR：9%）、短根瘤菌属（Bradyrhizobium，LRR：3.9%；LCR：4.1%）、贪噬菌属（Variovorax，LRR：0.7%；LCR：2%）、鞘氨醇杆菌属（Sphingobacterium，LRR：0.15%；LCR：1.5%）、噬几丁质菌属（Chitinophaga，LRR：1.9%；LCR：0.9%）、卡氏伯克霍尔德菌属（Burkholderia-Caballeronia-Paraburkholderia，LRR：1.4%；LCR：2.4%）、念珠菌属（Candidatus_Solibacter，LRR：2.5%；LCR：1.9%）、氨氧化菌属（Ellin6067，LRR：2.2%；LCR：2%）、芽单胞菌属（Gemmatimonas，LRR：2.1%；LCR：1.7%）、苔藓杆菌（Bryobacter，LRR：1.7%；LCR：1.5%）。相比于玉米-大豆轮作处理，大豆连作根际土壤鞘氨醇单胞菌属、短根瘤菌属、贪噬菌属、鞘氨醇杆菌属、卡氏伯克霍尔德菌属相对丰度较高。

由图1b可知，在属水平下（相对丰度>0.01%），玉米-大豆轮作和大豆连作根际土壤真菌的优势属均为散尾鬼笔属（Lysurus，LCR：29.5%；LRR：19.3%）、球盖菇属Stropharia，LCR：1.3%；LRR：15.9%）、镰刀菌属（Fusarium，LCR：20.5%；LRR：13.2%）、丽赤壳属（Calonectria，LCR：8.9%；LRR：0.2%）、拟松针藻属（Codinaeopsis，LCR：1.6%；LRR：4.2%）、亚隔孢壳属（Didymella，LCR：4.7%；LRR：1.3%）、球菌属（Powellomyces，LCR：1.3%；LRR：0.1%）、被孢霉属（Mortierella，LCR：2.4%；LRR：1.6%）、线虫捕捉菌属（Arthrobotrys，LCR：0.05%；LRR：1.2%）、鹅膏菌（Amanita，LCR：0.05%；LRR：0.8%）相比于玉米大豆轮作处理，大豆连作根际土壤镰刀菌属(Fusarium)、丽赤壳属(Calonectria)、亚隔孢壳属(Didymella)、被孢霉属(Mortierella)相对丰度较高。

图1 玉米-大豆轮作和大豆连作根际土壤微生物群落的优势属柱状图Figure 1 Dominant genus bar chart of rhizosphere soil microbial communities in maize soybean rotation and continuous soybean cropping

2.4 大豆连作根际土壤微生物群落特征

由表3 和表4 可知，对比轮作根际，连作根际土壤细菌富集OTU 数量显著降低，下降幅度为8.21%，连作根际土壤细菌香农指数、辛普森指数、chao1 指数和ace 指数均不同程度下降，下降幅度分别为15.33%、10.98%、6.28%和5.92。同样地，对比轮作根际，连作根际土壤真菌OTU数量、香农指数、chao1指数和ace指数均不同程度下降，但下降幅度较小。

表3 轮作和连作处理对大豆土壤细菌α多样性的影响Table 3 Rotation and continuous cropping treatments on soybean soil bacteriaα impact of diversity

表4 轮作和连作处理对大豆土壤真菌α多样性的影响Table 4 Effects of rotation and continuous cropping on fungi in soybean soil α impact of diversity

2.5 连作下土壤环境因子与土壤微生物的关系

冗余分析（RDA）表明，土壤环境因子与土壤微生物相对丰度具有较大的相关性（图4）。两处理处理下10个细菌属和真菌属表现出差异。细菌RDA 分布结构中，RDA 前两轴分别解释38.1%（RDA1）、23.5%（RDA2）。其中速效钾、速效磷、有机质、有机碳、全钾、蔗糖酶、过氧化氢酶、脲酶、亮氨酸氨肽酶及N-乙酰氨基葡萄糖苷酶是决定罗河杆菌属（Rhodanobacter）、慢生根瘤菌属（Bradyrhizobium）、鞘脂单胞菌属（Sphingomonas）、多嗜菌（Variovorax）、伯克霍尔德菌属（Burkholderia）、鞘脂杆菌属（Sphingobacterium）、芽单胞菌属（Gemmatimonas）和鞘氨醇菌属（Chitinophaga）相对丰度的主要土壤酶。其中，LRR 处理下，脲酶、亮氨酸氨肽酶、N-乙酰氨基葡萄糖苷酶、速效磷、有机质、有机碳、全钾与罗河杆菌属（Rhodanobacter）及芽单胞菌属（Gemmatimonas）表现出强烈相关性；LCR 处理下，速效钾、蔗糖酶、过氧化氢酶及慢生根瘤菌属（Bradyrhizobium）具有较强的正相关性（图3a、图3c）。真菌RDA 分布结构中，RDA 前两轴分别解释49.3%（RDA1）和28.8%（RDA2）。其中，速效钾、速效磷、有机质、有机碳、全钾、蔗糖酶、过氧化氢酶、脲酶、亮氨酸氨肽酶及N-乙酰氨基葡萄糖苷酶是决定散尾鬼笔属（Lysurus）、镰刀菌属（Fusarium）、球盖菇属（Stropharia）、被孢霉属（Mortierella）、丽赤壳属（Calonectria）、巨孢囊霉属（Gigaspora）、克鲁维酵母属（Kluyveromyces）、亚隔孢壳属（Didymella）、丝葚霉属（Papulaspora）和毛霉菌属（Mucoromycotina）相对丰度的主要土壤因子。LRR处理下，脲酶、亮氨酸氨肽酶、N-乙酰氨基葡萄糖苷酶、速效磷、有机质、有机碳、全钾与球盖菇属（Stropharia）、丝葚霉属（Papulaspora）及毛霉菌属（Mucoromycotina）表现出强烈相关性；LCR 处理下，速效钾、蔗糖酶、过氧化氢酶与镰刀菌属（Fusarium）、亚隔孢壳属（Didymella）、丽赤壳属（Calonectria）及被孢霉属（Mortierella）具有较强的正相关性（图3b、图3d）。

图3 轮作和连作处理根际环境因子与土壤微生物属水平的RDA分析Figure 3 RDA analysis of rhizosphere environmental factors and soil microbial genus levels in crop rotation and continuous cropping treatments

2.6 土壤微生物群落构建过程

图4 轮作和连作处理根际土壤微生物群落构建图Figure 4 Construction diagram of rhizosphere soil microbial community under rotation and continuous cropping treatments

3 讨论与结论

3.1 大豆连作、玉米-大豆轮作对根际土壤化学性质和酶活性的影响

大豆连作可造成土壤理化性质失衡[17]。玉米-大豆轮作可显著提高大豆根际土壤养分含量[10]。本试验研究结果显示，玉米-大豆轮作处理下，土壤全氮、溶解性有机碳、有机质等含量均高于连作处理。土壤酶活性随着土壤内部环境的变化而变化，如土壤过氧化氢酶与土壤逆境相关联。土壤脲酶活性变化可引起土壤养分的变化。亮氨酸氨肽酶是一种土壤氮获取酶，其活性变化与微生物氮利用效率呈正相关[18-20]。本研究结果表明，连作会提高土壤蔗糖酶和过氧化氢酶活性，降低土壤脲酶、亮氨酸氨肽酶活性和N-乙酰氨基葡萄糖苷酶活性，该结果与LIU等[18]研究结果一致。

3.2 大豆连作、玉米-大豆轮作对根际土壤微生物群落的影响

不同种植模式下，土壤微生物群落的组成、数量和多样性变化不尽相同。RAO 等[10]研究发现，与大豆连作相比，玉米-大豆轮作提高了大豆根际土壤微生物的数量和多样性指数。本研究中，相比于大豆连作，玉米-大豆轮作会提高大豆根际土壤细菌、真菌OTUs数量、ace指数和shannon指数。

相关研究表明，作物连作可造成其根际土壤微生物群落组成的变化[10]。本研究中，大豆连作土壤中鞘氨醇杆菌属（Sphingobacterium）和伯克霍尔德菌属（Burkholderia）相对丰度高于轮作处理，而轮作土壤中鞘氨醇菌属（Chitinophaga）相对丰度高于连作处理。鞘氨醇杆菌属（Sphingobacterium）和伯克霍尔德菌属（Burkholderia）是一种生防细菌，其耐受恶劣土壤生态环境[21-22]。本研究中，连作处理鞘氨醇杆菌属（Sphingobacterium）和伯克霍尔德菌属（Burkholderia）相对丰度升高，有可能是连作障碍激发了大豆根际鞘氨醇杆菌属（Sphingobacterium）和伯克霍尔德菌属（Burkholderia）生物防御效应以应对大豆连作障碍造成的。小麦连作障碍的研究表明，噬几丁质菌属（Chitinophaga）具有抑制或拮抗功能，可抑制病原菌在土壤中繁殖，其相对丰度的增加可进一步抑制土传病害的发生，进而促进植株健康生长[23]。在真菌群落中，有研究表明，球盖菇属（Stropharia）为有益真菌，球盖菇属（Stropharia）真菌相对丰度升高，可有效改善表层土壤理化性质，提升土壤有机质含量，减少土壤病虫害[24]。镰刀菌属（Fusarium）真菌为连作土壤典型病原菌，其相对丰度的升高导致党参、西瓜枯萎病的发生[25-26]。本研究中，与轮作处理相比，连作处理镰刀菌属（Fusarium）相对丰度显著增加，球盖菇属（Stropharia）真菌相对丰度显著降低。因此，连作改变土壤微生物群落组成，增加病原菌相对丰度，从而引起连作障碍的发生。

3.3 大豆连作、玉米-大豆轮作土壤因子与土壤微生物的关系

前人研究表明，土壤环境因子与土壤微生物相对丰度具有密切联系[27]。本研究中，轮作处理下，脲酶、亮氨酸氨肽酶、N-乙酰氨基葡萄糖苷酶、速效磷、有机质、有机碳及全钾显著影响罗河杆菌属（Rhodanobacter）和芽单胞菌属（Gemmatimonas）细菌；球盖菇属（Stropharia）、丝葚霉属（Papulaspora）和毛霉菌属（Mucoromycotina）真菌相对丰度；连作处理下，速效钾、蔗糖酶、过氧化氢酶显著影响细菌慢生根瘤菌属（Bradyrhizobium）、真菌镰刀菌（Fusarium）、亚隔孢壳属（Didymella）及丽赤壳属（Calonectria）、被孢霉属（Mortierella）相对丰度。镰刀菌（Fusarium）与土壤蔗糖酶具有显著正相关性，蔗糖酶进一步分解可为镰刀菌提供有利碳源[28]。本研究中，轮作处理蔗糖酶含量高于连作处理，轮作处理镰刀菌（Fusarium）相对丰度低于连作处理，这也说明，大豆连作条件下，土壤消耗了过量的蔗糖酶等土壤因子从而提升了镰刀菌（Fusarium）等致病真菌的相对丰度。

3.4 大豆连作、玉米-大豆轮作对根际土壤微生物群落组装过程的影响

为了更好地阐释土壤根际微生物群落构建对大豆连作的响应，研究了连作大豆、玉米-大豆轮作土壤根际细菌、真菌群落构建过程。本研究中，由于土壤化学性质（pH、土壤全氮、全碳）和土壤酶活性等变化，使得大豆连作根际土壤细菌群落构建主要由随机性过程转为确定性过程；玉米-大豆轮作根际土壤细菌群落构建主要由随机过程主导；大豆连作根际土壤真菌群落构建由随机过程主导[29]。有研究指出，西瓜连作时，pH 值和有机质含量等环境因素的选择会将细菌群落构建过程由随机性转换为确定性过程，而在连作年限增加后，细菌群落构建更容易受到随机性事件影响[14]。真菌群落构建过程中，随机性过程更为突出，这是由于一些土壤环境因子如根系分泌物、土壤理化指标、土壤酶等变化所造成的结果[30]。

连作会显著影响大豆根际土壤细菌和真菌群落结构，从而导致土壤有益菌（噬几丁质菌属Chitinophaga、球盖菇属Stropharia）相对丰度下降，病原菌（镰刀菌Fusarium、伯克霍尔德菌属Burkholderia）相对丰度升高，进而改变土壤微生物群落的组成和多样性，降低根际土壤应对环境变化的抵抗作用。玉米-大豆轮作模式下，土壤真菌相对丰度显著下降，根际土壤细菌群落组装也趋向于随机性过程。本研究结果表明，玉米-大豆轮作模式通过改变大豆根际土壤群落特征和提高土壤有益菌相对丰度等方式，缓解大豆连作障碍，进而促进大豆土壤环境健康，提升大豆产量。