西兰花残体还田对棉花黄萎病防治效果及其对不同生育时期土壤细菌群落的影响

赵卫松，郭庆港，李社增，王亚娇，鹿秀云，王培培，苏振贺，张晓云，马平

西兰花残体还田对棉花黄萎病防治效果及其对不同生育时期土壤细菌群落的影响

赵卫松，郭庆港，李社增，王亚娇，鹿秀云，王培培，苏振贺，张晓云，马平

（河北省农林科学院植物保护研究所/河北省农业有害生物综合防治工程技术研究中心/农业农村部华北北部作物有害生物综合治理重点实验室， 河北保定 071000）

【目的】研究土壤中添加西兰花残体（broccoli residues，BR）对棉花黄萎病及不同生育时期土壤细菌群落结构的影响，为棉花黄萎病的绿色生态防控和化学农药的减量施用提供新途径和思路。【方法】通过田间小区试验，分别以西兰花残体还田和不含有西兰花残体（CK）的耕层土壤中播种的棉花为试验对象，监测不同处理的棉花黄萎病的发生动态。采用实时荧光定量PCR（real-time PCR）和高通量测序（Illumina MiSeq）技术分别测定不同时期土壤中大丽轮枝菌（）DNA拷贝数量和土壤细菌群落，进而分析西兰花残体还田对土壤中病原菌数量变化和土壤细菌群落的影响。采用主成分分析探索西兰花残体还田对不同生育时期土壤细菌群落在门和属水平的结构变化特征与规律。【结果】西兰花残体还田后棉花黄萎病的发病率和病情指数均呈现出下降趋势，且发病高峰时间推迟，在发病高峰期对黄萎病的防治效果达到70.77%，建立了整个生育期内病情指数-时间病程发展曲线，平均防治效果为57.21%。与空白对照相比，西兰花残体还田处理在棉花生长季和非生长季的4个阶段土壤中大丽轮枝菌DNA拷贝数量分别降低了10.96%、11.11%、25.95%和11.25%。高通量测序分析表明，不同生育时期西兰花残体还田处理显著提高了土壤细菌多样性，丰富度指数Chao1均显著上升，ACE指数在播种前期、花铃期和吐絮期显著上升。与空白对照相比，在不同生育时期内，西兰花残体还田后土壤中的放线菌门（Actinobacteria）菌群的相对丰度呈现显著上升。同时研究表明，西兰花残体还田后优势菌群组成在不同生育时期存在差别，在播种前期和苗期，厚壁菌门（Firmicutes）菌群上升为优势菌群；在蕾期和花铃期，蓝藻细菌门（Cyanobacteria）菌群成为优势菌群。主成分分析表明，西兰花残体还田改变了苗期、蕾期、花铃期和吐絮期的细菌群落结构。进一步分析表明，西兰花残体还田显著增加了链霉菌属（）和芽孢杆菌属（）菌群的相对丰度。【结论】西兰花残体还田能够有效抑制棉花黄萎病的发生，降低土壤中大丽轮枝菌DNA拷贝的数量，改变土壤细菌群落结构，增加有益微生物的相对丰度，是一种有效控制棉花黄萎病的绿色措施。

西兰花残体；大丽轮枝菌；棉花黄萎病；细菌群落结构；高通量测序；生育时期

0 引言

【研究意义】棉花黄萎病（cotton verticillium wilt）是由大丽轮枝菌（）侵染引起的真菌性维管束病害，是导致棉花减产和纤维品质下降的主要病害[1-2]，已成为限制我国棉花产业可持续发展的重要瓶颈。目前，有关棉花黄萎病防治的研究主要集中在抗性品种、轮作及微生物菌剂和化学药剂应用等方面[1,3-5]。由于大丽轮枝菌具有致病力分化严重、变异速度快、传播途径多和在土壤中存活时间长等特点[2,6]，导致抗病品种出现抗性丧失的现象[7-8]；过去防治土传病害主要依靠化学药剂进行熏蒸处理，不仅降低了棉花的品质，还破坏了土壤的生态系统，影响环境。因此，亟需寻找环境友好型的黄萎病防控技术。【前人研究进展】利用植物及其残体是有效防治作物土传病害的途径之一[9-10]。西兰花（var.）残体（broccoli residues，BR）还田能够有效降低花椰菜[11-12]、茄子[13-14]、马铃薯[15-16]、草莓[17]等黄萎病的发生，Shetty等[11]研究表明，西兰花及其他部分十字花科作物体内含有硫代葡萄糖甙，其在芥子酶水解作用下降解为异硫氰酸酯，从而可控制土传病害的发生；许多学者研究发现，西兰花残体还田后降低了茄子、花椰菜和马铃薯黄萎病的发生，其防病机理主要是通过向环境中释放具有抑菌功能的次生代谢物质，改变土壤微生物多样性，降低土壤病原菌或微菌核数量、增加土壤中有益微生物菌群丰度[12-16]。近年来，随着测序技术的快速发展，高通量测序技术具有测序深度深和获得数据量大等优点，能更真实地揭示微生物群落的复杂性和多样性。【本研究切入点】关于西兰花残体还田对棉花黄萎病的防控研究尚少，笔者课题组前期研究发现，土壤中添加一定比例的西兰花残体对棉花黄萎病具有较好的防治效果[18]，但关于西兰花残体对棉花黄萎病及其不同生育时期土壤细菌群落结构的影响尚缺乏系统的研究。【拟解决的关键问题】以西兰花残体还田深翻后对棉花黄萎病的防治效果为研究主线，以棉花不同生育时期的土壤为研究对象，采用Illumina MiSeq高通量测序技术，分析西兰花残体还田后棉花不同生育时期土壤细菌多样性和群落组成的变化，揭示土壤微生物对西兰花残体处理的响应机制，阐释棉花黄萎病与土壤微生物间的互作关系，为深入研究西兰花残体防治棉花黄萎病的生态机制提供理论依据和技术指导，为棉花黄萎病的绿色防控和化学农药的减量施用提供新途径和思路。

1 材料与方法

1.1 试验材料

以西兰花品种‘炎秀’、棉花品种‘邯686’（耐黄萎病品种）为供试材料，由河北省邯郸市农业科学院提供，用于田间试验。

1.2 试验区域概况、试验设计与土壤样品采集

