首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
IBA诱导楸树嫩枝扦插不定根发育的转录组分析   总被引:1,自引:0,他引:1  
【目的】利用RNA-seq技术对楸树易生根品种‘豫楸1号’的嫩枝扦插生根的4个发育阶段进行转录组测序和分析,为阐明楸树不定根发育的分子机制奠定理论基础。【方法】首先利用浓度为2 g·L-1IBA处理‘豫楸1号’嫩枝,然后分别提取处理后0天(对照)、1天(激活期)、15天(愈伤形成期)的插穗基部的RNA及25天(不定根形成期)和35天(不定根伸长期)的根部的RNA进行转录组测序;接着测序得到raw reads,通过去除接头、重复序列、低质量的序列,得到高质量的clean reads,使用Trinity软件进行拼接获得contigs;随后,利用BLASTx、RSEM、Cytoscape及Me V4.9.0分子生物学软件对测序结果进行注释和差异基因鉴定;最后随机选取15个差异表达的RNA-Seq数据利用q PCR验证基因的表达。【结果】在62 955个unigenes中,31 646(50.26%)个unigenenes获得了注释;差异基因分析发现了11 100个差异基因,其中包括10 200个特异的和900个共同的差异基因;GO富集分析发现‘细胞骨架’仅在激活期显著富集,而‘DNA代谢过程’仅在愈伤组织形成期显著富集;KEGG富集分析显示参与丙烷合成、糖酵解和植物激素代谢等途径的基因对楸树不定根的形成具有重要的作用,并且随着楸树不定根发育进程的发展,参与糖酵解途径的基因数目逐渐减少而参与苯丙素合成的基因数量却在持续增加,表明在IBA刺激后,代谢通路的动态变化响应了不定根的发育进程。【结论】通过对楸树‘豫楸1号’不定根发育的4个阶段的转录组分析,发现细胞分裂素和乙烯间的互作促进楸树愈伤形成而生长素和油菜素内酯间的互作则促进了楸树不定根的伸长。虽然本研究不能完全解释‘豫楸1号’不定根形成的分子机制,但它可以作为进一步探索与这一复杂过程相关的候选途径和基因的有力工具,可为解析楸树不定根发育的分子机制和其他近缘物种的研究提供帮助。  相似文献   

2.
[目的]研究避光处理对麻竹笋苦涩味物质合成相关基因表达的影响。[方法]利用Illumina HiSeq~(TM)2500平台对自然生长(CK)和覆土处理(EP)2种类型的麻竹笋进行转录组测序,并对差异表达基因进行分析。[结果]转录组测序共获得36.45 Gb原始数据,经组装去冗余处理得到53 388个Unigene,将所得Unigene比对到Nr、Pfam、COG、Swissprot、GO和KEGG数据库中进行功能注释,发现共有31 462条Unigene与其它物种的基因具有同源性。对CK和EP处理所测得的Unigene进行表达量的比较分析,筛选出1 846个差异基因,其中,在EP处理中上调表达基因998个,下调表达基因848个。由KEGG代谢通路分析发现,差异基因显著地富集在32个代谢途径中,其中,包含与苦涩味物质合成相关的糖酵解途径、苯丙氨酸、酪氨酸和色氨酸合成途径等。进一步研究发现,参与麻竹笋芳香族氨基酸、单宁合成的PFK、ENO、PPY-AT/HPP-AT、LAR等酶基因在EP处理中表达量下降,荧光定量PCR验证结果与测序结果基本一致。[结论]避光处理抑制了苯丙氨酸、酪氨酸、单宁生物合成关键酶基因的表达,可能最终影响麻竹笋苦涩味物质的合成。  相似文献   

