基于高通量测序的阿尔泰蝠蛾(鳞翅目:蝙蝠蛾科)幼虫转录组分析

本研究采用Illumina HiSeq TM 2500测序平台对阿尔泰蝠蛾Hepialus altaicola Wang幼虫进行转录组测序及生物信息学分析.经序列拼接后共获得100133个Unigenes,总长度86319112 bp,平均长度862 bp,N50长度1628 bp.将Unigenes与NR、COG/KOG、Pfam、Swiss-Prot、GO、KEGG数据库比对,共获得38198条Unigenes,其中Nr数据库注释的Unigenes最多,为32381条,占32.34％.通过GO功能分类,共有13216个Unigenes在GO数据库中细胞组分、分子功能和生物学过程等3大类57个分支中找到注释;KEGG通路分析,共有15058条Unigenes被注释,归属于305条代谢通路.CDS预测发现54002条序列可被编码,占全部基因的53.93％.基因注释进一步获得311个与冷适应相关的代谢调节基因,并用FPKM值对基因表达量进行评估.本研究获得的转录组信息及分析结果,为进一步研究阿尔泰蝠蛾的基因功能及低温生态适应性奠定分子基础.     

在小麦和叶锈菌互作过程中TCTP的表达特征

翻译控制肿瘤蛋白是调节植物生长发育、参与植物抵抗逆境胁迫的重要调节因子。本研究以叶锈菌生理小种260与4种含有不同抗叶锈基因的近等基因系分别组成亲和程度不同的组合为对象,采用RT-qPCR和Western Blotting方法,从转录和翻译两个水平探讨TCTP在小麦与叶锈菌互作过程中的表达模式。研究结果表明,无论是转录水平还是翻译水平上,TCTP在不亲和组合中均表现为表达量显著上调,但在不同的Lr基因背景下其表达模式有所不同;而在亲和组合中并未表现出明显的变化趋势。基于TCTP mRNA和蛋白水平的结果,发现在小麦与叶锈菌互作体系中TCTP在转录水平和翻译水平均受到调控,推测其在小麦抵抗叶锈菌侵染的防卫反应中可能发挥正调控作用,这将为进一步研究TCTP的功能及其作用机制奠定了基础。     

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

采用新一代高通量测序技术平台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。这些研究极大地扩充了云南松的基因资源,将有助于云南松基因的发掘与利用、分子标记的开发及其种质资源遗传改良的研究等。     

基于RNA-Seq高通量测序技术的大鲵转录组分析

为发掘大鲵的功能基因,该研究以大鲵组织作为研究对象,利用RNA-seq技术对大鲵进行转录本测序和数据分析,经拼接组装共获得132 912条Unigenes,序列平均长度690 bp,N50为1 263 bp。另外从长度分布、GC含量与表达水平等方面对Unigenes进行评估,数据显示测序质量好,可信度高。此外,本研究也预测出132 416个能编码蛋白的Unigenes。使用9大数据库CDD、KOG、NR、NT、PFAM、Swiss-Prot、Tr EMBL、GO和KEGG注释大鲵转录组Unigenes,分别对应有24 049、18 406、36 711、15 858、20 500、27 515、36 705、28 879和10 958条Unigenes获得注释。其中,6 323条Unigenes在以上所有数据库中同时注释成功,39 672条Unigenes至少被一个数据库注释。KEGG分析结果显示:获得注释的10 958条Unigenes被划分到343个代谢通路中,参与信号转导通路的Unigenes数最多,共有2 644条(24.12%);其次是免疫系统通路的Unigenes有2 450条(15.29%)。最后本研究还检测到了29 790个SSR位点。通过对大鲵转录组测序,获得了大量的转录组信息,有助于进一步研究大鲵功能基因克隆、基因组学、遗传多样性分析和分子标记开发等。     

基于RNA-seq技术的楮头红转录组分析

