香蕉枯萎病是由真菌尖孢镰刀菌（Fusarium oxysporum）侵染引起的^[1]，又称为巴拿马病或黄叶病，于1904年在美国夏威夷首次被发现^[2]，是病虫害中对香蕉危害最为严重的一种。尖孢镰刀菌会破坏植株的维管束结构，使其木质部导管堵塞，从而导致植株假茎开裂、叶片发黄，随后植株枯萎死亡^[3]。

根际微生物组是抵御香蕉枯萎病病原菌（致病性尖孢镰刀菌）入侵并保护植物健康的关键屏障，其对病原菌的拮抗作用使得植物能够免受其侵染^[4]。然而，根际微生物群落在土传病害防治上的研究多集中于细菌和真菌^[5-6]，原生动物作为根际微生物组中另一个重要的微生物群体，受到的关注较少^[7-8]。原生动物在根际具有多种功能：部分原生动物可以直接对细菌及真菌进行捕食，扮演着根际土壤食物网中消费者的角色；部分原生动物寄生于植物；也有部分原生动物参与了分解有机质的过程，在碳氮循环和养分转化过程中发挥着重要作用等等^[9-11]。

研究表明，原生动物对于土壤环境变化的响应比细菌与真菌更强烈，患病植株中，根际原生动物的群落结构对病原菌种群数量的解释力度高于细菌与真菌^[12]，原生动物对植株的健康起到了预测作用^[13-14]。Foissner^[15]证实当土壤原生动物存在时，即使加入病原菌也不会影响作物的正常生长，而不接种原生动物的对照组则在几天之内发病枯萎。最新的研究结果表明吞噬型原生动物可能能够通过直接捕食作用来抑制病原青枯菌的数量，推动了对番茄青枯病的防治进展^[12]。而原生动物在香蕉土传枯萎病的防治上鲜有报道。因此，本研究通过采用近年来在探索微生物群落中使用较为广泛的高通量测序方法^[16]来探究连作蕉园中植株发病前后根际土壤原生动物群落特征以及其与病原菌的关系，以探明根际原生动物在香蕉枯萎病防控中的潜力。

试验于2018年10月至2019年8月在海南省澄迈县桥头镇（19°23'N，109°35'E）进行。该地区属热带季风气候，年平均气温23.8 ℃，降水量1 786 mm，且热雨同季，台风影响较弱。试验田土壤为砖红壤，土壤pH 5.53、有机质16.84 g·kg^–1、碱解氮93.47 mg·kg^–1、铵态氮14.19 mg·kg^–1、硝态氮10.81 mg·kg^–1、有效磷96.57 mg·kg^–1、速效钾372 mg·kg^–1。

试验田所种香蕉品种为巴西蕉（Mimosa nana Lour.），苗龄2个月，由海南大学儋州校区中国热带农业科学院试验研究基地提供。

2018年10月12日定植香蕉苗，香蕉行距2 m、株距2 m。2018年11月在蕉园中选取香蕉苗长势均匀相邻的两列，每列香蕉植株60棵作为样品采集对象，两列周边香蕉行设为保护行。水分管理：抽蕾前使田间持水量保持在85%左右，抽蕾后使田间持水量保持在75%左右。施肥管理：底肥每株施腐熟的羊粪有机肥5 kg，同时添加磷肥0.2 kg，均匀拌入施肥穴；前期（前3个月）每10~15天追施一次肥，第一个月每株每次淋施400倍尿素水溶液约4 kg，后两个月每株每次淋施200倍硫酸钾复合肥（15-15-15）水溶液约4 kg；中期每15~20天追施一次肥，每株每次淋施66 g尿素、60 g硫酸钾复合肥和166 g硫酸钾；后期每25~30天追施一次肥，每株每次淋施50 g尿素、80 g硫酸钾复合肥和110 g硫酸钾。

分别在抽蕾期前与抽蕾期采集土壤样品。1）抽蕾期前随机选取3棵健康香蕉植株，组成一个混合香蕉根际土。共采集15棵健康香蕉植株，组成5个混合重复样品，记为BH（Healthy plants before heading stage）。2）抽蕾期香蕉大量发病后选取发病率较高的一列，随机采集3棵发病香蕉植株，组成一个混合香蕉根际土。共采集15棵发病香蕉植株，组成5个混合重复样品，记为AD（Diseased plants at heading stage）。3）采集AD样品中每棵患病植株相邻的健康香蕉植株的根际土，同样三棵组成一个混合重复，共采集相邻健康香蕉植株15棵，组成5个混合样品，记为AH（Healthy plants at heading stage）。

