基于RNA-seq的能源植物芒转录组分析

芒(Miscanthus sinensis Anderss)是多年生C4草本植物,可为能量和纤维素产品生产提供高品质的木质纤维素材料,是一种理想的能源植物。采用Illumina Hi Seq?2000高通量测序技术,对芒花芽和叶芽进行转录组分析。经拼接组装共获得98 326个Unigene,序列平均长度822 bp,N50为1 337 bp。将Unigene序列与NR、NT、Swiss-Prot、KEGG、GO和COG数据库进行比对(Evalue1e-5),共有74 134条Unigene获得了基因注释,占总Unigene的75.40%。其中,通过GO功能分类,45 507个Unigene映射到GO不同的功能节点上;通过KEGG pathways分析,共有36 710个Unigene参与了128个代谢通路;比对到同源序列比例最高的物种分别为高粱(37 731,60.86%)、玉米(16 258,26.22%)、水稻(3 065,4.94%),共占所有同源序列的92.02%。此外,获得了芒C4关键酶相关基因24个。这些注释信息的完成为芒功能基因及相关候选基因的发掘提供了重要依据。     

朝仓花椒幼芽转录组测序数据组装及基因功能注释

采用2代Illumina Hi-Seq测序技术对朝仓花椒雌株结果期幼芽的转录组进行测序,构建了转录组数据库,获得49 672 080条Clean Reads数据,包含总长度为4 967 208 000 nt序列数据信息;经拼接组装,获得转录组基因信息长达55 902 243 nt的176 407个Contig片段,再通过进一步拼接,共获得平均长度为878 nt的96 475个Unigene片段。与Nt、Nr、Swiss-Prot、COG、GO、KEGG等数据库进行BLAST信息比对(E-value艽10-5),共获得68 315个注释基因,以及62 150个表达序列标签(EST)。与NR数据库比对,发现花椒转录组基因序列与同科柑橘属的甜橙(Citrus sinensis)和克莱门柚(Citrus clementina)具有较高的同源性,分别为43.96%及41.86%,与其他物种的同源性较低,均不足10%;可将花椒转录组的Unigene的功能通过与COG数据库进行注释比对划分为25类;根据GO数据库的注释,可将其分为生物过程、细胞组分和分子功能3大类共55分支,根据与KEGG数据库的比较分析发现朝仓花椒的转录组数据中含有代谢通路相关基因128类,其中包括多种化合物的代谢途径和次生代谢产物的生物合成途径,其中有较多的次生物质代谢途径,如萜类化合物生物合成途径、黄酮类生物合成代谢途径以及花青素的生物合成途径等。     

锥栗种仁转录组及淀粉和蔗糖代谢相关酶基因的表达分析

本研究采用高通量测序技术对淀粉积累高峰期的锥栗种仁进行转录组测序分析,对得到的Unigene进行功能注释、分类及代谢通路分析,并进一步通过实时定量PCR方法分析了7个与淀粉和蔗糖代谢相关酶基因在锥栗种仁发育过程中的表达特征。结果显示:de novo组装后共获得53629条Unigene序列,序列平均长度为746 bp;通过与其他核酸、蛋白质数据库的Blast搜索比对,共26739条Unigene序列获得基因注释,占All-Unigene的49.86%;与COG数据库比对后将其注释的14413条Unigene序列划分成25类;GO功能注释的33926条Unigene基因共分成细胞组分、分子功能和生物过程3大类58个分支;与KEGG数据库比对结果注释的5277条Unigene序列划分为116条代谢通路;7个与淀粉和蔗糖代谢相关酶基因表达量变化趋势中,CH.29636(淀粉合成酶,SS),CH.11971(淀粉分支酶,SBE)两个基因随着种仁的发育呈现逐渐上升的趋势,而其余5个基因(CH.31302、CH.33690、CH.19238、CH.30128、CH.13088)表达量随着种仁的发育呈先上升后下降的趋势,该结果与锥栗种仁发育期淀粉和蔗糖的变化规律一致。这些信息为锥栗果实品质形成相关功能基因的研究提供了重要依据。     

基于高通量测序的阿尔泰蝠蛾(鳞翅目:蝙蝠蛾科)幼虫转录组分析

