首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
采用高通量测序技术Illumlna HiSeq 2000对高木质素的象草品系eg7和低木质素的象草品系eg87(对照)茎组织进行转录组比较测序。测序获得了169630902个序列读取片段(reads),包含13788439920nt碱基信息。对reads进行序列组装,获得87641个单基因簇(unigene),平均长度580nt。从长度分布、GC含量等方面对unigene进行评估,数据显示测序质量好,可信度高。将获得的unigene与Nr、Nt、Swiss-Prot、COG、GO和KEGG数据库进行序列同源性比较和功能分析,62557个unigene与其他生物的已知基因具有不同程度的同源性,象草与高粱序列同源性最高。共鉴定出33323个差异表达基因,其中上调基因9704个(29.12%),下调基因23619个(70.88%);GO分析显示39968个unigene归为54个功能类别,大量unigene与细胞进程、代谢过程、催化活性等相关;KEGG pathway分析富集得到127条代谢通路,包括光合作用、betalain生物合成、苯丙烷类代谢、苯丙氨酸代谢等,苯丙烷类代谢途径差异基因富集程度高、差异基因数目最多,达285条,该途径中64条木质素单体合成酶基因表达上调,79条ClassⅢ型植物过氧化物酶基因表达下调、22条上调。挑选9个差异基因进行qRT-PCR验证,9个基因的表达趋势与高通量测序结果一致。为象草的分子生物学研究提供了宝贵的基因组数据,对于了解象草茎生物合成与木质素调控基因挖掘和多用途定向育种具有指导意义。  相似文献   

2.
基于高通量测序的海滨雀稗转录组学研究   总被引:1,自引:0,他引:1  
采用新一代高通量测序技术Illumina HiSeq 2000对海滨雀稗叶片转录组进行测序,结合生物信息学方法开展基因表达谱研究和功能基因预测。通过测序,获得了47520544个序列读取片段(reads),包含了4752054400个碱基序列(bp)信息。对reads进行序列组装,获得81220个单基因簇(unigene),平均长度1077 bp,序列信息达到了87542503 bp。另外从长度分布、GC含量、表达水平等方面对unigene进行评估,数据显示测序质量好,可信度高。数据库中的序列同源性比较表明,46169个unigene与其他生物的已知基因具有不同程度的同源性。海滨雀稗转录组中的unigene根据GO功能大致可分为细胞组分、分子功能和生物学过程三大类48个分支,其中有大量unigene与代谢进程、结合活性和细胞进程相关。将unigene与COG数据库进行比对,根据其功能大致可分为25类。KEGG 数据库作为参考,依据代谢途径可将unigene定位到112个代谢途径分支,包括苯丙氨酸代谢通路、植物与病原物互作、植物激素生物合成和信号转导、黄酮类化合物合成、萜类骨架生物合成、脂类代谢、RNA降解等。SSR位点查找发现,从81220个unigene中共找到22721个SSR位点。SSR不同重复基序类型中,出现频率最高的为A/T,其次是CCG/CGG和AGC/CTG。本研究首次对海滨雀稗转录组进行了分析,为草坪草的分子生物学研究提供了宝贵的基因组数据来源。  相似文献   

3.
18头肥育猪按饲粮均分为2组,添加红花籽油组和不添加油脂组(对照组)。试验期60d结束,全部屠宰,分别采集肝脏组织样品,进行高通量转录组测序,找出2处理组间的差异表达基因和差异代谢通路。利用Illumina HiSeqTM2500高通量RNA-seq测序技术对2组肥育猪肝脏RNA测序,使用TopHat2软件将测序得到的reads序列与猪(Sscrofa10.2)参考基因组序列比对,找出差异表达基因,在Nr、GO和KEGG数据库中进行功能注释、富集分析和聚类分析。结果显示,红花籽油组和对照组肝脏中差异表达基因共有1 114个,与对照组相比,红花籽油组表达上调的基因有579个,下调的基因有535个。GO功能分类注释到细胞组成、生物学过程和分子功能数据库中的差异表达基因分别有902,903,857个。注释到KEGG通路中差异表达基因414个,显著富集通路4个(P0.05)。  相似文献   

