肉牛miR-133a通过抑制PAX7促进成肌细胞分化研究

于 堃，刘瑞莉，刘贤勋，柏学进，董雅娟，3

(1.青岛农业大学 动物科技学院，山东 青岛 266109；2.山东省黑牛繁育工程技术研究中心，山东 青岛 266109；3.山东布莱凯特黑牛科技股份有限公司，山东 淄博 256306)

肉牛是畜牧业中重要的经济动物，骨骼肌的生长发育不仅会影响肉牛的产肉量，还会对牛肉品质有一定的影响，进而改变食用口感。因此，调节肌肉发育过程、增加肌肉量是肉牛遗传育种工作的主要研究内容。当前，越来越多的研究集中在动物骨骼肌生长发育的机理上，在这一过程中，microRNAs(miRNAs)的研究为热点。

miRNAs是一类进化上保守的内源性非编码小RNA分子，长度约为18～22个碱基，广泛存在于动物、植物以及病毒等生物的基因组中，通过在转录后水平调控基因的表达，成为许多生物过程的重要调节因子。研究证明，家畜骨骼肌发育的几乎所有方面都涉及miRNAs的调控，包括细胞迁移、增殖、分化和凋亡等[1]。已有研究表明，许多miRNAs以组织特异性或阶段特异性的方式表达，而最具特征的肌肉特异性miRNAs是miR-1、miR-206和miR-133家族，它们在心肌和骨骼肌中特异性表达，在骨骼肌卫星细胞的增殖和分化中发挥关键作用[2]。Safa等[3]研究报道，miR-1可通过靶向PAX7和HDAC4的3′UTR区，抑制PAX7和HDAC4基因的表达，促进牛骨骼肌卫星细胞的肌原性分化。过表达miR-1/miR-206，可以使TNF-α表达下调，降低其对成肌细胞分化的抑制程度[4]。Mok等[5]研究表明，miR-133靶向调控Hedgehog通路中Gli3基因，抑制miR-133表达会破坏细胞增殖、细胞外基质沉积和上皮化，进而损害肌节的形成和生长，调控胚胎肌肉发生。miR-133通过抑制山羊骨骼肌卫星细胞HDAC4和SRF基因的翻译，参与调控山羊骨骼肌细胞增殖和分化过程[6]，通过调控TGFBR1的表达，参与鸭成肌细胞的增殖分化[7]。miR-133a通过靶向BIRC5基因，抑制骨骼肌细胞凋亡，影响肌肉质量[8]。此外，有的研究报道，miR-133a可通过降解生长激素基因的转录抑制子-胰岛素样生长因子 1(IGF1)，进而造成肌肉肥大等异常现象[9]。

目前，研究虽然证实了bta-miR-133a参与调控骨骼肌的生长发育，但并未见bta-miR-133a调控靶基因参与成肌细胞增殖与分化过程的直接相关报道。本实验室(动物胚胎工程中心)选用布莱凯特黑牛和鲁西黄牛背最长肌构建转录组高通量测序的sRNA数据库[10-11]，布莱凯特黑牛是运用体细胞克隆技术、超数排卵技术、胚胎分割及胚胎冷冻技术以及杂交育种与分子标记辅助育种的方法培育而成的优良肉牛新种质[12-13]；鲁西黄牛是我国优秀的地方良种，素有“五花三层肉”和“山东膘牛”的美称。本试验基于布莱凯特黑牛和鲁西黄牛的数据库，利用生物信息学方法预测bta-miR-133a的靶基因，分析其保守性，并对靶基因进行GO富集和KEGG通路富集分析，筛选并确定目标miRNA(bta-miR-133a)；以C2C12细胞作为细胞模型进行功能验证，转染bta-miR-133a 模拟剂、模拟剂阴性对照、bta-miR-133a抑制剂和抑制剂阴性对照，2%马血清诱导分化，检测靶基因和增殖分化标志性基因的表达量变化，探究bta-miR-133a通过调控靶基因参与肉牛骨骼肌发育过程的分子机制，旨在为揭示bta-miR-133a调节肉牛骨骼肌生长发育的机制提供科学依据，为肉牛肉质性状的分子选育提供参考。

1 材料和方法

1.1 主要材料与试剂

