0.   引言

【研究意义】松材线虫病又称松树萎蔫病^[1]，被国际公认为最具危险性和毁灭性的松树疫病，是一种重大的检疫性森林病害，号称“松树的癌症”，具有蔓延迅速、发病时间短且致病力强、难防治、毁灭性强等特点^[2]。松材线虫（Bursaphelenchus xylophilus，PWN）是松材线虫病的病原体，主要以墨天牛属（Monochamus spp. ）昆虫为媒介进行传播^[3]，在我国南方地区以马尾松（Pinus massoniana）为主要寄主。探明松材线虫与寄主马尾松间的互作机理是探明松材线虫成灾机制并开发高效防治策略的重要前提。【前人研究进展】目前通过探索寄主、PWN和昆虫媒介之间的生态和遗传联系，开展了大量研究来揭示松材线虫病的致病机制^{[4, 5]}，现已提出了几个关键的致病因素，即寄主中活性氧（ROS）的异常调控、抗性基因的表达、代谢物的调控以及PWN-昆虫载体的共同进化，为松材线虫病的防治提供新的见解^[6−9]。ROS的异常调节是松材线虫病最重要的早期症状^[10]，ROS在生物新陈代谢中具有重要作用，其产生的代谢物质是生物体内一类重要的信号物质^[11]。植物抗病防御相关基因则在植物受侵染后的抗病能力调节中发挥作用^[12]，因此，对马尾松抗性基因的动态表达的研究至关重要，这有助于深入了解这些基因在马尾松抗病防御过程中的作用，并揭露其隐藏的机制。植物相关微生物群落为寄主植物提供了多种有益的优势，如病原体抗性^[13]。研究发现松树的内生细菌和真菌受到PWN入侵的影响，寄主内生菌群的多样性随着松材线虫病的严重程度而增加^{[14, 15]}。此外，一些产生次生代谢物的内生细菌已被证明在松材线虫病的防治中发挥作用，可生产抗生素并作为对抗PWN的生物控制剂^{[14, 16]}。根据分析寄主植物在受到病原物侵染后的代谢产物和代谢途径，能够深入了解植物如何应对松材线虫（PWN）的入侵^[17]。【本研究切入点】截至目前，大多数研究仅收集了单一的组学数据，缺乏多维度多变量的系统完整认知，这些数据仅提供了关于松材线虫病的知识片段^[13]，缺乏组学间的联动。【拟解决的关键问题】以4年生马尾松为研究材料，对有无PWN侵染的马尾松进行多组学协同分析，从多组学角度来探讨马尾松对PWN的防御机制，同步获取微生物群落、基因表达模式以及代谢物信息，建立不同类型数据间的动态关联，对松材线虫病传播多层次和多组学的整体调控具有较好的代表性，明晰PWN入侵前后马尾松体内的ROS变化、抗性反应特征、微生物的群落结构及代谢物变化，以期优化现有的整体防控体系，为松材线虫病的早期诊断提供依据，进而为开发长效防控机制提供理论依据。

1.   材料与方法

1.1   供试材料

1.1.1   马尾松的来源

供试4年生马尾松幼苗由福建省三明市沙县官庄国有林场（26.5603° N, 117.7455° E）提供。

1.1.2   PWN的来源

选择2龄且健康的PWN作为供试PWN^[11]，来源于福建省福州市林业科学研究院。

1.2   试验方法

1.2.1   马尾松的培养

马尾松在接种前先在福建农林大学海峡联合研究院植物生长室中培养2个月，生长室培养条件：光周期16L∶8D，光照度1 000 lx；温度28 ℃（光）/24 ℃（暗），湿度70%。

1.2.2   PWN的培养

首先在超净工作台内，挑取盘多毛孢菌（Pestalotia spp.）菌丝于马铃薯-葡萄糖-琼脂培养基（Potato Dextrose Agar，PDA）平板上，25 ℃下黑暗培养5 d，菌丝基本长满平板。从培养松材线虫的菌板即线虫生长培养基（Nematode Growth Medium，NGM）上挑取松材线虫到盘多毛孢菌的菌丝上，25 ℃黑暗培养，待盘多毛孢菌被线虫取食尽，再次将其转接到PDA上的盘多毛孢菌中扩繁，培养好的PWN可4 ℃保存^[18]。使用贝尔曼漏斗法^[19]从载有PWN的培养基中分离出PWN，24 h后收集于10 mL离心管中，3 000 r·min⁻¹离心5 min，去除上清液。将各离心管中的PWN汇集并混匀，配制成悬浊液。随机吸取1 mL PWN悬浊液于显微镜下观察计数，重复该步骤50次，计算每微升悬浊液中的PWN数量，最终制备成每毫升含5 000条PWN的悬浊液。

