不同海拔饲育的麦洼牦牛肺脏组织转录组学分析

为了探讨牦牛适应高海拔低氧环境的基因表达特征与规律,对在高海拔（3 560 m）和低海拔地区（478 m）饲育4个月的2.5~3岁健康雄性麦洼牦牛肺组织进行转录组测序。转录组测序采用Illumina高通量测序平台(HiSeqTM2500/4000)进行,并以qRT-PCR验证差异表达基因的表达量。结果显示,高海拔组牦牛肺脏转录组平均每个测序样本得到约5.76亿条Clean Reads,低海拔组牦牛中得到约6.10亿条Clean Reads,比对到参考基因组上的Reads数分别占91.74%和91.28%以上,共发现了2 047个新转录本。低海拔组与高海拔组牦牛肺脏组织之间共有199个差异表达基因,其中含89个差异上调表达基因和110个差异下调表达基因。所得差异表达基因富集在297个GO条目和146个KEGG通路中,包含62个低氧适应相关的GO条目和35个低氧适应相关代谢通路。其中低氧适应相关GO条目在生物过程、细胞组成和分子功能三种类别中占比最多的分别为细胞粘附、蛋白复合物和钙离子结合。低氧适应相关KEGG通路中占比最多的为肿瘤坏死因子(TNF)信号通路,其次为低氧诱导因子1(HIF-1)信号通路。qRT-PCR验证结果显示,Ⅱ类人类白细胞抗原α链(HLA-DOA、HLA-DRA)、补体因子 (C2)和甘露糖结合凝集素相关丝氨酸蛋白酶1（MASP1）基因的表达量变化与转录组测序结果相符。本研究为全局和深入理解牦牛肺组织转录本表达对高海拔低氧的响应提供了有价值的切入点。     

东方蜜蜂微孢子虫的基因结构优化及新基因鉴定

东方蜜蜂微孢子虫Nosema ceranae是一种寄生于蜜蜂中肠上皮细胞的单细胞真菌,对蜜蜂的健康危害严重,给世界各国的养蜂业造成较大损失。本研究基于前期获得的N.ceranae孢子的转录组数据对其已注释基因进行结构优化,并对未注释基因进行预测和分析。通过将测序得到的clean reads比对参考基因组和转录本重构,共对10个N.ceranae的已注释基因的5'端或3'端进行了延长。利用Cuffcompare软件将重构转录本与参考基因组进行比对,共鉴定出27个新基因,随机挑选9个新基因进行RT-PCR验证,均能扩增出符合预期的目的片段,表明预测出的新基因真实存在。有6个新基因能够注释到GO数据库和6个基因注释到KEGG数据库。进一步分析结果显示上述新基因注释到细胞等10个GO条目上,它们可能在N.ceranae的生命活动中具有重要功能。研究结果为N.ceranae的基因结构和功能注释信息的完善提供了有益补充,也为新基因的功能研究打下了基础。     

牦牛体外成熟卵母细胞差异转录组学研究

卵母细胞体外成熟环节是现代繁殖育种技术的基础环节,体外成熟的MⅡ期卵母细胞的质量直接影响了卵母细胞的后续受精及受精卵卵裂.由于卵母细胞的发育成熟是一个涉及大量基因转录调控的复杂过程.因此,转录组学研究是了解卵母细胞发生发育的关键.本研究以高海拔地区的特有牛属动物牦牛(Bos grunniens)卵母细胞作为研究对象,应用RNA-seq技术对牦牛GV期卵母细胞及体外成熟MⅡ期卵母进行转录组学测序及比对分析.经深度测序后,各获得一个包含51380686条过滤测序序列(cleanreads),4624261740碱基(bp)的GV期牦牛卵母细胞测序文库和一个包含50303412条过滤测序序列(clean reads),4527307080碱基(bp)的MⅡ期牦牛卵母细胞测序文库.基因覆盖率统计表明,GV期和MⅡ期文库中分别有16719条和16339条牦牛基因得到转录.通过比较分析MⅡ和GV期转录组数据,共筛选出4767个差异表达基因,其中1418个基因表达量上调和3349个基因表达量下调.GO功能分类注释显示,差异基因与"代谢过程"等生物学过程存在密切关联,同时结合分子功能相关类别等分子事件表现活跃.按表达量上调和下调分别对差异基因进行KEGG通路互作分析,结果表明,凋亡通路、内吞通路及代谢通路是上调的关键调控通路,而代谢通路、丙酮酸代谢、黏着斑及氧化磷酸化是下调的关键调控通路,其中代谢通路在上调及下调通路互作中均起关键调控作用.进一步对代谢通路进行分析,结果表明,组氨酸代谢、甘氨酸、丝氨酸和苏氨酸代谢、乙醛酸和二羧酸代谢及谷胱甘肽代谢是整个代谢通路的关键调控节点。此外,本研究发现,在牦牛卵母细胞体外成熟过程中,与免疫相关的白细胞介素家族及干扰素家族基因发生了显著差异表达.该研究结果为进一步解析牦牛卵母细胞发育的分子机制及全面了解牦牛繁殖的特异性提供了基础.     

