产蛋和停产武定鸡卵巢转录组SNP分析

为提高武定鸡生产性能,探索繁殖相关差异基因,对270日龄产蛋和停产一周以上的武定鸡卵巢组织进行转录组测序,将获得的数据进行SNP检测,利用GO富集和KEGG富集分析SNP所在基因。结果显示:产蛋武定鸡特有的SNP有116235个,停产武定鸡特有的SNP有156020个,两者共有的SNP有276467个;GO富集分析发现,SNP所在基因大量富集在生物进程中;KEGG富集分析发现,SNP所在基因显著富集在代谢途径、剪接、细胞周期、细胞凋亡、基础转录因子等通路;对GnRH信号通路中的5个基因SNP分析发现,有9个SNP位于外显子区,其中有4个SNP位于编码区并引起了氨基酸的改变。研究结果为武定鸡繁殖性能的相关研究提供了理论依据。     

基于转录组测序鉴定雌性鸡胚性腺不对称发育相关候选基因

本研究旨在通过转录组测序（RNA-seq）鉴定雌性鸡胚性腺不对称发育相关候选基因,为确定雌性鸡胚性腺不对称发育的分子机制提供参考。随机选择6胚龄（E6）和8胚龄（E8）的雌鸡胚胎各18～24只,分离并采集其两侧性腺,每6～8个同侧同期性腺组成一个样品,匀浆后提取总RNA经质控后送商业公司进行转录组测序。每个样品3个生物学重复。通过基因差异表达分析,GO、KEGG信号通路分析和基因集变异分析（GSVA）,筛选雌性鸡胚性腺不对称发育相关候选基因和生物学通路。经实时荧光定量PCR(qRT-PCR)对筛选出的差异表达基因进行验证。转录组差异表达分析结果显示雌性鸡胚左右两侧性腺在E6和E8分别存在195个和315个差异表达基因。GO功能富集和KEGG信号通路富集分析发现两侧性腺差异表达基因主要富集在卵子发生、细胞发育、piRNA代谢过程等生物学过程并筛选出Pou5f3、STK31、DMRTB1等与鸡胚性腺不对称发育相关的功能基因。本研究为进一步研究雌性鸡胚性腺不对称发育的分子调控机制提供了参考。     

低海拔与高海拔斑头雁肺脏的比较转录组学分析

为探究高低海拔斑头雁肺脏差异表达基因和差异代谢通路的变化,基于Illumina NovaSeq 6000测序平台对斑头雁肺脏组织进行转录组测序,分析差异表达基因,利用荧光定量PCR对部分差异表达基因进行验证,进一步对差异表达基因进行GO功能富集分析和KEGG信号通路富集分析。结果:与高海拔组斑头雁肺脏相比,筛选出差异表达基因84个,其中上调差异表达基因51个,下调差异表达基因33个。经GO功能注释后富集到21个显著性GO功能,"生物过程"显著富集到11个条目,涉及代谢过程以及内源性刺激的应答等。"细胞组成"显著富集到1个条目,涉及到细胞外区域。"分子功能"显著富集到9个条目,涉及到酶活性、激素、转录因子、钠离子跨膜转运等。注释到KEGG的显著性富集通路有3条,涉及到糖胺聚糖降解通路、丝裂原活化蛋白激酶信号通路以及丙氨酸、天冬氨酸和谷氨酸代谢通路。本研究结果丰富了斑头雁的基因资源,为该物种海拔适应相关基因的功能研究奠定了基础。     

肉鸡肌肉与脂肪组织基因组差异剪接基因分析