河北省邯郸市成安县位于冀南平原，东经114°29'—114°53'，北纬36°18'—36°30'，是全国棉花生产百强县、优质棉基地县。20世纪90年代初，该地区棉花黄萎病重病田发病率达89.9%，之后利用抗病品种、农业措施（如轮作）等手段，在一定程度上达到了防治黄萎病的目的，但该病原菌主要由土壤和种子传播，寄主范围广，流行性强，近年来常有发生的报道。

2017年3—11月在成安县道东堡乡大堤西村进行西兰花残体还田防治棉花黄萎病的试验（N：36º28′17″，E：114º45′55″）。选取地势平坦、肥力相对均匀、连续多年种植棉花的地块作试验田。2017年3月上旬，将西兰花幼苗（三叶一心）移栽到田间，行距70 cm，株距35 cm，种植密度约为2 720株/666.7 m2；2017年5月中旬收获西兰花，之后将剩余的植株地上部分和地下部分经粉碎机粉碎、旋耕机深翻25—30 cm作为处理区（T，施入西兰花残体约为3 500—3 800 kg/666.7 m2），未种植西兰花的地块作为空白对照区（CK）。2017年6月1日分别在处理区和对照区播种棉花，每个处理3次重复，每小区面积为45 m2（行长60 m，每行宽0.75 m），行距75 cm，株距35 cm。在棉花生长期间，按照当地生产技术进行水肥、植株虫草害管理。分别在播种前期、苗期、蕾期、花铃期和吐絮期，按照“W”形采用多点取样法，用取土器在棉花行间距离棉花根部10 cm采集0—20 cm耕作层的土壤，每小区采集约500 g土壤混匀装于自封袋中，于-80℃冰箱保存，用于后续土壤细菌群落分析。

1.3 测定项目及方法

1.3.1 黄萎病发生调查 在黄萎病发生时期，按照李社增等[19]分级标准于8月10日、8月21日、8月31日、9月12日和9月21日分别调查统计发病率和病情指数，分别记作DIT1、DIT2、DIT3、DIT4和DIT5。病情调查标准如下，0级：棉株健康，无病叶，生长正常；1级：棉株1/4以下叶片发病，变黄萎蔫；2级：棉株1/4以上，1/2以下叶片发病，变黄萎蔫；3级：棉株1/2以上，3/4以下叶片发病，变黄萎蔫；4级：棉株3/4以上叶片发病，或叶片全部脱落，棉株枯死。

为明确西兰花残体还田在整个棉花生育期内对黄萎病的防治效果，建立病情指数-时间的病程发展曲线图并计算病程发展曲线下面积（area under disease progress curve，AUDPC），进一步计算平均防治效果。

发病率（%）=（发病株数/调查总株数）×100；病情指数（DI）=∑（相应病级×各级病株数）/（最高病级×总株数）×100；AUDPC=1/2×[（DIT1+DIT2）/（T2-T1+1）+（DIT2+DIT3）/（T3-T2+1）][18]；平均防治效果（%）=（对照的AUDPC-处理的AUDPC）/对照的AUDPC×100。

1.3.2 土壤细菌群落结构的高通量测序 采用FastDNATMSPIN Kit for Soil（MP Biomedicals，Solon，OH，USA）试剂盒提取土壤基因组DNA，利用NanoDrop 2000分光光度计（Thermo Fisher Scientific Inc.，Waltham，MA，USA）检测DNA的纯度和浓度。DNA于-20℃保存。

PCR扩增：以上述基因组DNA为模板，土壤细菌16S rDNA的V3-4区扩增采用引物为338F（5′-ACT CCTACGGGAGGCAGCA-3′）和806R（5′-GGACTAC HVGGGTWTCTAAT-3′）[20]，根据序列中的保守区域设计相应引物，并添加样本特异性Barcode序列，采用NEB公司的Q5高保真DNA聚合酶进行PCR扩增。PCR扩增产物通过2%琼脂糖凝胶电泳进行检测，采用AXYGEN公司的凝胶回收试剂盒对目标片段进行切胶回收。

文库构建和测序：由上海派森诺生物科技有限公司完成。采用Illumina公司的TruSeq Nano DNA LT Library Prep Kit制备测序文库，采用Agilent High Sensitivity DNA Kit在Agilent Bioanalyzer上对文库进行质检，并采用Quant-iT PicoGreen dsDNA Assay Kit在Promega QuantiFluor荧光定量系统上对文库进行定量，合格后，使用MiSeq高通量测序技术进行上机测序。

1.3.3 土壤中大丽轮枝菌数量的测定 DNA的提取、重组质粒的制备和鉴定：采用CTAB法提取大丽轮枝菌的基因组DNA。利用大丽轮枝菌特异性引物DB19（CGGTGACATAATACTGAGAG）/DB20（GACGATGCGGATTGAACGAA）扩增ITS部分序列。将PCR产物与pMD19-T Vector（宝生物工程（大连）有限公司）16℃连接4 h后，转化于DH5感受态细胞中，37℃过夜培养，挑选白色菌落转接于含100 μg·mL-1氨苄青霉素的LB液体培养基中，37℃，200 r/min振荡培养12 h。提取重组质粒，利用NanoDrop 2000分光光度计测定浓度，计算单位体积质粒所含的DNA拷贝数。将重组质粒送生工生物工程（上海）股份有限进行插入片段的序列测定。测序结果显示，重组质粒中插入片段序列与大丽轮枝菌ITS序列同源性为100%。

质粒拷贝数=

实时荧光定量PCR反应及标准曲线的建立：将提取的质粒10倍梯度稀释，将其作为模板进行实时荧光定量PCR扩增。反应体系：2×Master Mix（德国DBI Bioscience公司）10 μL，大丽轮枝菌ITS上下游引物[21]（F：CCCGCCGGTCCATCAGTCTCTCTG；R：CGGGACTCCGATGCGAGCTGTAAC）各1 μL，DNA 1 μL，ROX 0.4 μL，ddH2O补至20 μL。扩增程序：95℃预变性30 s，95℃变性5 s，60℃退火并延伸30 s，扩增40个循环。以质粒拷贝数的对数值为横坐标，以循环阈值Ct值为纵坐标建立标准曲线。