利用RNA-seq技术对所构建的楮头红叶片的转录组进行测定,对原始reads进行过滤和组装,得到了51 305条质量较高的Unigenes,平均长度为921 nt,N50为1 490 nt。利用BLAST和BLAST2GO软件对这些从头组装的Unigenes进行注释。用NCBI蛋白质数据库(Nr)、非冗余核苷酸数据库(Nt)、基因本体论(GO)、直系同源基因簇(COG)和京都基因与基因组百科全书(KEGG)数据库做参考,共注释了40 532条Unigenes。注释到Nr、Nt、Swiss-Prot、KEGG、COG和GO库中的比例相对较高,分别为77.53%、56.18%、53.14%、46.58%、29.69%和60.72%。在蛋白质数据库中对所有的Unigenes进行blast以后,发现有39 302个CDS,用ESTscan预测了2 065个CDS。KEGG通路分析显示,参与次生代谢物生物合成的Unigenes有2 323条,占全部Unigenes的9.72%。其中有78条Unigenes编码了细胞色素P450家族蛋白,这些信息为药用植物次生代谢物生物合成关键基因的挖掘提供了理论参考。     

建兰花色形成的分子调控机理研究

为探究建兰花色形成的分子调控机理,该研究以同株建兰黄绿色、红色花瓣为实验材料,采用高通量Highseq测序技术进行转录组文库构建,从基因水平探索花色物质的合成代谢通路及关键调控基因的转录活性。结果表明:(1)转录组测序共获得131 110 030条过滤序列(Clean Reads)、106 479条单基因(Unigenes),通过NR、GO、COG、KEGG等公共数据库比对,获得了29 748个有注释信息的Unigenes。(2)建兰黄绿色花瓣与红色花瓣中包括20个类黄酮合成代谢相关差异表达基因,与黄绿色建兰花瓣相比,红色花瓣中776个基因转录表达量上升,589个基因转录表达量下降,在KEGG数据库被注释到93条代谢通路,涉及6条类黄酮合成代谢相关的途径,共20个差异表达基因(83Unigenes)。(3)QRT-PCR分析显示,所选20个基因在建兰2种颜色花瓣中的表达量比值趋势与转录组FPKM比值趋势一致,表明该研究获得的转录组数据具有较高的参考价值;其中分支酸、苯丙氨酸解氨酶、肉桂酸4-羟化酶、4-香豆酸:辅酶A连接酶基因表达上调,有利于类黄酮前体积累;查尔酮合成酶、二氢黄酮醇还原酶、花青素合成酶基因在黄绿色花瓣中几乎不表达,在红色花瓣中表达明显上调,可能与建兰花色形成相关。     

泡桐维管形成层及木质部区的转录组分析

采用Illumina HiSeqTM2500高通量测序技术,获得泡桐维管形成层及木质部区组织转录组的109021918条clean reads(16.35 Gb)。将clean reads从头组装得到104432个单基因簇(Unigene),平均长度662 nt。将组装得到的Unigenes与公共数据库进行序列比对,分别有40789(Nr:39.05%)、31675(NT:30.33%)、15539(COG:14.87%)、29168(GO:27.93%)、16316(KEGG:15.62%)、30499(SwissProt:29.20%)以及28828(Pfam:27.6%)个Unigenes获得功能注释。通过与GO数据库的比对分析,注释的29168个Unigenes归于生物过程、细胞组分及分子功能三大类的55个功能组;15539条Unigenes注释到COG数据库,被分为25个类别中;基于KEGG数据库可将16316个Unigenes归于130个代谢途径。此外,在泡桐的维管形成层及木质部区转录组中共检测出16118个简单序列重复(SSR)位点。为进一步挖掘泡桐重要功能基因提供了大量数据。     

NaCl胁迫下黑果枸杞转录组测序分析

研究黑果枸杞在不同浓度盐胁迫下基因表达谱变化情况,为进一步研究黑果枸杞抗盐分子机制奠定研究基础。对0(CK)、50、250 mmol/L NaCl溶液胁迫的黑果枸杞组培苗的根和叶在胁迫时间为0、1、12 h时分别取样,采用转录组测序(RNASeq)技术进行测序分析。结果表明,转录组测序共产生222.49 Gb原始数据,拼接出Unigenes 86 037条,注释到7大功能数据库(GO、KEGG、KOG、NR、Pfam、Swiss-Prot和egg NOG)上的Unigenes总数为46 594个,占总Unigenes的54.76%,还有38 929个Unigenes在这些数据库中没有得到注释。通过GO分类和KEGG Pathway富集性分析,分别归于51个GO类别和211条代谢途径。差异表达基因分析显示,黑果枸杞叶片和根的上调基因和下调基因数随着NaCl浓度的增大和处理时间的延长均呈增加趋势,叶片中的上调基因数(7 514)小于下调基因数(9 032),根中的上调基因数(12 347)大于下调基因数(11 559)。在黑果枸杞盐胁迫下转录组中发现28 325个SSR位点,最多的为单核苷酸SSR,占70.47%。综合分析表明,黑果枸杞对盐胁迫的反应是一个多基因参与、多个生物过程协同调控的过程,基因表达量的变化可能是基因调控的主要方式。     

