首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到17条相似文献,搜索用时 171 毫秒
1.
蜜蜂球囊菌基因结构优化及新基因鉴定   总被引:1,自引:1,他引:0  
为补充和完善蜜蜂球囊菌Ascosphaera apis基因组序列和功能注释信息,本研究利用RNA-seq技术对球囊菌菌丝和孢子进行深度测序,并基于高质量的转录组数据对已注释基因进行结构优化,对未注释基因进行预测、鉴定和分析。将测序得到的有效读段进行参考基因组比对和转录本重构;利用Cuffcompare软件将重构转录本与参考基因组进行比对。结果表明:共对101个已注释基因的5′端或3′端进行了延长;共鉴定出373个新基因,随机挑选10个新基因进行RT-PCR验证,其中8个能扩增出符合预期的目的片段,表明预测出的多数新基因真实存在;共计147个球囊菌新基因可注释到Nr和eggNOG数据库中,85个新基因注释到29个GO条目,66个新基因注释到33条代谢通路。本研究结果补充了球囊菌的基因结构和功能注释信息,也为新基因的功能研究打下了初步基础。  相似文献   

2.
西方蜜蜂的基因注释信息优化及新基因鉴定   总被引:1,自引:1,他引:0  
西方蜜蜂(Apis mellifera)全基因组序列早在2006年就已公布,但其基因注释信息并不完善。本研究利用前期获得的意蜂幼虫肠道RNA-seq数据对西方蜜蜂基因组的已知基因进行结构优化,对新基因进行预测、鉴定和分析。利用Cufflink软件进行转录本重构,进而与参考转录本序列进行比较,共对4 177个已公布西方蜜蜂基因的测序结果进行了优化,其中5′端延长、3′端延长、5′和3′端都延长的基因分别有2 383、2 392和2 384个。利用Cuffcompare软件将重构转录本与参考基因组进行比对,共预测出197个西方蜜蜂新基因,随机挑选20个新基因进行RT-PCR验证,有18对引物成功扩增出符合预期的目的片段。GO分类结果显示这些新基因涉及运输和代谢、能量代谢和信号分子和交互等21个GO terms。KEGG代谢通路(pathway)富集分析结果显示这些新基因富集于13个pathways,其中基因富集数最多的是细胞进程、单一有机体和结合。本研究为西方蜜蜂基因组信息提供了有益的补充和完善,也为新基因的功能研究打下了基础。  相似文献   

3.
【目的】利用已获得的纳米孔全长转录组数据对现有的东方蜜蜂微孢子虫(Nosema ceranae)参考基因组的基因序列和功能注释进行完善。【方法】采用TransDecoder软件预测东方蜜蜂微孢子虫基因的开放阅读框(open reading frame,ORF)及相应的氨基酸。利用gffcompare软件将全长转录本与参考基因组注释的转录本进行比较,对基因组注释基因的非编码区向上游或下游延伸,修正基因的边界。利用MISA软件鉴定长度在500 bp以上的全长转录本的简单重复序列(simple sequence repeat,SSR)位点,包括单核苷酸重复、双核苷酸重复、三核苷酸重复、四核苷酸重复、五核苷酸重复、六核苷酸重复、混合SSR等类型。通过Blast工具将鉴定到的新基因和新转录本比对Nr、KOG、eggNOG、GO和KEGG数据库,从而获得功能注释。【结果】共预测出2 353个完整ORF,其中长度分布在0—100个氨基酸的ORF最多,占总ORF数的72.12%。共对东方蜜蜂微孢子虫的2 340个基因进行了结构优化,其中5′端延长的基因有1 182个,3′端延长的基因有1 158个。共鉴定到1 658个SSR,其中单核苷酸重复、双核苷酸重复、三核苷酸重复、四核苷酸重复的数量分别为1 622、23、7和6个;单核苷酸重复类型的SSR密度最大,达到182.32个/Mb,其次为混合SSR、双核苷酸重复和三核苷酸重复,分别达到6.90、2.78和0.73个/Mb。共鉴定出954个新基因,其中分别有951、333、371、422和321个新基因可注释到Nr、KOG、eggNOG、GO和KEGG数据库。此外,还鉴定出6 164条新转录本,其中分别有6 141、2 808、2 932、3 196和2 585条新转录本可注释到Nr、KOG、eggNOG、GO和KEGG数据库。新基因和新转录本注释数量最多的物种均为东方蜜蜂微孢子虫,其次是蜜蜂微孢子虫(Nosema apis)。【结论】研究结果较好地完善了现有的东方蜜蜂微孢子虫参考基因组已注释基因的序列和功能注释,并补充和注释了大量参考基因组未注释的新基因和新转录本。  相似文献   