土壤中大丽轮枝菌数量的测定：在棉花的不同生育时期，按照1.2方法采集耕作层土壤，提取DNA，并将其作为模板进行实时荧光定量PCR反应，每个样品重复3次，得到不同时间样本的循环阈值，根据上述建立的标准曲线计算土壤中大丽轮枝菌的拷贝数。

1.4 数据分析

试验数据采用Microsoft Excel 2010进行整理，Origin 8.6进行作图，采用SPSS 17.0软件进行单因素方差分析、独立样本T检验和邓肯新复极差法（DMRT）多重比较统计分析。采用Canoco 4.5软件对样本微生物群落进行主成分分析（PCA），并利用CanoDraw软件进行作图；利用软件Mothur计算Alpha多样性指标，包括丰富度指数（ACE和Chao1）和均匀度指数（Shannon和Simpson）。

2 结果

2.1 西兰花残体还田对棉花黄萎病发生的影响

由图1可知，西兰花残体（BR）处理与空白对照（CK）相比，前者降低了黄萎病的发病率和病情指数。在8月底9月初，黄萎病的发生达到高峰，其中西兰花残体还田后发病率和病情指数分别为6.75%和1.80，而空白对照的发病率和病情指数分别为24.62%和6.76。空白对照7月30日开始发病并表现症状，而西兰花残体还田后开始发病时间和发病高峰推迟。

图1 棉花黄萎病发生的动态变化

2.2 全生育期内对棉花黄萎病的防治效果

西兰花残体还田后对棉花黄萎病的田间防治试验结果表明，在调查时间范围内（7月10日至9月21日）不同处理的棉花黄萎病病情指数呈上升趋势。从8月10日开始发病，在发病高峰期（8月31日）西兰花残体处理对黄萎病的防治效果达到70.77%，为了客观评价西兰花残体在整个生育期内对棉花黄萎病的防治效果，建立了病情指数-时间病程发展曲线，平均防治效果为57.21%（图2）。

图2 棉花黄萎病病情指数-时间的AUDPC

2.3 西兰花残体还田对土壤中病原菌数量的影响

利用实时荧光定量PCR建立了定量检测大丽轮枝菌的标准曲线，其方程为=-3.4309+46.756，相关系数2=0.9993，表明建立的标准曲线线性好、扩增效率高。西兰花残体处理后在棉花生长季和非生长季的不同时间土壤中病原菌数量的变化结果表明，在发病初期（7月11日），西兰花残体处理土壤样品中的大丽轮枝菌与空白对照差异不显著，在其他3个时间（8月21日、10月19日和11月28日），西兰花残体处理显著降低了土壤中大丽轮枝菌DNA拷贝数，与空白对照相比，土壤中病原菌的数量分别下降了11.11%、25.95%和11.25%（图3）。

图3 棉花土壤中大丽轮枝菌实时荧光定量PCR测定

2.4 西兰花残体还田对细菌群落多样性的影响

由表1可知，独立样本的t检验结果表明，与空白对照相比，西兰花残体处理后在棉花不同生育时期土壤细菌的Chao1指数均显著上升，ACE指数在播种前期、花铃期和吐絮期显著上升，而Shannon和Simpson指数差异不显著。

2.5 西兰花残体还田对细菌群落结构的影响

对不同生育时期样本的细菌群落结构进行主成分分析，由图4可知，所有样品在主成分坐标中位于不同的象限，说明不同生育时期样品的细菌群落结构存在差异。通过对不同生育时期的样品进一步分析，对于空白对照样品，分别位于4个象限，播种前期和苗期位于相同象限，蕾期、花铃期和吐絮期分别位于不同象限，说明细菌群落结构存在差异。对于西兰花残体处理样品，分别位于4个象限，其中播种前期和花铃期位于相同象限，苗期、蕾期和吐絮期分别位于不同象限。

门水平下西兰花残体还田在不同生育时期对土壤中细菌群落结构的影响结果表明，在播种前期（T0），空白对照处理（CKT0）获得30个菌群，其中有8个优势菌群。西兰花残体处理（BRT0）获得25个菌群，其中有9个优势菌群，相对丰度占总体的95%以上。变形菌门（Proteobacteria）、放线菌门（Actinobacteria）、酸杆菌门（Acidobacteria）、芽单孢菌门（Gemmatimonadetes）、绿弯菌门（Chloroflexi）、浮霉菌门（Planctomycetes）、拟杆菌门（Bacteroidetes）和硝化螺旋菌门（Nitrospirae）为两个处理共同的优势菌群，其相对丰度分别为30.31%和28.53%、21.96%和26.59%、12.99%和13.87%、13.54%和9.07%、7.36%和7.48%、3.96%和3.57%、2.90%和3.14%、2.77%和3.31%，其中西兰花残体处理后放线菌门、酸杆菌门、绿弯菌门、拟杆菌门和硝化螺旋菌门增幅分别为21.08%、6.77%、1.63%、8.28%和19.49%。此外，西兰花残体处理后厚壁菌门（Firmicutes）菌群上升为优势菌群，其相对丰度由0.55%升至2.17%（图5）。

表1 西兰花残体还田对棉花不同生育时期细菌测序量和群落多样性的影响

数据后不同小写字母表示差异显著（＜0.05）。下同

<

Different lowercases after the data indicate significant differences(＜0.05). The same as below

红色椭圆代表空白对照（1—15），绿色椭圆代表西兰花残体处理（16—30）The red ellipse represents the blank control (1-15), the green ellipse represents the broccoli residues treatment (16-30)。1—3、16—18：播种前期Pre-sowing stage；4—6、19—21：苗期Seedling stage；7—9、22—24：蕾期Bud stage；10—12、25—27：花铃期Flowering and boll stage；13—15、28—30：吐絮期Boll-opening stage。Proteobacteria：变形菌门；Actinobacteria：放线菌门；Acidobacteria：酸杆菌门；Gemmatimonadetes：芽单孢菌门；Chloroflexi：绿弯菌门；Planctomycetes：浮霉菌门；Nitrospirae：硝化螺旋菌门；Bacteroidetes：拟杆菌门；Firmicutes：厚壁菌门；Cyanobacteria：蓝藻细菌门；Others：其他

