基于转录组的葡萄花翅小卷蛾三类解毒酶基因家族分析

【目的】葡萄花翅小卷蛾是我国重要的检疫性害虫,一旦入侵我国,将会对我国葡萄产业和林果业生产造成严重损失,国外主要使用化学农药防治该虫。开展葡萄花翅小卷蛾转录组测序及体内细胞色素P450单加氧酶(cytochrome P450 monooxygenase,CYP)、羧酸酯酶(carboxylesterase,CarE)和谷胱甘肽S-转移酶(glutathione S-transferase,GST)三大解毒系的表达谱分析,可为葡萄花翅小卷蛾抗药性监测与抗药性治理研究提供依据。【方法】本研究利用BGISEQ-500平台对葡萄花翅小卷蛾进行转录组测序,通过从头组装获得unigene。【结果】经数据组装并去冗余后共获得44360条unigene;通过与多个公共数据库进行同源比对,共注释了28065条unigene,通过Nr数据库注释的unigene数量最多(26388条,59.49%),显示与棉铃虫的unigene同源性最高(24.11%)。基于基因序列分析表明,葡萄花翅小卷蛾有91条P450基因、31条CarE基因和22条GST基因。系统发育分析发现,葡萄花翅小卷蛾的P450基因主要集中在clade 3和clade 42个分支,与鳞翅目近源种的CYP基因聚集在一起;CarE主要包括外源化合物解毒相关的酯酶及脂质转运和代谢酯酶;GST基因聚类主要包含Delta、Omega、Epsilon、Theta和Microsomal等亚家族。【结论】本研究获得葡萄花翅小卷蛾转录组信息,鉴定了与葡萄花翅小卷蛾代谢抗性相关的基因,可为葡萄花翅小卷蛾杀虫剂和寄主植物的适应性机制研究提供数据和参考。     

草地贪夜蛾雄性成虫和5龄幼虫的转录组比较分析

【目的】草地贪夜蛾Spodoptera frugiperda是一种新近入侵我国的重要害虫。本研究旨在对草地贪夜蛾雄性成虫和5龄幼虫两个不同发育阶段的转录组进行比较分析。【方法】利用高通量测序技术对草地贪夜蛾雄性成虫和5龄幼虫进行转录组测序和数据组装,并对转录组数据进行功能注释和比较分析。【结果】经de novo组装共获得209 002条转录本,平均长度为687.55 bp,N50为982 bp。共有46 198条(57.43%) unigene在至少一个数据库中获得功能注释,其中1 713条(2.13%) unigene在所有数据库中均能获得注释。在GO数据库中获得205 269条unigene的注释,主要包括68个功能分类;在KEGG数据库中共有3 408条unigene得到注释,涉及277个代谢通路。共鉴定到424个嗅觉相关的基因,并且在雄性成虫和5龄幼虫之间的表达存在差异。通过比较转录组分析,在雄性成虫中鉴定到9 162个上调和6 399个下调差异表达基因(DEGs);功能富集分析发现在上调DEGs中涉及信息素以及信号转导的代谢通路显著富集,而下调DEGs中涉及解毒相关的通路显著富集。【结论】这些转录组数据为探究草地贪夜蛾的生长发育、嗅觉相关功能基因以及候选分子靶标提供资源信息。     

梨小食心虫幼虫中肠中高表达消化酶和解毒酶基因的筛选

