金茅黑鸡miR-206表达及靶基因预测分析

为了研究miR-206的组织表达分布及其在鸡骨骼肌生长发育中的表达规律,预测其可能的作用机制,本研究采集4周龄金茅黑鸡的心脏、肝脏、脾脏、肺脏、肾脏、胸肌、腿肌和腹脂8种组织及2、6、10、14、16周龄胸肌和腿肌组织,利用实时荧光定量PCR技术检测miR-206在8种组织及不同周龄胸肌和腿肌中的表达,利用生物信息学方法预测其靶基因并进行功能注释。结果显示,miR-206在公、母鸡胸肌和腿肌中表达量均极显著高于其他组织(P0.01),在腹脂、心脏、肝脏、脾脏、肺脏和肾脏中低表达,为鸡骨骼肌特异性表达miRNA;miR-206均在公、母鸡腿肌和胸肌组织生长早期(2周龄)表达最高,之后表达量逐渐下降。靶基因预测和功能分析发现,miR-206共有356个靶基因,其中21个与肌肉生成相关,多个靶基因已知参与调控肌肉生成,包括配对框7(paired box 7,Pax7)、胰岛素样生长因子1(insulin-like growth factor 1,IGF1)、脑衍生神经营养因子1(brain-derived neurotrophic factor 1,BDNF1)、视黄酸受体(retinoic acid receptor beta,RARB)和卷曲类受体7(frizzled class receptor 7,FZD7)基因。此外,发现miR-206靶基因富集到多条已知参与调控骨骼肌生成的信号通路,包括MAPK、Wnt、肌动蛋白骨架调控和黏着斑激酶信号通路。结合表达数据推测,miR-206可能通过靶向调控Pax7、IGF1、BDNF1、RARB、FZD7基因和MAPK、Wnt等信号通路参与调控鸡肌肉生成。本研究揭示了金茅黑鸡miR-206的组织表达分布及其在鸡骨骼肌生长过程中的表达规律,并预测了其调控骨骼肌生成的可能作用机理,为进一步研究miR-206在鸡肌肉生长发育中的作用机理和表达调控奠定了基础。     

鸡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信号通路调控肌肉细胞增殖分化,进而调控鸡肌肉生长发育。     

新西兰白兔MSTN基因克隆、生物信息学分析及表达谱构建

肌肉生长抑制素（MSTN）是骨骼肌发育的主要调控因子，为获得新西兰白兔MSTN基因序列及其表达规律，试验通过RT-PCR技术从新西兰白兔腿肌组织中扩增并克隆得到MSTN基因序列，利用在线网站对其进行生物信息学分析，并采用实时荧光定量PCR技术检测MSTN基因在新西兰白兔同一时期不同组织及不同时期腿肌组织中的表达量。结果显示，新西兰白兔MSTN基因CDS区序列长1 128 bp，编码375个氨基酸，其核苷酸序列与GenBank上公布的人、猪、马、小鼠、犬、牛和鸡的核苷酸序列相似性分别为95.9%、94.9%、94.9%、91.7%、91.5%、91.3%及84.2%。新西兰白兔MSTN蛋白属于酸性、亲水性分泌蛋白，且具有不稳定性，不含跨膜结构，含有1个信号肽，主要分布在内质网和线粒体中。MSTN蛋白二级结构和三级结构预测结果一致，以无规则卷曲为主，其次是α-螺旋、延伸链，为混合型蛋白。蛋白质互作关系预测发现，MSTN蛋白与PAX7、MYOG、IGF1、FST、Smad3及Smad2等与肌肉生长发育有关的蛋白之间存在相互作用。不同组织表达分析结果表明，新西兰白兔MSTN基因在心脏、肝脏、脾脏、肺脏、肾脏、腿肌、子宫和胃中均有表达，在腿肌中表达量最高，极显著高于其他组织（P<0.01），其次是肾脏，在肝脏中的表达量最低。不同时期腿肌组织中的表达量分析结果显示，新西兰白兔MSTN基因在180日龄腿肌组织中的表达量显著高于90和240日龄（P<0.05）。研究结果为深入探讨MSTN基因对家兔肌肉生长发育的调控机制奠定了基础。     