在苗期（T1），空白对照处理（CKT1）和西兰花残体处理（BRT1）分别获得27个菌群，分别得到9个和8个主要优势菌群，相对丰度占总体的97.09%和96.81%。变形菌门、放线菌门、酸杆菌门、芽单孢菌门、绿弯菌门、浮霉菌门、拟杆菌门和硝化螺旋菌门为共同的优势菌群，其相对丰度分别为31.66%和28.59%、20.07%和20.77%、13.20%和16.22%、12.77%和12.84%、7.11%和8.07%、3.02%和3.81%、3.49%和3.54%、4.76%和2.97%，其中西兰花残体处理后放线菌门、酸杆菌门、芽单孢菌门、绿弯菌门、浮霉菌门和拟杆菌门均呈现不同程度的上升，其增幅分别为3.49%、22.88%、0.55%、13.50%、26.16%和1.43%。此外，西兰花残体处理后厚壁菌门菌群上升为优势菌群，其相对丰度由0.51%升至0.79%，Saccharibacteria门菌群由1.01%降至0.34%（图5）。

在蕾期（T2），空白对照处理（CKT2）和西兰花残体处理（BRT2）分别获得25个和26个菌群，分别得到9个主要优势菌群，相对丰度占总体的96.01%和97.70%。变形菌门、放线菌门、酸杆菌门、芽单孢菌门、绿弯菌门、浮霉菌门、拟杆菌门、硝化螺旋菌门和蓝藻细菌门（Cyanobacteria）为共同的优势菌群，其相对丰度分别为24.13%和21.88%、21.12%和25.76%、15.79%和14.45%、11.87%和9.03%、11.80%和9.00%、3.49%和2.60%、2.68%和2.21%、3.78%和2.73%、1.36%和9.81%，其中西兰花残体处理后放线菌门和蓝藻细菌门增幅分别为21.97%和621.32%（图5）。

在花铃期（T3），空白对照处理（CKT3）和西兰花残体处理（BRT3）分别获得28个菌群，分别得到9个主要优势菌群，相对丰度占总体的96.53%和97.27%。变形菌门、放线菌门、酸杆菌门、芽单孢菌门、绿弯菌门、浮霉菌门、拟杆菌门和硝化螺旋菌门为两个处理共同的优势菌群，其相对丰度分别为26.44%和26.21%、25.76%和27.81%、12.62%和14.01%、12.97%和10.25%、9.63%和9.47%、2.24%和2.15%、2.76%和2.70%、3.08%和3.41%，其中西兰花残体处理后放线菌门、酸杆菌门和硝化螺旋菌门增幅分别为7.96%、11.01%和10.71%。此外，空白对照中的Saccharibacteria门成为优势菌群，而西兰花残体处理后蓝藻细菌门菌群成为优势菌群（图5）。

图5 西兰花残体在棉花不同生育时期对细菌群落结构的影响

在吐絮期（T4），空白对照处理（CKT4）和西兰花残体处理（BRT4）分别获得26个菌群，分别得到8个主要优势菌群，相对丰度分别占总体的97.22%和97.75%。变形菌门、放线菌门、酸杆菌门、芽单孢菌门、绿弯菌门、浮霉菌门、拟杆菌门和硝化螺旋菌门为共同的优势菌群，其相对丰度分别为32.23%和19.55%、20.18%和28.74%、15.31%和12.76%、9.62%和14.13%、10.51%和13.85%、4.11%和3.47%、1.83%和2.12%、3.42%和3.13%，其中西兰花残体处理后放线菌门、芽单孢菌门、绿弯菌门和拟杆菌门增幅分别为42.42%、46.88%、31.78%和15.85%（图5）。

细菌属水平研究结果表明（表2），在播种前期（T0），西兰花残体处理后（BRT0）土壤细菌组成占优势的前5类菌群为、鞘脂单胞菌属（）、芽孢杆菌属（）、和硝化螺菌属（），相对丰度分别为2.58%、1.57%、1.29%、1.16%和1.03%，与空白对照相比，鞘脂单胞菌属和菌群相对丰度降幅分别为27.98%和9.38%，而、芽孢杆菌属和硝化螺菌属菌群相对丰度分别为对照的9.92、3.91和1.27倍。

在苗期（T1），西兰花残体处理后（BRT1）土壤细菌组成占优势的前5类菌群为鞘脂单胞菌属、硝化螺菌属、气微菌属（）和芽单胞菌属（），与空白对照相比，鞘脂单胞菌属和硝化螺菌属菌群相对丰度差异不显著，而气微菌属和芽单胞菌属菌群相对丰度显著下降，降幅分别为26.85%、30.95%和18.52%。

在蕾期（T2），西兰花残体处理后（BRT2）土壤细菌组成占优势的前5类菌群为鞘脂单胞菌属、游动放线菌属（）、链霉菌属（）、和气微菌属，与空白对照相比，上述前5类菌群的相对丰度显著上升，其增幅分别为8.92%、159.70%、204.16%、76.47%和22.39%。

在花铃期（T3），西兰花残体处理后（BRT3）土壤细菌组成占优势的前5类菌群为鞘脂单胞菌属、类诺卡氏属（）、和芽单胞菌属，与空白对照相比，类诺卡氏属和菌群的相对丰度显著上升，其增幅分别为22.62%和23.17%，鞘脂单胞菌属、和芽单胞菌属的相对丰度显著下降，降幅分别为7.18%、11.57%和55.51%。

在吐絮期（T4），西兰花残体处理后（BRT4）土壤细菌组成占优势的前5类菌群为、鞘脂单胞菌属、芽单胞菌属、和玫瑰弯菌属（），与空白对照相比，芽单胞菌属和玫瑰弯菌属菌群相对丰度显著上升，其增幅分别为163.49%、96.49%和42.62%，鞘脂单胞菌属和菌群相对丰度显著下降，降幅分别为6.45%和41.22%。

3 讨论

3.1 作物残体还田对植物病害的影响