4.
为了明晰沙鞭(Psammochloa villosa)的转录组特征,本研究利用PacBio Sequel测序平台首次对其进行全长转录组测序和数据分析,结果共获得323 309个clean reads,自我矫正产生环形一致性序列(CCS)673 540个,预测得到蛋白编码区(CDS)序列28 447个;MISA软件共搜索得到沙鞭93 563个简单重复序列(SSR),分布于56 824条unigene上;转录本基因功能注释,共有166 541条序列得到NR注释,结果显示沙鞭与二穗短柄草(Brachypodium distachyon)亲缘最近;KOG数据库比对将97 892个unigene分为25个功能类别,其中注释较多的功能为翻译后修饰、蛋白质周转和分子伴侣;GO数据库比对共注释到87 930个unigene,分为生物过程、细胞组成和分子功能3大类及62个亚类分支;KEGG结果表明碳水化合物代谢、信号转导、能量代谢通路中注释基因较多。本研究结果丰富了沙鞭的遗传信息,为今后沙鞭关键耐旱基因的挖掘提供了理论依据。  相似文献   

5.
试验设计了椰子油组和对照组2个处理。采集育成猪肝脏组织,进行转录组高通量测序,找出2种处理间的差异表达基因及差异代谢通路。利用Illumina Hi SeqTM2 500高通量RNA-seq测序技术测序,使用Top Hat2软件将测序得到的reads序列与猪参考基因组(Sscrofa10.2)序列比对,找出差异表达基因,并在Nr、GO和KEGG数据库中进行功能注释、富集分析和聚类分析。结果表明:椰子油组和对照组肝脏中差异表达基因共有487个,与对照组相比,椰子油组表达上调的基因有222个,下调的基因有265个。在差异表达上/下调前15位基因中,有多个基因是参考基因组中没有的新基因和功能未知基因。GO功能分类注释到细胞组成、分子功能和生物学过程数据库中的差异表达基因数分别有398、368个和398个;注释到KEGG通路中差异表达基因数188个。  相似文献   

6.
7.
《中国兽医学报》2019,(12):2458-2466
旨在通过对奶牛干奶期和泌乳初期肝脏组织进行转录组测序和生物信息学分析,探明其不同生理阶段的差异表达基因和差异代谢通路。选取305 d产奶量接近、干奶天数一致、怀孕天数接近的中国荷斯坦奶牛3头,分别在干奶期(产前50 d)和泌乳初期(产后15 d)活体采集肝脏组织,利用Illumina Hiseq~(TM)2500高通量RNA-seq测序技术对干奶期和泌乳初期肝脏RNA测序,使用TopHat2软件将测序得到的reads序列与牛(UMD3.1.66)参考基因组序列比对,在GO和KEGG数据库中进行聚类分析和富集分析。结果显示,泌乳初期和干奶期肝脏中差异表达基因共有409个,与干奶期相比,泌乳初期表达上调的基因有235个,下调的基因有174个。GO功能分类注释到生物学过程、细胞组成和分子功能数据库中的GO条目分别有63,41,15条。注释到KEGG的显著富集通路有21条(P0.05)。该研究结果为挖掘奶牛产奶性状的候选基因提供了技术依据。  相似文献   

8.
旨在对美国短毛黑水貂快速生长过程中转录组差异表达进行分析。本试验将美国短毛黑水貂45、90日龄各3只健康公貂的胸肌组织作为试验材料,利用Illumina HiSeq~(TM)2500高通量测序平台建立转录组文库,对两个生长阶段差异显著性的基因进行GO功能富集和KEGG Pathway分析,寻找调控水貂快速生长期的相关候选基因和代谢通路,并利用qRT-PCR对转录组测序结果进行验证。结果表明,以P0.05,|log_2FC|1为条件筛选到279个显著差异基因,对这些差异基因进行GO和KEGG分析,筛选到与肌肉生长相关显著富集GO条目有66条、相关差异表达基因44个,其中上调20个,下调24个;与肌肉生长相关显著富集KEGG通路12条,如p53信号通路、PPAR信号通路、FOXO信号通路等。这些差异基因可能在水貂的快速生长过程中起到调控作用,可以作为后续研究水貂生长发育的候选基因。本研究结果为探究水貂生长发育调控机制及培育大体型水貂品种奠定基础。  相似文献   

9.
利用转录组测序方法分析了不同毛色(白色和海狸色)獭兔皮肤组织中基因的表达差异,旨在找到影响獭兔毛色的候选基因。采用Illumina Hi Seq TM2500测序平台对不同毛色獭兔耳组织中所有m RNA进行高通量测序,所得序列经质控、组装后比对到GO、COG、SWISS-PROT数据库中注释,并进行差异表达基因聚类分析和KEGG通路分析。结果表明:测序共获得Unigene 275 625个,差异表达基因12 408个,其中上调表达基因4 904个,下调表达基因7 504个。KEGG通路分析发现,这些差异基因显著地富集在13个代谢通路中,并在细胞色素代谢通路(metabolism of xenobiotics by cytochrome P450)中有8个差异表达的基因,其中CP2F1和CP2BB属于细胞色素P450家族成员,而GSTA4与黑素细胞分化有关,可能在獭兔毛色中发挥作用。  相似文献   

