简州大耳羊肌内脂肪细胞成脂分化差异表达基因的筛选与鉴定

旨在通过转录组测序技术获得简州大耳羊肌内前体脂肪细胞成脂分化前后的差异表达基因,并经生物信息学分析获得相关信号通路及可能发挥作用的关键功能候选基因。本研究以7日龄的健康简州大耳羊公羊为试验动物（n=3）,采用胶原酶消化法分离获得其肌内前体脂肪细胞;利用Illumina平台对肌内前体脂肪细胞和诱导分化5 d的肌内脂肪细胞cDNA样品（n=3）进行高通量测序;以|log₂fold change|>0和P<0.05为阈值筛选获得差异表达基因,并利用clusterProfiler R包对其进行GO功能和KEGG通路富集;最后利用实时荧光定量PCR（quantitative real-time PCR,qRT-PCR）技术检测功能候选基因在细胞分化前后的表达量变化。结果显示,共获得差异表达基因7 916个,其中4 143个为表达上调基因,主要富集到氧化磷酸化、核糖体合成和三羧酸循环通路等305条通路;3 773个为表达下调基因,且主要富集到甲状腺激素信号通路等303条通路;差异基因GO功能注释中52.8%为生物过程、13.4%是细胞组成和33.8%是分子功能。qRT-PCR结果表明,UCP3、ACACB、ACOT11、ACOX3、APOA1和WISP2在分化前后的山羊肌内脂肪细胞中的表达趋势与RNA-seq结果一致,提示这些基因适合作为下一步研究的功能候选基因。本研究筛选得到简州大耳羊肌内脂肪细胞成脂分化的差异基因,并确认了6个基因可作为功能候选基因。研究结果为阐明肉用山羊肌内脂肪细胞成脂分化的分子调控网络提供基础数据和系统资料。     

FST和FMOD基因在绵羊胎儿妊娠中后期的表达特征比较

本研究以骨骼肌表型差异明显的两个绵羊品种特克塞尔羊和乌珠穆沁羊为试验对象,通过对两个与生长发育相关的基因卵泡抑素(follistatin,FST)与纤维蛋白聚糖基因(fibromodulin,FMOD)在不同日龄、不同骨骼肌部位表达的比较分析,探讨绵羊骨骼肌早期发育规律;同时对绵羊胎儿的85、100、120和135日龄4个阶段,在背最长肌、半腱肌、股四头肌、半膜肌和股二头肌5个部位肌肉FST和FMOD基因的表达量进行了比较研究。实时荧光定量PCR结果表明,随着胎儿日龄的增加,FST基因在半腱肌、股四头肌中表达量有增加的趋势,同时在特克塞尔绵羊中的表达量高于乌珠穆沁羊。在100和120日龄半腱肌、120日龄半膜肌、120日龄股二头肌中FST基因表达量在两个品种间差异显著(P<0.05),在100和135日龄半膜肌、85日龄背最长肌、100日龄股四头肌、100日龄股二头肌中FST基因表达量在两个品种间差异极显著(P<0.01)。FMOD基因在股四头肌中表达量较高,而在其他4种肌肉组织中相对表达量较低,且在乌珠穆沁羊中的表达量高于特克塞尔羊中的表达量。FST和FMOD基因在绵羊早期骨骼肌发育中起着重要作用。     

基于WGCNA与GSEA方法挖掘调控中卫山羊羊毛弯曲相关基因

旨在探究中卫山羊羊毛弯曲随生长发生变化的分子机制，挖掘不同发育时期影响羊毛弯曲度的关键基因。试验采集了3只中卫山羊出生后45日龄（弯曲毛）和365日龄（直毛）的皮肤组织，进行转录组测序（RNA-seq），使用WGCNA与GSEA分析找出Hub基因。以P-value<0.05和|log₂FoldChange|≥1作为差异表达基因筛选的标准，45日龄为对照组，365日龄为试验组共筛选得到1 252个显著差异表达基因，包括812个表达上调基因和440个下调基因。基于WGCNA方法分析共得到14个模块，其中黄色和绿松石模块与被毛弯曲表型相关。利用String和Cytoscape构建网络，筛选到20个影响羊毛弯曲度的核心基因。从GSEA结果中筛选的被显著富集通路Hedgehog signaling pathway、JAK/STAT signaling pathway、TGF/BETA signaling pathway等参与了毛囊发育调控。综合KEGG、WGCNA、分子网络构建、GSEA分析结果，共同筛选得到基因CCL27、IL7、WNT2，推断这3个基因在调控毛囊发育中起到了关键的调控作用。本研究为进一步阐明动物毛发弯曲的分子机制提供了理论基础。     

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