红豆越橘幼叶转录组高通量测序及生物信息学分析

红豆越橘是一种新兴水果,具有较高的营养及食疗药用价值。本研究利用高通量测序技术对野生红豆越橘的幼叶转录组进行了测序,并对测序数据进行生物信息学分析。总共获得大约2 G的纯净数据,拼装了48 736条Unigenes。GO数据库注释到的Unigenes分别涉及到生物学过程、细胞成分及分子功能相关的共43种生理代谢功能;18 728条Unigenes能被KOG数据库注释,共涉及25条代谢通路;8 321条Unigenes能被KEGG数据库成功注释,共涉及到5个功能大类、20个功能中类、123条代谢通路;30 467条Unigenes可被Nr数据库注释;21 551条Unigenes可被Swissprot数据库注释;以上4个数据库共注释到30 545条Unigenes,占全部Unigenes的62.67%;被以上4个数据库均注释到的Unigene为6 075条,占全部Unigenes的12.47%。同时,生物信息学分析还显示:全部Unigenes中,有5 197条Unigenes编码转录因子,涉及57个家族;2 747条Unigenes编码抗性基因,涉及19个家族。共检测到8 099个SSR位点,其中2碱基重复的SSR位点达5 791个,占SSR位点总数的71.52%。以上研究结果对了解红果越橘的生长发育、生理生化等分子机制,对进一步进行红豆越橘分子方面的相关研究均具有较强的参考价值。     

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

该研究采用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基因与滇山茶的花叶形成密切相关。研究表明,滇山茶的花叶呈色是由于叶色相关代谢通路中大量的基因表达量发生变化引起的,并不是由单基因突变所致,病毒诱导的多个代谢途径基因沉默可能是滇山茶花叶呈色的主要原因。该研究丰富了滇山茶的转录组信息,并为滇山茶的彩叶品种培育提供了一定的参考。     

绿色杜氏藻转录组分析

为了深入了解绿色杜氏藻(Dunaliella viridis)基因信息及功能、耐盐相关通路(甘油脂代谢)及关键酶,本文首次通过Illumina HiSeq^TM 2000高通量测序技术对绿色杜氏藻转录组进行测序,利用Trinity软件将数据组装形成转录本,对所有转录本进行COG(Clusters of Orthologous Groups)、GO(Gene Ontology)和KEGG(Kyoto Encyclopedia of Genes and Genomes)分类和功能注释、Pathway注释以及蛋白编码区(Opening reading fragment,ORF)的预测,并对甘油脂代谢通路关键酶基因进行了分析。转录组测序共获得81 593个转录本,其中ORF共有77 117条,约占所有转录本的94.50%。COG分类结果表明,16 569条转录本被分为24个类别。GO分类结果表明,76 436条转录本被注释。在所有注释分类中,生物学过程转录本数量最多,为30 678条,占总转录本数的40.14%。KEGG分析结果表明,317个标准途径中包含26 428条转录本,含转录本最多的类别是代谢,为9949条(37.65%)。与代谢有关的途径为131条,占所有注释途径的41.32%。在甘油脂代谢通路中仅发现1条关键酶转录本(二羟丙酮激酶),该酶可能与绿色杜氏藻耐盐胁迫中甘油的合成有较大关系。本研究进一步完善了绿色杜氏藻的基因信息,为绿色杜氏藻代谢途径研究奠定了坚实的基础。     

基于转录组的八角挥发油合成相关基因挖掘与分析