3.
黑果枸杞果实发育过程中转录组测序分析   总被引:1,自引:0,他引:1  
【目的】研究黑果枸杞果实在不同发育阶段基因表达谱的变化情况,为进一步开展黑果枸杞果实发育过程中的分子生物学研究提供基础资料。【方法】以黑果枸杞果实发育过程中青果期、变色期和成熟期的果实为材料,利用新一代高通量测序技术平台Illumina Solexa进行转录组测序和数据de novo组装,并对得到的Unigene进行功能注释、分类及代谢通路分析。【结果】基于转录组测序数据得到N50值为1 743 bp及平均长度为1 262.65bp的黑果枸杞Unigene共有43 573条。在GO、KEGG、KOG、NR和SwissProt数据库上得到注释的Unigene总数为23 723条,占总Unigene的54.44%,有3 726条Unigene在这5个数据库中同时得以注释,但还有19 850个Unigene在这些数据库中没有得到注释。23 559条和17 212条Unigene分别在Nr和SwissProt数据库有同源比对信息。GO数据库中注释到的15 064个Unigene分为细胞组分、分子功能及生物学过程等3大类62个功能组。KOG数据库中注释到的13 128个Unigene功能系统分为25类,其中黄酮类代谢途径所属的Q类(次生代谢产物生物合成、运输和代谢)共获得了520个Unigene注释。以KEGG数据库为参考,将4 951个Unigene定位到215个代谢途径分支,其中有45个Unigene与类黄酮生物合成相关。在黑果枸杞果实转录组中发现16 815个SSR位点:最多的为单核苷酸SSR,占72.92%;最少的是6核苷酸,占0.04%。【结论】共得到43 573条Unigene,有23 723条Unigene在GO、KEGG、KOG、NR和SwissProt 5个数据库中得以注释,筛选出与类黄酮生物合成有关的Unigene 45个,这为与黑果枸杞果实品质有关基因筛选、克隆和功能分析等提供研究基础和理论依据。  相似文献   

4.
【目的】以自然突变的楸树雄性不育花芽为研究材料,分析与楸树雄性不育相关基因的表达模式,从基因的表达水平上揭示楸树雄性不育的分子机制,为楸树及木本植物的雄性不育研究提供有价值的参考。【方法】采用高通量测序技术对自然状态下楸树的雄性不育的花芽以及可育的花芽进行转录组测序,通过生物信息学对可育花芽与不育花芽的转录本进行比较分析,预测和筛选出与楸树雄性不育有关的基因。【结果】转录组测序共产生27.18 Gb数据,组装并去冗余后得到86 076个Unigene,然后将Unigene比对到7大功能数据库(NR,NT,GO,COG,KEGG,Swissprot,Interpro)进行功能注释,最终被7大数据库中任意一个数据库注释上的Unigene总数为64 600(75.05%)。将试验组(雄性不育花芽,SL)与对照组(可育花芽,FL)所测得的Unigene进行表达量的分析后,筛选出表达量有差异且可信度高的差异表达基因。在噪音分布中(差异表达倍数在2以上,可信度在0.8以上)共筛选出试验组中有6 915个基因上调表达,3 504个基因下调表达。在泊松分布中(差异表达倍数在2以上,错误发生率在0.001以下),SL-1 vs FL-1、SL-2 vs FL-2、SL-3 vs FL-3三个生物学重复中得出差异上调表达基因分别有13 979,13 513,13 055个,差异下调表达基因分别有12 170,13 807,10 411个。差异基因的GO功能分析表明,在生物过程中显著性富集的条目集群频率较高的是生殖过程、生殖发育过程、生殖系统发育、生殖结构发育,在分子功能中只有生长素流出跨膜转运蛋白活性显著性富集。差异基因KEGG通路分析中,差异基因映射到127个不同的生物途径,显著富集代谢通路中差异基因的生物途径主要为:代谢途径、次生代谢产物的生物合成、剪接体、RNA转运以及甘油磷脂及淀粉和蔗糖的代谢。在用差异基因与同源基因比对分析中,共比对出246个同源性高的Unigene,其COG功能分类多聚集在RNA加工和修饰,细胞周期控制、细胞分裂、染色体分区,转录等,同时将这些同源性高的差异表达基因映射到了丙酮酸代谢、植物激素信号转导途径中。【结论】楸树雄性不育的形成与生殖发育的多个过程、丙酮酸代谢途径、生长素流出跨膜转运蛋白活性以及油菜素内酯介导的信号转导通路密切相关。根据分析的结果和已经完成的相关细胞学观察,推测楸树雄性不育的形成可能是楸树体内丙酮酸代谢过程异常,抑制油菜素内酯的合成,使绒毡层发育异常并进一步影响到小孢子减数分裂,最终导致不育花粉的形成。  相似文献   

