马尾松转录组测序和分析

本研究首次构建了马尾松均一化cDNA文库,采用Illumina高通量测序技术对转录组进行了测序,利用生物信息学方法开展基因表达谱的研究、功能基因的预测。EST序列拼接获得83680个contig,其中33772个contig被注释为相应的331669对生物学功能,10647个contig被注释具有酶功能。根据KEGGpathway数据库,对马尾松转录组的contig进行Pathway生物学通路的注释和预测,共识别出10647个contig具有对应的1029种酶功能,并关联到135条生物学通路。SSR查找发现,从83680个contig中找到889个SSR位点,占contig总数的比例为1.06%。其中,三核苷酸重复所占比例最高,达到48.37%,其次是六核苷酸重复,为19.12%,比例最低的是四核苷酸重复,仅为4.72%,二核苷酸重复和五核苷酸重复基本相同,分别为14.62%和13.16%。SSR不同重复基元类型中,出现频率最高的为AT/AT,其次是AGC/CTG和AAG/CTT。     

基于高通量测序的野百合转录组分析

本研究基于新一代高通量测序技术平台Illumina HiSeq?4000对野百合进行转录组测序,对物种的转录组序列进行统计,并将得到的数据进行de novo组装,结果共获得47 605条Unigenes,总长度为33 972 306 bp,平均长度为713 bp,N50为 1 204 bp.将获得的 Unigenes ...     

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

为了获得甘薯耐盐转录组序列信息,挖掘差异表达基因及其相关代谢途径,本文以盐胁迫处理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 HiSeq 2000高通量测序平台对桐花树叶片转录组进行测序,经de novo组装后获得73 721个Unigenes,进一步利用7个公共数据库(NR, NT,GO, COG, KEGG, SwissProt和InterPro)进行比对,注释了50 338个Unigenes。结果表明,有34 991个Unigenes参与到136条KEGG通路上,其中在花青素和类黄酮合成途径中的Unigenes分别有49个和176个,这些基因可能参与了桐花树应对逆境的调控过程;此外,从转录组序列中搜索到34 509个SSR位点,其中二碱基重复出现次数最多,出现频率为65.75%;从转录组的数据中预测到了2 255个转录因子,分属于58个家族。本研究结果为桐花树的基因功能分析、分子标记开发和遗传多样性研究提供一定的理论参考。     

3个不同品种百香果转录组分析

以黄果、满天星和紫果3个不同品种百香果果皮和果肉为材料,采用2代Illumina高通量测序技术进行转录组测序,获得黄果果皮47 875 830条、果肉44 325 678条,满天星果皮47 019 772条、果肉45 472 434条,紫果果皮48 789 846条、果肉49 372 968条Clean reads片段,包含黄果果皮7 181 374 500 nt、果肉6 648 851 700 nt,满天星果皮7 052 965 800 nt、果肉6 820 865 100 nt,紫果果皮7 318 476 900 nt、果肉7 405 945 200 nt核苷酸序列信息;经拼接组装,分别获得平均长度为1 487 nt的黄果34 968条Unigene片段,平均长度1 361 nt的满天星41 896条Unigene片段和平均长度1 329 nt的紫果41 701条Unigene片段。GO富集分析表明,在生物过程、细胞组分和分子功能3大类别上,黄果、满天星和紫果间均分布有差异基因。KEGG富集分析表明,黄果与满天星在细胞胞吞作用过程中存在较大差异,黄果与紫果在玉米素生物合成和糖酵解途径中差异基因富集显著,满天星与紫果在缬氨酸、亮氨酸和异亮氨酸生物合成中存在较大差异。     

茶树花转录组微卫星分布特征

利用Perl语言,对茶树花转录组序列进行大通量SSR位点的发掘,发现含SSR的序列10 290条,共12 582个SSR,平均2.41 kb出现一个SSR。在茶树花的转录组中共发现340种碱基重复模式,所占比例最高的是(AG/CT)n(44.99%)。在49 586条注释成功的茶树花Unigene中,共发现10 490个SSR位点,其中位于编码区的1917个,其出现频率仅为0.102 SSR/1000 bp,而非编码区为3.072 SSR/1000 bp。在基因编码区中出现频率最高的是三碱基微卫星(1140,59.5%),其次是六碱基微卫星(524,27.3%)。茶树花转录组所含微卫星以重复长度小于20 bp的序列最多,大于20 bp的仅为25.2%。茶树花转录组中,含微卫星基因的平均表达水平显著低于不含微卫星基因,其中含复杂微卫星基因的平均基因表达水平最低。     

