土壤是植物生存和发展的基础。微生物作为土壤中最活跃的组成部分，由细菌、古菌、真菌等组成，数量庞大且种类多样，是土壤有机质和养分循环转化的主要推动力^[1]。不同的种植制度会改变土壤微生物组成，进而影响作物的生长发育。作物连作使土壤病原微生物富集，有益微生物减少，导致土壤微生态失衡，土传病害发生严重^[2-6]。轮作能够促进有益微生物活性，抑制有害微生物，在一定程度上减少连作障碍的发生，从而提高作物的产量和品质^[7-9]。

细菌是土壤微生物中数量最多、分布最广的组成部分^[10]，其群落结构与组成不仅能够反映土壤性状的变化趋势，而且能直接影响土壤生态功能的发挥^[11-12]。吕雅悠等^[13]研究发现根际细菌A21-4可以通过促进辣椒土壤酶的活性，从而提高土壤速效氮磷钾的含量。王娜^[14]的研究表明枯草芽孢杆菌DP10和链霉菌DK56对棉花土壤中大丽轮枝菌的抑制性最强，生防效果显著。目前，关于土壤微生物的研究已从传统的平板培养转变为分子生物学方法，尤其是高通量测序技术，其具有数据通量高、微生物检测结果多样丰富等优势而被迅速应用^[15]。近年来，该技术被广泛应用于植物土壤微生物的研究，但关于轮作下棉花根际土壤细菌变化的研究还鲜有报道。本研究利用高通量测序技术，通过分析花生、谷子轮作后棉花根际土壤细菌群落结构的演替规律，对选择合理种植措施和改善棉田土壤生态功能具有重要意义。

1 材料与方法 1.1 材料

棉花品种为新陆早61，谷子品种为豫18，花生品种为中花413。

1.2 方法 1.2.1 试验设计

试验地点位于新疆昌吉回族自治州的中国农业科学院西部农业研究中心。于2017年和2018年开展棉田轮作试验，设置棉花-棉花连作、花生-棉花轮作、谷子-棉花轮作、谷子-谷子连作、花生-花生连作共5种处理方式，每个处理3次重复，采用随机区组排列。各作物覆膜种植，一膜三行，水肥以膜下滴灌的方式施入。

1.2.2 棉花产量测定

2018年棉花吐絮后实收测产。产量数据通过SPSS(24.0)进行统计分析。

1.2.3 土壤样品采集

2018年9月10日，在各试验小区采用S形随机取样的方式，每小区取8个点混样。利用抖根法取20 cm深的粘附于作物根系的土壤鲜样，标记并立即置于冰上保存，带回实验室进行后续处理。

1.2.4 测序样品制备

采用CTAB方法提取土壤样品总DNA，利用细菌检测通用引物16S V4区引物(515F/806R)进行PCR扩增，产物经琼脂糖凝胶电泳检测并回收后进行浓度测定和等量混样，充分混匀后的产物再经2%琼脂糖胶凝胶电泳纯化，委托北京诺禾致源股份有限公司进行文库构建及上机测序。

1.2.5 测序数据分析

按照97%的一致性(Identity)，利用Uprase软件将所有样品的Clean Reads聚类成OTUs(Operational Taxonomic Units)，通过软件TBtools绘制韦恩图。使用Qiime软件计算Shannon指数、Simpson指数、Chao1指数、ACE指数等微生物α多样性指数；计算各样品间的Unifrac距离，构建UPGMA(Unweighted Pair-group Method with Arithmetic Mean)聚类树。通过稀释曲线绘制和主成分分析(PCA，Principal Component Analysis)对测序数据进行评价。

2 结果 2.1 不同轮作处理对棉花产量的影响

与棉花连作相比，花生-棉花轮作使棉花单铃重提高了17.2%，籽棉产量增加了32.19%；谷子-棉花轮作使棉花单铃重提高了14.5%，籽棉产量增加了18.13%(表 1)，表明轮作显著提高了棉花铃重和籽棉产量。

2.2 根际土壤细菌基因组测序评价

试验中所有样品经过质控平均得到80，167条有效数据（表 2），质控有效率达93.85%。随着测序量的不断增加，各处理稀释曲线趋向平坦（图 1），表明样品文库测序数据量合理，测序深度已基本覆盖到样品中的所有细菌。

2.3 不同种植方式下根际土壤细菌OTUs分布

以97%作为一致性标准，将所有序列聚类成为OTUs，5个处理共产生6 174个OTUs。各处理OTUs数量由表 3可知，轮作组（处理B和处理C）OUTs数量低于连作组（处理A、处理D、处理E），表明轮作处理降低了根际土壤细菌物种数量。

OTUs分布如图 2所示，所有处理共有的OTUs为2 760个，占全部OTUs的44.7%。3个连作处理的特有OTUs分布中，处理A与处理D之间，A特有OTUs为820个，D特有OTUs为982个；处理A与处理E之间，A特有OTUs为873个，E特有OTUs为1 086个。花生连作和谷子连作中特有OUTs数量明显多于棉花连作。2个轮作处理与棉花连作的特有OTUs分布中，处理A与处理B之间，A特有OTUs为821个，B特有OTUs为775个；处理A与处理C之间，A特有OTUs为805个，C特有OTUs为722个。同棉花连作相比，轮作处理降低了棉花根际土壤细菌中特有的OTUs数量。