目前国内外很多研究表明，作物残体或秸秆还田对病害具有积极的防治效果[12-13,22-24]。Subbarao等[12]研究发现，西兰花残体还田后能够抑制花椰菜黄萎病的发生，降低了土壤中大丽轮枝菌微菌核的数量；Inderbitzin等[13-14]研究表明，西兰花残体还田后降低了茄子黄萎病的发生和病原菌数量；Davis等[15]研究表明，甜玉米秸秆还田后能够降低马铃薯黄萎病的发生，防治效果达到60%—70%，进一步研究表明秸秆还田后土壤中大丽轮枝菌的数量没有直接受到影响，但减少了地上部和根部病原菌的定殖。国内学者发现，芥菜混入土壤中对茄子黄萎病的发生具有良好的防治效果[22]；玉米秸秆还田能够降低番茄枯萎病的发生[23]；小麦秸秆还田能够降低赤霉病菌（）的存活率，减少病害的发生[24]。然而某些植物残体或秸秆还田能够增加病害的发生[25-26]，如花生秸秆还田对棉花种子萌发具有抑制作用，且秸秆用量与棉花黄萎病的发生呈正相关关系，其能够为病原菌的繁殖提供场所，导致土壤中大丽轮枝菌微菌核数量的增加[25]；水稻秸秆还田后能够增加稻瘟病、鞘腐病和纹枯病的发生[26]。然而有关西兰花残体还田对棉花黄萎病的影响鲜有报道。

棉花黄萎病是由大丽轮枝菌引起的土传病害，防治较为困难。本研究以‘邯686’为供试棉花品种，系统研究了西兰花残体还田对黄萎病的影响，结果表明西兰花残体还田降低了黄萎病的发病率和病情指数，且能够使发病高峰期推迟，在高峰期对黄萎病的防治效果为70.77%，在全生育期平均防效为57.21%。同时，比较了西兰花残体处理对大丽轮枝菌数量的影响，发现西兰花残体还田降低了大丽轮枝菌数量，土壤中病原菌的数量与其发病程度相关，进一步验证了Subbarao等[12-13,22]的观点。由此推测其防治黄萎病的原因之一是能够降低病原菌的数量。

表2 属水平棉花不同生育时期细菌菌群相对丰度（10组中至少1组相对丰度＞1%）

3.2 西兰花残体还田对土壤细菌群落多样性、组成和结构的影响

微生物群落结构多样性是维持土壤生态系统稳定和可持续的重要保障，是影响土壤微生态环境中生理活性最重要的因子，也是衡量土壤微生态环境质量和土壤生产力的重要指标[27-30]。许多学者研究表明，作物秸秆还田、有机肥料或绿肥等处理措施不仅对作物生长具有促生作用，而且能够改变根际微生物数量、促使活性增强，对土壤生态环境起到修复作用[31-33]。

Inderbitzin等[13]研究表明，西兰花残体还田后土壤微生物区系由“真菌型”向“细菌型”占主导趋势的方向转化；刘建国等[34]研究发现，棉花秸秆还田后不仅降低了土壤有害真菌的数量，而且增加了可培养微生物的数量，提高了土壤微生物的多样性指数、丰富度指数和均匀度指数，土壤微生物区系向着细菌占主导的趋势发展。本研究结果表明，土壤中添加西兰花残体后在棉花不同生育时期与空白对照相比，土壤细菌多样性增加，其中Chao1指数均显著上升，ACE指数在播种前期、花铃期和吐絮期显著上升，Shannon和Simpson指数差异不显著。在细菌组成和结构方面，在不同生育时期，西兰花残体处理和空白对照中土壤细菌占优势的门水平菌群相似，其优势菌群为变形菌门、放线菌门、酸杆菌门、芽单孢菌门、绿弯菌门、浮霉菌门、拟杆菌门和硝化螺旋菌门。主成分分析表明，不同生育时期土壤细菌群落结构存在不同程度的差异，除了在播种前期外，在苗期、蕾期、花铃期和吐絮期西兰花残体处理与空白对照相比门水平菌群的微生物群落结构差异显著。由表2可知，属水平细菌随着生育时期变化，其优势种群发生改变。如西兰花残体还田后土壤中鞘脂单胞菌属菌群的相对丰度呈现出先上升后下降的趋势。此外，在相同生育时期，西兰花残体还田后土壤细菌优势菌群发生改变，其相对丰度呈现出上升或下降的趋势，如在吐絮期，和芽单胞菌属菌群相对丰度显著上升，其他菌群如鞘脂单胞菌属和相对丰度显著下降。

通过西兰花残体还田对棉花黄萎病防治效果的研究，结合其对不同生育时期土壤细菌群落的影响，得出西兰花残体还田后能够通过改变土壤细菌菌群的相对丰度，且不同生育时期主要优势菌群的相对丰度存在差别。

3.3 西兰花残体还田对土壤细菌有益种类的影响

链霉菌属、芽孢杆菌属、溶杆菌属菌群大多属于有益微生物，是重要的拮抗微生物种类，对植物病害的生物防治方面具有重要的作用[35-37]。本研究发现，在不同生育时期，西兰花残体还田后与空白对照相比，链霉菌属菌群的相对丰度呈显著上升趋势，在播种前期、苗期、蕾期、花铃期和吐絮期增幅分别为186.21%、52.17%、141.67%、41.03%和103.70%；对于芽孢杆菌属菌群而言，除了在吐絮期西兰花残体还田和空白对照之间差异不显著外，在其他时期菌群的相对丰度均呈现显著上升趋势，在播种前期、苗期、蕾期和花铃期增幅分别为290.91%、169.57%、32.50%、110.34%。此外，溶杆菌属菌群的相对丰度在蕾期、花铃期和吐絮期呈显著上升趋势。由此推测，西兰花残体还田后在一定程度上增加了有益微生物的相对丰度，进而对黄萎病起到防治作用。

已有研究表明，不同种植年限土壤微生物群落结构多样性具有明显差异，并且土壤微生物群落结构受土壤养分、pH、温度、水分和通气等环境条件的影响，因此能够通过改变土壤理化性质的管理措施，在一定程度上改变土壤微生物群落结构[38-39]，此外，微生物群落结构还与土壤类型、作物品种及类型等因素存在相关性[40-41]。本研究由于时间的局限性，仅从作物不同生育时期方面对土壤微生物群落结构进行了研究，西兰花残体处理与空白对照相比，在苗期、蕾期、花铃期和吐絮期门水平菌群的微生物群落结构差异显著。