选取的12月龄布莱凯特黑牛和鲁西黄牛各3头，采集背最长肌组织置于液氮保存。小鼠C2C12细胞系和bta-miR-133a转染载体，购自普健生物(武汉)科技有限公司；miRNA反转录试剂盒，购自TaKaRa Bio USA；Lipofectamine 3000试剂，购自赛默飞世尔科技；CCK-8检测试剂盒(CK04)，购自东仁化学科技(上海)有限公司。

1.2 小RNA文库的构建测序与质量控制

用TRIzol试剂(Invitrogen)提取总RNA，使用RIN≥7.0且A260/280的比值在0.8～2.0的高质量RNA样品，构建序列文库。采用Illumina Genome Analyzer对文库进行测序。原始测序数据由Fastp软件处理，删除包含适配器，ploy-N和低质量读取数据，获得干净原始读取。所有下游分析均基于高质量的数据。通过比对分析软件Bowtie 1.1.2(设置允许一个错配)将sRNA测序的干净读数定位到参考基因组(Bos_taurus.NCBI.795.1)进行比对。

1.3 bta-miR-133a的生物信息学分析

为了探索 bta-miR-133a 的生物学功能，从miRBase(http：//www.mirbase.org/)检索牛(Bostaurus，Bta)、人(Homosapiens，Hsa)、黑猩猩(Pantroglodytes，Ptr)、大鼠(Rattusnorvegicus，Rno)、小鼠(Musmusculus，Mmu)、马(Equuscaballus，Eca)的miR-133a成熟序列，分析其保守性。使用 TargetScan 7.2(http：//www.targetscan.org/vert_72/)软件对bta-miR-133a进行靶基因预测，以物种“cow”为背景，选择目的miR-133a，Total context++score数值大小，代表miRNA与该靶基因结合的概率；根据Aggregate PCT，代表靶基因的保守定位可能性，获得预测靶基因。使用DAVID 6.8(https：//david.ncifcrf.gov/)和KOBAS 3.0(http：//kobas.cbi.pku.edu.cn/kobas3/)在线软件对预测所有bta-miR-133a靶基因进行GO分析和KEGG 通路富集分析。使用DAVID 6.8上传靶基因信息，选择背景基因物种“Bostaurus”，使用“Functional Annotation Tool”选择“KEGG_PATHWAY” 进行通路富集分析。使用KOBAS 3.0中的“Gene-list Enrichment”，上传基因信息，选择“GO”数据库进行功能富集分析。

1.4 双荧光素酶报告试验

使用胶回收纯化试剂盒(Shanghai Generay Biotech Co.，Ltd.)对PCR产物进行回收，将胶回收产物连接在pSI-Check2载体上，连接产物经感受态细胞DH5α转化后，37 ℃ 250 r/min摇菌14 h，用菌液进行PCR鉴定，将阳性克隆菌液进行测序。构建靶基因PAX7的双荧光素酶突变载体。HEK-293T细胞在96孔板中避光接种培养，同时将DMEM、0.8 mg/mL转染试剂(普健生物科技有限公司)、0.16 μgPAX7质粒和0.5 μg bta-miR-133a或NC mimic共转染到细胞中，转染48 h后收集细胞检测。使用双荧光素酶检测试剂盒(Promega Dual-Luciferase System)检测双荧光素酶活性。

1.5 成肌细胞转染及诱导分化

正常培养的C2C12细胞，吸弃原培养液，PBS洗后加入胰酶消化2 min后终止消化，剩余细胞悬液用培养基重悬，4 000个/孔(96孔板)进行铺板(200 μL)，3个复孔，5% CO2，37 ℃过夜培养。当细胞密度达到90%时，用Lipofectamine 3000试剂转染浓度为20 μmol/L的空白组(Blank)、模拟剂阴性对照组(NC mimic)、bta-miR-133a模拟剂(bta-miR-133a mimic)、抑制剂阴性对照组(NC inhibitor)和bta-miR-133a抑制剂组(bta-miR-133a inhibitor)，转染6 h后，换GM培养基，12 h后用 2%马血清诱导分化，分别再培养24(1 d)，48(2 d)，120(5 d)，144 h(6 d)观察细胞生长状态，48 h后提取细胞总RNA及总蛋白。

1.6 Cell Counting Kit-8(CCK-8)细胞增殖率检测

转染后24，48，72 h分别从培养箱中取出1个96孔板，向每孔中加入20 μL CCK8 溶液，培养4 h后，检测450 nm处的OD值。

1.7 实时荧光定量PCR(qRT-PCR)

随机选取4个miRNAs(bta-miR-133a、bta-miR-196b、bta-miR-200b、bta-miR-11977)进行qRT-PCR检测验证。检测成肌细胞增殖分化标志基因PCNA、CCND1、Myh1、Myod1和靶基因PAX7的表达量。使用TRIzol法提取空白组(Blank)、模拟剂阴性对照组(NC mimic)、bta-miR-133a模拟剂(bta-miR-133a mimic)、抑制剂阴性对照组(NC inhibitor)和bta-miR-133a抑制剂组(bta-miR-133a inhibitor)5组转染48 h后细胞总RNA。将提取的细胞总RNA逆转录为cDNA，以cDNA为模板进行荧光定量PCR。荧光定量PCR反应体系为20 μL：RNA-free Water 7.8 μL，上游引物0.6 μL(10 μmol/L)，下游引物0.6 μL(10 μmol/L)，SYBRGreen mix 10 μL，cDNA 1 μL(80 ng/μL)。PCR反应程序：94 ℃预变性10 min；94 ℃变性30 s，60 ℃退火30 s，72 ℃延伸40 s，40个循环[14-15]。以GAPDH作为内参基因，对PAX7、PCNA、CCND1、Myh1、Myod1基因进行标准化；以U6作为内参，对bta-miR-196b、bta-miR-200b、bta-miR-11977、Novel_445、Novel_130、Novel_334进行标准化，采用2-ΔΔCT值法计算。实时荧光定量PCR引物序列见表1。

1.8 数据分析

采用 SPSS 20. 0 统计分析软件中的One-way ANOVA进行方差分析和显著性检验。

2 结果与分析

2.1 背最长肌组织中 miRNA的筛选

使用Illumina基因组分析仪系统测序后，构建sRNA文库。获得原始数据，去除低质量reads(质量值低于20、碱基数超过1或含有N的reads)、不带有3′接头的reads、含有5′接头的reads、插入片段缺失的reads和长度小于18 nt的reads以及含有polyA的reads(一条reads中大于70%的碱基都是A)后得到高质量片段(clean reads)。通过比对分析软件Bowtie(设置允许一个错配)将sRNA测序的clean reads定位到参考基因组上，共获得26 061 624个clean reads。根据比对上参考序列的reads，与miRBase数据库的牛物种的参考序列进行比较，匹配到已知miRNA上的Total和Unique clean reads的总数量分别为29 849 254，37 946。根据差异miRNA筛选标准：|log2(Fold_Change)|≥1且 Q≤0.05，获得差异表达的miRNA 71个，其中上调39个，下调32个，其结果见图1-A。根据与肉牛骨骼肌发育相关的功能数据分析，筛选出表达量(Transcripts per million，TPM)前10位的差异miRNA分别为bta-miR-1、bta-miR-133a、bta-miR-206、bta-miR-26a、bta-miR-27b、bta-miR-143、bta-miR-378、bta-miR-10b、bta-let-7g、bta-miR-499，其结果见图1-B。

2.2 miRNA 数据验证

为了验证sRNA测序鉴定的miRNA表达谱，随机选取sRNA文库中已知表达量的6个miRNA(bta-miR-196b、bta-miR-200b、bta-miR-11977、Novel_445、Novel_130、Novel_334)，通过qRT-PCR验证测序结果的准确性，其结果见图2。通过比较qRT-PCR与高通量测序的结果发现(表2)，这6个miRNA在RNA-Seq和qRT-PCR分析中表现出一致的结果，表明RNA-Seq的测序结果是真实可靠的。

表2 miRNA测序表达量数据Tab.2 miRNA sequencing expression data

2.3 bta-miR-133a的保守性分析

使用miRBase数据库检索出牛在内的6个物种miR-133a的成熟序列，其结果见表3，miR-133a(UUUGGUCCCCUUCAACCAGCUG)在各个物种中高度保守。

表3 不同物种 miR-133a 成熟序列Tab.3 The mature sequence of miR-133a in different species

2.4 bta-miR-133a靶基因预测和富集分析

通过GO富集和KEGG通路富集分析，筛选出与肉牛骨骼肌发育相关的bta-miR-133a，使用

TargetScan 7.2软件进行靶基因预测，其结果见表4，预测到bta-miR-133a的靶基因594个，其中排名前10位的靶基因分别为：GNB4、GPM6A、TMEM167A、PLEKHA3、RAP2C、PAX7、ERI2、MAML1、SGPP1、AFTPH。使用DAVID 6.8(https：//david.ncifcrf.gov/)和KOBAS 3.0(http：//kobas.cbi.pku.edu.cn/kobas3/)在线软件对预测到的bta-miR-133a靶基因进行GO富集和KEGG通路富集分析。通过GO富集分析可知(图3)，bta-miR-133a靶基因显著富集于细胞增殖的负调控、胚胎骨骼系统发育、肌动蛋白细胞骨架组织的调节、肌动蛋白丝结合、细胞连接等与骨骼肌发育相关的GO条目中，其中靶基因PAX7富集于肌动蛋白丝结合条目中。通过KEGG通路富集分析可知(表5)，bta-miR-133a靶基因显著富集于黏附连接、肌动蛋白细胞骨架的调节、PI3K-Akt、FoxO、MAPK、Ras、TGF-β等信号通路与骨骼肌生长发育密切相关。使用TargetScan 7.2软件预测到bta-miR-133a与靶基因PAX7的靶标位区2个，2个靶标位区相差683个碱基，详细结合靶标位点见表6。

表4 bta-miR-133a的部分预测靶基因Tab.4 Some predicted target genes of bta-miR-133a

表5 bta-miR-133a的靶基因在KEGG通路中富集结果Tab.5 The enrichment results of bta-miR-133a target genes in the KEGG pathway

表5(续)

2.5 bta-miR-133a与靶基因PAX7的靶标关系验证

通过荧光素酶报告基因系统检测试剂盒检测荧光素酶活性水平，其结果见图4，与转染空白载体组相比，转染bta-miR-133a 组极显著降低荧光素酶活性(P<0.01)，转染bta-miR-133a-mu组的荧光素酶活性无显著变化(P>0.05)。表明bta-miR-133a在PAX73′-UTR存在结合位点，PAX7是bta-miR-133a的靶基因。

2.6 bta-miR-133a促进肌管形成

由图5可知，诱导分化1 d，各组之间无显著差异(P>0.05)。诱导分化2 d，与空白组相比，bta-miR-133a过表达组出现细胞融合，形成多核细胞；bta-miR-133a抑制组细胞融合数量无明显变化，大部分依然是单核细胞。诱导分化5 d，与空白组相比，bta-miR-133a过表达组有大量肌管形成，肌管直径明显增大；bta-miR-133a抑制组开始出现少量肌管，肌管直径较短。诱导分化6 d，与空白组相比，bta-miR-133a过表达组肌管直径明显增粗，肌管长度增大，数量增多；bta-miR-133a抑制组肌管的直径较短，肌管长度较短，数量无显著差异(P>0.05)。

红色柱状图代表生物过程；蓝色柱状图代表分子功能；绿色柱状图代表细胞组成；1.RNA聚合酶Ⅱ启动子转录的正调控；2.转录的正调控；3.细胞增殖的负调控；4.胚胎骨骼系统发育；5.DNA修复的正调控；6.肌动蛋白细胞骨架组织的调节；7.经典Wnt信号通路的负调控；8.生长；9.血管生成；10.转化生长因子β受体信号通路；11.核苷酸结合；12.受体信号蛋白丝氨酸/苏氨酸激酶活性；13.ATP结合14.mRNA结合；15.蛋白质异二聚化活性；16.泛素蛋白连接酶结合；17.蛋白酪氨酸磷酸酶活性；18.肌动蛋白丝结合；19.核；20.神经元投射；21.核质；22.细胞质；23.细胞-细胞粘附连接；24.细胞连接。

The red bar graph represents biological process；The blue bar graph represents molecular function；The green bar graph represents cellular component；1.Positive regulation of transcription from RNA polymerase Ⅱ promoter；2.Positive regulation of transcription；3.Negative regulation of cell proliferation；4.Embryonic skeletal system development；5.Positive regulation of DNA repair；6.Regulation of actin cytoskeleton organization；7.Negative regulation of classical Wnt signaling pathway；8.Growth；9.Angiogenesis；10.Transforming growth factor beta receptor signaling pathway；11.Nucleotide binding；12.Receptor signaling protein serine/threonine kinase activity；13.ATP binding；14.mRNA binding；15.Protein heterodimerization activity；16.Ubiquitin protein ligase binding；17.Protein tyrosine phosphatase activity；18.Actin filament binding；19.Nucleus；20.Neuron projection；21.Nucleoplasm；22.Cytoplasm；23.Cell-cell adherens junction；24.Cell junction.

2.7 成肌细胞增殖率检测

转染后，分别在24，48，72 h进行CCK-8细胞增殖率检测。由图6可知，24 h，各处理组之间无显著差异(P>0.05)。48 h，与空白组相比，bta-miR-133a过表达组显著降低成肌细胞的增殖率(P<0.05)；bta-miR-133a抑制组无显著差异(P>0.05)。72 h，与空白组相比，bta-miR-133a过表达组极显著降低成肌细胞的增殖率(P<0.01)；bta-miR-133a抑制组显著提高成肌细胞的增殖率(P<0.05)。

2.8 过表达、敲除bta-miR-133a对靶基因PAX7和增殖分化相关基因mRNA表达水平的影响

转染后，48 h提取细胞RNA，使用qRT-PCR检测PAX7、PCNA、CCND1、Myh1、Myod1的表达量。由图7可知，bta-miR-133a过表达组极显著抑制靶基因PAX7的表达(P<0.01)，bta-miR-133a抑制组显著促进PAX7的表达(P<0.05)。在细胞增殖方面(图7-B)，bta-miR-133a过表达组极显著抑制PCNA的表达(P<0.01)，显著抑制CCND1的表达(P<0.05)，bta-miR-133a抑制组显著促进PCNA和CCND1的表达(P<0.05)；在细胞分化方面(图7-C)，bta-miR-133a过表达组极显著促进Myod1和Myh1的表达(P<0.01)，bta-miR-133a抑制组显著抑制Myod1的表达(P<0.01)，对Myh1的表达没有显著影响(P>0.05)。

3 讨论与结论

骨骼肌的生长发育过程大致分为成肌细胞分化形成阶段、成肌细胞迁移增殖阶段和肌管形成阶段，涉及多基因表达和信号途径调控，过程极其复杂，miRNAs则是该过程的重要调控因子[16-18]。研究表明，miRNAs参与了几乎所有细胞过程的调控，关于miRNAs在骨骼肌发育中的潜在机制和作用的研究也越来越多，大多采用高通量测序建库结合生物信息学分析的方法[19-21]。本研究以动物胚胎工程中心实验室建立的布莱凯特黑牛和鲁西黄牛背最长肌测序小RNA文库为数据背景，采用GO富集和KEGG通路富集分析筛选出与肉牛骨骼肌发育密切相关的差异miRNAs(bta-miR-1、bta-miR-133a、bta-miR-206、bta-miR-26a、bta-miR-27b、bta-miR-143、bta-miR-378、bta-miR-10b、bta-let-7g、bta-miR-499)，选取TPM表达量较高的bta-miR-133a作为研究对象。据报道，从人和小鼠肌肉组织中克隆鉴定出肌肉特异性miR-206，可以促进骨骼肌细胞分化，抑制小鼠C2C12细胞增殖[22]。同时在Ling等[23]的研究中发现，miR-27b可通过调控PAX3基因促进山羊骨骼肌卫星细胞分化而抑制细胞增殖。miR-378a的水平升高会抑制成肌抑制因子，进而正反馈调控Myod促进成肌分化[24]，敲除miR-378a会影响骨骼肌中脂肪酸和葡萄糖氧化的能量代谢，进而影响肌纤维类型的发育[25]。本研究通过保守性分析，比对不同物种miR-133a成熟序列(UUUGGUCCCCUUCAACCAGCUG)的同源性，发现miR-133a在各个物种中高度保守，表明miR-133a在生物进化过程中可能具有重要功能。通过生物信息学在线网站分析，预测到bta-miR-133a靶基因594个，发现预测结果排名前10的靶基因(GNB4、GPM6A、TMEM167A、PLEKHA3、RAP2C、PAX7、ERI2、MAML1、SGPP1、AFTPH)中，多与癌症相关，例如：GPM6A在肺癌组织中低表达会影响患者预后，GPM6A可作为预测肺癌发生、转移及判断预后的有效分子标记物[26]；RAP2C与脑胶质瘤[27]、喉鳞状细胞癌[28]、乳腺癌[29]等癌症发展相关；AFTPH被认为是乳腺浸润癌(BRCA)、肺鳞癌(LUSC)和胰腺癌(PAAD)的潜在治疗靶标和预后生物标志物[30]。靶基因PAX7主要参与心肌与骨骼肌的生长发育，在成年期和肌肉再生过程中起主要作用，对卫星细胞的生物发生和存活至关重要[31-32]。同时也有研究表明，PAX7可通过与肌肉调节因子家族成员之间的交叉抑制，维持着卫星细胞的自我更新，从而决定了骨骼肌维持和修复所必需的增殖、分化和自我更新[33-34]。以上研究分析表明，PAX7与肌肉生长发育调节机制的功能相关性较强，因此，选取PAX7作为bta-miR-133a的靶基因进行下一步的研究。此外，bta-miR-133a的预测靶基因PLEKHA3和ERI2在目前的研究中报道较少，功能尚待进一步研究。

基于本试验是探究不同品种间肉牛骨骼肌生长发育的分子机制，对bta-miR-133a靶基因进行GO富集和KEGG通路富集分析，推测这些靶基因的潜在生物学功能。发现bta-miR-133a靶基因主要显著富集在细胞增殖的负调控、肌动蛋白细胞骨架组织的调节、胚胎骨骼系统开发、肌动蛋白丝结合等GO条目中和MAPK、PI3K-Akt、FoxO、Ras、TGF-beta、肌动蛋白细胞骨架的调节等与骨骼肌发育密切相关的信号通路中。MAPK信号通路由4个亚家族组成，其中细胞外信号调节激酶(ERK)主要调控细胞生长和分化，且Ras信号通路与MAPK信号通路中的信号蛋白存在联系，共同参与调控肌肉的发育过程；IGF1通过介导PI3K-Akt信号通路促进成肌细胞的增殖和骨骼肌生长[35-36]；FoxO 信号通路中的FOXO1、FOXO3和FOXO32广泛参与肌肉的生长发育过程[37]；TGF-β信号通路中的TGFβ2、TGFβ3可抑制卫星细胞的分化，进而影响骨骼肌的生长发育[38]；靶基因PAX7富集在肌动蛋白丝结合的GO条目和肌动蛋白细胞骨架的KEGG通路中，表明bta-miR-133a与骨骼肌的发育过程密切相关，且可能通过靶基因PAX7参与调控该过程。已有研究表明，PAX7是维持静止的卫星细胞和肌肉再生所必需的基因[39-40]，可作为关键上游调节因子参与肌肉的生长发育过程[41]，发挥着维持骨骼肌再生和修复的作用。进一步研究表明，PAX7与骨骼肌的生长发育调节机制的功能相关性最强，因此，选取PAX7作为bta-miR-133a的靶基因。

基于上述分析讨论，预测bta-miR-133a的种子序列(CCUGGUU)与靶基因PAX7的靶标结合位区有2个，都位于3′ UTR区，中间相差683个碱基。进一步通过双荧光素酶报告试验发现，bta-miR-133a与PAX7存在靶向关系，证明了PAX7是bta-miR-133a的靶基因。在试验开展前考虑到bta-miR-133a运用的局限性，进行了bta-miR-133a的PCR扩增试验，测序结果发现该miRNA在小鼠和牛组织上的片段扩增不受影响，且在前人的相关研究中，已有相类似的细胞模型验证试验，Kong等[42-43]研究发现，miR-17在小鼠C2C12细胞和牛的骨骼肌成肌细胞中，皆发挥着促进成肌细胞分化的相同作用；在miR-143、miR-24-3p等miRNA上也有相类似的研究[44-46]。因此，本试验以小鼠C2C12细胞作为细胞模型，进行细胞功能验证试验。转染Blank、NC mimic、bta-miR-133a mimic、NC inhibitor和bta-miR-133a inhibitor，诱导分化后发现，诱导分化5 d，与空白组相比，bta-miR-133a过表达组的细胞有大量肌管形成，肌管长度较长，而转染bta-miR-133a抑制组形成的肌管较少，大部分成肌细胞还在融合中，表明bta-miR-133a过表达组的成肌细胞向肌管分化的进程加快，而bta-miR-133a抑制组被抑制。在成肌细胞分化的过程中，观察到bta-miR-133a对肌管形成前期的作用效果并不明显，但在诱导分化5 d后，观察到各个处理组之间的细胞形态差异最明显，表明bta-miR-133a对成肌细胞分化的作用明显；在诱导分化6 d后，观察到bta-miR-133a抑制组的肌管直径显著增大、长度增长，表明转染bta-miR-133a inhibitor导致肌管形成的高峰期延后。虽然bta-miR-133a抑制组的成肌细胞分化过程较慢，但各个处理组的成肌细胞均进入肌管形成过程，最终形成肌管，表明bta-miR-133a inhibitor可以在一定程度上抑制成肌细胞分化，减缓其分化速度，但无法直接阻断成肌细胞的融合分化。综上表明，bta-miR-133a促进肌管分化，且主要在肌管形成后期发挥作用，这与邢义珅[47]、秦瑞峰等[48]、陈永乐等[49]研究报道的关于成肌细胞增殖分化的过程相一致。采用CCK-8检测细胞增殖率发现bta-miR-133a过表达组抑制细胞增殖的速率在24～48 h最大，bta-miR-133a抑制组促进细胞增殖的速率也在24～48 h最大，表明bta-miR-133a抑制成肌细胞增殖，且主要在细胞增殖的前期发挥主要作用。

研究表明，关于miRNAs对成肌细胞增殖分化方面的研究多采用增殖标志性基因PCNA、CCND1和分化标志性基因Myod1、Myh1的表达量来反映成肌细胞增殖分化的程度[50-53]。通过qRT-PCR检测转染后48 h成肌细胞中增殖标志性基因PCNA、CCND1和分化标志性基因Myod1、Myh1的表达量，发现与空白组相比，转染bta-miR-133a minic，极显著抑制PCNA，显著抑制CCND1的表达，极显著促进Myod1和Myh1的表达；转染bta-miR-133a inhibitor，显著促进PCNA和CCND1的表达，显著抑制Myod1的表达。表明bta-miR-133a参与调控增殖细胞中DNA的复制过程[54]，抑制细胞从G1期到S期的过程；促进肌球蛋白的合成，促进肌卫星细胞向成肌细胞转化，最终融合形成肌管[55]。基于以上试验结果分析表明，bta-miR-133a抑制成肌细胞增殖，促进成肌细胞分化。通过qRT-PCR检测PAX7的表达量，发现bta-miR-133a抑制PAX7表达，与Nguyen等[56]在尼罗罗非鱼骨骼肌发现的miR-133a与假定靶基因PAX7的负调控关系相一致。有的研究报道PAX7基因的表达量下降会促进成肌细胞分化[41]，进一步表明bta-miR-133a促进成肌细胞分化。

本研究选取bta-miR-133a的目标靶基因时，结合了生物学信息分析以及相关功能性研究，虽然只选取了PAX7，但是已有研究表明bta-miR-133a对成肌细胞的增殖分化过程的调控可能不仅仅通过靶基因PAX7，也可能通过调控其他调节因子来参与骨骼肌的发育过程。在该过程中，miR-133可通过靶向调控血清反应因子(SRF)[57]、NPTB[58]或者UCP2[59]来参与肌细胞增殖分化过程。因此，bta-miR-133a与PAX7以及其他靶基因的分子调控网络还需要进行深入的探索。

bta-miR-133a可能通过抑制靶基因PAX7的表达，抑制成肌细胞的增殖，促进成肌细胞的分化，进而参与调节肉牛骨骼肌的生长发育。