4.
采用Blast工具将已鉴定到的意大利蜜蜂(Apis mellifera ligustica,简称意蜂)全长转录本比对Nr数据库,共鉴定到意蜂Hippo信号通路中相关的49个基因及其550条全长转录本。运用Gffcompare软件将鉴定到的Hippo信号通路相关全长转录本与西方蜜蜂参考基因组(Amel_HAv3.1)上注释的转录本进行比较,分别延长西方蜜蜂参考基因组注释到Hippo信号通路的9个基因的5'' UTR和7个基因的3'' UTR。通过Astalavista软件鉴定到Hippo信号通路相关的4个基因的7次可变剪切(AS)事件,包括2次外显子跳跃、2次内含子保留和3次可变5''端剪接。通过PCR验证了1个基因的AS事件真实性。利用TAPIS pipeline鉴定到意蜂Hippo信号通路相关的24个基因含有1个及以上可变多聚腺苷酸化(APA)位点,其中含有1个APA位点的基因数量最多。通过TBtool软件鉴定APA位点上游motif的一致性序列HVNCCNBNDYVDVVHVNDBVDYCNBBDNNNHNNVNNVMNN。  相似文献   

5.
根据中华蜜蜂(Apis cerana cerana)基因组序列,鉴定抗氧化酶系统的主要组分,并用RNA-Seq数据分析了工蜂行为发育中抗氧化基因在脑组织中的转录水平.研究表明,共有47个抗氧化基因,包括34个抗氧化酶基因和13个硫氧还蛋白基因.比较东方蜜蜂(Apis cerana)、西方蜜蜂(Apis mellifera)和黑腹果蝇(Drosophila melanogaster)的抗氧化基因,抗氧化酶系统的主要组分在进化上保守,但抗氧化酶基因家族的成员存在物种间差异.如东方蜜蜂和黑腹果蝇的SOD基因家族含3个成员,而西方蜜蜂有2个;果蝇有43个GST基因,而东方蜜蜂和西方蜜蜂仅有14和13个.RNA-seq分析显示,抗氧化基因的表达水平随着工蜂行为发育而发生规律性变化.工蜂脑表达的46个抗氧化基因,约有60%在由内勤蜂(哺育蜂和守卫蜂)转变为外勤蜂(采集蜂和侦察蜂)时明显表达上调,说明外勤蜂脑组织清除活性氧的酶促反应速率更高.  相似文献   