m6A甲基化在鸡肌肉生长发育中的表达研究

试验旨在研究RNA m6A修饰相关基因去甲基化酶Alk B同源蛋白5（Alk B homologue 5，ALKBH5）、去甲基化酶肥胖相关蛋白（fat mass and obesity-associated protein，FTO）、甲基转移酶样蛋白3（methyltransferase like 3，METTL3）、甲基转移酶样蛋白14（methyltransferase like 14，METTL14）和成肾细胞瘤1-结合蛋白（Wilms’tumor 1-associating protein，WTAP）在鸡骨骼肌发育过程中的表达，分析其与骨骼肌m6A甲基化水平的相关性。首先，利用实时荧光定量PCR技术检测m6A甲基化相关基因在金茅花鸡12（E12）、14（E14）、16（E16）、18（E18）胚龄和1日龄腿肌和胸肌组织中mRNA表达水平，以及其在鸡成肌细胞50%、100%增殖期和1、2、3、4、5 d分化期的mRNA表达水平；随后，利用m6A甲基化试剂盒检测金茅花鸡E12和1日龄腿肌和胸肌组织中m6A甲基化修饰水平，与m6A甲基化相关基因表达水平进行相关性分析。结果显示，m6A去甲基化基因ALKBH5和FTO mRNA表达水平在骨骼肌发育过程中显著上调（P<0.05），即在E12、E14低表达，E16、E18逐渐上调，1日龄达到最高。m6A甲基化写入基因METTL14、METTL3和WTAP mRNA表达水平在E12、E14、E16逐渐上升，E18下降，随后至1日龄表达量回升。在细胞增殖过程中，ALKBH5、FTO、METTL14、METTL3和WTAP基因表达均上调；在细胞分化过程中ALKBH5和FTO基因表达水平显著上调（P<0.05），在分化第5天达到最高。METTL14、METTL3和WTAP基因mRNA表达水平在细胞诱导分化的1、2、3、4 d表达量呈下降趋势，而在诱导分化的第5天有所回升。甲基化水平检测结果显示，腿肌和胸肌m6A甲基化水平变化趋势一致，均在胚胎发育过程中显著下降（P<0.05），至1日龄达到最低。相关性分析结果显示，鸡骨骼肌RNA m6A甲基化水平与m6A去甲基化修饰基因ALKBH5、FTO mRNA表达水平呈显著负相关（P<0.05）。综合以上试验结果，推测m6A甲基化修饰与鸡骨骼肌发育相关，而去甲基化基因ALKBH5、FTO可能通过调控RNA m6A甲基化水平，影响鸡骨骼肌发育。本研究结果为进一步研究m6A甲基化修饰调控鸡骨骼肌生长发育的功能和分子机制提供理论依据。     

PRELID3B基因mRNA在不同品种乌骨鸡肌肉组织中的表达规律研究

PRELID3B基因（PRELIP Domain Containing 3B）是调控乌骨鸡黑色素沉积的重要候选基因。为探讨PRELID3BmRNA在不同品种乌骨鸡肌肉发育中的表达情况，以丝羽乌骨鸡和余干乌骨鸡为材料，采用实时荧光定量PCR技术，检测出雏后不同发育时期胸肌和腿肌中PRELID3BmRNA的表达情况，并将基因表达量与肌肉色度L值进行相关性分析。结果发现，2个品种胸肌和腿肌的肌肉色度L值都随日龄的增长呈逐渐下降的趋势；2周龄后，2个品种的胸肌和腿肌L值出现差异（P<0.05）。丝羽乌骨鸡胸肌和腿肌PRELID3B mRNA表达变化趋势基本一致，都是前高后低；余干乌骨鸡胸肌和腿肌PRELID3BmRNA表达变化趋势不一致，其中胸肌在10周龄之后呈升高的趋势。14周龄之后，2个品种的胸肌和腿肌PRELID3BmRNA表达都出现差异（P<0.05）。品种间比较，14周龄之后胸肌PRELID3BmRNA表达在品种间出现差异（P<0.05）；腿肌则是在2周龄、6周龄和14周龄时品种间出现差异（P<0.05）。肌肉色度L值与基因表达的相关性结果显示，2个品种的腿肌...     

