基于高通量测序的文冠果转录组分析

应用Illumina Hi-seqTM2000高通量测序技术对文冠果花芽进行转录组分析。共获得N50为1 180bp、平均长度为686bp的unigene 58 311条。与公共数据库Nr和Swiss-Prot同源性比较后发现37 047条unigene获得基因注释,另有21 264条unigene未被注释。利用COG数据库将unigene分成25类。通过GO分类和KEGG Pathway富集性分析,将unigene分别归类于55个GO term和128个代谢途径。此外,在9 794条unigene中共搜索到12 213个SSR位点,单核苷酸重复基元出现频率最高(34.95%),其次分别为二核苷酸(32.74%)和三核苷酸(28.64%)。在获得的unigene中发掘出涉及4个开花调控途径(光周期途径、春化途径、GA途径和自主途径)多个基因的同源序列。研究结果可在一定程度上解析文冠果花芽形态分化的分子调控模式与机制。     

独行菜种子转录组的高通量测序及分析

独行菜种子为我国传统常用中药,从中已提取出多种药用活性成分,但目前尚不清楚其次级代谢过程中这些活性物质合成的遗传基础。采用Illumina Hiseq^TM 2000高通量测序平台对独行菜种子转录组进行测序,经de novo组装后获得40 303条unigene。进一步利用六大公共数据库进行同源比对,注释了27 935条unigene。研究发现,534个基因参与了独行菜次生物质的合成和代谢,其中在芥子苷、黄酮类和芪类化合物生物合成途径中的unigene分别有4个、19个和69个,在苯丙氨酸代谢途径中的unigene有92个。这些基因可能参与独行菜种子药性活性物质的生物合成,并分析获得了参与上述合成代谢途径的13个关键基因的同源序列。另外,从转录组序列中搜索到6 304个SSR位点,分布于5 306条unigene中,出现频率为15.64%。研究结果不仅为挖掘独行菜种子药用次生代谢物生物合成关键基因提供了基础数据信息,而且为独行菜遗传多样性研究和分子标记开发奠定了分子基础。     

火龙果转录组测序、基因表达与功能分析

利用高通量测序技术对火龙果(Hylocereus undulatus Britt)红肉品种‘大红二号’的花芽、果实和枝条不同发育阶段的基因表达进行研究。结果显示,转录组测序共获得468.68 Gb原始数据(Raw data),从头组装获得239 152条转录本和162 519条unigene,约53.74%的unigene得到注释。分别在43 506条和16 251条unigene中检测到600 283个SNP位点和56 147个SSR位点。基因表达分析结果表明,在火龙果不同组织Fl510、Fl513、Fl514、Fl518、F711、F715、S513、S419中分别有31、7、5、152、17、63、17、8个特异表达的unigene。通过GO和KEGG富集分析,发现了一些组织特异的GO条目和代谢通路,如在Fl510中富集的类萜骨架生物合成代谢通路等。本研究还对参与花发育的候选基因进行了鉴定和表达分析,他们包括COL基因、FT-like基因、分生组织决定基因和器官决定基因等。     

火龙果转录组测序、基因表达与功能分析

利用高通量测序技术对火龙果（Hylocereus undulatus Britt）红肉品种‘大红二号’的花芽、果实和枝条不同发育阶段的基因表达进行研究。结果显示,转录组测序共获得468.68 Gb原始数据（Raw data）,从头组装获得239 152条转录本和162 519条unigene,约53.74%的unigene得到注释。分别在43 506条和16 251条unigene中检测到600 283个SNP位点和56 147个SSR位点。基因表达分析结果表明,在火龙果不同组织Fl510、Fl513、Fl514、Fl518、F711、F715、S513、S419中分别有31、7、5、152、17、63、17、8个特异表达的unigene。通过GO和KEGG富集分析,发现了一些组织特异的GO条目和代谢通路,如在Fl510中富集的类萜骨架生物合成代谢通路等。本研究还对参与花发育的候选基因进行了鉴定和表达分析,他们包括COL基因、FT-like基因、分生组织决定基因和器官决定基因等。     

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

多效生长因子相关基因的生物信息学分析

