基于高通量测序技术的敖鲁古雅驯鹿鹿茸转录组分析

试验旨在通过高通量测序技术获得驯鹿的遗传信息和鹿茸组织的转录组特征,并进一步挖掘其中与鹿茸经济性状相关的关键基因信息。对驯鹿鹿茸组织提取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等,它们的功能可能与驯鹿生茸有关。本试验结果不仅对驯鹿的基因组数据进行了补充,而且还初步揭示了生长期驯鹿鹿茸的基因特征,为驯鹿分子遗传学研究、资源保护与利用及种群资源恢复提供了参考。     

基于高通量测序技术的敖鲁古雅驯鹿鹿茸转录组分析

试验旨在通过高通量测序技术获得驯鹿的遗传信息和鹿茸组织的转录组特征,并进一步挖掘其中与鹿茸经济性状相关的关键基因信息。对驯鹿鹿茸组织提取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等,它们的功能可能与驯鹿生茸有关。本试验结果不仅对驯鹿的基因组数据进行了补充,而且还初步揭示了生长期驯鹿鹿茸的基因特征,为驯鹿分子遗传学研究、资源保护与利用及种群资源恢复提供了参考。     

狍茸顶端组织转录组测序及生物信息学分析

为研究狍茸顶端组织的转录组信息,探讨与狍茸生长相关的基因结构及分子机理。使用Illumina Hiseq^(TM) 2000平台进行高通量测序,利用生物信息学方法进行分析。结果表明：通过测序、Trinity软件组装、比对拼接,获得了36 865个Unigenes,平均长932 nt。经过Nr、Swiss-Prot、KEGG、COG、GO等数据库进行比对,共22 983个Unigenes成功比对上Nr数据库,18 273个Unigenes被注释到GO功能分类中,8 668个Unigenes被归类到COG的25个功能类别中。采用FPKM法筛选与狍茸生长发育相关的高表达基因,其中表达量较高的生长因子有8个,转录因子有5个,细胞外基质有8个。利用高通量转录组测序技术(RNA-seq),成功建立了狍茸顶端组织转录组数据库,可为狍茸的发育分子机理研究提供数据支持。     

基于高通量转录组测序的红花籽油对肥育猪肝脏差异表达基因的影响

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)。     

椰子油对育成猪肝脏转录组高通量测序的影响

试验设计了椰子油组和对照组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个。     

梅花鹿茸重相关基因差异表达分析

通过分析不同茸重梅花鹿茸的差异表达基因,筛选可用于区分鹿茸重量的基因。以东丰县梅花鹿为试验对象,同一饲养条件下,采集不同梅花鹿的茸血,进行高通量测序,对差异表达基因进行GO、KEGG富集分析及qRT-PCR验证。结果表明：共有84个差异基因,其中,44个基因上调,40个基因下调。通过GO分析,这些基因共富集到318个条目,其中,分子功能118个、生物过程138个和细胞组分62个;通过KEGG分析,富集于细胞黏附分子(CAMs)、抗原加工递呈、Ⅰ型糖尿病等106条通路。通过P值筛选出8个基因进行qRT-PCR验证,趋势与测序数据相符,TRPC6基因在高产组的相对表达量高于低产组(P<0.01),有望成为筛选种公鹿的候选基因,为后续梅花鹿的良种选育提供理论基础。     

基于RNA-Seq技术的塔里木马鹿毛色相关差异表达基因筛选及初步分析