旨在鉴定和分析肉鸡肌肉和脂肪组织基因组差异剪接基因(differential splicing gene,DSG)。本研究采用Illumina Hiseq 2500测序平台对AA肉鸡胸肌和腹部脂肪组织各3个样品进行转录组测序,使用rMATS软件检测样品中可变剪接事件(alternative splicing,AS)和组织间的DSG,并对DSG进行GO功能富集和KEGG通路分析。结果表明,在肉鸡肌肉和脂肪组织中分别检测到(5 966.00±111.66)和(6 757.00±156.51)个可变剪接事件(P=0.002)。5种可变剪接类型中以外显子跳跃(skipped exon,SE)为主,分别占肌肉和脂肪组织中可变剪接事件总数量的54.92%和52.67%。肌肉和脂肪组织基因组间检测到513个显著的DSGs。GO基因富集分析发现,93个DSGs显著富集于细胞组分和分子功能中。信号通路分析发现,31个DSGs富集在肌动蛋白细胞骨架调节、焦点粘连、丙酮酸代谢和细胞内吞作用信号通路中,其中14个DSGs在肌肉和脂肪组织基因组中的表达水平未发生显著变化。本研究通过对肉鸡肌肉和脂肪组织基因组中DSG的鉴定与分析,为研究可变剪接及其调控机制、组织间差异剪接的多样性提供理论依据。     

敲除小RNA gcvB后沙门菌的转录组分析

为筛选与沙门菌小RNA gcvB相关的靶基因,采用高通量的转录组测序(RNA-seq)技术对野生型沙门菌和敲除小RNA gcvB沙门菌的mRNA进行对比分析,全面地预测GcvB的靶基因的数量与种类,通过GO和KEGG数据库对筛选基因进行比对注释,进一步分析筛选基因的主要功能及涉及的生物过程。并利用荧光定量PCR对部分差异表达基因进行验证。结果显示:差异表达基因共1 244个,其中基因表达上调678个,基因表达下调566个。GO富集分析显示,差异表达基因主要涉及氧化还原活性、铁离子结合、运动活性等分子功能,细菌鞭毛的细胞成分,细胞呼吸、钴胺素代谢过程和钴胺素生物合成过程等生物过程。KEGG数据库分析显示,差异表达基因主要参与细菌趋化性、卟啉和叶绿素代谢、不同环境中的微生物代谢、双组分系统、甲烷代谢等14个信号通路。对筛选出的10个差异表达基因进行荧光定量PCR验证,结果显示:差异表达基因相对表达量与RNA-seq测序结果基本一致。本研究为进一步探明小RNA gcvB与靶基因相互作用、小RNA的调控机制,以及沙门菌致病机制奠定了基础。     

肾上腺素对民猪脂肪细胞基因表达模式的影响

为了探究冷刺激对民猪脂肪细胞代谢的影响,试验以体外培养的前脂肪细胞为试验材料,通过添加终浓度为1 000 nmol/L的肾上腺素模拟冷刺激环境,同时设置未添加肾上腺素的对照组,诱导分化7 d后收集细胞进行转录组测序(RNA-seq)分析,筛选出表达水平发生变化的差异基因,对差异表达基因进行GO功能注释和KEGG信号通路富集分析,最后通过qRT-PCR验证试验检测RNA-seq结果的准确性。结果表明：肾上腺素会改变脂肪细胞中基因的表达模式,引起128个基因显著上调,93个基因显著下调,其中游离脂肪酸受体4(FFAR4)基因上调倍数最大,为4.63倍。对差异表达基因进行GO功能注释,在分子功能注释分析中有8个条目显著富集,在生物过程注释分析中有128个条目显著富集,在细胞组分注释分析中无显著富集的条目。对差异基因进行KEGG信号通路富集分析,仅神经活性配体-受体相互作用通路被富集。说明肾上腺素可以改变脂肪细胞原有的代谢模式,与脂肪酸代谢和糖原合成等相关的基因发生了显著变化,同时细胞内外的信号转导也发生了改变。     

美国短毛黑水貂快速生长过程中转录组差异表达分析

旨在对美国短毛黑水貂快速生长过程中转录组差异表达进行分析。本试验将美国短毛黑水貂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信号通路等。这些差异基因可能在水貂的快速生长过程中起到调控作用,可以作为后续研究水貂生长发育的候选基因。本研究结果为探究水貂生长发育调控机制及培育大体型水貂品种奠定基础。     

不同氧浓度下猪睾丸间质细胞转录组分析