本研究采用Illumina HiSeq TM 2500测序平台对阿尔泰蝠蛾Hepialus altaicola Wang幼虫进行转录组测序及生物信息学分析.经序列拼接后共获得100133个Unigenes,总长度86319112 bp,平均长度862 bp,N50长度1628 bp.将Unigenes与NR、COG/KOG、Pfam、Swiss-Prot、GO、KEGG数据库比对,共获得38198条Unigenes,其中Nr数据库注释的Unigenes最多,为32381条,占32.34％.通过GO功能分类,共有13216个Unigenes在GO数据库中细胞组分、分子功能和生物学过程等3大类57个分支中找到注释;KEGG通路分析,共有15058条Unigenes被注释,归属于305条代谢通路.CDS预测发现54002条序列可被编码,占全部基因的53.93％.基因注释进一步获得311个与冷适应相关的代谢调节基因,并用FPKM值对基因表达量进行评估.本研究获得的转录组信息及分析结果,为进一步研究阿尔泰蝠蛾的基因功能及低温生态适应性奠定分子基础.     

泡桐维管形成层及木质部区的转录组分析

采用Illumina HiSeqTM2500高通量测序技术,获得泡桐维管形成层及木质部区组织转录组的109021918条clean reads(16.35 Gb)。将clean reads从头组装得到104432个单基因簇(Unigene),平均长度662 nt。将组装得到的Unigenes与公共数据库进行序列比对,分别有40789(Nr:39.05%)、31675(NT:30.33%)、15539(COG:14.87%)、29168(GO:27.93%)、16316(KEGG:15.62%)、30499(SwissProt:29.20%)以及28828(Pfam:27.6%)个Unigenes获得功能注释。通过与GO数据库的比对分析,注释的29168个Unigenes归于生物过程、细胞组分及分子功能三大类的55个功能组;15539条Unigenes注释到COG数据库,被分为25个类别中;基于KEGG数据库可将16316个Unigenes归于130个代谢途径。此外,在泡桐的维管形成层及木质部区转录组中共检测出16118个简单序列重复(SSR)位点。为进一步挖掘泡桐重要功能基因提供了大量数据。     

樟树5种化学类型叶片转录组分析

樟树(Cinnamomum camphora )是樟科植物的一个代表种, 具有材用、药用、香料、油用和生态环境建设等多种用途。叶精油中富含利用价值极高的樟脑、芳樟醇、1,8-桉叶油素、异-橙花叔醇和右旋龙脑等萜类化合物。依据叶精油中主要成分的种类和含量, 可将樟树划分为脑樟、芳樟、油樟、异樟、龙脑樟5种化学类型。文章采用Illumina HiSeq™ 2000高通量测序技术, 对5种化学类型叶片转录组进行测序, 对测序得到的所有Unigene进行GO(Gene Ontology)、COG(Clusters of Orthologous Groups)和KEGG(Kyoto Encyclopedia of Genes and Genomes)分类, 给出功能注释和Pathway注释, 并预测Unigene蛋白编码区(Coding sequence, CDS)。De novo组装共获得156 278个Unigene, 序列平均长度584 bp, N50(覆盖50%所有核苷酸的最大Unigene长度)为1 023 bp。通过与其他核酸、蛋白数据库的Blast搜索比对, 共有55 955条Unigene获得了基因注释, 占所有Unigene的35.80%。其中, 有24 717条Unigene得到GO注释, 有21 806条Unigene得到COG注释。KEGG pathways分析结果表明, 共有3 350条基因(10.19%)注释到次生代谢生物合成途径, 其中参与单萜、二萜、倍半萜和萜类骨架合成的Unigene有424个。在单萜合成的代谢通路中, 有9条Unigene可能编码芳樟醇合成酶基因, 且表达分析结果显示, 芳樟醇合成酶基因在芳樟化学类型中优势表达, 在油樟化学类型中表达水平较低。这些注释信息的完成为樟树功能基因及相关候选基因的发掘提供了基础数据和重要依据。     

榆树叶片形成虫瘿过程中转录组学分析

