中华蜜蜂PLA2基因的生物信息学分析及响应高温胁迫的表达模式

磷脂酶A2 (PLA2)对细胞膜流动性、膜蛋白活性等细胞生理功能有重要调节作用。高温等环境压力对细胞的生理代谢有较大的影响,其在生物体抗逆过程中的调控作用一直备受关注。本研究以中华蜜蜂Hymenoptera: Apidae: Apis cerana cerana PLA2基因序列为基础,对其蛋白结构进行预测,分析该基因在不同温度和高温不同胁迫时间下的表达差异,以揭示该基因在中华蜜蜂抗高温过程中的生理功能。结果显示,中华蜜蜂PLA2基因包含745 bp的开放阅读框,编码169个氨基酸,蛋白分子量为19.3 kDa,无跨膜结构,不属于膜蛋白。氨基酸同源序列比对结果显示,中华蜜蜂PLA2序列与蜜蜂科昆虫的相似性最高,与其他膜翅目昆虫的相似性存在差异。qRT-PCR结果显示,PLA2在35℃、40℃和45℃处理下的表达量较高,较高的温度能诱导该基因表达的增加。同时,PLA2的表达受高温胁迫时间的影响,较长时间的高温胁迫导致PLA2基因的表达增加。本研究结果表明PLA2在中华蜜蜂应对高温胁迫时发挥重要生理功能。     

中华蜜蜂DNA甲基化转移酶Dnmt3基因克隆及表达谱分析

为探究中华蜜蜂Apis cerana cerana的DNA甲基化模式, 本研究采用RT PCR技术克隆了中华蜜蜂DNA甲基化转移酶3（Dnmt3）基因(GenBank登录号为JQ740768); 采用荧光定量PCR检测不同发育时期工蜂（4日龄蛹, 1, 7和30日龄成年蜂及产卵工蜂）和蜂王（4日龄蛹, 1日龄蜂王和产卵蜂王）头部的Dnmt3基因mRNA的表达量。结果表明： 该基因cDNA序列全长2 277 bp, 编码758个氨基酸残基, 预测的蛋白分子量为88.24 kD, 等电点为7.85。将中华蜜蜂与其他物种的Dnmt3基因的结构域进行比对, 同时将该基因推导的氨基酸序列与其他物种的Dnmt3氨基酸序列进行同源性比对和系统发育分析, 发现与西方蜜蜂的Dnmt3序列一致性高达99%。该基因在工蜂和蜂王不同发育时期均有表达, 1日龄工蜂与7日龄工蜂中没有显著差异(P>0.05), 30日龄工蜂中的表达量显著高于前两者 (P<0.05); 蜂王蛹中的表达量显著高于工蜂蛹 (P<0.05); 1日龄的蜂王中的表达量显著高于1日龄的工蜂(P<0.05); 产卵工蜂与产卵蜂王中的表达量没有差异（P>0.05）。这种表达情况提示其可能与工蜂劳动分工及蜜蜂卵巢发育有关。     

中华蜜蜂dynactin p62基因的克隆及不同发育阶段表达分析

为探究中华蜜蜂Apis cerana cerana dynactin p62基因的表达特性, 本研究克隆了中华蜜蜂dynactin p62的基因组DNA序列(GenBank登录号： JX101463) 和mRNA序列（GenBank登录号： JX101464）; 采用荧光定量PCR检测了中华蜜蜂dynactin p62在不同发育时期（3日龄和6日龄幼虫、 刚羽化出房蜜蜂）三型蜂中mRNA的表达量。结果表明： 该基因基因组DNA序列全长为2 403 bp, mRNA序列全长为1 491 bp, 编码496个氨基酸残基, 预测的蛋白分子量为56.49 kD, 等电点为8.31。系统发育分析表明中华蜜蜂dynactin p62与西方蜜蜂Apis mellifera dynactin p62聚成一支。该基因在不同发育时期均有表达, 在雌性蜜蜂（蜂王和工蜂）中, 刚羽化成虫期的表达量显著高于幼虫期(P<0.05), 并且同一发育时期相比, 工蜂的表达量显著高于蜂王(P<0.05); 而该基因在雄蜂中表达量没有明显的规律性。这些结果提示该基因可能与中华蜜蜂级型分化有关。     