10.
试验旨在通过高通量测序技术获得驯鹿的遗传信息和鹿茸组织的转录组特征,并进一步挖掘其中与鹿茸经济性状相关的关键基因信息。对驯鹿鹿茸组织提取RNA,检测RNA纯度、完整性及是否有污染,合格后构建文库,文库构建成功需要进行库检,库检合格后使用Illuminanovaseq平台进行转录组测序分析。利用Diamond、HMMER、KAAS、Blast2GO等软件对驯鹿鹿茸转录组序列进行了功能注释、基因丰度分析和其他分析。结果显示,测序获得47 818 202条raw reads,经过滤和质量检查得到了46 964 478条clean reads,组装后得到49 171个Unigenes,表明本次测序获得了高质量的鹿茸转录组。基于序列一致性分析,Unigenes与NR、Swiss-Prot、KOG、KEGG、Pfam、GO数据库比对注释成功总共占89.54%,共44 029条序列。Corset聚类后得到49 171个Unigenes,在这些Unigenes中共找到13 795个SNPs位点和18 446个SSRs位点。GO注释结果显示22 955个Unigenes得到注释,占总注释结果的46.68%。与KEGG数据库进行比对分析,发现49 171个Unigenes可能参与或涉及代谢途径,其中20 258个Unigenes共注释到32个KEGG代谢通路。此外,在转录组结果中发现了一些高丰度的基因,如OPNTMSB4、TMSB10等,它们的功能可能与驯鹿生茸有关。本试验结果不仅对驯鹿的基因组数据进行了补充,而且还初步揭示了生长期驯鹿鹿茸的基因特征,为驯鹿分子遗传学研究、资源保护与利用及种群资源恢复提供了参考。  相似文献   

11.
绵羊肺腺瘤病(OPA)是由绵羊肺腺瘤病毒(JSRV)引起的绵羊传染性肺肿瘤疾病。为探究JSRV对绵羊肺脏的致瘤机制,本研究从自然感染JSRV的羊(OPA患羊)肺肿瘤组织和健康羊肺组织中提取总RNA,构建二者的cDNA文库后采用Illumina Hi Seq 4000高通量测序平台进行转录组学测序(RNA-Seq),采用DESeqR筛选健康绵羊与患病绵羊肺组织中的转录差异基因,并以P<0.05和log2(Fold change)≥1筛选转录显著差异基因。通过GO和KEGG数据库对转录显著差异基因进行GO功能和KEGG信号通路的富集分析,并采用RT-qPCR对随机选择的10个转录显著差异基因进行验证。结果显示,与对照组相比,OPA羊肺肿瘤组织中共筛选到1 360个转录上调基因和783个转录下调基因,其中154个转录显著上调基因,212个转录显著下调基因。GO功能分析显示,转录显著差异基因显著富集在178个GO条目中,包括114个生物过程(BP)、19个细胞成分(CC)和45个分子功能(MF),主要涉及生长因子活性、复制后修复、NAD+二磷酸酶活性、核苷代谢过程和...  相似文献   

12.
试验旨在通过高通量测序技术获得驯鹿的遗传信息和鹿茸组织的转录组特征,并进一步挖掘其中与鹿茸经济性状相关的关键基因信息。对驯鹿鹿茸组织提取RNA,检测RNA纯度、完整性及是否有污染,合格后构建文库,文库构建成功需要进行库检,库检合格后使用Illuminanovaseq平台进行转录组测序分析。利用Diamond、HMMER、KAAS、Blast2GO等软件对驯鹿鹿茸转录组序列进行了功能注释、基因丰度分析和其他分析。结果显示,测序获得47 818 202条raw reads,经过滤和质量检查得到了46 964 478条clean reads,组装后得到49 171个Unigenes,表明本次测序获得了高质量的鹿茸转录组。基于序列一致性分析,Unigenes与NR、Swiss-Prot、KOG、KEGG、Pfam、GO数据库比对注释成功总共占89.54%,共44 029条序列。Corset聚类后得到49 171个Unigenes,在这些Unigenes中共找到13 795个SNPs位点和18 446个SSRs位点。GO注释结果显示22 955个Unigenes得到注释,占总注释结果的46.68%。与KEGG数据库进行比对分析,发现49 171个Unigenes可能参与或涉及代谢途径,其中20 258个Unigenes共注释到32个KEGG代谢通路。此外,在转录组结果中发现了一些高丰度的基因,如OPN、TMSB4、TMSB10等,它们的功能可能与驯鹿生茸有关。本试验结果不仅对驯鹿的基因组数据进行了补充,而且还初步揭示了生长期驯鹿鹿茸的基因特征,为驯鹿分子遗传学研究、资源保护与利用及种群资源恢复提供了参考。  相似文献   