此外，土壤微生物不仅是土壤有机质和养分转化与循环的动力，本身也是土壤养分的储存库，对土壤养分的转化与供应起着重要作用，是土壤肥力水平的有效指标[39]。研究表明，植物秸秆或残体还田能够改变土壤养分，与秸秆或残体类型、还田量、还田时间以及还田土壤类型有关[42]。开展西兰花残体还田对土壤理化性质的变化及其相关性研究，可为进一步衡量和评价西兰花残体还田对棉田土壤肥力和土壤质量的影响提供依据。

4 结论

西兰花残体还田是一种操作简便且节约成本的棉花黄萎病防治方法，既有效解决了西兰花残体资源浪费的问题，又杜绝了西兰花残体随意处置所造成的环境问题。通过利用西兰花残体还田，可改善土壤微生物群落结构，降低土壤中病原菌数量，增加有益微生物的相对丰度，从而达到对黄萎病菌的抑制作用，是一种有效的控制棉花黄萎病的措施。

[1] 王娜, 陶伟, 陈双林, 闫淑珍. 植物内生细菌在棉花体内的定殖动态及对棉花黄萎病的生物防治效果. 植物保护学报, 2016, 43(2): 207-214.

Wang N, Tao W, Chen S L, Yan S Z. Colonization trends of endophytic bacteria in cotton and their biological control effect on cotton Verticillium wilt., 2016, 43(2): 207-214. (in Chinese)

[2] Fradin E F, Thomma B P. Physiology and molecular aspects of Verticillium wilt diseases caused byand, 2006, 7(2): 71-86.

[3] Erdogan O, Benlioglu K. Biological control of Verticillium wilt on cotton by the use of fluorescentspp. under field conditions., 2010, 53(1): 39-45.

[4] Kurt S, Dervis S, Sahinler S. Sensitivity ofto prochloraz and prochloraz-manganese complex and control of Verticillium wilt of cotton in the field., 2003, 22(1): 51-55.

[5] Zheng Y, Xue Q Y, Xu L L, Xu Q, Lu S, Gu C, Guo J H. A screening strategy of fungal biocontrol agents towards Verticillium wilt of cotton., 2011, 56(3): 209-216.

[6] 周兆华, 张天真, 潘家驹, 陈兆夏, 蒋伟民. 大丽轮枝菌在棉花品种上的致病力分化研究. 中国农业科学, 2000, 33(2): 51-57.

Zhou Z H, Zhang T Z, Pan J J, Chen Z X, Jiang W M. Pathogenicity differentiation ofon cotton varieties., 2000, 33(2): 51-57. (in Chinese)

[7] 朱荷琴, 简桂良, 宋晓轩. 棉田黄萎病菌致病型群落结构研究. 棉花学报, 2004, 16(3): 147-151.

Zhu H Q, Jian G L, Song X X. Community structure of pathogenic type ofKleb. in cotton field., 2004, 16(3): 147-151. (in Chinese)

[8] 简桂良, 卢美光, 朱荷琴. 棉田黄萎病菌致病型结构初步研究. 植物保护学报, 2005, 32(1): 109-110.

Jian G L, Lu M G, Zhu H Q. A preliminary study on the pathogenic structure ofin cotton field., 2005, 32(1): 109-110. (in Chinese)

[9] 巩彪, 张丽丽, 隋申利, 王秀峰, 魏珉, 史庆华, 杨凤娟, 李岩. 大蒜秸秆对番茄根结线虫病及根际微生态的影响. 中国农业科学, 2016, 49(5): 933-941.

Gong B, Zhang L L, Sui S L, Wang X F, Wei M, Shi Q H, Yang F J, Li Y. Effects of garlic straw application on controlling tomato root-knot nematode disease and rhizospheric microecology., 2016, 49(5): 933-941. (in Chinese)

[10] 连宾, 王进军, 陆玲. 植物与微生物的化感作用. 南京师大学报(自然科学版), 2007, 30(1): 88-94.

Lian B, Wang J J, Lu L. Research progress of allelopathy on microorganism to plant., 2007, 30(1): 88-94. (in Chinese)

[11] Shetty K G, Subbarao K V, Huisman O C, Hubbard J C. Mechanism of broccoli-mediated Verticillium wilt reduction in cauliflower., 2000, 90(3): 305-310.

[12] Subbarao K V, Hubbard J C, Koike S T. Evaluation of broccoli residue incorporation into field soil for Verticillium wilt control in cauliflower., 1999, 83(2): 124-129.

[13] Inderbitzin P, Ward J, Barbella A, Solares N, Izyumin D, Burman P, Chellemi D O, Subbarao K V. Soil microbiomes associated with Verticillium wilt-suppressive broccoli and chitin amendments are enriched with potential biocontrol agents., 2018, 108(1): 31-43.

[14] Ikeda K, Banno S, Furusawa A, Shibata S, Nakaho K, Fujimura M. Crop rotation with broccoli suppresses verticillium wilt of eggplant., 2015, 81(1): 77-82.

[15] Davis J R, Huisman O C, Everson D O, NoltE P, Sorenson L H, Schneider A T. The suppression of Verticillium wilt of potato using corn as a green manure crop., 2010, 87(2): 195-208.

[16] Larkin R P, Honeycutt C W, Olanya M O. Management of Verticillium wilt of potato with disease-suppressive green manures and as affected by previous cropping history., 2011, 95(5): 568-576.

[17] Njoroge S M C, Kabir Z, Martin F N, Koike S T, Subbarao K V. Comparison of crop rotation for Verticillium wilt management and effect onspecies in conventional and organic strawberry production., 2009, 93(5): 519-527.

[18] 赵卫松, 李社增, 鹿秀云, 郭庆港, 王亚娇, 王培培, 张晓云, 董丽红, 苏振贺, 马平. 西兰花植株残体还田对棉花黄萎病的防治效果及其安全性评价. 中国生物防治学报, 2019, 35(3): 449-455.

Zhao W S, Li S Z, Lu X Y, Guo Q G, Wang Y J, Wang P P, Zhang X Y, Dong L H, Su Z H, Ma P. Safety evaluation of broccoli residues returning on cotton growth and its control effect on Verticillium wilt., 2019, 35(3): 449-455. (in Chinese)

[19] 李社增, 马平, 刘杏忠, Huang H C, 陈新华. 利用拮抗细菌防治棉花黄萎病. 华中农业大学学报, 2001, 20(5): 422-425.

