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

本研究基于新一代高通量测序技术平台Illumina Hi SeqTM4000对野百合进行转录组测序,对物种的转录组序列进行统计,并将得到的数据进行de novo组装,结果共获得47 605条Unigenes,总长度为33 972 306 bp,平均长度为713 bp,N50为1 204 bp。将获得的Unigenes与Nr、Swiss-Prot、KEGG以及KOG数据库进行比对,结果显示,分别有28 104、17 739、11 984及14 682条Unigenes成功注释。通过与KOG数据库进行比对,可分为25个不同的功能注释。与GO数据库进行比对,结果显示,共有33 254条Unigene获得注释,这些功能注释分为三大类50个功能亚类。其中,生物过程最多。以KEGG数据库参考,共有11 984条Unigenes参与133条代谢途径分支,以代谢相关的通路较为集中,找到了与花青素合成关键酶的Unigenes。本研究极大地丰富了野百合的基因资源,为进一步开展野百合功能基因及分子标记育种等方面的研究提供了一定理论支持与依据。     

四川道地中药金龙胆草叶片转录组特性研究

本研究以道地药材金龙胆草叶片为研究对象,利用Illumina Hi Seq 2500高通量测序技术构建了金龙胆草转录组数据库,获得了42 903 527条Reads数据,通过与Nr、GO、COG、KOG、KEGG等数据库比对,最终获得了38 199个具有注释信息的Unigenes。其中以KEGG数据库为参考,依据代谢通路将Unigenes分成117类,包括萜类骨干合成、二萜类合成、黄酮类化合物生物合成及聚糖生物合成等路径,分别有99,23,161及70个Unigenes映射到上述途径,这些Unigenes可能参与金龙胆草主要活性成分三萜皂苷、特征成分苦蒿素、金龙胆草黄酮及金龙胆草多糖的生物合成。金龙胆草转录组测序工作的完成,极大地扩充了金龙胆草的基因资源,为药用功能基因的发掘与利用、遗传改良及有效成分含量的提高等研究奠定基础。     

白鲜种子休眠解除过程中的转录组分析

以解除休眠过程不同阶段的种子为材料,通过高通量测序技术对样品转录组进行测序分析。共获得135.98 Gb的原始数据,302 453条Unigenes,注释到七大功能数据库(GO、KEGG、KOG、NR、TrEMBL、Pfam、Swiss-Prot),分别有139 276(NR:46.05%)、128 449(TrEMBL:42.47%)、86 498(Swiss-Prot:28.60%)、90 564(KOG:29.94%)、103 481(KEGG:34.21%)、108 540(GO:35.89%)以及90 410(Pfam:29.89%)条Unigenes得到功能注释。通过GO分类和KEGG Pathway富集性分析,分别归于59个GO类别和145条代谢途径。同时,有6 081个Unigenes被注释为转录因子。利用高通量测序技术获得白鲜转录组信息,有助于从分子水平对白鲜进行深入研究。     

药食同源植物甘葛藤的全长转录组分析

为了深入了解甘葛藤转录组的整体水平及黄酮类生物合成通路基因。利用高通量测序PacBio Sequel平台,以甘葛藤根、茎、叶的混合样品为材料,使用单分子长读数测序技术(SMRT)对甘葛藤进行全长转录组测序及分析。平台共获得10 994 967个高质量reads和384 072条全长非嵌合序列(FLNC),测序数据经质控后获得90 856个转录本;获得的所有转录本经NR、SwissProt、KOG、KEGG、GO数据库进行注释和功能分类,结果有85 239个单基因被注释,NR注释数量最多为84 675个,占93.2%;KEGG注释的基因最少,22 330个基因被注释到132条途径,代谢途径分布的基因较多(9 368,41.95%)。预测到3 507个转录因子,bHLH转录因子家族的基因最多。14 127个基因被分配到17个R基因类别,主要为RLP类。检测到33 660个SSR序列,多为AG/CT类型。分析黄酮类生物合成途径,发现与黄酮类合成相关的基因110个,其中,26个编码HCT,3个编码CHS,7个编码CHI。PacBio测序平台能获得更长的转录本,SMRT技术能够深入挖掘甘葛藤转录数据,比第二代测序技术能够获得更高的转录本注释率。在高通量全长转录组水平对甘葛藤进行了研究,为甘葛藤的分子生物学研究提供了较可靠、全面的转录组数据,为进一步开发甘葛藤的分子标记和挖掘优良基因提供了科学依据。     

