榆树叶片虫瘿形成过程中蛋白质组学分析

为探究榆树虫瘿叶片形成的分子机制,以不同时期榆瘿蚜取食诱导的榆树叶片为试材,利用iTRAQ技术分析榆树虫瘿叶片形成过程中蛋白表达丰度的差异。通过质谱鉴定,共得到2 689个蛋白,与KEEG数据库进行比对,发现2 145个蛋白被注释到126个不同的代谢通路中,和本研究相关的蛋白有12条,涉及代谢途径的蛋白数量最多,为813(37.9%)个。未被榆瘿蚜取食的叶片与榆瘿蚜取食诱导叶片形成虫瘿的前期、中期、后期相比,差异蛋白分别有418个、390个、244个,筛选出虫瘿发育的共有差异蛋白29个,这些蛋白在榆瘿蚜取食形成虫瘿的初始形成期和成长分化期持续发挥作用的有过氧化物酶、过氧化氢酶等氧化还原酶类;在前期持续发挥作用的有翻译控制肿瘤蛋白TCTP和肌动蛋白;在榆瘿蚜取食形成虫瘿的开裂期起到重要抗性作用的有溶质的抗坏血酸过氧化物酶、70 kD热激蛋白等。在榆树虫瘿叶片发展过程中,数个蛋白基因涉及了应激防御反应、氧化还原、免疫系统过程和光合作用等生理反应过程,建议进一步进行榆树虫瘿叶片形成的分子机制研究。     

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

本研究采用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值对基因表达量进行评估.本研究获得的转录组信息及分析结果,为进一步研究阿尔泰蝠蛾的基因功能及低温生态适应性奠定分子基础.     

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

红豆越橘是一种新兴水果,具有较高的营养及食疗药用价值。本研究利用高通量测序技术对野生红豆越橘的幼叶转录组进行了测序,并对测序数据进行生物信息学分析。总共获得大约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%。以上研究结果对了解红果越橘的生长发育、生理生化等分子机制,对进一步进行红豆越橘分子方面的相关研究均具有较强的参考价值。     

基于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位点。通过对大鲵转录组测序,获得了大量的转录组信息,有助于进一步研究大鲵功能基因克隆、基因组学、遗传多样性分析和分子标记开发等。     

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

为了增强对资源植物山莨菪的深入了解,本研究采用高通量测序技术对山莨菪进行转录组测序分析,经过处理得到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%。该结果丰富了山莨菪的转录组信息数据,为该物种分子生物学方面的研究奠定了基础,有助于进一步开展对山莨菪的合理保护及开发利用工作。     

基于高通量测序的麦红吸浆虫转录组分析

麦红吸浆虫Sitodiplosis mosellana(Gehin)是一种世界性的小麦害虫。为获得其转录组信息,本研究采用新一代高通量测序技术Illumina HiSeq TM 2000对麦红吸浆虫成虫转录组进行测序。共获得转录组样本数据量为27.88 G,经分析共获得59257个Unigenes,总长度49861164 bp,最短20 bp,最长29282 bp,平均长度841 bp。将Unigenes序列与NR、NT、Swiss-Prot、KEGG、GO和KOG数据库进行比对(e≤10-10),共获得95029个结果。通过GO功能分类,共有19584个Unigenes在GO数据库中细胞组分、分子功能和生物学过程等3大类50个功能组中找到对应。与KOG数据库进行比对,共有11279个麦红吸浆虫Unigenes被注释,按功能大致可分为26类。通过KEGG pathways分析,共有9110个麦红吸浆虫Unigenes被注释,分别归属于细胞进程、环境信息进程、遗传信息进程、新陈代谢和有机体系统5大类代谢途径,主要包括细胞生长与死亡、细胞运动、信号转导、能量代谢等32类代谢途径。CDS预测发现30088条序列可被编码,占全部基因的50.78%。SSR位点查找发现,在59257个Unigenes中共找到36323个SSR位点,发生率为61.30%。本研究获得的巨大的麦红吸浆虫转录组信息,为麦红吸浆虫的功能基因挖掘提供了重要的信息资源。     

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

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

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

采用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)位点。为进一步挖掘泡桐重要功能基因提供了大量数据。     

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

为探究建兰花色形成的分子调控机理,该研究以同株建兰黄绿色、红色花瓣为实验材料,采用高通量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连接酶基因表达上调,有利于类黄酮前体积累;查尔酮合成酶、二氢黄酮醇还原酶、花青素合成酶基因在黄绿色花瓣中几乎不表达,在红色花瓣中表达明显上调,可能与建兰花色形成相关。     

感染球孢白僵菌的烟粉虱若虫免疫应答转录组分析