转录组学技术在茶树抗逆性的研究进展

茶树(Camellia sinensis)属山茶科山茶属,是一种重要的木本经济作物,其叶片富含多种活性物质,其中尤以儿茶素、咖啡因、茶氨酸及其他游离氨基酸等活性成分对茶叶品质有重要影响。逆境胁迫是影响茶树生长发育、产量及品质的重要因素,人们对茶树应答逆境胁迫的分子机理越来越关注。转录组学作为功能基因组学的重要组成部分,它从RNA水平上全面研究物种的基因表达情况。利用该技术研究逆境胁迫下茶树基因表达的情况,有利于阐明茶树对逆境胁迫的应答机理。本综述主要介绍应用转录组学技术研究茶树应答生物胁迫和非生物胁迫(干旱,温度,盐碱,重金属等)的进展,以期为茶树抗逆性分子机理的研究提供理论参考。     

扁秆荆三棱不同组织混合样品的转录组分析

扁秆荆三棱是一种常见的稻田杂草,旨在补充其转录组信息,为相关防治工作提供支持。基于高通量测序技术,在Illumina Solexa HiSeq 2000平台上对扁秆荆三棱的茎、叶、根茎和球茎的混合样品进行转录组分析。经拼接组装共获得了59 788个Unigene,序列的平均长度842 bp,N50为1 402 bp。将获得的Unigene与5个通用公共数据库(NR、 Swiss-Prot、 KEGG、 GO、 KOG)进行比对(Evalue1e-5),35 221条Unigene获得了基因注释,占总Unigene的58.91%。在12 823个Unigene中共搜索到9 698个SSR位点,其中二核苷酸和三核苷酸的重复类型占所有SSR位点的95.93%。通过KEGG pathways分析,共有13 141个Unigene参与了291个代谢通路,获得了扁秆荆三棱淀粉合成功能相关Unigene 60个,根茎生长功能相关Unigene 14个。其中赤霉素相关基因、核糖体代谢通路等反映了地下根茎网络系统的扩张趋势以及活跃的能量需求,淀粉合成基因则说明其球茎中的营养储备开始于地下系统扩张的早期,这些发现为深入研究地下根茎网络调控机制,或杂草防治或湿地保护等实践工作提供了参考。     

长瓣兜兰花2个不同时期转录组分析

为了更好地认识长瓣兜兰,并开发其园艺价值,以长瓣兜兰花器官为材料,利用RNA-seq技术对长瓣兜兰花蕾和花朵进行转录组测序.结果表明,共获得95659条unigene.将unigene比对到NR、KOG、Swissprot、KEGG等数据库进行注释,共发现有61629条unigene得到注释,占全部unigene的64...     

桂花(Osmanthus fragrans)转录组SSR特征分析

为解析桂花转录组数据中SSR的特点及潜在利用价值,本研究基于‘早银桂’和‘橙红丹桂’转录组数据信息,利用MISA软件从桂花的117 595个Unigenes中进行SSR位点的查找、分析,共检测到16 015个SSR位点,平均密度是1/5.2 kb。桂花SSR基序97.73%的重复次数在5～10次之间,重复次数最多的是6次(28.12%);SSR长度分布在12～36 bp之间,其中12～16 bp的最多,占66.8%,其总体变化趋势是随着重复次数增加重复次数急速下降。基于已知的SSR位点,利用primer3设计引物并从中挑选出196对引物,其中92对引物(占46.94%)已被注释到参与信号的转录翻译、苯丙氨酸代谢、植物激素信号转导等功能。综上所述‘,早银桂’和‘橙红丹桂’转录组SSR位点出现频率高,类型较丰富、多态性潜能较高,为今后桂花分子育种提供了理论依据。     