根际土样放入冰盒内带回实验室，抖落根系表层土壤，抖落的根围土壤过2 mm筛后风干，用于土壤理化性质的测定；然后用无菌剪刀将根系剪成5 cm左右小段，装入到含有500 mL无菌水的组培瓶内，然后以170 r·min^–1的速度振荡30 min。振荡完成后，用纱布过滤到灭菌后的锥形瓶内，然后将瓶内悬浊液转移到50 ml离心管内，配平后用冷冻离心机4 ℃，10 000 r·min^–1离心3 min，弃去上清液后底部土壤即为根际土^[17]，取约100 g放置到–80 ℃冰箱，用于提取DNA，随后进行高通量测序。

土壤理化性质分析：pH采用电位法测量，水土比为5︰1；通过油浴法测定土壤有机质；0.03 mol·L^–1氟化铵和0.025 mol·L^–1盐酸浸提—钼锑抗比色法测定土壤有效磷；1 mol·L^–1中性乙酸铵溶液浸提—火焰光度计法测定土壤速效钾；2 mol·L^–1氯化钾溶液浸提—紫外分光光度计法测硝态氮；2 mol·L^–1氯化钾溶液浸提-靛酚蓝比色法测铵态氮；具体方法详见《土壤农化分析》^[18]。

土壤可培养病原菌尖孢镰刀菌、细菌、真菌的数量均通过稀释涂布计数法测定^[19]，可培养尖孢镰刀菌计数使用K2培养基，28 ℃培养72 h；细菌采用LB培养基，30 ℃培养24 h；真菌采用马丁氏培养基，28 ℃培养72 h。将培养后计数平板上形成的菌落数转换成每克干土形成的菌落数（Colony forming unit，CFU），以CFU·g^–1干土表示；原生动物数量采用改进的“3级10倍”环式稀释法计数分析^[20]，称取2 g风干土样后加入198 mL蒸馏水，充分震荡1 h，得到稀释度为10²的土壤悬浮液；吸取2 mL悬浮液，加入18 mL无菌水，充分摇匀，此时稀释度为10³；按此方法配制稀释度为10⁴的土壤悬浮液；每级各取1 mL接种于各自的培养环内，每级稀释液接种2个培养皿10个环，将培养皿置于25 ℃光照培养箱中培养，分别在第4、7、11天于显微镜下镜检，并依据“3级10倍”环式稀释法原生动物密度换算表计算每克烘干土壤中原生动物的数量。

取0.80 g混合土壤样品，用土壤DNA提取试剂盒（PowerSoil DNA Isolation Kit）按照说明书上的要求提取每个混合土壤样品的DNA提取。提取好的DNA用核酸定量仪（NanoDrop 2000，Thermo Scientific，USA）测定其核酸浓度，然后存于–20 ℃冰箱备用。每个重复分别吸取20 μL样品送至美格基因公司进行测序，以基因组DNA为模板，根据测序区域的选择，使用带barcode的特异引物及TaKaRa Premix Taq® Version 2.0（TaKaRa Biotechnology Co.，Dalian，China）进行PCR扩增。原生动物18SV4区域采用引物528F（5'-GCGGTAA TTCCAGCTCCAA-3'）和706R（5'-AATCCRAGAAT TTCACCTCT-3'）^[21]；真菌测序ITS1区域采用引物ITS1F（5'-CTTGGTCATTTAGAGGAAGTAA-3'）和ITS2R（5'-GCTGCGTTCTTCATCGATGC-3'）^[22]。PCR反应条件：98 ℃预变性2 min，进入热循环：98 ℃变性15 s，55 ℃，退火30 s，72 ℃延伸30 s，共30个循环，最后72 ℃延伸5 min。每个样本进行3个重复，并将同一样本的PCR产物进行混合。用1%琼脂糖凝胶电泳检测PCR产物的片段长度和浓度。按照NEBNext® Ultra™ II DNA Library Prep Kit for Illumina®（New England Biolabs，USA）标准流程进行建库操作。使用Illumina Nova 6000平台对构建的扩增子文库进行PE250测序。原始数据上传至NCBI Sequence Read Archive（SRA）中，序列号为：PRJNA762227。

