基于转录组测序挖掘广西麻鸡生长发育相关基因

【目的】通过对1和60日龄广西麻鸡腿肌进行转录组测序(RNA-Seq),筛选出广西麻鸡生长发育关键基因。【方法】选择1和60日龄健康母鸡各3只,分别采集腿肌组织,利用Illumina Novaseq^TM 6000平台进行mRNA转录组测序,利用DESeq2软件进行差异表达基因的分析,并对差异表达基因进行GO功能和KEGG通路富集分析,进一步筛选出生长发育相关基因,并用实时荧光定量PCR验证转录组测序结果。【结果】转录组测序结果显示,6个样本分别获得34 296 198~41 647 642条clean reads。与1日龄腿肌相比,60日龄腿肌共获得2 304个差异表达基因,其中998个基因上调,1 306个基因下调。GO功能富集分析发现,共获得富集条目572条,其中生物过程381条,细胞组分70条,分子功能121条。KEGG通路富集分析发现,共获得34个显著富集的信号通路,其中心肌收缩、MAPK信号通路、紧密连接、肌动蛋白细胞骨架信号通路、心肌细胞中的肾上腺素能信号传导等与生长发育相关。通过GO功能和KEGG通路富集分析获得5个生长发育相关基因,分别为肌球蛋白重链10(MYH10)、肌球蛋白链15(MYH15)、成纤维细胞生长因子10(FGF10)、成纤维细胞生长因子16(FGF16)、肌肉生长抑制素(GDF8)基因,其中MYH10、MYH15、FGF10为上调差异表达基因,FGF16和GDF8为下调差异表达基因。所选差异表达基因实时荧光定量PCR结果与转录组测序的基因表达水平一致,说明转录组测序结果可靠。【结论】本试验通过对广西麻鸡1和60日龄2个不同生长阶段进行转录组测序分析,获得MYH10、MYH15、FGF10、FGF16、GDF8共5个与生长发育相关的关键基因,为后续进一步探讨广西麻鸡生长发育的分子调控机制提供参考。     

基于RNA-Seq数据分析葡萄糖诱导绵羊卵泡颗粒细胞衰老的差异表达基因

【目的】 探讨葡萄糖对绵羊卵泡颗粒细胞衰老及其基因表达的影响。【方法】 将原代分离培养的绵羊卵泡颗粒细胞分别用含17.5(H组)和2 mmol/L(L组)葡萄糖的培养基进行体外培养, 细胞处理72 h后, 利用β-半乳糖苷酶染色检测细胞衰老, 采用RNA-Seq技术进行转录组测序, 对差异表达基因进行GO功能富集分析和KEGG通路富集分析, 并用实时荧光定量PCR验证锚定蛋白重复域蛋白1(ANKRD1)、白细胞介素8(IL-8)、催产素(OXT)、烟酰胺腺嘌呤二核苷酸磷酸氧化酶4(NOX4)基因的表达量。【结果】 H组β-半乳糖苷酶染色阳性细胞率极显著高于L组(P<0.01)。转录组测序结果显示, 两组间存在401个差异表达基因, 其中上调基因153个, 下调基因248个, 高糖诱导的细胞异常相关基因9个, 细胞周期相关基因6个。GO功能富集分析发现, 差异表达基因参与了细胞过程、生物调节、代谢过程、多细胞生物过程及细胞运动等生物过程, 在细胞成分上主要富集在胞外区、膜、突触、细胞器及超分子复合体等, 在分子功能主要富集在催化活性、整合、转运活性、分子功能调节剂及结构分子活性等。KEGG信号通路富集分析发现, 差异表达基因主要富集在糖尿病并发症中的晚期糖基化产物(AGE)-晚期糖基化终末产物受体(RAGE)信号通路、肿瘤坏死因子(TNF)信号通路、过氧化物酶体增殖物激活受体(PPAR)信号通路等。实时荧光定量PCR结果表明, 与L组相比, H组IL-8基因表达水平极显著上调(P<0.01), ANKRD1基因表达水平显著上调(P<0.05), OXT基因表达水平极显著下调(P<0.01), NOX4基因表达量呈上升趋势, 但差异不显著(P>0.05), 实时荧光定量PCR结果与转录组测序结果一致。【结论】 17.5 mmol/L葡萄糖可诱导绵羊卵泡颗粒细胞衰老及衰老相关基因表达变化, 为葡萄糖诱导颗粒细胞衰老的功能和分子机制提供了依据。     

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