旨在基于RNA-Seq技术对塔里木马鹿毛色相关基因进行筛选及分析。采用Illumina Hi Seq TM2000测序平台对塔里木马鹿和天山马鹿的皮肤组织进行转录组测序,所得序列经质控、组装后比对到NR、Swiss-Prot、COG、KOG、KEGG、GO和Pfam数据库中注释,并对差异表达基因进行筛选、功能注释和富集分析。结果表明,测序获得25 038个有注释信息的Unigenes,比对分析显示,塔里木马鹿与天山马鹿有922个差异表达基因,其中上调表达基因495个,下调表达基因427个;GO功能富集分析结果显示,568个差异表达基因富集到61个GO条目上,分别参与了生物学过程、细胞组分及分子功能;KEGG代谢通路富集分析发现,在差异表达基因中富集最显著的代谢通路是ECM-受体相互作用。利用实时荧光定量PCR（qRT-PCR）方法分析与塔里木马鹿毛色相关的7个候选基因的转录水平变化来验证转录组测序结果的准确性和可靠性,这些基因的表达趋势与转录组测序结果相一致。ECM-受体相互作用、蛋白质消化与吸收、PI3K-Akt信号通路及与黑色素合成相关的酪氨酸等通路可能与塔里木马鹿的毛色有关;候选基因MITF、Ggt1、VDR、PTPRF、CⅡTA、ARPC5L、POMC等可能在塔里木马鹿毛色形成过程中发挥重要作用。本研究结果为今后塔里木马鹿毛色相关基因的分子调控机制方面及挖掘潜在的新基因提供了丰富的试验数据。     

福清山羊和戴云山羊夏季发情期卵巢组织比较转录组分析

试验旨在从分子遗传学角度探讨戴云山羊夏季维持较高繁殖力的遗传机制,从中发掘有用的遗传变异信息。利用转录组测序方法对福清山羊和戴云山羊夏季发情期卵巢组织进行研究,筛选品种间的差异表达基因（differentially expressed genes,DEGs）,利用基因本体（gene ontology,GO）、蛋白相邻类的聚簇（cluster of orthologous groups of proteins,COG）和京都基因与基因组百科全书（kyoto encyclopedia of genes and genomes,KEGG）数据库对筛选的差异表基因进行功能注释,并进行差异表达基因聚类分析,通过实时荧光定量PCR验证测序结果。结果显示,福清山羊和戴云山羊发情期卵巢的DEGs共357个,其中上调基因221个,下调基因136个。357个DEGs中的287个基因能够被GO数据库注释,78个DEGs能够被COG数据库注释,219个DEGs能够被KEGG数据库注释。KEGG数据库注释结果显示,这些差异基因显著富集的信号通路为蛋白消化吸收（protein digestion and absorption）、吞噬体（phagosome）、肺结核（tuberculosis）、胞外基质受体互作（ECM-receptor interaction）、金黄色葡萄球菌感染（Staphylococcus aureus infection）、同种异体移植物排斥（allograft rejection）等通路。经实时荧光定量PCR验证,所选基因（转录本）表达变化模式与转录组测序结果一致,表明测序结果可靠。本试验利用高通量测序获得大量山羊卵巢转录组信息,有助于从分子水平对山羊卵巢进行深入研究。     

福清山羊和戴云山羊夏季发情期卵巢组织比较转录组分析

试验旨在从分子遗传学角度探讨戴云山羊夏季维持较高繁殖力的遗传机制,从中发掘有用的遗传变异信息。利用转录组测序方法对福清山羊和戴云山羊夏季发情期卵巢组织进行研究,筛选品种间的差异表达基因(differentially expressed genes,DEGs),利用基因本体(gene ontology,GO)、蛋白相邻类的聚簇(cluster of orthologous groups of proteins,COG)和京都基因与基因组百科全书(kyoto encyclopedia of genes and genomes,KEGG)数据库对筛选的差异表基因进行功能注释,并进行差异表达基因聚类分析,通过实时荧光定量PCR验证测序结果。结果显示,福清山羊和戴云山羊发情期卵巢的DEGs共357个,其中上调基因221个,下调基因136个。357个DEGs中的287个基因能够被GO数据库注释,78个DEGs能够被COG数据库注释,219个DEGs能够被KEGG数据库注释。KEGG数据库注释结果显示,这些差异基因显著富集的信号通路为蛋白消化吸收(protein digestion and absorption)、吞噬体(phagosome)、肺结核(tuberculosis)、胞外基质受体互作(ECM-receptor interaction)、金黄色葡萄球菌感染(Staphylococcus aureus infection)、同种异体移植物排斥(allograft rejection)等通路。经实时荧光定量PCR验证,所选基因(转录本)表达变化模式与转录组测序结果一致,表明测序结果可靠。本试验利用高通量测序获得大量山羊卵巢转录组信息,有助于从分子水平对山羊卵巢进行深入研究。     

基于高通量测序的水貂胚胎滞育期和激活期卵巢转录组分析