【目的】挖掘梨小食心虫Grapholita molesta幼虫中肠中高表达消化酶和解毒酶基因,为今后研究以肠道为靶标的新型农药和转基因作物提供理论依据。【方法】基于梨小食心虫4龄幼虫中肠转录组高通量测序数据的FPKM值,筛选高表达基因,进行GO功能注释和KEGG通路富集分析,并使用BLAST软件进行比对筛选高表达的消化酶和解毒酶基因,利用MEGA对这些高表达的消化酶和解毒酶及其他鳞翅目昆虫的同源蛋白进行系统发育分析。利用qRT-PCR技术对梨小食心虫幼虫不同龄期中肠中的高表达代表性消化酶和解毒酶基因表达量进行定量分析和验证。【结果】在GO数据库中注释了103 677个在梨小食心虫4龄幼虫中肠中高表达基因,包括细胞组分、分子功能和生物学进程三大类功能共41个分支。KEGG通路分析表明,10 846个高表达基因参与了5类生化代谢通路。筛选到具有完整开放阅读框的消化酶基因17个［5个胰蛋白酶(trypsin, TRY)基因、3个氨肽酶(aminopeptidase, APN)基因和9个羧肽酶(carboxypeptidase, CP)基因］和解毒酶基因32个［11个谷胱甘肽S-转移酶(glutathione S-transferase, GST)基因、13个细胞色素P450(cytochrome P450, CYP450)基因和8个羧酸脂酶(carboxylesterase, CarE)基因］。系统发育分析结果表明,梨小食心虫的消化酶同源聚类分支较为分散,GSTs和CYP450s分支聚类较为集中,但都至少与1个鳞翅目昆虫同源蛋白聚在一支。qRT-PCR验证结果表明,消化酶和解毒酶基因在不同龄期梨小食心虫幼虫中肠中的表达量差异显著,表达量均在4龄幼虫期最高。【结论】本研究成功筛选和验证部分梨小食心虫幼虫中肠中高表达的消化酶和解毒酶基因,明确其与鳞翅目其他昆虫同源蛋白的进化关系。研究结果为鳞翅目其他近缘昆虫的转录组分析和以肠道为靶标的害虫防治提供了参考。     

烟粉虱不同发育阶段解毒代谢酶基因的特异性表达

【目的】基于烟粉虱Bemisia tabaci转录组数据,系统分析了烟粉虱解毒代谢酶基因在噻虫嗪抗性品系中的表达模式,探讨了这些基因在烟粉虱不同发育阶段差异表达的生物学意义。【方法】分别收集室内长期饲养的烟粉虱噻虫嗪抗性和敏感品系的卵、4龄若虫和刚羽化1 d的雌成虫,在烟粉虱转录组数据库中挑选8 394条解毒代谢相关基因设计探针,通过探针杂交,得到烟粉虱噻虫嗪抗性品系表达谱芯片,比较了这些基因在抗性烟粉虱3个不同发育阶段的表达情况。并随机挑选了9个基因,在抗感品系间3个不同发育阶段进行了荧光定量PCR验证。【结果】在抗性烟粉虱的卵和4龄若虫发育阶段,共有3 424个差异表达基因,其中489个是编码3类解毒代谢酶(羧酸酯酶、谷胱甘肽-S-转移酶和细胞色素P450多功能氧化酶)的基因;有14个基因在4龄若虫发育阶段过量表达,其中10个为P450基因,4个为GST家族基因。在抗性烟粉虱的4龄若虫和雌成虫发育阶段,总共有1 273个差异表达基因,193个为3类解毒代谢酶家族的基因,其中有9个P450家族基因在雌成虫期的表达量超过4龄若虫期的10倍。此外,表达谱芯片分析还筛选到了一些候选抗性基因。qRT-PCR验证显示在这些候选基因中,与敏感品系相比,9个基因在抗性烟粉虱的3个不同发育阶段表达上调,其中GST基因家族的p_06027和P450基因p_06013在抗性品系的卵和4龄若虫中过量表达;p_05885和p_07806和编码CYP6家族蛋白的p_00988在抗性品系的4龄若虫期的表达量上调;p_05916和p_00478在抗性品系卵和4龄若虫期表达量很低,而在成虫期过量表达;而p_00059和p_00428在抗性品系雌成虫发育阶段表达量显著上调,其中编码CYP4C1的p_00059的差异表达倍数在雌成虫期约为15.15倍。【结论】表达谱芯片分析结果提示,CYP6和CYP4C1基因的过量表达可能会是烟粉虱抗性产生的机制之一。解毒代谢酶基因在烟粉虱不同发育阶段的特异性表达,可能与其在抵御杀虫剂胁迫时体内能量的分布及有效利用率有关,也可能是害虫在环境选择压下的一种适应机制。     

