NaCl胁迫下柳树的转录组测序与初步分析

从组学水平分析盐胁迫下柳树内在分子机制,为柳树耐盐研究及耐盐基因的挖掘利用提供理论依据。本研究通过转录组测序技术对‘盐柳1号’(Salix psammophila’Yanliu No.1’)和‘渤海柳1号’(Salix matsudana’Bohailiu No.1’)经150 mmol/L NaCl胁迫处理后的叶片和正常叶片(对照)进行高通量转录组测序,并对获得的unigene进行从头组装和注释分析。结果表明:转录组测序共获得183987条Unigenes,平均长度为1080.18 bp,分别有120130条、140813条、98066条、88425条、47982条、83732条Unigenes被注释到NT、NR、COG、SwissProt、KEGG、COG和GO数据库,共有149864条(81.45%)Unigenes得到注释。在GO功能注释中,共得到55个GO功能小类,在KEGG代谢通路分析时,获得了135条KEGG通路。该转录组测序数据质量高,结果覆盖面广,为柳树耐盐基因挖掘和研究提供了一定的理论参考。     

滇白前种子转录组测序及功能注释

采用高通量测序技术平台Illumina Novaseq 6000对滇白前种子进行转录组测序,共获得平均长度为753.9 bp的43 663个Unigenes。注释到六大功能数据库(NR,COG,Pfam, Swiss-Prot, GO, KEGG)上的Unigene总数为29 177个。滇白前Unigenes比对到NR数据库共有26 688条,与甜菜、藜麦、栓皮栎、菠菜有较高同源性;25 657条Unigenes在COG数据库得到26 500个注释,分为23类;19 942条Unigenes在GO数据库得到79 481个注释,按功能分为细胞组分、分子功能和生物过程三大类,分别有13、16、23个亚类,其中参与的生物过程较多;11 527条Unigene富集在KEGG数据库的101条代谢通路中,代谢相关的通路占比最大。同时,有527个Unigenes被注释为转录因子;有3 995条SSR标记被挖掘。利用高通量测序获得滇白前转录组信息,有助于从分子水平对滇白前进行深入研究。     

虎眼万年青鳞茎转录组特性分析

虎眼万年青具有很高的药用价值和经济价值,鳞茎是其主要的药用和繁殖器官,然而虎眼万年青鳞茎的分子生物学信息严重缺乏成为制约其研究生产的瓶颈。为获得虎眼万年青鳞茎的转录组信息,应用高通量测序技术对虎眼万年青鳞茎进行测序分析。采用His4000测序平台共获得8.2 G原始数据,经过整理分析后,获得7.6 G有效数据和121 223条Unigene;经过与NR、GO、COG及KEGG等数据库进行比较分析后,56 374条Unigene获得NR数据库的注释,39 263条Unigene获得GO数据库的注释,12 388条Unigene获得COG功能注释,5 405条Unigene参与了新陈代谢途径,2 477条Unigene参与了次生代谢途径,166条Unigene参与了萜类的合成,195条Unigene参与了生物碱的合成,295条Unigene作为转录因子参与新陈代谢过程调控。研究结果为采用生物技术手段提高虎眼万年青的药用成分的含量及改善其园艺性状提供了分子数据。     

白羊草转录组和差异性表达分析

[目的]获得白羊草(Bothriochloa ischaemum)转录组数据库和差异表达基因。[方法]选取对照组和干旱胁迫组作为受试材料,采用Illumina HiSeq 4000进行转录组测序,并进行系统的生物信息学分析。[结果]共获得25.23 Gb数据(118580条Unigenes)。将所得到的Unigenes与基因库数据进行比较,对Unigenes进行功能注释,共50070条Unigenes获得功能注释。白羊草转录组中的Unigene根据GO功能大致可分为细胞组分、分子功能和生物学过程三大类共54个分支。将Unigene与KOG数据库进行比对,根据其功能大致可分为25类。KEGG数据库作为参考,依据代谢途径可将Unigene定位到116个代谢途径分支。通过差异性分析发现,差异性表达基因共3987条,其中上调基因占39.45%,下调基因占60.55%,相差较大。[结论]本研究首次对白羊草转录组进行了分析,为白羊草分子生物学研究提供了宝贵的基因组数据库资源。     