目的:用生物信息学方法分析多效生长因子(PTN)潜在的分子功能。方法:利用由美国亚利桑那癌中心提供的生物信息学数据库,对前期用小鼠全基因组表达谱芯片检测到的Ptn相关基因进行生物信息学分析,通过GO Terms分析这些基因所属的功能群体,用Pathway Miner分析这些基因参与调控的信号通路。结果:370个由芯片检测得到的Ptn相关基因中,在GO Terms数据库中找到231个基因,其中参与细胞成分构成的基因占31.83%,具有分子功能的基因占35.34%,而参与生物学过程的基因占32.83%;在Pathway Miner数据库中找到105个基因。这些基因相关的信号通路有230条,分别属于细胞和调控过程通路以及代谢通路。结论:PTN是一个重要的细胞因子,可能参与机体的免疫与防御反应、炎症反应,以及细胞的增殖、凋亡调控等。     

基于转录组的苹小卷叶蛾杀虫剂靶标及解毒代谢相关基因分析

【目的】建立苹小卷叶蛾Adoxophyes orana转录组数据库,挖掘杀虫剂靶标及解毒代谢相关基因。【方法】采用Illumina HiSeq~(TM) 2000高通量测序技术对苹小卷叶蛾进行转录组测序,挖掘并分析杀虫剂靶标基因;利用qPCR检测6个杀虫剂靶标基因在苹小卷叶蛾卵、幼虫、蛹和成虫各不同发育阶段的表达;挖掘并分析苹小卷叶蛾转录组中解毒代谢相关基因的代谢通路及进化关系。【结果】通过组装有效序列共获得48 610条unigene(GenBank登录号:GGMW00000000)。挖掘鉴定到155个杀虫剂靶标unigene;qPCR结果显示,1个蜕皮激素受体(ecdysone receptor, ECR)、2个乙酰胆碱酯酶(acetylcholinesterase, AChE)、1个氯离子通道蛋白(chloride channel, CLC)、1个几丁质酶(chitinase, CS)和1个鱼尼丁受体(ryanodine receptor, RyR)基因在苹小卷叶蛾不同发育阶段均存在表达差异。挖掘鉴定到69个羧酸酯酶(carboxylesterase, CarE)unigene、66个谷胱甘肽S-转移酶(glutathion S-transferase, GST)unigene和205个细胞色素P450(cytochrome P450)unigene等解毒代谢相关基因,共鉴定20个CarE unigene, 32个GST unigene和30个P450 unigene与有毒物质代谢相关的通路有关。基于氨基酸序列对具有完整ORF的unigene聚类分析结果显示:12个CarEs中9个为G类,即鳞翅目保幼激素类;分别有10个AoGSTs属于Delta和Epsilon亚家族;18个P450全部聚到CYP3集团。【结论】该研究有助于苹小卷叶蛾杀虫剂靶标基因的挖掘及抗药性的研究。     

东亚特有珍稀蕨类植物岩穴蕨(碗蕨科)高通量转录组测序及分析

岩穴蕨(Monachosorum maximowiczii)隶属于碗蕨科稀子蕨属, 是东亚中高海拔地区所特有的珍稀濒危植物。为了在分子水平对该物种有进一步的认识, 本文首次利用二代高通量测序技术(RNA-seq)对岩穴蕨进行了转录组测序分析。通过Illumina Hiseq 2500测序平台, 共计获得4.95 Gb原始数据(raw data), 经过滤后得到4.83 Gb有效数据(clean reads), 并进行从头组装得到了101,448条unigene。其中, 54,106条unigene预测到完整的开放阅读框。我们利用目前已知的51个植物基因组数据, 对岩穴蕨的unigene进行了详尽的功能注释, 并通过GO、COG、KEGG注释进一步了解了这些编码基因的作用方式、特征以及所参与的代谢通路。同时, 转录因子分析结果也为岩穴蕨的环境适应机制研究提供了初步线索。     

基于高通量测序的辽东栎转录组学研究

应用Illumina Solexa Hiseq 2000高通量测序技术对辽东栎的芽、花、叶及果实的混合样品进行转录组测序,结果共获得3.8 Gb的有效数据。应用Trinity软件对有效序列从头拼接去重复后,共获得95 800条unigene,总长度为73.57 Mb,最大长度、平均长度和N50分别为11 284 bp、768 bp和1 373 bp。利用Blastx与公共数据库Nr和Swiss-Prot的同源性比较(E值1×10-5)发现,38 163条unigene未发现与公共数据库中的序列具有同源性。通过KEGG数据库中参与淀粉合成与代谢的pathway分析,共发掘出67条参与淀粉合成的unigene,编码9个关键酶。此外,在13 380条unigene中共搜索到15 901个SSR位点,其中二核苷酸和三核苷酸的重复类型占所有SSR位点的98.16%。     

