基于RNA-seq数据鉴定苏姜猪肉质性状相关基因

试验通过对苏姜猪背最长肌和腿肌进行转录组分析,旨在挖掘影响苏姜猪肌肉肉质性状的候选基因。选取90和180日龄苏姜猪各3头,利用RNA-seq技术对90日龄苏姜猪背最长肌、180日龄苏姜猪背最长肌和腿肌进行测序,对差异表达基因进行GO功能和KEGG Pathway显著性富集分析。结果显示,180日龄苏姜猪背最长肌和90日龄苏姜猪背最长肌中有1 655个基因差异表达,其中474个基因表达上调,1 181个基因表达下调,GO功能和KEGG Pathway显著性富集分析发现,474个表达上调基因主要参与肌肉纤维、转录调节、钙信号通路、MAPK信号通路、胰岛素信号通路、FOXO信号通路等生物学功能,1 181个表达下调基因主要参与免疫、造血、溶酶体、绑定、酶活性等生物学功能;180日龄苏姜猪背最长肌和腿肌中有383个基因差异表达,其中70个基因表达上调,313个基因表达下调,GO功能和KEGG Pathway显著性富集分析发现,70个表达上调基因未能显著富集,313个表达下调基因主要参与细胞分化增殖、肌肉发育、细胞外基质、cGMP-PKG信号通路等生物学功能。本试验获得了苏姜猪背最长肌和腿肌组织的转录组信息,筛选出6个与苏姜猪肉质性状相关的候选基因：FOXO1、FOXO3、FOXO4、MYF6、A-FABP和H-FABP,为深入研究苏姜猪肉质性状的分子机理奠定基础。     

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

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

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