【目的】研究不同氧浓度下猪睾丸间质细胞表达谱差异,筛选缺氧导致猪睾丸间质细胞损伤的关键基因。【方法】将猪睾丸间质细胞分为缺氧组(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...     

云南半细毛羊精子冷冻前后转录组表达差异分析

【目的】 本研究旨在对云南半细毛羊精子冷冻前后差异表达的基因与精子冷冻损伤的关系进行研究。【方法】 采用电刺激法采集3只云南半细毛羊精液,各分成2份,一份直接用于提取精子RNA,另一份经过冷冻保存7 d后再进行精子RNA提取。利用获得的精子RNA,采用单细胞转录组测序技术分别检测新鲜精子和冻融精子的mRNA转录组,构建mRNA转录组文库,用Hiseq2500测序仪进行测序,对测序结果做进一步过滤处理,然后对差异表达基因进行GO功能注释和KEGG通路分析。【结果】 6份精液样本共产生623 399 754条原始测序序列,过滤处理后得到514 313 802条（82.50%）高质量的Clean reads。生物信息学分析结果表明,冷冻前后差异表达显著的基因共1 213个。其中,解冻后表达量下调的有1 208个,上调的有5个。通过GO功能注释和KEGG通路富集分析,差异表达基因分别归类于56个GO分类,参与40个信号通路。其中与精子功能密切相关的主要包括精子细胞过程、代谢过程、应激反应、生殖、质膜和氧化还原反应;精子冷冻前后差异表达基因参与的信号通路主要有炎症、细胞因子-细胞因子受体相互作用和细胞凋亡等通路。【结论】 云南半细毛羊精子的冷冻损伤可能与获得的差异表达基因参与的生物过程有关。此外,获得的差异表达基因也可作为分子标记物应用于绵羊冻精品质的评价。     

绵羊超数排卵卵巢差异表达基因及生物信息学分析

【目的】筛选出影响经超数排卵处理的多浪羊卵巢卵泡发育的关键基因、生物过程及信号通路,为进一步了解多浪羊超数排卵状态下卵巢生长发育机制提供依据。【方法】对15只2～3岁身体状况良好的多浪羊进行超数排卵处理后,采集超数排卵中早期(第11天)、中期(第13天)、晚期(第15天)卵巢组织,利用高通量测序技术进行转录组测序,将测序得到的reads序列与参考基因组进行序列比对,筛选出差异表达基因;通过STRING数据库对差异表达基因进行蛋白互作网络分析,采用Cytoscape软件进行可视化编辑处理,利用eggNOG-mapper软件进行GO功能和KEGG通路富集分析。【结果】超数排卵卵巢第11和13天共筛选出223个差异表达基因,获得259条互作关系,差异表达基因显著注释在细胞生长、核糖核蛋白复合等条目中,其中显著差异表达的基因为ITGB3和SNAI1,未发现显著富集通路。超数排卵卵巢第13和15天共筛选出694个差异表达基因,获得785条互作关系,差异表达基因显著注释在DNA复制的启动、DNA复制、DNA新陈代谢过程等条目中,其中显著差异表达的基因主要有CASP3、NOTCH1、WNT2和RUN...     

野生盘羊×巴什拜羊杂交二代尾部脂肪差异表达lncRNA的鉴定与分析

【目的】 分析不同脂尾型绵羊尾部脂肪组织中的长链非编码RNA(long non-coding RNA,lncRNA)的表达谱,探究lncRNA与绵羊尾部脂肪沉积机制的关系,为揭示绵羊尾脂沉积机制提供理论依据。【方法】 选择新疆地方绵羊巴什拜羊(肥尾型)和野生盘羊×巴什拜羊杂交二代(小尾型)作为研究对象,采集2个群体尾部脂肪组织,利用转录组测序技术筛选差异表达lncRNA并预测靶基因,对其进行GO功能和KEGG通路富集分析;利用实时荧光定量PCR技术验证转录组测序的表达结果。【结果】 通过差异表达分析共筛选出728个差异表达lncRNAs,其中270个表达下调,458个表达上调;通过差异表达lncRNA靶基因的GO功能和KEGG通路富集分析,筛选到634个靶基因参与疾病、免疫、蛋白修饰、细胞代谢等相关功能分类,共涉及76条信号通路,部分lncRNA介导的靶基因SCD、GPAM、THRSP、FASN等与绵羊尾部脂肪沉积相关。实时荧光定量PCR与转录组测序结果趋势相同。【结论】 lncRNA在绵羊脂尾进化中对尾部脂肪沉积具有重要的调控作用,为从lncRNA的角度分析绵羊脂肪沉积提供了理论依据。     

A型塞内卡病毒感染PK-15细胞环状RNA表达谱分析

【目的】分析A型塞内卡病毒(Seneca virus A,SVA)感染猪肾上皮细胞(PK-15)后环状RNA(circular RNA,circRNA)表达谱的变化,探究circRNAs在SVA感染过程中发挥的潜在调控作用。【方法】本研究基于Illumina HiSeq 2000平台对SVA感染及对照PK-15细胞circRNAs转录组进行测序,并对差异表达circRNAs的来源基因进行GO功能及KEGG通路富集分析,通过实时荧光定量PCR对新发现的circRNAs进行表达水平鉴定。【结果】转录组测序结果显示,SVA感染组、PK-15细胞对照组中共有432种circRNAs被发现。在SVA感染PK-15细胞中,87种circRNAs表达水平相对对照组PK-15细胞显著上调,74种circRNAs表达水平显著下调。对差异表达circRNAs的来源基因进行GO功能注释分析表明,SVA感染组差异表达circRNAs来源基因主要富集于核仁、细胞器膜和细胞大分子代谢、发育等生理过程。KEGG通路富集结果显示,SVA感染组差异表达circRNAs来源基因主要富集于胞吞过程、cAMP信号通路、Rap...     

湖羊卵巢不同发育阶段的miRNA鉴定与分析

旨在分析不同发育阶段绵羊卵巢组织中miRNA的表达水平,了解miRNA表达及其调控机理。本研究采用转录组测序技术鉴定1和8月龄湖羊卵巢组织中的miRNA,运用生物信息学方法对差异表达的miRNA进行靶基因预测,并对靶基因进行GO和KEGG功能注释,筛选参与卵巢发育的miRNAs。结果显示,6个样品中检测到的已知miRNA个数分别为2 186、2 091、2 217、2 136、2 176、2 088,新miRNA的个数分别为613、567、668、640、662、450。筛选到701个差异表达novel miRNAs,38个差异表达已知miRNAs。对差异表达miRNAs的靶基因进行GO和KEGG分析,发现主要富集在血管内皮生长因子通路和卵巢类固醇生成等过程。对3个差异miRNAs进行qRT-PCR验证,与测序结果一致。oar-let-7b、oar-miR-29a、oar-miR-134-3p、oar-miR-541-3p等的靶基因富集在TGF-beta通路、卵巢卵泡生长等与卵巢发育相关的生物过程。oar-miR-148a、oar-miR-136的靶基因调控颗粒细胞增殖和类固醇生成,可能参与调控绵羊卵巢发育。     

伊犁鹅产蛋前后各期卵巢转录谱的构建及卵泡发育相关基因的分析

【目的】构建产蛋前后各期伊犁鹅卵巢转录组文库,结合生物信息学分析揭示不同产蛋期伊犁鹅卵巢组织差异表达基因,并鉴定出影响鹅卵巢发育的关键基因,为伊犁鹅的繁殖调控提供理论参考。【方法】选择处于开产期（KL组）、产蛋期（CL组）及休产期（XL组）的伊犁鹅各4只,屠宰后取其卵巢组织,用于构建卵巢转录组文库。通过新基因挖掘、基因差异表达、基因注释以及蛋白互作网络分析筛选出与鹅卵泡发育相关的候选基因;随机挑选8个差异表达基因,应用实时荧光定量PCR验证其表达情况。【结果】伊犁鹅卵巢组织学结果表明,在开产期时,伊犁鹅卵巢表面有大量初级卵泡,而产蛋期卵巢则显示出卵泡的层级性,在休产期卵巢中可观察到卵泡出现向内凹陷、闭锁的现象。通过转录组测序（RNA-Seq）,从构建的12个伊犁鹅卵巢cDNA文库中获得有效读数57 811 186~85 328 377条,Q30值均>93.38%,每个样品所产的测序读数于鹅参考基因组上的比对率在82.79%~89.24%。在卵巢中注释的新基因共有1 112个,单核苷酸多态性（SNP）位点总数为1 642 273~2 425 069个;SNP和插入缺失（InDel）均主要注释于内含子区域;各时期的可变剪切类型主要集中为最后1个外显子可变剪切（TSS）及第1个外显子可变剪切（TTS）。在KL vs CL、XL vs CL及XL vs KL组中分别有337、1 136和525个差异表达基因,共有差异表达基因为α2A肾上腺素能受体（ADRA2A）、表皮蛋白（CP）、非转移性黑色素瘤糖蛋白B （GPNMB）及α-1-抗胰蛋白酶样（LOC106033756）。GO功能富集分析发现,差异表达基因主要富集在肽基酪氨酸磷酸化的正调控、细胞黏附及质膜外侧等过程。KEGG通路富集分析发现,差异表达基因主要显著富集于神经活性配体-受体相互作用、ECM-受体相互作用及类固醇生物合成等。结合蛋白互作网络分析,筛选到与鹅卵巢发育相关的潜在调控因子Bruton’s酪氨酸激酶（BTK）、血小板源性生长因子受体α（PDGFRA）、整合素β3（ITGB3）等。实时荧光定量PCR结果显示,RNA-Seq结果准确可靠。【结论】本研究揭示了产蛋前后各期伊犁鹅卵巢组织中的基因表达差异,筛选到影响伊犁鹅卵巢发育的神经活性配体-受体相互作用、ECM-受体相互作用、类固醇生物合成等重要通路与BTK、PDGFRA和ITGB3等关键候选基因,为了解鹅卵巢组织调控产蛋性能的分子机制提供了理论支撑。     

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

通过对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)的表达频率最高,说明这些信号通路可能也参与下丘脑-垂体-卵巢轴对母牛繁殖活动的调节。     

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