基于高通量转录组测序的草鱼雌雄性腺差异表达基因分析

草鱼雌性个体生长优势明显高于雄性,挖掘精巢和卵巢差异表达的功能基因对草鱼生殖调控具有重要的价值及应用前景。采用Illumina HiSep 2500转录组测序技术分别对12月龄的3尾雌性草鱼和3尾雄性草鱼的性腺组织的mRNA进行测序分析,精巢和卵巢共得到8 363个差异表达基因。与卵巢相比,精巢上调基因6 446个(77%),下调基因1 917个(23%),差异明显。将上述差异基因在Nr、GO和KEGG数据库进行比对,获得48种Go功能注释分类和313条KEGG代谢通路。分析发现,Go功能分类中"绑定"功能富集差异基因的数量最多,占总数的15.8%;KEGG代谢通路中"信号转导"通路富集差异基因的数量最多,占总数的10.2%。结合转录本数据,挑选与草鱼性腺发育相关的关键差异基因Dmrt1、amh、foxl2、cyp19a1a进行分析,研究显示,与卵巢相比,Dmrt1和amh基因在草鱼精巢中极显著高表达(p0.01);与精巢相比,cyp19a1a和foxl2基因在草鱼卵巢中显著高表达(p0.05),与转录组数据分析结论一致,说明Dmrt1和amh基因与早期草鱼精巢发育调控相关,cyp19a1a和foxl2基因与早期草鱼卵巢发育调控相关。本实验为进一步了解草鱼性别基因的调控机理提供数据依据。     

颜氏大疣蛛毒腺转录组学研究

颜氏大疣蛛Macrothele yani个体大、毒性强、易取毒，是农林生态系统中害虫的重要天敌，也是天然药物研究工作者进行新药挖掘和毒液开发利用的重要选择对象。本研究利用高通量测序平台Illumina Hiseq 2000对颜氏大疣蛛雌雄成体的毒腺进行转录组测序和生物信息学分析，共获得转录组样本数据37.8 G，总计301 024条unigenes，总长度170 512 372 bp，最短201 bp，最长36 105 bp，平均长度566 bp，GC含量平均值为38.17%，N50长度为705 bp。将unigenes序列与NR（非冗余蛋白序列数据库）、String（蛋白质相互作用数据库）、Swiss-prot（蛋白质序列数据库）、Pfam（蛋白质家族数据库）、GO（基因本体论）、KOG（真核生物蛋白质直系同源数据库）、KEGG（京都基因与基因组百科全书）数据库进行比对（e≤10^-10），分别获得31 409、6 549、11 817、8 216、11 480、7 804条unigenes注释。通过生物信息学对雌雄颜氏大疣蛛unigenes进行表达量分析，筛选高表达量且高可信度的差异表达基因，并进行对这些差异基因进行GO和KEGG功能富集分析；与雄蛛相比，雌成体有205条unigenes表达量上调，113条unigenes表达量下调。最后在所有的转录本数据中共比对到68条毒素相关序列。本研究获得的颜氏大疣蛛转录组信息，为颜氏大疣蛛的功能基因挖掘提供了重要的信息资源。     

荔枝蒂蛀虫转录组及嗅觉相关基因分析

