松树蜂毒腺转录组分析

Analysis of the venom gland transcriptome of Sirex noctilio (Hymenoptera: Siricidae)

【目的】构建入侵种松树蜂Sirex noctilio毒腺转录组数据库,筛选并分析松树蜂毒腺基因数据。【方法】采用新一代高通量测序平台Illumina HiSeqTM 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)。【结论】本研究获得的松树蜂毒腺转录组数据为松树蜂毒液组分的鉴定和生物学功能的研究奠定了一定的理论基础。

【Aim】 To establish the venom gland transcriptome database of the wood wasp, Sirex noctilio, and to screen and analyze the genetic data of its venom gland.【Methods】 The venom gland transcriptome of S. noctilio female adults was sequenced using an Illumina HiSeqTM 4000 platform, assembled and subjected to bioinformatics analysis.【Results】 A total of 12.7 Gb clean data of the venom gland transcriptome of S. noctilio female adults was obtained, and the clean reads were assembled into 37 098 unigenes with a mean length of 968 bp and an N50 length of 2 364 bp. A total of 13 515 unigenes were annotated using BlastX searches against the databases. Most unigenes(11 108 unigenes) were annotated to the NR database(accounting for 29.94% of the total), showing the highest homology with the unigenes of Nasonia vitripennis in the NR database by 815 similar unigenes(accounting for 7.29% of the total). A total of 5 726 unigenes were annotated in the GO database and divided into three categories(biological processes, cellular component and molecular function) and 63 subclasses based on their biological functions. In KEGG database, 7 602 unigenes were assigned to the 357 known metabolic pathways. By further screening and identification, 43 olfaction-related genes of S. noctilio including 25 odorant receptor(OR) genes, 10 chemosensory protein(CSP) genes, five ionotropic receptor(IR) genes, and three odorant binding protein(OBP) genes were obtained. In addition, 11 laccase genes including five laccase-1 genes, four laccase-2 genes, one laccase-4 gene and one laccase-9 gene were screened. Among the annotated genes, a LAC2 gene had the highest expression level, with the FPKM value of 21 126.【Conclusion】 The venom gland transcriptome data of S. noctilio obtained in this study provide a theoretical foundation for the component identification and biological function study of S. noctilio venom.