QDPR基因克隆及其在乌蒙凤鸡不同生长阶段组织表达研究

旨在克隆乌蒙凤鸡醌型二氢生物喋呤还原酶（quinoid dihydropteridine reductase，QDPR）基因，并检测其在不同性别乌蒙凤鸡不同生长阶段各组织中的表达差异，以研究QDPR基因在乌蒙凤鸡生长发育过程中的调控作用。本试验对乌蒙凤鸡QDPR基因CDS区进行PCR扩增和克隆，并对得到的序列进行生物信息学分析，通过实时荧光定量PCR检测QDPR基因在4、12、20和28周龄乌蒙凤鸡公鸡及母鸡的心、肝、脾、肺、肾、胸肌和腿肌组织中的相对表达量。结果表明，QDPR基因开放阅读框长度为417 bp，共编码139个氨基酸，编码蛋白为酸性不稳定蛋白，乌蒙凤鸡的QDPR基因序列与原鸡、日本鹌鹑和珠鸡的QDPR基因同源性较高；QDPR基因在不同性别乌蒙凤鸡不同生长阶段各组织中均有表达，其中，公鸡4、12和20周龄肺组织中QDPR基因表达量最高，极显著高于同一周龄其他组织（P<0.01），母鸡4周龄肝、脾、肺中QDPR基因相对表达量为所有时期最高，极显著高于12、20周龄（P<0.01），公鸡大多数组织中QDPR基因表达量随时间增长总体呈现出先升后降的趋势，母鸡为先降后升。结果显示，QDPR基因在乌蒙凤鸡内脏组织中的表达均高于肌肉组织，且公母鸡不同时期表达情况有所差异，试验结果可为进一步研究QDPR基因在鸡生长发育中的调控作用提供参考。     

miR-495-3p对山羊骨骼肌细胞增殖分化的影响

【目的】 研究miR-495-3p对山羊骨骼肌细胞增殖分化的影响及其在不同组织中的表达情况,为探究miRNA在肌肉发育中的调控机制提供理论基础。【方法】 利用生物信息学方法预测miR-495-3p的靶基因;实时荧光定量PCR检测miR-495-3p及其靶基因在心脏、肝脏、脾脏、肺脏、肾脏、背最长肌、腿肌等组织中的表达水平。构建miR-495-3p的过表达(miR-495-3p mimic、miR-495-3p mimic NC)及抑制物(miR-495-3p inhibitor、miR-495-3p inhibitor NC),在山羊骨骼肌细胞汇合度达60%~70%时进行转染,并用含2%马血清的培养基进行诱导分化,用CCK-8检测细胞增殖活力,实时荧光定量PCR检测过表达和干扰效率及增殖分化相关基因配对盒基因(paired-box 7,Pax7)、细胞蛋白周期E(Cyclin E)、肌细胞生成素(myogenin,myoG)、生肌因子5(recombinant myogenic factor 5,Myf5)的表达;构建miR-495-3p靶基因的野生型和突变型载体并转染至293T细胞,用双荧光素酶测定miR-495-3p与其靶基因的结合情况。【结果】 生物信息学预测结果表明,miR-495-3p的靶基因为心肌肌动蛋白α1(alpha cardiacactin 1,ACTC1),且miR-495-3p和ACTC1在1和10月龄山羊背最长肌中表达量均差异极显著(P<0.01)。miR-495-3p和ACTC1在心脏、肝脏、脾脏、肺脏、肾脏、背最长肌、腿肌6种组织中均有表达,在背最长肌中表达量较高。细胞转染结果表明,miR-495-3p mimic极显著促进miR-495-3p的表达(P<0.01),miR-495-3p inhibitor极显著抑制miR-495-3p的表达(P<0.01),说明miR-495-3p mimic和miR-495-3p inhibitor可用于后续试验。实时荧光定量PCR结果表明,过表达miR-495-3p促进骨骼肌细胞分化相关标志基因MyoG、Myf5的表达,抑制miR-495-3p则降低MyoG、Myf5基因的表达;过表达、干扰miR-495-3p对Pax7、Cyclin E基因的表达及细胞增殖活力均无显著影响(P>0.05)。双荧光素酶报告结果表明,miR-495-3p与ACTC1基因存在靶向关系。【结论】 miR-495-3p对山羊骨骼肌细胞的增殖无显著影响,但可通过靶向ACTC1基因促进山羊骨骼肌细胞分化。     