【目的】研究不同氧浓度下猪睾丸间质细胞表达谱差异,筛选缺氧导致猪睾丸间质细胞损伤的关键基因。【方法】将猪睾丸间质细胞分为缺氧组(1%氧气,H24组)、正常组(15.75%氧气,N24组),每组3个重复,在相应条件下分别处理24 h,用CCK8法检测细胞存活率;利用转录组测序技术筛选出缺氧组和正常组睾丸间质细胞中的差异表达基因,利用DESeq 2.0软件进行差异表达基因分析,并对差异表达基因进行GO功能和KEGG通路富集分析,进一步筛选与氧含量相关的调控基因。【结果】缺氧刺激24 h后,与正常组相比,缺氧组猪睾丸间质细胞存活率极显著降低(P<0.01);转录组测序共获得1 654个差异表达基因,其中1 242个基因上调,412个基因下调。GO功能富集分析发现,共获得富集条目8条,其中生物过程1条,分子功能7条;KEGG通路富集分析发现,共获得10个显著富集的信号通路,包括HIF-1信号通路、JAK-STAT信号通路和糖尿病并发症中的AGE-RAGE信号通路等。通过GO功能和KEGG通路富集分析共获得6个与缺氧应激相关的基因,分别为信号传导及转录激活蛋白5(signal transd...     

应用RNA-seq数据开展鸭基因组可变剪接的鉴定与分析

可变剪接是畜禽增加转录组和蛋白质结构与功能多样性的一种重要机制。研究应用前期测定的不同发育时期(2、4、6周龄)北京鸭胸肌和皮脂RNA-seq数据,挖掘北京鸭胸肌和皮脂内可变剪接及其功能注释。结果显示:1 6个样品中共鉴定出30 465个可变剪接,这些事件发生在7 587个基因上,可变剪切类型以外显子跳跃、内含子滞留、可变的5′剪接位点和可变的3′剪接位点等为主;2可变剪接的发生具有组织和发育阶段特异性;3共有可变剪接的基因主要富集于物质合成、物质结合及酶活性相关的GO条目中,胸肌特有可变剪接基因显著富集到骨骼肌发育相关的GO条目中,而皮脂特有的可变剪接基因富集到油脂合成的GO条目中。研究系统鉴定了北京鸭胸肌和皮脂内的可变剪接和基因,为北京鸭分子育种提供理论基础。     

初情期母牛下丘脑-垂体-卵巢转录组分析

通过对3头18月龄安格斯牛下丘脑-垂体-卵巢组织进行转录组测序和生物信息学比较分析,筛选和挖掘与肉牛繁殖活动相关候选基因和信号通路。通过转录组测序,进行表达基因分析、GO和KEGG富集分析以及信号通路筛选。结果显示,共获得61.92 Gb Clean Data,各样品Clean Data均达到6.88 Gb,Q30碱基百分比在92.42%及以上。有15154个基因在所有组织中都有表达,有4780个基因只在2种组织器官中共同表达,有5809个基因只在其中的一个组织器官中表达。下丘脑-垂体-卵巢各组织所表达的基因中,在GO数据库得到注释的分别为2168,2101,13777个,在KEGG数据库中得到注释分别有12477,12446,12176个。在所有KEGG通路中,Notch信号通路(notch signaling pathway)、神经营养因子信号通路(neurotrophin signaling pathway)和血管内皮生长因子信号通路(VEGF signaling pathway)的表达频率最高,说明这些信号通路可能也参与下丘脑-垂体-卵巢轴对母牛繁殖活动的调节。     

富营养饮食重塑代谢性疾病易感猪关键组织mRNA剪接模式研究

旨在探讨富营养饮食对代谢性疾病易感猪关键组织可变剪接（alternative splicing，AS）和剪接因子（splicing factors，SFs）的影响。本研究以健康状态、月龄和体重相同的12头代谢性疾病易感公猪为材料，随机将其分为两组，其中8头猪饲喂高脂高能量饲料2个月作为富营养饮食组（TG-HFHSD），4头猪饲喂正常饲料2个月作为对照饮食组（TG-CD）。对TG-HFHSD和TG-CD的脂肪、肝、肌肉和下丘脑组织进行RNA-seq测序，利用rMATS软件鉴定组间差异AS，对差异剪接基因（differential splicing genes，DSGs）进行GO和KEGG富集分析，并利用DESeq2软件和皮尔逊相关性系数分析肝组织中的SFs。结果显示，脂肪、肝、肌肉和下丘脑中分别鉴定到48 279、39 536、47 888和52 210个AS事件，均以SE类型为主，占比在77.09%～78.29%之间。脂肪、肝、肌肉和下丘脑中的DSGs数分别为528、1 070、570和600个，其中肝增加的DSGs部分主要是RI和MXE两种类型，分别为177和471个，而其它组织仅为12～28个和99～110个。GO和KEGG富集显示，肝DSGs除富集到细胞结构相关的过程外，还显著富集到可变剪接及mRNA代谢、脂质代谢、DNA损伤及修复、补体与凝血级联等过程。肝中筛选到56个差异表达的SFs，其中CLK1、PGC-1ɑ、USP4等与糖脂代谢和代谢性疾病密切相关。本研究表明，富营养饮食显著增加了肝DSGs，通过剪接调控、糖脂代谢以及DNA损伤来影响肝代谢性疾病的发生，将为肝代谢性疾病早期发生机制研究提供参考。     

撒坝猪与大白猪背脂差异circRNA筛选与调控网络分析

为探究撒坝猪背脂中circRNAs的表达模式与背脂沉积的关系，本研究采用RNA-seq和生物信息学方法筛选撒坝猪和大白猪背部皮下脂肪组织差异表达circRNA，使用MiRanda预测circRNA与miRNA的靶向关系，TargetScan和miRDB预测miRNA的靶基因，对靶基因进行GO和KEGG富集分析，构建与脂质代谢相关的circRNA-miRNA-mRNA互作网络。结果显示：以|log2fold change|≥1且Padj<0.05筛选到27个差异表达circRNA，在撒坝猪中上调7个，下调20个；GO和KEGG富集分析显示，靶基因主要富集在PI3K-Akt信号通路等与脂质代谢相关的过程；circRNA-miRNA-mRNA互作网络显示，novel＿circ＿0003253上调，靶向miR-148a-3p和miR-182下调，靶基因ROCK1和SNAP23表达上调，促进脂肪生成相关基因表达，可能导致脂肪沉积速率增加，脂肪生成增多，脂滴和脂肪细胞直径增大；novel＿circ＿0002314下调，靶向miR-29b上调，靶基因SPARC下调，可能促进皮下脂肪细胞增殖。随...     

MSTN基因编辑牛肌肉转录组测序及生物信息学分析

为探究肌生长抑制素（myostatin,MSTN）基因对牛肌肉发育的具体调控机制,本研究选取同一牛场健康鲁西黄牛10头,其中通过转基因技术得到的基因编辑牛MSTN^-/-和同种非转基因野生型牛各5头。分别采集两组牛腿臀肌肉样品,利用IIlumina HiSeq高通量测序技术进行转录组测序分析,通过生物信息学方法比较两组样本间的差异表达基因,并进行GO和KEGG富集分析,最后利用实时荧光定量PCR验证转录组测序数据。结果显示,基因编辑型牛和野生型牛之间共检测到18 071个基因。在log₂|FoldChange|≥ 1.48条件下,筛选出406个差异表达基因,其中347个显著上调,59个显著下调。GO功能富集分析显示,MSTN基因编辑后显著影响915个功能类别（P<0.05）,差异基因主要参与结合、生物系统调节、免疫系统等相关功能。KEGG通路富集分析结果共涉及211个通路,差异基因主要富集在细胞黏附分子、趋化因子信号通路、细胞因子互作等信号通路上,进一步从中筛选出可能参与细胞生长、肌肉发育的差异基因（CD14、KIT、CSF1R、FBP1、DUSP4、ULBP21、PRKCB、SPN、CHAD、SRC）。实时荧光定量PCR检测结果显示,所选差异基因表达水平与转录组表达水平一致,证明测序结果的可靠性。本研究结果表明,MSTN基因发挥作用后可以介导多个下游基因表达,从而影响相关信号通路及生物学过程;同时,所筛选出的差异表达基因可作为进一步研究骨骼肌调控机制的候选靶标。     

基于转录组测序筛选高盐培养下双峰驼肾皮质细胞差异表达基因

为探讨双峰驼的高盐适应机制并筛选出受高盐调控的基因，本试验采用RNA-Seq技术对双峰驼肾皮质细胞进行转录组测序分析。试验分为2个组，等渗培养基处理的肾皮质细胞为对照组（Control），高渗培养基处理的肾皮质细胞为高渗组（HS），进行转录组测序，从而筛选出一系列差异表达基因（DEGs），并用GO富集和KEGG富集分析对DEGs进行基因功能和途径的注释。结果显示，在高盐环境下双峰驼肾皮质细胞中共有4 854个显著DEGs;GO富集分析显示，DEGs显著富集在G蛋白偶联受体活性、受体结合、核酸结合转录因子活性、质膜、离子通道的活动和免疫反应等功能中；KEGG通路富集显示，DEGs显著富集在信号传导、辅助因子和维生素代谢、信号分子相互作用、运输和分解代谢、免疫系统等通路中；采用实时荧光定量PCR技术随机检测了参与高盐代谢的3个差异表达基因，其表达趋势与RNA-Seq一致，可以验证RNA-Seq技术的准确性。因此，双峰驼肾皮质细胞的高盐耐受性可能与这些途径和通路中重要基因的调控有关。本试验结果将为进一步揭示双峰驼耐盐性的分子调控机制提供重要的参考依据。     

共轭亚油酸影响边鸡胸肌脂类代谢的相关差异基因的筛选

旨在通过转录组测序技术筛选出共轭亚油酸(conjugated linoleic acids,CLA)影响边鸡胸肌脂类代谢的相关基因,为日粮中CLA调控边鸡脂类代谢的作用机制提供理论依据。本试验选取90日龄健康边鸡180只,随机分为5组,每组3个重复,每组日粮中分别添加不同比例的CLA 0%(对照组)、0.5%、1.0%、1.5%和2.0%,预饲期1周,正饲期6周。饲养试验结束后采集胸肌组织进行转录组测序,使用DESeq2软件对测序数据进行差异表达分析,对差异表达基因进行GO功能和KEGG通路富集分析,筛选出与胸肌脂类代谢相关的差异表达基因,利用qPCR对差异表达基因进行验证。通过转录组测序数据分析,共筛选出1 229个差异表达基因,其中594个基因上调,635个基因下调。随机选取9个差异表达基因进行qPCR验证,其相对表达量变化趋势与测序结果一致。GO功能分析发现,差异表达基因主要集中在生物学过程中的细胞过程、单一生物过程、生物调节和代谢过程。KEGG通路富集分析发现,各试验组中差异表达基因显著富集在不同的信号通路中,主要富集在ECM-受体相互作用、黏着斑、吞噬体、细胞凋亡、肌动蛋白细胞骨架调节和细胞因子-细胞因子受体相互作用等通路中。通过GO和KEGG功能注释分析,共筛选出MCAT、APOA1、PTGDS、ALDH3A2、PLTP、FABP3、FOXO1、UCP3和FABP4等18个主要的参与脂类代谢相关的差异表达基因,可能在调控边鸡胸肌脂类代谢中发挥重要的作用。本研究通过转录组测序筛选出CLA影响边鸡胸肌脂类代谢过程中的主要调控基因,发现了日粮中CLA影响差异表达基因作用的主要生物学过程和信号通路,为今后研究CLA影响边鸡胸肌脂类代谢的分子作用机制奠定了基础。     

利用RNA-seq技术筛选影响文昌鸡腹脂沉积的重要候选基因

本研究利用RNA-seq技术筛选影响文昌鸡腹脂沉积的重要候选基因，旨在为低脂文昌鸡的分子育种提供理论依据。本研究结果显示，共有15 996个基因在文昌鸡腹部脂肪组织中表达，16 926个基因在儋州鸡腹部脂肪组织中表达，15 580个基因在文昌鸡和儋州鸡腹部脂肪组织中均有表达，95个基因在文昌鸡和儋州鸡腹部脂肪组织中差异表达。GO富集分析结果显示差异表达基因显著富集在细胞过程、分子活性物质运输等过程（P<0.05）。KEGG富集分析结果显示，差异表达基因在PPAR信号通路、TGF-β信号通路、不饱和脂肪酸的生物合成等途径富集。通过筛选与验证，最后确定SCD、SMOC1、IGFBP2、MSMB基因为影响文昌鸡腹脂沉积的重要候选基因。     

基于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等可能在塔里木马鹿毛色形成过程中发挥重要作用。本研究结果为今后塔里木马鹿毛色相关基因的分子调控机制方面及挖掘潜在的新基因提供了丰富的试验数据。     

捻转血矛线虫阿苯达唑耐药株给药前后比较转录组学分析

为探明捻转血矛线虫阿苯达唑耐药株给药前后在转录组水平上的差异,本研究采用Illumina Hiseq4000对给药前后的捻转血矛线虫进行转录组测序,筛选得到差异表达基因并通过GO和KEGG数据库对其进行功能注释及富集性分析,并利用荧光定量PCR验证部分差异表达基因。结果显示,共筛选获得851个显著差异表达基因,其中包括584个上调基因和267个下调基因。通过GO功能富集分析显示,有458、418和367个基因分别注释到生物学过程、分子功能和细胞组分三大类;对差异表达基因进行KEGG富集分析显示,有173个差异表达基因参与到75条KEGG通路中,显著富集在核糖体合成、细胞凋亡、过氧化物酶增殖体激活受体(PPAR)等信号通路。本试验初步筛选了阿苯达唑耐药株给药前后的差异表达基因,为深入探索捻转血矛线虫耐药性分子机制、筛选对耐药性检测的分子标记及耐药性早期鉴别诊断方法的建立提供重要数据。     

通过比较转录组分析鉴定山麻鸭开产调控基因

蛋鸭开产与卵巢发育有紧密联系，本实验旨在探索山麻鸭开产前后卵巢形态变化及基因表达的差异，鉴别参与启动开产的关键基因，为后期蛋鸭开产相关的分子遗传机制研究奠定基础。在山麻鸭开产前后的S1期（81日龄）和S2期（137日龄）分别选取3个个体的卵巢基质进行转录组测序，对转录组数据开展差异表达分析、功能富集分析及关键基因筛选，同时对关键基因进行荧光定量PCR(qRT-PCR)验证。在差异倍数log2转换的绝对值（|log2(FoldChange)|）大于1.5、矫正P值<0.05的筛选条件下鉴别出382个差异表达基因。基于382个差异表达基因的GO功能注释发现与卵泡生长和卵巢功能调节相关的通路：细胞外基质的组织、细胞外基质结构成分和含有胶原蛋白的细胞外基质相关途径；KEGG富集分析发现上述基因显著富集于ECM受体相互作用和PI3K-Akt信号通路。在涉及上述重要富集通路的基因中有9个基因交集：CD44、COL1A1、COL6A1、CREB3L1、ITGB2、THBS2、VWF、VEGFC、WNT5B，通过qRT-PCR验证了这9个基因在转录组数据差异表达方面的可靠性。在这9个基因中，CO...     

牦牛发情期卵巢比较转录组学研究

旨在进一步了解牦牛发情期卵巢的分子机制,解析牦牛繁殖的特殊性。本研究应用RNA-seq技术对牦牛和平原黄牛发情期卵巢进行转录组高通量测序及全基因组差异表达模式比对分析。通过比较分析牦牛和黄牛卵巢转录组数据,共筛选出1 307个差异表达基因,其中661个基因表达量上调和646个基因表达量下调。进一步功能分析表明,这些差异基因涉及多种GO分类及KEGG通路。其中GO分类注释显示,差异基因与细胞粘附、激素调控等生物学过程存在密切关联,同时钙离子结合、阳离子跨膜转运等分子事件表现活跃。KEGG通路分析显示,补体和凝血级联通路的富集水平最高,其次为细胞色素P450相关通路。昼夜节律等一些新型通路,尽管与生殖功能没有明显关联,但也表现出显著富集。本研究首次对比牦牛和黄牛发情期卵巢转录组数据,筛选并分析相关差异基因。该研究结果为进一步阐述牦牛卵巢的基本分子机理提供基础,同时也为全面理解牦牛繁殖特异性提供新的思路。