【目的】 通过分析马身猪和大白猪背最长肌转录组测序数据,筛选肌内脂肪沉积相关基因,并对其表达特性进行分析。【方法】 采集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个与猪肌内脂肪沉积相关的候选基因,可为后续进一步探讨猪肌内脂肪沉积调控机制提供一定的参考。     

藏羊MSTN基因克隆、生物信息学及组织表达分析

【目的】 对藏羊肌肉生长抑制素（myostatin，MSTN）基因进行克隆和生物信息学分析，检测其在藏羊不同组织中的表达，为探究MSTN基因在藏羊中的生物学功能提供参考。【方法】 以藏羊背最长肌组织cDNA为模板，克隆藏羊MSTN基因完整CDS区序列并测序，用SeqMan程序对测序结果进行拼接，并用BLAST在线程序对组装后的序列进行分析鉴定。用生物信息学软件进行相似性比对、系统进化树构建及生物信息学分析，用实时荧光定量PCR检测MSTN基因在藏羊不同组织中的表达量。【结果】 藏羊MSTN基因CDS全长为1 128 bp，编码375个氨基酸。藏羊MSTN基因氨基酸序列与绵羊、牦牛、牛、猪、恒河猴、人、黑猩猩、犬、鸡及斑马的相似性依次为100.0%、93.4%、93.4%、95.2%、94.4%、94.2%、94.4%、93.1%、87.8%和87.5%；系统进化树分析结果表明，藏羊与绵羊的亲缘关系最近，与斑马和鸡的亲缘关系最远。藏羊MSTN蛋白属于亲水性分泌蛋白，且具有不稳定性，不含跨膜结构，含1个信号肽，存在31个潜在的磷酸化位点、2个N-糖基化修饰位点，主要分布在线粒体和细胞质中；MSTN蛋白二级结构以无规卷曲为主，其次是α-螺旋、延伸链和β-转角；三级结构预测结果与二级结构一致。实时荧光定量PCR结果显示，MSTN基因在藏羊不同组织中均有表达，其中在臂三头肌和半腱肌的表达量显著高于肺脏、下丘脑、心脏、肝脏、十二指肠、瘤胃和肾脏（P<0.05）。【结论】 成功克隆了藏羊MSTN基因；该基因在藏羊肌肉中的表达量高于内脏组织。该结果为进一步研究MSTN基因对藏羊肌肉生长发育的调控机制奠定了基础。     

基于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组）与小尾寒羊（X组）皮下脂肪组织中的特异性表达基因,并研究其潜在的作用,为理解绵羊脂肪组织发育规律以及对脂代谢相关疾病的预防和治疗研究提供依据。本研究选取脂肪沉积能力存在差异、健康无病、体况良好、种内个体体重相近（约50 kg）的雌性成年多浪羊和小尾寒羊为试验材料,分为D组（试验组）和X组（对照组）,每组3个重复,采集位于背最长肌的皮下脂肪组织,应用RNA-Seq技术和生物信息学方法进行转录组测序并对结果进行分析。以|Fold change|≥2,P adjust≤0.05为标准筛选差异表达基因,通过对差异表达基因进行GO功能注释和KEGG通路富集分析,得到与脂肪沉积和脂代谢有关的差异基因。为了验证测序数据的可靠性,本研究随机选取了6个差异表达基因进行qRT-PCR验证。结果显示,在6个样本中共检测到38 672个已知的mRNAs,新的mRNAs为1 606个,在两组中共有839个差异表达基因,其中有320个差异基因在多浪羊组中上调表达,有519个差异基因在多浪羊组中下调表达。通过GO功能注释分析发现,差异表达基因主要参与脂质分解代谢过程、脂质生物合成过程、脂质分解代谢负调控过程、MAPK级联反应调控、对甘油三酯的反应等生物学过程。KEGG通路富集结果显示,差异表达基因显著富集到了PI3K-Akt、MAPK、胰岛素以及PPAR等信号通路中。qRT-PCR结果与测序结果一致,表明测序结果可靠。通过对多浪羊和小尾寒羊皮下脂肪组织进行转录组测序以及生物信息学分析,筛选到与脂肪沉积和脂代谢相关的差异表达基因,这些基因主要参与脂质生物合成、脂质代谢等过程,其中COL1A1、AKT2、SCD、LPL、PCK1与PPP2R5A可能在多浪羊与小尾寒羊的皮下脂肪组织的沉积与代谢中发挥重要作用。     

绵羊FGF10基因克隆、序列分析及其在毛囊发育中的表达分析