5.
油茶花发育转录组测序及相关基因表达分析   总被引:1,自引:0,他引:1  
通过对油茶成花过程的转录组测序及其成花相关基因的表达分析,结果表明:转录组测序总共获得28448847个 reads,5742023480 bp 数据量,GC 含量为46.52%;拼接成大于200 bp 以上的 Unigenes 有94476条,N50长度为806 bp,其中1 kbp以上的Unigenes 共12643条,占Unigene总数的13.38%; Unigenes 在各数据库中功能注释数目,在 COG中有9095条,在 GO中有27201条,在 KEGG中有6431条,在 Swissprot中有24534条,在TrEMBL中有36393条,在Nr中有36400条,在Nt中有30858条;茎尖中FLC,FCA和FT基因表达量较 AP1,AP2和 PI 基因低,FT基因( ID:Unigene60063)是油茶成花的关键基因,PI 基因( ID:Unigene56059)与雄蕊发育关系紧密。  相似文献   

6.
[目的]采用高通量测序技术Illumina Mi Seq250获得蕨类植物芒萁的孢子体转录组数据,以期为芒萁的生长、发育、代谢调控、微进化机制分析等提供重要的分子信息。[方法]应用生物信息学方法对测序获得的大量单基因簇(Unigene)进行基因功能注释、代谢途径及微卫星分析等。[结果]本研究共获得18 463 296条序列读取片段(reads),总碱基数为4.62Gbp序列信息,经序列组装最终得到63 169个Unigene,平均单条Unigene长度为863 bp,N50为1 587 bp,其中分布在200 500 bp长度区间Unigene占总数的55.4%。数据库中的序列同源性比较表明,26 826个Unigene与其他物种的已知基因具有不同程度的同源性。芒萁转录组中的Unigene根据GO功能大致可分为细胞组成、分子功能和生物学过程3大类47个分支,其中有大量的Unigene与细胞进程、绑定活性、代谢过程和催化活性相关。将Unigene与COG数据库进行比对,根据其功能大致可分为26类。以KEGG数据库作为参考,依据代谢途径可将Unigene定位到257个代谢途径分支。此外,利用MISA软件检索2 6碱基微卫星,共找到13 286个SSR。在不同长度微卫星中,三核苷重复数量最多,占总数的40.41%。在各重复基序类型中出现频率最高的为AG/CT(14.45%)与AAG/CTT(12.39%)。利用重复基序开发的多态性SSR标记,可应用于芒萁不同个体的基因型分型鉴定。[结论]本研究获得了较高质量的芒萁转录组数据库,揭示了芒萁孢子体生长发育过程中表达基因的功能总体特征,可为芒萁进一步的功能基因挖掘和分子标记规模化开发奠定基础。  相似文献   

7.
【目的】对白腐菌偏肿革裥菌在木质和非木质环境下的转录组进行测序,为白腐菌降解木材的机制研究提供支持。【方法】采用高通量测序技术对木屑和非木屑处理条件下的菌丝样本进行转录组测序。利用eggNOG、GO、KEGG等数据库对转录本进行注释和比较分析,预测和筛选出偏肿革裥菌与木材降解有关的基因。应用荧光定量PCR(qPCR)检测mnp2等11个与木质素降解相关的差异基因在木屑处理5天条件下的表达量,结合转录组数据对这些基因的表达量变化趋势进行分析。【结果】偏肿革裥菌转录组测序共得到38.9 Gb Clean Data,各样品Clean Data平均达到6.49 Gb,各样品的Reads与参考基因组的比对效率在71.23%~74.25%之间。使用edgeR软件进行差异表达分析,得到差异表达基因898个,其中上调351个、下调547个。差异基因被注释到GO数据库的有251个,注释到KEGG数据库的有223个,注释到eggNOG数据库的有704个。eggNOG分析表明,差异基因表达多聚集在能量产生和转换、转录后修饰、蛋白质代谢、伴侣关系、次生代谢物的生物合成、碳水化合物的运输和分解代谢等功能分类下。GO分析表明,显著性富集与频率较高的生物过程是丙酮酸代谢过程、类异戊二烯生物合成过程、天冬氨酸家族氨基酸代谢过程和木质素分解过程。在KEGG通路分析中,差异基因分布到86个不同的生物途径,差异基因显著富集的代谢通路主要为:碳代谢、柠檬酸循环、芳香化合物降解、乙醛酸代谢。qPCR结果表明在木屑处理条件下基因mnp2、mnp3、lip9、lip2、laccase1和nadp的表达量明显上调;mnp10 s、mrp、cyp450-1、GroES1和GroES2的表达量明显下调,qPCR与转录组测序结果在表达量差异倍数上存在较小偏差,但整体趋势一致,可以证明转录组测序结果是正确可靠的。【结论】偏肿革裥菌降解木材与碳代谢、芳香化合物降解等通路密切相关;与木质素分解等生物过程密切相关。根据基因功能注释的结果得到11个与白腐菌降解木质素相关的重要差异表达基因。  相似文献   