试验旨在从分子遗传学角度探讨戴云山羊夏季维持较高繁殖力的遗传机制,从中发掘有用的遗传变异信息。利用转录组测序方法对福清山羊和戴云山羊夏季发情期卵巢组织进行研究,筛选品种间的差异表达基因（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验证,所选基因(转录本)表达变化模式与转录组测序结果一致,表明测序结果可靠。本试验利用高通量测序获得大量山羊卵巢转录组信息,有助于从分子水平对山羊卵巢进行深入研究。     

初生及成年牦牛大脑皮质的转录组分析

旨在分析初生和成年牦牛大脑皮质转录组表达谱,筛选影响大脑皮质发育的候选基因。本研究选取初生(1～7日龄(NYB),n=3)和成年(3～4岁(AYB),n=3)牦牛,放血处死,采集大脑皮质,利用Illumina HiSeq平台对转录组进行测序与分析,筛选得到差异表达基因,在此基础上对差异基因进行GO功能注释、KEGG富集分析和实时荧光定量PCR(qRT-PCR)方法验证。测序结果共分析得到了4 790个显著差异表达的基因(P<0.05),其中包括2 670个上调表达的差异基因,2 120个下调表达的差异基因。随机选取了9个差异基因进行qRT-PCR验证,表达趋势与转录组测序结果一致,表明测序结果可靠。GO功能注释发现其中占比较多的条目为与蛋白合成及ATP利用相关的转录调控、RNA聚合酶II正调控转录、蛋白结合及ATP结合;KEGG富集分析发现占比较多与大脑发育相关的通路为调控轴突发育及重塑的轴突导向通路,与神经胶质细胞增殖相关的MAPK信号通路,与突触传递相关的神经活性配体-受体相互作用、 PI3K-Akt、钙信号及磷脂酶D信号通路,与免疫防御相关的Toll样受体通路,与低氧适应调...     

