木薯（Manihot esculenta Crantz）是全球9亿多人的主粮,同时也是重要的饲料作物、工业原料作物和生物能源[1]。木薯作为“先锋作物”,在我国农业外交和对外科技援助中发挥着举足轻重的作用。我国生产的木薯绝大部分用于制备乙醇,在玉米等粮食作物作为乙醇生产原料供给不足导致燃料乙醇产业遇到瓶颈的形势下,发展“非粮”燃料乙醇已成为世界燃料乙醇技术的发展趋势[2]。2018年8月22日国务院常务会议明确提出坚持加快建设木薯燃料乙醇的项目,木薯产业发展一直受到国家各级管理部门的高度重视[3]。二斑叶螨（Tetranychus urticae Koch）是世界危险性害螨,其寄主多达50科200余种[4,5]。二斑叶螨是为害木薯最严重的四大有害生物之一,严重暴发成灾主要集中在木薯种植后6~8个月左右,此时的木薯地已封行,化学药剂防治困难[6,7]。
培育与利用抗虫品种是国内外公认的最经济、最有效、最简便的一项带有方向性的害虫防治途径[8,9],而阐明作物的抗虫性机理能够为选育抗虫品种提供理论依据。作物抗虫性是多个基因、酶、蛋白和代谢物共同参与调控的结果,具有典型的数量遗传防御效应特征[10],因此,作物抗虫性的形成原因及可能的调控网络机制需要通过多维度、多因素、多级联关系进行考量,而转录组测序分析是开展作物抗虫性机理研究的强有力工具,已广泛应用于多种作物抗虫基因的筛选、挖掘和功能验证。
水杨酸和茉莉酸介导的防御信号途径在作物抗虫性形成过程中发挥着重要作用,在植物受到害虫为害时能激活体内相应的防御基因表达,可以诱导植物产生影响昆虫消化的防御蛋白或对昆虫有毒害作用的次生代谢物质的直接防御[11,12,13],也可以诱导产生植物释放对昆虫有趋避作用或吸引天敌的挥发物的间接防御[14,15]。当前关于水杨酸和茉莉酸信号途径在拟南芥[16]、烟草[17]等模式植物,水稻[18]、小麦[19]、棉花[20]等大宗经济作物抗虫性的作用机制研究已有较多报道。
当前尚未见害虫为害后抗、感木薯品种转录组差异分析相关报道。因此,本研究通过比较二斑叶螨为害后抗、感螨木薯品种转录组差异,筛选水杨酸、茉莉酸差异表达基因,并通过qPCR验证这两个信号途径关键基因的差异表达,以期初步阐明木薯抗螨性机理,为抗螨木薯品种选育和分子设计育种提供理论基础。
1 材料与方法
1.1 材料
1.1.1 供试木薯 本研究选择遗传稳定、健康无病虫的抗螨参照标准木薯品种C1115和感螨参照标准木薯品种BRA900作为供试木薯[3]。C1115和BRA900均由农业农村部国家木薯种质资源圃提供。
1.1.2 供试虫源 二斑叶螨为本实验室以感螨参照标准木薯品种BRA900植株长期继代饲养的室内试验种群。饲养温度为(26±1) ℃,相对湿度为80%。
1.1.3 供试试剂 多糖多酚植物总RNA提取试剂盒购自于天根生化科技（北京）有限公司,cDNA合成试剂盒、去基因组DNA试剂盒购自宝日医生物技术（北京）有限公司,qPCR试剂盒购自于铂尚生物技术（上海）有限公司。琼脂糖、乙醇、氯仿等均为国产分析纯试剂。
1.2 方法
1.2.1 二斑叶螨接种和取样 二斑叶螨接种和取样参照实验室前期建立的方法进行[21],为保证试验材料的一致性和方法的准确性。选取温室种植3个月、长势一致、无病虫的C1115和BRA900健康植株,发育一致的二斑叶螨雌成螨按50头/叶分别接种于抗、感木薯植株中下部3张叶片背面,再用浸润甘油的棉絮包裹住木薯叶柄基部以防止二斑叶螨在叶片间迁移。分别在二斑叶螨接种为害前（0 d,对照）、接种为害后1、8 d采集叶片用于水杨酸、茉莉酸含量测定及RNA提取。每品种设3个重复,每重复取中下部3张叶片（1个重复3张叶片）。
1.2.2 二斑叶螨为害后转录组差异分析 采用1.2.1的方法进行二斑叶螨的接种和木薯叶片取样。取样后的木薯叶片马上放入液氮中粉碎研磨,并取不同处理各0.5 g样品委托北京诺禾致源科技股份有限公司进行转录组数据测定。去除原始数据的接头序列和低质量reads得到clean data。计算Phred数值大于30的碱基占总体碱基的百分比Q30,当这2个数值均大于90%时表示测序的质量良好,能够满足建库的要求。将clean data与木薯参考基因（ http://www.cassava-genome.cn/index.jsp）进行序列比对,获得mapped data。对样品组中的基因表达量进行差异表达、GO（Gene Ontology）功能注释和KEGG（ https://www.kegg.jp/kegg/pathway.html）通路富集分析。在差异表达基因筛选方面,从差异倍数和显著水平2个水平进行评估,阈值设定一般为|log2(Fold Change)|>1且qvalue<0.005。差异基因GO富集分析方面,选择直观反映生物过程（biological process）、细胞组分（cellular component）和分子功能（molecular function）富集的GO条目上差异基因的个数分布情况,最终挑选富集最显著的30个GO条目作图展示。KEGG富集程度通过rich factor、Q value和富集到此通路上的基因个数来衡量。其中rich factor指差异表达的基因中位于该pathway条目的基因数目与所有有注释基因中位于该pathway条目的基因总数的比值。rich factor越大,表示富集的程度越大。Q value是做过多重假设检验校正之后的P value,Q value的取值范围为[0,1],越接近于零,表示富集越显著,最终挑选富集最显著的20条pathway条目作图展示。
1.2.3 差异表达基因的qRT-PCR验证 以1.2.2转录组分析中筛选出的差异表达的水杨酸、茉莉酸信号途径基因为样本,并结合文献报道,进一步筛选具有指示和标记作用的水杨酸和茉莉酸途径关键差异表达基因进行qPCR验证。采用1.2.1的方法进行二斑叶螨的接种和取样。采用天根公司（TIANGEN）多糖多酚植物总RNA提取试剂盒提取螨害前后的木薯叶片RNA,用Nanodrop和1%琼脂糖凝胶电泳检测RNA的完整性。取去除gDNA后的RNA样品1.0 μg用于cDNA的合成。cDNA样品经无核酸酶水5倍稀释后作为qRT-PCR的模板,分别以木薯Metub为内参基因,相关基因的引物信息如表1所示。qPCR反应条件为：95 ℃预温育1 min后,以40个循环完成如下程序：95 ℃变性15 s,56 ℃退火30 s,72 ℃延伸20 s。分别以来源于未经二斑叶螨为害的木薯叶片cDNA为对照,螨害后的木薯水杨酸、茉莉酸途径基因表达量变化情况以为害前的相对倍数表示,根据Livak等[22]的2-ΔΔCT方法计算而得,每个处理均设置3个重复。
表1 木薯水杨酸、茉莉酸信号途径基因表达量qPCR分析引物信息Tab. 1 Primer information of qPCR analysis of cassava salicylic acid and jasmonic acid signaling pathway gene expression |
| 基因类别 Gene type | 基因名称 Gene name | 基因号 Gene ID | 引物序列（5°-3°） Primer sequence (5°-3°) | 温度 Tm/℃ | 产物大小 Product size/bp |
|---|---|---|---|---|---|
| 水杨酸信号途径基因 | PAL2 | MANES_07G098700 | F:AGAGAGCTCATCAGGTTCTTG | 53.9 | 190 |
| R:ACATGGGGTGATGTTGTTGTTG | 55.8 | ||||
| 4CL3 | MANES_09G127000 | F:TGCTGTTGTCCCGCAAAAAG | 56.7 | 216 | |
| R:GCCAGCTTGGCTTTGATGTC | 57.6 | ||||
| WRKY7 | MANES_05G176800 | F:TTAACAGGGGAAACCGACAG | 54.6 | 163 | |
| R:AAGAGGAACCGCACTTGAGA | 56.1 | ||||
| NPR3 | MANES_14G155900 | F:ACCTGGAAAACAGAGTGGCA | 56.7 | 77 | |
| R:AGCATGTGCGATATCCCTGG | 57.7 | ||||
| 茉莉酸信号途径基因 | JAR1 | MANES_05G025500 | F:TCATCTACTACGCCCAAGTTTAAGA | 58.4 | 118 |
| R:CCATCTCAGCGCCCTTTGTA | 58.5 | ||||
| LOX2 | MANES_10G148800 | F:GCCATGGAAACACTTCTTGCT | 55.9 | 297 | |
| R:CTAGCCCCGTTTTGGGATCA | 57.6 | ||||
| OPR11 | MANES_18G120300 | F:GTTCACACCTCCAAGGCGA | 56.9 | 103 | |
| R:AACTCCATCAAACCCAGCTT | 58.0 | ||||
| 内参基因 | Tub | MANES_18G123607 | F:ATGTTGAGTTCAATATGCAGATCC | 58.1 | 154 |
| R:CCCAGAACTTGGCTCCGATT | 61.7 |
1.3 数据处理
采用SPSS软件进行数据分析,多组数据的显著性差异分析采用One-Way ANOVA-Fisher中的LSD方法,所有数据均为3个生物学重复的平均值,显著性水平为α=0.05。
2 结果与分析
2.1 二斑叶螨为害后转录组差异分析
2.1.1 转录组数据质量评估 由表2可看出,6个处理的18个木薯样品的转录组测序数据总共得到6.93 G的clean data。质量分析结果显示,18个木薯样品Q30碱基百分比分别均超过90.00%,表明18个木薯转录组测序结果的数量与质量都己合格,可进行下一步分析。
表2 转录组数据产出质量评估表Tab. 2 Quality evaluation of transcriptome data output |
| 样品名 Sample name | 原始测序数据Raw reads | 过滤的数据 Clean reads | 过滤的碱基 Clean bases | Q30碱基百分比 Q30 base percentage/% | GC含量 GC content/% | 总体映射（占比） Total mapped (Percentage/%) |
|---|---|---|---|---|---|---|
| 抗-0d-1 | 9 158 864 | 9 075 399 | 0.45G | 97.03 | 43.43 | 8 270 745 (91.13) |
| 抗-0d-2 | 8 416 585 | 8 303 435 | 0.42G | 96.90 | 43.45 | 7 565 295 (91.11) |
| 抗-0d-3 | 9 274 111 | 9 184 437 | 0.46G | 96.97 | 45.12 | 8 371 538 (91.15) |
| 抗-1d-1 | 8 365 650 | 8 289 396 | 0.41G | 97.76 | 43.49 | 7 537 968 (90.94) |
| 抗-1d-2 | 8 022 944 | 7 591 327 | 0.38G | 93.77 | 43.00 | 6 797 283 (89.54) |
| 抗-1d-3 | 8 510 062 | 8 325 590 | 0.42G | 96.91 | 42.94 | 7 466 261 (89.68) |
| 抗-8d-1 | 6 331 530 | 6 170 882 | 0.30G | 93.62 | 43.94 | 5 492 702 (89.01) |
| 抗-8d-2 | 6 406 610 | 6 263 472 | 0.31G | 93.84 | 42.57 | 5 600 715 (89.42) |
| 抗-8d-3 | 6 940 368 | 6 482 502 | 0.32G | 93.75 | 43.74 | 5 729 021 (88.38) |
| 感-0d-1 | 9 329 966 | 8 850 070 | 0.44G | 97.76 | 43.78 | 7 965 958 (90.01) |
| 感-0d-2 | 6 706 724 | 6 475 586 | 0.32G | 93.82 | 43.11 | 5 768 537 (89.08) |
| 感-0d-3 | 6 398 849 | 6 201 887 | 0.31G | 91.66 | 43.94 | 4 812 664 (78.64) |
| 感-1d-1 | 8 714 151 | 8 289 426 | 0.41G | 96.94 | 43.05 | 7 469 774 (90.11) |
| 感-1d-2 | 9 095 073 | 8 584 397 | 0.43G | 96.91 | 43.05 | 7 666 189 (89.30) |
| 感-1d-3 | 6 801 897 | 6 553 439 | 0.33G | 93.93 | 43.22 | 5 119 819 (78.12) |
| 感-8d-1 | 7 227 219 | 6 768 748 | 0.34G | 93.76 | 43.43 | 5 977 438 (88.31) |
| 感-8d-2 | 10 069 575 | 9 816 740 | 0.49G | 96.97 | 41.42 | 8 672 946 (88.35) |
| 感-8d-3 | 7 879 163 | 7 768 360 | 0.39G | 96.85 | 41.30 | 6 782 986 (89.32) |
2.1.2 差异表达分析 用于转录组分析的抗、感螨参照标准木薯品种受螨害1、8 d后的症状如图1所示,随螨害时间的延长,感螨参照标准木薯品种BRA900叶片枯黄褪绿的症状越发明显,而抗螨参照标准木薯品种C1115叶片无显著螨害症状。图2结果表明,与为害前（0 d）相比,受二斑叶螨为害1 d后,抗螨木薯C1115共有1023个基因差异表达,其中540个基因上调表达,483个基因下调表达（图2A）;为害8 d后共有1684个基因差异表达,其中757个基因上调表达,927个基因下调表达（图2B）;与为害前（0 d）相比,受二斑叶螨为害1 d后,感螨木薯BRA900共有2137个基因差异表达,其中1001个基因上调表达,1136个基因下调表达（图2C）;为害8 d后共有1751个基因差异表达,其中907个基因上调表达,844个基因下调表达（图2D）。进一步对抗感木薯品种间的转录组差异表达基因分析表明,受二斑叶螨为害1 d后,抗螨木薯相较于感螨木薯品种共有882个基因差异表达,其中358个基因上调表达,524个基因下调表达（图2E）;受二斑叶螨为害8 d后,抗螨木薯相较于感螨木薯品种共有251个基因差异表达,其中102个基因上调表达,149个基因下调表达（图2F）。二斑叶螨为害前（0 d）、为害1、8 d后,共同差异表达的基因有91个,其中48个上调表达,43个下调表达（图3）。
图1 用于转录组测序分析的抗、感螨参照标准木薯品种及螨害症状Fig. 1 Standard mite-resistant and mite-susceptible cassava cultivars used for transcriptome sequencing and the symptoms caused by mite infestation |
图2 二斑叶螨为害后抗、感螨木薯品种差异表达基因火山图Fig. 2 Volcano plot of differentially expressed genes of mite-resistant and mite-susceptible cassava cultivars |
2.1.3 差异表达基因的GO及KEGG富集分析 图4结果表明,对于不同样本间的比较,GO的功能注释均集中在生物学过程（biological process）和分子功能（molecular function）这两大类。其中受二斑叶螨为害1、8 d后,抗、感螨木薯品种间在生物学过程中富集差异基因个数最多的均为氧化还原过程（oxidation-reduction process）、单有机体代谢过程（single-organism metabolic process）（1 d）、代谢过程（metabolic process）和碳水化合物代谢过程（carbohydrate metabolic process）（8 d）,而为害1、8 d后在分子功能方面富集差异基因个数最多的为氧化还原酶活性（oxidoreductase activity）和催化活性（catalytic activity）。KEGG富集分析方面（图5）,抗、感螨木薯品种差异表达基因分别富集到125条和132条,选取的显著富集的20个条目中,与作物抗虫性关系最密切的通路有次生代谢物质生物合成（biosynthesis of secondary metabolites）、苯丙烷生物合成（phenylpropanoid biosynthesis）、类黄酮生物合成（flavonoid biosynthesis）。
图4 二斑叶螨为害后抗、感螨木薯品种差异表达基因GO功能注释分类统计图Fig. 4 GO functional annotation of differentially expressed genes of mite-resistant and mite-susceptible cassava cultivars infested by T. urticae |
2.2 抗、感螨木薯品种被二斑叶螨为害后差异表达基因的qPCR验证
2.2.1 用于qPCR验证的水杨酸、茉莉酸信号途径基因的筛选 在通过转录组分析获得的二斑叶螨为害后木薯水杨酸、茉莉酸信号途径差异表达基因的基础上,进一步参考前人研究已经明确的具有指示和标记作用的2个信号途径关键基因,最终筛选了木薯水杨酸信号途径基因PAL2、4CL3、WRKY7、NPR3和茉莉酸信号途径JAR1、LOX2和OPR11作为qPCR验证的基因（表2）。
表2 用于木薯水杨酸、茉莉酸信号途径qPCR验证的基因Tab. 2 Genes used for qPCR verification of salicylic acid and jasmonic acid signaling pathway in cassava |
| 序号 No. | 害虫 Pest | 植物 Plant | 基因名称Gene name | 基因功能 Gene function | 参考文献 References | 转录组分析差异表达情况 C1115 (R) vs BRA900 (S) Difference of gene expression as analyzed by transcriptome C1115 (R) vs BRA900 (S) |
|---|---|---|---|---|---|---|
| 1 | 二斑叶螨 | 番茄 | PAL | 启动或加强苯丙烷类代谢,刺激木质素积累 | [23] | 1 d上调,8 d上调 |
| 2 | 茶小绿叶蝉 | 茶树 | 4CL | 促进防御次生代谢物合成,刺激木质素积累 | [24] | 1 d上调,8 d下调 |
| 3 | 大豆蚜 | 大豆 | WRKY7 | 转录调控激活水杨酸途径 | [25] | 1 d上调,8 d下调 |
| 4 | 烟草天蛾 | 烟草 | NPR | 调节水杨酸诱导的免疫反应 | [26] | 1 d上调,8 d下调 |
| 5 | 根结线虫 | 拟南芥 | JAR1 | 催化茉莉酸-氨基酸缀合物合成 | [27] | 1 d上调,8 d上调 |
| 6 | 桃蚜 | 辣椒 | LOX | 植物间诱导防御通讯信号物质,调控茉莉酸合成 | [28] | 1 d上调,8 d上调 |
| 7 | 烟粉虱 | 辣椒 | OPR | 植物防御反应还原性蛋白 | [29] | 1 d上调,8 d上调 |
2.2.2 二斑叶螨为害后水杨酸信号途径差异表达基因qPCR验证 图6A结果表明,抗螨参照标准木薯品种C1115水杨酸信号途径基因PAL2、4CL3、WRKY7和NPR3的表达量表现出随螨害时间延长而显著降低的趋势,其中螨害1 d后表达量达到峰值,4个基因表达量分别为螨害前（0 d）的5.11、3.74、2.40、1.92倍。而感螨参照标准木薯品种BRA900上述4个基因的表达量在螨害后均维持在显著高于螨害前的水平（P<0.05）（图6B）。在同一螨害时间内,相对于BRA900,C1115中PAL2、4CL3和NPR3基因的表达量在螨害1 d时显著上调,而在螨害8 d时显著下调（图6C）,这与转录组分析的基因差异表达结果具有一致性。
图6 二斑叶螨为害前后水杨酸途径基因在抗、感木薯品种中的相对表达量A：C1115;B：BRA900;C：C1115 vs BRA900;不同小写字母表示同一基因的表达量具有显著差异（P<0.05）。 Fig. 6 Changes of salicylic acid pathway gene transcription levels in mite-resistant and mite-susceptible cassava cultivars A: C1115; B: BRA900; C: C1115 vs BRA900; Different lowercase letters indicate significant difference of transcription of identical gene (P<0.05). |
2.2.3 二斑叶螨为害后茉莉酸信号途径差异表达基因qPCR验证 图7A结果表明,抗、感螨参照标准木薯品种被二斑叶螨为害后,茉莉酸信号途径基因的表达量变化趋势并不完全一致。抗螨木薯C1115受螨害后,茉莉酸信号途径基因JAR1、LOX2和OPR11的表达量与螨害前相比均显著提
高（P<0.05）,表达量的最高值均出现在螨害后1 d,分别达到了螨害前的7.16、7.28、8.50倍,随螨害时间的延长,表达量均表现出下降的趋势。而感螨木薯BRA900受螨害后,上述3个基因的表达量也随着螨害时间的延长逐渐降低,并且螨害8 d后JAR1和LOX2的表达量显著降低至低于螨害前的水平（P<0.05）（图7B）。此外,螨害前,抗螨品种C1115中JAR1、LOX2和OPR11的表达量仅分别为感螨品种BRA900的1.15倍、1.14倍和1.57倍,而受螨害1 d和8 d后,上述3个基因在C1115中的表达量均显著高于BRA900（P<0.05）（图7C）。上述茉莉酸信号途径基因qPCR验证的结果与转录组分析的差异表达情况基本一致。
图7 二斑叶螨为害前后茉莉酸途径基因在抗、感木薯品种中的相对表达量A：C1115;B：BRA900;C：C1115 vs BRA900;不同小写字母表示同一基因的表达量具有显著差异（P<0.05）。 Fig. 7 Changes of jasmonic acid pathway gene transcription levels in mite-resistant and mite-susceptible cassava cultivars A: C1115; B: BRA900; C: C1115 vs BRA900; Different lowercase letters indicate significant difference of transcription of identical gene (P<0.05). |
3 讨论
转录组分析结果表明,与未受螨害之前相比,随着螨害时间的延长,抗螨木薯品种C1115的差异表达基因数基本保持不变,而感螨木薯品种BRA900的差异表达基因数则显著下降,并且抗螨木薯相对于感螨木薯的差异表达基因也呈现出随螨害时间延长而显著下降的趋势。推断抗螨木薯品种在螨害期间能够稳定激活自身的防御反应,而感螨木薯品种因螨害后植株严重受损,体内的多个生化反应可能停止或被抑制。这在烟粉虱为害抗感棉花品种[30]、枕叶盲蝽为害抗感水牛草品种[31]的转录组分析中也能观察到与本研究相类似的结果。差异表达基因的GO和KEGG富集分析结果表明,差异表达基因主要富集在次生代谢物质合成、苯丙烷生物合成、类黄酮生物合成和氧化还原相关反应方面,而这些生物学过程又与作物抗虫性息息相关,尤其是苯丙烷合成通路的激活是启动水杨酸信号途径的关键[32],而植物内源茉莉酸亦能够激活类黄酮生物合成途径[33]。因此本研究转录组分析的结果可以总体上大致反映抗、感螨木薯品种抗性差异的分子基础。
水杨酸和茉莉酸途径是植物抵御害虫、病原菌和病毒等有害生物入侵过程中研究最广泛的2条防御信号途径。因此,本研究从转录组分析获得的差异表达基因中选择了具有代表性的指示和标记基因PAL2、4CL3、WRKY7、NPR3（水杨酸途径）和JAR1、LOX2、OPR11（茉莉酸途径）做进一步的qPCR验证。结果表明,上述基因在转录组分析中的差异表达情况和qPCR验证的结果基本一致,并且qPCR分析结果还显示2个途径基因的表达量变化取食并不一致。许多研究表明,双子叶植物中JA和SA防御信号途径既有相互拮抗[34,35],也存在相互协同[36,37],这2个信号途径之间表现出何种相互作用模式和害虫和作物的种类有关。一些研究表明刺吸式昆虫为害后可诱导水杨酸信号途径的激活,而抑制茉莉酸信号途径。如烟粉虱为害利马豆和拟南芥后,植物内源水杨酸含量及其调控的防御基因的表达量均增加,同时抑制内源茉莉酸含量及其调控的防御基因的表达[38]。B型烟粉虱能够诱导拟南芥的水杨酸途径基因并抑制茉莉酸途径基因的表达[39]。烟粉虱为害番茄后,主要诱导水杨酸防御反应,同时抑制茉莉酸防御反应,如何瑜晨[34]研究测定了茉莉酸信号路径上LOX-D基因的表达情况,结果表明,LOX-D在烟粉虱为害12 h时其表达量显著增加,但随着为害时间延长,其表达量逐渐降低,并恢复到对照水平。也有研究表明,刺吸式害虫为害能同时激活植物的水杨酸和茉莉酸途径。如姚新建[40]研究发现抗性小麦品系35-E4和敏感品系35-A20被禾谷缢管蚜为害后,其水杨酸和茉莉酸信号转导途径中相关基因（PAL、LOX、NPR）的表达量均显著提高,但是表达模式存在差异,抗性品系35-E4中相关基因的表达量高于敏感品系35-A20。张羽宇等[41]研究发现西花蓟马为害能显著诱导菜豆植株茉莉酸和水杨酸信号转导途径相关防御酶基因LOX,PAL和PR-2的表达,并在菜豆植株不同部位叶片产生系统抗性。赵丽艳[42]研究发现,蚜虫为害既激活了水杨酸信号途径基因如PAL、4CL的表达,也激活了茉莉酸信号途径基因LOX的表达。上述研究表明不同刺吸式口器害虫为害不同寄主植物后,SA和JA介导的防御反应却显著不同。本研究结果表明,二斑叶螨为害后抗螨木薯C1115和感螨木薯BRA900后水杨酸信号途径基因PAL2、4CL3、WRKY7和NPR3均能够被显著诱导,而抗、感木薯中,茉莉酸信号途径基因JAR1、LOX2和OPR11的表达量在为害早期（1 d）均较为害前显著提高。然而随着为害时间的延长,上述3个基因的表达量虽然在抗螨木薯C1115中逐渐降低,但依然维持在显著高于螨害前的水平,而在感螨木薯BRA900中这3个基因的表达量降低至显著低于螨害前的水平,一方面说明二斑叶螨可能采取了和烟粉虱相同的为害策略,即通过诱导水杨酸途径,抑制茉莉酸途径达到为害感螨木薯的目的,而抗螨木薯的茉莉酸途径均能被诱导产生防御反应免于受害。
本研究从转录组分析、水杨酸、茉莉酸途径关键差异表达基因的qPCR验证角度初步解析了抗、感木薯品种存在抗性差异的可能原因,但仅从基因表达层面仍不能充分阐明木薯的抗螨性机理,后续研究可采用多组学联合分析的深入探讨基因-蛋白-代谢物的共表达调控网络,挖掘并验证关键抗虫基因的功能,以期系统全面阐明木薯品种的抗螨性机理。
责任编辑：黄东杰