8.
为研究杉木Cunninghamialanceolata杂交组合内子代生长性状分离的分子机制,以杉木杂种和亲本为研究对象,从基因表达水平揭示超亲杂种优于低亲杂种的形成分子机制,采用IlluminaHiseq4000高通量测序技术对不同生长势的杉木杂种(龙15×1339)的超高亲子代(HF1)和超低亲子代(LF2)及其亲本(龙15和1339)进行无参转录组测序和差异比较。结果表明,所有样本经转录组测序共产生Clean reads 5.8E+08条,总拼接长度为49 803 726nt,将Clean reads在6个数据库(Nr,Swiss-prot,KOG,KEGG,Pfam,GO)进行BLASTX分析,比对结果产生80 171个基因;杉木同一组合内的HF1比LF2的生产力高,是因为HF1在不同的GO和KEGG terms内和terms间两个层面上的差异表达基因的分布处于不均匀、不平衡状态;HF1的GO terms和KEEG terms的差异基因的分布不均匀、不平衡;LF2 GO terms和KEEG terms的差异基因的分布趋于均匀、平衡状态。利用分子生物学技术对育种群体展开遗传多样性研究,选择表达基因互补的亲本进行杂交组配,这样可以降低杂种间生长性状的分离,以达到收获期大体一致的目的,是一条可行的技术路线。  相似文献   

9.
【目的】白花泡桐Paulownia fortunei (Seem.) Hemsl.为中国南方重要的尾矿造林树种。为了揭示泡桐根系受重金属胁迫后的转录调控机制,明确泡桐根系响应重金属胁迫的相关基因和重要通路。【方法】以生长在铅锌矿、南方红壤的白花泡桐根为实验材料,进行高通量测序,将组装的基因序列在Nr、Nt、Pfam、KOG/COG、SwissProt、KEGG和GO数据库中进行基因功能注释,同时对在红壤、铅锌矿中差异表达的泡桐根部基因进行了分析。【结果】1)经过测序共得到366 701 960条reads序列,总碱基数为53.90 Gb,GC平均含量达44.58%。将测序所得reads再经Trinity软件组装后得到126 061条基因序列,其中96 133条序列获得了注释信息,占总数的76.25%。2)注释到Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO数据库的基因序列数分别为63 938、64 685、68 840、29 886、66 433、28 747、68 840条。3)相对于红壤生长条件,铅锌矿胁迫时泡桐根部分别有154个基因上调和4 143个基因下调表达。4)差异表达基因的GO富集性分析发现,这些基因的功能主要集中于结合(Binding)、细胞器(Organelle)、细胞蛋白代谢(Cellular protein metabolic process)。将差异表达基因在KEGG数据库中比对发现,差异基因主要富集于核糖体、内质网蛋白加工、剪接体、胞吞作用、RNA转运、吞噬代谢、mRNA监测等7个代谢通路。其中,核糖体是最主要的代谢途径。q RT-PCR所检测的6个差异基因表达模式与RNA-seq结果保持一致,证实了RNA-seq结果的可靠性。【结论】本研究提供了泡桐根部在铅锌矿和红壤条件下的转录组信息,并分析了泡桐根部转录组变化情况,为寻找泡桐根部响应重金属胁迫的基因打好基础。  相似文献   