2.4 不同种植方式下根际土壤细菌α多样性分析

为进一步分析不同种植方式下根际土壤细菌的多样性和丰富度，采用α多样性分析方法，分析结果见表 4，反映物种多样性的Shannon指数、Simpson指数显示，花生-棉花轮作的Shannon指数低于棉花连作，Simpson指数无差异；谷子-棉花轮作的Shannon指数高于棉花连作，Simpson指数无差异。反映物种丰富度的Chao1指数、ACE指数显示，2个轮作处理分别与棉花连作相比，其各项指数均降低；花生-棉花轮作与花生连作相比，其各项指数均显著降低；谷子-棉花轮作与谷子连作相比，其各项指数均显著降低。同连作处理相比，轮作处理降低了根际土壤细菌丰富度，物种数量减少。

2.5 不同种植方式下根际土壤细菌群落组成及相对丰度分析

在门的水平上(图 3)，5种种植方式共发现10个主要菌门(> 1%)，分别是变形菌门(Proteobacteria，29.69%-33.38%)、放线菌门(Actinobacteria，15.09%-19.39%)、拟杆菌门(Bacteroidetes，7.63%- 11.31%)、芽单胞菌门(Gemmatimonadetes，7.70%- 10.23%)、酸杆菌门(Acidobacteria，6.70%-9.03%)、厚壁菌门(Firmicutes，1.84%-4.55%)、绿弯菌门(Choroflexi，5.54%-6.96%)、浮霉菌门(Planctomycetes，2.83%-4.75%)、未辨识菌门(unidentified- Bacteria，0.78%-2.94%)、疣微菌门(Verrucomicrobia，1.61%-2.43%)。这些优势菌门中的变形菌门在谷子连作中占比33.38%，花生连作中占比32.15%，棉花连作中占比29.69%，同棉花连作相比，轮作处理提高了变形菌门的相对丰度，花生-棉花轮作提高了1.58%，谷子-棉花轮作提高了1.29%；浮霉菌门在棉花连作中相对丰度为2.83%，是5个处理中丰度最低的，但在谷子连作中的相对丰度高达3.97%，花生连作中的相对丰度高达4.74%，同棉花连作相比，花生-棉花和谷子-棉花两个轮作处理的浮霉菌门相对丰度分别提高了18.2%和13.7%。拟杆菌门在3个连作处理(棉花连作、谷子连作、花生连作)中的相对丰度分别9.19%、9.54%、7.63%，在2个轮作处理(花生-棉花轮作和谷子-棉花轮作)中的相对丰度分别为11.31%和11.07%。与连作处理相比，轮作处理提高了拟杆菌门的相对丰度。

在属的水平上(图 4)，所有处理相对丰度排名前十的菌属分别有鞘氨醇单胞菌属(Sphingomonas，2.98%-3.95%)、未知细菌属(unidentified_Bacteria，0.64%-2.67%)、未知酸杆菌属(unidentified_ Acidobacteria，1.79%-2.82%)、Pontibacter(0.13%- 1.83%)、溶杆菌属(Lysobacter，0.62%-1.93%)、拟杆菌属(Bacteroides，0.01%-0.70%)、莱茵海默氏菌属(Rheinheimera，0.007%-0.88%)、Haliangium (1.02%-1.54%)、类诺卡氏菌属(Nocardioides，1.05%-1.81%)、小梨形菌属(Pirellula，0.59%- 1.32%)。其中鞘氨醇单胞菌属(Sphingomonas)在5个处理中的相对丰度最大(≥ 3%)，该菌属在花生连作和谷子连作中相对丰度最高，分别为3.55%、3.96%，在棉花连作中丰度最低为2.98%，通过轮作处理后，该菌属在棉花根际土壤细菌中的相对丰度均提高，在花生-棉花轮作中的丰度提高到3.93%，在谷子-棉花轮作中的丰度提高到3.30%。该结果表明轮作处理提高了棉花根际土壤中鞘氨醇单胞菌属的相对丰度。

2.6 土壤细菌群落组成PCA分析及丰度UPGMA聚类分析

基于OTU水平的PCA分析结果如图 5所示，主成分1（PCA1）和主成分2（PCA2）对样品的贡献值分别为42.69%和21.41%。处理A（棉花-棉花）与处理B（花生-棉花）、处理C（谷子-棉花）两个处理的距离较近，而与处理D（谷子-谷子）、处理E（花生-花生）两个处理的距离较远。

UPGMA聚类分析如图 6所示，棉花-棉花、花生-棉花和谷子-棉花3个处理因其根际细菌群落结构相似度高被聚类在一起，谷子-谷子、花生-花生2个处理因菌群结构相似被聚类为一个分支，这与PCA分析的结果一致。

3 讨论