灰飞虱唾液腺三大解毒酶家族的转录组分析

灰飞虱Laodelphax striatellus （Fallén）是危害多种禾本科经济作物的重要刺吸式害虫。唾液腺对刺吸式口器昆虫取食植物尤其重要, 其分泌的唾液可以帮助刺吸式口器昆虫刺穿植物、 消化食物、 解毒植物的次生物质。细胞色素P450单加氧酶（cytochrome P450 monooxygenase, P450）、 谷胱甘肽S 转移酶（glutathione S transferase, GST）和羧酸酯酶（carboxylesterase, CarE）是昆虫主要的解毒酶系。为了分析解毒酶基因在灰飞虱唾液腺中的表达谱, 本研究对灰飞虱成虫唾液腺进行转录组测序、 重头组装和注释, 并与豌豆蚜Acyrthosiphon pisum和西方蜜蜂Apis mellifera的同源蛋白进行系统发育分析。发现有9个谷胱甘肽S 转移酶（glutathione S transferase, GST）基因、 22个羧酸酯酶（carboxylesterase, CarE）基因和39个细胞色素P450单加氧酶（cytochrome P450 monooxygenase, P450）基因在灰飞虱唾液腺中表达。通过对同源蛋白进行系统发育分析, 发现灰飞虱唾液腺大部分的CarE是参与消化/解毒和激素/信息素的加工, 而参与神经/发育的CarE很少; 灰飞虱唾液腺表达的P450基因远远少于豌豆蚜和西方蜜蜂基因组的P450基因数, 且只有CYP6和CYP4家族的成员; GST家族在3种昆虫的保守性最高。研究结果为灰飞虱对寄主植物和杀虫剂的适应性研究奠定了基础。     

新疆和云南番茄潜叶蛾种群对六种杀虫剂的敏感性及其与解毒酶活性的关系

【目的】番茄潜叶蛾Tuta absoluta 是新入侵我国的对番茄具有毁灭性危害的入侵害虫,目前入侵我国的番茄潜叶蛾种群对杀虫剂的抗性尚无报道。本研究旨在明确新疆和云南番茄潜叶蛾田间种群对6种常用杀虫剂的敏感性及其与解毒酶活性的关系。【方法】采用浸叶法测定6种常用杀虫剂对番茄潜叶蛾新疆和云南种群2龄幼虫的室内毒力。通过对2龄幼虫的生物测定确定3种增效剂［CYP450抑制剂胡椒基丁醚(PBO)、酯酶抑制剂磷酸三苯酯(TPP)和GST抑制剂丁烯二酸二乙酯(DEM)］对氯虫苯甲酰胺的增效作用。采用酶活性分析测定室内敏感种群和田间抗性种群(新疆种群) 2龄幼虫体内解毒酶［细胞色素P450酶(CYP450)、谷胱甘肽S-转移酶(GST)和羧酸酯酶(CarE)］活性,以确定杀虫剂抗性与解毒酶活性的关系。【结果】番茄潜叶蛾云南种群对6种杀虫剂的敏感性由高到低依次为甲维盐、溴虫腈、多杀菌素、茚虫威、氯虫苯甲酰胺和高效氯氰菊酯。新疆种群对6种杀虫剂的敏感性由高到低依次为甲维盐、溴虫腈、氯虫苯甲酰胺、多杀菌素、茚虫威和高效氯氰菊酯。与室内敏感种群相比,云南和新疆种群对氯虫苯甲酰胺的抗性水平最高,抗性倍数分别为212.7和169.3倍。生物测定结果表明,3种增效剂PBO, TPP和DEM均对氯虫苯甲酰胺无明显增效作用。酶活性测定结果表明,番茄潜叶蛾室内敏感种群和田间抗性种群之间2龄幼虫中CYP450, GST和CarE活性无显著差异。【结论】番茄潜叶蛾新疆和云南种群对测试的6种杀虫剂产生不同程度的抗性,对氯虫苯甲酰胺的抗性最高,番茄潜叶蛾对杀虫剂的抗性与解毒酶活性无关。本研究的结果对番茄潜叶蛾的田间防治和杀虫剂抗性治理具有指导意义。     

