绿豆象触角转录组及嗅觉相关基因的分析

【目的】建立绿豆象Callosobruchus chinensis触角转录组数据库,挖掘绿豆象的基因数据。【方法】采用高通量测序平台(Illumina Hiseq)对绿豆象成虫触角进行转录组测序、序列组装及生物信息学分析。【结果】共获得51.10 Gb的clean data(NCBI SRA数据库登录号:SRP119884)及83 535条unigenes,长度在1 kb以上的unigenes有15 075条,unigenes的N50为1 492 bp。使用BLAST软件将绿豆象触角Unigene序列与公共数据库比对,共注释到22 148条unigenes,其中NR数据库注释的unigenes最多,为18 744条,且与赤拟谷盗Tribolium castaneum的相似基因最多,达39.57%。通过GO数据库注释,将unigenes功能分为细胞组分、分子功能与生物学过程三大类54个分支,其中参与代谢进程以及结合与催化活性的unigenes比例较大。KEGG代谢途径分析表明,9 084条unigenes形成289条代谢通路,其中核糖体代谢通路所含unigene最多,为516条。进一步基因注释分析筛选到140个嗅觉相关基因,包括12个气味结合蛋白(odorant binding protein,OBP)基因,4个化学感受蛋白(chemosensory protein,CSP)基因,116个气味受体(odorant receptor,Or)基因,1个离子型受体(ionotropic receptor,IR)基因和7个感觉神经元膜蛋白(sensory neuron membrane protein,SNMP)基因,并用FPKM值对基因表达量进行评估。【结论】本研究获得了绿豆象触角转录组数据,为进一步研究绿豆象的基因功能分析及嗅觉感受机制奠定分子基础。     

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

【目的】瓢虫科食性高度分化。本研究旨在通过建立七星瓢虫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%)。【结论】触角嗅觉相关基因数目不是昆虫食性分化的主因。本研究获得了七星瓢虫触角转录组学资源,初步探讨了嗅觉相关基因同瓢虫食性分化的关系,为了解瓢虫乃至昆虫食性分化的分子基础提供了信息。     

大蜡螟触角转录组及嗅觉相关基因分析

【目的】大蜡螟Galleria mellonella Linnaeus是世界范围内蜂群中普遍存在的害虫,给世界养蜂业带来了巨大的危害。本研究旨在建立大蜡螟触角转录组数据库,挖掘大蜡螟嗅觉相关基因,为今后利用嗅觉基因靶标防治大蜡螟提供分子生物学信息数据。【方法】利用高通量测序平台(Illumina HiSeqTM 2500)对大蜡螟触角进行转录组测序和组装,通过筛选转录组数据库,鉴定出嗅觉相关基因。【结果】成功构建了大蜡螟触角转录组,经序列拼接获得188 278条unigenes,平均长度为781 bp,N50为1 161 bp。使用BLAST软件将大蜡螟触角unigenes序列与7大公共数据库比对,共注释到108 047条unigenes,其中NR数据库注释的unigenes最多,为78 746条,占41.82%,且NR注释的大蜡螟unigenes中与家蚕Bombyx mori类似基因最多,为18.4%。将unigenes与GO数据库比对发现,共有69 794条unigenes注释到GO数据库中,按照其所参与的生物过程、细胞组分和分子功能对其进行GO功能分类,共含有55个亚类,其中参与结合和催化功能及代谢过程的unigenes比例较大;KEGG代谢途径分析表明,14699条unigenes形成231条代谢通路,其中被注释在信号传导通路中的unigenes最多,为2 401条,占16.33%。进一步基因注释分析,鉴定得到226个候选嗅觉相关基因,包括52个气味结合蛋白(Odorant binding proteins,OBPs)基因,36个化学感受蛋白(Chemosensory proteins,CSPs)基因,80个气味受体(Odorant receptors,ORs)基因,56个离子型受体(Ionotropic receptors,IRs)基因和2个感觉神经元膜蛋白(Sensory neuron membrane proteins,SNMPs)基因;通过比较雌、雄大蜡螟转录组鉴定出114个差异表达基因,包括5个嗅觉相关基因(OBP8、OBP17、CSP3、CSP34和OR15);114个差异表达基因中,66个基因在雌蛾触角中高表达,48个基因在雄蛾触角中高表达。【结论】本研究获得了大蜡螟触角转录组数据,并鉴定出候选嗅觉相关基因,结果为进一步研究大蜡螟的基因功能及嗅觉感受机制奠定分子基础。     

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