1.2.3   马尾松接种松材线虫

取株高、生长状况一致的健康马尾松幼苗用于接种松材线虫，采用人工皮接法^[20]，在距离马尾松基部10 cm的主茎的一侧用无菌刀片斜切出长度约1 cm的伤口，伤口深至木质部，剖开树皮，继而在伤口内固定一个经灭菌处理的脱脂棉球。用封口膜将棉球固定形成斜漏斗状，将树皮和棉球包围其中。用封口膜缠紧漏斗下部，用镊子夹紧接口处和下部，确保不漏水；最后用微量移液器往棉球中注入1 mL的上述松材线虫悬浊液（5 000 条·mL⁻¹），以接种无菌水（ddH₂O）的马尾松作为对照，分别记作PWN（+）组和PWN（−）组。收集处理14 d后的马尾松接种伤口以上10 cm处的茎干及松针，立即放入液氮中并储存在−80 ℃ 冰箱中，用于后续试验。

为验证松材线虫接种效率，采用MoBio PowerSoil DNA分离试剂盒（货号12855-50，MoBio，美国）提取每株马尾松茎干的总基因组DNA。采用NanoDrop 2000光度计（Thermo Fisher Scientific，USA）对DNA的浓度和质量进行检测。完整性分析通过1%（质量分数）的琼脂糖凝胶电泳来完成。提取的DNA在−20 ℃下储存备用。将上述马尾松幼苗分离出的DNA作为模板，通过PCR反应进一步验证，扩增出条带代表有效接种。以Internal control作为内参引物，以Bx-cathepsin作为扩增的目的引物^[21]并列于表1。

表  1  特异引物

Table  1.  Specific primers

引物名称 Primer name	正向引物 Forward primer（5′-3′）	反向引物 Reverse primer（5′-3′）
内对照 BxPP2A	AACGTCATTTCTAGCCGCCA	TCAGCCCTACAAACCCCTCT
Bx-cathepsin	TTGCATTCTACGGCCAGTCC	ACTGACTTTCGATGGCTCCG
c60547.graph_c0	TAAATTCCAAGTGCCCCGCA	ACCGTGATACACATTTCAGA
c64867.graph_c0	AAATCGTGTGTGTCCCTGCA	GGTTGCAATGATAACGGCCC
c68789.graph_c0	CGCCCGAATCTCTGCACTTA	TCGATGGTCTTGGTGATGGC
c81022.graph_c0	TTGGCTGTACAGATTCCCGT	ACCTATGGATGTCTGCTCCA
c82953.graph_c0	ACTGTTAACCTGGCTCACGG	CTACGCAAATTCACCGCCAC

下载: 导出CSV 

| 显示表格

1.2.4   寄主马尾松中氧化应激指标

取PWN（+） 组和PWN（−） 组的马尾松针叶研磨成粉末，ROS定量参照苏州格锐思生物科技有限公司的ROS含量试剂盒（货号ROS-1-Y，科铭生物技术有限公司，中国苏州）。H₂O₂定量试验参照苏州格锐思生物科技有限公司的过氧化氢含量试剂盒（货号G0112W，微板法96样）。采用SPSS 26.0进行单因素方差分析（ANOVA Tukey）来确定各组之间的差异。

1.2.5   寄主马尾松中抗性基因转录水平的定量分析

用Trizol试剂（Invitrogen，美国）从PWN（+）组和PWN（−）组马尾松针叶的总RNA，通过PrimeScriptTM RT试剂盒（TaKaRa，日本）逆转录为cDNA。马尾松体内抗松材线虫基因参照蔡守平等^[21]研究中筛选出来的抗性基因（编码应激反应通路的c60547.graph_c0、c82953.graph_c0，萜类生物合成途径通路的c64867.graph_c0、c68789.graph_c0及合胞体形成通路的c81022.graph_c0）合成相应的检测引物，根据 HieffTM qPCR SYBR Green Master Mix（Low Rox Plus，叶森，中国上海）操作步骤，使用QuantStudio 6 Flex PCR（ABI，MA，美国）扩增仪进行实时荧光定量聚合酶链反应（qRT-PCR）分析，特异性引物和参考基因列于表1中。应用2^−ΔΔCT方法计算PST相对表达水平^[22]。设置3个生物学重复，3个技术性重复，通过单因素方差分析（ANOVA Tukey检验）确定各组之间差异的显著性水平，使用R软件（版本4.1.1）绘制聚类分组热图。