红光熊蜂卵黄原蛋白基因的cDNA全长序列克隆和表达分析

红光熊蜂Bombus ignitus Smith是许多经济作物和野生植物的重要授粉昆虫之一。卵黄原蛋白基因(vitellogenin,Vg)在昆虫的生殖调控中和行为方面起到重要的作用,本试验对Vg基因全长cDNA的克隆和测序及在蜂王、工蜂和雄性蜂三型蜂中的表达分析得出:Vg基因的全长cDNA为5 481 bp,GenBank中的登录号为FJ913883,有一个完整的开放阅读框(ORF),编码1 772个氨基酸,N-末端的前16个氨基酸为一个信号肽。接近C-末端区域存在保守的GL/ICG基元,其后含有9个半胱氨酸,而且DGXR位于GL/ICG基元上游18个氨基酸残基处。其氨基酸序列与韩国的熊蜂B.ignitus和B.hypocrita相似性高达95%,与西方蜜蜂Apis mellifera的相似性达到51%。Vg的mRNA首先在蜂王蛹期的白眼蛹(Pw)时期出现,其表达量在蜂王整个蛹期发育过程中呈上升趋势,且在黑眼蛹(Pbd)时期达到最高,在成年蜂的脂肪体中的表达量仍在升高,甚至更高。Vg也在工蜂蛹期的白眼蛹(Pw)时期被检测到,然后在整个蛹期发育过程中呈现上升趋势,在刚羽化出房时达到高峰,Vg的mRNA水平随着成年蜂日龄的增加而增加,到15日龄时达到最高,然后呈现下降趋势。对于雄性蜂,Vg的mRNA虽然卵黄原蛋白基因的mRNA水平几乎在整个蛹期发育阶段都表达,但是表达水平非常低,只有在刚羽化出房时期表达水平较高。     

Cloning and expression profiling of the DNA methyltransferase Dnmt3 gene in the Chinese honeybee, Apis cerana cerana (Hymenoptera:Apidae )

为探究中华蜜蜂Apis cerana cerana的DNA甲基化模式,本研究采用RT-PCR技术克隆了中华蜜蜂DNA甲基化转移酶3( Dnmt3)基因(GenBank登录号为JQ740768);采用荧光定量PCR检测不同发育时期工蜂(4日龄蛹,1,7和30日龄成年蜂及产卵工蜂)和蜂王(4日龄蛹,1日龄蜂王和产卵蜂王)头部的Dnmt3基因mRNA的表达量.结果表明:该基因cDNA序列全长2 277 bp,编码758个氨基酸残基,预测的蛋白分子量为88.24 kD,等电点为7.85.将中华蜜蜂与其他物种的Dnmt3基因的结构域进行比对,同时将该基因推导的氨基酸序列与其他物种的Dnmt3氨基酸序列进行同源性比对和系统发育分析,发现与西方蜜蜂的Dnmt3序列一致性高达99％.该基因在工蜂和蜂王不同发育时期均有表达,1日龄工蜂与7日龄工蜂中没有显著差异(P＞0.05),30日龄工蜂中的表达量显著高于前两者(P＜0.05);蜂王蛹中的表达量显著高于工蜂蛹(P＜0.05);1日龄的蜂王中的表达量显著高于1日龄的工蜂(P＜0.05);产卵工蜂与产卵蜂王中的表达量没有差异(P＞0.05).这种表达情况提示其可能与工蜂劳动分工及蜜蜂卵巢发育有关.     

西方蜜蜂发育基因AmWnt1的克隆鉴定及时空表达分析