致倦库蚊GSTe1和CYP9J40的时空表达特点

【目的】细胞色素P450单加氧酶(P450)和谷胱甘肽-S-转移酶(GST)是生物降解外源有毒物质的重要解毒酶家族,在昆虫适应杀虫剂选择的过程中发挥重要的作用。以往在致倦库蚊Culex pipiens quinquefasciatus对硫磷抗性研究中,发现了一个Epsilon GST基因(GSTe1)和一个P450 CYP9家族基因(CYP9J40)在抗性品系中高表达。研究它们在蚊虫中的时空表达模式对理解它们的生物学功能有重要意义。【方法】本研究利用实时定量PCR和原位杂交的方法探讨了这两个解毒酶基因在蚊虫中的时空表达。【结果】定量PCR显示GSTe1和CYP9J40在对硫磷抗性蚊虫中的转录水平分别是敏感蚊虫的4.2倍和1.9倍,它们在蚊虫整个生活周期的表达模式相似,均为幼虫期表达量相对较低,从蛹期开始表达量显著升高,到成虫期雌虫达到顶峰,而雄虫与幼虫期表达量相当。原位杂交结果表明这两个基因在肌肉、消化道和卵巢中都有表达,且CYP9J40的表达更广泛,在脑和胸部神经节也有明显表达。【结论】基因表达的发育阶段特异性和广泛的组织分布说明GSTe1和CYP9J40除了参与蚊虫对杀虫剂的抗性,还可能具有其它生理功能,如对雌虫发育有重要的影响。     

亚致死浓度溴虫腈和毒死蜱对等钳蠊螨生长繁殖和解毒酶的影响

【目的】害虫综合治理中,化学防治对天敌资源也具有一定的杀伤力。为缓解这一矛盾,本研究旨在探究溴虫腈和毒死蜱亚致死浓度处理对等钳蠊螨Blattisocius dentriticus生长繁殖及解毒酶活性和基因表达量的影响。【方法】利用药膜法处理2-3日龄等钳蠊螨成螨24 h,测定溴虫腈和毒死蜱的亚致死剂量;测定和比较LC_(10)和LC_(30)剂量这两种药剂处理后F_0和F_1代的产卵量、产卵期、卵孵化率和雌成螨寿命等生物学特性的变化;通过酶活力分析和RT-PCR分别测定LC_(10), LC_(30)和LC_(50)剂量下这两种药剂处理后钳蠊螨成螨体内谷胱甘肽-S-转移酶(glutathione S-transferase, GST)、细胞色素P450(cytochrome P450, CYP450)、羧酸酯酶(carboxylesterase,CarE)酶活力及其基因表达量变化。【结果】溴虫腈和毒死蜱处理24 h对等钳蠊螨成螨的LC_(50)分别为42.56 mg/L和72.42 mg/L。溴虫腈和毒死蜱LC_(10)和LC_(30)剂量处理等钳蠊螨雌成螨后,与对照(清水处理)相比,仅F_0代雌成螨寿命和产卵期显著缩短(P0.05),而产卵量和卵孵化率无明变化。酶活力测定结果发现,GST, CYP450和CarE的活力在溴虫腈和毒死蜱LC_(10)和LC_(30)剂量处理后无明显变化,而LC_(50)剂量下,上述3种酶的活力均显著增加(P0.05)。基因表达结果表明,溴虫腈LC_(10), LC_(30)和LC_(50)剂量处理下2个GST基因(BdGST3和BdGST6)、3个CYP450基因(BdCYP2-4)和5个CarE基因(BdCarE1-5)表达均显著上调;在毒死蜱这3个剂量处理下3个GST基因(BdGST1,BdGST3和BdGST4)、3个CYP450基因(BdCYP2,BdCYP5和BdCYP6)和2个CarE基因(BdCarE1和BdCarE2)表达量均显著上调。【结论】结果表明,LC_(10)和LC_(30)剂量的溴虫腈和毒死蜱亚致死剂量会抑制F_0代雌成螨的生长繁殖;LC_(10), LC_(30)和LC_(50)剂量下这两种药剂可诱导等钳蠊螨GST, CarE和CYP450基因表达;LC_(50)剂量能明显诱导等钳蠊螨体内GST, CarE和CYP450活性上升。该研究为等钳蠊螨抗性品系的筛选及田间应用提供了理论依据。     