1.2.6   寄主马尾松样品的宏基因组测序与分析

将PWN（+） 组的马尾松针叶取5个样品等量混合作为1个试验样本，PWN（−）组作相同处理，将混合好的2组样本在液氮中粉碎。使用 MoBio PowerSoil DNA 分离试剂盒（12855-50，MoBio，United States）从每个样品中提取总基因组 DNA， 采用NanoDrop 2000分光光度计（Thermo Fisher Scientific，United States）测量DNA的数量和质量。然后，通过超声破碎机将DNA破碎成300 bp左右的片段，作为测序文库。而后，所有DNA测序文库在北京奥维森基因科技有限公司的Illumina HiSeq平台上进行高通量测序。测序过程中，使用Illumina Analysis Pipeline v2进行图像分析、碱基读取和误差评估。测序完成后使用Trimmomatic 软件对原始数据进行质控，包括去除接头序列和低质量数据。保留的高质量测序数据与NR数据库进行比较，得到样品中微生物群落的组成信息；此外，通过计算Alpha多样性、香农 （Shannon-Wiener index）和辛普森 （Simpson index）多样性指数表示寄主马尾松中微生物的多样性；以上分类和多样性计算均使用DIAMOND 中的相关软件包完成^[23]。

1.2.7   代谢组测序和分析

代谢组测序使用的样品与宏基因组测序一致，样品的准备按Jia 等^[23]的方法进行。样品经Agilent 1290 Infinity LC 超高效液相色谱系统（UHPLC）分离后，用 Triple TOF 6600 质谱仪进行质谱分析，分别采用电喷雾电离正离子和负离子模式进行检测。设置流动相的流速为0.35 mL·min⁻¹，流动相为0.1%甲酸水溶液（A）（A117-50，Thermo Fisher Scientific，Waltham，MA，USA）和0.1%甲酸的乙腈溶液（B）（A998-4，Thermo Fisher Scientific，Waltham，MA，USA），柱温为45 ℃，自动进样器温度为4 ℃，进样体积为5 μL。而后，按照Jia等 ^[23] 的方法，从测序数据中对代谢物进行定量和主成分分析，并筛选差异代谢物。

2.   结果与分析

2.1   松材线虫接种效果的评估

松材线虫的接种效果如图1所示，PWN（+）_1、PWN（+）_2、PWN（+）_3（接种松材线虫的马尾松）扩增出松材线虫目的基因的特异性条带，且条带大小与阳性对照PWN（松材线虫）扩增出的特异性条带大小相同，而PWN（−）（注射ddH₂O的马尾松）和阴性对照双蒸水（ddH₂O）没有扩增出特异性条带，说明本试验使用的人工皮接法可以使松材线虫成功接入到马尾松幼苗中。

图  1  松材线虫接种效果的验证

以PWN作为阳性对照，以ddH₂O作为阴性对照。

Figure  1.  Validation of PWN inoculation efficiency

PWN was used as a positive control and ddH₂O was used as a negative control.

下载: 全尺寸图片 幻灯片

2.2   松材线虫入侵引起寄主马尾松氧化应激指标的变化

ROS作为植物抗性反应中一种十分重要的信号物质^[24]，H₂O₂是植物光呼吸过程中的代谢产物，具有较高的跨膜通透性且化学性质较为稳定，当植物受到生物或非生物刺激时，它能够诱导植物体内抗氧化保护酶活性的改变，是ROS中主要的信号分子，在防卫反应中扮演重要角色^[11]。本研究在接种第14天测定ROS含量以量化松材线虫的致病性，结果表明，PWN（+）（接种松材线虫的马尾松）体内的ROS含量和H₂O₂含量均显著高于PWN（−）（接种ddH₂O的马尾松）（P＜0.05），其中，ROS含量提高了3.2倍，H₂O₂含量提高了1.7倍（图2）。

图  2  有无接种松材线虫对寄主马尾松ROS（a）和H₂O₂（b）含量的影响