根据每个样品所特有的barcode信息，利用QIIME软件将序列分开形成各自的fastaq文件。利用FLASH V1.2.7软件将分开的fastaq文件合并，再利用QIIME软件将低质量的序列（质量分数小于0.5或长度低于200bp或单序列）滤掉^[23]。用USEARCH软件按照UPARSE pipeline分析优质序列，以97%的相似度标注生成OTU（Operational taxonomic units）并筛选出代表序列^[24]。利用RDP classifier（https://pyro.cme.msu.edu）对OTU代表序列进行比对，病原菌使用UNITE Fungal ITS数据库^[25]，原生动物利用PR2数据库（Protist Ribosomal Reference database）进行比对^[26]，删去红藻门（Rhodophyta）、链藻门（Streptophyta）、后生动物门（Metazoa）、真菌门（Fungi）和未分类的后角藻门序列（unclassified Opisthokonta sequences），以生成一个原生动物的OTU分类学综合信息表，用于之后的分析^[11]。根据原生动物行使的功能将其分为：吞噬型（Phagotrophic），光养型（Phototropic），植物致病型（Plant pathogenic），寄生型（Parasitic），腐生型（Saprophytic）和兼养型（Mixotrophic）^[27-29]。

通过R语言的relaimpo包进行可培养微生物与原生动物数量对病原菌种群数量的相对重要性分析，采用IBM SPSS Statistics 20统计分析软件，通过独立样本t检验进行数据比较，检验处理间的差异显著性（P < 0.05）。将USEARCH程序中生成的OTU表格格式进行转化后，用Mothur对转化后的文件进行Alpha多样性的分析，包括Shannon指数，Chao1指数^[30]。基于贝叶斯距离（Bray-distance）进行各处理间的主坐标分析（Principal Coordinates Analysis，PCoA）。相似性分析（ANOSIM，Analysis of similarity）与置换多元方差分析（PERMANOVA，Permutational multivariate analysis of variance）用于判断差异是否到达显著水平（P < 0.05），通过R语言的vegan包完成。原生动物各个功能类型在不同生长期不同处理的相对丰度采用Microsoft Excel 2017软件进行数据处理与作图。不同处理根际土壤中相对丰度前50的原生动物OTU热图通过R语言的gplots包完成。采用R语言中ggcorrplot和corrplot包进行Spearman相关性分析^[31]。

表 1为抽蕾前后可培养细菌、真菌和原生动物数量及其对病原菌种群数量的影响力。抽蕾前可培养细菌、真菌的数量均低于抽蕾后，病原菌在发病植株中的数量显著高于健康植株（P < 0.05），抽蕾后原生动物在发病植株中的数量也高于健康植株。对比可培养细菌与真菌，原生动物的数量对病原菌的种群数量变化的影响力更高。

表 1 Table 1

表 1 土壤可培养微生物数量及对病原菌种群数量的影响力 Table 1 Number of culturable soil microbes and the relative importance for pathogenic Fusarium number 可培养微生物 Culturable microbes 细菌 Bacteria /（log（CFU·g–1 dry soil）） 真菌 Fungi /（log（CFU·g–1 dry soil）） 原生动物 Protist /（log（ind·g–1 dry soil）） 病原菌 Fusarium /（log（CFU·g–1 dry soil）） BH 7.59±4.73b 5.34±3.52c 4.60±3.55a 3.93±2.49b AD 7.81±5.35a 5.81±3.94a 4.64±3.94a 4.20±2.01a AH 7.81±5.38a 5.80±3.91b 4.32±3.67b 3.93±2.03b 相对重要性 Relative importance 20.45% 21.56% 47.19% — 注：BH：抽蕾前健康植株；AD：抽蕾后发病植株；AH：抽蕾后健康植株。不同字母表示不同处理间差异显著（P < 0.05）。下同。Note：BH：Healthy plants before heading stage；AD：Diseased plants at heading stage；AH：Healthy plants at heading stage. Different letters indicated significant differences between treatments（P < 0.05）. The same below.

表 1 土壤可培养微生物数量及对病原菌种群数量的影响力 Table 1 Number of culturable soil microbes and the relative importance for pathogenic Fusarium number

表 2可以看出，抽蕾后香蕉植株根围土壤的pH显著高于抽蕾前（P < 0.05），抽蕾后铵态氮、硝态氮、速效钾均显著低于抽蕾前（P < 0.05），抽蕾后的健康植株根围土中pH、硝态氮、有效磷、有机质显著高于发病植株（P < 0.05），抽蕾后发病植株的铵态氮、速效钾都显著高于发病植株（P < 0.05）。

表 2 Table 2