为更好地挖掘八角(Illicium verum)挥发油合成相关基因，该文对挥发油性状差异显著的优良无性系桂角69号及普通品种砧01号叶片进行了转录组测序及组装注释，并对差异表达基因进行了GO分类和KEGG通路分析。结果表明：(1)转录本经组装后获得84 182条序列，使用NR、NT、Swiss-Prot、KEGG、KOG、GO和Pfam数据库进行序列比对，共注释了59 161条序列，筛选出30 572个差异表达基因。与砧01号相比，桂角69号叶片中上调基因有15 025个，下调基因有15 547个。(2)GO分类结果显示共有20 287个差异基因被注释。KEGG分析结果表明，有21 600个差异基因被注释到133条KEGG通路上，其中挥发油合成相关的单萜生物合成通路、萜类骨架生物合成通路、苯丙素合成通路中的芳樟醇合酶、月桂烯合酶、香叶基香叶基焦磷酸合酶、肉桂酰辅酶A还原酶、咖啡酸3-O-甲基转移酶、肉桂醇脱氢酶等关键酶基因呈差异表达。(3)转录因子分析发现差异表达基因分布于31个转录因子家族，其中MYB家族序列数量最多。该文利用转录组测序技术分析八角优良无性系与普通品种叶片的差异基因及...     

基于高通量测序的极小种群野生植物长梗杜鹃转录组分析

为加强特有濒危植物长梗杜鹃资源的评价、保护和鉴定工作,及为今后遗传育种和改善其农艺性状提供有益参考。本研究采用新一代高通量测序平台Illumina HiSeq 4000对其转录组测序,得到的数据过滤后进行de novo组装并聚类去冗余,获得74 092个Unigenes,平均长度、N₅₀、Q₂₀、Q₃₀以及GC含量分别为938 nt、1 616 nt、98.22%、95.20%和43.24%,其中1 Kb以上的Unigenes有23 879条。通过与七大功能数据库比对,分别有39 876（NR：53.82%）、38 065（NT：51.38%）、27 384（Swissprot：36.96%）、16 099（COG：21.73%）、30 401（KEGG：41.03%）、17 518（GO：23.64%）以及29 676（Interpro：40.05%）条Unigenes获得功能注释。长梗杜鹃转录组中的Unigenes根据GO功能大致可分为生物学过程、细胞组分和分子功能3大类56亚类,其中执行生物学过程的基因最多,占41.53%;与COG数据库比对,根据其功能大致可分为25类;以KEGG数据库为参考,可归为6个代谢通路大类、32类代谢途径,并发掘出176条与人类疾病相关的Unigenes,包括内分泌及代谢疾病（167条）和抗药性（9条）。根据注释结果共检测出39 418个CDS,未注释上的Unigenes使用ESTScan预测后获得3 194个CDS。同时,预出到1 488个编码TF的Unigenes,以及检测到57 927个SNP多态位点。该转录组分析为今后长梗杜鹃乃至杜鹃属植物功能基因挖掘与利用、基因克隆、抗性机理分析、遗传资源分类和进化、分子标记开发以及分子辅助育种等研究提供了基础数据和重要参考。     

刺猎蝽属种间转录组比较分析