本试验旨在获得中国美利奴羊成纤维细胞生长因子10（fibroblast growth factor 10，FGF10）基因的编码区（CDS）全长序列并进行生物信息学分析，随后对FGF10基因在中国美利奴羊毛囊发育过程中的表达特征进行分析，明确其在中国美利奴羊毛囊发育过程中的表达模式，为进一步研究FGF10 mRNA表达水平与中国美利奴羊毛囊生长发育的表达调控机制奠定理论基础。采用PCR扩增获得中国美利奴羊FGF10基因CDS，并克隆到zero PCR@^TM-Blunt进行测序验证；利用实时荧光定量PCR技术检测FGF10在中国美利奴羊毛囊发育过程中的表达差异。结果表明，绵羊FGF10基因CDS长度为696 bp（序列上传GenBank，获得登录号：MT872422），编码231个氨基酸，与牛和山羊的氨基酸序列同源性达100%，存在1个信号肽和1个跨膜结构域，其为分泌通路信号蛋白；实时荧光定量PCR分析表明，FGF10基因在中国美利奴羊毛囊发育过程中均表达，在毛囊发育第85天表达最高，显著高于其他毛囊发育时期（P<0.05）。本研究获得中国美利奴羊FGF10基因完整的编码区序列和毛囊发育过程中的表达特征，生物信息学分析发现，FGF10基因编码区序列具有物种间的保守性，同时FGF10在绵羊毛囊不同发育阶段的皮肤组织中表达，由此表明，FGF10基因可能在绵羊毛囊的生长发育过程中发挥重要的生物学作用。     

简州大耳羊肌内脂肪细胞成脂分化差异表达基因的筛选与鉴定

旨在通过转录组测序技术获得简州大耳羊肌内前体脂肪细胞成脂分化前后的差异表达基因,并经生物信息学分析获得相关信号通路及可能发挥作用的关键功能候选基因。本研究以7日龄的健康简州大耳羊公羊为试验动物(n=3),采用胶原酶消化法分离获得其肌内前体脂肪细胞;利用Illumina平台对肌内前体脂肪细胞和诱导分化5 d的肌内脂肪细胞cDNA样品(n=3)进行高通量测序;以|log_2fold change|0和P0.05为阈值筛选获得差异表达基因,并利用clusterProfiler R包对其进行GO功能和KEGG通路富集;最后利用实时荧光定量PCR(quantitative real-time PCR, qRT-PCR)技术检测功能候选基因在细胞分化前后的表达量变化。结果显示,共获得差异表达基因7 916个,其中4 143个为表达上调基因,主要富集到氧化磷酸化、核糖体合成和三羧酸循环通路等305条通路;3 773个为表达下调基因,且主要富集到甲状腺激素信号通路等303条通路;差异基因GO功能注释中52.8%为生物过程、13.4%是细胞组成和33.8%是分子功能。qRT-PCR结果表明,UCP3、ACACB、ACOT11、ACOX3、APOA1和WISP2在分化前后的山羊肌内脂肪细胞中的表达趋势与RNA-seq结果一致,提示这些基因适合作为下一步研究的功能候选基因。本研究筛选得到简州大耳羊肌内脂肪细胞成脂分化的差异基因,并确认了6个基因可作为功能候选基因。研究结果为阐明肉用山羊肌内脂肪细胞成脂分化的分子调控网络提供基础数据和系统资料。     

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基因发挥作用后可以介导多个下游基因表达,从而影响相关信号通路及生物学过程;同时,所筛选出的差异表达基因可作为进一步研究骨骼肌调控机制的候选靶标。     

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

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

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

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

ZBED6基因对巴马香猪脾发育的作用机制分析