本研究旨在克隆鉴定西方蜜蜂Apis mellifera发育相关基因AmWnt1,分析其在不同发育时期和刚出房工蜂不同组织的表达特征,为进一步研究Wnt1基因功能提供理论参考。根据NCBI中AmWnt1基因序列信息,利用Primer 6.0设计引物,RT-PCR扩增AmWnt1基因完整的CDS序列,进行生物信息学预测,用推导的氨基酸序列构建系统进化树;利用荧光定量PCR检测该基因在卵（1日龄、2日龄和3日龄）、幼虫（1日龄、3日龄和5日龄）、预蛹（1日龄和3日龄）、蛹（0日龄、2日龄、4日龄、6日龄和8日龄）、刚出房工蜂、哺育蜂和采集蜂以及刚出房工蜂8个组织中相对表达量。克隆获得西方蜜蜂的Wnt1基因CDS序列,命名为AmWnt1,上传NCBI,获得GenBank登录号MT993937。全长1 239 bp,编码412个氨基酸,预测等电点为9.48,相对分子质量为46.40313 kDa。序列比对和系统进化树结果表明：AmWnt1蛋白与其它膜翅目昆虫聚为一类,其中和东方蜜蜂Apis cerana亲缘关系最近,序列相似度为99.50%。时空表达谱结果表明：AmWnt1基因在各个时期中均有表达,且在胚胎后期表达量最高,预蛹期和蛹前期表达量相对较高,其它时期表达量相对较低;AmWnt1基因在刚出房工蜂头、胸、触角表达量高于其它组织。AmWnt1可能参与西方蜜蜂胚胎晚期的神经系统发育和化蛹过程中的四肢发育等关键历程,为进一步研究AmWnt1功能提供了理论参考。     

平颏海蛇毒磷脂酶A2基因的多样性

通过对平颏海蛇毒腺cDNA文库的随机测序和PCR筛选,分离到3种未曾报道的磷脂酶A2（PLA3）cDNA,分别命名为PLA2－8,PLA2－9,PLA2－384。这些cDNA均含有完整的读码框（长度分别为456bp、438bp和438bp）,编码含27个氨基酸残基的信号肽和长度分别为4．8、7．8和8．4。氨基酸序列分析和二级结构预测表明,3种PLA2均属I类磷酯酶A2家族,其中PLA2－8与来自澳大利亚虎蛇Notechis scutatus scutatus的PLA2有81％的同源性;而PLA2－9和PLA2－384则与来自海蛇Enhydrina schistosa的肌毒PLA2有90％的同源性,这不仅预示着他们在功能上可能存在某些相似性,也反映出平颏海蛇PLA2的结构多样性。这组海蛇PLA2同功酶基因的成功克隆为进一步揭示PLA2的结构与功能关系及的分子机制提供了新的信息。     

西方蜜蜂磁受体基因AmMagR的克隆及表达分析

本研究旨在克隆西方蜜蜂Apis mellifera磁受体基因AmMagR的序列,并分析该基因在工蜂不同日龄、不同组织中的表达特征,以期为该基因的功能研究提供理论基础。以西方蜜蜂工蜂为材料,提取其总RNA,通过RT-PCR技术克隆西方蜜蜂MagR基因的序列;利用多种生物信息学软件对其核酸及氨基酸序列分析;采用实时荧光定量PCR技术(qPCR)分析在西方蜜蜂工蜂不同日龄(1日龄、5日龄和21日龄)和不同组织(头部、胸部和腹部)中MagR基因的相对表达量。克隆到西方蜜蜂MagR基因,并命名为AmMagR(GenBank登录号:MH635411),其序列长度为748 bp,编码区长390 bp,编码130个氨基酸,其蛋白质分子量为14.142 kDa,理论等电点为9.06,无信号肽,无跨膜结构,且在第24～126位氨基酸之间得到一个结构域Fe-S_biosyn。系统发育树显示,西方蜜蜂AmMagR与小蜜蜂Apis florae AfIscA(IscA别称MagR)、中华蜜蜂Apis cerana cerana AcIscA基因聚成一支。荧光定量PCR结果表明,AmMagR基因在西方蜜蜂不同日龄的工蜂头部、胸部和腹部均有表达。5日龄工蜂的相对表达量最高,5日龄工蜂头部和胸部的AmMagR基因表达量显著高于1日龄(P0.05)和21日龄(P0.05),腹部AmMagR基因表达量显著高于21日龄(P0.05),但与1日龄的比较差异不显著(P0.05)。1日龄、5日龄和21日龄工蜂胸部的AmMagR基因表达量均显著高于头部和腹部(P0.05);1日龄工蜂头部AmMagR基因表达量均高于腹部(P0.05),但5日龄和21日龄工蜂AmMagR基因在头部与腹部的表达量无显著差异(P0.05)。明确了该基因在西方蜜蜂工蜂不同日龄和不同组织中的表达模式,并推测AmMagR可能参与了西方蜜蜂的定位、归巢过程。     