旨在通过转录组测序技术筛选出共轭亚油酸(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影响边鸡胸肌脂类代谢的分子作用机制奠定了基础。     

猪肌内脂肪沉积相关基因的筛选及其表达特性分析

【目的】 通过分析马身猪和大白猪背最长肌转录组测序数据,筛选肌内脂肪沉积相关基因,并对其表达特性进行分析。【方法】 采集180日龄马身猪和大白猪背最长肌,采用RNA-Seq技术对其进行转录组测序,筛选差异表达基因并进行GO功能和KEGG通路富集分析;采用GeneCards在线查询基因功能,进一步筛选与肌内脂肪沉积相关的差异基因;采用实时荧光定量PCR技术检测差异基因在2个品种猪不同组织以及肌内脂肪细胞成脂分化过程中的表达变化。【结果】 获得2个品种猪背最长肌差异表达基因共280个,其中128个表达显著上调,152个表达显著下调。GO功能富集分析注释到46个条目,其中生物过程24条,分子功能8条,细胞组分14条。KEGG通路分析显著富集到PPAR信号通路、MAPK信号通路和脂肪细胞因子信号等脂肪沉积相关通路。GeneCards功能查询获得TRAF2、DUSP1、ACOT4、NR4A1、SLC27A6、PLIN5共6个与脂肪沉积相关的基因。实时荧光定量PCR分析表明,马身猪和大白猪背最长肌中NR4A1、DUSP1、PLIN5基因表达量差异显著或极显著(P<0.05;P<0.01),TRAF2、ACOT4和SLC27A6基因表达量差异不显著(P>0.05),但表达变化趋势与测序结果一致;马身猪腹部皮下脂肪组织中NR4A1和PLIN5基因表达量极显著高于背部皮下脂肪组织(P<0.01),而大白猪背部皮下脂肪组织中NR4A1和DUSP1基因表达量极显著高于腹部皮下脂肪组织(P<0.01),PLIN5基因表达量显著高于腹部皮下脂肪组织(P<0.05)。在背部皮下脂肪组织中,NR4A1、DUSP1和PLIN5基因在大白猪中的表达量极显著高于马身猪(P<0.01);在腹部皮下脂肪组织中,NR4A1和DUSP1基因在大白猪中的表达量显著高于马身猪(P<0.05),PLIN5基因表达量在2个品种间无显著差异(P>0.05)。体外培养的猪肌内脂肪细胞分析表明,随着成脂天数的增加,NR4A1基因表达量呈下降趋势,DUSP1和PLIN5基因表达量均呈上升趋势。【结论】 本研究获得了NR4A1、DUSP1、PLIN5 3个与猪肌内脂肪沉积相关的候选基因,可为后续进一步探讨猪肌内脂肪沉积调控机制提供一定的参考。     

基于RNA-seq数据鉴定苏姜猪肉质性状相关基因

试验通过对苏姜猪背最长肌和腿肌进行转录组分析,旨在挖掘影响苏姜猪肌肉肉质性状的候选基因。选取90和180日龄苏姜猪各3头,利用RNA-seq技术对90日龄苏姜猪背最长肌、180日龄苏姜猪背最长肌和腿肌进行测序,对差异表达基因进行GO功能和KEGG Pathway显著性富集分析。结果显示,180日龄苏姜猪背最长肌和90日龄苏姜猪背最长肌中有1 655个基因差异表达,其中474个基因表达上调,1 181个基因表达下调,GO功能和KEGG Pathway显著性富集分析发现,474个表达上调基因主要参与肌肉纤维、转录调节、钙信号通路、MAPK信号通路、胰岛素信号通路、FOXO信号通路等生物学功能,1 181个表达下调基因主要参与免疫、造血、溶酶体、绑定、酶活性等生物学功能;180日龄苏姜猪背最长肌和腿肌中有383个基因差异表达,其中70个基因表达上调,313个基因表达下调,GO功能和KEGG Pathway显著性富集分析发现,70个表达上调基因未能显著富集,313个表达下调基因主要参与细胞分化增殖、肌肉发育、细胞外基质、cGMP-PKG信号通路等生物学功能。本试验获得了苏姜猪背最长肌和腿肌组织的转录组信息,筛选出6个与苏姜猪肉质性状相关的候选基因:FOXO1、FOXO3、FOXO4、MYF6、A-FABP和H-FABP,为深入研究苏姜猪肉质性状的分子机理奠定基础。     

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

【目的】解析番鸭回肠组织中与腹脂沉积相关的主效基因。【方法】从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等基因与番鸭腹部脂肪沉积相关。     

藏鸡与白羽肉鸡垂体转录组差异表达研究

试验旨在通过对藏鸡和白羽肉鸡垂体转录组测序和生物信息学分析,筛选出与鸡生长发育相关的差异表达基因（differentially expression genes,DEGs）及信号通路。本研究选用42日龄健康藏鸡和白羽肉鸡各3只,分别采集垂体组织利用Illumina HiSeq 2000平台进行转录组mRNA测序,对筛选到的DEGs筛选DEGs,GO和KEGG数据库功能注释和富集分析,实时荧光定量PCR验证随机挑选的DEGs表达水平。结果显示,藏鸡和白羽肉鸡分别获得126 094 302和125 666 442条clean reads。与白羽肉鸡相比,藏鸡垂体组织中共有DEGs 392个,其中196个为上调基因,196个为下调基因（P<0.05）,其中包括生长激素（GH）基因、生长激素释放激素受体（GHRHR）基因和胰岛素样生长因子2 mRNA结合蛋白1（IGF2BP1）基因。GO和KEGG富集分析发现,DEGs显著富集于细胞黏附分子信号通路和神经活性配体-受体相互作用信号通路等细胞通讯相关的信号通路,GH、GHRHR和IGF2BP1基因也显著富集于神经活性配体-受体相互作用信号通路。实时荧光定量PCR结果显示,转录组测序结果准确可靠。综上研究,本研究初步揭示了影响藏鸡和白羽肉鸡生长发育速度差异的关键候选基因和信号通路,为了解鸡垂体组织调节生长发育的分子机制提供了理论依据。     

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

基于单细胞转录组测序技术筛选牛体内囊胚性别特异性基因

【目的】探究牛体内雌性和雄性囊胚在mRNA水平上的特异性差异,挖掘雌性和雄性囊胚发育差异相关候选基因。【方法】参考GenBank中牛牙釉质基因（AMEL）序列（AMELX基因,登录号：NM_001014984.1;AMELY基因,登录号：NM_174240.2）设计引物,以胚胎内细胞DNA为模板,采用巢式PCR扩增技术鉴定单个牛体内囊胚的性别。选用确定性别的单个雌性和雄性囊胚为试验材料,采用Smart-Seq2扩增技术构建单个胚胎测序文库,应用Illumina HiSeq Xten高通量测序平台对单个胚胎进行单细胞转录组测序（scRNA-Seq）,并进行了基因差异表达、GO功能和KEGG通路分析。【结果】通过巢式PCR扩增胚胎内细胞中的AMEL基因,确定了胚胎性别;基因表达分析筛选出两组之间差异表达基因（DEGs）6 160个,其中雌性特异性基因675个,雄性特异性基因305个;GO功能注释发现雌性特异性基因显著注释在核苷酸结合、信号转导调控、多细胞生物发育、细胞骨架等条目,雄性特异性基因显著注释在氧化磷酸化、线粒体、线粒体内膜、核糖体等条目;KEGG通路富集分析发现7条与雌性和雄性囊胚发育差异相关的通路,分别为代谢途径、糖酵解、磷酸戊糖途径、细胞衰老、氧化磷酸化、调节干细胞多能性的信号通路和Wnt信号通路;并利用GO和KEGG结果筛选出5个在牛体内囊胚期胚胎发育过程中具有直接或间接作用的基因：FBP1、GADD45G、FHL2、FOSB和WNT2B。【结论】牛体内雌性和雄性囊胚之间存在广泛的转录差异,性别特异性基因FBP1、GADD45G、FHL2、FOSB和WNT2B可能是直接或间接影响胚胎发育的关键基因。     

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

【目的】获取青年梅花鹿和老年梅花鹿卵巢组织的转录组信息,为深入研究梅花鹿卵巢衰老的潜在机制提供依据。【方法】以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损伤修复、炎症反应、凋亡等生物事件与卵巢衰老有关...     

野生盘羊×巴什拜羊杂交二代尾部脂肪差异表达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的角度分析绵羊脂肪沉积提供了理论依据。     

基于猪小肠上皮细胞炎症模型分析产肠毒素大肠杆菌诱导的转录组表达变化

【目的】分析产肠毒素大肠杆菌（enterotoxigenic Escherichia coli，ETEC）诱导猪小肠上皮细胞（IPEC-J2）炎症反应的主要宿主调控基因及相关分子通路，为进一步揭示ETEC感染引发炎症反应的作用机制提供理论依据。【方法】用感染复数（multiplicity of infection，MOI）为100的ETEC菌液感染IPEC-J2细胞，18 h后收集细胞上清，通过ELISA方法检测炎性因子白细胞介素1β（interleukin 1 β，IL-1β）的分泌水平。由Illumina PE150测序平台进行高通量转录组测序，通过edgeR v 3.12.1软件对测序结果进行差异表达分析，利用GOseq和KOBAS 2.0软件对得到的差异表达mRNA进行GO功能和KEGG通路富集分析，随机挑选5条差异表达mRNA，利用实时荧光定量PCR对其验证。【结果】试验成功建立了ETEC感染诱导的IPEC-J2细胞炎症模型。与对照组相比，炎症模型组共发现529条差异表达mRNA，其中236条表达上调，293条表达下调，包括HSP90AA1、CGNL1、CADPS2、EHBP1等免疫相关基因。GO功能分析表明，ETEC感染后差异表达mRNA显著富集于细胞组分和生物过程，包括细胞内成分（intracellular part）、细胞质（cytoplasm）、核成分（nuclear part）、胞内膜结合细胞器（intracellular membrance-bounded organelle）、膜结合细胞器（membrance-bounded orgabelle）及细胞进程（cellular process）等。KEGG通路分析表明，炎症模型中差异表达基因大多数被注释到疾病相关通路及Toll样受体信号通路、mTOR信号通路、细胞周期、黏附连接等免疫相关信号通路。利用实时荧光定量PCR验证随机挑选的5条差异表达mRNA，结果与测序报告中差异表达mRNA变化趋势一致，证明测序结果可信。【结论】HSP90AA1、CGNL1、CADPS2和EHBP1等基因可能在ETEC黏附IPEC-J2细胞并诱导炎症反应中发挥关键作用，并通过Toll样受体、mTOR等信号通路进行免疫调控。研究结果为进一步揭示ETEC感染诱导炎症反应的分子机制及相关调控网络提供参考依据。     

ZBED6基因敲除猪肾脏转录组分析

【目的】 探究锌指BED结构域结合蛋白6(zinc finger BED domain-containing protein 6,ZBED6)基因敲除对猪肾脏组织基因转录表达的影响,并解析ZBED6基因调控猪肾脏代谢的靶基因及其通路。【方法】 对8月龄ZBED6基因敲除巴马香猪(ZBED6 KO)和同月龄野生型巴马香猪(WT)的肾脏组织(n=3)进行总RNA提取,利用实时荧光定量PCR检测胰岛素样生长因子2(insulin-like growth factor 2,IGF2)和差异基因的表达量;以Illumina Hiseq高通量测序技术对ZBED6 KO和WT猪肾脏组织mRNA进行转录组测序(RNA-Seq),用猪Sus scrofa 11.1作为参考基因组序列,通过生物信息学软件TopHat、Cufflink、Cuffmerge和Cuffdiff对测序数据进行分析,筛选ZBED6 KO和WT猪肾脏组织中的差异表达基因,对差异表达基因进行层次聚类和KEGG通路富集分析,并通过实时荧光定量PCR验证RNA-Seq结果中差异表达基因的可靠性。【结果】 实时荧光定量PCR结果显示,ZBED6 KO猪肾脏IGF2基因mRNA表达量显著高于WT猪(P<0.05)。测序共获得78 G数据量,每个样本比对率在87.3%以上,测序质量良好;ZBED6 KO和WT猪肾脏组织之间共检测到25 213个基因,以调整后P<0.05和log₂|FoldChange|>2为筛选标准,得到299个差异表达基因,其中上调的基因103个,下调的基因199个。热图和主成分分析(PCA)结果显示,ZBED6 KO和WT组猪肾脏基因组内表达模式相似,组间表达模式不同;KEGG通路富集分析显示,差异基因主要参与视黄醇及疾病代谢相关通路。ZBED6基因敲除后,其中有9个差异表达基因(CYP2C42、AOX1、ENSSSCG00000036274、RDH16、CYP2A19、ENSSSCG00000022724、CYP26B1、CYP1A1、RDH5)被富集到视黄醇代谢通路中,可能参与调控猪肾脏代谢平衡。实时荧光定量PCR检测发现,7个差异表达基因(CYP2C42、AOX1、RDH16、CYP2A19、CYP26B1、CYP1A1、RDH5)的表达模式与RNA-Seq分析结果一致,证实RNA-Seq结果的可靠性。【结论】 本试验利用RNA-Seq技术分析了ZBED6基因敲除对巴马香猪肾脏代谢的影响,其通过调控肾脏多个下游基因表达,从而影响其代谢相关信号通路,试验结果为阐明ZBED6基因功能提供了材料。     

鸡胸腺组织中免疫相关circRNA鉴定

【目的】 筛选鸡胸腺中免疫相关的环状RNA （circular RNA，circRNA），研究circRNA在鸡免疫调节中的作用。【方法】 采用高通量测序平台Illumina PE150对黄羽肉鸡和山地乌骨鸡的胸腺组织进行转录组测序，用生物信息学分析软件find_circ和CIRI识别circRNA分子，鉴别其基本结构特征。通过差异显著性分析筛选差异表达circRNA，再经过GO功能和KEGG通路富集分析筛选出免疫相关circRNA，用IRESfinder和PFAM软件进行circRNA编码潜能预测，并用实时荧光定量PCR进行验证。【结果】 circRNA在鸡胸腺中有丰富的表达，共识别到8 015个circRNAs，其中115个circRNAs在胸腺组织中的表达量呈显著差异（log₂|FoldChange|>1，P<0.05），在黄羽肉鸡中上调表达基因54个，下调表达基因61个（以山地乌骨鸡为对照）。GO功能富集结果显示有104个GO条目与免疫相关，涉及23个差异表达cricRNAs，占比高达20%（23/115）。KEGG通路分析结果显示，差异表达circRNA来源基因集中在细胞因子-细胞因子受体相互作用、RIG-Ⅰ样受体信号通路、Toll样受体信号通路等7个免疫相关通路中，共涉及7个差异表达circRNAs。有5个circRNAs在GO功能和KEGG通路富集分析中重复出现：circ_0002478（IL-1R1）、circ_0003208（MAPK11）、circ_0001674（STX8）、circ_0005105（RIPK2）和circ_0003590（BRAF），推测它们可能在免疫系统调节中扮演着重要角色。实时荧光定量PCR结果显示，5个circRNAs在2个鸡品种胸腺组织中的表达差异与测序结果一致，circ_0002478（IL-1R1）不仅在两品种间差异极显著（P<0.01），而且具有编码多肽潜能，它可能以多种方式参与免疫调节。【结论】 鸡胸腺组织中有丰富的circRNA，在与免疫相关的circRNA中circ_0002478（IL-1R1）可能通过多种方式参与免疫调节，可作为鸡免疫力研究的重点关注靶点。     

鸡miR-10b-5p生物信息学及组织表达分析

【目的】 通过对苏禽3号肉鸡的miR-10b-5p进行靶基因预测、生物信息学分析及其在鸡不同生长时期组织表达变化规律研究,探索其在鸡肌肉生长发育中作用。【方法】 通过文献和miRBase检索脊椎动物的miR-10b-5p序列,利用Ensembl数据库查询并确定miR-10b-5p在基因组中的位置,根据前体序列构建系统进化树;使用miRDB、microT和TargetScan网站预测miR-10b-5p靶基因,并对靶基因进行GO功能和KEGG通路分析,进而揭示miR-10b-5p的基因功能。采用实时荧光定量PCR技术检测不同生长时期鸡miR-10b-5p组织表达量,探索其在大脑、小脑、垂体、下丘脑、胸肌、腿肌、卵巢、心脏、肝脏、脾脏、肺脏和肾脏组织中表达变化。【结果】 miRBase检索共获得59种动物59条miR-10b-5p序列,即每一种动物都只有一种miR-10b-5p形式。基因定位分析发现,鸡miR-10b-5p位于2号染色体上的HOX基因家族中。常见物种miR-10b-5p成熟序列比对分析结果表明,miR-10b-5p的成熟序列相似性较高,物种间较为保守。系统进化树分析发现,miR-10b在哺乳类中的灵长目、啮齿目、鸟类、鱼类中各自先聚为一支,然后再共同聚为一类,这表明miR-10b在进化过程中是保守的。靶基因预测发现,miR-10b-5p共有300个靶基因。GO功能分析发现,靶基因主要富集到染色质组装和基因表达的转录后调控、组蛋白甲基转移酶复合物、转录因子复合物等功能。KEGG通路分析表明,靶基因主要富集到Wnt和p53信号通路上。组织表达分析显示,miR-10b-5p在检测的各个组织中均有表达;3日龄时大脑、小脑、垂体、胸肌和腿肌中miR-10b-5p表达量显著高于心脏、下丘脑、肾脏和肺脏等组织(P<0.05);90日龄时,垂体、胸肌、腿肌和大脑中miR-10b-5p表达量显著高于其他组织(P<0.05)。与3日龄雏鸡相比,90日龄时垂体、胸肌和腿肌中miR-10b-5p表达量均显著上升(P<0.05)。【结论】 鸡miR-10b-5p是组织广泛表达的miRNA,miR-10b在进化过程中是保守的,其可能通过Wnt及p53信号通路调控肌肉细胞增殖分化,进而调控鸡肌肉生长发育。     

鸡繁殖性状近交衰退相关lncRNA的筛选

【目的】 筛选与鸡繁殖性状近交衰退相关的长链非编码RNA （long non-coding RNA，lncRNA），为深入探究lncRNA对鸡繁殖性能近交衰退调控机制提供参考。【方法】 以狼山鸡高、低近交群体为试验素材，通过转录组测序技术分析狼山鸡下丘脑和卵巢中lncRNA表达情况，筛选高、低近交组间差异表达lncRNA，并对其进行顺式调控靶基因预测及功能分析。【结果】 狼山鸡下丘脑和卵巢中共鉴定出4 222个lncRNAs，高、低近交组间比较发现，下丘脑中差异表达lncRNAs 35个，卵巢中差异表达lncRNAs 215个。下丘脑中差异lncRNAs中预测到顺式调控靶基因98个，这些靶基因显著富集于出生或孵化时胚胎发育终止、胚胎心导管发育、视黄酸应答等繁殖相关生物过程（P<0.05），涉及lncRNA MSTRG.9196.4、MSTRG.9195.2、MSTRG.6254.2以及相应靶基因DNAAF2和FKBP1B等胚胎发育相关基因；卵巢中差异lncRNAs预测到顺式调控靶基因414个，这些靶基因富集到卵母细胞减数分裂、MAPK、叶酸合成等信号通路，包括MSTRG.7683.1、MSTRG.13604.4、MSTRG.16570.1、MSTRG.8330.5、MSTRG.8330.4等神经内分泌调节及配子生成相关的lncRNAs。【结论】 本研究在下丘脑和卵巢中筛选到了一系列鸡胚胎发育及配子生成调控相关差异lncRNAs，这些lncRNAs可作为鸡繁殖性能近交衰退候选lncRNAs，为进一步揭示鸡繁殖性能近交衰退调控机制提供线索。     

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

【目的】 本研究旨在对云南半细毛羊精子冷冻前后差异表达的基因与精子冷冻损伤的关系进行研究。【方法】 采用电刺激法采集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个信号通路。其中与精子功能密切相关的主要包括精子细胞过程、代谢过程、应激反应、生殖、质膜和氧化还原反应;精子冷冻前后差异表达基因参与的信号通路主要有炎症、细胞因子-细胞因子受体相互作用和细胞凋亡等通路。【结论】 云南半细毛羊精子的冷冻损伤可能与获得的差异表达基因参与的生物过程有关。此外,获得的差异表达基因也可作为分子标记物应用于绵羊冻精品质的评价。