榆瘿蚜取食侵染榆树叶片形成了榆树虫瘿,本研究采用新一代的高通量Illumina Hi Seq^TM 2000技术测序平台对榆瘿蚜取食刺激的榆树叶片进行转录组测序和功能注释,利用生物学方法对基因表达和功能进行研究。测序获得23.19 Gb碱基序列信息,通过对测序数据进行序列过滤、拼接和去冗余,共获得102 017个Unigenes,通过NR与BLAST等数据库比对,其中有37 899个(37.15%) Unigense被注释。利用KOG、GO、KEGG等数据库对榆树虫瘿叶片的Unigense进行比对,按功其能将匹配的Unigenes基因划分25大类;GO注释将信息归纳为基因的3大主类,57个亚类;以KEGG数据库为参考,将Unigene定位到110个不同的代谢通路,包括氧化应激防御、植物激素信号转导、碳水化合物以及次生物代谢等代谢相关的Unigenes通路。本研究通过二代高通量转录组测序技术研究榆瘿蚜侵染下榆树虫瘿的相关基因,为今后研究榆瘿蚜侵染榆树叶片形成虫瘿的分子机理提供了基础资料。     

文冠果果实转录组测序及分析

为了研究文冠果果实转录组中功能基因表达情况,采样RNA-seq技术对文冠果不同发育时期果实进行测序分析。结果显示:两个时期样品材料测序组装后共获得68 298个Unigene序列,有24 691个Unigene得到注释,占36.15%;GO数据库注释显示14 766条Unigene分为细胞组分、分子功能及生物学过程3大类54个功能组;KOG数据库中注释到的13 994条Unigene功能系统分为25类;以KEGG代谢途径数据库为依据,可将8 284个文冠果果实转录组Unigene分为127个代谢通路;在文冠果果实转录组中发现11 732个SSR位点,其中最多的为单核苷酸SSR,占60.85%。本研究为文冠果果实分子生物学研究提供了一定的基础和参考。     

多年宿根甘蔗的转录组分析及差异基因挖掘

甘蔗(Saccharum hybirds)宿根性直接关系到甘蔗生产成本和种植效益,品种、环境和栽培措施均会影响甘蔗宿根能力,但品种的种性是影响宿根性最关键因素。国内外从分子生物学层面解释甘蔗宿根性差异的文献鲜见报道,从转录水平分析不同宿根年限GR2号和ROC22号基因表达的差异。结果表明：转录组测序共得到100 558条转录本和25 582条Unigene。获得53 790条Unigene的注释结果,分别在NR、Swiss-Prot、KEGG、COG、KOG、GO和Pfam数据库进行比对。GO功能注释共分成三大类,51小类,6年宿根蔗注释到1 029个差异基因,3年宿根蔗注释到3 391个差异基因。主要KEGG代谢通路有8条,分别涉及植物激素信号转导,淀粉和蔗糖代谢,甘氨酸、丝氨酸和苏氨酸代谢,苯丙素生物合成,同源重组,DNA复制,错配修复及植物病原体相互作用。筛选出脱落酸(ABA)相关差异基因6个,分别为脱落酸不敏感蛋白2基因(ABI2)、bZIP转录因子超家族蛋白、碱性亮氨酸拉链型转录因子基因(ABI5)、aba响应元件结合因子1基因(ABF1)、G-box结合因子基因(GBFs)...     

基于转录组测序的扎龙湿地野大麦耐盐性分析

通过转录组测序技术对扎龙湿地野大麦幼苗在200 mmol/L NaCl胁迫2 d后的叶片进行分析。对测序结果进行De novo拼接后,将差异表达基因在GO、KEGG数据库中进行比对注释。结果表明:NaCl胁迫2 d后,对照组与处理组分别获得Unigene序列61038个和46754个,检测到差异表达基因25465个,差异基因GO功能注释到3个大类的55个功能组。差异表达基因被注释到135个Pathway上,直观地显示出NaCl胁迫下野大麦幼苗体内发生调节及改变的代谢过程和信号通路,其中主要涉及光合作用、脯氨酸代谢、叶绿素代谢等途径。通过对野大麦幼苗叶片转录组分析,发掘相关耐盐基因,有助于培育野大麦及其近缘作物新品种,并对揭示植物耐盐性分子机制及相关代谢途径具有重要意义。     

基于RNA-Seq技术的蜡蚧轮枝菌转录组分析（英文）