籽鹅和莱茵鹅MSTN和MyoG基因表达及其与屠宰性状的相关分析

本研究旨在探讨肌肉生成抑制素(myostatin,MSTN)和肌细胞生成素(myogenin,MyoG)基因对鹅骨骼肌生长发育的影响。试验以莱茵鹅和籽鹅为研究对象,采用实时荧光定量PCR方法测定MSTN、MyoG基因在籽鹅和莱茵鹅胸肌、腿肌中的差异表达情况,运用统计软件对基因表达情况与屠宰性状间的相关性进行分析。结果显示,莱茵鹅胸肌重和胸肌率极显著高于籽鹅(P<0.01),籽鹅胸肌MSTN、MyoG mRNA表达量极显著高于莱茵鹅(P<0.01);在腿肌中,籽鹅MSTN mRNA表达量极显著低于莱茵鹅(P<0.01),籽鹅与莱茵鹅MyoG mRNA表达量间无显著差异(P>0.05)。同一品种MSTN mRNA表达水平也存在差异,在籽鹅胸肌和腿肌中表达量存在极显著差异(P<0.01);在莱茵鹅胸肌和腿肌中表达量存在显著差异(P<0.05)。就MyoG mRNA而言,在籽鹅和莱茵鹅胸肌和腿肌中表达量存在极显著差异(P<0.01)。MSTN基因表达与屠宰性能相关性分析表明,胸肌中MSTN mRNA表达量与活重、胸肌重和胸肌率呈极显著负相关(P<0.01);胸肌和腿肌中MyoG mRNA表达量与活重呈极显著正相关(P<0.01)。说明MSTN和MyoG基因可能对肌肉生长分别起正负调控作用。     

Myf6基因在不同地方鸡种早期生长发育中的表达分析

为了探究Myf6基因在京海黄鸡和金茅黄鸡体内的差异表达规律，研究其在两个品种鸡中可能发挥的作用机制，试验采用实时荧光定量PCR技术对Myf6基因在不同品种鸡不同组织中的表达情况进行了检测。结果表明：京海黄鸡和金茅黄鸡公母鸡胸肌和腿肌中的Myf6基因相对表达量均比较高，且极显著高于脾脏(P<0.01);京海黄鸡心脏组织中的Myf6基因相对表达量显著高于脾脏(P<0.05),金茅黄鸡心脏组织中的相对表达量极显著高于脾脏(P<0.01);两个品种鸡的其余组织中Myf6基因相对表达量与脾脏相比差异不显著(P>0.05)。京海黄鸡公母鸡胸肌和腿肌组织中Myf6基因的表达规律基本一致，仅在4周龄胸肌组织中有稍微下降的趋势，但总体上呈现上升的趋势；腿肌组织中的相对表达量比较稳定，从0周龄到16周龄呈现上升趋势，16周龄时达到最高水平。金茅花鸡公母鸡胸肌和腿肌组织中Myf6基因的相对表达量从0周龄到4周龄呈上升趋势，8周龄有稍微下降的趋势，之后又呈现上升的趋势。在两个品种中，无论是公鸡还是母鸡，胸肌中Myf6基因的相对表达量始终显著高于腿肌(P<0.05)。说明Myf6基...     

鸡miR-15a对胸肌肌内脂肪沉积的调控作用