不同小写字母代表差异显著（P＜0.05）。

Figure  2.  Effects of PWN inoculation on ROS (a) and H₂O₂ (b) contents in P. massoniana

Different lowercase letters represent significant difference at P＜0.05.

下载: 全尺寸图片 幻灯片

2.3   松材线虫入侵对寄主马尾松抗性基因表达量的变化

通过对2组马尾松体内抗松材线虫基因表达量进行定量，结果表明c60547.graph_c0、c82953.graph_c0在PWN（+）组中的表达水平显著高于PWN（−）组（P＜0.05），c64867.graph_c0、c68789.graph_c0、c81022.graph_c0在PWN（+）组中的表达水平显著低于PWN（−）组（P＜0.05）（图3）。

图  3  松材线虫入侵前后马尾松体内抗性基因的变化

Figure  3.  Changes of resistance genes in P. massoniana before and after PWN inoculation

下载: 全尺寸图片 幻灯片

2.4   松材线虫入侵对寄主马尾松宏基因组的变化

Alpha多样性主要用于衡量一定生境内或样品内生物种类菌群多样性情况，PWN （−）组的Shannon指数（图4a）、Simpson指数（图4b）、Invsimpson指数（图4c）均显著高于PWN （+）组，说明松材线虫入侵后会显著降低马尾松体内微生物的多样性和物种集中性。

图  4  不同处理条件下马尾松微生物群的多样性和结构功能

不同处理组微生物群落的多样性Shannon（a）、Simpson（b）、Invsimpson（c）指数差异，不同小写字母代表差异显著（P＜0.05）；（d）代表样本在目水平上的聚类柱状图；（e）样本微生物在KEGG level3水平上具有显著性差异的功能预测。

Figure  4.  Diversity and structural functions of microflora in P. massoniana under varied conditions

Diversity of microbial communities in different sample groups differed in Shannon (a), Simpson (b), and Invsimpson (c) indexes; different lowercase letters represented significant differences at P＜0.05; d: clustering histogram of samples at order level; e: functional prediction of microorganisms at KEGG level 3.

下载: 全尺寸图片 幻灯片

根据宏基因组测序分析得到的有效序列注释结果，选取各样品在目水平上相对丰度高于1%且前10的物种，生成物种累积相对丰度图，结果表明PWN （+） 组和PWN （−）组的优势目（图4d）均为肉座菌目（Hypocreales）＞肠杆菌目（Enterobacterales）＞伞菌目（Agaricales）＞伯克氏菌目（Burkholderiales），但单胞菌目（Hyphomicrobiales）的相对丰度在PWN （+） 组＜PWN （−）组。

如第三层级KEGG功能注释显著差异功能堆积柱状图（图4e）所示，PWN （+）组有3条通路丰度高于PWN （−）组：复制修复通路（Replication and repair）、DNA复制通路[DNA replication （PATH:ko030300）]、DNA复制蛋白通路[DNA replication proteins（BR:ko03032）]。

2.5   松材线虫入侵对寄主马尾松代谢组的变化

对2组马尾松代谢物的化学成分进行主成分分析，结果表明两组样本得分差异显著，样本全部处于95%置信区间内，第一主成分（PC1）的贡献率为15.6%，正交主成分（PC2）得分为81.6%，两个主成分的贡献率之和为97.2%（图5a）。在两个处理组的马尾松中共检测到2 072个代谢物。将筛选差异代谢物的结果以火山图（volcano plot）的形式进行可视化（图5b），共筛选到716个差异代谢物，其中PWN （+）组相对于PWN （−）组上调的差异代谢物有365个，占比51%，下调的差异代谢物有351个，占比49%。不同代谢物在PLS-DA模型第一主成分的变量投影重要度（Variable Importance in the Projection, VIP）显示前20个代谢物，说明两个处理组的样本间差异较大（图5c），其中3-Hydroxy-1-(3-hydroxy-4-methoxyphenyl)-2-{4-[(1E)-3-hydroxy-1-propen-1-yl]-2-methoxyphe-noxy}propyl-beta-D-glucopyranosid、根皮素（Phloretin）、熊去氧胆酸（Ursodiol）、鼠李蒽酚（Cascaroside D）、舍吲哚（Sertindole）、促微丝聚合剂（Jasplakinolide）、17-trifluoromethylphenyl trinor PGF2alpha isopropyl ester、羧苄青霉素（Carbenicillin）、前列腺素C2（Prostaglandin C2）这9个代谢物在PWN （+）组的含量中显著高于PWN （−）组中。KEGG通路的气泡图（只展示top 20的结果）如图5d所示，其中差异代谢物主要富集在ABC转运蛋白（ABC transporters）通路、花生四烯酸代谢（Arachidonic acid metabolism）通路、类黄酮生物合成（Flavonoid biosynthesis）通路、甘油磷脂代谢（Glycerophospholipid metabolism）通路等。聚类热图分析显示PWN （+）和PWN （−）两个处理组间诱导的代谢物组成不同（图5e）。