5种杀虫剂对西花蓟马和花蓟马的毒力及其生理酶活性的影响

【目的】探究不同杀虫剂对重要入侵害虫西花蓟马及其本地近缘种花蓟马的毒力及对保护酶和解毒酶活性的影响,为进一步研究2种害虫的抗性管理提供依据。【方法】采用浸渍法测定5种田间常用杀虫剂对西花蓟马和花蓟马的毒力,并测定杀虫剂亚致死浓度(LC25)下2种蓟马体内保护酶和解毒酶活性的差异。【结果】不同杀虫剂对2种蓟马的毒力依次为:乙基多杀菌素甲维盐阿维菌素吡虫啉噻虫嗪,乙基多杀菌素对西花蓟马和花蓟马的LC_(50)分别为0.28和0.03 mg·L~(-1)。不同药剂的亚致死剂量(LC_(25))对西花蓟马和花蓟马体内保护酶和解毒酶活性普遍具有诱导作用。其中,阿维菌素对西花蓟马超氧化物歧化酶(SOD)活性诱导作用最强,为326.40 U·mg~(-1),是对照的9.37倍,而乙基多杀菌素对花蓟马SOD活性诱导作用最强,为245.35 U·mg~(-1),是对照的9.32倍;吡虫啉对西花蓟马和花蓟马过氧化物酶(POD)诱导作用最强,分别为298.67和246.79 U·mg~(-1),是对照的37.10和20.57倍;阿维菌素对西花蓟马和花蓟马过氧化氢酶(CAT)和羧酸酯酶(CarE)诱导作用最强,分别为298.67、246.79 U·mg~(-1)(CAT活性)和12.53、11.99 U·mg~(-1)(CarE活性);乙基多杀菌素对西花蓟马和花蓟马谷胱甘肽转移酶(GST)和乙酰胆碱酯酶(AChE)诱导作用最强,分别为77527.59、66927.39 U·mg~(-1)(GST活性)和2.34、2.22 U·mg~(-1)(AChE活性)。【结论】5种杀虫剂中,乙基多杀菌素对2种蓟马的毒力最强;西花蓟马对杀虫剂的解毒代谢能力强于花蓟马。     

七星瓢虫触角转录组及嗅觉相关基因分析

【目的】瓢虫科食性高度分化。本研究旨在通过建立七星瓢虫Coccinella septempunctata触角转录组数据库,探讨其触角嗅觉相关基因与食性分化的关系。【方法】采用Illumina HiSeq 4000高通量测序平台对七星瓢虫成虫触角转录组进行测序、组装、注释,挖掘嗅觉相关基因,并与已发表的茄二十八星瓢虫Henosepilachna viginyioctopunctata转录组进行比较。【结果】共获得七星瓢虫触角转录组31 775条unigenes,其中69.71%的序列得以注释,NR数据库中注释最多,为20 539条。据注释信息,挖掘到27个嗅觉相关基因,包括1个气味结合蛋白(odorant binding protein, OBP)基因,13个化学感受蛋白(chemosensory protein, CSP)基因,4个气味受体(odorant receptor, OR)基因,7个味觉受体(gustatory receptor, GR)基因和2个感觉神经元膜蛋白(sensory neuron membrane protein, SNMP)基因。相对应地,在植食性茄二十八星瓢虫转录组中鉴定到38个嗅觉相关基因,包括七星瓢虫中未发现的1个离子型受体(ionotropic receptor, IR)基因。在各类型嗅觉相关基因中,茄二十八星瓢虫转录组的OBP基因比例(13.16%)高于七星瓢虫触角转录组的(3.70%),而七星瓢虫触角转录组的GR基因比例(25.93%)则高于茄二十八星瓢虫转录组的(13.16%)。【结论】触角嗅觉相关基因数目不是昆虫食性分化的主因。本研究获得了七星瓢虫触角转录组学资源,初步探讨了嗅觉相关基因同瓢虫食性分化的关系,为了解瓢虫乃至昆虫食性分化的分子基础提供了信息。     