【目的】筛选烟粉虱Bemisia tabaci(Gennadius)若虫应对白僵菌Beauveria bassiana侵染的应答基因,以进一步研究烟粉虱免疫反应的分子机制。【方法】采用新一代高通量测序技术对感染和非感染白僵菌的烟粉虱4龄若虫进行了测序分析,并筛选了差异表达基因;利用生物信息学工具对转录组测序得到的基因进行了功能注释、分类以及参与的信号通路展示。【结果】组装得到非冗余Unigene 232 554个,其N50和N90分别为1 153 bp和260 bp,平均长度为67 424 bp。对所有的基因进行差异表达分析发现:以P<0.05为标准筛选得到了1 166个差异表达基因,其中上调表达基因有474个,下调表达基因有692个,其中,与免疫相关的基因有405个;GO富集分析发现:有416个GO term有富集现象,包括156个生物学过程(66 402个Unigenes),89个细胞组分(27 645个Unigenes)和154个分子功能(73 417个Unigenes);KEGG代谢通路分析将烟粉虱转录组1 145个DEG匹配到309个通路上,其中76个通路得到了富集。【结论】对405个可能参与到烟粉虱若虫对白僵菌侵染的免疫识别和防御的基因进行了测序。研究结果为从分子水平上开展应用虫生真菌防治烟粉虱的研究奠定信息学基础。     

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位点,最多的为单核苷...     

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

为更好地挖掘八角(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多态位点。该转录组分析为今后长梗杜鹃乃至杜鹃属植物功能基因挖掘与利用、基因克隆、抗性机理分析、遗传资源分类和进化、分子标记开发以及分子辅助育种等研究提供了基础数据和重要参考。     

闽楠转录组分析及基因功能注释

采用第二代Illumina HiSeq测序技术对闽楠的木质部、韧皮部、叶片进行转录组测序，分别获得Clean Reads片段41 383 707条、43 343 922条、44 191 586条，经转录本拼接后得到序列总长度达120 535 288 bp的383 331条Conting片段，进一步组装得到平均长度为542 bp的151 729条Unigenes。将闽楠转录组Unigenes进行基因功能注释，与NR数据库比对发现，其与葡萄的相似序列最多（34%），与黄瓜、野草莓、大豆的同源性较低（各占3%）；进行GO功能注释，可将其划分为生物过程、细胞成分、分子功能3大类共计52个分支，与eggNOG数据库比对可分为25类，通过KEGG功能注释可知转录组中涉及的基因共参与了176条代谢通路，其中核糖体和碳代谢获得的注释较多。另外通过MISA软件分析，共获得35 972个SSR位点。其中，单核苷酸、二核苷酸和三核苷酸为优势重复类型，SSR位点数分别为21 762（60.50%），8 931（24.83%），4 924（13.69%）。闽楠转录组分析及基因功能注释为深入开展闽楠遗传育种及分子生物学相关研究奠定基础。     

绿色杜氏藻转录组分析

为了深入了解绿色杜氏藻(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条关键酶转录本(二羟丙酮激酶),该酶可能与绿色杜氏藻耐盐胁迫中甘油的合成有较大关系。本研究进一步完善了绿色杜氏藻的基因信息,为绿色杜氏藻代谢途径研究奠定了坚实的基础。     

基于转录组测序的铁皮石斛植物甾醇生物合成相关基因分析

为了解铁皮石斛(Dendrobium officinale)植物甾醇的生物合成途径,利用Illumina HiSeq 4000高通量测序技术对茎和叶进行转录组测序,比较植物甾醇生物合成关键酶基因的表达。结果表明,共获得43 085条Unigenes,其中24 459条在Nr、Swiss-prot、KOG和KEGG数据库获得注释,7 333条获得共同注释。KEGG代谢途径分析表明,铁皮石斛植物甾醇生物合成分为3个阶段,共有50个Unigenes (30种酶)参与。表达分析表明,DXR和HMED的表达量明显高于MK和MVD;成熟期茎、叶SMT1的表达量比生长期高,SMT2则生长期高于成熟期;同一时期,SMT1和SMT2在叶的表达量都比茎高。这为铁皮石斛植物甾醇的开发利用和调控植物甾醇生物合成研究提供了科学依据。     

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 HiSeq 2500平台进行全转录组测序,干净序列经组装后得单一基因(Unigene),将Unigene在非冗余蛋白/核酸数据库(nonredundant protein database,NR)、核酸序列数据库(nucleotide sequence database,NT)、基因本体论数据库(gene ontology,GO)、蛋白质真核同源数据库(clusters of eukaryotic orthologous groups,COG)、京都基因与基因组百科全书(Kyoto encyclopedia of genes and genomes,KEGG)、蛋白质序列数据库SwissProt和Interpro进行生物信息学分析。结果表明:共获得干净序列6.67 Gb,组装后得到58 646个Unigene,平均长度为1 023 bp,Unigene在上述数据库的总体注释率为69.25%。其中,在GO数据库中,20 762个基因被注释到生物功能、细胞组分和分子功能3个本体的52个功能组,COG注释了20 633个基因并将其划分为25个...