图  5  松材线虫入侵前后诱导的马尾松代谢物组成

（a） 不同处理的OPLS-DA图；（b）不同处理间差异代谢物筛选火山图；（c）不同处理的差异代谢物VIP散点图。显示基于VIP值的前20种差异代谢物；(d)不同处理的代谢通路分析气泡图。（e）不同处理的所有差异代谢物的层次聚类分析热图。

Figure  5.  Metabolites in P. massoniana before and after PWN inoculation

(a) OPLS-DA plots for different treatments; (b) differential metabolite screening volcano plot between different groups; (c)VIP scatter plots of differentially differentiated metabolites in different groups; top 20 differential metabolites based on VIP values are displayed; (d)Bubble map of metabolic pathway analysis in different groups;(e) Hierarchical cluster analysis heat map of all differential metabolites of groups.

下载: 全尺寸图片 幻灯片

3.   讨论与结论

松材线虫是当前林业发展面临的主要威胁之一，松树在受到松材线虫侵染后，其表型组、转录组、微生物组及代谢组均会产生明显变化，共同激发寄主马尾松的抗逆反应^{[1, 3]}。在众多寄主植物应对病原物的防卫反应中，作为信号分子的H₂O₂、ROS等会在病原物入侵后迅速产生并参与寄主植物对病原物的过敏性反应^[25]。叶建仁^[2]曾指出当松材线虫入侵寄主后，寄主会通过产生大量的活性氧（ROS）和萜烯类物质等化学防御物质来抵抗或抑制松材线虫的入侵。俞禄珍等^[26]的研究发现黑松受松材线虫侵染后树体内ROS活性和NO含量均升高，胡龙娇等^[27]研究表明H₂O₂含量的积累对松树抗病性具有重要的影响，何龙喜等^{[11, 28]}的研究表明不同抗性松树体内H₂O₂、ROS含量等活性的变化差异均与其抗性密切相关。陈玉惠等^[29]研究表明在不同抗性松树与松材线虫互作早期，仅有感病松树接种松材线虫后H₂O₂不断上升，意味着活性氧及其相关酶产生特异变化在松材线虫病的发生发展中起着相当关键的作用。本研究通过松材线虫侵入寄主马尾松后的表型组学数据进行测定，发现PWN（+）组中的ROS及H₂O₂含量均会显著增高（P＜0.05），这与前人的研究结果一致，说明松材线虫侵入马尾松体内后，寄主植物马尾松会产生防御反应，ROS和H₂O₂含量是马尾松受松材线虫侵染后介导抗病信号转导的重要信号物质。

基因表达分析为深入了解基因功能提供了关键的信息，已经在多个领域得到了广泛应用，其深入研究探索了疾病、非生物胁迫等相关基因表达和调控机制，揭示了生命活动的变化规律^[30]。本研究利用qPCR检测松材线虫侵染下的马尾松中相关抗性基因的动态表达情况，结果发现PWN（+）组中c60547.graph_c0和c82953.graph_c0这两个基因的表达水平均显著（P＜0.05）高于PWN（-）组，说明属于编码应激反应通路的这两个抗性基因在马尾松抗松材线虫病中与ROS呈正相关；而萜类生物合成途径通路的c64867.graph_c0和c68789.graph_c0以及合胞体形成通路c81022.graph_c0在PWN（+）组的表达水平则显著（P＜0.05）低于PWN（-）组，说明这3个抗性基因在马尾松抗松材线虫病中与ROS呈负相关，将不利于马尾松的解毒作用而导致植株系统抗性的不断减弱，可见本试验中马尾松的基因表达量会受松材线虫的影响。这与Liu等^[30]的研究结果一致。同时Han等^[31]的研究也表明，松树会响应松材线虫的抗性基因发生表达，意味着松材线虫的入侵能够激发马尾松幼苗启动一系列复杂的早期防御机制。谢婉凤等^[32]也研究了松材线虫侵染对马尾松基因表达的影响，发现松材线虫侵染致使马尾松中的信号传递与抗病防御相关的基因表达受到抑制，破坏其系统抗性。