10.
【目的】研究美洲黑杨杂种F1的基因表达模式,从基因表达水平揭示杂种优势形成机制,为杨树杂种优势的深度开发利用提供有益参考。【方法】采用Illumina HiSeq 2000高通量测序技术对不同生长势美洲黑杨杂种F1(3个生长势超过双亲、2个生长势低于亲本)及其亲本进行转录组测序和差异比较。【结果】转录组测序共产生171 154 127条Reads(平均长度200 bp,约31.32 Gb),将处理后Reads与毛果杨参考基因组数据库比对,有61.89%的Reads可以与参考基因组匹配。超亲杂种F1 Vs亲本中筛选出342个基因显著差异表达(上调表达87个,下调表达255个);低亲杂种F1 Vs亲本共筛选出577个基因差异表达(上调表达146个,下调表达431个);超亲杂种F1 Vs低亲杂种F1中共筛选出486个基因显著差异表达(上调表达200个,下调表达286个)。筛选出377个对杂种优势产生具有重要作用的差异表达基因,其中,有4个差异基因在超亲F1 Vs亲本、低亲F1 Vs亲本及超亲F1 Vs低亲F1中为下调表达;72个下调表达差异基因在超亲F1 Vs亲本及低亲F1 Vs亲本2个对比组中均显著差异表达;有129个差异基因仅在超亲F1 Vs亲本及超亲F1 Vs低亲F1 2个对比组中显著表达,包括19个差异基因显著上调表达,110个差异基因显著下调表达;有172个基因仅在低亲F1 Vs亲本及超亲F1 Vs低亲F1中显著差异表达。功能注释及代谢途径分析结果表明,超亲F1 Vs亲本、超亲F1 Vs低亲及低亲F1 Vs亲本3个对比组中分别有167个、233个及280个差异显著基因能被功能注释,分别涉及46个、45个及51个生物学功能,在分子功能、生物学过程、细胞组分3种GO分类中,其中差异基因主要富集在"催化活性"、"胺的代谢过程"及"氧化还原酶活性"等功能中,参与碳水化合物代谢、氨基酸代谢、能量代谢途径及外来物质的降解和代谢途径。【结论】杂种优势的形成,可能是由于相关基因的显著差异表达,调控光合作用、物质代谢吸收等与生长紧密联系的代谢活动,进而促进生长优势的产生。  相似文献   

11.
【目的】油桐种仁中产出的桐油经济利用价值极高。桐油中α-桐酸的含量高达70%,然而植物体内α-桐酸代谢通路的研究还未见报道,这对直接筛选油桐α-桐酸代谢通路相关酶基因造成一定的困难。α-亚麻酸作为α-桐酸的同分异构体,其代谢通路的研究则较为深入,能为α-桐酸代谢通路的解析提供参考。因此,本研究期望在油桐种仁转录组数据的基础上,解析油桐的α-亚麻酸代谢途径,为油桐α-桐酸代谢机理的阐明提供理论参考。此外,通过调控这些基因的表达模式以及开发与之紧密连锁的分子标记,可大大加快油桐遗传改良和分子育种的进程。【方法】采用 RNA-Seq技术对油桐种仁3个不同油脂合成期的转录组进行比较,获得大量差异表达的Unigene,并将这些Unigene归类于128个代谢途径。在此基础上,通过GO分类和Pathway富集性分析,解析油桐α-亚麻酸代谢通路并分析通路中相关酶基因在油脂合成期的表达变化规律。【结果】通过对Ⅰ,Ⅱ,Ⅲ期的油桐种仁RNA测序,共获得长度为200~3000个核苷酸的非冗余Unigene序列58439条,其中能够比对到公共数据库中已知基因序列的 Unigene共有41059条,占所有非冗余基因的70.3%。不同长度的非冗余 Unigene序列与数据库中序列匹配的效率不同,越长的序列匹配效率越高。序列长度大于2000 bp的序列匹配效率达到98.28%,而500~1000 bp和100~500 bp的序列分别只有78.86%和48.99%的匹配效率。3个种仁油脂合成期的转录组数据中共有105个 Unigene可被富集于α-亚麻酸代谢途径,占所有非冗余 Unigene的0.47%。从3个转录组数据的两两比较中鉴别出一些差异表达Unigene,其中也有一些可被富集于α-亚麻酸代谢途径。通过在KEGG数据库中进行检索后发现,105个 Unigene序列分别对应于14个α-亚麻酸代谢途径关键酶基因,这些基因在其他物种中都有同源基因与之对应。通过基因表达模式分析发现,整体上与合成代谢相关的基因在油脂合成期呈现上调的表达模式,而与分解代谢相关的基因则呈现下调的表达模式。【结论】在油桐种仁转录组数据的基础上,解析油桐α-亚麻酸代谢途径,获得与α-亚麻酸代谢相关的重要酶基因并分析它们在油脂合成期的表达模式,这对后续研究具有重要的启示作用。  相似文献   