中华蜜蜂感觉神经元膜蛋白基因克隆、组织表达分析及原核表达

为明确中华蜜蜂Apis cerana cerana嗅觉形成中重要功能因子的信号转导通路, 本研究利用RT-PCR方法, 克隆了中华蜜蜂感觉神经元膜蛋白 (sensory neuron membrane protein, SNMP) 基因编码区, GenBank登录号为KC012595, 命名为AccSNMP1。序列分析表明, 该编码区开放阅读框长1 563 bp, 编码520个氨基酸, 推测的编码蛋白的相对分子量和等电点分别为58.02 kD和5.83。同源性比较发现, 中华蜜蜂AccSNMP1与其他昆虫感觉神经元膜蛋白基因同源性差异很大, 在氨基酸水平上与西方蜜蜂Apis mellifera SNMP基因一致性达99.2%, 与熊蜂Bombus impatiens SNMP基因一致性达90.9%, 而与赤拟谷盗Tribolium castaneum SNMP基因一致性仅为22.7%。系统发育树显示中华蜜蜂与西方蜜蜂遗传距离最近。实时荧光定量PCR结果分析表明, AccSNMP1在触角中表达量最高, 在足中表达量较高, 与胸、 腹、 头（去除触角和喙）、 喙中表达量相比差异显著（P<0.05）。构建原核表达载体pEASY-E1-AccSNMP1, 经IPTG诱导, 中华蜜蜂感觉神经元膜蛋白在大肠杆菌Escherichia coli BL21 （DE3）中高效表达。结果为进一步研究AccSNMP1在中华蜜蜂体内的作用机理奠定了基础。     

细胞色素b基因序列与11种熊蜂的系统进化