【目的】荔枝蒂蛀虫Conopomorpha sinensis Bradley是专一性危害我国华南地区荔枝和龙眼的重要害虫,隐蔽性强,防治困难,基因组信息缺乏。本研究的目的是获得荔枝蒂蛀虫的基因数据,寻求有效控制害虫的分子靶标。【方法】采用新一代高通量测序技术Illumina Hi SeqTM4000对荔枝蒂蛀虫进行转录组测序和生物信息学分析。【结果】经序列拼接获得68 996条unigenes。进一步利用七大公共数据库进行同源比对,注释了22 348 unigenes。注释到Nr数据库的unigenes数量最多,达27.01%,其中Nr注释的荔枝蒂蛀虫unigenes中与小菜蛾Plutella xylostella unigenes同源性最高,达34.1%。将unigenes与GO数据库比对发现,15 585条unigenes根据其功能大致可分为3类47亚类。KEGG pathways分析表明,7 272条unigenes定位为267个代谢通路。基因注释进一步筛选鉴定获得100个荔枝蒂蛀虫嗅觉相关基因;与鳞翅目相关气味结合蛋白基因联合分析发现,与荔枝蒂蛀虫气味结合蛋白基因直系同源的有18组,部分基因形成独立一簇。【结论】本研究首次获得了荔枝蒂蛀虫的转录组数据,研究结果为生物控制荔枝蒂蛀虫提供了重要的基础数据和候选分子靶标。荔枝蒂蛀虫独有的气味结合蛋白基因可能与其生境中特有的化学物质相关。     

虹鳟肝组织新转录本分析及基因结构优化

目的虹鳟热应激下肝RNA-seq数据中新转录本的分析及已注释基因结构优化。方法以虹鳟肝为材料提取总RNA,构建cDNA文库,并利用Illumina双端测序Hiseq 2500平台进行测序。运用Cufflinks软件对测序数据进行组装,将其与虹鳟参考基因组进行序列比对。结果发掘新转录本6555个,其中30个新转录本在热应激前后差异表达(P<0.05)。与GO数据库比对对新转录本进行功能注释,获得3097个新转录本的注释。与KEGG数据库比对,共有3617个新转录本注释到284条代谢通路中。对19 424个已注释基因的结构进行优化,延伸了14 719个基因的5′端和14 796个基因的3′端。结论通过对发掘的6555个新转录本分析,并对19 424个已注释基因结构优化,为虹鳟基因组注释信息的完善提供了有力的借鉴,并为进一步了解虹鳟热应激的机制提供更有力的理论基础。     

RNA-Seq高通量测序技术在果树功能基因组学研究的应用进展

RNA-Seq又称为"转录组测序技术",它能够从整体水平研究物种基因功能以及基因结构,揭示特定生物学过程。应用该项技术能有效地解决因果树基因组庞大等问题,有效地揭示和阐明植物体内次生代谢生物合成途径及代谢机制。主要论述了RNA-Seq技术的应用方法与流程,通过对测序数据的处理分析与比对等方法,分析了RNA-Seq技术在果树新基因挖掘、次生代谢产物合成途径的应用进展,基于RNA-Seq的高通量测序技术对于发现和揭示具有重要生物学功能和经济价值的次生代谢产物合成关键基因及其调控机制提供了可能。     

发菜RNA-seq转录组分析及耐旱相关基因筛选

发菜是一种陆生耐旱固氮蓝藻,具有重要的生态和资源价值,但目前发菜的基因组和转录组信息还比较缺乏。本研究采用Illumina Hi Seq TM 2500高通量测序平台对充分吸水和失水后的发菜藻体进行转录组分析,经过拼接组装共获得30 013 560条clean reads,3 100个Unigene,序列平均长度854 nt。将Unigene序列与NR、NT、Swiss-Port、PFAM、GO、COG、KEGG(E-value1.0E~(-5))等数据库进行比对,共有2 874个Unigene获得了注释。通过GO功能分类,917个Unigene映射到GO不同功能节点上;通过Unigene与COG数据库的比对,可将有COG功能的514个Unigene分为20类;通过KEGG pathway分析,有1 220个编码基因参与了135条已知的代谢通路。提示谷胱甘肽、类胡萝卜素、淀粉和蔗糖、ABC转运体等代谢途径的部分基因在失水藻体上调表达,说明这些通路与发菜耐旱性关系密切。因此这些注释信息为今后发菜耐旱基因的挖掘提供了重要依据,同时也为发菜分子生物学的进一步研究提供了试验基础。     

