类志贺邻单胞菌感染杂交鲟肠道组织转录组分析

为探究病原菌感染杂交鲟肠道转录组变化情况,本试验以常见病原菌类志贺邻单胞菌接种约360 d杂交鲟,于接种24 h后,采集杂交鲟肠道组织样本,提取组织总RNA,采用Illumina HiSeqTM 2000进行转录组测序,筛选杂交鲟肠道差异表达免疫应答基因并进行GO功能分类和KEGG信号通路分析,结果显示：与正常对照组比较,试验感染组杂交鲟肠道差异表达基因为13 542个,其中上调基因为9 774个,下调基因为3 768个;GO分析发现,杂交鲟肠道差异表达基因显著性富集到生物学过程的主要有免疫系统过程、免疫效应过程、刺激反应等;显著性富集到分子功能的主要有DNA结合、信号受体活性、核酸转录因子活性等;显著性富集到细胞组分的主要是胞外区、胞外区组成部分、细胞膜等。KEGG分析发现,杂交鲟肠道差异表达基因参与的免疫信号通路有RIG-Ⅰ样受体信号通路、Toll样受体信号通路和细胞溶质DNA传感途径。这些研究结果为深入探究杂交鲟对肠道病原微生物感染的防御分子机制奠定了良好的基础。     

基于RNA-Seq筛选ALV-J感染汶上芦花鸡肝差异表达致病相关基因

为挖掘感染ALV-J汶上芦花鸡的肝差异表达致病相关基因。本研究以汶上芦花鸡为试验素材,分别在42和300日龄进行ALV血液病毒分离检测筛选ALV阴性和阳性个体,对ALV阳性个体有病理变化的肝组织和阴性个体肝组织进行PCR检测,分别选取3只ALV阳性（G1组）和阴性（G2组）个体的肝组织,并对G1组PCR产物测序、聚类分析确定ALV亚型,利用RNA-Seq测序技术筛选差异表达基因,并进行功能分析,利用荧光定量qRT-PCR对部分差异表达基因进行验证。结果表明,G1组样品经PCR检测、PCR扩增产物测序和聚类分析确定ALV为J亚型,转录组分析发现,共有42个差异基因在GO和KEGG中富集,其中,上调基因18个,下调基因24个。随机选取的5个差异表达基因qRT-PCR验证结果与RNA-Seq测序结果相一致。GO分析显示,差异表达基因主要涉及细胞过程、代谢过程、对刺激的反应、免疫系统等;KEGG分析显示,差异表达基因主要富集在细胞过程、信号传导、疾病、新陈代谢等信号通路。本研究通过转录组分析发现了影响ALV-J致病性的多个基因和关键信号通路,为深入了解ALV-J致病的分子机制提供了思路。     

C型产气荚膜梭菌外毒素致小鼠肠道损伤的转录组分析

旨在探索产气荚膜梭菌外毒素造成机体炎性损伤及免疫调控紊乱的毒性机制,为产气荚膜梭菌病的致病机理研究提供理论基础。将C型产气荚膜梭菌强毒株C59-2培养上清腹腔注射BALB/c小鼠,采集小肠样本进行转录组测序,筛选差异表达基因,并对其进行GO功能注释和KEGG通路富集分析。结果显示,共获得40.99 Gb有效碱基,筛选后共得到795个差异表达基因,其中,229个基因表达上调,566个基因表达下调,对随机选取的10个基因进行q-PCR验证,其相对表达量与转录表达谱一致。GO功能注释主要涉及G蛋白偶联核苷酸受体活性和G蛋白偶联嘌呤核苷酸受体活性等。KEGG通路富集分析发现,主要富集在TNF信号通路、IL-17信号通路、p53信号通路、FOXO信号通路、Toll样受体信号通路、NF-κB信号通路等。产气荚膜梭菌外泌毒素侵入机体,会激活TNF等炎性信号通路,进而造成肠道发生炎性损伤甚至坏死。     