【目的】建立花椒窄吉丁 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类化学感受相关基因表达量没有差异。【结论】本研究首次获得了花椒窄吉丁转录组数据,为进一步研究花椒窄吉丁的基因功能奠定了分子基础。     

黑须污蝇不同发育阶段转录组及嗅觉相关基因分析

【目的】建立黑须污蝇Wohlfahrtia magnifica转录组数据库,挖掘黑须污蝇基因数据,为更深入地研究双峰驼阴道蝇蛆病提供理论依据。【方法】采用高通量测序平台Illumina HiseqTM4000对黑须污蝇幼虫、蛹和成虫进行转录组测序,并进行生物信息学分析。利用黑须污蝇幼虫、蛹和成虫转录组数据,通过基因差异表达分析以及GO功能显著性富集分析的方法筛选出嗅觉相关基因,并通过荧光定量PCR技术对黑须污蝇幼虫、蛹和成虫中OBP99b,OBP56a和OBP99a 3个嗅觉相关基因表达水平进行验证。【结果】结果显示,每个黑须污蝇样本的转录组数据量在4.96 Gb以上,G+C含量在35.35%以上,Q20含量在97%以上。与NCBInr, GO, KEGG, Pfam, Swiss-Prot和eggNOG数据库进行对比,共注释到73 303条unigenes。其中eggNOG数据库注释到的unigenes最多,为35 066条,分布在23个蛋白功能中,其中的翻译修饰、蛋白质周转的unigenes所占比例较大;GO数据库注释到的unigenes数为29 193条,功能包括分子功能、细胞组分、生物学过程三大类50个分支,其中参与生物学过程和氧化还原过程的unigenes比例较大;KEGG数据库注释到的unigenes数为15 068条,其中参与信号转导过程的unigenes比例较大。上述数据表明该转录组测序结果较好,为后续研究奠定了基础。依据黑须污蝇幼虫、蛹、成虫转录组数据筛选到30个嗅觉相关基因,其中包含9个气味结合蛋白(OBP)基因,进一步基因注释发现OBP99b,OBP56a和OBP99a基因在黑须污蝇不同发育阶段表达差异显著(|log2FC|>1,P<0.05),这在荧光定量PCR分析结果中得到了验证。【结论】本研究获得了黑须污蝇幼虫、蛹和成虫转录组数据,筛选出黑须污蝇各发育阶段嗅觉相关基因,并验证了OBP99b,OBP56a和OBP99a基因在黑须污蝇幼虫、蛹和成虫3个发育阶段差异表达,为防治双峰驼阴道蝇蛆病提供了新思路。     

荔枝蒂蛀虫转录组及嗅觉相关基因分析

【目的】荔枝蒂蛀虫Conopomorpha sinensis Bradley是专一性危害我国华南地区荔枝和龙眼的重要害虫,隐蔽性强,防治困难,基因组信息缺乏。本研究的目的是获得荔枝蒂蛀虫的基因数据,寻求有效控制害虫的分子靶标。【方法】采用新一代高通量测序技术Illumina Hi SeqTM4000对荔枝蒂蛀虫进行转录组测序和生物信息学分析。【结果】经序列拼接获得68 996条unigenes。进一步利用七大公共数据库进行同源比对,注释了22 348 unigenes。注释到Nr数据库的unigenes数量最多,达27.01%,其中Nr注释的荔枝蒂蛀虫unigenes中与小菜蛾Plutella xylostella unigenes同源性最高,达34.1%。将unigenes与GO数据库比对发现,15 585条unigenes根据其功能大致可分为3类47亚类。KEGG pathways分析表明,7 272条unigenes定位为267个代谢通路。基因注释进一步筛选鉴定获得100个荔枝蒂蛀虫嗅觉相关基因;与鳞翅目相关气味结合蛋白基因联合分析发现,与荔枝蒂蛀虫气味结合蛋白基因直系同源的有18组,部分基因形成独立一簇。【结论】本研究首次获得了荔枝蒂蛀虫的转录组数据,研究结果为生物控制荔枝蒂蛀虫提供了重要的基础数据和候选分子靶标。荔枝蒂蛀虫独有的气味结合蛋白基因可能与其生境中特有的化学物质相关。     