通过比较一段357b的细胞色素b的序列,分析了熊蜂属5亚属11种熊蜂系统发育关系。5亚属包括Bombus( 5种）,Thoracobombus(3种）,Mendacibombus(1种）,Fervidobombus(1种）和Pyrobombus(1种）,该序列有65个单变异态位点和71个简约信息多态位点,翻译成119个氨基酸序列后有45个氨基酸变异位点。根据P－距离构建的邻接树（NJ tree)和最大简约树（MP tree)都显示同样的结果：Mendacibombus (B.avinovielllus)分化最早;Fervidobombus(B.pensylvanicus)次之;Pyrobombus(B.impatiens)和Bombus形成姊妹群;Bombus亚属是单系群,其中B.ignitus在所研究的5个种中分化最早。     

Molecular cloning and characterization of a venom phospholipase A2 from the bumblebee Bombus ignitus

Phospholipase A₂ (PLA₂) is one of the main components of bee venom. Here, we identify a venom PLA₂ from the bumblebee, Bombus ignitus. Bumblebee venom PLA₂ (Bi-PLA₂) cDNA, which was identified by searching B. ignitus venom gland expressed sequence tags, encodes a 180 amino acid protein. Comparison of the genomic sequence with the cDNA sequence revealed the presence of four exons and three introns in the Bi-PLA₂ gene. Bi-PLA₂ is an 18-kDa glycoprotein. It is expressed in the venom gland, cleaved between the residues Arg44 and Ile45, and then stored in the venom sac. Comparative analysis revealed that the mature Bi-PLA₂ (136 amino acids) possesses features consistent with other bee PLA₂s, including ten conserved cysteine residues, as well as a highly conserved Ca²⁺-binding site and active site. Phylogenetic analysis of bee PLA₂s separated the bumblebee and honeybee PLA₂ proteins into two groups. The mature Bi-PLA₂ purified from the venom of B. ignitus worker bees hydrolyzed DBPC, a known substrate of PLA₂. Immunofluorescence staining of Bi-PLA₂-treated insect Sf9 cells revealed that Bi-PLA₂ binds at the cell membrane and induces apoptotic cell death.     

Sensitivities of three bumblebee species to four pesticides applied commonly in greenhouses in China

Abstract  This study was conducted to evaluate the effects of four pesticides (Mosplian, Kingbo, Score and Lvrtong) applied commonly in greenhouses in China, on three bumblebee species ( Bombus hypocrita , Bombus ignitus and Bombus patagiatus ). The study used a contact experiment and oral toxicity LD₅₀ values. The results showed that the mortality for B. hypocrita after contacting the four pesticides was significantly lower than B. patagiatus and B. ignitus , but there was no significant difference between B. patagiatus and B. ignitus . The oral toxicity median lethal dose (LD₅₀) value of Mosplian to B. hypocrita (0.0028 μg active ingredient/bee) was significantly higher than that to B. ignitus (0.0023 μg active ingredient/bee) and B. patagiatus (0.0021 μg active ingredient/bee). Of the bee species, it can be concluded that B. hypocrita was the least susceptible to the four pesticides. The mortality rates of each bumblebee species after contact with Mosplian were significantly higher than for the other three pesticides and the control group. For Kingbo, the rates were significantly higher than the control group, but Score and Lvrtong exposed groups showed no significant increase in mortality relative to the control group. It can therefore be concluded that the pesticides differ in their negative influences on bumblebees, and that Mosplian is the most harmful.     

不同蜂在温室草莓园的访花行为和传粉生态学比较

应用意大利蜜蜂、本土小峰熊蜂和引自欧洲的地熊蜂3种蜂在北京为温室草莓传粉,对其传粉行为和生态特性进行研究。结果表明:意大利蜜蜂的出巢温度、工作起点温度、平均单花停留时间和访花间隔时间均显著高于小峰熊蜂和地熊蜂(P<0.01),而这些行为在2种熊蜂之间差异不显著(P>0.05);3种蜂携带花粉回巢的比例依次为意大利蜜蜂>小峰熊蜂>地熊蜂,而其携带的花粉活性依次为意大利蜜蜂>小峰熊蜂=地熊蜂;2种熊蜂开始活动的时间比意大利蜜蜂早,日活动时间比意大利蜜蜂长;地熊蜂在一天中的活动频率最高(P<0.01);经3种蜂授粉的草莓果实营养品质差异不显著(P>0.05),与对照相比也无差异,而经过3种蜂授粉的草莓果实畸形率降低、单果质量增大;经2种熊蜂授粉的草莓果实畸形率低于意大利蜜蜂(P<0.05)。本土小峰熊蜂能够有效替代意大利蜜蜂或地熊蜂为温室草莓授粉。     

管氏肿腿蜂毒液过敏原3基因的克隆及表达分析

为研究管氏肿腿蜂Scleroderma guani毒液过敏原3基因SgA3的功能，通过RT-PCR克隆SgA3基因，采用荧光定量PCR分析其在不同发育阶段和组织中的表达特征，并利用载体pSUMO-Mut对其进行原核表达。克隆获得SgA3基因开放阅读框长699 bp，编码233个氨基酸，其中信号肽位于N端第1~23位氨基酸。多序列比对分析发现，SgA3与丽蝇蛹集金小蜂Nasonia vitripennis、粉蝶盘绒茧蜂Cotesia glomerata、红火蚁Solenopsis invicta和常见黄胡蜂Vespula vulgaris过敏原3的氨基酸序列一致性分别为50.44%、50%、48.89%和47.83%。荧光定量PCR结果表明，SgA3基因在雌成虫中高表达，且在毒液器官中的表达量显著高于雌成虫其他组织和雄成虫中的表达量。利用pSUMO-Mut表达载体，成功表达并纯化得到高纯度的SgA3重组蛋白。该研究结果为进一步研究SgA3基因的功能奠定了基础。     

梨小食心虫普通气味受体基因GmolOR20的克隆及表达分析

【目的】通过克隆梨小食心虫Grapholita molesta触角中的普通气味受体(odorant receptor, OR)OR20基因,明确其在不同发育期及成虫不同组织中的表达特征,为进一步研究其功能提供理论依据。【方法】根据梨小食心虫雌虫触角转录组数据,利用RT-PCR克隆梨小食心虫OR20基因的完整开放阅读框;采用qRT-PCR检测该基因在不同发育期(卵、1-5龄幼虫、蛹和雌雄成虫)、成虫不同组织(触角、去除触角的头、胸、腹、足、翅)以及不同日龄(1, 3, 5和7日龄)成虫触角中的表达量。【结果】克隆获得梨小食心虫GmolOR20基因cDNA序列(GenBank登录号:MH898864)。该基因完整开放阅读框为1 284 bp,编码427个氨基酸,预测蛋白分子量为49.83 kD,理论等电点为8.57,具有7个跨膜结构域。序列比对和系统进化树结果表明,梨小食心虫GmolOR20与苹果蠹蛾Cydia pomonella CpomOR15和豆荚小卷蛾Cydia nigricana CnigOR15亲缘关系较近,氨基酸序列一致性分别为87%和84%。发育表达模式结果显示,GmolOR20在不同发育时期均有表达,雌雄成虫中的表达量显著高于其他发育期的表达量(P<0.05),但雌、雄虫间的表达量差异不显著。组织表达模式结果表明,GmolOR20主要在成虫触角中高丰度表达,且雌虫触角中的表达量极显著高于雄虫触角中的表达量(P<0.01);GmolOR20在不同日龄成虫的触角中均有表达,且在1和3日龄成虫触角中的表达水平显著高于其他日龄(P<0.05)。【结论】根据GmolOR20基因的表达谱分析结果,推测GmolOR20可能参与梨小食心虫对植物挥发物和性信息素的识别。     

意大利蜜蜂和小峰熊蜂在温室桃园的传粉行为及其影响因素

意大利蜜蜂Apis mellifera ligustica和小峰熊蜂Bombus hypocrita是我国设施农业中的两种主要授粉昆虫, 二者的传粉方式和传粉效率因不同作物而各有不同。为了选用最理想的蜂种为温室作物授粉, 提高作物的授粉效率, 我们于2008-2010年连续3年在北京的温室桃园进行了意大利蜜蜂和小峰熊蜂的传粉行为及其影响因素的研究。结果表明, 两种蜂的日活动规律不同, 和意大利蜜蜂相比, 小峰熊蜂的活动起点温度低、 日工作时间长、 单花访问时间长, 采集蜂在温室内距离蜂巢不同距离之间扩散比较均匀。而意大利蜜蜂采集蜂的数量随授粉半径的增大而明显减少。温室内环境因子对两种蜂传粉活动的影响基本一致, 温度对两种蜂的传粉活动影响最大, 其次是湿度、 单花花蜜浓度和光照强度, 而单花花蜜体积对蜂活动无明显影响。本研究认为对于温室桃授粉应优先选用授粉性能稳定的小峰熊蜂, 并且适当调节温室内环境条件, 以提高其授粉效率。     

Molecular cloning and fibrin(ogen)olytic activity of a bumblebee (Bombus hypocrita sapporoensis) venom serine protease

Although several bee venom serine protease genes have been previously described, fibrin(ogen)olytic activity of these serine proteases has been reported for only two bumblebees to date, Bombus ignitus and B. terrestris. Here, we cloned venom serine proteases from the other bumblebee species, B. hypocrita sapporoensis and B. ardens ardens. The venom serine protease genes of B. h. sapporoensis and B. a. ardens consist of 358 amino acids and 357 amino acids, respectively. We compared the predicted mature protein sequences of these serine protease genes to those previously reported for other bees. A phylogenetic analysis shows that B. h. sapporoensis venom serine protease is further immediately close to B. ignitus and B. terrestris venom serine proteases, excluding the venom serine protease of B. a. ardens. Using B. h. sapporoensis venom serine protease (Bs-VSP), we identified that Bs-VSP acts as a fibrin(ogen)olytic enzyme. We also found that Bs-VSP activates prothrombin and directly degrades fibrinogen into fibrin degradation products. Our results further define roles for bumblebee venom serine proteases as fibrin(ogen)olytic agents.     

松树蜂毒腺转录组分析

【目的】构建入侵种松树蜂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)。【结论】本研究获得的松树蜂毒腺转录组数据为松树蜂毒液组分的鉴定和生物学功能的研究奠定了一定的理论基础。