长途运输影响安格斯育肥牛肝脏功能关键基因的筛选

【目的】研究长途运输对育肥牛体重、屠宰率及肝脏mRNA表达谱的影响,筛选长途运输影响肝脏功能的关键基因。【方法】选用体重相近的安格斯育肥公牛22头,随机分为2组,对照组9头就地屠宰,试验组13头经长途运输后屠宰,观察长途运输对体重和屠宰率的影响;此外,两组分别采集3头体重相近牛的肝脏组织,进行转录组测序(RNA-Seq),筛选差异表达基因,并对其进行GO功能和KEGG通路富集分析。【结果】长途运输后牛活体重降低4.5%,屠宰率较对照组无显著差异(P>0.05);RNA-Seq结果发现,长途运输组肝脏组织有182个基因显著上调表达,91个基因显著下调表达。GO功能与KEGG通路富集分析显示,182个上调表达基因显著富集在刺激反应的调节、过氧化物酶体和氧化还原酶活性,以及cAMP、MAPK和长寿调控通路等,与能量代谢及肝细胞损伤密切相关。根据GO功能与KEGG通路富集分析结果,筛选出4个与育肥牛长途运输后肝脏能量代谢变化、肝脏细胞损伤及损伤后适应性保护机制密切相关的基因：ACOX3、FOS、HSPA1A、HSPA6。【结论】安格斯育肥牛经长途运输后,活体重较运输前降低4.5%,筛选得...     