物种多样性是生物多样性的核心，反映了生物群落结构的复杂性，是体现群落结构发展阶段及稳定性的一种表现形式^[33]。Shannon和Simpson指数分别反映物种多样性大小和集中程度，是评价微生物群落多样性的重要指标^[34]。本研究通过对松材线虫入侵寄主植物马尾松进行宏基因组分析发现，PWN（−） 组的微生物多样性均高于PWN（+）组，这表明由于松材线虫的侵染降低了马尾松微生物的多样性，温晓健等^[35]通过对松材线虫侵染前后马尾松树体内微生物多样性进行分析，发现健康马尾松树体内的真菌多样性高于枯死马尾松，但健康和枯死马尾松的细菌群落具有一定的相似性，这与本研究结果一致。松树体内的微生物种群与松材线虫病的发生和发展密切相关，本研究结果表明优势真菌是肉座菌目（Hypocreales）且相对丰度PWN（+）组＞PWN（−）组，Al Khoury等^[36]研究发现肉座菌目（Hypocreales）长期以来一直被认为是害虫防治剂的重要微生物，朱兆香等^[37]研究表明肉座菌目（Hypocreales）具有较强的杀伤能力，因而作为生物防治制剂广泛使用。可见PWN侵染后马尾松Hypocreales相对丰度的增加，可能是马尾松应对PWN的生物防治策略。

通过对两组马尾松代谢物进行分析，发现PWN（+）组高于PWN（−）组的差异代谢物主要有根皮素、熊去氧胆酸、羧苄青霉素等。根皮素是一种二氢查耳酮类黄酮，它可提高机体的抗氧化能力，抵御细胞毒性^[38]，抑制体内自由基异常积累^[39]，本研究结果中PWN（+）组中氧化性指标物质ROS及H₂O₂含量显著高于PWN（−）组，从而说明根皮素的升高可保护植物免受氧化损伤。熊去氧胆酸对脂质代谢至关重要，可减少细胞周期调控蛋白的表达，具有降脂和免疫调节等重要特性^[40]。羧苄青霉素是一种半合成青霉素抗生素，对革兰氏阴性菌具有广谱的抗菌活性，可以干扰细菌的细胞壁合成，同时对植物组织显示出低毒性^[41]。差异显著的代谢物显著富集于ABC转运蛋白通路、花生四烯酸代谢通路、类黄酮生物合成通路、甘油磷脂代谢通路。其中，ABC转运蛋白通路功能为主动向胞外排放种类繁多的内源或外源毒素，包括真菌自身产生的侵染寄主植物所需的真菌毒素、其他拮抗微生物所产生的抗生素和植物防御性物质等^[42]。花生四烯酸调节多种脂质信号介质的生理反应和病理过程，包括细胞增殖、凋亡、代谢和迁移，具有促炎症和促细胞凋亡作用^[43]。植物体中类黄酮合成途径是目前研究的最为清楚的次生代谢途径，周科等^[44]表明具有抗病性的类黄酮生物合成通路，与植物抗逆境胁迫直接相关。赵琳儒等^[45]曾在研究中表明甘油磷脂代谢通路在细胞膜融合、胞饮作用、膜运转作用及膜中酶的催化活性等方面起着重要作用，同时也在维持细胞膜结构和功能过程中发挥着优势。这说明松材线虫的入侵，激发了寄主马尾松产生与抗病和抗逆有关的防御性次生代谢物来对抗各种胁迫。

综上，本研究结果表明当松材线虫入侵马尾松时，寄主会启动一系列复杂的防御反应，这些反应并非孤立存在，而是通过多种机制协同作用，共同激发寄主马尾松的抗性反应，应对松材线虫的侵染。这对于全面深入认识松材线虫病的致病机理和寻找新的防治途径有着十分重要的意义。在未来的工作中，我们将通过筛选一系列能够在松树体内能够轻易定植并大量繁殖的功能性微生物或次生代谢物，通过人工外源导入的方式实现植物抗性的提升却不占有植物生长发育的资源，从而提供更有效的大规模虫害综合治理策略。