摘要

【目的】通过对大豆受豆卷叶螟幼虫胁迫下的转录组和蛋白质组结果进行联合分析,筛选出一些与大豆抗豆卷叶螟相关的候选基因,为深入认识大豆抗豆卷叶螟的分子调控机制奠定基础。【方法】以高抗材料赶泰-2-2(HR)和高感材料皖82-178(HS)为研究对象,运用RNA-Seq技术和i TRAQ技术鉴定出豆卷叶螟幼虫取食诱导0和48 h时样品间的差异表达基因(DEGs)和差异表达蛋白(DEPs),将在蛋白水平和转录水平关联到的所有可定量数据进行关联分析,计算蛋白水平和转录水平间的相关系数。【结果】蛋白质鉴定结果表明,HR48/HR0、HS48/HS0、HR0/HS0和HR48/HS48比较组中分别鉴定出236、250、213、211个DEPs;转录组鉴定结果表明,HR48/HR0、HS48/HS0、HR0/HS0和HR48/HS48比较组中分别鉴定出1 064、680、605、468个DEGs。定量关联分析结果表明,HR48/HR0、HS48/HS0、HR0/HS0和HR48/HS48比较组的相关系数r分别为0.156、0.2687、0.1149和0.035;HR48/HR0、HS48/HS0、HR0/HS0和HR48/HS48比较组中分别关联到11、9、3和4个DEPs/DEGs,其中蛋白和m RNA表达趋势相同的基因分别有11、9、2和3个,蛋白和m RNA表达趋势相反的基因分别有0、0、1和1个。生物信息学分析结果表明,这些差异蛋白功能涉及代谢途径、核糖体、类黄酮生物合成、亚油酸代谢、氨基糖-核苷酸代谢、氨酰生物合成、亚麻酸代谢、次生代谢产物的生物合成、苯基丙酸生物合成、RNA运转、谷胱甘肽代谢、抗坏血酸盐和aldarate代谢等。功能关联分析结果表明,HR48/HR0比较组中有27条Pathway通路共关联到101个DEPs和23个DEGs,HS48/HS0比较组中有15条Pathway通路共关联到147个DEPs和16个DEGs,HR0/HS0比较组中有18条Pathway通路共关联到82个DEPs和10个DEGs,HR48/HS48比较组中仅有1条Pathway通路关联到71个DEPs和2个DEGs。同时关联发现胰蛋白酶抑制剂A、K型胰蛋白酶抑制剂、查尔酮异构酶4、脂氧合酶9、α-加氧合酶1、9S-脂氧合酶、植物凝集素前体、POD12、应激蛋白SAM22和c APX1等可能是大豆抗豆卷叶螟潜在的靶标蛋白(基因)。q RT-PCR表达分析结果表明,5个DEPs/DEGs在HR48/HR0比较组的表达趋势与它们的RNA-Seq和i TRAQ结果基本一致。【结论】确定了一些与豆卷叶螟的抗性相关蛋白(基因)和代谢通路,这些差异蛋白可能直接或间接参与大豆对豆卷叶螟的抗虫反应。