松树蜂毒腺转录组分析

【目的】构建入侵种松树蜂Sirex noctilio毒腺转录组数据库,筛选并分析松树蜂毒腺基因数据。【方法】采用新一代高通量测序平台Illumina HiSeq^TM 4000对松树蜂雌成虫毒腺进行转录组测序、数据组装和生物信息学分析。【结果】共获得12.7 Gb松树蜂雌成虫毒腺有效转录组数据,并组装到37 098条unigenes,平均长度968 bp,N50长度为2 364 bp。将所得的unigenes数据使用BlastX与各大数据库比对,共注释到13 515条unigenes,并且在NR数据库中注释的unigenes最多,共11 108条(占总数的29.94%),其中相似基因占比最高的物种为丽蝇蛹集金小蜂Nasonia vitripennis,达815条(占总数的7.29%)。在GO数据库中注释到5 726条unigenes,根据功能被分为生物学进程、细胞组分和分子功能3大类63个亚类。KEGG代谢通路分析表明,7 602条unigenes注释到357个代谢通路。根据基因注释信息进一步筛选到43条嗅觉相关基因,包括嗅觉受体(odorant receptor, Or)基因25条、化学感受蛋白(chemosensory protein, CSP)基因10条、离子型受体(ionotropic receptor, IR)基因5条和气味结合蛋白(odorant binding protein, OBP)基因3条。此外,还筛选出11条漆酶基因,包括漆酶1(laccase1, LAC1)基因5条、漆酶2(laccase2, LAC2)基因4条、漆酶4(laccase4, LAC4)基因1条和漆酶9(laccase9, LAC9)基因1条,且其中1条LAC2基因在所有被注释的基因中表达量最高(FPKM值=21 126)。【结论】本研究获得的松树蜂毒腺转录组数据为松树蜂毒液组分的鉴定和生物学功能的研究奠定了一定的理论基础。     

大垫尖翅蝗转录组分析

【目的】大垫尖翅蝗Epacromius coerulipes(Ivanov)是广泛分布的草原蝗虫之一,但其基因资源缺乏。为了获得大垫尖翅蝗的基因数据,对其进行了转录组测序和分析。【方法】利用Illumina公司paired-end转录组的测序技术进行从头组装。【结果】总计获得了63 033条unigenes,平均长度为772 bp,N50为1 589 bp。通过BLAST搜索,确定有25 132条(39.87%)unigenes与NCBI数据库已知的蛋白质相匹配,其中有24 841,16 490,11 558和8 013条unigenes成功注释到Nr,Swiss-Prot,GO和COG数据库中。KEGG数据库中,7 218条unigenes形成218条代谢或信息通路。其中,189条unigenes参与外源性物质或药物的代谢通路。进一步分析显示,213条unigenes被确认为可能参与外源性物质的解毒作用,29条unigenes被确定为编码杀虫剂的目标蛋白。此外,检测到5 696条简单重复序列。【结论】该转录组测序分析将为进一步研究大垫尖翅蝗的基因功能分析及杀虫剂的抗药性机制分析奠定分子基础。     

我国棉花主产区棉蚜对吡虫啉的抗性监测及抗性机理