试验旨在获取胚胎滞育期与激活期水貂卵巢组织的转录组信息，挖掘滞育的胚胎激活前后水貂卵巢差异表达基因（DEGs）及其功能信息，为探讨卵巢信号调控水貂胚胎滞育的分子机制提供参考。随机采集8只健康雌性水貂（胚胎滞育期和激活期各4只）的卵巢组织为样本，利用Illumina HiSeq平台RNA-Seq技术对其进行转录组测序，筛选在水貂胚胎滞育期和激活期卵巢中的DEGs并进行生物信息学分析。结果显示，测序获得661 195 568个raw reads，经过滤后获得650 834 900个clean reads，组装后得到389 895个unigenes，通过与Nr、GO、KOG和KEGG数据库比对，对unigenes进行注释，其中Nr数据库注释了156 419个unigenes；GO数据库注释了122 657个unigenes；KOG数据库注释了58 320个unigenes；KEGG数据库注释了72 653个unigenes。滞育期和激活期卵巢中有1 797个DEGs，与胚胎滞育期水貂卵巢相比，激活期卵巢有1 298个DEGs显著上调，499个DEGs显著下调。GO功能分析发现，DEGs显著富集的生物学过程主要有跨膜信号受体活性、信号受体活性、G蛋白偶联受体活性、酶联受体蛋白信号通路、细胞表面受体信号通路、细胞周期阻滞、芳香酯酶活性。KEGG通路分析发现，水貂滞育期和激活期卵巢中的DEGs显著富集于神经活动配体-受体相互作用信号通路。本研究利用高通量测序技术获得水貂胚胎滞育期与胚胎激活期卵巢的转录组信息，为深入探究水貂卵巢调节胚胎滞育的分子机制奠定了基础。     

青稞响应燕麦镰孢菌侵染的转录组分析

为了明确青稞与燕麦镰孢菌在转录水平上的互作,以抗病青稞NQK-01-03和感病青稞甘青2号为材料,以Illumina HiSeq Xten为平台,对抗感病青稞在接菌前、接菌20d和40d后茎基部的转录组进行测序。结果表明:经过测序质量控制,抗感病青稞茎基部共得到112.97Gb Clean Data。DEGs基因表达注释共获得差异表达基因4 280个,其中,上调基因2 037个,下调基因2 243个。通过BLAST软件与COG、GO、KEGG、KOG、NCBI-NR、Pfam、Swiss-Prot和eggNOG数据库进行比对,28 869条Unigene获得注释信息,NCBI-NR数据库注释到的Unigene数量最多达到3 981条,占全部注释基因的13.79%。差异表达基因的KEGG通路富集Q-value≤0.05主要有5条,分别为苯丙氨酸代谢、异喹啉类生物碱的生物合成、酪氨酸代谢、苯丙基类生物合成和光合作用天线蛋白。抗病青稞NQK-01-03和感病青稞甘青2号之间GO注释系统包含生物学过程、细胞组分和分子功能3个主要分支。eggNOG和NCBI-NR注释的新基因数目最多分别达到1 833和2 396个,所占比例分别为59.22%和77.42%。     

青稞响应燕麦镰孢菌侵染的转录组分析CSCD

为了明确青稞与燕麦镰孢菌在转录水平上的互作,以抗病青稞NQK-01-03和感病青稞甘青2号为材料,以Illumina HiSeq Xten为平台,对抗感病青稞在接菌前、接菌20d和40d后茎基部的转录组进行测序。结果表明:经过测序质量控制,抗感病青稞茎基部共得到112.97Gb Clean Data。DEGs基因表达注释共获得差异表达基因4 280个,其中,上调基因2 037个,下调基因2 243个。通过BLAST软件与COG、GO、KEGG、KOG、NCBI-NR、Pfam、Swiss-Prot和eggNOG数据库进行比对,28 869条Unigene获得注释信息,NCBI-NR数据库注释到的Unigene数量最多达到3 981条,占全部注释基因的13.79%。差异表达基因的KEGG通路富集Q-value≤0.05主要有5条,分别为苯丙氨酸代谢、异喹啉类生物碱的生物合成、酪氨酸代谢、苯丙基类生物合成和光合作用天线蛋白。抗病青稞NQK-01-03和感病青稞甘青2号之间GO注释系统包含生物学过程、细胞组分和分子功能3个主要分支。eggNOG和NCBI-NR注释的新基因数目最多分别达到1 833和2 396个,所占比例分别为59.22%和77.42%。     

