云斑天牛成虫触角转录组及嗅觉相关基因分析

Analysis of antennal transcriptome and olfaction-related genes of adult Batocera horsfieldi (Hope)

【目的】 建立云斑天牛Batocera horsfieldi (Hope)成虫触角转录组数据库,深入挖掘云斑天牛的基因数据信息。【方法】 采用高通量测序平台(Illumina HiSeq)对云斑天牛成虫触角进行转录组测序、序列组装及生物信息学分析。【结果】 云斑天牛成虫触角转录组共获得137 485条Transcript序列和69 214条Unigene序列;其中,Transcript序列平均长度1 142 bp,Unigene序列平均长度1 983 bp。将Unigene分别比对到NR、NT、SwissProt、KO、PFAM、GO、KOG数据库进行基因功能注释,NR注释41 636条,NT注释14 895条,KO注释19 287条,SwissProt注释33 442条,PFAM注释34 687条,GO注释35 321条,KOG注释20 582条。NR注释表明,72.0%的云斑天牛Unigene与赤拟谷盗Tribolium castaneum和中欧山松大小蠹Dendroctonus ponderosae具有相似性。基因功能注释分类表明,云斑天牛成虫触角转录组在GO数据库三大类中包含5个最主要功能,分别是细胞过程、代谢过程、单有机体过程、结合和催化活性,分别占20 912、19 086、17 202、21 477和15 823条Unigenes;云斑天牛转录组在KOG数据库26个功能目录共注释20 585条Unigenes,其中,翻译后修饰、蛋白质转换和伴侣共有1 977条,一般功能预测3 285条(最多),信号传导机制3 053条,合计8 315条占全部Unigenes 40.39%;总共19 287条Unigenes分至5个KEGG功能类别,其中细胞过程6 793条,环境信息处理6 255条,遗传信息处理3 038条,代谢3 852条,有机系统3 508条。进一步基因功能注释分析筛选得到161个嗅觉相关基因,包含96个气味结合蛋白(Odorant binding protein,OBP),34个化学感受蛋白(Chemosensory protein,CSP)和31个气味受体(Odorant receptor,OR)。【结论】 本研究获得了云斑天牛成虫触角转录组数据库,为进一步研究云斑天牛的基因功能分析及嗅觉感受机制奠定了分子基础。

[Objectives] To obtain genetic data on, and establish an antenna transcriptome database for, the adult white-striped longhorn beetle, Batocera horsfieldi (Hope). [Methods] The antenna transcriptome of adult B. horsfieldi was sequenced using an Illumina HiSeq 4000 platform and analyzed with bioinformatic software. [Results] The antennal transcriptome database of adult B. horsfieldi contained 137 485 transcripts and 69 214 Unigenes. The mean number of transcripts was 1 142 bp, and the mean number of Unigenes was 1 983 bp. Of the Unigenes found, 41 636 were annotated in the NR database, 14 895 were annotated in the NT database, 19 287 were annotated in the KO database, 33 442 were annotated in the SwissProt database, 34 687 were annotated in the PFAM database, 35 321 were annotated in the GO database, and 20 582 were annotated in the KOG database. The results of NR comment showed that 72.0% of B. horsfieldi Unigenes were similar to those of Tribolium castaneum and Dendroctonus ponderosae. Transcriptome gene function classification indicated five main functional classifications annotated in the GO database;these were cellular process, metabolic process, single-organism process, binding and catalytic activity, each classification containing 20 912, 19 086, 17 202, 21 477 and 15 823 unigenes, respectively. 20 585 unigenes were annotated in 26 groups of the KOG database, 1 977 of which were involved in post-translational modification, protein turnover and chaperoning. Only 3 285 unigenes were involved in general function prediction, and 3 053 in signal transduction mechanisms, together comprising 40.39% of the total. 19 287 unigenes were annotated in 5 functional classifications in the KEGG database;cellular process (6 793), environmental information processing (6 255), genetic information processing (3 038), metabolism (3 852) and organismal system (3 508). By further screening and identification, 161 olfaction-related genes, including 96 odorant binding protein (OBP) genes, 34 chemosensory protein (CSP) genes and 31 odorant receptor (OR) genes, were identified. [Conclusion] The transcriptome database of the antenna of adult B. horsfieldi was successfully obtained. The results provide a foundation for the further study of gene function and the molecular mechanisms responsible for olfactory reception in B. horsfieldi.