为研究微小核糖核酸miR-15a在固始鸡肌内脂肪沉积过程中的作用,利用茎-环荧光定量PCR(stem-loop quantitative RT-PCR)检测了miR-15a在固始鸡4个发育阶段脂肪代谢相关组织(腹脂、皮质、胸肌、肝脏)中的时序表达情况,测定了固始鸡胸肌肌内脂肪含量的动态变化,评价并验证了miR-15a与其预测靶基因在胸肌中的互作表达。结果表明:miR-15a的表达水平随胸肌发育逐渐升高,并与胸肌肌内脂肪含量呈明显的正相关;生物信息学分析显示,miR-15a预测靶基因在脂肪酸代谢、不饱和脂肪酸生物合成、脂肪酸降解等信号通路中显著富集。在胸肌中,miR-15a与其潜在靶基因酰基辅酶A氧化酶(ACOX1)、甾醇载体蛋白2(SCP2)、ATP结合盒子(ABCD3)的表达呈负相关,且在鸡成纤维细胞中过表达miR-15a可显著抑制这些靶基因的表达。上述结果说明,miR-15a可能通过其潜在靶基因ACOX1、SCP2、ABCD3参与鸡胸肌肌内脂肪沉积的调控。     

Expression of MSTN and MyoG Genes in Zi and Rhine Goose and their Correlation with Carcass Traits

The aim of this study was to investigate the effects of MSTN and MyoG genes on goose skeletal muscle growth.In this study,MSTN and MyoG genes expression were detected in breast and leg muscle of Zi and Rhine goose by Real-time PCR,and the correlations between genes expression levels and carcass traits were investigated.The results showed that the breast muscle weight and breast muscle rate of Rhine goose were extremely significant higher than Zi goose (P<0.01).MSTN and MyoG mRNA expression in breast muscle of Zi goose were significantly higher than that of Rhine goose,and the mRNA level of MSTN in leg muscle of Rhine was extremely significant higher than that of Zi goose (P<0.01),there was no significant difference of MyoG mRNA between Zi goose and Rhine goose (P>0.05).There was extremely significant difference between MSTN mRNA expression in breast muscle and leg muscle of Zi goose (P<0.01).MSTN mRNA expression in leg muscle was significantly higher than that of breast muscle of Rhine goose (P<0.05).There was extremely significant difference between MyoG mRNA expression in breast muscle and leg muscle of Zi goose and Rhine goose (P<0.01).There was a extremely significant negative correlation between MSTN mRNA expression in breast muscle with body weight,breast muscle weight and breast muscle percentage (P<0.01).There was a extremely significant positive correlation between MyoG mRNA expression in breast muscle and leg muscle with body weight (P<0.01).MSTN and MyoG gene might have positive and negative regulation effect on muscle growth.     

金茅黑鸡肉品质性状的关联转录组分析

为深入了解家禽肉品质性状形成的分子机制,本实验测定金茅黑鸡胸肌的肉品质,利用RNA-seq技术进行转录组测序,分析转录组基因表达量与肉品质性状间的关联性。结果显示:与金茅黑鸡14周龄肉色显著相关的基因5个,pH57个,剪切力33个,系水力17个;PLEKHS1、MUC等基因为影响肉色的候选基因,LHX9、SPDEF、GRXCR1等基因为系水力的候选基因,NPY、POMC、FGF20、MGAT4C和GUK1基因为影响肌肉pH的候选基因,EDAR基因和HS6ST3基因为肌肉剪切力的主要候选基因;功能分析发现这些基因在调节肌肉纤维、肌间脂肪沉积、脂类代谢、糖代谢过程和磷酸二酯酶水解等方面发挥重要作用,进而影响胸肌肌肉品质的形成。     

外源注射miR-499-5p对肉鸡腓肠肌转录组的影响