表 2 香蕉根围土壤的化学性质 Table 2 Soil chemical properties of banana in soil 处理 Treatment pH 有机质 Organic matter /（g·kg–1） 铵态氮 ${\text{NH}}_{\text{4}}^{\text{ + }} - {\text{N}}$ /（mg·kg–1） 硝态氮 ${\text{NO}}_{\text{3}}^ - - {\text{N}}$ /（mg·kg–1） 有效磷 Available phosphorus /（mg·kg–1） 速效钾 Available potassium /（mg·kg–1） BH 6.40±0.02c 14.52±0.99b 15.54±0.33a 12.52±0.63a 114.75±1.79b 80.53±0.77a AD 6.51±0.02b 14.41±0.82b 14.66±0.23b 5.45±0.31c 99.13±2.83c 77.19±1.02b AH 6.74±0.02a 16.47±0.38a 13.56±0.20c 7.92±0.43b 128.33±2.30a 68.97±1.39c

表 2 香蕉根围土壤的化学性质 Table 2 Soil chemical properties of banana in soil

不同样品中香蕉根际土壤的原生动物Alpha多样性指数如图 1所示。抽蕾前的健康香蕉根际中原生动物的香农指数和Chao1指数均显著高于抽蕾期（P < 0.05）。抽蕾期患病植株香农指数显著低于临近的健康植株（P < 0.05），而Chao1指数值与健康植株无显著差异。

图 1 原生动物的Alpha多样性指数：a）香农指数；b）Chao1指数 Fig. 1 Alpha Diversity Indexes of protists in different treatments: a)Shannon index; b)Chao1 index

图 2主坐标分析表明，不同样品之间原生动物群落结构存在显著差异（P < 0.05）。抽蕾前的健康植株与抽蕾后根际原生动物的群落结构在第一主成分上显著分离（P < 0.05），在抽蕾期发病与健康植株的原生动物群落结构在第二主成分上显著分离（P < 0.05）。第一主成分和第二主成分分别占总贡献率的50.02%和30.86%。

图 2 不同处理间原生动物群落主坐标分析 Fig. 2 Principal coordinate analysis of protist communities between treatments

不同样品土壤原生动物各个功能类型相对丰度如图 3所示。无论是在抽蕾前后、发病与健康植株中吞噬型原生动物（Phagotrophic protist）相对丰度均显著高于其他功能类型的原生动物。吞噬型原生动物在抽蕾后相对丰度显著升高（P < 0.05），且在发病植株中最高（59.47%）。光氧型原生动物（Phototropic protist）在抽蕾前相对丰度较高，在抽蕾后相对丰度显著降低（P < 0.05）。植物致病型原生动物（Plant pathogenic protist）在抽蕾后的相对丰度显著升高（P < 0.05），且在发病植株中达到最高（14.97%）。

注：不同字母代表同一功能类型的原生动物在不同处理之间的差异显著（P < 0.05）。  Note：Different letters indicated significant differences between treatments of the same functional taxonomic protist（P < 0.05）. 图 3 各个功能类型在不同生长期不同处理的相对丰度 Fig. 3 Relative abundance of different functional taxonomic protist rhizosphere soil different treatments with different growth periods

在抽蕾前后、发病与健康植株之间原生动物群落组成差异显著（图 4）。所有样品中，吞噬型原生动物的相对丰度都较高，且多属于丝足虫（Cercozoa），其次为光氧型原生动物。抽蕾前健康植株中，相对丰度最高的为光氧型原生动物OTU18（Ochrophyta，Bacillariophyta_X_unclassified），到抽蕾期，无论是健康或是发病植株，OTU18的相对丰度都降低。抽蕾期的发病植株中，相对丰度最高的为一种植物致病型原生动物OTU5（Pseudofungi，Pythium），而抽蕾期健康植株中，相对丰度最高的为吞噬型原生动物OTU16（Cercozoa，Group-Te），且OTU16（Cercozoa，Group-Te）在健康植株中相对丰度显著高于发病植株（P < 0.05）。

图 4 不同处理根际土壤中相对丰度前50的原生动物OTU热图 Fig. 4 Heat map of the top 50 protist OTUs in the rhizosphere soil samples collected from different treatments

不同时期病原菌的相对丰度具有显著差异（P < 0.05），抽蕾前病原菌的相对丰度较高；抽蕾后病原菌在发病植株中的相对丰度高于健康植株，差异具有显著性（P < 0.05，图 5a）。

图 5 不同处理下病原菌的相对丰度（a）及相对丰度前20个原生动物OTU与病原菌的Spearman相关性（b） Fig. 5 The relative abundance of Fusarium in different treatments(a)and the Spearman correlations between the relative abundance of top 20indicator protistan taxa(OTUs)and Fusarium(b)

相对丰度前20个原生动物OTU与病原菌的相关性分析显示，与病原菌呈显著的负相关关系的OTU有9个（P < 0.05），其中吞噬型原生动物有7个，2个属于纤毛虫（Ciliophora），1个属于变形虫（Lobosa），4个属于丝足虫门（Cercozoa），包含了Group-Te、Cercomonas和Filosa-Sarcomonadea属。