ZBED6基因作为一个转录因子,能够与IGF2结合进而调节肌肉的生长发育。但其在脾生长发育中的作用尚不清楚,本研究利用RNA-seq测序技术比较ZBED6基因敲除巴马香猪（ZBED6-KO）和同日龄野生型巴马香猪（WT）的脾组织转录组,探究ZBED6基因对巴马香猪脾组织的发育影响。利用t检验对ZBED6-KO猪和WT猪脾组织大小的表型差异及ZBED6直接调控的靶基因IGF2的表达量进行显著性分析。提取ZBED6-KO猪和WT猪脾组织的总RNA,在Illumina Hiseq 2500平台进行RNA-seq分析。以猪Sus scrofa11.1为参考基因组,用转录组分析的标准流程筛选ZBED6-KO猪和WT猪脾组织中的差异表达基因。用DAVID在线网站对差异表达基因进行GO和KEGG富集分析。然后,随机选取7个差异表达的基因,利用实时荧光定量PCR技术验证测序结果的准确性与可靠性。结果显示,与WT猪相比,ZBED6-KO猪脾的重量和IGF2基因表达量均显著增加（P<0.05）,表明ZBED6基因的敲除对巴马香猪脾组织的生长发育有一定的促进作用。测序结果显示,各样本至少获得4G的数据量,每个样本的Clean Ratio及Q30比率均在90%以上,83.94%以上的reads可比对在猪的参考基因组上,表明测序数据质量良好,真实可靠;对测序数据进行转录组分析,共筛选到161个差异表达基因,其中上调基因90个,下调基因71个;差异表达基因的层次聚类分析显示,ZBED6-KO猪的3个个体（spleen1、spleen3、spleen6）的表达模式相似,WT的3个个体（spleen2、spleen4、spleen5）的表达模式相似,进一步证明测序数据的准确可靠;GO和KEGG富集分析中,富集到10条显著的GO条目以及5条显著的信号通路,2条与肌肉发育相关的通路;实时荧光定量PCR试验随机检测7个差异表达基因的表达模式与RNA-seq分析结果相一致,证实了测序数据的可靠性。以巴马香猪为模型,利用RNA-seq技术研究ZBED6基因的敲除对中国地方猪脾发育的作用影响,为挖掘ZBED6基因的更多功能提供了基础。     

基于瘤胃转录组探究双峰驼沙漠适应性分子机制

旨在比较胚胎期和成年期双峰驼瘤胃在组织形态、编码基因表达上的变化,挖掘影响双峰驼瘤胃发育的关键基因和通路,从消化系统探究双峰驼的沙漠适应性机制。本研究选取3峰10~12岁健康状况良好的成年期阿拉善双峰驼,以及3峰9~10月龄健康状况良好的胚胎期阿拉善双峰驼,采集瘤胃组织样品,制作石蜡组织切片,观察成年期与胚胎期双峰驼瘤胃,比较其组织结构之间的差异。分别提取成年期和胚胎期双峰驼的瘤胃组织总RNA,利用Illumination Hiseq 2000测序平台分别进行转录组测序,对转录组数据进行质控、比对、差异基因筛选、GO、KEGG富集分析。随机选择6个差异表达基因进行荧光定量试验,关联转录组数据与荧光定量试验的表达趋势。结果表明,胚胎期瘤胃组织中上皮细胞和肌纤维细胞清晰可见,分布密集,在成年期的瘤胃组织中,可见明显的肌纤维,肌纤维直径较宽,肌纤维间的空隙较大。转录组测序结果显示,每个样本获得至少10G的数据量,各样本的质控率都在90%以上,Q30数据都在88%以上。对测序数据进行分析,以胚胎期为对照组,成年期为试验组,筛选到1 207个差异表达基因,其中上调基因456个,下调基因751个;对差异表达基因进行层次聚类分析,结果显示,成年期的3个个体（M1、M2、M3）表达模式接近,胚胎期的3个个体（T1、T2、T3）表达模式接近。对差异表达基因进行GO和KEGG富集分析,结果显示,上调差异表达基因显著富集到62条显著的GO条目,下调差异表达基因显著富集到366条显著的GO条目,73条显著的KEGG通路。差异表达基因主要富集在代谢过程的负调控、RNA生物合成过程的负调控、基因表达的负调控等GO条目中,KEGG主要富集到MAPK信号通路、PI3K-Akt信号通路、糖尿病并发症中的AGE-RAGE信号、胰岛素信号通路、醛固酮的合成与分泌等通路中。同时筛选到MAPK12、MAPK13、FABP5、PPARγ、CaMK1等与双峰驼沙漠适应性相关的基因。荧光定量结果显示,差异基因在成年期和胚胎期瘤胃组织中的表达模式与RNA-Seq的结果一致。上述结果表明,双峰驼的瘤胃组织可能与脂肪细胞分化有关。在沙漠环境中,双峰驼可能降低瘤胃组织代谢效率、通过胰岛素抵抗作用提高血糖耐受性以及促进醛固酮的合成与分泌以调节血压平衡,以更好地在沙漠干燥缺水的环境中生存。     

系统分析多组织转录组鉴定影响猪脂肪沉积的关键基因

