-
绒毛白蜡为木犀科(Oleaceae)落叶乔木,原产于北美,现在我国华北、内蒙古和长江中下游地区均有栽培,是重要的木材和观赏物种。白蜡窄吉丁(Agrilus planipennis Fairmaire)属鞘翅目吉丁总科吉丁科(Coleoptera: Buprestidae),是木犀科梣属(Fraxinus)树木的一种重要蛀干害虫[1],其幼虫在树木的韧皮部、形成层和木质部浅层蛀食为害,严重危害我国北方地区引种栽植的各个白蜡树种,如洋白蜡(Fraxinus pennsylvanica Marsh.)和绒毛白蜡(Fraxinus velutina Torr)[2]。白蜡窄吉丁幼虫的危害极具隐蔽性,其在树皮内蛀食为害时不将碎木和虫粪推出坑道,从树皮外很难发现树木受害。成虫的飞翔能力较弱,新羽化的成虫常常在同一棵树上或附近树上继续产卵为害[3]。因此,成虫一旦侵染树干成功,危害会逐年加重,通常1~3年即可导致树木死亡[4]。由于白蜡窄吉丁幼虫生活的隐蔽性,常规的防治方法很难奏效,所以目前对于白蜡窄吉丁的防治主要集中在成虫期,一般采用诱捕器、粘虫板、阻虫网以及药剂喷雾等。目前国内针对白蜡吉丁幼虫天敌的研究较多,主要是为利用天敌昆虫的生物防治服务[5-6]。
植物在长期进化过程中,形成了通过物理结构和有毒次生代谢产物抵御植食性昆虫危害的机制[7]。次级代谢产物的应答策略可以通过分子水平上的变化来解析其产生的分子机制。目前,国内外对于白蜡树抗白蜡窄吉丁的相关研究相对较少,尤其是针对白蜡树响应害虫侵害时的韧皮部组织中转录组水平变化的研究更是鲜见报道。鉴于此,作者对绒毛白蜡与白蜡窄吉丁幼虫互作过程中韧皮部的分子表达谱进行了研究。通过对健康与受害白蜡树的韧皮部材料进行转录组水平测序、组装和注释,获得了寄主树与幼虫互作过程中被激活的各类转录因子家族,并分别对差异表达基因进行GO (Gene Ontology)、KEGG(Kyoto Encyclopedia of Genes and Genomes)注释和富集分析。研究结果将为绒毛白蜡树抗虫分子机制研究奠定理论基础,同时也为丰富木犀科植物抗虫相关基因库,促进白蜡抗逆分子育种和优良品系培育提供重要科学依据。
-
分别取健康(Non-infested 1-3)和受害(ASF 1-3)的白蜡树韧皮部,对样品文库进行了RNA-Seq测序,分别获得Raw Reads 49287036、44918932、50780898、49722100、45082112和45802214条,去除有接头(Adapter related)、N含量超过10% (Containing N)以及低质量(Low quality)的reads,最终得到(Clean reads)分别为48778052、44592026、50393962、49349330、44708080和45432400条;去除Phred值小于20的碱基片段,样本高于Q20测序结果比例分别为98.11%、98.34%、98.3%、98.29%、98.3%和98.25%;碱基G和C的数量总和占总的碱基数量的百分比分别为44.75%、43.58%、43.3%、43.57%、43.4%和43.37% (表1)。结果表明测序质量合格,可进行下一步的生物学分析。
-
对健康与受害白蜡树韧皮部进行样品间相关性及基因表达差异分析,结果显示,健康与受害白蜡树韧皮部之间相关性较低,均分布在0.7~0.8之间,而各自生物学重复间的相关性相对较高,均高于0.9,说明我们前期分别采集健康与受害白蜡树韧皮部并设置3组重复的实验设计较为合理(图1A)。进一步对健康与受害白蜡树韧皮部的表达差异基因通过火山图进行分析,显著上调的基因以红色表示,显著下调的基因以绿色表示,无显著性差异的基因以灰色表示(图1B)。共得到了3388个差异表达基因(DEGs),其中受害相对于健康白蜡树韧皮部表达下调的DEGs有1247个,表达上调的DEGs有2141个,上调的基因数明显高于下调基因数。结果表明,虫害胁迫下白蜡树韧皮部的大部分基因受到激活,少数基因受到抑制。
表 1 转录组优化组装结果评估
Table 1. Evaluation of transcriptome optimized assembly results
样品名称
Sample原始序列数据
Raw reads过滤后测序数据
Clean reads碱基错误率
Error rate/%Phred
> 20的百分比/%GC百分比
GC content/%ASF-1 49722100 49349330 0.02 98.29 43.57 ASF-2 45082112 44708080 0.02 98.3 43.4 ASF-3 45802214 45432400 0.02 98.25 43.37 Non-infested1 49287036 48778052 0.02 98.11 44.75 Non-infested2 44918932 44592026 0.02 98.34 43.58 Non-infested3 50780898 50393962 0.02 98.3 43.3 -
为明确健康与受害绒毛白蜡树韧皮部差异表达基因的生物学功能,对这些基因进行GO功能注释分析,将差异表达基因分为生物学过程(biological process)、细胞组成(cellular component)和分子功能(molecular function)三大类。差异表达基因可主要归纳于9个生物学过程、7个细胞组分以及4个分子功能(图2)。在生物学过程中以细胞过程(cellular process)以及代谢进程(metabolic process)差异表达基因数量居多;在分子功能中差异表达基因集中于催化活性(catalytic activity)和结合元件(binding)两个方面,另外,转运活性(transporter activity)以及核酸结合转录因子活性(nucleic acid binding transcription factoractivity)也出现了基因差异表达(表2)。表明白蜡树韧皮部在白蜡窄吉丁胁迫下主要通过以上途径的差异表达基因来进行相关应答及调控。
表 2 差异基因GO功能注释分类统计
Table 2. Statistical table of functional annotation classification of differential gene GO
GO 术语
GO Term差异基因表达数
Seq Number信号 signaling 97 细胞成分组织或生物合成
Cellular component organization or biogenesis184 刺激响应 Response to stimulus 214 本地化 Localization 220 生物过程调节 Regulation of biological process 315 生物调节 Biological regulation 328 单个有机体过程 Single-organism process 751 代谢过程 Metabolic process 1 202 细胞过程 Cellular process 1 217 高分子络合物 Macromolecular complex 231 细胞器部件 Organelle part 245 细胞器 Organelle 741 单元部件 Cell part 1 031 细胞 Cell 1 058 膜部件 Membrane part 1 225 膜 Membrane 1 274 核酸结合转录因子活性
Nucleic acid binding transcription factor activity143 转运活性 Transporter activity 247 结合元件 Binding 1 953 催化活性 Catalytic activity 1 963 -
为了系统分析健康与受害白蜡树韧皮部差异表达基因的基因功能、联系基因组信息和功能信息,作者将差异表达基因按照参与的pathway通路或行使的功能进行分类,共建立了20条代谢通路(pathway)。图3表示健康与受害白蜡树韧皮部差异基因KEGG代谢途径,由图可知,这些差异基因参与六大类代谢,分别为代谢(Metabolism)、遗传信息处理(Genetic Information Processing)、环境信息处理(Environmental Information Processing)、细胞过程(Cellular Processes)、生物体系统(Organismal Systems)和人类疾病(Human Diseases)。这些差异基因主要参与代谢通路较多,其中碳水化合物代谢(Carbohydrate metabolism)的差异表达基因最多,达到6397个,其次是氨基酸代谢 (Amino acid metabolism),达到3632个,另外这些差异基因参与的其他代谢通路有能量代谢、脂质代谢以及其他次生代谢物的生物合成(表3)。结果表明,这些基因在多个方面较多的参与受害白蜡树响应白蜡窄吉丁胁迫的生理过程,而次级代谢产物极有可能在此过程中发挥重要作用。
表 3 差异基因KEGG代谢途径分类统计
Table 3. Classification and statistics of metabolic pathway of differential gene KEGG
KEGG术语
KEGG Term差异基因表达数
Seq Number核苷酸代谢 Nucleotide metabolism 1 655 其他次生代谢物的生物合成
Biosynthesis of other secondary metabolites985 糖的生物合成和代谢 Glycan biosynthesis and metabolism 877 萜类化合物和聚酮类化合物的代谢
Metabolism of terpenoids anf polyketides1 052 脂质代谢 Lipid metabolism 2 594 能量代谢 Energy metabolism 2 880 其他氨基酸的代谢 Metabolite of other amino acids 1 716 代谢的辅助因子和维生素
Metabolite of cofactors and vitamins1 610 氨基酸代谢 Amino acid metabolism 3 632 碳水化合物代谢 Carbohydrate metabolism 6 397 转录 Transcription 2 465 翻译 translation 7 022 复制和修复 Replication and repair 884 折叠、分类和降解 Folding,sorting and degradation 5 179 信号转导 Signal transduction 2 093 膜运输 Membrane transport 525 运输和分解代谢 Transport and catabolism 3 561 环境适应 Environmental adaptation 1 476 耐药性:抗菌 Drug resistance:antimicrobial 8 内分泌及代谢性疾病 Endocrine and metabolic diseases 209 -
为获取健康与受害白蜡树韧皮部差异表达基因的主要参与代谢通路,作者对集中的差异基因进行KEGG功能富集分析,发现健康与受害白蜡树韧皮部差异基因主要在122条通路(pathway)中显著富集。图4表示q-value值最小(0.0~0.2)的前20条pathway,差异表达基因数量以点的大小来表示,不同的q-value范围以点的颜色区分。由图可知,内质网蛋白加工(Protein processing in endoplasmic reticulum)的差异表达基因最多,达到78个;其次是内吞作用 (Endocytosis),达到75个(表4)。另一方面,植物激素信号转导(Plant hormone signal transduction)、植物-病原互作(Plant-pathogen interaction)、MAPK信号通路-植物(MAPK signaling pathway-plant)、RNA运输(RNA transport)、氨基糖和核苷酸糖代谢(Amino sugar and nucleotide sugar metabolism)、淀粉和蔗糖代谢(Starch and sucrose metabolism)、嘌呤代谢(Purine metabolism)的q-value值最小,说明这些通路富集程度最高。结果表明,虫害胁迫时以上通路中相关基因的差异表达最为活跃。
图 3 差异基因KEGG代谢途径分类统计柱状图
Figure 3. Histogram of metabolic pathway classification of differential gene KEGG
表 4 绒毛白蜡受白蜡窄吉丁危害下差异基因富集程度排名前20的pathway条目
Table 4. The top 20 pathway entries of differential gene enrichment under pest stress of ash tree
KEGG术语
KEGG Term差异基因表达数
Seq NumberKEGG编号
KEGG ID内质网蛋白加工 Protein processing in endoplasmic reticulum 78 map04141 内吞作用 Endocytosis 75 map04144 植物激素信号转导 Plant hormone signal transduction 70 map04075 植物-病原互作 Plant-pathogen interaction 46 map04626 MAPK信号通路-植物 MAPK signaling transduction 44 map04016 RNA运输 RNA transport 39 map03013 氨基糖和核苷酸糖代谢 Amino sugar and nucleotide sugar metabolism 31 map00520 淀粉和蔗糖代谢 Starch and sucrose metabolism 29 map00500 嘌呤代谢 Purine metabolism 28 map00230 mRNA监测通路 MRNA monitoring pathway 26 map03015 RNA降解 RNA degradation 25 map03018 苯丙素的生物合成 Biosynthesis of phenylpropanoid 23 map00940 卟啉和叶绿素代谢 Porphyrin and chlorophyll metabolism 21 map00860 果糖和甘露糖代谢 Fructose and mannose metabolism 21 map00051 氧化磷酸化 Oxidative phosphorylation 21 map00190 丙酮酸代谢 Pyruvate metabolism 20 map00620 泛素介导的蛋白水解作用 Ubiquitin mediated proteolysis 20 map04120 萜类化合物生物合成支柱 sesquiterpenoid and triterpenoid biosynthesis 20 map00900 光合作用 Photosynthesis 19 map00195 半乳糖代谢 Galactose metabolism 18 map00052 -
通过对健康与受害白蜡树韧皮部进行转录因子家族分析,共发现了20个转录因子家族,分别为MYB超级家族、bZIP、C2H2、C2C2、AP2/ERF、bHLH、C3H、NAC、WRKY、LBD、B3超级家族、GRAS、HSF、FAR1、MADS、LOB、TCP、SBP、NF-Y和CAMTA转录因子家族(图5)。其中,C3H转录因子家族差异显著基因最多为73个,上调48个,下调25个;其次为BHLH转录因子家族,共58个基因差异显著,上调40个,下调18个;另外,NAC、MYB、B3、GRAS、SBP等转录因子家族基因表达量达到显著差异(表5)。CAMTA转录因子家族没有差异显著基因,推测此转录因子家族对白蜡窄吉丁胁迫无明显响应作用。综上所述,以上20个转录因子家族中除CAMTA转录因子家族外均可能参与白蜡树响应白蜡窄吉丁胁迫的过程,其中C3H、BHLH、NAC、MYB和B3转录因子家族极有可能在白蜡树抵御白蜡窄吉丁受害的过程中发挥重要作用,该结果为本试验后续研究奠定理论基础。
图 5 健康与受害白蜡树韧皮部转录因子家族统计
Figure 5. Family statistics of transcription factors in phloem of Fraxinus velutina
表 5 转录因子家族基因表达量差异统计
Table 5. Statistical table of gene expression difference of transcription factor family
转录因子家族名称
Transcription factor familyUp Down C3H 48 25 BHLH 40 18 NAC 25 32 MYB 28 22 B3 20 25 GRAs 16 9 SBP 15 8 WRKY 20 2 LOB 8 10 HSF 7 5 bZIP 0 11 NF-Y 1 8 TCP 1 5 C2H2 1 4 AP2/ERF 2 0 LBD 1 1 FAR1 0 2 MADS 0 2 CAMTA 0 0
绒毛白蜡(Fraxinus velutina Torr)韧皮部响应白蜡窄吉丁(Agrilus planipennis Fairmaire)危害的转录组变化
Transcriptome Changes in the Phloem of Fraxinus velutina Torr Response to Infection of Agrilus planipennis Fairmaire
-
摘要:
目的 为深入研究绒毛白蜡树响应白蜡窄吉丁危害的分子机制。 方法 以健康与受害的绒毛白蜡树树干韧皮部为试验材料,利用RNA-Seq测序技术对绒毛白蜡树韧皮部的转录组表达变化进行分析。 结果 通过对比健康与受害韧皮部的转录组数据,共鉴定出3388个DEGs,其中受害韧皮相对于健康白蜡树韧皮表达上调的DEGs有2141个,表达下调的DEGs为1247个;通过GO功能注释,将差异基因划分为20个功能类别,包括细胞过程、代谢过程、催化活性、结合元件、转运活性以及核酸结合转录因子活性等;对差异表达基因分别进行KEGG代谢途径及功能富集分析,其中代谢途径共20个,包括碳水化合物代谢,氨基酸代谢,能量代谢和脂质代谢等,另外差异表达基因分别在122条通路中均有富集,包括植物-病原体互作、代谢进程、碳水化合物结合、核酸结合转录因子活性等;通过转录组家族统计,发现健康与受害白蜡树韧皮部共有20个转录因子家族,其中C3H、BHLH、NAC、MYB、B3、GRAS和SBP等转录因子家族基因表达量均达到显著差异。 结论 研究结果为揭示白蜡树应对虫害胁迫反应的分子机制提供分子与理论依据。 Abstract:Objective To study the molecular response mechanism of Fraxinus velutina Torr to the infection of Agrilus planipennis Fairmaire. Method The phloem of non-infested and infested F. velutina stem were used as the experimental materials to analyze the transcriptome expression changes of the phloem using RNA-Seq sequencing technology. Result By comparing the transcriptome data of non-infested and infested phloem, a total of 3388 DEGs were identified, among which 2141 were up-regulated and 1247 were down regulated. Through GO functional annotation, the differential genes were divided into 20 functional categories, including cellular processes, metabolic processes, catalytic activities, binding elements, transport activities and nucleic acid binding transcription factor activities. The metabolic pathways and functional enrichment of KEGG were analyzed. There were 20 metabolic pathways, including carbohydrate metabolism, amino acid metabolism, energy metabolism and lipid metabolism. In addition, differentially expressed genes were enriched in 122 pathways, including plant-pathogen interactions, metabolic processes, carbohydrate binding, nucleic acid-binding transcription factor activities. According to transcriptome family statistics, it was found that there were 20 transcription factor families in phloem of non-infested and infested F. velutina, among which C3H, bHLH, NAC, MYB, B3, GRAS and SBP were significantly different. Conclusion The study will provide theoretical basis for revealing the molecular response mechanism of F. velutina to insect stress. -
Key words:
- Fraxinus velutina
- / Agrilus planipennis
- / phloem
- / transcriptome
-
表 1 转录组优化组装结果评估
Table 1. Evaluation of transcriptome optimized assembly results
样品名称
Sample原始序列数据
Raw reads过滤后测序数据
Clean reads碱基错误率
Error rate/%Phred
> 20的百分比/%GC百分比
GC content/%ASF-1 49722100 49349330 0.02 98.29 43.57 ASF-2 45082112 44708080 0.02 98.3 43.4 ASF-3 45802214 45432400 0.02 98.25 43.37 Non-infested1 49287036 48778052 0.02 98.11 44.75 Non-infested2 44918932 44592026 0.02 98.34 43.58 Non-infested3 50780898 50393962 0.02 98.3 43.3 表 2 差异基因GO功能注释分类统计
Table 2. Statistical table of functional annotation classification of differential gene GO
GO 术语
GO Term差异基因表达数
Seq Number信号 signaling 97 细胞成分组织或生物合成
Cellular component organization or biogenesis184 刺激响应 Response to stimulus 214 本地化 Localization 220 生物过程调节 Regulation of biological process 315 生物调节 Biological regulation 328 单个有机体过程 Single-organism process 751 代谢过程 Metabolic process 1 202 细胞过程 Cellular process 1 217 高分子络合物 Macromolecular complex 231 细胞器部件 Organelle part 245 细胞器 Organelle 741 单元部件 Cell part 1 031 细胞 Cell 1 058 膜部件 Membrane part 1 225 膜 Membrane 1 274 核酸结合转录因子活性
Nucleic acid binding transcription factor activity143 转运活性 Transporter activity 247 结合元件 Binding 1 953 催化活性 Catalytic activity 1 963 表 3 差异基因KEGG代谢途径分类统计
Table 3. Classification and statistics of metabolic pathway of differential gene KEGG
KEGG术语
KEGG Term差异基因表达数
Seq Number核苷酸代谢 Nucleotide metabolism 1 655 其他次生代谢物的生物合成
Biosynthesis of other secondary metabolites985 糖的生物合成和代谢 Glycan biosynthesis and metabolism 877 萜类化合物和聚酮类化合物的代谢
Metabolism of terpenoids anf polyketides1 052 脂质代谢 Lipid metabolism 2 594 能量代谢 Energy metabolism 2 880 其他氨基酸的代谢 Metabolite of other amino acids 1 716 代谢的辅助因子和维生素
Metabolite of cofactors and vitamins1 610 氨基酸代谢 Amino acid metabolism 3 632 碳水化合物代谢 Carbohydrate metabolism 6 397 转录 Transcription 2 465 翻译 translation 7 022 复制和修复 Replication and repair 884 折叠、分类和降解 Folding,sorting and degradation 5 179 信号转导 Signal transduction 2 093 膜运输 Membrane transport 525 运输和分解代谢 Transport and catabolism 3 561 环境适应 Environmental adaptation 1 476 耐药性:抗菌 Drug resistance:antimicrobial 8 内分泌及代谢性疾病 Endocrine and metabolic diseases 209 表 4 绒毛白蜡受白蜡窄吉丁危害下差异基因富集程度排名前20的pathway条目
Table 4. The top 20 pathway entries of differential gene enrichment under pest stress of ash tree
KEGG术语
KEGG Term差异基因表达数
Seq NumberKEGG编号
KEGG ID内质网蛋白加工 Protein processing in endoplasmic reticulum 78 map04141 内吞作用 Endocytosis 75 map04144 植物激素信号转导 Plant hormone signal transduction 70 map04075 植物-病原互作 Plant-pathogen interaction 46 map04626 MAPK信号通路-植物 MAPK signaling transduction 44 map04016 RNA运输 RNA transport 39 map03013 氨基糖和核苷酸糖代谢 Amino sugar and nucleotide sugar metabolism 31 map00520 淀粉和蔗糖代谢 Starch and sucrose metabolism 29 map00500 嘌呤代谢 Purine metabolism 28 map00230 mRNA监测通路 MRNA monitoring pathway 26 map03015 RNA降解 RNA degradation 25 map03018 苯丙素的生物合成 Biosynthesis of phenylpropanoid 23 map00940 卟啉和叶绿素代谢 Porphyrin and chlorophyll metabolism 21 map00860 果糖和甘露糖代谢 Fructose and mannose metabolism 21 map00051 氧化磷酸化 Oxidative phosphorylation 21 map00190 丙酮酸代谢 Pyruvate metabolism 20 map00620 泛素介导的蛋白水解作用 Ubiquitin mediated proteolysis 20 map04120 萜类化合物生物合成支柱 sesquiterpenoid and triterpenoid biosynthesis 20 map00900 光合作用 Photosynthesis 19 map00195 半乳糖代谢 Galactose metabolism 18 map00052 表 5 转录因子家族基因表达量差异统计
Table 5. Statistical table of gene expression difference of transcription factor family
转录因子家族名称
Transcription factor familyUp Down C3H 48 25 BHLH 40 18 NAC 25 32 MYB 28 22 B3 20 25 GRAs 16 9 SBP 15 8 WRKY 20 2 LOB 8 10 HSF 7 5 bZIP 0 11 NF-Y 1 8 TCP 1 5 C2H2 1 4 AP2/ERF 2 0 LBD 1 1 FAR1 0 2 MADS 0 2 CAMTA 0 0 -
[1] Burr S J, Mccullough D G, Poland T M. Density of emerald ash borer (Coleoptera: Buprestidae) adults and larvae at three stages of the invasion wave[J]. Environmental Entomology, 2018, 47(1): 121-132. doi: 10.1093/ee/nvx200 [2] Davis J C, Shannon J P, Bolton N W, et al. Vegetation responses to simulated emerald ash borer infestation in Fraxinus nigra dominated wetlands of Upper Michigan, USA[J]. Canadian Journal of Forest Research, 2017, 47(3): 319-330. doi: 10.1139/cjfr-2016-0105 [3] 王小艺, 杨忠岐, 魏 可. 白蜡窄吉丁(鞘翅目: 吉丁甲科)的生物防治研究进展[J]. 中国生物防治学报, 2015, 31(5):666-678. [4] Limback C K. Tree vigor and its relation to emerald ash borer (Agrilus planipennis Fairmaire) adult host preference and larval development on green and white ash trees[J]. Dissertations & Theses - Gradworks, 2010: 1487177. [5] 孙 元, 王皙玮, 段建军, 等. 白蜡窄吉丁综合防控研究进展[J]. 中国农学通报, 2012, 29(18):158-161. [6] Hai Z T, Tong G R, Liu H, et al. Host range of emerald ash borer, Agrilus planipennis Fairmaire, its damage and the countermeasures[J]. Acta Entomologica Sinica, 2005, 48: 594-599. [7] Togola A, Boukar O, Belko N, et al. Host plant resistance to insect pests of cowpea (Vigna unguiculata L. Walp.): achievements and future prospects[J]. Euphytica, 2017, 213(11): 239. doi: 10.1007/s10681-017-2030-1 [8] 窦京海. 白蜡树害虫调查与几种新白蜡害虫的发生规律及防治[D]. 泰安, 山东农业大学, 2016. [9] Chamorro M L, Volkovitsh M G, Poland T M, et al. Preimaginal stages of the emerald ash borer, Agrilus planipennis Fairmaire (Coleoptera: Buprestidae): an invasive pest on ash trees (Fraxinus)[J]. Plos One, 2012, 7(3): e33185. doi: 10.1371/journal.pone.0033185 [10] Orlova-Bienkowskaja M J, Belokobylskij S A. Discovery of the first European parasitoid of the emerald ash borer Agrilus planipennis Fairmaire (Coleoptera: Buprestidae)[J]. European Journal of Entomology, 2014, 111(4): 594-596. doi: 10.14411/eje.2014.061 [11] 路纪芳, 王小艺, 杨忠岐. 中国白蜡窄吉丁研究进展[J]. 应用昆虫学报, 2012, 49(3):785-792. doi: 10.7679/j.issn.2095-1353.2012.118 [12] Zogli P, Libault M. Plant response to biotic stress: is there a common epigenetic response during plant-pathogenic and symbiotic interactions?[J]. Plant Science, 2017, 263: 89-93. doi: 10.1016/j.plantsci.2017.07.008 [13] Nalam V J. 9-lipoxygenase oxylipin pathway in plant response to biotic stress[J]. Dissertations & Theses-Gradworks, 2012: 3533629. [14] Melvin P, Bankapalli K, D’Silva P, et al. Methylglyoxal detoxification by a DJ-1 family protein provides dual abiotic and biotic stress tolerance in transgenic plants[J]. Plant Molecular Biology, 2017, 94: 381-397. doi: 10.1007/s11103-017-0613-9 [15] Misra P, Pandey A, Tiwari M, et al. Modulation of transcriptome and metabolome of tobacco by arabidopsis transcription factor, AtMYB12, leads to insect resistance[J]. Plant Physiology, 2010, 152(4): 2258-2268.