基于Solexa高通量测序的黄曲条跳甲转录组学研究

黄曲条跳甲Phyllotreta striolata (Fabricius) 是十字花科蔬菜的重要害虫。为深入了解其遗传信息, 本研究应用新一代高通量测序技术Illumina’s Solexa平台对黄曲条跳甲成虫的转录组进行测序, 并结合SOAPdenovo拼接聚类等分析软件, 获取大量的EST和挖掘功能基因。本文最终获得了4 924条序列重叠群（contig）, 其中包含2 209种与黑腹果蝇Drosophila melanogaster蛋白基因具直系同源的独立基因（unigene）和610种黄曲条跳甲物种特有的unigene。结合Gene Ontology （GO）数据库进行分析, 发现大部分的unigene具结合能力（binding capability）和催化活性（catalytic activity）; 上百种unigene可聚类于生物学过程分类中的配子发生、 生殖腺发育和交配行为等重要功能。另外, 结合KEGG Pathway数据库分析发现, 共有363种unigene参与或涉及了40种代谢路径, 其中生物钟调控路径和植物次生代谢物路径等相关基因的发现, 有助于深入研究黄曲条跳甲行为发生的内在机理。Solexa高通量测序技术作为昆虫功能基因组研究的重要手段, 为发掘黄曲条跳甲功能基因发挥了重要作用, 也为在分子水平上研发黄曲条跳甲的防治新策略提供了更翔实的基因信息。     

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

【目的】建立绿豆象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值对基因表达量进行评估。【结论】本研究获得了绿豆象触角转录组数据,为进一步研究绿豆象的基因功能分析及嗅觉感受机制奠定分子基础。     

牦牛体外成熟卵母细胞差异转录组学研究

卵母细胞体外成熟环节是现代繁殖育种技术的基础环节,体外成熟的MⅡ期卵母细胞的质量直接影响了卵母细胞的后续受精及受精卵卵裂.由于卵母细胞的发育成熟是一个涉及大量基因转录调控的复杂过程.因此,转录组学研究是了解卵母细胞发生发育的关键.本研究以高海拔地区的特有牛属动物牦牛(Bos grunniens)卵母细胞作为研究对象,应用RNA-seq技术对牦牛GV期卵母细胞及体外成熟MⅡ期卵母进行转录组学测序及比对分析.经深度测序后,各获得一个包含51380686条过滤测序序列(cleanreads),4624261740碱基(bp)的GV期牦牛卵母细胞测序文库和一个包含50303412条过滤测序序列(clean reads),4527307080碱基(bp)的MⅡ期牦牛卵母细胞测序文库.基因覆盖率统计表明,GV期和MⅡ期文库中分别有16719条和16339条牦牛基因得到转录.通过比较分析MⅡ和GV期转录组数据,共筛选出4767个差异表达基因,其中1418个基因表达量上调和3349个基因表达量下调.GO功能分类注释显示,差异基因与"代谢过程"等生物学过程存在密切关联,同时结合分子功能相关类别等分子事件表现活跃.按表达量上调和下调分别对差异基因进行KEGG通路互作分析,结果表明,凋亡通路、内吞通路及代谢通路是上调的关键调控通路,而代谢通路、丙酮酸代谢、黏着斑及氧化磷酸化是下调的关键调控通路,其中代谢通路在上调及下调通路互作中均起关键调控作用.进一步对代谢通路进行分析,结果表明,组氨酸代谢、甘氨酸、丝氨酸和苏氨酸代谢、乙醛酸和二羧酸代谢及谷胱甘肽代谢是整个代谢通路的关键调控节点。此外,本研究发现,在牦牛卵母细胞体外成熟过程中,与免疫相关的白细胞介素家族及干扰素家族基因发生了显著差异表达.该研究结果为进一步解析牦牛卵母细胞发育的分子机制及全面了解牦牛繁殖的特异性提供了基础.     

基于高通量测序的云南松转录组分析

采用新一代高通量测序技术平台Illumina Hiseq 2 000对云南松转录组测序,得到的数据进行de novo组装,获得80 000条Unigenes,N50为1 881 nt、平均890 nt。与公共数据库进行比对,注释到NR、NT、Swiss-Prot数据库的Unigenes分别为43 434、46 415、29 418条。将Unigenes与COG数据库比对,有14 792条Unigenes成功注释,根据功能大致分成25类;与GO数据库比对,有26 743条Unigenes获得注释,按功能分为细胞组分、分子功能和生物过程3大类55亚类,其中参与的生物过程较多;以KEGG数据库参考,有25 873条Unigenes参与128条代谢途径分支,以代谢相关的通路较为集中,并找到与木质素合成关键酶的Unigenes。这些研究极大地扩充了云南松的基因资源,将有助于云南松基因的发掘与利用、分子标记的开发及其种质资源遗传改良的研究等。     