旨在探讨外源注射miR-499-5p对肉鸡腓肠肌肌纤维直径及基因表达的影响，完善miR-499-5p对肉鸡骨骼肌发育的分子调控机制。本研究选用14日龄健康大恒肉鸡公鸡16只，随机分为2组，每组8个重复，均正常饲喂日粮，分别在腓肠肌部位肌肉注射agomiR-NC（对照组，5 nmol）和agomiR-miR-499-5p（注射组，5 nmol），为保证试验效果，每3 d注射1次，共注射6次。注射完成后将16只肉鸡屠宰，采集腓肠肌组织进行肌纤维直径测定和转录组学（RNA-seq）分析。结果显示，注射miR-499-5p可导致肉鸡腓肠肌肌纤维直径显著减小（P<0.05）。转录组测序结果显示，注射组和对照组分别获得308 263 164和316 696 152条clean reads，与鸡参考基因组的比对率分别为90.29%和89.57%。与对照组相比，注射组共有差异表达基因18个，其中上调基因12个，下调基因6个（P<0.05）。其中差异表达最大的两个基因分别是原癌基因FOS（fos proto oncogene，FOS）和早期生长应答蛋白基因（EGR1）。GO富集分析发现，差异表达基因显著富集于跨膜受体蛋白丝氨酸/苏氨酸激酶信号通路、转录调控、脂质合成和骨骼肌再生等过程。KEGG富集分析结果显示，差异表达基因显著富集于丝裂原活化蛋白激酶（MAPK）信号通路和Wnt信号通路。利用生物信息学软件预测到miR-499-5p的靶基因共78个，GO富集分析发现，靶基因也显著富集于丝氨酸/苏氨酸激酶活性的负调控和转录因子活性等过程。本研究表明，注射miR-499-5p可改变肉鸡腓肠肌肌纤维直径从而影响肌纤维特性，差异表达基因中的FOS、EGR1、CYR61和WNT7A基因可能在miR-499-5p的作用下通过上述通路参与了脂肪沉积和肌肉发育的调控。     

猪miR-206前体的克隆、表达及生物信息学分析

研究骨骼肌中表达的miRNAs对深入探明骨骼肌生长发育的调控机制具有十分重要的意义。采用比较基因组学与生物信息学相结合的方法,从猪骨骼肌中克隆了猪miR-206的前体（Precursor）,对该miRNA前体的二级结构及与其它物种序列的同源性进行了分析。RT-PCR结果显示,该前体miRNA在猪大多数组织中均有表达,在骨骼肌中表达量最高。生物信息学分析表明miR-206前体序列在不同物种间具有较高的保守性。     

乌头驴与三粉驴骨骼肌中miRNAs的表达鉴定与分析

试验旨在分析microRNAs （miRNAs）在德州驴骨骼肌中的表达情况,探讨小RNA （small RNA）在驴骨骼肌发育中的调控机制。采集德州驴2个品系（乌头驴和三粉驴）的背最长肌组织,构建乌头驴背最长肌（WB）和三粉驴背最长肌（SB）的small RNA测序文库并进行测序,运用生物信息学技术对WB和SB之间差异表达的miRNAs进行分析,预测和解析可能影响驴产肉性状的miRNAs,从中挑选了6个表达量较高的miRNAs,利用实时荧光定量PCR技术来验证测序结果的准确性。结果表明,在WB和SB文库中分别鉴定出保守miRNAs有982和1 149个,新miRNAs有106和143个。在WB与SB比较组中基于条件|log₂（fold_change）|≥2、P≤0.05获得17个差异表达的miRNAs,其中表达量上调的有5个,下调的有12个。通过GO注释,发现差异表达的miRNAs显著富集于调控细胞迁移分化、细胞分裂和骨骼肌细胞收缩等生物学过程（P<0.05）,并且挖掘到与驴骨骼肌发育相关的bta-miR-378d_L+1R-1_1ss22CA、bta-miR-378d_R-2,其靶向NPY1R、GPR152、BEST1、KCNH8、SPAG9等多个基因。KEGG富集结果表明,miRNAs富集于Wnt、MAPK、泛素蛋白水解系统等与骨骼肌发育相关的信号通路中,其中在Wnt信号通路中共涉及到17个差异表达miRNAs的207个靶基因,尤其以PC-5p-84483_31、hsa-miR-186-3p_R-3和cfa-miR-329b_1ss19CT的靶基因最多。实时荧光定量PCR结果与miRNAs测序结果一致,并发现除eca-miR-181b外,其他5个miRNAs在SB中的表达水平均高于WB。本研究首次对驴骨骼肌miRNAs进行转录组学分析,揭示了乌头驴和三粉驴miRNAs的表达差异,研究结果为阐明德州驴骨骼肌生长发育的分子调控机制提供参考。