不同浓度乙酸处理北美冬青(Ilex verticillata)叶片的转录组分析

为探究乙酸提高北美冬青抗旱能力的分子机制,通过对CK (清水)、T1 (20 mmol/L乙酸溶液)、T2(40 mmol/L乙酸溶液)处理后的北美冬青叶片进行转录组测序及分析。通过组装共得到120 421条Unigene序列,共计73.75%的序列在相关数据库中得到注释;通过3个组差异表达基因比较得到以下结果:与CK相比,T1处理后的北美冬青叶片有10 229个基因表达量上调,559个基因表达量下调;T2分别有10 959个上调和938个下调;T1相比较T2有3 524个上调,4 287个下调。GO富集分析显示差异基因主要涉及结合、催化活性、转运活性、细胞成分的组织与生物发生、细胞定位、生物调节细胞、细胞分离及细胞器等。KEGG分析表明在使用乙酸溶液处理后,北美冬青植株内的乙烯、水杨酸、油菜素内酯等植物激素信号增强。综合以上分析结果,北美冬青可能通过代谢调节、物质能量准备及有害自由基的清除等途径提高处理植株应对干旱胁迫的能力,其中激素通路信号提高倍数最强的乙烯途径值得深入研究。     

特色茶树品种在生态观光茶园景观设计中的应用

以我国丰富的特色茶树品种资源为核心,提出如何科学合理地将这些特色品种资源应用在生态观光茶园总体布局设计、立体设计、季相设计、色彩设计和文化设计中,为我国生态观光茶园建设提供参考。     

基于高通量测序的芍药属杂交种胚转录组分析

为了研究牡丹、芍药远缘杂交胚败育的分子机理,本研究使用BGISEQ-500平台对牡丹、芍药杂交种胚转录组进行测序,通过测序,共得到了92.02 Gb数据。通过拼接组装去冗余后得到了86 195个Unigene,平均长度为1 189 bp。86 195个Unigene在七大功能数据库中进行注释,最终分别有49 172 (NR:57.05%)、38 352(NT:44.49%)、36 477 (Swiss Prot:42.32%)、38 905 (KOG:45.14%)、37 993 (KEGG:44.08%)、26 832 (GO:31.13%)以及37 758 (Pfam:43.81%)个Unigene获得功能注释。同时还检测出21 998个SSR分布于17 567个Unigene中,其中二核苷酸和三核苷酸这两种重复类型分别有5 628个和2 906个。在17 567个Unigene中预测出1 671个编码转录因子。不同发育时期的胚转录组数据中共筛选出13 974个差异表达基因,其中上调基因有8 647个,下调基因有5 327个。该转录组测序分析为寻找牡丹、芍药远缘杂交胚败育相关基因提供了一定的理论参考。     

休眠期后板栗混合花芽的转录组分析

板栗是中国重要的木本粮食作物,具有极高的经济价值,但板栗雄花比例远远高于雌花,是限制板栗产量的主要因素。为了从分子水平上探讨板栗雌雄花的发生机制,本研究利用高通量测序技术对休眠期后板栗的两种混合花芽的转录组进行比较分析。总计获得了3 772个差异表达基因,其中2 173个基因上调,1 599个基因下调。这些差异基因经过GO分类和KEGG pathway富集分析分别归类于47个GO term和251条代谢途径。从转录组数据中获得了春化、光周期、赤霉素(GA)等途径的开花相关基因,还鉴定多个开花相关内源激素途径的基因。将基因VIN3,GAI,ELF3,GID1的qPCR检测结果与其转录组测序结果进行比对,发现虽然两种检测结果在差异倍数上存在差别,但4个基因的表达量的总体变化趋势一致,说明转录组测序结果有较高的可信性。本研究结果为揭示板栗雌雄花性别比例调控机制和花期人工调控提供了一定的科学参考。     

基于RNA-seq的铁线莲转录组信息分析