药用资源植物山莨菪的转录组信息分析

为了增强对资源植物山莨菪的深入了解,本研究采用高通量测序技术对山莨菪进行转录组测序分析,经过处理得到71 463个Unigenes。通过与多个数据库进行比对,对基因进行分类和分析注释,最终成功获得注释的基因有47 624条。将Unigenes比对到KOG蛋白质库中,有13 110个基因被注释,共有26个子类;比对到NR库中后有39 621个Unigenes被注释;转录本与Swissprot、Tr EMBL的比对结果得到GO功能注释信息,注释得到的29 309个Unigenes可被分为分子功能、生物学过程和细胞组分3个大类,62个子类;以KEGG数据库为参考,3 679条基因被注释,参与的代谢通路可归为4个大类,分别是代谢相关的通路、遗传信息处理、细胞过程、环境信息处理,其中与代谢相关的通路最多,约占所有代谢通路的一半。对山莨菪的药用活性成分的代谢通路及相关Unigenes数量和类型的统计结果表明,与生物碱相关的代谢通路最多,萜类和苯丙素类所对应的Unigenes数量最多。另外,结果还检测到31 382个SNP位点,6种SSR重复类型,其中单碱基重复类型所占的比例最高,每百万碱基中出现的单碱基重复的SSR个数有56. 52个,占45. 30%。该结果丰富了山莨菪的转录组信息数据,为该物种分子生物学方面的研究奠定了基础,有助于进一步开展对山莨菪的合理保护及开发利用工作。     

基于转录组测序的滇山茶花叶呈色机理分析

该研究采用Illumina Hi-Seq2500高通量测序技术,对滇山茶的花叶(绿叶区和金叶区)的cDNA进行了转录组测序,分析滇山茶的花叶呈色机理。测序结果得到了228 862条单基因序列,筛选得到了4 851条差异表达基因。2 291条差异表达基因在GO数据库中有功能注释,1 826条差异表达基因被注释到COG数据库,这些差异表达基因被分为23个主要功能类别,其中大部分基因与滇山茶的生长和代谢有关。1 363条差异表达基因注释到了KEGG数据库,进一步筛选出了差异表达基因注释序列最多的15条代谢通路,其中卟啉和叶绿素代谢途径中的GluRS、δ-ALAD、ALAS基因和类黄酮合成途径中ANS、LAR、CHS基因与滇山茶的花叶形成密切相关。研究表明,滇山茶的花叶呈色是由于叶色相关代谢通路中大量的基因表达量发生变化引起的,并不是由单基因突变所致,病毒诱导的多个代谢途径基因沉默可能是滇山茶花叶呈色的主要原因。该研究丰富了滇山茶的转录组信息,并为滇山茶的彩叶品种培育提供了一定的参考。     

沙葱萤叶甲卵滞育的转录组学分析

【目的】建立沙葱萤叶甲Galeruca daurica滞育卵转录组数据库,挖掘卵滞育相关的基因以及代谢和信号通路,在转录组水平探讨卵滞育的分子机制。【方法】采用Illumina NovaSeq6000高通量测序平台对沙葱萤叶甲滞育卵与解除滞育卵进行转录组测序,并进行生物信息学分析;利用DESeq软件分析沙葱萤叶甲滞育卵与解除滞育卵中的差异表达基因,对差异表达基因进行KEGG通路富集分析;利用qRT PCR技术对10个差异表达基因的表达模式进行验证。【结果】基于沙葱萤叶甲滞育卵与解除滞育卵转录组测序结果,共获得53 389个unigene,其中差异表达基因2 145个,24个差异表达基因与保幼激素信号及脂肪酸生物合成和降解相关。与解除滞育卵相比,滞育卵转录组中1 297个基因上调表达,富集于124条KEGG通路,其中核糖体通路显著富集;848个基因下调表达,富集于73条KEGG通路,其中MAPK信号通路和糖胺聚糖生物合成通路显著富集。qRT-PCR结果表明,随机选取的10个差异表达基因的表达趋势与RNA-Seq转录组测序结果完全一致。【结论】保幼激素,脂肪酸生物合成和降解,核糖体,MAPK信号及糖胺聚糖生物合成等通路可能在沙葱萤叶甲卵滞育调控中起着重要的作用。     