【目的】建立刺猎蝽属Sclomina 3种的转录组数据库,分析近缘种间转录组表达差异,探讨在转录组水平的种间趋异情况。【方法】采用Illumina HiSeq^TM4000高通量测序平台进行刺猎蝽属齿缘刺猎蝽S.erinacea、兴仁刺猎蝽S.xingrensis和广西刺猎蝽S.guangxiensis 3种转录组测序;利用Nr, Swiss-Prot, COG/KOG和KEGG数据库进行基因功能注释;对种间差异表达基因(differentially expressed genes, DEGs)进行GO功能注释和KEGG通路富集分析。【结果】齿缘刺猎蝽、兴仁刺猎蝽和广西刺猎蝽转录组测序组装获得42 215个unigenes。21 117个基因在上述4个数据库中得到注释(各数据库注释基因数为：Nr, 20 522; Swiss-Prot, 15 550; COG/KOG, 13 969; KEGG,10 850)。兴仁刺猎蝽vs齿缘刺猎蝽转录组有3 390个DEGs (803个上调,2 587个下调),兴仁刺猎蝽vs广西刺猎蝽转录组有12 543个DEGs (6 63...     

小麦条锈菌PsCdc2基因的克隆及在条锈菌侵染小麦后的转录表达分析

【目的】克隆小麦条锈菌细胞分裂基因PsCdc2,分析该基因在条锈菌接种小麦后不同时间点的表达特征。【方法】利用PCR和RT-PCR技术克隆PsCdc2的cDNA序列和基因组序列,采用生物信息学技术预测分析该基因编码蛋白的保守结构域及基本特性,对该蛋白进行系统发育分析,构建进化树;运用实时荧光定量RT-PCR技术,以PsCdc2在夏孢子时期的表达情况为对照,分析该基因在亲和及非亲和互作中不同时间点的表达特征。【结果】PsCdc2基因组序列长2279 bp,由11个外显子和10个内含子构成,开放阅读框为885 bp,编码294个氨基酸,分子量为33.14 kDa,等电点为6.26。编码蛋白含两个保守的激酶特征位点,一个跨膜螺旋区域。PsCdc2基因编码蛋白与小麦秆锈菌、新型隐球菌、玉米瘤黑粉菌等多种真菌的Cdc2高度相似,其中与小麦秆锈菌的Cdc2亲缘关系最近,序列相似性达73.1%。实时荧光定量RT-PCR结果表明,在亲和组合中,该基因在条锈菌接种小麦的前期上调表达,其中接种后12 h时表达量最高,约为夏孢子中表达量的1.62倍,接种后24-268 h,基因表达基本呈下调趋势,其中96 h基因表达量最低,仅为夏孢子时期的0.07倍。在非亲和组合中,该基因表达基本呈下调趋势,在接种后各个时间点的表达量均低于在夏孢子中的表达量,其中接种后12 h时表达量最高,但仅为夏孢子中表达量的0.34倍;接种后96 h表达量最低,为夏孢子中表达量的0.02倍。【结论】PsCdc2可能通过调控条锈菌的细胞周期循环参与了侵染前期初生菌丝生长和吸器母细胞的形成,与条锈菌的致病性相关。本文首次报道了小麦条锈菌的Cdc2基因,为进一步揭示条锈菌细胞周期调控的本质及研究开发靶向Cdc2的新型农药,以及实现对小麦条锈病的新型药剂防治提供了理论基础。     

基于RNA-Seq技术的蜡蚧轮枝菌转录组分析（英文）

【目的】蜡蚧轮枝菌是一种防治植物病原物和害虫的生防菌,弄清其致病的分子基础具有重要的意义。【方法】利用Illumina HiS eq技术测序蜡蚧轮枝菌的cD NA文库,并进行RNA-Seq序列拼接和功能注释。【结果】共获得1634670条高品质reads,碱基GC含量48.50%。14856条Unigene经组装拼接后,在Nr、Swiss-Prot、COG和KEGG数据库中分别比对得到1467、9379、6518和4188条Unigenes。采用GO分类工具将21403条Unigenes分成24个功能组,主要包含在细胞组成(1369)、分子功能(1679)和生物学过程(1928)3个大类里,在COG数据库注释的基础上,把6518条Unigenes分成38个功能组,通过GO和COG丰度识别分类,聚类的主要是一般功能预测、次生代谢产物的合成、运输和分解代谢以及信号转导机制起作用的独立基因。能与KEGG路径匹配的108条Unigenes中大部分和新陈代谢途径及次生代谢产物的合成有关。【结论】这些结果将有助于找寻蜡蚧轮枝菌的致病因子,从而丰富其致病机理。     

发菜RNA-seq转录组分析及耐旱相关基因筛选

发菜是一种陆生耐旱固氮蓝藻,具有重要的生态和资源价值,但目前发菜的基因组和转录组信息还比较缺乏。本研究采用Illumina Hi Seq TM 2500高通量测序平台对充分吸水和失水后的发菜藻体进行转录组分析,经过拼接组装共获得30 013 560条clean reads,3 100个Unigene,序列平均长度854 nt。将Unigene序列与NR、NT、Swiss-Port、PFAM、GO、COG、KEGG(E-value1.0E~(-5))等数据库进行比对,共有2 874个Unigene获得了注释。通过GO功能分类,917个Unigene映射到GO不同功能节点上;通过Unigene与COG数据库的比对,可将有COG功能的514个Unigene分为20类;通过KEGG pathway分析,有1 220个编码基因参与了135条已知的代谢通路。提示谷胱甘肽、类胡萝卜素、淀粉和蔗糖、ABC转运体等代谢途径的部分基因在失水藻体上调表达,说明这些通路与发菜耐旱性关系密切。因此这些注释信息为今后发菜耐旱基因的挖掘提供了重要依据,同时也为发菜分子生物学的进一步研究提供了试验基础。     

Transcriptome analysis of Verticillium lecanii based on RNA-Seq technology

[目的]蜡蚧轮枝菌是一种防治植物病原物和害虫的生防菌,弄清其致病的分子基础具有重要的意义。[方法]利用Illumina HiSeq技术测序蜡蚧轮枝菌的cDNA文库,并进行RNA-Seq序列拼接和功能注释。[结果]共获得1634670条高品质reads,碱基GC含量48.50%。14856条Unigene经组装拼接后,在Nr、Swiss-Prot、COG和KEGG数据库中分别比对得到1467、9379、6518和4188条Unigenes。采用GO分类工具将21403条Unigenes分成24个功能组,主要包含在细胞组成（1369）、分子功能（1679）和生物学过程（1928）3个大类里,在COG数据库注释的基础上,把6518条Unigenes分成38个功能组,通过GO和COG丰度识别分类,聚类的主要是一般功能预测、次生代谢产物的合成、运输和分解代谢以及信号转导机制起作用的独立基因。能与KEGG路径匹配的108条Unigenes中大部分和新陈代谢途径及次生代谢产物的合成有关。[结论]这些结果将有助于找寻蜡蚧轮枝菌的致病因子,从而丰富其致病机理。     

钩苞大丁草高通量转录组测序及差异表达分析

该研究采用Illumina Hi-Seq2500高通量测序技术,对叶背有纤维和无纤维发育的两组钩苞大丁草(Gerbera delavayi Franch.)叶片样品的cDNA进行转录组测序,分析其叶背毡毛纤维发育机理。测序结果得到了108 694条单基因序列,进一步筛选得到了1 605条差异表达基因,838条差异表达基因在GO数据库具有功能定义,512条差异表达基因在COG分类体系中具有详细的蛋白功能释义,315条差异表达基因注释到了KEGG数据库中。其中,氨基糖和核苷酸糖代谢途径中控制纤维素合成的相关基因,苯丙烷类生物合成途径中控制木质素合成的相关基因,以及激素信号转导途径中控制生长素信号转导的相关基因表达量下调;激素信号转导途径中控制细胞分裂素、脱落酸信号转导的相关基因表达量上调。研究结果在一定程度上丰富了钩苞大丁草的基因信息,并为后续的遗传改良提供了基础数据。     

榆树叶片形成虫瘿过程中转录组学分析

榆瘿蚜取食侵染榆树叶片形成了榆树虫瘿,本研究采用新一代的高通量Illumina Hi Seq^TM 2000技术测序平台对榆瘿蚜取食刺激的榆树叶片进行转录组测序和功能注释,利用生物学方法对基因表达和功能进行研究。测序获得23.19 Gb碱基序列信息,通过对测序数据进行序列过滤、拼接和去冗余,共获得102 017个Unigenes,通过NR与BLAST等数据库比对,其中有37 899个(37.15%) Unigense被注释。利用KOG、GO、KEGG等数据库对榆树虫瘿叶片的Unigense进行比对,按功其能将匹配的Unigenes基因划分25大类;GO注释将信息归纳为基因的3大主类,57个亚类;以KEGG数据库为参考,将Unigene定位到110个不同的代谢通路,包括氧化应激防御、植物激素信号转导、碳水化合物以及次生物代谢等代谢相关的Unigenes通路。本研究通过二代高通量转录组测序技术研究榆瘿蚜侵染下榆树虫瘿的相关基因,为今后研究榆瘿蚜侵染榆树叶片形成虫瘿的分子机理提供了基础资料。