6.
【目的】利用RNA seq技术对中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫肠道参考转录组进行de novo组装,并进行功能及代谢通路注释,进而利用该转录组数据进行中蜂幼虫的SSR分子标记鉴定。【方法】实验室条件下饲养中蜂幼虫,将纯化的蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)孢子饲喂3日龄幼虫,剖取4、5和6日龄幼虫肠道,液氮速冻。将健康幼虫肠道与感染球囊菌的幼虫肠道同时进行Illumina测序。通过对raw reads的过滤得到clean reads,利用Trinity软件组装得到unigenes。通过BLASTx(E-value10-5)比对NCBI Nr、Swiss-Prot、KOG和KEGG数据库,对unigenes进行功能和代谢通路注释。利用MISA软件对所有unigenes进行SSR搜索,并利用Primer Premier 5软件设计特异性SSR引物,通过常规PCR对来源于北京、辽宁兴城和四川成都的中蜂幼虫肠道样本进行SSR位点鉴定。【结果】中蜂幼虫肠道的RNA seq共得到35 670 000条reads,de novo组装得到43 557个unigenes,平均长度为898 nt。共有18 225个unigenes可被注释到上述公共蛋白数据库,单独注释到NCBI Nr、Swiss-Prot、KOG和KEGG数据库的unigenes数分别为3 899、443、37和10个。KOG注释结果显示,11 442条unigenes分布于25个基因家族,其中注释到RNA加工和修饰家族的基因数最多,达1 249个。9 679个unigenes的GO分类结果显示,在生物学进程分类中,注释到细胞进程的基因最多,达4 201个,在分子功能和细胞组分类中,注释到结合与细胞的基因数最多,分别为4 935和2 900个。4 517个unigenes可注释到KEGG数据库中的216个代谢通路,注释到核糖体的基因数最多,达385个。利用MISA软件,可在7 763个unigenes搜索到13 448个SSR位点,随机选取20对SSR引物对国内3个不同来源的中蜂幼虫肠道样本的SSR位点进行扩增,有6对引物可鉴定出SSR分子标记。【结论】成功组装并注释了中蜂幼虫肠道参考转录组,可为中蜂及其幼虫的分子生物学及组学研究提供重要的参考信息,也可用于补充、丰富和检验东方蜜蜂的参考基因组,基于此转录组数据开发出6个中蜂的SSR分子标记,可应用于中蜂的基因图谱构建、基因多样性分析、基因定位等研究,也说明利用转录组数据开发非模式生物SSRs的方法可行。  相似文献   

7.
杨冠煌等(1986)发表蜜蜂阿坝亚种Apis cerana abansis,其亚种学名并没有拉丁化,应为Apis cerana abaensis,因多种原因未得到认可,但形态测定和分子数据分析认可中国阿坝地理种群Aba cerana。笔者介绍蜜蜂种下分类时,测评了38个形态学指标的测量方法、阿坝地理种群研究历史和背景、青藏高原分化出阿坝地理种群的环境条件,以及该种群个体形态特征、生物学和分布;中国是世界东方蜜蜂多样性热点地区之一,提出了中国东方蜜蜂(工蜂)中华亚种Apis cerana cerana、印度亚种Apis cerana indica、海南地理种群Hainan cerana及阿坝地理种群Aba cerana的形态鉴定表,肯定中国东方蜜蜂的多样性;从中国蜂业生产实际出发,首次较好地解决了东方蜜蜂种和种下阶元(亚种和地理种群)的关系,提出保护东方蜜蜂阿坝地理种群的措施。  相似文献   

8.
杨冠煌等(1986)发表蜜蜂阿坝亚种Apis cerana abansis,其亚种学名并没有拉丁化,应为Apis cerana abaensis,因多种原因未得到认可,但形态测定和分子数据分析认可中国阿坝地理种群Aba cerana。笔者介绍蜜蜂种下分类时,测评了38个形态学指标的测量方法、阿坝地理种群研究历史和背景、青藏高原分化出阿坝地理种群的环境条件,以及该种群个体形态特征、生物学和分布;中国是世界东方蜜蜂多样性热点地区之一,提出了中国东方蜜蜂(工蜂)中华亚种Apis cerana cerana、印度亚种Apis cerana indica、海南地理种群Hainan cerana及阿坝地理种群Aba cerana的形态鉴定表,肯定中国东方蜜蜂的多样性;从中国蜂业生产实际出发,首次较好地解决了东方蜜蜂种和种下阶元(亚种和地理种群)的关系,提出保护东方蜜蜂阿坝地理种群的措施。  相似文献   