连作导致土壤病原菌富集，土壤质量恶化，从而使得作物产量降低。研究表明，轮作能协调作物与土壤之间的关系，改善根际微生态环境，有助于提高作物产量，是防治作物发生连作障碍的有效途径之一^{[3, 16-18]}。本研究中对比棉花连作、轮作处理中棉花的单铃重和籽棉产量均显著提高，表明轮作处理对棉花产量有积极作用。

植物根际为微生物提供能源物质，有研究证明轮作制度通过改变作物的种类和生物量的方式增加了土壤微生物的食物来源，进而导致土壤微生物多样性的改变^[19]。本研究通过高通量测序技术以97%的一致性将序列聚类出共6 174个OTUs，根据各处理OTUs分布绘制韦恩图，用来分析不同样本（组）之间共有、特有的OTUs，其中棉花连作、谷子连作和花生连作3个连作处理中共有OTUs为3 224个，棉花连作中特有OTUs为463个；棉花参与的处理（棉花连作、花生-棉花轮作和谷子-棉花轮作）中共有OTUs为3 236个，棉花连作中特有OTUs为408个，这表明轮作处理后的棉花根际土壤中共有细菌OUTs增多，特有OTUs减少。α多样性指数用来解释物种的多样性和丰富度，细菌α多样性指数显示，代表根际土壤细菌多样性的Shannon指数和Simpson指数无显著差异，但代表样品群落丰富度的ACE指数和Chao1指数显示，轮作处理的数值低于连作处理。通过对本试验根际土壤细菌OTUs分布和细菌α多样性指数的分析可知，轮作处理使棉花根际土壤细菌物种数量减少，丰富度明显降低。有研究表明连作可以导致病原微生物富集，降低土壤微生物活性及数量^[20-22]，所以本试验存在轮作措施促进了棉花根际土壤中某些有益细菌的活性，例如拟杆菌门、变形菌门、浮霉菌门和鞘氨醇单胞菌属的富集，从而抑制了病原菌存活的现象。

土壤中庞大的微生物数量意味着其群落结构的复杂性，群落结构的变化趋势以群落组分变化为基础^[23]，不同的作物类型，对应其土壤细菌群落的组成成分不同，这就意味着生长作物的变化为某些特异性细菌提供了生存环境。为此，菌群针对作物的改变重新进行了选择性适应。王悦等^[24]通过研究不同种植模式下丹参根际土壤微生物群落结构发现轮作模式下芽单胞菌门相对丰度低于连作和套作，轮作和套作时浮霉菌门和拟杆菌门的相对丰度则高于连作。试验通过分析不同种植方式下的根际土壤细菌群落组成及相对丰度发现，门水平上的变形菌门、浮霉菌门和属水平上的鞘氨醇单胞菌属，其相对丰度在花生连作和谷子连作中高，而在棉花连作中低，通过轮作处理后，其菌门（属）在棉花根际土壤细菌中的相对丰度均明显提高；另有拟杆菌门在轮作中的相对丰度均高于连作。这表明了轮作处理中前茬作物的存在给予了棉花根际土壤中这些菌门（属）适宜的生存环境，促进了它们的活性。土壤细菌在植物生长过程中发挥着重要作用，孙建光等^[25]从玉米根际土壤中分离出鞘氨醇单胞菌属菌株GD542，该菌株具有较高的固氮酶活性，可以提高土壤中的全氮含量；王秀红等^[26]研究发现变形菌门、浮霉菌门、拟杆菌门等是好氧堆肥过程中共享的优势菌门，它们的存在提高了植物对土壤碳源的利用效率。通过前人的研究成果可以推断，试验中受轮作处理影响的细菌菌门（属）可以为棉花的生长带来积极地影响。

主成分分析^[27]是一种基于欧式距离的应用 方差分解，对多维数据进行降维，从而提取出数据中最主要的元素和结构的方法。各样品在PCA图中的距离越接近，则它们的群落组成越相似。UPGMA^[28]是为了研究不同样本间相似性的聚类分析方法。基于OTU水平的PCA分析和UPGMA聚类分析发现，花生-棉花轮作处理、谷子-棉花轮作处理和棉花连作处理的根际土壤细菌组成相似性更强，这表明虽然棉花不同的种植方式会影响其菌群相对丰度的高低，但其群落结构组成总体上不会发生显著变化。这与Li等^[29]发现根系分泌物虽可改变细菌群落中某些主要种群的丰富度，但是不足以改变细菌的群落结构的研究结论一致。

本研究使用的高通量测序技术仅为二代测序技术，不能精确到种的水平，而且试验只检测了棉花根际土壤中细菌的变化。未来应使用精细到种水平上的三代测序技术，来监测棉田根际土壤中细菌、真菌的群落结构变化，这对进一步探究棉花根际土壤微生态功能具有重要意义。

4 结论

本研究发现对棉田进行花生-棉花轮作处理可以使棉花的增产效果更显著。通过高通量测序技术，检测各处理的根际土壤细菌发现，轮作处理中的前茬作物（花生、谷子）影响了棉田根际土壤菌群的丰度，它们的存在提高了棉田中拟杆菌门、变形菌门、浮霉菌门、鞘氨醇单胞菌属的相对丰度，这些菌门（属）的含量与棉花的生长发育密切相关。