中华大仰蝽转录组学研究及SSR新标记开发

【目的】中华大仰蝽Notonecta chinensis为中国和日本冲绳分布的重要水生天敌昆虫,可用于蚊虫的生物防治。本研究旨在建立中华大仰蝽转录组数据库,挖掘其基因信息。【方法】采用高通量测序平台Illumina NextSeq500对中华大仰蝽进行转录组测序、de novo组装及生物信息学分析;利用MISA软件基于转录组unigenes数据进行SSR新分子标记筛选。毛细管电泳检测SSR多态性。【结果】总计获得34782282条clean reads(NCBI SRA数据库登录号:SRR13259254),组装成37801条unigenes,N50为913 bp。将unigenes与已知数据库比对进行基因功能注释,分别有36474,32470,27781,35079和5638条序列注释到nr,Swiss-Prot,GO,eggNOG和KEGG数据库。通过GO数据库注释,unigenes的功能可分为生物学过程、细胞组分和分子功能三大类,其中参与细胞、细胞部分及结合功能的unigenes比例较大。eggNOG数据库注释结果显示,37801条unigenes归到25个基因家族,注释到未知功能的最多。KEGG代谢通路富集分析显示,5638条unigenes注释到245个代谢通路,注释到核糖体的数目最多。此外,用MISA软件在转录组测序数据中的37801条unigenes中搜索到3124个SSR位点(占总unigenes的8.26%),发生频率为7.07%。通过PCR筛选出16个SSR位点。7个中华大仰蝽地理种群3个位点NcCF/NcCR,NcKF/NcKR和NcLF/NcLR的多态信息含量(PIC)分别为0.870,0.902和0.857,具高度多态性。【结论】本研究成功获得了中华大仰蝽转录组数据,为其基因功能分析提供了分子理论基础;SSR新标记的开发为中华大仰蝽遗传多样性分析、隐存种鉴定及基因图谱构建提供了更丰富的候选分子标记。     

梨小食心虫幼虫中肠转录组及SSR分子标记分析

【目的】梨小食心虫Grapholitha molesta是一种重要的世界性果树害虫。昆虫中肠在食物消化、免疫反应、生长发育等方面发挥着重要作用。本研究旨在建立梨小食心虫幼虫中肠转录组数据库,挖掘其基因信息。【方法】利用高通量测序平台(Illumina HiSeq X Ten)对梨小食心虫幼虫中肠进行转录组测序、组装及生物信息学分析;进而利用转录组数据进行梨小食心虫幼虫中肠SSR分子标记鉴定。【结果】总计获得梨小食心虫幼虫中肠转录组96 419条unigenes,与公共数据库比对,共注释到57 300条unigenes。通过GO数据库注释,将unigenes的功能分为生物学进程、细胞组分和分子功能三大类55个功能区,其中参与细胞进程、细胞和细胞组分及结合功能的unigenes比例较大。KOG结果显示,10 090条unigenes归到25个基因家族,注释到一般功能预测的最多。KEGG数据库中,10 250条序列注释到232个代谢通路,注释到核糖体的基因数最多。此外,用MISH软件,在10 600条unigenes中搜索到12 690个SSR位点,通过PCR开发出6个SSR位点。【结论】本研究成功组装梨小食心虫幼虫中肠参考转录组,为以中肠为靶标的害虫防治提供理论依据。     