意大利蝗转录组及性腺发育相关基因分析

【目的】意大利蝗Calliptamus italicus是新疆荒漠、半荒漠草原优势危害种类之一,其生殖相关基因信息缺乏。本研究对意大利蝗精巢和卵巢组织进行高通量转录组测序,旨在揭示意大利蝗转录组特征,并获得与性腺发育相关的基因数据。【方法】利用Illumina高通量测序平台(Hi Seq/Mi Seq)对意大利蝗进行了转录组测序,并进行序列分析。【结果】获得718 872条unigenes,平均长度631 bp,N50为1 186 bp。通过同源性比对,成功注释254 597条unigenes,其中,Nr数据库注释占28.87%,比例最高。Nr数据库注释的unigenes与黑褐草蚁Lasius niger的相似基因序列最多,为10.80%。与GO数据库比对发现,意大利蝗转录组unigene可分为3大类:涉及分子功能的unigenes有31 392条,涉及细胞组分的unigenes有20 586条,与生物学过程有关的unigenes有57 014条。KEGG pathways分析表明,23 666条unigenes归属于251个通路。涉及生殖系统发育的通路有卵母细胞减数分裂、胰岛素信号通路、Wnt信号通路、促性腺激素释放激素信号通路、转化生长因子β信号通路以及泛素介导的蛋白水解、黄体酮介导的卵母细胞成熟、昆虫激素生物合成等。进一步筛选得到部分与性腺发育相关的基因,包括vg,vgr,sox 21b,testis-specific serine/threonine-protein kinase,spermatogenesis-associated protein和vasa-like gene等。【结论】本研究获得了意大利蝗转录组数据及性腺发育相关基因,为进一步研究意大利蝗生殖调控机理奠定了分子基础。     

野生矮牡丹高通量转录组测序分析

为探究野生稷山矮牡丹(Paeonia jishanensis)花青素等抗氧化活性物质生物合成的遗传基础,本文采用高通量测序技术Illumina NextSeqTM500对其幼嫩叶片、嫩茎和花芽混合池样本进行转录组测序分析。共获得39450个unigene,其中3563个unigene中有4106个SSR位点。在Swiss Prot数据库中注释到29420个unigene,其中11个与花青素生物合成过程相关,22个与类黄酮生物合成过程相关,15个与异黄酮生物合成过程相关;NR数据库中注释高达100%,有12个与类黄酮生物合成相关,7个与花青素生物合成相关;GO数据库中注释到24291个unigene,分为生物学过程、细胞组分及分子功能等三大类53个功能组,有21个unigene与黄酮、黄酮醇及异黄酮生物合成有关,4个unigene与花青素合成有关;eggNOG数据库中注释到38238个unigene,涉及植物生长发育过程绝大多数生命活动过程;KEGG数据库注释到5709个unigene,分别定位到6大类别42个代谢途径,其中25个unigene与花青素生物合成相关。研究表明,在稷山矮牡丹转录组序列中,共鉴定和发掘出多个涉及类黄酮化合物生物合成、花青素生物合成过程的相关基因序列及SSR位点,为下一步开展相关代谢合成途径研究奠定了基础。     

沙葱萤叶甲对蜕皮激素响应的转录组分析

为建立沙葱萤叶甲Gauleruca daurica成虫对蜕皮激素(20-hydroxyecdysone, 20E)响应的转录组数据库,挖掘对20E响应的基因以及代谢和信号通路,并在转录组水平探讨20E调控生殖滞育的分子机制,本研究采用Illumina HiSeqTM4000高通量测序平台对20E及二甲基亚砜(dimethyl sulfoxide, DMSO)处理后的沙葱萤叶甲成虫进行了转录组测序,共获得80 313个unigene;与阴性对照DMSO相比,共获得201个差异表达基因,其中106个上调、95个下调。GO和KEGG富集分析表明,多数差异表达功能基因富集于各种代谢通路,其中核黄素代谢(riboflavin metabolism)、溶酶体(lysosome)和泛酸与乙酰辅酶A合成(pantothenate and CoA biosynthesis)通路显著富集(q<0.05)。结果表明,20E可能通过影响多种代谢通路调控沙葱萤叶甲的生殖滞育。     

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

叉角厉蝽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个气味结合蛋白基因(Odor...