携带与未携带松材线虫的松墨天牛转录组差异表达基因分析

【目的】本研究旨在建立携带与未携带松材线虫Bursaphelenchus xylophilus的松墨天牛Monochamus alternatus成虫各组织的转录组数据库,揭示松墨天牛对松材线虫响应的转录组整体表达特征。【方法】以携带和未携带松材线虫的松墨天牛成虫的表皮、气管和脂肪体为材料,采用Illumina HiSeq TM 2000测序平台开展转录组测序,利用Trinity软件对RNA-Seq数据进行从头组装,利用NCBI数据库进行基因注释。通过DEG-Seq分析携带和未携带松材线虫的松墨天牛成虫中的差异表达基因,对上调表达基因进行GO和KEGG代谢途径富集分析,并通过实时荧光定量PCR技术对部分差异表达的基因包括ENV,CDK1,HSP70-C,HSP70,HSP75,DUO,PABPN1和IGFP基因的表达水平进行验证。【结果】测序过滤后共获得55059条单基因簇(unigene),平均长度为1536 bp。将unigene与数据库中的序列进行BLASTX比对,成功注释24354条unigenes,其中4022条unigene在GO数据库中获得注释,6098条unigene在KEGG数据库中获得注释。经过GO与KEGG富集分析显示,松墨天牛成虫携带松材线虫后,差异表达基因主要为与其气管中压力调节、组织与DNA修复、激素响应方面相关的基因。实时荧光定量PCR分析结果显示CDK1,HSP70,HSP75,DUO,PABPN1和IGFP基因的表达量在携带松材线虫的松墨天牛气管中相比未携带松材线虫的松墨天牛气管中的显著上调。【结论】本研究初步阐明松墨天牛携带松材线虫后转录组的整体表达模式,为进一步研究松墨天牛的基因功能及对松材线虫的胁迫响应过程奠定了基础。     

大蜡螟触角转录组及嗅觉相关基因分析

【目的】大蜡螟Galleria mellonella Linnaeus是世界范围内蜂群中普遍存在的害虫,给世界养蜂业带来了巨大的危害。本研究旨在建立大蜡螟触角转录组数据库,挖掘大蜡螟嗅觉相关基因,为今后利用嗅觉基因靶标防治大蜡螟提供分子生物学信息数据。【方法】利用高通量测序平台(Illumina HiSeqTM 2500)对大蜡螟触角进行转录组测序和组装,通过筛选转录组数据库,鉴定出嗅觉相关基因。【结果】成功构建了大蜡螟触角转录组,经序列拼接获得188 278条unigenes,平均长度为781 bp,N50为1 161 bp。使用BLAST软件将大蜡螟触角unigenes序列与7大公共数据库比对,共注释到108 047条unigenes,其中NR数据库注释的unigenes最多,为78 746条,占41.82%,且NR注释的大蜡螟unigenes中与家蚕Bombyx mori类似基因最多,为18.4%。将unigenes与GO数据库比对发现,共有69 794条unigenes注释到GO数据库中,按照其所参与的生物过程、细胞组分和分子功能对其进行GO功能分类,共含有55个亚类,其中参与结合和催化功能及代谢过程的unigenes比例较大;KEGG代谢途径分析表明,14699条unigenes形成231条代谢通路,其中被注释在信号传导通路中的unigenes最多,为2 401条,占16.33%。进一步基因注释分析,鉴定得到226个候选嗅觉相关基因,包括52个气味结合蛋白(Odorant binding proteins,OBPs)基因,36个化学感受蛋白(Chemosensory proteins,CSPs)基因,80个气味受体(Odorant receptors,ORs)基因,56个离子型受体(Ionotropic receptors,IRs)基因和2个感觉神经元膜蛋白(Sensory neuron membrane proteins,SNMPs)基因;通过比较雌、雄大蜡螟转录组鉴定出114个差异表达基因,包括5个嗅觉相关基因(OBP8、OBP17、CSP3、CSP34和OR15);114个差异表达基因中,66个基因在雌蛾触角中高表达,48个基因在雄蛾触角中高表达。【结论】本研究获得了大蜡螟触角转录组数据,并鉴定出候选嗅觉相关基因,结果为进一步研究大蜡螟的基因功能及嗅觉感受机制奠定分子基础。     

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

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