基于转录组数据挖掘牦牛皱胃发育代谢的关键候选基因

旨在比较不同年龄段牦牛皱胃指数测定、转录组表达谱的变化,挖掘影响牦牛皱胃发育代谢的信号通路和关键基因,为进一步探讨牦牛皱胃发育机制提供理论基础。本研究对1日龄、20日龄、60日龄、15月龄与3岁的牦牛皱胃组织进行转录组测序,对转录组数据进行质控、比对、差异表达基因筛选,并对差异表达基因进行GO功能注释和KEGG富集分析。为进一步验证测序数据的可靠性,随机选取5个差异表达基因进行qRT-PCR验证。结果显示,在发育过程中,瘤胃重的增长最快,其次是瓣胃、网胃和皱胃。转录组结果显示,以1日龄组为对照,在20日龄、60日龄、15月龄和3岁组中分别鉴定到1 310、1 715、1 931和2 199个差异表达基因,4个组共有基因565个;以前一个时间点为对照,20日龄、60日龄、15月龄和3岁组分别鉴定到1 310、861、569和597个差异表达基因,4个对比组共有基因9个。GO功能注释发现,以1日龄组为对照,20日龄、60日龄、15月龄和3岁组分别富集到1 191、2 578、1 117和2 835条显著条目,不同年龄组的前10条显著性GO条目中有共有的,也有各年龄阶段特有的条目。KEGG富...     

不同腹脂率番鸭回肠转录组差异分析

【目的】解析番鸭回肠组织中与腹脂沉积相关的主效基因。【方法】从2 000只番鸭中随机选取200只70日龄的番鸭进行屠宰,采集腹脂并称重,计算腹脂率。按腹脂率大小排序,选取腹脂率最高的5只和最低的5只番鸭的回肠组织进行转录组测序分析。利用Illumina NovaSeq 6000高通量RNA-Seq测序技术对2种腹脂率番鸭的回肠进行转录组测序,使用RSEM软件将测序得到的reads序列与Trinity拼接的参考基因组进行序列比对,筛选出差异表达基因,利用GO和KEGG数据库进行功能注释、富集分析和聚类分析,并利用皮尔森相关性分析方法分析了腹脂重和腹脂率与差异表达基因的相关性。【结果】高腹脂率和低腹脂率番鸭的回肠差异表达基因共有602个;与低腹脂率番鸭相比,高腹脂率番鸭回肠中有285个基因上调,317个基因下调。GO和KEGG富集分析显示,有5个与脂肪代谢相关的GO条目(脂质应答、类固醇激素介导的信号通路、类固醇激素应答、类固醇激素刺激的细胞应答和脂质的细胞应答)和4个KEGG通路(PPAR信号通路、初级胆汁酸的生物合成、花生四烯酸代谢和胆汁分泌),并筛选出与脂肪代谢显著相关的6个基因(P＜0.05),包括NR5A2、ACBP、ACOX2、FABP2、FABP4和FABP5,其中,ACOX2和FABP5基因与腹脂重和腹脂率呈显著正相关(P＜0.05),ACBP基因与腹脂重和腹脂率呈极显著正相关(P＜0.01)。【结论】不同腹脂率番鸭回肠转录组存在差异,发现5个GO条目和4个KEGG通路,NR5A2、ACBP、ACOX2、FABP2、FABP4和FABP5等基因与番鸭腹部脂肪沉积相关。     

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