鸭瘟病毒感染鸭十二指肠黏膜炎症相关基因及调控通路的筛选

试验旨在分析鸭瘟病毒（duck plague virus，DPV）感染樱桃谷鸭后十二指肠黏膜转录组水平的变化，筛选出DPV感染后参与炎症的相关基因和调控通路。DEV-GZ株经腿部肌肉接种35日龄鸭后，于接种后24、48和72 h采集鸭十二指肠黏膜组织样本，提取总RNA进行转录组测序，对数据进行质控与注释，筛选与炎症相关的差异表达基因及调控通路，并应用实时荧光定量PCR对部分差异表达基因进行验证。结果显示，对照组与试验组总计得到21.39 Gb的Clean Bases，测序样本有效序列占原始序列的99%以上，映射率高于83%。鸭接毒DPV后24 h十二指肠黏膜的差异表达基因有221个，参与炎症相关生物学过程的有3个，均为下调基因；接毒DPV后48 h差异表达基因有499个，参与炎症相关生物学过程的有17个，均为下调基因；接毒DPV后72 h差异表达基因有798个，参与炎症相关生物学过程的有20个，其中上调基因13个、下调基因7个。GO数据库分析显示，参与炎症相关的差异表达基因主要是CCR8、CCL19、CCL4、IL-17、IRF1、IFN-γ、SLC11A1等，涉及急性炎症反应、趋化因子受体活性、急性期反应、炎症反应和炎症反应的正调节等生物学过程。KEGG数据库分析显示，差异表达炎症相关基因主要富集在Toll样受体信号通路、NF-κB信号通路、IL-17信号通路、TNF信号通路和炎症性肠病等；实时荧光定量PCR结果显示，IFN-γ和MALT基因的相对表达倍数与转录组结果基本一致。本试验完成了DPV感染樱桃谷鸭十二指肠的转录组测序分析，获取了差异表达基因的功能注释信息，初步揭示了参与DPV感染的炎症相关通路，为深入探究DPV的致病机制提供了理论参考。     

鸡白痢沙门菌感染雏鸡的骨髓miRNA表达谱分析

旨在探究鸡白痢沙门菌(Salmonella enterica serovar Pullorum,S.Pullorum)感染雏鸡后骨髓miRNA的差异表达特征,为深入了解鸡白痢沙门菌的致病机制提供理论基础。将7日龄SPF雏鸡随机分为两组,分别口服鸡白痢沙门菌和PBS,于感染24 h后采集骨髓进行miRNA测序,筛选差异倍数≥2且P值≤0.05的差异表达miRNAs进行靶基因预测以及GO、KEGG富集分析,随机选取6个miRNAs进行qRT-PCR验证,构建与免疫过程相关的miRNA-mRNA靶点网络。通过miRNA测序,共获得20个已知的差异表达miRNAs,其中11个上调,9个下调。qRT-PCR结果表明,miRNA变化趋势与测序结果一致。GO分析结果表明,差异表达基因主要富集在膜运输、信号转导、免疫系统、碳水化合物代谢、糖类的生物合成和代谢等,KEGG的信号通路主要富集在Notch信号通路、Hedgehog信号通路、PPAR信号通路、AMPK信号通路、Hippo信号通路等,miRNA-mRNA网络互作发现,gga-miR-1466和gga-miR-6643-5p可能是参与免疫过程相关的关键miRNA。本研究解析了鸡白痢沙门菌感染雏鸡的骨髓miRNA表达谱特征,为了解鸡白痢沙门菌和鸡之间相互作用的复杂分子致病机制提供了新的见解,为防控鸡白痢沙门菌感染提供了新策略。     

非洲猪瘟病毒感染后猪外周血淋巴细胞Toll样受体信号通路的lncRNA-mRNA分析