12.
油松雌雄球花高通量基因表达谱芯片分析   总被引:1,自引:0,他引:1  
生殖调控是针叶树种遗传改良与良种生产过程中的重要环节,已经有许多实践经验积累,然而对其中雌雄球花发育及调控的分子生物学基础研究仍相当匮乏.以油松多个体多组织样本均一化cDNA文库为材料,通过转录组测序,获得46 584个油松单基因簇;基于这些基因信息,设计并定制油松表达谱芯片;通过对雌雄球花基因表达芯片分析,筛选到2 460个在球花发育后期显著差异表达的基因,其中1 309个在雄球花优势表达,1 151个在雌球花优势表达.分别从雄球花与雌球花优势表达的基因中各挑选1个表达水平差异超过800倍的基因,由于其功能未知,分别命名为MCE1(Male Cone Expressed 1)与FCE1(Female Cone Expressed 1),并进一步分析它们在雌雄球花不同发育阶段的表达规律.结果表明:MCE1只在雄球花各发育阶段高水平表达,而在雌球花各发育阶段表达水平极低,推测其参与雄球花的整个发育过程;FCE1在早期雄球花中大量表达,但随着雄球花的发育其表达量会急剧降低,而在雌球花中却始终稳定地高水平表达,推测其在雌雄球花不同发育阶段发挥不同的作用.研究结果可为油松及近缘树种雌雄球花发育的分子生物学研究提供基础和新的思路.  相似文献   

13.
[目的]利用高通量转录组测序技术,在硝态氮或铵态氮处理条件下,对杨树根尖差异表达基因进行了筛选和研究,同时分析和描述了差异表达基因对杨树根尖生长发育的影响,为后续开发高氮素吸收利用效率的杨树新种质提供科学依据。[方法]以灰杨幼苗根尖为材料,用0.5mmol·L-1硝态氮(NO3-)和0.5mmol·L-1铵态氮(NH4+)对幼苗处理10 d,并对植株根尖进行转录组测序以及生物信息学分析。[结果]硝态氮处理下的主根长度几乎是铵态氮处理条件下的一倍。从两种不同氮形态处理杨树根尖转录组文库中,筛选到2 207个差异表达基因。通过差异基因GO和KEGG功能聚类分析,分别获得50个GO功能聚类和20条KEGG通路。进一步利用MapMan分析,筛选出36个氮代谢通路过程、各类氨基酸的生物合成以及代谢过程相关的差异表达基因。对这些差异基因互作调控网络分析发现,硝酸还原酶(Potri.005G172400)基因通过响应不同氮形态,在影响杨树根尖生长发育过程中发挥了重要作用...  相似文献   

14.
【目的】深入研究云南金花茶优良性状、基因功能和基因分布特征,可为云南金花茶乃至山茶属植物功能基因的探索以及该物种的遗传保护提供基础生物信息学理论参考。【方法】采用Illumina Hi Seq 2000技术,对云南金花茶叶片进行转录组测序及相关数据分析,在对测序数据整理、de novo组装后,获得Unigene序列,继而利用相关生物信息学数据库进行序列比对。【结果】组装后,共获得95 979条Unigenes,其中N50(拼接转录本不小于总长50%的长度)为1 660 nt,平均长度为1 124 nt,Q20和Q30序列(处理后质量高于20和30的碱基)分别占96.39%和91.28%。经比对,在Nr、Nt、Swiss-Prot中能得到注释的Unigene分别为58 830、43 623、44 315条。在GO中,有41 905条(43.66%)Unigenes注释224 129个GO功能,将其分为3个大类和56亚类,所占比例最多的为生物过程这一类功能。在KOG中,有23 499条(24.48%)Unigenes注释26 430个KOG功能信息,将其分为26个基因功能大类,其中表达量较高的分别是一般功能基因和翻译后修饰、蛋白质转化、伴侣功能的相关基因。另外,共有23 214条(24.18%)Unigenes在KEGG中得到注释,根据其涉及的相关通路可将其归为19个亚类,其中以代谢通路较为丰富。此外,利用相关数据库和ESTScan软件对云南金花茶转录组Unigene进行CDS比对和预测,共预测到26 428条CDS,其长度集中在100~500 nt,占总CDS的82.10%。【结论】云南金花茶所含基因信息丰富,利用分析得到的所有注释信息可以更深层次地探索其基因组信息和基因分布情况。  相似文献   