13.
为探讨屎肠球菌自溶状态下基因表达变化,利用高通量测序技术进行RNA-Seq分析。差异基因分析显示,自溶状态下共有783个基因表达出现差异,其中在自溶状态下上调和下调的基因分别为689和94个。GO功能聚类分析显示,差异基因主要富集在细胞、细胞组分、生物合成、有机物质生物合成几个方面。KEGG富集分析发现,差异基因主要参与了核糖体、氨基酰基-tRNA生物合成、嘌呤代谢、肽聚糖生物合成、嘧啶代谢、同源重组、氨基酸合成等通路,且下调的差异基因主要富集在碳代谢、糖酵解/糖异生、丙酮酸代谢等能量代谢相关通路。说明屎肠球菌自溶状态下代谢出现了明显改变,对碳源的利用不足和合成能力的下降可能是导致屎肠球菌自溶发生的重要原因。通过分析屎肠球菌自溶状态下基因的表达情况,得到了充分的基因数据,为进一步研究屎肠球菌的自溶发生机制提供了基础。 [关键词] 屎肠球菌|自溶|转录组|差异表达基因|GO功能  相似文献   

14.
木质素是影响生物质降解利用机制的重要组分,但其生物降解机制尚不清晰。利用白腐真菌糙皮侧耳菌株降解玉米秸秆过程中的基因表达变化,尝试挖掘木质素降解的关键生物代谢途径。通过发酵过程中5个木质素含量变化显著关键时间点的转录组测序发现,发酵过程中共有139个上调基因和308个下调基因。基因功能分类数据库(GO)功能分析结果表明,上调基因差异显著性较高的功能为囊泡运输、细胞定位、胞内定位、细胞过程、组蛋白甲基化、转移酶活性、异构酶活性和磷酸葡糖变位酶。基因产物在细胞中代谢途径数据库(KEGG)注释结果表明,上调基因差异显著性较高代谢通路为赖氨酸代谢、β-丙氨酸代谢、淀粉蔗糖代谢、半乳糖代谢和戊糖磷酸。13条与木质素降解直接有关的路径为苯甲酸代谢、二噁英代谢、酚代谢、多环芳香烃(PAH)代谢、氨基苯甲酸代谢、萘代谢、甲苯代谢、脂肪酸代谢、醚酯代谢、二甲苯代谢、碳代谢、过氧化物酶体和溶酶体。  相似文献   

15.
试验旨在探索产肠毒素大肠杆菌(enterotoxigenic Escherichia coli,ETEC)感染猪小肠上皮细胞(IPEC-J2)诱导的microRNA (miRNA)表达谱变化,为解析宿主miRNA在ETEC感染过程中的调控作用提供理论基础。利用Illumina 6000 Novoseq SE50测序平台分别对ETEC感染前后的IPEC-J2进行高通量测序,用Bowtie与参考基因组比对,用DESeq R Package进行miRNA差异性分析。通过miRanda和RNAhybid共同预测差异表达miRNA的靶基因,对差异表达miRNA靶基因进行GO功能和KEGG通路分析。随机选取5个miRNAs,对测序结果进行实时荧光定量PCR验证。结果显示,IPEC-J2在感染前后的sRNA文库经过滤分别得到12 889 260和11 203 056条clean reads。感染前后文库中,miRNA所占比例最高,分别为73.16%和54.10%;分别有97.98%和69.83%长度为18~40 nt的sRNA可比对到参考基因组,表明测序质控良好。长度在22~24 nt的序列大部分首位碱基偏向U,2~8位点出现频率最高的碱基分别为AGCUUAU。共发现311个已知miRNAs,128个新miRNAs。在2个文库中,长度为23 nt的miRNA序列占比最高,分别为41.42%和23.56%。感染后共筛选到140个差异表达miRNAs,其中74个表达上调,66个表达下调。GO分析表明,miRNA靶基因显著富集于代谢过程、正向调节代谢过程、细胞成分或生物合成、免疫系统、细胞内部分和细胞器等功能。KEGG分析表明,差异表达miRNA靶基因显著富集于赖氨酸降解、生产IgA的肠道免疫网络、NF-κB信号通路和T细胞受体信号通路等。实时荧光定量PCR验证结果表明,随机选取的5个miRNAs表达趋势与测序结果一致,表明测序准确可靠。综上所述,IPEC-J2的miRNAs参与了ETEC感染过程,为进一步揭示调控ETEC感染的关键miRNA及其作用机制提供科学依据。  相似文献   