【目的】蜡蚧轮枝菌是一种防治植物病原物和害虫的生防菌,弄清其致病的分子基础具有重要的意义。【方法】利用Illumina HiS eq技术测序蜡蚧轮枝菌的cD NA文库,并进行RNA-Seq序列拼接和功能注释。【结果】共获得1634670条高品质reads,碱基GC含量48.50%。14856条Unigene经组装拼接后,在Nr、Swiss-Prot、COG和KEGG数据库中分别比对得到1467、9379、6518和4188条Unigenes。采用GO分类工具将21403条Unigenes分成24个功能组,主要包含在细胞组成(1369)、分子功能(1679)和生物学过程(1928)3个大类里,在COG数据库注释的基础上,把6518条Unigenes分成38个功能组,通过GO和COG丰度识别分类,聚类的主要是一般功能预测、次生代谢产物的合成、运输和分解代谢以及信号转导机制起作用的独立基因。能与KEGG路径匹配的108条Unigenes中大部分和新陈代谢途径及次生代谢产物的合成有关。【结论】这些结果将有助于找寻蜡蚧轮枝菌的致病因子,从而丰富其致病机理。     

基于高通量测序的玉米象转录组分析

玉米象Sitophilus zeamais是一种世界性的储粮害虫,但其基因信息尚不完善.本研究利用高通量测序平台Illumina HiSeq TM2000对玉米象成虫进行转录组测序,总计获得64358条unigenes,总长度39481752 bp,最短201 bp,最长29046 bp,平均长度613 bp.将unigenes序列在7大数据库Nr、GO、Swissprot、KOG、Pfam、KEGG、TrEMBL中进行注释,分别有25271、20747、16477、12060、9471、3370、25500条unigenes获得注释.通过GO功能分类,共有20747条unigenes在GO数据库中生物学过程、细胞组分、分子功能3大类68个亚类功能组中找到对应.KOG结果显示,12606条unigenes归到25个基因家族.通过KEGG代谢通路分析,共有3370条unigenes被注释,分别归属于细胞进程、环境信息进程、遗传信息进程、新陈代谢和有机体系统5大类代谢途径,共33个亚类274个功能通路.进一步的数据分析,共鉴定出146081个SNP位点和5002个简单序列重复(SSR)位点.本研究获得的玉米象转录组信息,为玉米象的功能基因挖掘提供了重要的信息资源.     

Transcriptome analysis of Verticillium lecanii based on RNA-Seq technology

[目的]蜡蚧轮枝菌是一种防治植物病原物和害虫的生防菌,弄清其致病的分子基础具有重要的意义。[方法]利用Illumina HiSeq技术测序蜡蚧轮枝菌的cDNA文库,并进行RNA-Seq序列拼接和功能注释。[结果]共获得1634670条高品质reads,碱基GC含量48.50%。14856条Unigene经组装拼接后,在Nr、Swiss-Prot、COG和KEGG数据库中分别比对得到1467、9379、6518和4188条Unigenes。采用GO分类工具将21403条Unigenes分成24个功能组,主要包含在细胞组成（1369）、分子功能（1679）和生物学过程（1928）3个大类里,在COG数据库注释的基础上,把6518条Unigenes分成38个功能组,通过GO和COG丰度识别分类,聚类的主要是一般功能预测、次生代谢产物的合成、运输和分解代谢以及信号转导机制起作用的独立基因。能与KEGG路径匹配的108条Unigenes中大部分和新陈代谢途径及次生代谢产物的合成有关。[结论]这些结果将有助于找寻蜡蚧轮枝菌的致病因子,从而丰富其致病机理。     

圆口铜鱼亲本与子代肝脏转录组特征及对比分析

通过Illumina HiseqTM 2000测序平台首次开展了圆口铜鱼(Coreius guichenoti)亲本与子代肝脏转录组测序并对比分析其测序结果。对转录组进行拼接和组装, 共获得80688个Unigene, Unigene的长度主要分布在401—600 bp, 占总Unigene的43.56%。与Nr、GO、COG、KEGG等公共数据库比对并进行功能注释, 经分析圆口铜鱼亲本与子代差异表达基因共2701个, 其中上调1240个, 下调1461个。差异基因表达模式聚类热图显示, 同一样本不同重复基因表达相似。将差异基因映射到代谢通路KEGG数据库进行富集分析, 并对37个KEGG显著富集的代谢通路(P<0.05)构建其基因共表达网络, 结果显示丙酸盐代谢、丙酮酸代谢、原核生物固碳途径、乙醛和二羧酸代谢、脂肪酸代谢、萜类骨架生物合成是KEGG差异表达基因的核心代谢途径, 且这些基因在圆口铜鱼子代中表达量均显著上调。该结果丰富了圆口铜鱼基因组数据, 并首次从分子水平比较了圆口铜鱼亲本和子代的代谢差异, 同时解释了在同一循环水系统中子代不易患病的现象。     