15.
为研究牡丹花型的分子机理,挖掘调控牡丹雌蕊瓣化的候选基因,揭示牡丹不同花型的分子机理,以牡丹‘黑海撒金’(tree peony ‘HEIHAISAJIN’)和‘璎珞宝珠’(‘YINGLUOBAOZHU’)为材料,采用Illumina Hiseq测序技术进行花瓣的转录组测序,通过序列的拼接、组装和过滤,得到牡丹的转录组;比较‘黑海撒金’和‘璎珞宝珠’中转录本特定基因的表达量,获得差异表达基因;将差异表达基因进行GO功能注释和KEGG代谢通路分类,获得与花发育相关的功能基因。通过de novo的拼接和组装,获得了49 191个unigenes;差异显著性分析得到2 735个差异表达基因(DEGs),其中1 082个为上调基因,1 653个差异表达基因下调。GO功能注释发现,DEGs主要分布在25个功能区,最显著富集的为90 S preribosome。KEGG代谢通路中显著富集的通路有13条,最显著富集的通路为Ribosome。对功能基因进行挖掘和比对,发现与开花相关的ABCDEF模型基因有3个,分别为agamous-like MADS-box protein 65(AGL65), floral homeotic protein APETALA 2(AP2)和EMBRYO SAC DEVELOPMENT ARREST 3(EDA3)。认为利用转录组测序技术,可以分析雌雄蕊发育正常和雌蕊瓣化的不同花型的牡丹花瓣的功能基因。3个ABCDEF模型基因中,AGL65和EDA3为首次在牡丹中发现的花器官发育基因。说明这3个基因很可能是调控牡丹雌蕊瓣化的关键基因,参与牡丹花型的发育。该研究系统地挖掘与花型相关的基因,并发现了新的可能调控牡丹花型基因的成员,不仅为调控牡丹花器官发育挖掘了新的重要功能基因,也为揭示牡丹花型发育的分子机理奠定了重要的基础。  相似文献   

16.
17.
[目的]了解马尾松干旱胁迫转录组序列的功能分布和SSR位点分布特征,并探索与干旱关联的SSR位点。[方法]对马尾松幼苗进行持续干旱胁迫,选取干旱10、15、25 d及正常供水对照的马尾松针叶样品,通过总RNA提取、Illumina测序、raw reads去冗、Trinity拼接获得unigenes。利用Blast比对,对Unigene序列进行GO、KOG、KEGG功能注释及分类;利用Misa软件进行SSR位点批量搜寻,primer 3. 0软件进行规模化SSR引物设计;利用GOSeq(1. 10. 0)、KOBAS(v2. 0. 12)软件对差异表达的含SSR位点Unigene进行GO、KEGG显著性富集分析。[结果]马尾松转录组194 821个Unigene中,101 806个Unigene获得注释,包含64 943个GO功能注释、35 880个KOG功能注释及30 882个KEGG注释。搜寻到6 728个SSR位点,分布于6 367个Unigene中,SSR出现频率为3. 45%;重复类型以单、三、二核苷酸为主,分别占总SSR的35. 82%、33. 03%和25. 22%;重复基序以A/T、AT/AT、AG/CT、AGC/CTG、AAG/CTT为主;基序长度以10 20 bp的短序列SSR为主;基序重复次数以5~10次重复占优势;批量设计13 338对SSR引物。422个含SSR位点Unigene具有差异表达; KEGG富集分析发现,有11个含SSR位点的差异Unigene参与了光合作用、植物激素信号传导及类胡萝卜素合成等3个与干旱响应相关的代谢途径。[结论]从较高质量的马尾松转录组中获得101 806个具有注释的Unigene;从6 367个Unigene中挖掘出6 728个SSR位点; 422个含SSR位点的差异表达Unigene中筛选出11个与干旱关联的SSR功能位点,为抗旱功能基因定位及马尾松抗旱分子机制研究奠定基础。  相似文献   