基于高通量测序的南酸枣转录组分析

以南酸枣雌株叶片和成熟果实为材料,利用Illumina Hi Seq TM 4 000测序平台对其转录组进行测序,共获得14.15 G有效数据,通过序列拼接组装得到46 936条Unigene,平均长度为1 287 bp。36 299条Unigene(77.33%)在7大数据库(NR,NT,KO,Swiss-Prot,PFAM,GO和KOG)中得到注释,其中与GO数据库比对上的26 014条Unigene可分为生物过程、细胞组分和分子功能3大类53分支,注释到KOG数据库中的8 459条Unigene依据功能可分为25类。与KEGG数据库比对,6 225条Unigene分属5大类246条代谢通路中。南酸枣叶片和果实转录组中共发现5 631个差异表达基因,包括1 930个上调基因和3 701个下调基因。本研究获得的转录组数据将有助于开展南酸枣功能基因挖掘与利用、分子辅助育种和其种质资源遗传改良等方面的研究。     

盐胁迫下不同甘薯品种的转录组测序分析

为了获得甘薯耐盐转录组序列信息,挖掘差异表达基因及其相关代谢途径,本文以盐胁迫处理0 d、3 d和6d的耐盐甘薯品种‘榕薯819’以及不耐盐甘薯品种‘榕薯910’的叶片为材料,借助高通量测序技术进行转录组测序分析。结果表明, 2个品种共获得157,252条Unigenes,平均组装长度为576 bp。其中有83,264条Unigenes在七大数据库中得到注释,占总数的52.95%。NR注释分类结果显示,在牵牛花(Ipomoea nil)中比对到同源序列的Unigenes最多,共43,620条,占总数的57.05%。Unigenes在KOG数据库中的注释主要富集在普通功能预测(8752个)、信号转导机制(5067个)以及翻译后修饰、蛋白转换、分子伴侣(4471个)中。差异表达分析显示,在‘榕薯819’中,盐处理3 d和6d的样品差异表达基因数分别为323个和3752个,共参与了33个GO功能分类项和302条KEGG代谢通路。在‘榕薯910’中,差异表达基因数则分别为5554个和7395个,共参与了50个GO功能分类项,涉及了329条KEGG代谢通路。以部分差异表达基因的转录组数据为基础,...     

外源因子处理下蓝莓不同发育阶段果实的转录组分析

采用Illumina测序技术对在醋酸钙、硫酸铵和蔗糖处理后蓝莓不同发育阶段的果实进行转录组测序,获得Clean Reads 2723731442条,经组装得到平均长度为753.65 nt的87608条Unigene。将转录组Unigene进行基因功能注释,其中39867条Unigene能被NR数据库注释,与葡萄同源序列最多,占8.58%;与GO数据库比对发现,有29661条Unigene获得注释,分别匹配到生物过程、细胞组成和分子功能三大类共59个分支;与KOG数据库进行比对,发现有21992条Unigene具有功能信息,分别涉及25类;根据KEGG数据库的注释信息进行Pathway注释,参与的代谢通路共有246条;共检测到8704个SSR位点,其中双碱基重复的SSR占78.57%。本研究为探索外源物质调控蓝莓果实生长发育、生理代谢的分子机理提供了理论基础。     

基于高通量测序的大花序桉顶芽转录组生物信息学分析