大垫尖翅蝗转录组分析

【目的】大垫尖翅蝗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条简单重复序列。【结论】该转录组测序分析将为进一步研究大垫尖翅蝗的基因功能分析及杀虫剂的抗药性机制分析奠定分子基础。     

羽化和性成熟时中华蜜蜂蜂王和雄蜂转录组分析

【目的】为了系统了解中华蜜蜂Apis cerana cerana蜂王和雄蜂转录组特征,丰富蜜蜂转录组数据信息。【方法】本研究利用高通量测序的方法分别检测中华蜜蜂蜂王和雄蜂刚出房、性成熟时期以及性成熟期蜂王生殖系统和雄蜂生殖系统之间转录组表达差异。【结果】经过测序获得质量值不低于20的碱基比例（Q20）均高于90%;所有reads组装成90 839个unigenes,平均长度1 549 bp;基于5个数据库(NR,Swiss Prot,GO,COG和KEGG)进行比对,共有45 112个unigenes被注释。差异基因表达分析发现,与刚出房时相比,性成熟的蜂王和雄蜂均在表皮蛋白、细胞色素P450、气味结合蛋白等家族基因上存在显著差异表达,而且这些差异表达基因与蜜蜂生长发育和性成熟过程中蜜蜂骨骼发育、生殖系统发育、嗅觉发育等方面有关;性成熟蜂王与性成熟雄蜂之间以及它们生殖系统之间在气味结合蛋白基因方面存在显著差异。【结论】结果表明,中华蜜蜂在性成熟过程中,体内大量基因的表达发生了变化。这些结果揭示了中华蜜蜂性成熟发育的整体基因表达特征,在得到大量转录组unigene序列的同时,获得了一批与蜜蜂性成熟有关的基因序列,为深入开展中华蜜蜂生长发育与繁殖研究提供了丰富的数据资源。     

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

【目的】建立花椒窄吉丁Agrilus zanthoxylumi成虫触角转录组数据库,挖掘嗅觉相关基因,为今后研究其触角的化学感受机制及生物防控提供理论支撑。【方法】采用高通量测序平台IlluminaNovaSeq 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 perkilobase per million mapped reads)值对嗅觉相关基因表达量进行分析。【结果】花椒窄吉丁雌雄成虫触角转录组测序共获得36 209条基因和90 982条转录本,其N₅₀分别为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值表明,嗅觉相关基因AzanOBP1和AzanOBP2在雌成虫触角中不表达,在雄成虫触角中微量表达;AzanOBP3在雄成虫触角中高丰度表达。【结论】首次获得了花椒窄吉丁成虫触角转录组数据,筛选到了花椒窄吉丁OBP, CSP, Or, IR和SNMP等的嗅觉相关基因。推测触角中高丰度表达的OBPs对雄成虫识别同类异性释放的信息素或寄主植物释放的挥发物起关键作用。研究结果可为花椒窄吉丁化学感受基因功能分析及嗅觉感受机制研究奠定分子基础。     

草地贪夜蛾雄性成虫和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中涉及解毒相关的通路显著富集。【结论】这些转录组数据为探究草地贪夜蛾的生长发育、嗅觉相关功能基因以及候选分子靶标提供资源信息。     

中华大仰蝽转录组学研究及SSR新标记开发

[目的]中华大仰蝽Notonecta chinensis为中国和日本冲绳分布的重要水生天敌昆虫,可用于蚊虫的生物防治.本研究旨在建立中华大仰蝽转录组数据库,挖掘其基因信息.[方法]采用高通量测序平台Illumina NextSeq500对中华大仰蝽进行转录组测序、de novo组装及生物信息学分析;利用MISA软件基于...     

携带与未携带松材线虫的松墨天牛转录组差异表达基因分析