山东栒子全长转录组测序数据组装及功能注释

山东栒子是中国山东省的特有种,现已处于极度濒危状态。目前,山东栒子的分子生物学研究较少,基因数据库资源极度缺乏,急需探究其生物学遗传信息,以加快对山东栒子的保护遗传学工作。本研究以山东栒子叶片、花、成熟果实为实验材料,利用PacBio Sequel测序平台对其转录组进行全长转录组测序。共得到高质量去冗余的转录本53932个,作为最终转录本序列。预测到的CDS区共52490个;对其SSR位点进行分析,对测序得到Unigenes进行单核苷酸至六核苷酸重复的SSR位点搜索,共搜索到26796个SSR位点。对非冗余转录本利用BLAST软件与NR、Swissprot、GO、COG、KOG、KEGG 6个数据库进行比对,一共成功注释了53319个Unigenes,其中与NR数据库进行比对中注释为苹果的的相关基因数量最多,其次是白梨和桃;在GO数据库中有33305条山东栒子Unigenes被注释分类,由生物学过程、细胞成分和分子功能三部分组成;在与COG数据库比对中,共有23910条比对到了同源序列,且一共被分为25类;在与KEGG数据库的一系列比对中,可将Unigenes映射到126条代谢通路中。本研究在高通量全长转录组水平对山东栒子进行了系统研究,这为进一步开展山东栒子的分子标记开发和挖掘优良基因提供了科学依据,从而推动山东栒子的保护与利用。     

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

采用高通量测序技术平台Illumina Novaseq 6000对滇白前种子进行转录组测序,共获得平均长度为753.9 bp的43663个Unigenes.注释到六大功能数据库(NR,COG,Pfam,Swiss-Prot,GO,KEGG)上的Unigene总数为29177个.滇白前Unigenes比对到NR数据库共有...     

广藿香幼叶与成熟叶片转录组数据组装及基因功能注释

为了进一步阐明广藿香中药用活性成分生物合成的分子机制,本研究以海南广藿香幼叶及成熟叶片为材料,采用BGISEQ-500高通量测序平台进行转录组测序,分别获得了63 751 826条和65 949 390条clean reads,平均读长为90 nt。De novo组装后将All-unigene分别注释到Nr、KOG、GO、KEGG、Swiss-Prot、Inter Pro数据库,对每个数据库注释的Unigene数目进行统计,共有162 509条Unigene有对应的功能信息,其中105 430条Unigene被注释到Nr数据库,显示与芝麻有69.87%的相似度;有83 369条Unigene被注释到KOG数据库,根据功能将其分为25类;有12 261条Unigene与GO数据库中的基因具有相似性,将其归为3大类中49个功能组;有79 053条Unigene被注释到KEGG的代谢通路中,分属于124类代谢通路,包括次生代谢物质生物合成、倍半萜和三萜类化合物生物合成、黄酮和黄酮醇生物合成、花青素生物合成等。该研究结果对广藿香药用活性成分生物合成与代谢、关键酶基因克隆以及分子标记开发等研究有一定的帮助。     