为了获得珍贵用材树种大花序桉顶芽转录组数据及预测关键基因功能,本研究基于Illumina HiSeq X Ten测序技术获得大花序桉顶芽转录组原始数据,经Trinity组装拼接获得高质量Unigene,并与NR、Swiss-Prot、GO、KOG、egg NOG和KEGG等生物信息数据库进行序列比对和功能注释,利用MISA软件进行SSR位点搜索和分析。从大花序桉顶芽中共获得26 587条高质量Unigene,平均长度为1 279.69 bp;共有22 099条Unigene至少在一个数据库中被成功注释,其中,11 507条Unigene被注释到KOG数据库中25个功能类别,以参与一般功能基因的数量最多;GO数据库中,所注释到的14 105条Unigene分别匹配到生物功能、细胞组分和分子功能3大类50个功能基因区,其中执行生物过程所占比例最多;KEGG功能注释共发现有7 117个Unigene参与127条代谢通路,以代谢相关的基因最丰富;共有1 021条Unigene注释到转录因子数据库,分布于65个家族,其中比例最大的是bHLH和MYB家族;3 274条Unigene注释到植物抗性基因数据库,分布于13个类别,相匹配基因数量最大的是RLP和TNL。MISA软件共检测到12 366个SSR位点,分布密度为1/2.75 kb,重复基元类型丰富,标记开发潜力大。本研究利用高通量测序获得丰富的顶芽转录组信息,可以为大花序桉分子辅助育种提供丰富的资源。     

RNA-seq揭示碱性盐(NaHCO3)对细叶百合鳞茎基因表达的影响

本研究旨在揭示细叶百合在盐碱逆境中的基因表达表达情况,为合理利用盐碱地和大面积种植百合提供理论基础。以一年生的细叶百合鳞茎为材料,经过20 mmol/L NaHCO₃处理24 h后,用Illumina HiSeqTM2000测序平台进行转录组测序。共测得56828个mRNA的注释信息,其中,55433个Unigene被注释到Nr数据库,26973条Unigene被注释到KOG数据库,23610条Unigene被注释到GO数据库,13142条Unigene被注释到KEGG数据库的五大类中。共鉴定出390个差异表达基因。选取了9个与盐碱胁迫密切相关的基因进行qPCR验证,9个基因的表达结果与转录组的结果基本趋于一致。通过细叶百合在碳酸盐(NaHCO₃)逆境下的转录组对比数据,提供了许多的基因表达通路和表达量的差异,为合理利用盐碱地和大面积种植百合提供理论基础。     

百合未授粉成熟雌蕊转录组测序研究

本研究旨在大量发掘百合雌蕊特异表达基因,为研究百合杂交障碍的分子机制提供依据。利用Illumina高通量测序平台,进行百合未授粉成熟雌蕊转录组测序和数据组装,并对获得的Unigene进行功能注释、分类和代谢通路分析。“干柱头”雌蕊LoP1和“湿柱头”雌蕊LoP2分别获得40792条和39708条Unigene。2个样品的Unigene通过序列拼接,去冗余处理,以及同源转录本聚类,得到44175条All-Unigene。基于NR、NT、Swiss-Prot、GO、KEGG和COG等6个数据库进行相似性比对（E值≤10-5）,依次分别有30343、21647、21281、19274、12964和22227个基因被注释。KEGG分析中,19274条AllUnigene被注释到128个代谢通路上。基因分析显示,百合未授粉成熟雌蕊已为授粉受精做好了物质、能量、信号转导以及抗病原等准备工作。     

干旱胁迫下枇杷叶片的转录组分析