绿眼赛茧蜂转录组分析及嗅觉相关蛋白基因的鉴定

[目的]绿眼赛茧蜂Zele chlorophthalmus是草地螟Loxoste gesticticalis幼虫重要天敌之一.本研究通过构建绿眼赛茧蜂触角转录组数据库,挖掘其与嗅觉相关的蛋白基因,为更好的利用绿眼赛茧蜂防治草地螟发挥其生防潜能提供理论依据.[方法]以Illumina Novaseq 6000高通量测序平台为基础,将绿眼赛茧蜂触角的基因进行转录组测序、组装序列,以及完成其生物信息学的研究分析,并对绿眼茧赛蜂触角的相关嗅觉基因做鉴定分析.[结果]成功构建绿眼赛茧蜂转录组数据库,数据库中的unigenes序列为65228条,N50为3882 bp.使用BLAST软件将绿眼赛茧蜂触角unigenes序列各自和Pfam、Swiss-Prot、NR、COG、KEGG、GO权威数据库进行对比,并完成基因功能相关注释,共注释基因数为18662条,占总数的28.61％.其中,NR数据库获得的注释最多,占总数的24.61％,为15863条,KEGG数据库获得的注释最少,为9612条(14.91％),其他依次为Pfam数据库注释数据库12164条(18.86％)、COG数据库注释15584条(24.17％)、GO数据库注释11634条(18.05％),Swiss-Prot数据库注释达到11634条,为总数的18.86％.借助GO数据库对unigenes的注释,其功能供分为三大类,可以细分49个分支,主要包括分子功能、细胞组分以及生物学过程.通过注释基因功能对嗅觉相关基因进行筛选,共发现151条和嗅觉有关的蛋白基因,包括3个感觉神经元膜蛋白基因、22个离子型受体基因、23个味觉受体基因、83个气味受体基因、6个化学感受蛋白基因、14个气味结合蛋白基因.[结论]成功收集了绿眼赛茧蜂触角转录组相关数据,并对与嗅觉相关的蛋白进行鉴定分析,为深入研究基因功能及嗅觉分子机制提供理论依据.     

意大利蝗转录组及性腺发育相关基因分析

【目的】意大利蝗Calliptamus italicus是新疆荒漠、半荒漠草原优势危害种类之一,其生殖相关基因信息缺乏。本研究对意大利蝗精巢和卵巢组织进行高通量转录组测序,旨在揭示意大利蝗转录组特征,并获得与性腺发育相关的基因数据。【方法】利用Illumina高通量测序平台(Hi Seq/Mi Seq)对意大利蝗进行了转录组测序,并进行序列分析。【结果】获得718 872条unigenes,平均长度631 bp,N50为1 186 bp。通过同源性比对,成功注释254 597条unigenes,其中,Nr数据库注释占28.87%,比例最高。Nr数据库注释的unigenes与黑褐草蚁Lasius niger的相似基因序列最多,为10.80%。与GO数据库比对发现,意大利蝗转录组unigene可分为3大类:涉及分子功能的unigenes有31 392条,涉及细胞组分的unigenes有20 586条,与生物学过程有关的unigenes有57 014条。KEGG pathways分析表明,23 666条unigenes归属于251个通路。涉及生殖系统发育的通路有卵母细胞减数分裂、胰岛素信号通路、Wnt信号通路、促性腺激素释放激素信号通路、转化生长因子β信号通路以及泛素介导的蛋白水解、黄体酮介导的卵母细胞成熟、昆虫激素生物合成等。进一步筛选得到部分与性腺发育相关的基因,包括vg,vgr,sox 21b,testis-specific serine/threonine-protein kinase,spermatogenesis-associated protein和vasa-like gene等。【结论】本研究获得了意大利蝗转录组数据及性腺发育相关基因,为进一步研究意大利蝗生殖调控机理奠定了分子基础。     

花椒窄吉丁转录组及化学感受相关基因的分析

