花椒窄吉丁触角转录组及嗅觉相关基因分析

摘要

[目的]建立花椒窄吉丁Agrilus zanthoxylumi成虫触角转录组数据库,挖掘嗅觉相关基因,为今后研究其触角的化学感受机制及生物防控提供理论支撑.[方法]采用高通量测序平台Illumina NovaSeq 6000对花椒窄吉丁雌雄成虫触角进行转录组测序,用Trinity软件对获得的高质量reads进行序列拼接与组装;使用BLAST软件将触角转录组数据比对NR,NT,Swiss-Prot,GO,KEGG,BLASTX,eggNOG,Pfam,TmHMM,SignalP,KO,Map,BLASTP和RNAMMER 公共数据库;基于初步筛选到的花椒窄吉丁候选气味结合蛋白(odorant binding protein,OBP)和化学感受蛋白(chemosensory protein,CSP)以及其他鞘翅目昆虫的同源蛋白的核苷酸序列,利用MEGA软件进行系统进化分析.运用RPKM(reads per kilobase per million mapped reads)值对嗅觉相关基因表达量进行分析.[结果]花椒窄吉丁雌雄成虫触角转录组测序共获得36 209条基因和90 982条转录本,其N50分别为2 103和2 523 bp,组装完整性较高.注释到NR数据库的基因最多(41.62％),其中与赤拟谷盗Tribolium castaneum相似基因所占比例最高(19％).在GO数据库中比对到11 614个基因,按功能分为细胞组分、分子功能与生物学进程三大类57个分支,其中分子功能大类中的结合(70.57％)与催化活性(45.51％)相关基因占比最多;KEGG代谢途径分析表明,7 427条基因参与了5类代谢通路,其中涉及信号转导的基因最多,为815条;筛选到7个候选OBP基因和5个具有全长开放阅读框的CSP基因,其编码蛋白均具有化学感受蛋白家族的典型特征.系统进化分析表明,花椒窄吉丁 OBPs和CSPs分别与苹果小吉丁A.mali的OBPs和CSPs氨基酸序列一致性最高.RPKM值表明,嗅觉相关基因AzanOBPl和AzanOBP2在雌成虫触角中不表达,在雄成虫触角中微量表达;AzanOBP3在雄成虫触角中高丰度表达.[结论]首次获得了花椒窄吉丁成虫触角转录组数据,筛选到了花椒窄吉丁OBP,CSP,Or,IR和SNMP等的嗅觉相关基因.推测触角中高丰度表达的OBPs对雄成虫识别同类异性释放的信息素或寄主植物释放的挥发物起关键作用.研究结果可为花椒窄吉丁化学感受基因功能分析及嗅觉感受机制研究奠定分子基础.