分析干旱胁迫下枇杷叶片的转录组,挖掘功能基因并对其差异表达基因进行筛选和分析,为枇杷抗旱提供理论依据。利用新一代高通量测序技术测序,对测序结果进行de novo拼接、功能注释和ORF预测,将差异表达基因在COG、GO和KEGG数据库中进行比对注释。测序结果表明,获得转录本共88 530个,平均长度为740.64 bp,ORF41 748条。COG、GO和KEGG数据库将转录本分别划分为24,54个功能类别及291条代谢通路中。25 197个差异表达的基因在30条代谢通路中显著富集,与酶活性、激素合成代谢和信号转导等相关的差异基因积极响应枇杷干旱,其中双萜类、油菜素类固醇和类胡萝卜素3条生物合成途径中的差异基因呈现出较为一致的表达。采用实时荧光定量(qRT-RCR)对选取的差异基因进行验证,其中过氧化物酶、蛋白激酶byr2、丝氨酸苏氨酸蛋白激酶、脱落酸8'-羟化酶、吲哚-3-乙酰乙酸合酶相关基因上调表达,细胞色素P450 734A1基因下调表达。为枇杷提供了较为全面的基因信息和代谢途径数据,为干旱胁迫下枇杷分子调控机制的深入研究奠定了基础。     

碰碰香叶片的转录组测序分析

为探讨碰碰香(Plectranthus tomentosa)叶片香气的分子机制,本研究对其进行了转录组测序分析,挖掘挥发性物质合成的关键基因,为植物次生代谢和花卉的分子育种提供借鉴和参考。分析结果共得到了4.2 G有效数据,通过序列拼接得到56 883条Unigene。其中,有29 617条Unigene得到GO注释,18 704条Unigene得到COG注释。KEGG pathways分析结果表明,共有1 446条Unigene注释到新陈代谢途径,625条Unigene注释到次生代谢生物合成途径,其中参与单萜、二萜、倍半萜和萜类骨架合成的Unigene有59条,生物碱类29条,至少369条Unigene参与了甲苯、松香油、酯类和烷烃等挥发性物质的代谢途径。本研究为探讨碰碰香挥发性物质的分子机制,挖掘重要功能基因,及其研究开发应用提供了基础数据。     

猕猴桃抗溃疡病转录组分析和基因功能注释

为了发掘中国野生猕猴桃资源抗溃疡病基因,利用Illumina HiSeq测序平台对不同时间接种溃疡病菌的毛花猕猴桃进行mRNA高通量测序。结果表明,共获得68731个Unigene,其中有48414个基因注释到Nr、SwissProt、KEGG和COG/KOG数据库。KOG注释显示,通用功能预测、转录后修饰、蛋白代谢、伴侣蛋白、信号转导机制,翻译、核糖体结构和生物合成所占比例最高;GO注释表明参与生物过程的差异表达基因数目最多,细胞组分和分子功能次之;KEGG通路分析表明差异基因的富集以生物合成和代谢为主。接种后不同时间点共获得63050个差异表达基因。通过SSR位点分析,从68731个Unigene中鉴定出13652个SSRs位点;同时获得了各类型转录因子1580个;R基因3727个。本研究结果对猕猴桃抗溃疡病基因发掘和抗溃疡病新品种选育具有重要的理论指导和育种实践意义。     

梁山慈竹高通量转录组测序及差异表达基因分析

分析梁山慈竹及其体细胞突变体的转录组,挖掘功能基因并对其差异表达基因进行筛选和分析,为梁山慈竹遗传改良提供理论依据。利用RNA-Seq技术进行转录组测序,对测序结果进行de novo拼接和功能注释;并对差异表达基因进行筛选及COG、GO、KEGG数据库中进行比对注释,此外,基于Swiss-Prot功能注释结果,分析纤维素和木质素相关功能基因的表达量差异。测序结果表明,共获得86 575 631条reads,de novo组装得到84 741条unigenes,共有49 829条被Nr、COG、GO、KEGG、Swiss-Prot注释。从梁山慈竹实生植株(对照)和体细胞突变体No.30这2个测序样本中,筛选出3 572条差异表达unigenes,757条差异表达unigenes在COG分类体系中具有详细的蛋白功能释义,2 213条差异表达unigenes在GO数据库具有功能定义,385条unigenes被注释到94条KEGG Pathways中。纤维素合成相关纤维素合酶、过氧化物酶、泛素连接酶和热休克蛋白在梁山慈竹体细胞突变体No.30中表达量升高,木质素合成相关MYB4、4-香豆酸CoA连接酶、肉桂醇脱氢酶、肉桂酰-CoA还原酶和漆酶在突变体中表达量降低。提供了全面的梁山慈竹转录组信息,获得了一批在梁山慈竹纤维素和木质素生物合成过程中有重要功能的基因序列。     