18.
[目的]研究马尾松抗松材线虫病相关基因,为解析马尾松抗性机理和高抗马尾松早期选择提供理论依据。[方法]对高抗和易感2种马尾松基因型接种松材线虫,在接种后第1、15和30天取样进行转录组高通量测序,通过对高抗和易感马尾松差异基因识别以及富集分析,以筛选与抗松材线虫病相关的组成型差异基因。[结果]对高抗和易感马尾松接种松材线虫后的转录组进行比较,在接种后第1、15和30天分别获得2 866、679和1 657个差异基因,且在接种松材线虫后不同时间点间,共有差异基因相对较少。对差异基因进行GO富集分析,结果显示:在接种后第1、15和30天,最显著的生物学进程是氧化-还原过程,而GO项刺激反应、转录调节以及香叶基二磷酸代谢过程也显著富集。对这几类GO项相关基因进一步分析,发现接种后3个时间点GO项刺激反应中共包括26个R基因,除了2个R基因外,其它R基因表达皆为组成型,在接种松材线虫和对照间差异不显著,而GO项转录调节中仅2个差异基因被nr注释为ERF转录因子,其余为未知基因。接种松材线虫后,ERF转录因子在高抗马尾松上的表达量始终高于易感马尾松上的表达量,且在接种松材线虫和对照间表达量变化不显著,也为组成型基因表达。GO项香叶基二磷酸代谢过程中,GGPPS 3个基因也在高抗马尾松上具有更高表达量,为组成型基因表达。这些基因中TIR-NBS-LRR基因(c65785.graph_c0)、ERF转录因子(c78073.graph_c0)和3个GGPPS基因在高抗马尾松中表达量较高,而在易感马尾松中表达量较低,甚至为0,可作为开发分子标记的候选基因用于鉴定抗性马尾松。[结论]高抗和易感马尾松在接种松材线虫后第1天其表达量差异达到显著水平的基因最多,其中,LRR基因、ERF转录因子和GGPPS与马尾松的抗性相关,部分TIR-NBS-LRR基因、ERF转录因子和GGPPS有望被开发为鉴定高抗马尾松的分子标记。  相似文献   

19.
[目的]为了挖掘比拉底白刺耐盐相关基因,对其盐胁迫下差异表达基因进行筛选分析。[方法]以比拉底白刺幼苗为材料,用200 mmol·L^-1 NaCl对幼苗处理7 d,并对胁迫处理和对照植株叶片进行转录组测序及生物信息学分析。[结果]有效序列组装共得到应答盐胁迫的168463条unigenes和196个差异表达基因。通过差异基因GO和KEGG功能聚类,分别获得64个GO功能小类和25条KEGG通路。进一步基因相互作用网络分析发现,转录调控、氧化还原以及抗逆相关基因在比拉底白刺应答盐胁迫中发挥重要作用,其中,筛选到3个重要的节点基因,分别是热激同源蛋白基因、L型凝集素类受体激酶基因和Win类蛋白基因。[结论]本研究获得了盐胁迫下比拉底白刺的差异表达基因及功能注释信息,有助于理解其耐盐的分子机制,为后续开发耐盐分子标记及通过基因编辑改良植物耐盐特性提供了科学依据。  相似文献   

20.
以杜仲果实和叶片转录组为基础,研究了MEP途径系列基因的表达差异,为杜仲基因的克隆、功能鉴定及重要成分的分子积累机理研究提供重要的依据。结果表明,杜仲转录组中共42条Unigene被MEP途径注释。DXS的DXS2和DXS3在幼果和叶片中表达量具有显著差异,且DXS2在幼果和叶片中的表达量最高;该酶在幼果和成熟果实中有7条Unigene被注释,其中DXS6在成熟果实中特异表达,的表达量具有显著差异。DXR在幼果和叶片中有4条Unigene被注释,其中DXR1、DXR2、DXR3的表达量具有显著差异;该酶在幼果和成熟果实中有3条Unigene被注释,其中DXR3在成熟果实中特异表达,3条Unigene的表达量均无显著差异,DXR1在幼果和成熟果实中表达量均最高。MCT在幼果和叶片中仅1条被注释,表达量差异显著;该酶在幼果和成熟果实中仅1条被注释,其表达量无显著差异。CMK在幼果和叶片中仅CMK1被注释,且表达量无显著差异;该酶在幼和成熟果实中仅CMK1被注释,且表达量无显著差异。MDS在幼果和叶片中仅MDS1被注释,其表达量差异显著;该酶在幼果和成熟果实中有2条Unigene被注释,且表达量无显著差异。HDS在幼果和叶片中仅HDS1被注释,表达量具有显著差异;该酶在幼果和成熟果实中仅HDS1被注释,表达量有显著差异。HDR在幼果和叶片中有5条Unigene被注释,其中HDR1和HDR5的表达量具有显著差异;该酶在幼果和成熟果实中有7条Unigene被注释,其中HDR3、HDR4、HDR5、HDR6的表达量有显著差异。IPI在幼果和叶片中仅IPI1被注释,其表达量无显著差异;该酶在幼果和成熟果实中有2条Unigene被注释,其中IPI1的表达量有显著差异。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号