试验旨在获取胚胎滞育期与激活期水貂卵巢组织的转录组信息，挖掘滞育的胚胎激活前后水貂卵巢差异表达基因（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显著富集于神经活动配体-受体相互作用信号通路。本研究利用高通量测序技术获得水貂胚胎滞育期与胚胎激活期卵巢的转录组信息，为深入探究水貂卵巢调节胚胎滞育的分子机制奠定了基础。     

高温胁迫下罗非鱼源无乳链球菌的比较转录组分析

【目的】从分子水平了解高温季节罗非鱼无乳链球菌病容易暴发的原因。【方法】在35和28℃分别培养无乳链球菌,提取总RNA并进行转录组测序;利用华大基因Dr.Tom自主分析平台进行差异表达基因分析,并针对差异基因进行了KEGG通路和GO功能的分类和富集;通过实时荧光定量PCR随机检测10个基因的差异表达,确定转录组数据的可靠性。【结果】共得到35和28℃培养条件下的显著性差异表达基因357个,其中上调基因229个,下调基因128个;GO数据库富集能得到多个具有显著性的功能集,多集中在碳源代谢相关的途径,包括碳水化合物分解代谢过程(carbohydrate catabolic process)、糖原代谢过程(glycogen metabolic process)、细胞葡聚糖代谢过程(cellular glucan metabolic process)等;差异基因中发现了多个与无乳链球菌致病性相关的基因,如烷基过氧化氢还原酶亚基、cAMP因子、C5a肽酶等编码基因;实时荧光定量PCR与转录组数据间皮尔逊相关系数R值为0.979(P<0.001),说明转录组数据可靠有效。【结论】高温能够显...     

基于多组学解析布鲁氏菌BvfA对睾丸支持细胞线粒体NADH-泛醌氧化还原酶复合体的影响

【目的】明确布鲁氏菌BvfA对宿主睾丸支持细胞(TM4细胞)的损伤作用,以揭示布鲁氏菌对宿主的致病机制。【方法】利用RNA-Seq双端测序、Label-free蛋白组学技术和LC-MS非靶代谢组学技术分析感染布鲁氏菌S19△bvfA与S19菌株的TM4细胞的组学差异,采用KEGG数据库对差异表达蛋白及差异代谢物通路进行分析,并利用平行反应监测(PRM)对蛋白表达量进行验证。【结果】感染S19△bvfA和S19菌株的TM4细胞差异表达基因共400个,差异表达蛋白共422个,差异代谢产物共271种。多组学联合分析发现,S19△bvfA和S19菌株感染TM4细胞差异表达基因、差异表达蛋白和差异代谢产物共同富集的KEGG通路为内源性大麻素信号通路,在该通路中发现S19△bvfA菌株感染TM4细胞的线粒体NADH-泛醌氧化还原酶复合体(ComplexⅠ)中有20个蛋白表达量下调。PRM验证结果显示,ComplexⅠ中Ndufb11、Ndufa9、Ndufv1、Ndufv2、Ndufs8和Ndufs2蛋白的表达趋势与Label-free蛋白组学技术测定结果一致。【结论】BvfA在布鲁氏菌感染TM4...     

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

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