水稻苗期低温应答转录组分析

温度对水稻生长发育及产量建成很重要,为了进一步明确低温应答对水稻苗期生长发育的调节机制,以粳稻品种中花-11为试验材料,在17℃低温处理的条件下,利用新一代高通量测序手段-转录组测序(RNA-Seq)开展水稻苗期低温应答转录组分析。通过转录组分析,分别获得干净的高通量序列为46 384 074条和52 483 052条,鉴定并筛选出2 044个差异表达基因,其中,1 252个基因表达上调,792个基因表达下调; GO注释分析将全部的差异基因分为分子功能、生物学过程和细胞组分三大类,各自占比为63. 29%,8. 81%,27. 90%。利用KEGG数据库具体分析低温调节的2个通路光合作用通路和苯丙氨酸代谢通路中差异表达基因的数量及位置,并进一步预测了新基因的功能及与光合作用和苯丙氨酸代谢相关基因的互作蛋白基因。结果表明,低温胁迫对水稻的影响主要体现在生物学过程中,在光合作用、次生代谢的生物合成和苯丙氨酸代谢均有明显作用。另外,差异表达的WRKY转录因子为耐冷机制进一步研究提供方向与依据。     

白芨转录组特性分析

白芨(Bletilla striata)具有较高的药用、经济和观赏价值,但是其基因组和转录组序列未知,严重影响了其的研究开发和利用。本研究采用His4000测序平台对白芨的全株进行了转录组测序分析,共获得原始数据6.8 G,有效数据6.7 G,243 410条Unigene,经过与NR、GO、KOG及KEGG等数据库进行比较分析后,83 541条Unigene被注释到NR数据库,50 178条Unigene被注释到GO数据库,10 007条Unigene在KOG数据库获得注释,43 637条Unigene在Swissprot数据库获得注释,15 321条被注释到KEGG代谢途径中,2 021条Unigene参与了糖类代谢,1 309条Unigene参与了氨基酸合成和代谢,120条Unigene参与了萜类合成,106条转录因子与代谢相关;微卫星位点有31 958个,其中单核苷酸最多,15 709个,占49.16%,其次为二核苷酸和三核苷酸,分别有9 145个和7 104个,占28.62%和22.23%。本研究为白芨的重要功能基因挖掘、遗传育种及其研究开发提供了参考和依据。     

基于Illumina HiSeq 2000测序技术对大花序桉根的转录组特性分析

大花序桉(Eucalyptus cloeziana)是中国重要用材林树种,但目前对其生物信息学研究缓慢,在一定程度上限制了大花序桉分子育种以及品种改良。采用Illumina Hi Seq TM 2000高通量测序平台对1年生大花序桉根系转录组测序和de novo组装,将得到的Unigene与公共数据库比对,同时进行转录组分析。结果显示,组装获得53 433条Unigene,其平均长度890 bp,N50长度为1 587 bp;有34 700条Unigene获得注释信息,占全部Unigene的64.94%,其中19 327条Unigene注释到KOG数据库,被分到25个类别中,共得到32 302个KOG功能注释信息;有11 197条Unigene注释到GO数据库中,共获得了54 971个GO注释功能,归于细胞组分、分子功能和生物学过程三大类;有12 181条Unigene得到KEGG注释,其中6 493条Unigene归入128条代谢途径,发现有57条Unigene参与氮代谢途径。通过软件查找获得13 290个SSR位点,二核苷酸重复类型占的频率最高,其次是三核苷酸、四核苷酸和六核苷酸,五核苷酸重复频率最低。本研究首次对大花序桉转录组进行分析,为深入开展大花序桉分子生物学研究提供基础数据来源。     