基于转录组测序的新城疫疫苗免疫鸡胸腺差异表达基因分析

为研究新城疫(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疫苗后机体免疫应答反应提供了研究基础。     

牦牛皱胃全长转录组测序分析

【目的】获得牦牛(Bos grunniens)皱胃全长转录组数据库,深入挖掘牦牛皱胃功能基因。【方法】采用PacBio Sequel高通量测序系统,使用单分子实时(single molecular real time, SMRT)测序技术对成年牦牛皱胃全长转录组进行测序,对原始数据进行质控和去冗余分析,再比对参考基因组获取过滤后的非重复序列(Unigenes),使用多种生物信息软件对牦牛皱胃全长转录本数据进行功能注释、转录因子注释、编码区预测、简单重复序列(simple sequence repeats, SSR)分析及可变剪接分析。【结果】通过测序共获得14 467 420条子序列,平均子序列长度为3 344.23 bp,质控得到循环一致性序列(CCS)有296 840条,全长非嵌合(FLNC)序列有277 402条,过滤和去冗余后对比参考基因组最终获得8 556条Unigenes。通过与NR、Swiss-Prot、KEGG、KOG、eggNOG、GO和Pfam数据库比对,对Unigenes进行注释,其中NR数据库注释了8 544条Unigenes; Swiss-Prot数据库注释了...     

基于转录组测序的奶牛泌乳初期与干奶期肝脏差异表达基因分析

旨在通过对奶牛干奶期和泌乳初期肝脏组织进行转录组测序和生物信息学分析,探明其不同生理阶段的差异表达基因和差异代谢通路。选取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)。该研究结果为挖掘奶牛产奶性状的候选基因提供了技术依据。     

大型迪庆藏猪不同生长阶段背脂与腹脂脂质代谢差异基因及调控网络分析

旨在筛选大型迪庆藏猪不同生长阶段、不同部位脂质代谢差异的关键基因。本研究选择胎次相同、出生日期相近、体重10 kg左右的大型迪庆藏猪36头,随机分为3组,相同条件育肥,分别在平均体重达40、80和120 kg时屠宰,测定胴体性能,每组采集3头猪的背脂和腹脂进行高通量转录组测序,测序数据经拼接、比对,筛选与脂质代谢相关的显著差异基因并进行GO、KEGG分析、基因互作网络分析。结果表明,40、80和120 kg大型迪庆藏猪腹脂vs.背脂分别筛到486、765和339个差异表达显著基因,随机挑选的EGR2、SOD3等5个差异表达显著基因的qPCR结果与转录组测序结果一致,差异表达显著基因主要富集在肌肉收缩、细胞黏附、间充质细胞增殖正调节等GO条目,心肌收缩、PI3K-Akt信号通路、Hippo信号通路等KEGG通路;40 kg组EGR2、RARRES2、TMOD4和SFRP2基因位于网络核心,EGR2、RARRES2、SFRP2在腹脂上调,TMOD4下调;80 kg组THBS1、PPARA、NRIP1和LPL基因位于网络核心,4个核心基因均在腹脂上调;120 kg组HTRA1、TSHR、LR...     

基于RNA-Seq高通量测序技术对肉牛下丘脑基因表达分析及新基因挖掘研究

通过对9头不同月龄的安格斯母牛下丘脑组织进行转录组测序和生物信息学比较分析,筛选和挖掘与肉牛生长、繁殖活动相关候选基因和信号通路。通过转录组测序,进行表达基因分析、GO和KEGG富集分析以及信号通路筛选。结果显示,共获得59.74 Gb clean data,各组的样品clean data均达到6.64 Gb, Q30碱基百分比平均在94.60%及以上。有18 979个基因在所有组织中都有表达,有1 706个基因只在2个月龄段下丘脑组织中共同表达,有2 149个基因只在其中的1个月龄段下丘脑组织中表达。在所有下丘脑所表达的基因中,GO、KOG和KEGG数据库得到注释的分别为2 621,12 930和13 453个。在所有KEGG通路中,血管内皮生长因子信号通路、Notch信号通路和Hippo信号通路在各年龄段下丘脑组织中检出频率最高,说明这些信号通路可能通过下丘脑和性腺轴对母牛生殖活动产生影响。     