9.
对前期获得的家蚕中肠组织转录组数据进行进一步生物信息学分析,利用Cufflinks软件对Mapped Reads进行组装,并与参考基因组注释信息进行比对,共发掘新基因788个。利用Blast软件将发掘的新基因与NR,Swiss-Prot数据库比对,其中746个新基因得到功能注释。这些新发掘基因在家蚕28个连锁群上都有分布,其中在第15连锁群上最多。有198个新基因注释到KEGG数据库中,分布于85条已知的通路中,将新基因与COG数据进行比对,并进行功能注释与分类,总共有258个新基因得到了注释,被分为22个COG类别,部分新基因的克隆分析,与预期结果一致,这些发现为以后进一步研究新基因的功能提供了基础信息。  相似文献   

10.
东方蜜蜂Apis Cerana Fab.在我国的分布及其亚种分化   总被引:5,自引:1,他引:4  
 本文论述了东方蜜蜂在中国的地理分布范围:除新疆和内蒙古自治区北部以外的广大地区;确定了东方蜜蜂在我国有五个亚种分化;除指名亚种Apis cerana cerana Fabricu外,肯定了西藏亚种Apis ceranaskorikovi Maa,发现了我国蜜蜂亚种的新纪录印度亚种Apis cerana indica Fabricus,提出了阿坝亚种Apis cernaabansis Yang et Kuang.subsp.now.和海南亚种Apis cerana hainana Yanget Krang sbusp.now.两个新亚种。  相似文献   

11.
旨在分析一株牦牛源植物乳杆菌SWUN5815的全基因组序列测序,并对该菌株的功能特性基因进行挖掘。基于PacBio RSII测序平台对乳杆菌SWUN5815进行全基因组测序分析和GO、KEGG、COG和NR数据库进行基因组基本功能注释。结果表明:1)植物乳杆菌SWUN5815的基因组大小为3.27 Mb,GC含量44.59%;2)预测到3 258个编码基因,编码基因总长度为2 814 147bp,平均长度为864bp;3)编码基因通过功能数据库对该菌株基因组基本功能注释和代谢通路基因信息注释得知,SWUN5815基因组中包括5个基因参与抗氧化活性过程,2个基因参与免疫过程;基因组中包含了抗生素、生物降解等代谢过程;可调控免疫和炎症的相关通路;基因组中含有18种已知细菌素基因;有4个合成ABC型细菌素转运系统的功能序列。综上,SWUN5815全基因组测序及分析挖掘其特异性功能基因有利于研究其益生特性提供基础,为菌株作为抗生素替代品的益生菌和饲料添加剂的应用研究提供了重要的参考数据。  相似文献   

12.
为探究中华蜜蜂(Apis cerana cerana,以下简称中蜂)与意大利蜜蜂(Apis mellifera ligustica,以下简称意蜂)幼虫的球囊菌(Ascospahaera apis)抗性差异产生的原因,利用Venn分析、GO分类分析、KEGG代谢通路分析以及Ka/Ks分析对二者的转录组进行比较研究。Venn分析结果显示,中蜂与意蜂的同源基因有17 656个,归属于8 111个基因家族,二者的特有基因分别有1 158和241个,分别归属于468和86个基因家族。GO分类结果显示,中蜂和意蜂的特有基因分别富集在38和28个GO条目。KEGG代谢通路富集分析结果显示,中蜂和意蜂的特有基因分别富集于96和21个代谢通路;进一步分析发现中蜂的特有基因富集在内吞作用、溶酶体、泛素介导的蛋白水解和黑化作用等4个细胞免疫通路,以及MAPK信号通路和Toll-like受体信号通路等2个体液免疫通路,而意蜂仅有1个特有基因富集在MAPK信号通路。Ka/Ks分析结果显示受到强烈正向选择、弱正向选择、中性选择和负向选择的单拷贝同源基因分别有37、281、2 630和3 269个。对受到强烈正向选择的基因进行GO分类和KEGG代谢通路富集分析,GO结果显示中蜂和意蜂的单拷贝同源基因富集的GO条目类型相同;KEGG结果显示中蜂的单拷贝基因仅富集在氧化磷酸化。上述结果表明免疫相关基因数量的差异是二者的球囊菌抗性差异产生的重要原因之一,中蜂在与球囊菌的协同进化过程中可能通过提高能量利用从而限制球囊菌的增殖。本研究结果可为阐明中蜂及意蜂幼虫的球囊菌抗性差异产生的分子机制提供参考。  相似文献   