16.
以重组自交构建的燕麦(Avena L.)近等基因系为研究对象,通过构建转录组测序文库并进行Illumina RNA-seq双端测序及生物信息学分析,了解雄蕊生长发育过程,特别是花粉发育基因调控过程。共获得原始数据23.13 Gb,有效数据20.61 Gb。通过 de novo拼接,获得的 unigene 总长度约为 95.95 Mb,均contig长度为716 nt,N50长度为1.07 kb。转录物注释结果显示,82992条转录物序列具有同源比对信息,COG和GO 功能分类工具分别将已注释转录物序列划分为 25个和56 个功能类别。通过对2个材料间转录物序列及表达水平比较,初步确定在雄蕊形成过程中3643条序列表达量上调,50条序列表达量上调10倍以上;12312 条序列表达量下调,2041条序列表达量下调10倍以上。  相似文献   

17.
为了寻找不同品种马骨髓之间的差异表达基因,本研究构建了蒙古马-骨髓和纯血马-骨髓基因表达谱文库,经Illumina Miseq高通量测序,分别获得了7 219 956和7 116 170个有效reads;两个文库映射到11、8和1号染色体的reads最多,而映射到29、30和31号染色体的reads最少;两个文库共筛出318个差异表达基因,其中与免疫相关的基因有21个。GO功能富集分析结果表明差异表达基因与转运、细胞间信号传导、生物学过程和细胞质等有关;KEGG通路富集分析结果表明,差异最大的显著富集KEGG条目为神经系统。这些结果为研究马匹抗病力的遗传和免疫机理等提供基础数据。  相似文献   

18.
19.
为研究新城疫(ND)疫苗接种的青脚麻鸡胸腺组织相关基因与调控通路变化,对无特定病原体(SPF)的青脚麻鸡人工接种ND疫苗48 h后,取胸腺组织采用Illumina HiSeq 2500高通量测序技术进行转录组测序,使用TopHat2软件对测序得到的reads序列与鸡(Gallus_gallus-5.0)的参考基因组比对,筛选出差异表达基因后,在GO和KEGG数据库中进行功能注释、富集分析和聚类分析。结果显示:接种ND疫苗后,共筛选到1 514个差异表达基因,其中1 016个基因上调,498个基因下调;GO功能分类注释到生物学过程、细胞组成和分子功能数据库中的差异表达基因分别有823、815和985个,注释到KEGG通路中差异表达基因1 274个,显著富集通路52个(P0.05);在差异基因中筛选出70个与鸡先天免疫相关的基因,并使用STRING构建蛋白质相互作用网络。研究提示:ND疫苗免疫接种后激活了青脚麻鸡先天免疫相关基因或通路,改善了体内的抗体水平,对探索鸡的宿主-病原体相互作用具有重要意义,为进一步探究青脚麻鸡接种ND疫苗后机体免疫应答反应提供了研究基础。  相似文献   

20.
本文旨在通过转录组测序和生物学信息分析探究在鸡的成骨发育过程中可能发挥重要作用的调控基因和生物学通路。从18胚龄明星肉鸡的胫骨获取骨髓间充质干细胞(BMSC),并进行体外培养与染色鉴定。试验分为试验组和对照组,每组3个生物学重复,试验组将BMSC进行体外成骨诱导分化21 d,对照组不进行处理。分别提取试验组和对照组细胞样品的总RNA,进行转录组分析。结果表明:每个样品平均获得了23 M clean reads数,大约20 M clean reads被比对到鸡的参考基因组上。相比于对照组细胞,试验组筛选到2 715个差异表达基因,包含1 437个上调的差异基因,1 278个下调的差异基因。GO富集分析表明,差异基因主要富集在骨骼系统发育、骨矿化、骨化以及整合素结合等相关功能;KEGG通路主要富集在细胞的焦点黏着、物质合成及能量代谢以及细胞外基质受体相互作用等相关通路。在20条骨发育相关的GO条目或KEGG通路中总共富集到182个基因,其中上调的基因81个,下调的基因101个。通过查找NCBI、Web of Science等相关数据库以及文献,对BMSC成骨分化的过程中发挥重要作用的基因进...  相似文献   

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

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