基于高通量测序的梅花鹿卵巢衰老的转录组分析

【目的】获取青年梅花鹿和老年梅花鹿卵巢组织的转录组信息,为深入研究梅花鹿卵巢衰老的潜在机制提供依据。【方法】以4岁龄青年梅花鹿和10岁龄老年梅花鹿卵巢作为试验样本,对样本进行石蜡切片、苏木精-伊红染色和卵泡计数;用Trizol法提取卵巢样本RNA构建转录组文库,应用Illumina HiSeq^TM2000测序平台测序,对测序结果进行差异表达基因、GO功能、KEGG通路和蛋白互作网络分析,并挑选部分差异表达基因进行实时荧光定量PCR验证。通过RIPA裂解液提取卵巢样本蛋白并对组蛋白H3赖氨酸4三甲基化(H3K4me3)蛋白在青年和老年梅花鹿中的表达情况进行Western blotting检测。【结果】组织学分析结果表明,老年梅花鹿的卵巢卵泡数显著低于青年梅花鹿(P<0.05),而闭锁卵泡数极显著高于青年梅花鹿(P<0.01)。转录组学分析结果表明,老年和青年梅花鹿卵巢间共有1 477个差异表达基因,其中表达上调基因503个,表达下调基因974个。GO功能和KEGG通路富集结果表明,免疫系统过程、转录、DNA损伤修复、炎症反应、凋亡等生物事件与卵巢衰老有关...     

牦牛肺转录组图谱绘制及特有正选择遗传进化基因探究

旨在利用RNA-Seq的高通量测序技术从基因组转录水平解析牦牛低氧适应及相关遗传机制。本研究以牦牛和黄牛肺为样本,经Illumina HiSeq2 500平台测序并进行生物信息学分析。结果表明,在黄牛转录测序文库中共获得了54 720 968条序列,包含4 924 882 170bp,在牦牛中共获得了51 641 282条序列,包含4 647 715 380bp。牦牛基本转录组分析发现,8 123条mRNA的5′-或3′-末端在原有基础上发生了延伸,同时还发现了7 059个新转录本,预测其中2 795个新转录本具有编码蛋白的能力;GO分析表明,共有14 764个基因得到分类注释,涉及细胞组分、生物学过程及分子功能3个大类59个小类,其中细胞、细胞过程及绑定类别最为富集;KEGG分析表明,14 485个基因涉及到258个通路,代谢途径通路最为富集,其次为粘着斑通路及阿米巴病通路。将牦牛肺转录组与黄牛进行比较,Ka/Ks分析筛选出39个正选择基因,这些基因的GO分类注释结果显示,在细胞组分TOP10中与核糖体相关的类别所占比例最大(4/10),生物学过程TOP10中与免疫细胞相关的类别所占比例最大(7/10),分子功能TOP10中与酶活性相关的类别所占比例最大(5/10);KEGG注释结果表明,这些基因共涉及56个通路,最为富集的前10个通路中涉及糖类能量、维生素、核酸等相关代谢的通路占到了7/10。同时发现,最为富集前10个通路中有3个涉及到疾病,推测这与牦牛的自我保护机制有关。本研究通过RNA-Seq转录组测序技术获得了牦牛肺正常转录组数据库,描绘出了牦牛肺正常转录组图谱,为进一步的完善牦牛基因组数据库提供了有价值的数据。同时通过与黄牛转录组的比较,筛选出了相关的正选择基因,为进一步在分子水平揭示牦牛高原低氧适应的独特进化过程及遗传分子机制提供了参考依据。     

利用转录组测序方法研究獭兔毛色相关基因

利用转录组测序方法分析了不同毛色(白色和海狸色)獭兔皮肤组织中基因的表达差异,旨在找到影响獭兔毛色的候选基因。采用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与黑素细胞分化有关,可能在獭兔毛色中发挥作用。