Li S Z, Ma P, Liu X Z, Huang H C, Chen X H. Biological control of cotton Verticillium wilt by antagonistic bacteria., 2001, 20(5): 422-425. (in Chinese)

[20] Lu R, Li Y P, Li W X, Xie Z S, Fan C L, Liu P X, Deng S X. Bacterial community structure in atmospheric particulate matters of different sizes during the haze days in xi’an, China., 2018, 637/638: 244-252.

[21] 肖蕊, 余真真, ElsharawyA A, 魏锋, 杨家荣. 土壤中棉花黄萎病菌SYBR Green Ⅰ荧光RT-PCR定量检测技术研究. 菌物学报, 2011, 30(4): 598-603.

Xiao R, Yu Z Z, ElsharawyA A, Wei F, Yang J R. SYBR Green I real time RT-PCR assay for quantitatively detectingthe occurrence ofof cotton in naturally infestedsoil., 2011, 30(4): 598-603.(in Chinese)

[22] 李淑敏, 郑成彧, 张润芝, 杨自超, 曲红云, 刘彤彤, 袁睿, 姚小桐, 王雪蓉, 许宁, 张春怡. 生物熏蒸对大棚连作茄子产量和黄萎病发病率影响. 东北农业大学学报, 2017, 48(5): 35-41.

Li S M, Zheng C Y, Zhang R Z, Yang Z C, Qu H Y, Liu T T, Yuan R, Yao X T, Wang X R, Xu N, Zhang C Y. Effect of biofumigation on yield and Verticillium wilt incidence of continuous eggplant in greenhouse., 2017, 48(5): 35-41. (in Chinese)

[23] 乔俊卿, 刘邮洲, 余翔, 梁雪杰, 陈志谊, 刘永锋, 张英. 集成生物防治和秸秆还田技术对设施番茄增产及土传病害防控效果研究. 中国生物防治学报, 2013, 29(4): 547-554.

Qiao J Q, Liu Y Z, Yu X, Liang X J, Chen Z Y, Liu Y F, Zhang Y. Evaluation of yield increasing and control efficiency of tomato soil-borne diseases under the integrated application of straw returning and biocontrol agent., 2013, 29(4): 547-554. (in Chinese)

[24] 张政文, 孙秀娟, 胡乃娟, 黄卉, 朱利群, 卞新民. 麦秸掩埋还田对赤霉菌(sehw.)存活率和水稻带菌的影响. 植物病理学报, 2014, 44(4): 443-448.

Zhang Z W, Sun X J, Hu N J, Huang H, Zhu L Q, Bian X M. Effect of straw ditch-buried return depth onsurvival rate and carrier rate of rice., 2014, 44(4): 443-448. (in Chinese)

[25] Chawla S, Woodward J E, Wheeler T A. Influence ofinfested peanut residue on wilt development in subsequent cotton., 2012, 2012: Article ID 212075.

[26] 李洪林, 刘凤艳, 龚振平, 靳学慧, 黄世文. 稻秆还田对水稻主要病害发生的影响. 作物研究, 2012, 26(1): 7-10.

Li H L, Liu F Y, Gong Z P, Jin X H, Huang S W. Effect of straw return back to paddy field on occurrence of rice major diseases., 2012, 26(1): 7-10. (in Chinese)

[27] Enwall K, Nyberg K, Bertilsson S, Cederlund H, Stenström J, HalliN S. Long-term impact of fertilization on activity and composition of bacterial communities and metabolic guilds in agricultural soil., 2007, 39(1): 106-115.

[28] Schutter M E, Sandeno J M, Dick R P. Seasonal, soil type, and alternative management influences on microbial communities of vegetable cropping systems., 2001, 34(6): 397-410.

[29] Gómez-Sagasti M T, Alkorta I, Becerril J M, Epelde L, Anza M, Garbisu C. Microbial monitoring of the recovery of soil quality during heavy metal phytoremediation., 2012, 223(6): 3249-3262.

[30] 吴林坤, 林向民, 林文雄. 根系分泌物介导下植物-土壤-微生物互作关系研究进展与展望. 植物生态学报, 2014, 38(3): 298-310.

Wu L K, Lin X M, Lin W X. Advances and perspective in research on plant-soil-microbe interactions mediated by root exudates., 2014, 38(3): 298-310. (in Chinese)

[31] Srivastava R, Khalid A, Singh U S, Sharma A K. Evaluation of arbuscular mycorrhizal fungus, fluorescentandformulation againstf. sp.for the management of tomato wilt., 2010, 53(1): 24-31.

[32] 袁英英, 李敏清, 胡伟, 张静, 赵兰凤, 李华兴. 生物有机肥对番茄青枯病的防效及对土壤微生物的影响. 农业环境科学学报, 2011, 30(7): 1344-1350.

Yuan Y Y, Li M Q, Hu W, Zhang J, Zhao L F, Li H X. Effect of biological organic fertilizer on tomato bacterial wilt and soil microorganism., 2011, 30(7): 1344-1350. (in Chinese)

[33] 蒋岁寒, 刘艳霞, 孟琳, 朱春波, 李想, 沈标, 石俊雄, 杨兴明. 生物有机肥对烟草青枯病的田间防效及根际土壤微生物的影响. 南京农业大学学报, 2016, 39(5): 784-790.

Jiang S H, Liu Y X, Meng L, Zhu C B, Li X, Shen B, Shi J X, Yang X M. Effects of biological organic fertilizer on field control efficiency of tobaccowilt and soil microbial in rhizosphere soil., 2016, 39(5): 784-790. (in Chinese)

[34] 刘建国, 卞新民, 李彦斌, 张伟, 李崧. 长期连作和秸秆还田对棉田土壤微生物活性的影响. 应用生态学报, 2008, 19(5):1027-1032.

Liu J G, Bian X M, Li Y B, Zhang W, Li S. Effects of long-term continuous cropping of cotton and returning cotton stalk into field on soil biological activities., 2008, 19(5):1027-1032. (in Chinese)

[35] Alam M, Dharni S, Abdul-Khaliq, Srivastava S K, Samad A, Gupta M K. A promising strain ofsp. with agricultural traits for growth promotion and disease management., 2012, 50(8):559-568.

[36] Goudjal Y, Zamoum M, Sabaou N, Mathieu F, Zitouni A. Potential of endophyticspp. for biocontrol of Fusarium root rot disease and growth promotion of tomato seedlings., 2016, 26(12):1691-1705.