【目的】丰富非洲猪瘟病毒（African swine fever virus,ASFV）感染后猪外周血淋巴细胞长链非编码RNA（long non-coding RNA,lncRNA）表达谱,并进一步挖掘影响Toll样受体信号通路的调控网络。【方法】试验动物感染ASFV,于第7天采集外周血并分离得到外周血淋巴细胞,运用Illumina高通量组学测序对外周血淋巴细胞中lncRNA进行测序,原始数据经处理后筛选获得差异表达的lncRNA,并进行靶基因预测,利用生物信息学方法对靶基因进行GO功能和KEGG信号通路富集分析,初步绘制与Toll样受体信号通路相关的lncRNA-mRNA调控网络,并对lncRNA-ENSSSCG00000041959在内的4个lncRNAs进行实时荧光定量RT-PCR验证。【结果】共筛选到73个差异表达的lncRNAs,其中上调表达lncRNAs 38个,下调表达lncRNAs 35个。GO功能分析结果显示,靶基因显著富集在调节免疫系统过程、防御反应、生物刺激、病毒反应和先天性免疫;KEGG信号通路富集分析显示,大部分靶基因与细胞循环、疾病及免疫应答有关,其中免疫相关信号通路有Toll样受体信号通路、TNF信号通路、产生IgA的肠道免疫网络等。进一步挖掘出lncRNA-ENSSSCG00000041959-RIPK1和lncRNA-ENSSSCG00000041959-IRAK1可能是影响Toll样受体信号通路的重要调控网络,实时荧光定量RT-PCR与测序结果一致。【结论】本研究初步鉴定出lncRNA-ENSSSCG00000041959-RIPK1和lncRNA-ENSSSCG00000041959-IRAK1可能是影响Toll样受体信号通路的lncRNA-mRNA调控网络,为进一步探索lncRNA调控ASFV感染机体免疫反应奠定了理论基础。     

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

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

基于转录组测序分析研究枸芪多糖对育肥猪肠黏膜免疫功能的作用机理

本试验旨在研究枸芪多糖对育肥猪肠黏膜免疫调控的关键基因,阐明其对肠黏膜免疫功能的可能作用机理。选用80~90日龄的健康育肥母猪180头,随机分为2组,每组6个重复,每个重复15头猪。对照组饲喂基础饲粮,试验组在基础饲粮中添加0.1%的枸芪多糖,试验期90 d。采集育肥猪的空肠肠段,提取RNA后进行转录组测序分析差异免疫基因的表达。结果显示:1)与对照组相比,饲粮中添加枸芪多糖差异表达的免疫基因有75个,其中有26个上调基因和49个下调基因。2)基因本体(GO)富集分析显示差异表达基因主要参与免疫调节、免疫应答和应激反应,信号通路分析显示差异表达基因主要富集在以白细胞介素-17(IL-17)信号通路、NOD样受体信号通路为主的免疫信号通路。3)通过对肠黏膜差异表达的关键基因原癌基因Jun(JUN)和原癌基因Fos(FOS)的探究可见,与对照组相比,试验组白细胞介素-18(IL-18)、补体5(C5)、JUN和FOS的基因表达受到抑制,Jun和Fos的蛋白表达显著降低(P<0.05)。综上所述,在饲粮中添加枸芪多糖可提高肥育猪的肠黏膜免疫功能。枸芪多糖通过调节JUN、FOS、IL-18等免疫相关基因的表达,影响IL-17和NOD样受体信号通路,以此来抑制炎症反应。     

鸭肠炎病毒感染鸭肝脏组织miRNA表达谱差异分析