【目的】建立花椒窄吉丁 Agrilus zanthoxylumi 转录组数据库,挖掘其基因组数据。【方法】采用高通量测序平台Illumina NovaSeq6000对花椒窄吉丁成虫不同组织(触角、头、胸、腹、足、翅)进行转录组测序、序列组装,使用BLAST软件将花椒窄吉丁unigenes序列与公共数据库进行比对,BLAST同源性搜索的方法从中筛选出与其化学感受相关的基因,并与其他已研究发表的鞘翅目昆虫化学感受相关基因的核酸序列比对,进行系统发育分析。利用RPKM值对这些基因在花椒窄吉丁成虫不同组织中的表达量进行分析。【结果】经测序及序列拼接后共得到80 320条unigenes和 169 398 条contigs,其G+C比例分别是37.63%和39.18%,平均长度分别为828.75和1 084.33 bp;unigenes长度主要分布在200~400 bp的共有37 374条,占全部unigenes的46.53%。在NR数据库中注释的与花椒窄吉丁转录组相似基因所属物种分布为赤拟谷盗 Tribolium castaneum 所占比例26.80%,其后依次是小家鼠 Mus musculus (13.87%),蛀犀金龟 Oryctes borbonicus (5.26%),山松甲虫 Dendroctonus ponderosae (4.10%),黑蚁 Lasius niger (2.36%),其他物种所占比例是47.61%。在GO数据库中比对到27 488个unigenes,可分为分子功能、细胞组分和生物学进程三大类共59个过程;利用KOG数据库功能注释分为25类,注释到unigenes 最多功能类别的是普通功能,共 4 666个,最少的是与细胞移动相关,共40个;将80 320个unigenes映射到KEGG数据库中,21 104条unigens注释到代谢通路35个,占26.27%,注释最多的代谢途径是转录途径,共涉及到2 037条unigenes。从NR数据库中共注释到8个气味受体(odorant receptor, OR)基因、7个离子型受体(ionotropic receptor, IR)基因和6个味觉受体(gustatory receptor, GR)基因,各组织中的基因表达量分析发现,这3类化学感受相关基因在成虫触角中的表达量显著高于其他组织中的表达量,雌、雄触角中表达量差异最显著的为 AzanOR 1,而在雌、雄成虫足中这3类化学感受相关基因表达量没有差异。【结论】本研究首次获得了花椒窄吉丁转录组数据,为进一步研究花椒窄吉丁的基因功能奠定了分子基础。     

昆虫RNA-Seq数据的分析流程

随着高通量RNA测序(RNA-Seq)技术的发展和测序成本迅速下降,RNA-Seq技术已经成为生物学研究的重要工具,为生物学家全面地了解和研究转录组提供了机遇。高通量测序具有读长短、存在一定比例的测序错误、数据量大等特点,因此RNA-Seq数据分析与基因组分析和传统的EST数据分析有所不同。本文通过介绍不同的测序平台、原始数据产生和低质量数据过滤的计算流程,对短序列比对、转录组拼接、功能注释、以及差异表达分析进行了研究和分析,最后对RNA-Seq在昆虫学研究中的应用进行了综述,并对RNA-Seq技术进行了总结和展望。     

基于转录组测序的羽衣甘蓝叶色相关基因分析

叶色是羽衣甘蓝重要的观赏性状之一。本研究采用Illumina Hi-Seq2500高通量测序技术,对基因型纯合的紫叶和白叶羽衣甘蓝叶片进行转录组测序,筛选差异基因并与GO和KEGG数据库比对进行注释分析,分析羽衣甘蓝叶色形成相关基因。结果显示,获得高质量短读序共104 608 770条,筛选出紫叶相对白叶的差异表达基因1 993个,其中上调表达基因1 094个,下调表达基因899个。根据GO功能分类可分为生物过程、细胞组分和分子功能3大类64功能组。根据KEGG代谢通路分析可以分为171类,在叶色相关的类黄酮生物合成途径中黄酮醇合成酶(FLS)、二氢黄酮醇4-还原酶(DFR)上调以及类胡萝卜素生物合成途径中的类胡萝卜素β-环化酶上调与紫叶形成关系密切。本研究丰富了羽衣甘蓝的转录组信息,获得了一些差异表达基因,为进一步研究羽衣甘蓝叶色形成的遗传机制提供了有价值的信息。     

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

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

绿色杜氏藻转录组分析

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