[37] 李社增, 鹿秀云, 马平, 高胜国, 刘杏忠, 刘干. 防治棉花黄萎病的生防细菌NCD-2的田间效果评价及其鉴定. 植物病理学报, 2005, 35(5):451-455.

Li SZ, LU X Y, Ma P, Gao S G, Liu X Z, Liu G. Evaluation of biocontrol potential of a bacterial strain NCD-2 against cotton Verticillium wilt in field trials., 2005, 35(5): 451-455. (in Chinese)

[38] 王光华, 金剑, 徐美娜, 刘晓冰. 植物、土壤及土壤管理对土壤微生物群落结构的影响. 生态学杂志, 2006, 25(5): 550-556.

Wang G H, Jin J, Xu M N, Liu X B. Effects of plant, soil and soil management on soil microbial community diversity., 2006, 25(5): 550-556. (in Chinese)

[39] 宋蒙亚, 李忠佩, 吴萌, 刘明, 江春玉. 不同种植年限设施菜地土壤微生物量和群落结构的差异. 中国农业科学, 2015, 48(18): 3635-3644.

Song M Y, Li Z P, Wu M, Liu M, Jiang C Y. Changes in soil microbial biomass and community structure with cultivation chronosequence of greenhouse vegetables., 2015, 48(18): 3635-3644.(in Chinese)

[40] 武爱莲, 丁玉川, 焦晓燕, 王劲松, 董二伟, 郭珺, 王浩. 玉米秸秆生物炭对褐土微生物功能多样性及细菌群落的影响. 中国生态农业学报, 2016, 24(6): 736-743.

Wu A L, Ding Y C, Jiao X Y, Wang J S, Dong E W, Guo J,WangH. Effects of corn-stalk biochar on soil microbial functionaldiversity and bacterial community in cinnamon soils., 2016, 24(6): 736-743. (in Chinese)

[41] 侯建伟,邢存芳,卢志宏,陈芬,余高. 不同秸秆生物炭对贵州黄壤细菌群落的影响. 中国农业科学, 2018, 51(23): 4485-4495.

Hou J W, Xing C F, Lu Z H, Chen F, Yu G. Effects of the different crop straw biochars on soil bacterial community of yellow soil in Guizhou., 2018, 51(23): 4485-4495.(in Chinese)

[42] Senesi N, Plaza C, Brunetti G, Polo A. A comparative survey of recentresults on humic-like fractions in organic amendments and effects onnative soil humic substances., 2007, 39(6):1244-1262.

Control Efficacy of Broccoli Residues on Cotton Verticillium Wilt and Its Effect on Soil Bacterial Community at Different Growth Stages

ZHAO Weisong, GUO Qinggang, LI Shezeng, WANG Yajiao, LU Xiuyun, WANG Peipei, SU Zhenhe, ZHANG Xiaoyun, MA Ping

(Plant Protection Institute of Hebei Academy of Agricultural and Forestry Sciences/IPM Centre of Hebei Province/Key Laboratory of IPM on Crops in Northern Region of North China, Ministry of Agriculture and Rural Affairs, Baoding 071000, Hebei)

【Objective】The objective of this study is to research the effects of broccoli residues (BR) on cotton verticillium wilt and soil bacterial community structure at different growth stages, and to provide new ways and ideas for green ecological control of cotton verticillium wilt and reduction of chemical pesticides.【Method】In the field plot experiment, the soil planted with cotton was treated with broccoli residues and no broccoli residues (CK). the incidence dynamics of cotton verticillium wilt in different treatments were monitored at different growth stages. The number of DNA copies ofwas determined by real-time quantitative PCR, and the structure of soil bacterial communities was measured by high-throughput sequencing (Illumina MiSeq). The effects of broccoli residues returning to soil on the quantity of pathogen and the community of soil bacteria were analyzed. Principal component analysis (PCA) was used to explore the characteristics and rules of structural changes of soil bacterial community at the level of phylum and genus in different growth stages.【Result】The incidence and disease index of cotton verticillium wilt showed a downward trend after the return of broccoli residues to soil, respectively. Moreover, the peak period of verticillium wilt was delayed, and 70.77% control efficacy was achieved at the peak period. the development curve of disease index and time course of disease in the whole growth period was established, and the average control efficacy was 57.21%. Compared to the blank control, the number of DNA copies ofwas decreased by 10.96%, 11.11%, 25.95% and 11.25% at the 4 stages of cotton growing and non-growing season, respectively, after the return of broccoli residues to soil. Illumina MiSeq analysis showed that broccoli residues treatment significantly increased soil bacterial diversity. The richness index Chao1 was increased significantly at all growth stages, and the ACE index was significantly increased at the pre-sowing stage, flowering and boll stage and boll-opening stage.Compared to the blank control, the relative abundance of Actinobacteria was increased significantly at different growth stages after the return of broccoli residues to soil. Meanwhile, the results showed that there were differences in the composition of dominant microorganisms at different growth stages after the return of broccoli residues to soil. For example, Firmicutes increased to the dominant microflora at the pre-sowing stage and seedling stage, and Cyanobacteria became the dominant microflora at the bud stage and flowering and boll stage. Principal component analysis showed that broccoli residues to soil changed the structure of bacterial community at seedling stage, bud stage, flowering and boll stage and boll-opening stage. Further analysis showed that the relative abundance ofandwas significantly increased after the return of broccoli residues to soil.【Conclusion】Broccoli residues returning to soil can effectively inhibit the occurrence of cotton verticillium wilt, reduce the number of DNA copies ofin soil, change the structure of soil bacterial community, and increase the relative abundance of beneficial microorganisms, which is an effective green measure to control cotton verticillium wilt.

broccoli residues;; cotton verticillium wilt; bacterial community structure; Illumina MiSeq; growth stage

2019-06-03；

2019-06-28

国家重点研发计划（2017YFD0200601）、国家公益性行业（农业）科研专项（201503109）、国家棉花产业技术体系（CARS-15-17）

赵卫松，Tel：0312-5927076；E-mail：zhaoweisong1985@163.com。通信作者马平，Tel：0312-5915678；E-mail：pingma88@126.com

（责任编辑 岳梅）