为探究鸭感染鸭肠炎病毒（Duck enteritis virus,DEV）后肝脏miRNA表达谱的差异,试验以DEV-GZ株经腿部肌肉接种30日龄麻鸭,于感染后66、90和114 h采集鸭肝脏组织样本,提取组织总RNA,经质检合格后,采用高通量测序技术对对照组和试验组样品进行miRNA测序,筛选出DEV感染鸭肝脏组织的差异表达miRNA,对其进行生物信息学GO功能分类和KEGG信号通路分析,并随机选取部分差异表达miRNA进行实时荧光定量PCR验证。结果显示,鸭感染DEV后66、90和114 h,肝脏组织差异表达miRNAs数量分别为227、225和231个。GO功能注释显示,感染鸭肝脏差异表达miRNA在生物过程分类中主要为细胞过程、单有机体过程和代谢过程类别;在细胞成分分类中主要是细胞、细胞部分和细胞器类别;在分子功能分类中主要是绑定分子功能和催化活性功能类别。KEGG通路富集显示,差异表达miRNA主要涉及PI3K-Akt、JAK-STAT、磷脂酰肌醇信号通路系统、ECM-受体相互作用、MAPK、Wnt、Toll样受体、IL-17、脂质代谢、钙离子信号通路和cAMP等信号通路,其中感染66 h差异表达miRNA主要在生物系统及神经系统中发挥作用;感染90 h主要在内分泌系统及消化系统中发挥作用;感染114 h主要在全身生物、免疫和消化系统等中发挥作用。选取10个差异表达miRNAs进行实时荧光定量PCR验证,结果与高通量测序结果一致。表明DEV感染对鸭肝脏组织miRNA表达具有显著影响,为从宿主miRNA角度揭示DEV致病机制提供了参考依据。     

新型鹅细小病毒感染鸭的肝、胸腺、回肠转录组差异表达分析

旨在探究新型鹅细小病毒(novel goose parvovirus,NGPV)与宿主细胞相互作用的分子致病机制.本研究采用转录组测序技术对感染NGPV SD株14 d雏鸭的肝、胸腺和回肠进行分析.结果 显示:感染组回肠上调基因有78个,下调基因有31个;感染组胸腺上调基因有111个,下调基因有78个;感染组肝中上调基...     

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

【目的】通过对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个与生长发育相关的关键基因,为后续进一步探讨广西麻鸡生长发育的分子调控机制提供参考。     

猪A型塞内卡病毒感染诱导猪肾细胞(PK-15)长链非编码RNA差异表达谱的分析

为研究猪A型塞内卡病毒(SVA)感染猪肾细胞(PK-15)后对于PK-15 lncRNA表达谱的影响,本研究进行了高通量测序以检测SVA感染诱导PK-15产生的差异表达lncRNA。结果显示,在SVA感染后24 h,共诱导1043个lncRNA差异表达,其中已知的lncRNA 420个,未知的lncRNA 623个。GO富集和KEGG通路分析显示,lncRNA靶基因富集于p53信号通路、RIG-I样受体信号通路、MAPK信号通路、IgA生成的肠道免疫网络通路等。根据高通量测序结果选择了8个差异表达的lncRNA采用荧光定量PCR验证,验证结果与高通量测序结果一致。本研究初步表明lncRNA参与了SVA对PK-15细胞的感染过程,为SVA的分子机制研究提供了新的方向。     

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

试验旨在通过对藏鸡和白羽肉鸡垂体转录组测序和生物信息学分析,筛选出与鸡生长发育相关的差异表达基因（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结果显示,转录组测序结果准确可靠。综上研究,本研究初步揭示了影响藏鸡和白羽肉鸡生长发育速度差异的关键候选基因和信号通路,为了解鸡垂体组织调节生长发育的分子机制提供了理论依据。     

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

【目的】分析产肠毒素大肠杆菌（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感染诱导炎症反应的分子机制及相关调控网络提供参考依据。     

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

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

基于转录组数据分析环形泰勒虫对诱导宿主细胞转化相关基因的影响

