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

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

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

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

基于高通量测序的夏黑葡萄叶片衰老转录组分析

为探究夏黑葡萄叶片衰老相关基因的表达情况及代谢通路,本研究以叶果比12:1修剪的夏黑葡萄为材料,采用高通量测序技术对采集花后45、65、85和105 d的叶片进行转录组分析。结果显示,花后45 d的样品中得到50 855 290条clean reads,GC含量为46.49%;花后65 d的样品中得到50 432 588条clean reads,GC含量为46.13%;花后85 d的样品中得到52 169 372条clean reads,GC含量为47.10%;花后105 d的样品中得到50 988 842条clean reads,GC含量为47.01%。与参考基因组进行序列比对,发现各样品74.20%以上的reads比对到参考基因组,绝大部分分布在外显子区,其中有72.43%的reads比对到参考基因组的唯一位置。各样品间基因差异表达分析得到共有差异基因2 583个。GO功能注释分析发现这些差异基因主要参与生物过程中的光合作用,以及与核酸结合转录因子活性等分子功能相关;SNP位点分析结果发现在每个样品中平均有417 490个变异位点,主要分布在外显子区、内含子区、基因上下游区域。本研究明确了夏黑葡萄叶片衰老过程中相关基因的变化,为进一步挖掘夏黑葡萄叶片衰老过程中关键基因的功能奠定了理论基础。     

基于转录组测序分析的黄精种子休眠解除相关差异基因研究

旨在获得黄精转录组数据库并挖掘参与其种子发育和休眠解除相关基因,以休眠解除前后的黄精种胚为试材,利用新一代高通量测序手段对供试样品进行转录组测序,并进行系统的生物信息学分析.黄精种子休眠解除前后样品中共得到79716个差异表达基因,上调的表达基因有60074个,下调的表达基因有19642个.休眠解除前后的黄精种胚中共有...     

马尾松转录组测序和分析

本研究首次构建了马尾松均一化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。     

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

为了获得甘薯耐盐转录组序列信息,挖掘差异表达基因及其相关代谢途径,本文以盐胁迫处理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代谢通路。以部分差异表达基因的转录组数据为基础,...     

基于高通量测序的露地菊(Chysanthemum×grandiflora)盐胁迫转录组分析

为揭示露地菊生长发育及耐盐胁迫的应答机制和分子基础,本研究以盐胁迫处理的露地菊及其对照为材料,使用Illumina Hiseq2500高通量测序平台对转录组进行测序,分别获得了60370448和71415448条Clean reads,通过序列拼接组装得到45591条Unigene,平均长度724 bp。有37675条Unigene在七大数据库(COG,GO,KEGG,KOG,Pfam,Swiss-Prot,NR)中得到注释。通过比对露地菊盐胁迫处理组和对照组样品间Unigene的表达量及在各数据库中的注释情况,统计得到:有4143条差异表达基因获得注释;有2441条差异表达基因在GO数据库中获得功能注释;注释到COG数据库中的2281条差异表达基因依据功能可分为25类;有1062条差异基因映射到KEGGPathway数据库中,涉及了199个代谢通路,包括核糖体途径、植物激素信号传导途径、淀粉蔗糖代谢、碳代谢、氨基酸的生物合成等。本研究获得的转录组数据将有助于揭示露地菊生长发育及耐盐胁迫的应答机制和分子基础,及相关抗性基因的挖掘和分子辅助育种等方面的研究。     

基于高通量测序技术的杭白芷(Angelica dahurica)根转录组数据分析

为获得杭白芷转录组信息特征,本研究利用Illumina HiSeqⅩTen测序平台对杭白芷根进行高通量转录组测序,获得高质量序列(Clean reads)47742445条,Trinity denovo组装后得到47044条Unigenes,平均长度1164.20 nt。BLAST分析显示分别有32208(68.46%)、23049(48.99%)、10479(22.27%)、17883(38.01%)、28201(59.95%)、20731(44.07%)、55(0.12%)条Unigenes在数据库NR、Swiss-Prot、KEGG、KOG、eggNOG、GO、Pfam中获得注释,可归为GO分类的生物过程、细胞组分和分子功能3大类57分支,涉及205个KEGG代谢通路,其中包括27个次生代谢通路。蛋白编码框序列32303个,高等植物转录因子58个家族,借助MISA软件发现10020个SSR,其中二碱基重复最丰富,有4336个,出现频率为43.27%;五碱基重复SSR最少仅占0.37%。本研究获得了大量基因序列信息以及SSR信息,为今后开展相关分子机制研究提供了数据资源和理论基础。     

桐花树叶片的转录组分析

为了进一步了解桐花树耐盐的分子机制,本研究采用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个家族。本研究结果为桐花树的基因功能分析、分子标记开发和遗传多样性研究提供一定的理论参考。     

基于RNA-seq牡丹SSR标记开发及通用性分析

从牡丹花色候选基因中开发一套SSR标记,为牡丹花色分子遗传学研究和分子标记辅助育种提供帮助.本研究基于高通量转录组测序技术(RNA-seq)获得了278条涉及调控牡丹花色形成的Unigenes序列,利用SSRIT在128条Unigenes中检测到215个SSR位点,出现频率为46.04％.其中优势重复基序为二核苷酸、三...     