13.
为研究蒙古马高负荷运动训练前后miRNA差异表达,利用二代测序技术对6匹经4个月高负荷运动训练的蒙古马进行了调查,筛选训练前后发生差异表达的相关基因及其所参与的生物学过程。在训练前与训练后两个时期分别采集肌肉样品,建立蒙古马肌肉miRNA文库,对训练前后表达量发生显著变化的差异miRNA进行靶基因预测,并进行GO功能富集与KEGG Pathway分析。结果表明:1)在蒙古马高负荷训练前后两个文库中,分别筛选得到366和346个miRNA,共表达miRNA达248种,分别预测得到7 620和7 535个靶基因。其中高丰度表达的miRNA有miR-1、miR-378、miR-21、miR-101、miR-133a等;2)靶基因参与调控的生物学过程包括:发育过程、解剖结构发育、蛋白质结合、酶结合等;3)KEGG Pathway结果显示,在210条被注释的信号通路中与运动代谢相关的有:调节肌动蛋白细胞骨架,钙信号通路、心肌收缩等。本研究结果可为探究蒙古马强耐力的基本分子机理提供理论基础,同时也为了解蒙古马运动学特性提供新思路。  相似文献   

14.
为探明骆驼斯氏副柔线虫雌虫和雄虫的转录组差异,挖掘出与性别决定相关的功能基因,采用Illumina HiSeqTM2000分别对其雌虫和雄虫样本进行转录组测序,构建cDNA文库,通过建库质量评估和组装效率评估后对其进行功能注释及相关生物信息学分析。结果显示,比较雌虫和雄虫的转录本,获得6 430个差异表达基因,其中2 131个基因上调;4 299个基因下调。对差异表达基因进行GO功能分类,有341 328和714个Unigenes分别注释到生物学过程、细胞组成和分子功能三大类41个分支中。KEGG数据库分析,有1 116个差异表达基因参与到198条KEGG通路中,显著富集在机体免疫调节,卵母细胞成熟等通路。功能聚类分析中,雌虫样本在生殖发育、产卵、生殖行为中高表达;雄虫在主要精子蛋白、丝氨酸/苏氨酸特异性蛋白磷酸酶中高表达。此外,在雄虫和雌虫中筛选出Tab-3,Tab-5、Fem-1和Mog-3等10种线虫性别决定相关基因,以及8种雄性性别特异表达基因和4种雌性性别特异表达基因,这些基因在虫体性别决定和性别特异性表达中发挥重要作用。本研究利用生物信息学技术预测出雌虫和雄虫差异表达基因在功能分类和代谢通路等方面的生物学特征,挖掘出性别决定基因和性别特异性表达相关基因,为进一步开展斯氏副柔线虫的功能基因组学研究、药物靶点预测和疫苗候选基因筛选提供有价值的数据资源。  相似文献   