与病原菌呈显著的正相关关系的OTU同样有9个（P < 0.05），其中有5个吞噬型原生动物，均属于丝足虫门（Cercozoa），包含了Flectomonas、Neoheteromita和Sandonidae_unclassified属。

长期连续种植同一作物与不合理的施肥方式是香蕉枯萎病发生的重要原因，同时，土壤环境的问题例如土壤板结、土壤容重下降与土壤微生物群落逐渐失衡，这些均会导致土传病原菌的加速繁殖，而根际土壤的微生物群落对作物的生长和健康有重要的影响。为解决这些问题，改善连作蕉园根际土壤微生物群落成为至关重要的一步。原生动物作为土壤微生物群落的重要组成部分，关于其在土壤功能方面的研究一直较少。本研究中，土壤可培养微生物和原生动物对病原菌影响的相对重要性分析结果表明，对比可培养细菌与真菌，原生动物的数量对病原菌种群数量变化的解释率更高，与之前的研究一致^[32]。因此找到某种与病原菌有着负相关关系甚至是抑制作用的原生动物可能会成为一种防治土传病害的有效手段。

之前的研究发现健康与发病的番茄植株中原生动物的群落结构有显著分离^[12]，这与本研究结果一致。本研究中，在根际土壤中原生动物不仅在不同时期的群落结构显著分离，在同一时期，健康与发病植株中原生动物的群落结构也同样差异显著。随着香蕉植株的生长，根际土壤中原生动物的多样性与丰度均呈下降趋势，健康植株中原生动物群落的多样性较发病植株中的更高，这与Ceja-Navarro等^[33]的研究结果相似，他们发现健康柳枝稷根际土壤中的原生动物多样性与丰富度随植物生长而降低。这可能是由于香蕉到达抽蕾期，根际土壤环境如细菌和真菌群落结构等发生改变，原生动物群落内部对细菌真菌等食物的竞争加剧，而某些原生动物的竞争力可能大于其他类群的原生动物，成为优势种，甚至有可能直接捕食其他类群的原生动物^[34]。

本研究中，香蕉抽蕾前后，根际原生动物群落组成发生明显变化，并且抽蕾期发病植株和健康植株根际原生动物群落组成差异明显。无论是在健康还是发病植株中，吞噬型原生动物的相对丰度均显著高于其他功能类群的原生动物。相比于健康植株，发病植株中吞噬型原生动物相对丰度较高。同时，与病原菌具有显著相关关系的原生动物功能类型大部分为吞噬型原生动物，且主要属于丝足虫门，这与Guo等^[35]的研究一致，他们认为丝足虫可能是土壤原生动物群落中最重要的组成。同时Guo等^[35]还指出，原生动物可能是通过与根际中某些微生物群落的互作从而对植物生长起到积极的作用。因此在本研究中，吞噬型原生动物尤其是丝足虫门在健康和发病植株之间的差异，可能与发病植株中的病原菌数量的变化相关^[36]。此外，Xiong等^[12]发现在患病植株中，植物致病型原生动物Pythium属与番茄青枯病病原菌呈显著正相关，与本研究结果相似，说明患病植株可能感染了不同的病原菌，在其中形成了致病型微生物群落，从而使植株更易患病。

本研究中，丝足虫门的部分高丰度原生动物属与病原菌呈现显著的相关关系。丝足虫门中，Group-Te属与病原菌呈显著负相关，在抽蕾期的健康植株中相对丰度最高，显著高于发病植株，Cercomonas属与病原菌也呈显著负相关，且在健康植株中的相对丰度较高，与最新的研究结果相似^[35]。由此可见，丝足虫门部分高丰度属的变化可能会对病原菌的数量产生负面的影响，从而促进植物健康，可能会在防控香蕉枯萎病上产生一定潜力。随后的研究中，有必要将重点放在分离与纯化和病原菌呈负相关的吞噬型原生动物上，有针对性地探索吞噬型原生动物与病原菌的作用机理，以及明确是否需要搭配其他生防菌，才能更好地提高防治香蕉枯萎病的效率。

香蕉根际土壤中原生动物的群落结构和组成在抽蕾前后、发病与健康植株中差异显著。香蕉抽蕾后根际土壤中原生动物群落多样性与丰富度均降低，且发病植株低于健康植株。吞噬型原生动物尤其是丝足虫门在健康和发病植株之间的差异，可能与发病植株中的病原菌数量的变化相关。与病原菌相关的吞噬型原生动物可能在防控香蕉枯萎病上具有一定潜力。