20个芍药切花品种遗传多样性的AFLP分析

为了探索不同芍药切花品种之间的亲缘关系特征,利用荧光标记AFLP分子标记技术,采用9 对PstI /MseI 引物组合对20 个芍药切花品种进行了DNA遗传多样性分析。结果表明：9 对引物均显示多态性,共产生谱带733 条,其中多态性条带699 条,多态性位点平均为95.36%,Nei 基因多样性指数为1.4163~1.5552,平均1.4740;Shannon 信息指数为0.4197~0.4929,平均0.4490。经UPGMA聚类分析后,20 份芍药品种的Dice 相似系数0.458~0.741,在相似系数0.52 处,可将供试的20 个切花芍药品种分为5类。研究认为,药芍品种间的遗传差异较大,具有丰富的遗传多样性,为芍药品种分类及亲缘关系分析提供参考依据。     

基于转录组测序的苦瓜热激蛋白相关基因分析

为获得苦瓜叶片高温胁迫状态下表达的相关基因,本研究以高温胁迫下苦瓜叶片为实验材料,采用Illumina HiSeq 2500高通量测序技术对正常温度(26℃)和高温胁迫(38℃)下苦瓜叶片转录测序后进行组装,再利用KEGG、Pfam、Swiss-Prot和NR数据库平台对热激蛋白的相关功能基因筛选,共获得24个热激蛋白相关基因,包括得20个热激蛋白基因、2个位于叶绿体中的小热激蛋白以及2个热激转录因子。聚类分析研究表明,涉及到热激蛋白相关的20个Unigene广泛分布于3大类中,其中c42567和c43094这两个亲缘关系最近;两个小热激蛋白c42963和c12133的亲缘关系次之;c35683和c3040两个热激转录因子亲缘关系最远。通过本次研究中获取的24个热激蛋白相关Unigene信息,有助于深入研究苦瓜耐高温胁迫机理等提供数据信息。     

芍药分生组织决定基因APETALA2(AP2)的克隆及生物信息学分析

花分生组织决定基因APETALA2(AP2)属于植物ABCDE模型基因,在花器官发育过程中起着重要的调控作用。为进一步了解芍药AP2基因的生物学功能,利用RACE扩增和测序技术克隆plAP2基因序列,利用生物信息学在线程序对其序列特征、蛋白结构及功能、亚细胞定位进行预测,并利用MEGA 5.0构建不同植物AP2分子进化树,最后利用qPCR检测其在内外花瓣中的差异表达情况。结果显示,克隆获得芍药AP2基因(plAP2)cDNA序列全长1 935 bp,其ORF全长为1 578 bp,编码525个氨基酸。蛋白结构与功能分析表明,plAP2蛋白为亲水性不稳定蛋白,无跨膜结构和信号肽,表明为非分泌蛋白;核定位信号位于氨基酸序列139-147(KKSRRGPRS);二级结构包括α-螺旋(24%)、β-折叠(19%)、β-转角(28%)和无规则卷曲(28%);plAP2蛋白存在8个糖基化位点和64个磷酸化位点,plAP2蛋白包含2个相同的保守结构域:AP2/(Ethylene-Responsive factors,ERF)(151-213aa和243-306aa)。亚细胞定位主要在细胞质(45.0%)中,少量分布于微体、线粒体基质间隙和溶酶体;进化树分析表明,芍药AP2基因与牡丹高度同源且亲缘关系最近;qPCR检测显示外瓣AP2表达量均极显著高于内瓣(P0.01)。克隆出芍药AP2全长cDNA序列,系统地揭示了plAP2蛋白基本结构、功能位点区域、细胞定位以及组织表达情况,为今后深入研究plAP2基因功能提供基础素材和理论参考。     

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通路。该转录组测序数据质量高,结果覆盖面广,为柳树耐盐基因挖掘和研究提供了一定的理论参考。     

基于转录组测序的向日葵对干旱胁迫及复水差异表达分析

为研究向日葵在不同程度干旱胁迫下基因的表达调控模式,对抗旱自交系(17062)和非抗旱自交系(167072)分别进行干旱胁迫0、4、8、12、16 d以及复水处理,本研究设置干旱胁迫第8天、第16天以及复水处理与各梯度的对照组进行转录组测序分析。共设有6个差异分组,经过测序共得到277.31G的Clean data,Q30碱基>91%,平均GC含量为45.37%。样品间进行差异基因筛选,得到的DEG数量分别为2365、2571、3149、2692、4593、4980。对各差异基因进行GO和KEGG功能注释,GO分析表明,6个比较组参与生物过程差异基因较多,参与细胞组分的次之,参与分子功能的较少;KEGG分析表明干旱胁迫对向日葵的植物激素信号转导、蔗糖与淀粉代谢、氨基酸代谢有较大的影响,对生理特性的测定显示可溶性糖和可溶性蛋白与正常供水相比具有上升趋势,表明在逆境胁迫下,植物会通过对基因的表达调控的调节来改变相关生理生化代谢途径。本研究通过对向日葵叶片进行转录组分析,对揭示抗旱分子机制及相关代谢途径具有重要意义。     

桂花(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位点出现频率高,类型较丰富、多态性潜能较高,为今后桂花分子育种提供了理论依据。