【目的】本研究旨在建立携带与未携带松材线虫Bursaphelenchus xylophilus的松墨天牛Monochamus alternatus成虫各组织的转录组数据库,揭示松墨天牛对松材线虫响应的转录组整体表达特征。【方法】以携带和未携带松材线虫的松墨天牛成虫的表皮、气管和脂肪体为材料,采用Illumina HiSeq TM 2000测序平台开展转录组测序,利用Trinity软件对RNA-Seq数据进行从头组装,利用NCBI数据库进行基因注释。通过DEG-Seq分析携带和未携带松材线虫的松墨天牛成虫中的差异表达基因,对上调表达基因进行GO和KEGG代谢途径富集分析,并通过实时荧光定量PCR技术对部分差异表达的基因包括ENV,CDK1,HSP70-C,HSP70,HSP75,DUO,PABPN1和IGFP基因的表达水平进行验证。【结果】测序过滤后共获得55059条单基因簇(unigene),平均长度为1536 bp。将unigene与数据库中的序列进行BLASTX比对,成功注释24354条unigenes,其中4022条unigene在GO数据库中获得注释,6098条unigene在KEGG数据库中获得注释。经过GO与KEGG富集分析显示,松墨天牛成虫携带松材线虫后,差异表达基因主要为与其气管中压力调节、组织与DNA修复、激素响应方面相关的基因。实时荧光定量PCR分析结果显示CDK1,HSP70,HSP75,DUO,PABPN1和IGFP基因的表达量在携带松材线虫的松墨天牛气管中相比未携带松材线虫的松墨天牛气管中的显著上调。【结论】本研究初步阐明松墨天牛携带松材线虫后转录组的整体表达模式,为进一步研究松墨天牛的基因功能及对松材线虫的胁迫响应过程奠定了基础。     

叉角厉蝽触角转录组及嗅觉相关基因分析

叉角厉蝽Eocanthecona furcellata是一种在生物防治方面有重要作用和潜力的捕食性天敌昆虫,触角是昆虫进行信息交换的重要器官,而嗅觉相关基因则是调控天敌昆虫捕食行为的重要分子基础。为获得叉角厉蝽触角转录组数据库,挖掘叉角厉蝽嗅觉相关基因,本研究利用Illumina高通量测序平台对叉角厉蝽雌雄成虫与5龄若虫触角进行转录组测序。成功构建了叉角厉蝽触角转录组,获得了67 843条unigenes,N50长度为2 300 bp。与七大公共数据库比对注释到27 686条unigenes,其中NR数据库注释最多(33.33%),且与茶翅蝽Halyomorpha halys相似度最高(64.20%)。14 258条注释到GO数据库中,分为生物过程、细胞组分和分子功能3个大类42个亚类;KEGG代谢途径分析表明,7 703条形成282条代谢通路,其中被注释在信号传导通路中的unigenes最多(11.50%)。进一步基因注释分析,鉴定得到134个候选嗅觉相关基因,32个化学感应蛋白基因(Chemosensory protein genes,CSP),10个气味结合蛋白基因(Odorant binding protein genes,OBP),21个气味受体基因(Gustatory receptor genes,GR),48个嗅觉受体基因(Olfactory receptor genes,OR),17个离子型受体基因(Ionotropic receptor genes,IR),6个感觉神经元膜蛋白基因(Sensory neuron membrane protein genes,SNMP)。通过比较叉角厉蝽5龄若虫、雌雄成虫触角转录组,共筛选出7 324个差异表达基因,分析发现5龄若虫与成虫间的差异表达基因较多,雌雄成虫之间差异表达基因较少;利用FPKM值对OBP与CSP基因进行表达量分析发现,大部分OBP与CSP的表达量在雌雄间差异不大,而若虫与成虫相比基因表达量差异较大,除少部分基因在5龄若虫中高表达外,其余大部分基因均在雌雄成虫间高表达。本研究获得了叉角厉蝽触角转录组数据,并鉴定出候选嗅觉相关基因,结果为进一步研究叉角厉蝽的嗅觉感受机制及昆虫化学生态奠定分子基础。