本研究首次采用高通量测序技术对同一株上三种不同花被类型的铁线莲'薇安'进行了转录组测序,构建了铁线莲的cDNA文库,并对转录组数据进行过滤及拼接、功能注释、基因的差异表达分析等。通过转录组测序,一共获得了 146 036条序列,99 652个Unigene基因,长度分布在201-14 638 bp之间,平均长度977 bp。对聚类后得到的Unigene进行同源搜索和功能注释,共9 133个Unigene与七个数据库(NR, GO,NT, Pfam, KEGG, KOG, Swissprot)均具有同源性,占总数目的9.16%。经SSR位点查找与分析,从检测到的99 652个序列总数中识别出23 035个SSR位点,包含序列的SSR的数量有18 243个,总长为129 654 975 bp 该研究结果可为铁线莲'薇安'不同花被类型的差异基因表达及筛选和探究重被花形成的分子机制提供了科学依据。     

基于高通量测序的天目铁木转录组分析

天目铁木为中国特有的珍稀濒危植物,5株原生株分布在浙江西天目山.随着环境条件的不断恶化,对濒危植物的保护变得尤为重要,从20世纪80年代初开始对天目铁木开展研究,然而到目前为止对于天目铁木分子生物学和基因功能方面的研究仍旧较少.为了获得天目铁木的基因数据和初步了解基因功能,本研究利用Illumina Hiseq高通量测...     

柳杉维管形成层及木质部区转录组序列及基因表达分析

柳杉是中国南方重要的针叶用材树种,具有树形高大,纹理直等特点,已广泛用于板材和建筑生产中,研究柳杉木材形成过程中维管形成层及木质部区转录组特征,为木材形成主要过程的分子机理,木材形成有关基因调控,培育优良材质的林木提供理论参考。本研究以柳杉形成层组织为材料,通过Illumina Hiseq测序平台对4个不同发育阶段的柳杉维管形成层进行转录组测序;对Unigene进行了蛋白功能注释、分类及KEGG代谢通路分析等。测序数据参照无参转录组测序流程进行分析,一共产生105.9 Gb的数据,过滤后Clean reads均达到90%以上。Clean reads经Trinity软件进行组装,一共产生64 969个Unigene,对Unigene进行拼接,拼接后的Contigs序列有29 381条、Singletons有35 588条,总长度为83 003 836 bp。将Unigen与七大功能数据库(NR, NT, GO, COG, KEGG, Swissprot, Interpro)进行比对,共有42 836个Unigene得到注释,占总Unigene的66.05%。GO注释共有89 644个转录本得到注释。对GO注释转录本进行分类,其中有38 432个转录本(42.87%)注释为生物过程,有32 749个转录本(36.53%)注释为细胞组分,有18 463个转录本(20.6%)注释为分子功能。KEGG注释有31 580个Unigene得到注释。分为6大类为:细胞代谢过程、环境信息处理、遗传信息处理、人类疾病、代谢作用、生物体系统。其中代谢作用涉及的Unigene最多达18 580个,占KEGG总注释Unigene的58.83%。另有37 762个Unigene得到COG注释,被分为25类,其中注释数目最多的类型是仅预测一般功能的Unigene,占总数的16.9%;其次是转录,比例为8.8%;再次是复制、重组及修复,占总数的7.8%。本研究最后分析了可能参与木材形成重要功能基因,为探讨木材发育的分子机制及进行重要功能基因的挖掘提供了提供了依据。     

扁豆荚皮中花青素合成的转录组分析

扁豆在中国广泛栽培,紫色扁豆品种的豆荚中的花青素含量高于绿色品种,但扁豆中花青素合成的调控机制尚待研究。鉴于扁豆基因组测序工作尚未完成,为了解析扁豆花青素合成的调控机制,本研究利用Illumina转录组测序技术,获得了假定Unigene 94982个。与紫色扁豆相比,绿色品种的荚皮中上调Unigene1065个,下调3086个。通过GO、KEGG基因富集分析以及20个表达差异最大的Unigene功能预测,可以发现,绿色扁豆品种抗锈病能力较强,可能与抗氧化能力、氧化还原反应、40S核糖体小亚基(与翻译有关)蛋白高于紫色品种有关,而扁豆豆荚颜色由紫变绿的原因与生长素类AUX蛋白、类衰老特定的半胱氨酸蛋白酶等编码基因的表达水平下调有关。本研究将为解析扁豆花青素合成的调控机制提供借鉴。