绿色杜氏藻转录组分析

为了深入了解绿色杜氏藻(Dunaliella viridis)基因信息及功能、耐盐相关通路(甘油脂代谢)及关键酶,本文首次通过Illumina HiSeq^TM 2000高通量测序技术对绿色杜氏藻转录组进行测序,利用Trinity软件将数据组装形成转录本,对所有转录本进行COG(Clusters of Orthologous Groups)、GO(Gene Ontology)和KEGG(Kyoto Encyclopedia of Genes and Genomes)分类和功能注释、Pathway注释以及蛋白编码区(Opening reading fragment,ORF)的预测,并对甘油脂代谢通路关键酶基因进行了分析。转录组测序共获得81 593个转录本,其中ORF共有77 117条,约占所有转录本的94.50%。COG分类结果表明,16 569条转录本被分为24个类别。GO分类结果表明,76 436条转录本被注释。在所有注释分类中,生物学过程转录本数量最多,为30 678条,占总转录本数的40.14%。KEGG分析结果表明,317个标准途径中包含26 428条转录本,含转录本最多的类别是代谢,为9949条(37.65%)。与代谢有关的途径为131条,占所有注释途径的41.32%。在甘油脂代谢通路中仅发现1条关键酶转录本(二羟丙酮激酶),该酶可能与绿色杜氏藻耐盐胁迫中甘油的合成有较大关系。本研究进一步完善了绿色杜氏藻的基因信息,为绿色杜氏藻代谢途径研究奠定了坚实的基础。     

基于RNA-Seq的甘蔗主茎和分蘖茎转录组建立及初步分析

本研究以甘蔗与斑割复合体杂交F1代中1个株系在主茎和分蘖茎的节间长度存在明显差异的植株为实验材料,利用RNA-seq测序技术对甘蔗主茎(节间短)和分蘖茎(节间长)叶片样品RNAs进行了转录组测序并进行从头组装分析。结果表明,共获得测序数据量11.16Gb,其中甘蔗主茎(BG-stalk)为5.67Gb,分蘖茎(BG-tiller)为5.49Gb。去冗余后,denovo组装得到69204条Unigene,并对这些Unigene进行七大功能数据库(NR,NT,GO,COG,KEGG,Swissprot和Interpro)注释,结果发现,有83.73%(57942条)的Unigene得到注释,16.27%(11262条)的Unigene未被注释。所有的Unigene共有9156个SSR(simple sequence repeat)位点,其中三核苷酸重复最多、四核苷酸重复最少,且CCG/CGG出现的频率最高。差异表达基因分析显示,分蘖茎(BG-tiller)样品表达的上调基因有1842个,下调基因的有2663个。进一步对这些差异表达基因进行GO和Pathway功能分类分析,分别获得57个功能小组和19条生物通路。本研究对甘蔗主茎和分蘖茎叶片转录组信息进行初步分析,获得了分蘖茎表达上调和下调的差异基因,为后面进一步分析调控甘蔗节间长短差异的基因及挖掘甘蔗分蘖生育调控相关基因提供数据参考。     

基于转录组测序的滇山茶花叶呈色机理分析

该研究采用Illumina Hi-Seq2500高通量测序技术,对滇山茶的花叶(绿叶区和金叶区)的cDNA进行了转录组测序,分析滇山茶的花叶呈色机理。测序结果得到了228 862条单基因序列,筛选得到了4 851条差异表达基因。2 291条差异表达基因在GO数据库中有功能注释,1 826条差异表达基因被注释到COG数据库,这些差异表达基因被分为23个主要功能类别,其中大部分基因与滇山茶的生长和代谢有关。1 363条差异表达基因注释到了KEGG数据库,进一步筛选出了差异表达基因注释序列最多的15条代谢通路,其中卟啉和叶绿素代谢途径中的GluRS、δ-ALAD、ALAS基因和类黄酮合成途径中ANS、LAR、CHS基因与滇山茶的花叶形成密切相关。研究表明,滇山茶的花叶呈色是由于叶色相关代谢通路中大量的基因表达量发生变化引起的,并不是由单基因突变所致,病毒诱导的多个代谢途径基因沉默可能是滇山茶花叶呈色的主要原因。该研究丰富了滇山茶的转录组信息,并为滇山茶的彩叶品种培育提供了一定的参考。