本试验主要通过研究布帕伐醌(BW720)处理前后,环形泰勒虫转化牛单核细胞(TaNM Ⅰ)转录组数据的变化情况,为进一步研究环形泰勒虫诱导宿主细胞转化的分子机制奠定基础。以环形泰勒虫感染的牛单核细胞为试验材料,设置对照组(二甲基亚砜,DMSO)和试验组(BW720),利用Illumina HiSeq4000测序平台进行高通量测序,将测序得到的原始数据(Raw reads)与GenBank和Rfam数据库进行比对过滤,通过质控获得最终数据(Clean read)。利用Venn分析软件筛选出DMSO组和BW720组中差异表达明显的基因,进行功能注释(GO)以及信号通路(KEGG)分析。随机选择10个基因,利用荧光定量PCR (qPCR)检测基因在不同处理条件下细胞中的表达量。在验证测序结果的准确性后,选出差异表达明显的基因,分析其在TaNM Ⅰ细胞的凋亡、增殖等信号通路方面的作用。结果显示：BW720组和DMSO组中分别得到6 854 019 704和6 627 265 854条Raw reads,过滤、质控后分别得到22 925 606和22 171 427条Clean reads。BW720组和DMSO组中筛选出差异表达明显的基因共计4 054个(P＜0.05),其中,2 146个基因显著上调,1 908个基因显著下调;进一步筛选出共同差异表达的基因367个,其中,显著上调的196个,显著下调的171个。同时,qPCR验证随机选择10个基因的表达趋势与转录组测序结果一致,表明测序结果可靠。然后对差异表达基因进行GO功能注释,筛选到与细胞增殖相关的20个生物学过程、15个细胞组分以及11个分子功能条目。KEGG富集分析结果显示,在Top20的信号通路中筛选出与环形泰勒虫转化细胞有关的一些信号通路,如：癌症信号通路、PI3K-Akt信号通路、MAPK信号通路等。对这些信号通路进一步分析发现：PI3KR3、FOXO1、IL23A、FZD3、AKT、MMP9等基因在环形泰勒虫转化细胞的研究中也有相关文献的报道。本研究表明,在TaNM Ⅰ细胞中DAPK1、FZD3、FOXO3、PI3KR3等可能在感染细胞无限增殖的过程中发挥作用,为后续研究TaNM Ⅰ细胞无限增殖机制奠定了基础。     

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

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

基于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,为深入研究苏姜猪肉质性状的分子机理奠定基础。     

羊源D型多杀性巴氏杆菌感染小鼠脾脏转录组分析

试验旨在探索羊源D型多杀性巴氏杆菌入侵宿主脾脏组织中引发的免疫应答途径。首先利用羊源D型多杀性巴氏杆菌(HN01菌株)感染小鼠,建立羊源D型多杀性巴氏杆菌感染的动物模型;之后利用转录组测序技术获得感染小鼠与正常小鼠的脾脏转录组数据,并使用COG、KOG、eggNOG、GO、KEGG数据库对测序结果中的差异表达基因(differentially expressed genes,DEGs)进行功能注释与分析,同时对于显著富集到关键免疫通路的差异表达基因使用STRING软件和KEGG mapper进行蛋白互作分析,筛选出核心通路中起关键作用的基因;最后选取关键的10个基因进行实时荧光定量RT-PCR验证。转录组分析结果显示,与正常组相比,感染组中筛选出3 380个差异表达基因(P<0.01,log₂|FoldChange|≥0.5),其中1 691个基因上调,1 689个下调。基因功能富集分析结果表明,感染组脾脏中的差异表达基因主要发挥信号转导的功能,其主要参与的生物途径包括细胞因子与细胞因子受体互作通路、趋化因子信号通路、HIF-1信号通路、TNF信号通路。蛋白互作分析筛选出约28个核心差异表达基因,结合实时荧光定量RT-PCR验证后,其中9个基因的表达结果与测序一致,分别是C3、Cd4、Cxcl13、Lck、Gnai1、Grap2、IL-6、Cxcr6及Serping1基因。本研究初步证明了脾脏在抵抗多杀性巴氏杆菌入侵中参与了一系列免疫应答反应,为进一步研究羊源D型多杀性巴氏杆菌与宿主之间相互作用的分子机制奠定理论基础。