四川山鹧鸪的转录组组装和注释

四川山鹧鸪Arborophila rufipectus是中国特有的珍稀濒危鸟类.本研究对1只成年雄性四川山鹧鸪个体的心脏、肝脏和肾脏进行了转录组测序、组装和注释.其原始序列过滤后分别产生了5.70 G、4.60 G和5.16 G数据.286661条转录本经过Trinity组装并去掉冗余后共得到234488个基因.BUS...     

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

玉米象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)位点.本研究获得的玉米象转录组信息,为玉米象的功能基因挖掘提供了重要的信息资源.     

基于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个。这些注释信息的完成为芒功能基因及相关候选基因的发掘提供了重要依据。     

绿色杜氏藻转录组分析

为了深入了解绿色杜氏藻(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条关键酶转录本(二羟丙酮激酶),该酶可能与绿色杜氏藻耐盐胁迫中甘油的合成有较大关系。本研究进一步完善了绿色杜氏藻的基因信息,为绿色杜氏藻代谢途径研究奠定了坚实的基础。     

闽楠转录组分析及基因功能注释

采用第二代Illumina HiSeq测序技术对闽楠的木质部、韧皮部、叶片进行转录组测序,分别获得Clean Reads片段41383707条、43343922条、44191586条,经转录本拼接后得到序列总长度达120535288 bp的383331条Conting片段,进一步组装得到平均长度为542 bp的151729条Unigenes。将闽楠转录组Unigenes进行基因功能注释,与NR数据库比对发现,其与葡萄的相似序列最多(34%),与黄瓜、野草莓、大豆的同源性较低(各占3%);进行GO功能注释,可将其划分为生物过程、细胞成分、分子功能3大类共计52个分支,与eggNOG数据库比对可分为25类,通过KEGG功能注释可知转录组中涉及的基因共参与了176条代谢通路,其中核糖体和碳代谢获得的注释较多。另外通过MISA软件分析,共获得35972个SSR位点。其中,单核苷酸、二核苷酸和三核苷酸为优势重复类型,SSR位点数分别为21762(60.50%),8931(24.83%),4924(13.69%)。闽楠转录组分析及基因功能注释为深入开展闽楠遗传育种及分子生物学相关研究奠定基础。     

基于高通量测序的慈竹笋转录组分析与基因功能注释

慈竹是我国四川当地的优势丛生竹种之一,其纤维长度和质量较优异,是造纸、纺织等工业的良好原料。本文利用Illumina Hi SeqTM 2000平台,对10、50、100和150 cm高的慈竹笋进行转录组分析,共得到69.28 M条读长(Reads),经从头拼接、组装和聚类后得到111 137条非重复序列基因Unigene,其中共有63 094条注释到COG、GO、KEGG、Swiss-Prot和Nr数据库中。这些Unigene不仅具有一般的功能,如转录和信号转导等,还涉及到蔗糖转运与代谢、次级代谢产物及细胞壁的生物合成等方面。不同高度慈竹笋的纤维素合成酶基因存在差异表达,发现了可能调控慈竹生长发育以及纤维素和木质素生物合成的相关基因,为慈竹品种改良提供一定的理论基础。     

野生矮牡丹高通量转录组测序分析

为探究野生稷山矮牡丹(Paeonia jishanensis)花青素等抗氧化活性物质生物合成的遗传基础,本文采用高通量测序技术Illumina NextSeqTM500对其幼嫩叶片、嫩茎和花芽混合池样本进行转录组测序分析。共获得39450个unigene,其中3563个unigene中有4106个SSR位点。在Swiss Prot数据库中注释到29420个unigene,其中11个与花青素生物合成过程相关,22个与类黄酮生物合成过程相关,15个与异黄酮生物合成过程相关;NR数据库中注释高达100%,有12个与类黄酮生物合成相关,7个与花青素生物合成相关;GO数据库中注释到24291个unigene,分为生物学过程、细胞组分及分子功能等三大类53个功能组,有21个unigene与黄酮、黄酮醇及异黄酮生物合成有关,4个unigene与花青素合成有关;eggNOG数据库中注释到38238个unigene,涉及植物生长发育过程绝大多数生命活动过程;KEGG数据库注释到5709个unigene,分别定位到6大类别42个代谢途径,其中25个unigene与花青素生物合成相关。研究表明,在稷山矮牡丹转录组序列中,共鉴定和发掘出多个涉及类黄酮化合物生物合成、花青素生物合成过程的相关基因序列及SSR位点,为下一步开展相关代谢合成途径研究奠定了基础。     

茶银尺蠖触角转录组及相关嗅觉基因的分析

本研究利用Illumina Hiseq 4000高通量测序平台对茶银尺蠖Scopula subpunctaria(Herrich-Schaeffer)雌雄成虫触角进行转录组测序,共获得93 335条unigenes,N50长度为1 151 bp。通过基因功能注释,有23 581条unigenes注释到六大公共数据库中,其中注释到NR数据库的unigenes最多,占39.96%,且与家蚕Bombyx mori相似基因数目最多,占20.32%;10 718条unigenes注释到GO数据库的生物过程、细胞组分和分子功能3个大类46个亚类;5 464条unigenes注释到KEGG数据库的38个小类代谢通路,其中注释在信号传导通路的unigenes最多,占10.64%。根据注释信息,共筛选得到60个候选嗅觉基因,包含17个气味结合蛋白基因（Odorant binding protein genes,OBPs）、13个化学感受蛋白基因（Chemosensory protein genes, CSPs）、14个嗅觉受体基因（Olfactory receptor genes,ORs）、13个离子型受体基因（Ionotropic receptor genes,IRs）和3个感觉神经膜蛋白基因（Sensory neuron membrane protein genes, SNMPs）。利用FPKM值对SsubOBPs和SsubCSPs基因进行表达量分析发现,SsubOBP2和SsubOBP3在茶银尺蠖雌雄成虫触角中均高丰度表达（FPKM>2000）, SsubCSP4、 SsubCSP7和SsubCSP9在茶银尺蠖雌雄成虫触角中均高丰度表达（FPKM>900）。本研究首次获得了茶银尺蠖触角转录组数据,并鉴定出候选嗅觉相关基因,结果为进一步研究茶银尺蠖的嗅觉感受机制提供分子基础。     

黄秋葵花和果荚转录组测序及类黄酮代谢差异表达分析

采用Illumina Hi seq 2500转录组高通量测序,构建黄秋葵花和果荚转录组文库,并利用测序评估、差异基因功能注释等生物信息学方法进行分析。研究结果表明：(1)分别获得黄秋葵花和果荚有效数据7.12 Gb和8.14 Gb,碱基百分比(Q30)均达到91.0%以上。(2)获得差异表达基因(DEGs)1 336个,其中上调基因319个,下调基因1 017个。(3)获得功能注释的基因有1 131个,GO将455注释基因归成41个功能小类,主要涉及代谢过程、催化活性、单一生物过程和细胞过程等过程;KOG注释了472个DEGs,涉及23个功能分类,其中与次生代谢直接相关过程O和Q类别获得111个注释结果;KEGG将372个DEGs注释到80条代谢通路中,获得F3H、F3′5′H、DFR、ANR、ANS、GT、LAR共10个关键差异基因,其中F3H、DFR在黄秋葵花中表现上调效应,F3′5′H、ANR、LAR在黄秋葵果荚中表现显著上调效应,ANS、GT则分别在花和果荚中均有上调或下调效应。(4)实时定量PCR分析显示,其中7个差异表达基因得到的相对表达量与转录组表达谱分析趋势完全一致。(5)类黄酮代谢途径分析表明,黄秋葵花通过F3H、DFR、ANS、GT途径将NAR催化为生成花青素苷;果荚则将NAR通过F3′5′H将催化为DHM,后在FLS催化下生成黄酮醇类物质等;部分NAR在F3H、 DFR催化下,生成无色飞燕草素苷元,再分别在ANS、LAR作用下,进入花青素苷元和原花青素合成途径。该研究结果丰富了黄秋葵转录组信息,为黄秋葵类黄酮物质纯化和功能利用提供参考依据。     

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

本研究采用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值对基因表达量进行评估.本研究获得的转录组信息及分析结果,为进一步研究阿尔泰蝠蛾的基因功能及低温生态适应性奠定分子基础.