15.
旨在通过miRNA测序技术对捻转血矛线虫阿苯达唑给药前后的耐药虫株差异表达miRNA及网络调控分析,获取其转录本中全miRNA表达谱及耐药相关差异表达基因和miRNA,探究miRNA在捻转血矛线虫阿苯达唑耐药虫株给药前后的差异。采用Illumina Hiseq2 000平台对捻转血矛线虫阿苯达唑耐药虫株给药前和给药后的虫体进行Small RNA-seq测序,将Raw Data与参考基因组和miRBase数据库进行比对分析,使用RNAhybrid、Miranda和TargetScan对miRNA靶基因进行预测,筛选捻转血矛线虫耐药虫株给药前后表达量显著差异的miRNA并通过Stem-loop qRT-PCR验证,同时对差异显著的miRNA-mRNA调控关系采用Cytoscape软件进行可视化分析。将靶基因向KEGG和GO两大权威数据库映射,预测和分析靶基因功能注释情况和参与调控的通路情况。结果显示:1)捻转血矛线虫阿苯达唑耐药虫株给药前和给药后共筛选出显著差异表达的miRNA 188个,其中125个上调,63个下调,且qRT-PCR结果与转录组测序结果表达情况一致。2)显著差异表达的144个miRNA共关联到353个miRNA-mRNA负调控关系,注释到46条GO功能条目以及代谢,运输和降解等327条通路。3)筛选出的候选靶基因主要富集在ABC转运蛋白,HIF-1信号通路,cGMP - PKG 信号通路以及P450药物代谢通路,证实在药物胁迫的影响下捻转血矛线虫主要通过代谢和转运的方式参与耐药性调控。本研究为捻转血矛线虫耐药性诊断辅助分子标记的筛选和虫体耐药性调控因素及作用机制的研究提供参考。  相似文献   

16.
为筛选出策勒黑羊种质资源繁殖和抗逆性状分子标记和相关基因,本研究以100只无亲缘关系的策勒黑羊为研究对象,基于Illumina ovine SNP 50 K bead chip进行基因分型;用Plink1.90对基因组数据进行质量控制;对有效SNPs计算哈代-温伯格值,选择从小到大排列后前1%的分子标记,参考绵羊基因组(Oar_v4.0)进行注释,对注释到的基因进行GO和KEGG分析。结果表明:1)共获得480个SNPs, 417个候选基因,其中,繁殖性状基因13个、抗逆性状相关基因30个;2)繁殖性状相关的基因有SOX9、GTF2A1、GNAQ等,抗逆性相关的基因有TMEM154、ZNF70、CREB3L2等。综上,本研究解析了策勒黑羊繁殖和抗逆性状基因的遗传规律,特别是对抗逆性相关候选基因研究,发现策勒黑羊群体携带抗肺炎相关的分子标记,为策勒黑羊保种和改良提供了分子水平的参考。  相似文献   

17.
目的 预测、挖掘和分析涉及赤桉Eucalyptus camaldulensis低温胁迫应答的miRNA,为研究其调控赤桉低温胁迫应答的分子网络奠定基础。方法 采用高通量测序对低温处理组和对照(CK)组的赤桉组培苗茎尖进行小RNA测序。以miRBase21.0、Rfam14.1和巨桉E. grandis基因组为参考数据库,利用Bowtie、miREAP和miRDeep2等软件进行miRNA预测,使用RNAfold对预测到的miRNA前体进行二级结构的折叠;采用psRNATarget预测靶基因,通过DEGSeq包分析差异表达的miRNA,并对它们进行GO注释和KEGG富集分析。结果 在赤桉中,共预测到隶属于54个家族的392个已知miRNA和97个新miRNA;其中,CK组共预测到282个已知miRNA,65个新miRNA;低温处理组共预测到329个已知miRNA,51个新miRNA。挖掘到80个在低温处理下显著差异表达的miRNA,包括55个上调和25个下调。GO基因功能注释和KEGG富集分析的结果表明,这些差异表达miRNA可能通过参与代谢通路、次生代谢物的生物合成、细胞膜的改变、信号转导和生物调节等响应低温胁迫。此外,还挖掘到25个可能与ICE1-CBFs-COR通路有关的miRNA。结论 借助高通量测序和生物信息学软件初步得到了低温胁迫下差异表达的赤桉miRNA,为进一步分析这些miRNA在赤桉低温胁迫中的分子功能提供一些参考。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号