基于RNA-Seq技术的盐胁迫向日葵转录组信息分析

本研究采用RNA-Seq技术,进行了120 mmol/L Na Cl胁迫下向日葵耐盐品种P50和盐敏感品种P29转录组测序和De novo组装。P50和P29分别获得106 275和119 350条unigenes,平均长度分别为716 bp和685 bp。2个样品的长序列聚类后,获得了110 751条All-unigenes,总长度为96 970 017 nt,平均长度为876 nt。对获得的All-unigenes进行功能注释,注释到NR、NT、Swiss-Prot、KEGG、COG和GO数据库的All-unigenes分别是71 610、52 569、50 827、46 676、32 215和52 406个,所有注释上的All-unigenes是77 536个。基于NR和KEGG的注释结果,对All-unigenes进行GO和KEGG功能分类,分别获得55个功能小类和128个Pathways注释。研究结果为进一步进行P50和P29转录组差异表达基因分析及耐盐相关基因挖掘奠定了基础。     

基于高通量测序的柚木边材转录组分析

本研究利用Illumina HiseqTM4000测序平台对柚木(Tectona grandis L. F.)边材组织进行转录组测序,获得39.65 Gb的数据。拼接组装共得到90 843个Unigene,平均长度、N50以及GC含量分别为1 415 bp,2 208 bp和41.28%。将获得的Unigene与七大功能数据库进行比对,分别有64 416 (NR:70.91%)、69 281 (NT:76.26%)、28 777 (COG:31.68%)、18 630 (GO:20.51%)、49 594 (KEGG:54.59%)、44 707 (Swissprot:49.21%)以及50 938 (Interpro:56.07%)个Unigene获得功能注释。经过GO数据库的比对分析,18 630个Unigene被注释到生物过程、细胞组分和分子功能3大类别55个亚类。与COG数据库进行比对分析,28 777个注释Unigene按功能被划分为25类。基于KEGG数据库,44 595个Unigene序列注释到6大类,21个亚类代谢通路中。根据注释结果预测出2 772个编码转录因子的Unigene,检测出26 773个SSR位点,以及39 856个SNP位点。本研究为柚木分子育种工作的开展提供数据和参考。     

低温胁迫对冰菜转录组水平响应分析

利用二代高通量测序技术对低温胁迫处理的冰菜进行测序,构建冰菜转录组数据库.分别得到24.13 Gb有效数据和24045条Unigene的注释,得到DEGs 1902个(T0 vs T1)和2134个(T0 vs T2).T0 vs T1组和T0 vs T2组分别有40和41个功能小类化归GO数据库;分别有20和24个功能分类注释到KOG数据库;155和272条基因注释到KEGG数据库,并分别富集在74和105条代谢通路.T0 vs T1组DEGs主要注释到植物信号转导等4个代谢通路;T0 vs T2组DEGs注释到苯丙醇类生物合成等11个代谢通路,其中正向影响代谢途径:丙酮醇类生物合成、嘌呤代谢、谷胱甘肽代谢、脂肪酸代谢、类黄酮代谢、氨基酸的生物合成代谢途径等;负向影响代谢途径:植物-病原互作、植物激素信号转导、淀粉和蔗糖代谢等途径.通过对淀粉和蔗糖代谢途径关键基因分析表明:低温胁迫1 h (T1),海藻糖6-磷酸合酶、海藻糖-6-磷酸酯酶、β-淀粉酶、葡萄糖-1-磷酸腺苷酰转移酶、糖原磷酸化酶等5个关键基因表现为上调表达,未见下调表达基因;低温胁迫36 h (T2),海藻糖-6-磷酸合酶、己糖激酶、β-淀粉酶等3个关键基因上调表达,葡萄糖内酯-1,3-β-葡萄糖苷酶基因下调表达.选取淀粉和蔗糖代谢途径中8个DEGs,经RT-qPCR分析,8个DEGs的相对表达量与转录组表达水平相符.