旨在挖掘影响松辽黑猪脂肪沉积的关键基因及脂肪、肝和肌肉在体内的功能。本研究选择体重100 kg左右健康且背膘厚差异显著的6头(高、低各3头)松辽黑猪为试验动物,利用高通量转录组测序技术检测其脂肪、肝和背最长肌组织中基因的表达水平,鉴定不同脂肪沉积猪和不同组织中的差异表达基因,并分析差异表达基因的生物学功能。结果表明,在不同分组的猪中发现135个差异表达基因,其中部分参与了PPAR信号通路、AMPK信号通路、代谢通路、脂肪酸代谢和甘油代谢等通路。经生物学功能分析发现,EHHADH、ME1、SCD、OLR1、PHGDH、ACLY、LEP和CYP超家族基因等基因为影响猪脂肪沉积的关键基因。在不同组织的差异表达表达基因中,脂肪组织中高表达的基因显著富集在胰岛素信号通路、MAPK信号通路、三羧酸循环、氧化磷酸化等通路;肝中高表达基因显著富集在多种物质的代谢、脂肪酸的降解、氨基酸的合成等通路;背最长肌中高表达的基因主要参与了蛋白质的降解、PI3K-Akt信号通路、氧化磷酸化通路、Wnt信号通路、磷脂酰肌醇信号通路等通路。不同组间差异表达基因分析结果提示,EHHADH、ME1、SCD、OLR1、PHGDH、ACLY、LEP和CYP超家族基因等基因是影响脂肪沉积的关键基因;不同组织间差异表达基因表明,脂肪组织是脂肪合成的主要部位,而肝和肌肉组织主要涉及脂肪酸的降解。本研究结果对脂肪性状的遗传改良、机理解析有一定意义。     

Cutaneous smooth muscle tumors in the dog and cat

Cutaneous smooth muscle tumors may arise from arrector pili muscles and from smooth muscles of the dermal vasculature. This report describes histologic and immunohistochemical features of eight arrector pili hamartomas in 8 dogs, 15 piloleiomyomas in 10 dogs and 3 cats, 10 piloleiomyosarcomas in 9 dogs and 1 cat, 1 angioleiomyoma in 1 cat, and 9 angioleiomyosarcomas in 6 dogs and 3 cats. Hamartomas and tumors arising from arrector pili muscles preferentially originated from the dorsal trunk. 5/5 (100%) arrector pili hamartomas, 10/12 (83%) piloleiomyomas, 4/5 (80%) piloleiomyosarcomas, 1/1 (100%) angioleiomyoma, and 6/7 (86%) angioleiomyosarcomas were positive for smooth muscle actin. 5/5 (100%) arrector pili hamartomas, 10/12 (83%) piloleiomyomas, 4/5 (80%) piloleiomyosarcomas, 1/1 (100%) angioleiomyoma, and 1/7 (14%) angioleiomyosarcomas were positive for desmin. Two incompletely excised canine angioleiomyosarcomas recurred locally. Metastases were not reported.     

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

试验旨在分析鸭瘟病毒（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的致病机制提供了理论参考。     

金黄色葡萄球菌在生物被膜态与浮游态的转录组差异表达分析

金黄色葡萄球菌是引起奶牛细菌性乳腺炎的主要原因，生物被膜的形成是金黄色葡萄球菌在不利环境条件下持久性存在的关键因素。探索同一株菌在生物被膜态与浮游态生长状态下的耐药性与其生长状态的相关性，可为进一步探究金黄色葡萄球菌的耐药性机制奠定基础。本研究培养了金黄色葡萄球菌的生物被膜，使用光学显微镜和扫描电镜观察其形成过程。测定并比较了9种抗菌药物对32株金黄色葡萄球菌在生物被膜态和浮游态的最小抑菌浓度，并对两种状态下的金黄色葡萄球菌进行转录组学测序，筛选出具有显著性差异的细胞信号通路和表达基因，同时对主要差异表达的基因进行RT-qPCR验证。结果发现，在生物被膜形成前期，随着培养时间的延长，显微镜下观察到的生物被膜态菌聚集面积越来越大，结构也越来越紧密，培养至72 h后，生物被膜逐渐开始分散。MIC测定结果显示浮游态菌的抑菌浓度低于生物被膜态菌。转录组结果显示两种状态菌的差异表达基因共1 512个，其中，生物被膜态菌中上调基因760个，下调752个。GO与KEGG富集分析显示，相比于浮游态菌，生物被膜态菌中与代谢相关的通路显著富集，其次为氨基酸的生物合成和ABC转运蛋白通路。与生物被膜形成相关的基因，如编码ABC转运蛋白的基因表达上调，而与代谢途径相关的基因下调。RT-qPCR验证了10个主要差异基因，其表达差异趋势与转录组测序结果一致。这些差异可能对金黄色葡萄球菌生物被膜态的高耐药性和细菌毒力的研究有所帮助。