【目的】由于生长周期短、繁殖率高,棉蚜Aphis gossypii容易对杀虫剂产生抗药性。本研究旨在明确我国棉花主产区棉蚜对吡虫啉的抗性水平及抗性机理。【方法】采用浸叶法测定了北京海淀,河北廊坊和邯郸,山东德州,河南许昌,以及新疆奎屯和阿克苏地区棉蚜对吡虫啉的抗性水平;测定了不同种群棉蚜3种解毒酶（多功能氧化酶、羧酸酯酶、谷胱甘肽S-转移酶）及乙酰胆碱酯酶的活性;并对靶标基因烟碱型乙酰胆碱受体（nAChR）β1亚基基因进行了突变检测。【结果】北京海淀、河南许昌和河北邯郸的棉蚜对吡虫啉敏感;河北廊坊、新疆阿克苏、山东德州及新疆奎屯地区的棉蚜对吡虫啉的抗性倍数(resistance ratio, RR)分别为22.6, 26.3,53.5和61.1倍,为中等水平抗性。酶活力对比研究发现,阿克苏和奎屯地区的棉蚜多功能氧化酶的比活力分别是敏感种群(北京种群)的1.7和1.8倍,羧酸酯酶的比活力分别是敏感种群的1.6和1.7倍,谷胱甘肽S-转移酶的比活力均是敏感种群的1.5倍,但是乙酰胆碱酯酶比活力在棉蚜种群间差异不显著。靶标基因突变检测表明,河北廊坊、新疆阿克苏、山东德州及新疆奎屯棉蚜种群nAChR β1亚基均存在与吡虫啉抗性相关的精氨酸到苏氨酸（R81T）突变。【结论】结果提示,多功能氧化酶、羧酸酯酶和谷胱甘肽S-转移酶活力升高以及nAChR β1亚基R81T突变与棉蚜对吡虫啉的抗性形成相关。     

花椒窄吉丁转录组及化学感受相关基因的分析

【目的】建立花椒窄吉丁 Agrilus zanthoxylumi 转录组数据库,挖掘其基因组数据。【方法】采用高通量测序平台Illumina NovaSeq6000对花椒窄吉丁成虫不同组织(触角、头、胸、腹、足、翅)进行转录组测序、序列组装,使用BLAST软件将花椒窄吉丁unigenes序列与公共数据库进行比对,BLAST同源性搜索的方法从中筛选出与其化学感受相关的基因,并与其他已研究发表的鞘翅目昆虫化学感受相关基因的核酸序列比对,进行系统发育分析。利用RPKM值对这些基因在花椒窄吉丁成虫不同组织中的表达量进行分析。【结果】经测序及序列拼接后共得到80 320条unigenes和 169 398 条contigs,其G+C比例分别是37.63%和39.18%,平均长度分别为828.75和1 084.33 bp;unigenes长度主要分布在200~400 bp的共有37 374条,占全部unigenes的46.53%。在NR数据库中注释的与花椒窄吉丁转录组相似基因所属物种分布为赤拟谷盗 Tribolium castaneum 所占比例26.80%,其后依次是小家鼠 Mus musculus (13.87%),蛀犀金龟 Oryctes borbonicus (5.26%),山松甲虫 Dendroctonus ponderosae (4.10%),黑蚁 Lasius niger (2.36%),其他物种所占比例是47.61%。在GO数据库中比对到27 488个unigenes,可分为分子功能、细胞组分和生物学进程三大类共59个过程;利用KOG数据库功能注释分为25类,注释到unigenes 最多功能类别的是普通功能,共 4 666个,最少的是与细胞移动相关,共40个;将80 320个unigenes映射到KEGG数据库中,21 104条unigens注释到代谢通路35个,占26.27%,注释最多的代谢途径是转录途径,共涉及到2 037条unigenes。从NR数据库中共注释到8个气味受体(odorant receptor, OR)基因、7个离子型受体(ionotropic receptor, IR)基因和6个味觉受体(gustatory receptor, GR)基因,各组织中的基因表达量分析发现,这3类化学感受相关基因在成虫触角中的表达量显著高于其他组织中的表达量,雌、雄触角中表达量差异最显著的为 AzanOR 1,而在雌、雄成虫足中这3类化学感受相关基因表达量没有差异。【结论】本研究首次获得了花椒窄吉丁转录组数据,为进一步研究花椒窄吉丁的基因功能奠定了分子基础。