基于高通量测序技术的杭白芷(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高通量测序技术进行转录组测序,并对得到的测序数据进行de novo组装,获得97 503条Unigenes,N50为1 010 bp、平均长度852 bp。与公共数据库进行比对,注释到NR、Swiss-Prot数据库的Unigenes分别为60 926、44 249条。将Unigenes与COG数据库比对,有42 344条Unigenes成功注释,根据功能大致分成25类;与GO数据库比对,有15 947条Unigenes获得注释,按功能可分为生物学过程、细胞组分和分子功能3大类55亚类,其中参与的生物学过程较多;以KEGG数据库参考,有26 267条Unigenes参与330条代谢途径分支,以代谢相关的途径较为集中;差异基因表达分析我们得出,根与茎的差异表达基因为2 826个,根与叶之间的差异表达基因为3 230个,这两组差异基因相对校多,除此之外茎与叶差异基因相对较少;差异基因的GO分类、GO富集以及pathway富集分析表明,主要在细胞组分、蛋白质活性、生物合成以及代谢过程等所占比例较高。这些研究结果极大地扩充了火力楠的基因资源,将有助于火力楠基因的发掘与利用、分子标记的开发及其种质资源遗传改良的研究等。     

云南火焰兰高通量转录组测序分析

云南火焰兰(Renanthera imschootiana Rolfe)是濒临灭绝的附生植物,关于其分子水平的研究少有报道。为了获得云南火焰兰相关基因组基础信息,本研究通过Illumina HiSeqTM4000技术对云南火焰兰叶片进行转录组测序,并利用相关生物信息学数据库对其进行序列比对和数据分析。分析结果表明:共获得77 888条功能基因(Unigenes),其N50为1 348 bp,平均长度为805 bp,Q20和Q30序列分别占98.04%和94.32%。经比对,在Nr、KEGG、KOG、Swiss-prot、GO等数据库中能得到注释的Unigene分别为31 298、29 356、14 852、15 697、17 645条。在GO数据库中,共得到24 473个GO功能,将其分为3个大类和47个亚类,以生物过程这一类别所注释的基因数最多;KOG功能注释信息有17 743个,将其分为25个基因功能大类,其中基因数目较多的分别是一般功能基因和信号传导机制的相关基因;在KEGG中涉及的相关通路可将其归为19个亚类,其中以代谢通路较为突出。此外,对云南火焰兰转录组Unigene进行CDS比对和预测,共预测到2 317条CDS,其大小主要集中在200～500 bp,占总CDS的79.41%。本研究可为云南火焰兰乃至火焰兰属植物功能基因的挖掘、利用以及遗传育种等方面的工作提供一定的科学依据。     

姜黄根茎的高通量转录组测序与生物信息学分析

为获得姜黄(Curcuma longa)的转录组特征信息,本研究采用Illumina HiSeqΧTen高通量测序平台对姜黄根茎进行高通量转录组测序并进行系统的生物信息学分析。共获得7.18Gb Clean数据,组装了50194条unigenes,平均长度961.3 bp,N50为1 339 bp。数据库比对显示,姜黄根茎转录组unigenes在NR、Swiss-Prot、KEGG、KOG、eggNOG、GO、Pfam数据库中分别注释到38 802条(77.30%)、27 869条(55.52%)、14 725条(29.34%)、22 225条(44.28%)、37 317条(74.35%)、25 863条(51.53%)、26 137条(52.07%)。注释结果显示,姜黄与野生型马来西亚蕉的同源序列最多,unigenes在GO数据库中注释到参与生物过程、细胞组分和分子功能3个大类50小类,KOG功能分类获得25个不同的功能群,涉及128个KEGG代谢通路,其中包括21个次生代谢通路。在植物抗性基因(PRG)数据库中分别注释到3 718条unigenes;借助MISA软件发现7 183...     

不同种质油茶的比较转录组分析

本研究采用高通量Illumina转录组测序技术对‘万海4号’(WH4)、‘博白大果’(BBDG)、‘长林40’(CL40)、‘岑软’(CR)和‘红花油茶’(HH)五个不同种质的油茶进行转录组测序并进行差异比较分析。结果表明,WH4、BBDG、CL40、CR以及HH通过组装分别获得了71 852、111 693、75 436、15 899和96 205个Unigenes,将所有Unigenes与公共数据库(Nr, Swiss-Prot, KOG, KEGG)进行相似性比对并注释,五种油茶注释基因所占百分比依次为65.51%、46.41%、60.57%、91.94%、54.98%;KEGG富集分析表明WH4、BBDG、CL40和HH中基因丰度最高的都是核糖体途径,而CR中基因丰度最高的是碳代谢途径;GO富集分析显示,WH4、BBDG、CL40、CR和HH基因分布结果一致。结合KEGG、GO和KOG富集分析,绘制出了油茶三萜皂苷生物合成途径通路图,筛选出该途径中的关键酶及相关基因,并对其进行了差异分析。该研究结果完善了油茶三萜皂苷生物合成途径相关基因的表达特征,可为油茶三萜皂苷的生物合成提供...     

白刺花种子转录组分析

为探索白刺花(Sophora viciifolia)硬实形成的相关机制,采用Illumina Hi Seq TM 2000高通量测序技术对白刺花种子转录组进行测序,利用Trinity软件将数据组装形成转录本,对所有转录本进行Nr(NCBI nonredundant protein sequences)、Nt(NCBI nucleotide sequences)、Pfam(protein family)、KOG/COG(eu Karyotic ortholog groups/clusters of orthologous groups)、Swiss-Prot(A manually annotated and reviewed protein sequence database)、KEGG(kyoto encyclopedia of genes and genomes)和GO(gene ontology)分类和功能注释、Pathway注释,并对种子形成的代谢通路中的相关基因进行了分析。转录组共获得了333 339 724条初始序列,总长为335 557 bp,初始序列组装获得序列片段的平均长度与N50值分别为282 bp和537 bp;与KOG功能注释、GO分类及KEGG代谢通路分析后,获得了44 840个GO功能注释、46 126个KOG功能注释以及89 494个PFAM注释:并从KEGG通路中找到有色氨酸代谢、半胱氨酸和甲硫氨酸的代谢途径的编码基因片段分别有66和37个。     

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

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

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

猪毛菜不同部位转录组特征分析

本研究以猪毛菜的地上部分和地下部分为实验材料,通过转录组测序分析探讨猪毛菜不同组织部位的基因表达特征和代谢物质形成有关的基因调控。测序数据分析结果显示获得192 777条转录本序列。差异表达筛选得到14 424个差异基因,其中上调表达基因6 137个,下调表达基因8 287个;GO功能富集分析得到14 424个差异基因。差异基因与KEGG数据库进行比对,共有5 169个差异基因注释在138个通路上,其中参与次生代谢途径共11条通路呈显著性富集。进一步分析黄酮类化合物合成和异喹啉生物碱生物合成通路中的差异基因和其对应的相关酶。随机选取了黄酮类化合物合成通路上的8个基因,通过RT-qPCR分析验证了转录组数据的准确性。本研究利用高通量测序技术和生物信息分析获得猪毛菜不同部位的转录组信息特征,为猪毛菜的药用成分研究和基因功能鉴定提供了科学依据。     

药用植物紫红獐牙菜的转录组信息分析

紫红獐牙菜是重要的民族药,本研究采用高通量测序技术对紫红獐牙菜进行转录组测序并分析。结果显示,完成的15个样品的转录组测序,获得112.17 Gb Clean data,各样品Clean data均达到6.18 Gb,Q30碱基百分比在94.96%及以上。组装后共获得47 106条unigene,其中长度在1 kb以上的unigene有18 971条。通过与KEGG、GO、KOG、COG等多个数据库进行比对,对unigene进行功能注释,共获得35 375条unigene。GO数据库中注释到的26 570条unigene,可分为生物学过程、细胞组分和分子功能三大类共44个亚类。以KEGG数据库为参考,21 321条基因被注释,参与的代谢通路分为5大类,分别是细胞过程、环境信息处理、遗传信息处理、代谢和有机系统,其中与代谢相关的通路最多,约占所用通路的72.5%;KEGG代谢通路分析发现紫红獐牙菜中1 147条unigene参与到22个次生代谢标准通路中,有77条unigene参与编码环烯醚萜苷类合成通路中的25个关键酶。此外,紫红獐牙菜转录组中找到9 027个SSR重复位点,6种SSR...     

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

为了获得珍贵用材树种大花序桉顶芽转录组数据及预测关键基因功能,本研究基于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,重复基元类型丰富,标记开发潜力大。本研究利用高通量测序获得丰富的顶芽转录组信息,可以为大花序桉分子辅助育种提供丰富的资源。     

基于高通量测序的露地菊(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个代谢通路,包括核糖体途径、植物激素信号传导途径、淀粉蔗糖代谢、碳代谢、氨基酸的生物合成等。本研究获得的转录组数据将有助于揭示露地菊生长发育及耐盐胁迫的应答机制和分子基础,及相关抗性基因的挖掘和分子辅助育种等方面的研究。     

桐花树叶片的转录组分析

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