巴什拜羊不同发育阶段骨骼肌中miR-486靶基因的研究

2021-07-30张伟王世银高莉杨力伟邓双义刘晓娜石国庆甘尚权

中国农业科学 2021年14期

张伟，王世银，高莉，杨力伟，邓双义，刘晓娜，石国庆，甘尚权

张伟1,2，王世银1，高莉1，杨力伟1，邓双义1，刘晓娜1，石国庆2，甘尚权2

1新疆农业职业技术学院，新疆昌吉 831100；2新疆农垦科学院畜牧兽医研究所/省部共建绵羊遗传改良与健康养殖重点实验室，新疆石河子 832000

【】为了深入挖掘巴什拜羊不同发育阶段骨骼肌中miR-486的靶基因，为最终阐明巴什拜羊优良肉质性状形成的分子调控机制奠定基础，同时为巴什拜羊的持续选育提供理论依据。分别采集了巴什拜羊胎儿期40，50，60，80，100和120 d，以及出生当天，出生后30，60和90 d共计10个不同发育阶段的骨骼肌组织，分别利用胎儿期6个阶段的混合mRNA和出生后4个阶段的混合mRNA构建了胎儿期和出生后骨骼肌中miR-486调控靶基因的cDNA文库，并利用高通量测序技术对cDNA文库中的靶基因信息进行了深入挖掘。在对所得候选靶基因功能分析的基础上，选择10个候选靶基因，使用qRT-PCR检测其在巴什拜羊上述10个不同发育阶段骨骼肌中表达量的变化，结合课题组前期获得的miR-486在巴什拜羊不同发育阶段骨骼肌中的表达规律，对其靶向调控关系及生物学功能进行初步验证和分析。进而选择其中的4个候选靶基因使用荧光素酶报告载体试验和巴什拜羊骨骼肌卫星细胞中的靶向调控试验最终确认其靶向调控关系。通过对两个cDNA文库高通量测序数据的分析，胎儿期和出生后分别获得了123个和118个miR-486的靶基因，因其靶向调控关系未经实验验证，故暂称为候选靶基因。其中有96个为两个阶段共表达的候选靶基因，其余27个和22个分别为胎儿期和出生后骨骼肌中特异表达的候选靶基因。GO和KEGG分析结果表明所得候选靶基因显著富集于与肌肉分化和肌纤维发育相关的代谢和信号转导通路，诸如PI3k-Akt、MAPK、Wnt、Adherens junction和Regulation of actin cytoskeleton等。qRT-PCR检测结果表明所选10个候选靶基因在不同发育阶段巴什拜羊骨骼肌中均有表达，但其表达变化规律有所不同。、、、、、和在胎儿期巴什拜羊骨骼肌中呈高表达，出生后其表达量显著下调，而的表达量则呈相反的变化趋势，和在胎儿期和出生后巴什拜羊骨骼肌中的表达量未见显著变化。对其中4个候选靶基因的进一步试验验证，双荧光素酶报告载体试验结果表明，miR-486可以高效结合在其候选靶基因、、和的靶位点上，进而极显著抑制其后萤火虫荧光素酶的活性；巴什拜羊骨骼肌卫星细胞中的靶向调控试验结果亦表明miR-486可以显著下调上述4个候选靶基因的mRNA水平，抑制其生物学功能。所以可以确证在巴什拜羊骨骼肌中miR-486可以靶向调控、、和的表达。对巴什拜羊不同发育阶段骨骼肌中miR-486的靶基因进行了深入挖掘，全面解析了miR-486及其靶基因在巴什拜羊骨骼肌发育过程中参与的生物学过程和信号通路，所得候选靶基因数据可靠性高，可为阐明巴什拜羊优良肉质性状的分子调控机制奠定基础。

巴什拜羊；miR-486；靶基因；骨骼肌；cDNA文库

0 引言

【研究意义】巴什拜羊（Ovis aries）是我国一个优良的肉脂兼用型地方绵羊品种，以生产优质羔羊肉而著称，原产于新疆塔城地区裕民县，现主要分布在塔城地区的裕民县、额敏县和托里县等地，2014年存栏量约150万只[1]，该品种羊的养殖一直是当地牧民收入的主要来源。巴什拜羊素有贵族羊之称，其优良的肉质性状在国内外绵羊品种中亦不多见[2]。但是，近年来随着人们健康意识的提高以及肉食品消费观念的转变，低脂肪羊肉的需求量在逐渐增加。巴什拜羊由于胴体脂肪含量较高[3]，已无法满足市场的需求，亟需进行遗传改良[4]。本研究对巴什拜羊不同发育阶段骨骼肌中miR-486的靶基因进行深入研究，将为最终阐明巴什拜羊优良肉质性状形成的分子调控机制，进而指导其科学选育奠定基础。【前人研究进展】miRNA（microRNA）是一类小分子非编码RNA，约为22个核苷酸大小，在骨骼肌的发育过程中发挥着重要的调控作用[5-6]。miR-486首先是从人胚胎肝组织中鉴定出来，在心肌和骨骼肌中高表达，Alexander等发现miR-486作为PTEN/AKT 途径中重要的调节因子在抗肌萎缩蛋白缺陷型肌肉和杜氏肌肉营养不良症中（duchenne muscular dystrophy, DMD）发挥着重要的作用[7]。miR-486在敲除小鼠的骨骼肌中显著上调，体外过表达miR-486 可以诱导肌管肥大，抑制miR-486可以降低AKT 的活性，说明miR-486是骨骼肌中连接MSTN信号和IGF-1/Akt/ mTOR途径重要的中间调控因子[8]。目前对miR-486的研究主要集中在人类肌肉相关疾病方面，在绵羊等家畜肌肉发育中的相关研究鲜有报道。【本研究切入点】本课题组前期利用高通量测序技术对巴什拜羊不同发育阶段骨骼肌中miRNA的差异表达进行了研究，发现miR-486在出生后巴什拜羊骨骼肌中的表达量极显著上调，为胎儿期的100多倍[9]。进一步研究表明miR-486在巴什拜羊100 d胎儿骨骼肌和心肌中的表达量显著高于胃肠组织（＜0.05），而极显著高于肝、脾和肺等其他组织（＜0.01）[10]，提示miR-486在巴什拜羊骨骼肌发育中可能发挥着重要作用。miRNA通过调控其靶基因而发挥其生物学功能，但miR-486在巴什拜羊骨骼肌中的靶基因尚不清楚。【拟解决的关键问题】本研究通过构建miR-486靶基因cDNA文库结合高通量测序，对巴什拜羊胎儿期和出生后骨骼肌中miR-486的靶基因进行深入研究，并对部分靶基因进行试验验证，将为阐明miR-486在巴什拜羊骨骼肌发育过程中的调控机制奠定基础。

1 材料与方法

1.1 实验动物及骨骼肌样品采集

试验用巴什拜羊购自裕民县巴什拜种羊繁育中心。试验羊的饲养、前期处理及样品采集工作于2017年8月至2018年2月在新疆雪羚生物科技有限责任公司种羊场完成。

选择20只体况良好，胎次在2—3胎，来自不同家系的巴什拜羊（Ovis aries）母羊，统一进行饲养管理。15 d后进行同期发情和人工授精处理。从输精后的第15天开始每天分别在早上10点和下午6点观察母羊的发情状态，剔除返情的母羊，结合B超检查最终选择10只妊娠母羊进入后续实验。分别在妊娠期的第40，50，60，80，100和120天屠宰1头妊娠母羊，采集胎儿的四肢骨骼肌组织并投入液氮速冻，然后转入-80 ℃冰箱保存备用。其余4只母羊继续饲养至羔羊出生，并于出生当天，及出生后30，60和90 d屠宰羔羊，按照上述方法采集和保存骨骼肌组织。

1.2 RNA提取及基因组DNA消化

参照TRIzol（Invitrogen, USA）说明从上述骨骼肌样品中提取总RNA，使用Bioanalyzer 2100（Agilent, Germany）评估总RNA完整性。参照Oligotex®mRNA纯化试剂盒（Promega, USA）的说明，从总RNA中分离纯化mRNA，然后使用RNase-free DNase I（Ambion Inc., USA）消化所得mRNA，以去除微量基因组DNA污染，避免其对cDNA文库质量的影响。

1.3 cDNA文库的构建和高通量测序

将胚胎期6个阶段骨骼肌mRNA混合为一个样品，出生后4个阶段骨骼肌mRNA混合为另一个样品，参照甘尚权等[11]建立的方法分别构建胎儿期和出生后巴什拜羊骨骼肌miR-486的靶基因cDNA文库。现以胎儿期cDNA文库为例，简要说明建库的过程。将胎儿期骨骼肌mRNA混合样品分为两份，分别用反转录引物A（5¢-3¢: CTCGACTGAGTTGCCG，斜体加粗碱基为巢式引物的上游引物序列，下划线所示碱基为miR-486与其靶基因结合的种子序列）和B（5¢-3¢: AACGCG+oligodT20，下划线所示碱基为巢式引物的下游引物互补序列），参照PrimeScript™ Double Strand cDNA Synthesis Kit（TAKARA）试剂盒操作说明合成cDNA双链，分别称之为A库和B库。由于用于构建A库的反转录引物3¢末端的10个碱基为miR-486与其靶基因结合的种子序列，只有该序列与mRNA上的序列互补结合，才能完成反转录，所以A库中的cDNA应主要为miR-486的靶基因序列；用于构建B库的反转录引物3¢末端是20个oligodT，所以B库中的cDNA应是相应骨骼肌组织中所有表达功能基因序列的集合。使用酚/氯仿法抽提纯化A库和B库的cDNA双链后将二者等体积混合，95 ℃变性5 min，然后50 ℃保温5 min形成A库和B库的杂交cDNA文库。最后以杂交的cDNA为模板，使用巢式引物（上游引物5¢-3¢: TGAGTCGGCAACTCAGTCGAG和下游引物5¢-3¢: AGCGCAGCTCATCATCTGCAT）进行巢式PCR扩增。此时，只有A库中miR-486的靶基因链与其在B库中对应的互补链形成的杂交双链同时带有上下游巢式引物，所以在巢式PCR反应中可以被指数级扩增。而其他形式的杂交双链只带有上游或者下游巢式引物，在巢式PCR反应中只能被线性扩增。所以，经巢式PCR扩增后的产物，绝大部分为miR-486的靶基因序列，本研究中称之为miR-486靶基因cDNA文库。将构建的两个cDNA文库分别采用Illumina HiSeq 3000测序平台进行双向高通量测序。

1.4 高通量测序数据处理

测序原始数据使用SOAPnuke过滤，除去接头序列以及未知碱基超过10%和低质量碱基（Q≤10）超过50%的低质量reads。获得的clean reads使用SOAPaligner和SOAP2以牛（v3.1）和绵羊（v3.1）基因组序列为参考进行拼接组装和注释。

1.5 候选靶基因序列筛选及分析

按照靶基因cDNA文库建库原理，所得miR-486靶基因序列应该含有反转录引物A的末端10碱基序列（GTCATGTCCT）。因此以该10碱基序列为种子序列对上述分析获得的基因序列进行比对筛选，含有该序列的基因序列即为miR-486靶基因序列，由于所得靶基因尚未经过实验验证，存在一定数量的假阳性结果，所以暂时称为候选靶基因。为了深入了解所得候选靶基因的生物学功能，利用David在线软件进行基因本体论[12]（gene ontology, GO）（http://geneontology. org）分析，从细胞组分、生物学过程和分子功能3方面对其进行功能注释，然后利用京都基因和基因组百科全书[13]（Kyoto encyclopedia of genes and genome, KEGG）数据库（http://www.kegg.jp/kegg ）分析所得miR-486候选靶基因参与调控的信号通路，设置≤0.05为显著性阈值。

1.6 部分候选靶基因在巴什拜羊不同发育阶段骨骼肌中表达量检测

在以上靶基因cDNA文库测序结果分析的基础上，查阅近年来骨骼肌发育相关信号通路的文献报道，从靶基因cDNA文库中选择10个与骨骼肌生长发育调控密切相关的候选靶基因，以-actin为内参基因，以出生后第90天的表达量为对照，使用qRT-PCR对其在巴什拜羊不同发育阶段骨骼肌中的表达量进行相对定量，并参照课题组前期检测的miR-486的表达规律[10]，对10个候选靶基因与巴什拜羊骨骼肌发育的关联性，及其与miR-486的靶向调控关系进行初步分析。10个候选靶基因及内参基因在Ensemble数据库（http://asia.ensembl.org/index.html）中的登录号及qRT-PCR引物信息见表1。引物使用oligo6.0设计，由生工生物工程（上海）有限公司合成。mRNA的反转录参照PrimeScript Ⅳ 1st strand cDNA Synthesis Mix（TAKARA）完成，qRT-PCR扩增采用20 μL体系：TB Green Premix ExⅡTli RNaseH Plus）（2×）（TAKARA, 大连）10 μL，上、下游引物（10 μmol·L-1）各 0.8 μL，cDNA模板 2 μL，ddH2O 6.4 μL。扩增条件：95 ℃预变性30 s；95 ℃ 10 s，58 ℃ 20 s，72 ℃ 20 s，40个循环。扩增反应结束后以 0.1 ℃/s 的速度进行熔解曲线分析。qRT-PCR检测由LightCycler 96（Roche, Germany）完成。

表1 10个候选靶基因qRT-PCR引物信息表

*登录号指基因在Ensemble数据库（http://asia.ensembl.org/index.html）中的Transcript ID

*Accession number means the transcript ID of gene in Ensemble database (http://asia.ensembl.org/index.html)

1.7 miR-486与部分候选靶基因靶向调控关系的验证

由于通过构建靶基因文库获取的靶基因存在一定比例的假阳性结果，为了验证上述靶基因文库的可靠性，结合上述10个候选靶基因在巴什拜羊不同发育阶段表达量的变化规律，以及课题组前期检测到的miR-486的表达规律[10]，从上述10个候选靶基因中选择、、和4个表达量与miR-486的表达量呈明显负相关的候选靶基因作进一步验证。

1.7.1 双荧光素酶报告载体试验在候选靶基因序列上miR-486结合靶位点上下游约100 bp片段范围设计PCR引物，以巴什拜羊基因组为模板PCR扩增获得靶位点区域约200 bp大小的野生型片段，同时根据种子序列设计突变引物，与上述PCR引物进行重叠PCR对靶位点进行点突变，获得靶位点处4个碱基突变的突变型片段。将野生型片段和突变型片段经I和dIII双酶切后连入PGL4双荧光素酶报告载体（Promega, USA），并使用双酶切和测序确认载体构建的正确性。上述载体构建引物信息见表2。

表2 荧光素酶报告载体构建引物

用于载体构建引物序列中小写字母为限制性内切酶I和dIII的切割序列及保护碱基；用于靶位点突变引物中小写字母为致突变碱基

The bases in lowercase of vector construction primer were restriction enzyme cutting sites (italics and bold) and protective bases ofI anddIII restriction enzyme cutting sites and protective bases, the ones of mutating bind site primer were mutated bases of the miR-27b target site in MSTN 3-UTR

参照Lipofectamine 2000（Invitrogen, USA）的操作说明，将构建好的报告载体及miR-486 mimics（表3）分别转入Hela细胞，培养48 h后裂解细胞，按照Dual-Luriferase Reporter Assay System（Promega, USA）的操作说明检测萤火虫荧光素酶和海肾荧光素酶的活性。

1.7.2 骨骼肌卫星细胞中miR-486干扰表达试验检索miRBase数据库，根据其收录的bta-miR-486序列（MIMAT0009329）设计其模拟物miR-486 mimics，抑制物miR-486 inhibitor及其阴性对照negative control（表3）并由上海吉玛基因合成。扩大培养巴什拜羊骨骼肌卫星细胞（通过原代培养获得），参照Lipofectamine 2000（Invitrogen, USA）的操作说明将miR-486 mimics，miR-486 mimics negative control，miR-486 inhibitor和miR-486 inhibitor negative control分别转入骨骼肌卫星细胞，以不转染任何试剂的细胞组作为空白对照。转染后的细胞继续培养至24、48和72 h时分别消化收集细胞，按照1.2所述方法提取总RNA，并按照1.6所述方法及引物检测不同处理组骨骼肌卫星细胞中4个候选靶基因mRNA的水平。

表3 miR-486 mimics、miR-486 inhibitor及其阴性对照序列

1.8 数据分析

qRT-PCR结果采用2-△△Ct法进行计算，双荧光素酶检测结果采用萤火虫荧光素酶活性/海肾荧光素酶活性计算。试验均重复3次，结果用“平均值±标准差”表示，使用SPSS13.0进行单因素方差分析，<0.05表示差异显著，<0.01表示差异极显著。

2 结果

2.1 靶基因cDNA文库高通量测序及分析

胎儿期和出生后miR-486靶基因cDNA文库高通量测序分别得到6.21M和5.94M原始数据，剔除接头序列以及低质量reads后分别获得5.73M和5.42M的clean reads。进一步拼接组装、剔除冗余序列后分别得到283条和276条Unigenes，平均序列长度分别为673 bp和681 bp。将胎儿期和出生后Unigenes序列进行blastn比对，发现分别有68和61个为胎儿期和出生后特异表达Unigenes，有215个为两个时期的共表达Unigenes。

2.2 候选靶基因注释及功能分析

按照靶基因cDNA文库建库原理，所得miR-486候选靶基因序列应该含有反转录引物A的末端10碱基序列（GTCATGTCCT）。因此以该10碱基序列为种子序列对上述分析获得的Unigenes序列进行比对筛选，胎儿期和出生后分别获得123个和118个符合条件的Unigenes，其中96个为两个阶段共表达Unigenes，其余27个和22个分别为胎儿期和出生后骨骼肌中特异表达的Unigenes。由于miR-486与所得Unigenes的靶向调控关系尚未经过试验验证，所以暂时称之为miR-486的候选靶基因。

GO分析结果表明，胎儿期和出生后分别有96个和87个候选靶基因得到注释信息，胎儿期的候选靶基因显著富集于26个细胞组分条目，32个生物学过程条目和18个分子功能条目。出生后的候选靶基因显著富集于20个细胞组分条目，36个生物学过程条目和23个分子功能条目，富集基因数前10的GO条目及富集于相应条目的候选靶基因占该阶段检测到候选靶基因的比例见图1。

KEGG通路富集分析结果表明，胎儿期和出生后候选靶基因显著富集于多个与肌肉发育相关的信号通路，如PI3k-Akt、MAPK、Wnt和Adherens junction等，提示miR-486可能通过调控一些与肌肉发育相关的基因及信号通路，在巴什拜羊骨骼肌发育调控中发挥着重要的作用。胎儿期和出生后候选靶基因显著富集的前20个信号通路见图2。

2.3 miR-486及其部分候选靶基因在巴什拜羊不同发育阶段骨骼肌中的表达规律

qRT-PCR检测结果表明所选10个候选靶基因在巴什拜羊骨骼肌中均有表达（图3），但在不同发育阶段的表达量存在差异。表达量在胎儿期呈持续下降趋势，在胎儿期第40和50 天表达量最高，极显著高于其他时期（＜0.01），胎儿期第60天开始逐步下降，至出生后降至最低且保持相对稳定（图3-A）。表达量从胎儿期至出生后呈持续下降趋势，在胎儿期第40、50和60天表达量最高且差异不显著（＞0.05），但极显著高于其他时期（＜0.01）；胎儿期第80、100和120天时表达量下调，但仍极显著高于出生后各阶段（＜0.01）；出生当天至出生后第30天表达量进一步下降，这2个时期差异不显著（＞0.05），但极显著高于出生后第60和90天的表达量（＜0.01），至出生后第90天降至最低（图3-B）。和在胎儿期前期高表达，胎儿期第40、50、60和80天均极显著高于其他各阶段（＜0.01），但至胎儿期第100天迅速下调，至出生后第90天一直保持较低的表达水平（图3-C和D）。表达量在胎儿期第40、50和60天表达量最高且差异不显著（＞0.05），至胎儿期第80和100天时表达量极显著下调（＜0.01），至胎儿期第120天时又出现极显著上调（＜0.01），但至出生后即出现极显著下调（＜0.01），直至出生后第90天降至最低水平（图3-G）。在胎儿期第40和50天时呈高表达，至胎儿期第60和100天时出现2次极显著下调（＜0.01），而胎儿期第80和120天时出现2次极显著上调（＜0.01），至出生后一直呈现持续下调趋势（图3-H）。基因在胎儿期第40、50、60和80天表达量最高且差异不显著（＞0.05），胎儿期第100 d时表达量极显著下调（＜0.01），胎儿期第120 d时又出现极显著上调（＜0.01），至出生后再次出现极显著下调（＜0.01），并在出生当天及出生后30、60和90 d保持稳定表达（图3-I）。上述7个基因在巴什拜羊不同发育阶段骨骼肌中的表达变化存在一定差异，但是总体而言呈胎儿期高表达，出生后低表达的趋势。已有的研究结果表明miRNA对其靶基因的调控一般为负调控，而课题组前期研究发现miR-486在巴什拜羊胎儿期总体呈低表达，出生后表达量逐渐上调[10]。由此可以推断，在上述7个候选靶基因与miR-486巴什拜羊骨骼肌中存在靶向调控关系的可能性很大。

A和B分别为胎儿期和出生后候选靶基因GO分析结果，Y轴为GO条目名称，X轴为相应条目富集的候选靶基因占注释于该类型候选靶基因的比例

A和B分别为胎儿期和出生后候选靶基因参与通路分析结果，Y轴为pathway名称，X轴为富集因子，气泡大小反映了富集到某个信号通路中的候选靶基因的数量占该通路已发现参与基因的比例，颜色变化反映了某个通路富集的显著性水平

A-J分别为10个候选靶基因在巴什拜羊不同发育阶段骨骼肌中相对表达量的qRT-PCR检测结果，横轴为巴什拜羊骨骼肌不同发育阶段，1-10分别指胎儿期第40、50、60、80、100和120天，以及出生当天，出生后30、60和90 d共10个不同的发育阶段，纵轴为相应基因的相对表达量，图中相同字母或无字母表示差异不显著（P＞0.05），不同小写字母表示差异显著（P＜0.05），不同大写字母表示差异极显著（P＜0.01）。下同

和在胎儿期和出生后巴什拜羊骨骼肌中的表达量未见显著变化（图3-E和3-J），推测这2个基因的表达可能是巴什拜羊骨骼肌生长发育所必需，故其表达量保持相对稳定；而胎儿期呈低表达，但表达量逐步上调，至出生后第90 天达到较高的表达水平（图3-F）。结合miRNA与其靶基因靶向调控的规律，推测这3个候选靶基因与miR-486存在靶向调控关系的可能性不大。

2.4 荧光素酶报告载体验证

为了进一步验证miR-486与候选靶基因的靶向调控关系，从上述10个候选靶基因中选择、、和4个候选靶基因进行双荧光素酶报告载体实验验证。将4个候选靶基因野生型（WT）和突变型（Mut）靶位点序列插入到双荧光素酶报告载体海肾荧光素酶和萤火虫荧光素酶之间（图4-A），荧光素酶活性检测结果表明miR-486与插入野生型靶位点序列报告载体共转染细胞中萤火虫荧光素酶的活性均极显著低于对照组（＜0.01），而miR-486与插入突变型靶位点序列报告载体共转染细胞中萤火虫荧光素酶的活性与对照组无显著差异（＞0.05，图4-B），说明miR-486可与候选靶基因靶位点序列高效结合，从而抑制其后萤火虫荧光素酶的表达，由此也初步证明miR-486与、、和4个基因存在靶向调控关系。

2.5 骨骼肌卫星细胞中靶向调控关系验证

巴什拜羊骨骼肌卫星细胞中的验证结果表明，对于和，当转入miR-486 mimics造成细胞中miR-486的表达水平呈上调状态时，继续培养至24、48和72 h的细胞中均检测到和的mRNA水平随着培养时间的延长呈显著下调趋势（＜0.05），空白对照组、miR-486 inhibitor组和两个阴性对照组细胞中和的mRNA水平未见显著变化（＞0.05）（图5-A和C）；的mRNA水平在培养至24 h细胞中未见显著变化（＞0.05），48 和72 h的mRNA水平变化与和相同（图5-D）；mRNA水平在miR-486 inhibitor组24、48和72 h的细胞中均呈显著上调趋势（＜0.05），在miR-486 mimics组48和72 h的细胞中均呈显著下调趋势（＜0.05），其余处理组之间无显著差异（＞0.05）（图5-B）。

A为候选靶基因片段插入双荧光素酶报告载体位置及4个候选靶基因野生型（WT）和突变型（Mut）靶位点序列，序列中带下划线的碱基为通过重叠PCR获得的突变碱基；B为插入野生型和突变型靶位点片段荧光素酶报告载体的荧光素酶相对活性检测结果，NC指miR-486的阴性对照（Negative control）

综合荧光素酶报告载体实验和骨骼肌卫星细胞中靶向调控关系验证实验结果，可以确认在巴什拜羊骨骼肌中miR-486可以靶向调控、、和4个基因。

A、B、C和D分别为候选靶基因PTEN、Foxo1、IGF1R和 PIK3R1在转染细胞中24、48和72 h的相对表达量

3 讨论

动物的骨骼肌是由单核的成肌细胞相互融合而成的多核肌纤维构成的[14]。成肌细胞来源于中胚层，迁移至躯干和四肢等骨骼肌发育的位置后，成肌细胞先增殖达到一定的数量，然后相互融合为肌管，并最终发育成肌纤维[15-16]，这一过程一直持续至动物出生前才结束。在不同的动物中，肌管的形成可以明显的区分为两至三个时期。就大部分动物而言，第一波和第二波肌管分别来自于胚胎和胎儿成肌细胞，并分别发育成初级和次级肌纤维[15]。但在一些大型动物中，例如牛、绵羊、猪和人，在妊娠后期还会有第三波肌管形成[17-20]，这可能与这些动物具有发达的骨骼肌有一定的关系。利用光学显微镜定量研究绵羊骨骼肌的超微结构，发现胎儿期的40—80 d和100—120 d为骨骼肌生长发育的关键时期[17]。李雪娇等对中国美利奴羊骨骼肌发育进行了系统研究，发现胎儿期第75天前后肌纤维的形成基本完成，至第105天时肌纤维的数量基本恒定，此后骨骼肌的发育主要表现为肌纤维直径和长度的增加[21-22]。石田培等利用全转录组测序对中国美利奴羊胎儿期骨骼肌中circRNA的表达情况进行研究，发现胎儿期第85、105和135天这3个阶段显著富集了大量与肌肉发育相关的基因和通路[23]。提示这些阶段骨骼肌中的基因表达模式发生了变化，以支持骨骼肌细胞不同生长发育状态的转换和维持。鉴于此，本研究进一步细化采样阶段，对巴什拜羊胎儿期第40、50、60、80、100和120天，以及出生当天，出生后第30，60和90天骨骼肌中miR-486的靶基因进行深入挖掘，发现有96个候选靶基因在胎儿期和出生后均有表达，有27个和22个候选靶基因分别在胎儿期和出生后特异表达。对其中10个候选靶基因mRNA水平的定量研究结果表明，即使是共表达的候选靶基因，在巴什拜羊不同发育阶段骨骼肌中的表达量亦有很大的变化，、和等7个候选靶基因在胎儿期巴什拜羊骨骼肌中呈高表达，出生后其表达量显著下调，而在胎儿期低表达，出生后表达量出现显著上调，仅有和两个候选靶基因的表达量未见显著变化。结合前人研究结果，推测这些候选靶基因可能在维持不同发育阶段骨骼肌细胞增殖和分化过程中发挥着重要的作用。

miR-486由（）的第40个内含子编码[24]，第40—42外显子编码的一个小分子Ank1蛋白仅在肌肉中表达，可将肌小节连接在肌质网上[25]。到目前为止未发现miR-486尚有其他家族成员，且仅在哺乳动物中表达，其序列高度保守，在鸟类、鱼类、两栖类和更低等的后生动物中均未检测到miR-486的表达[24]。miR-486在巴什拜羊的骨骼肌和心肌中高表达[10]，且其表达量在出生后巴什拜羊骨骼肌中极显著上调[9]，推测miR-486在巴什拜羊骨骼肌发育过程中可能发挥着重要的作用，本研究中发现的miR-486部分候选靶基因已证明与骨骼肌发育调控和机能维持密切相关，初步支持这一推测的可靠性。（phosphatase and tensin homolog）是一个重要的肿瘤抑制基因，可以通过抑制PI3K的活性以及降低Akt的磷酸化水平影响下游mTOR信号通路，最终影响细胞的生长和存活[26]。在成年人的骨骼肌中，骨骼肌卫星细胞处于静止、激活和分化的动态平衡状态，但是当被敲除，大量的骨骼肌卫星细胞将不经过增殖即进入激活、提前成熟分化状态，进而造成骨骼肌中卫星细胞的消耗而最终影响骨骼肌损伤的修复和再生[27]。可以促进小鼠肌细胞融合，同时参与骨骼肌纤维分化的调控，在富含快肌纤维的骨骼肌组织中的表达量显著高于富含慢肌纤维的骨骼肌组织[28-29]。在小鼠肌细胞中，miR-486抑制其靶基因和的表达进而激活PI3K/Akt信号通路参与肌肉的生长的调控和内稳态的维持[24]。在杜氏肌应用不良症（duchenne muscular dystrophy, DMD）患者的肌肉组织中miR-486的表达量显著下调，同时在DMD动物模型骨骼肌中过表达miR-486，可降低DOCK3和PTEN的水平，进而提高AKT的磷酸化水平并使肌营养不良症状得到明显缓解[30]。慢性肾病（chronic kidney disease, CKD）患者由于E3连接酶的激活，泛素蛋白酶体系统被活化，进而引起肌肉中蛋白质的分解，但是若上调肌肉中miR-486的水平可显著抑制E3连接酶的表达，使患者的肌肉量明显增多，进一步研究表明这一过程是miR-486通过抑制的转录，同时增加Foxo1磷酸化水平，降低PTEN磷酸化水平，最终使Akt的磷酸化水平显著升高进而激活相关信号通路而实现[31-32]。本研究已证实，在巴什拜羊骨骼肌中miR-486可以靶向调控、、和4个基因的表达，但其在巴什拜羊骨骼肌发育过程中的生物学功能尚有待后续试验阐明。

KEGG分析结果表明，miR-486通过调控其候选靶基因的表达同时参与多个信号通路的调控。例如同时参与mTOR signaling pathway、PI3K-Akt signaling pathway、FoxO signaling pathway、MAPK signaling pathway、Ras signaling pathway 和Adherens junction等信号通路的调控，同时参与Phosphatidylinositol signaling system、PI3K-Akt signaling pathway、FoxO signaling pathway、mTOR signaling pathway和Insulin signaling pathway等信号通路的调控，很多信号通路已证明与骨骼肌生长发育的调控密切相关[33-37]，在小鼠中的研究结果表明，其骨骼肌中miR-486的表达受的调控，通过miR-486转录水平的调控可将myostatin信号通路和IGF-1/Akt/ mTOR信号通路有效整合[8]。在人类骨骼肌中，miR-486可通过调控PI3K-Akt、ubiquitin proteasome、FoxO和mTOR信号通路，使骨骼肌组织很好地适应有氧运动训练[38-39]，提示miR-486可能在沟通和协调巴什拜羊骨骼肌内相关信号通路中发挥着重要作用。

在细胞质中miRNA首先与Argonaute蛋白组装成为miRNA诱导的沉默复合体（miRISC），然后miRNA通过其种子序列与靶基因mRNA互补结合引导miRISC结合到其靶基因mRNA的靶位点上，进而调控其靶基因的翻译[40-41]，但在动物中miRNA是通过抑制其靶基因mRNA的翻译还是诱导其降解实现其对靶基因的调控尚存在争议，有研究结果表明在动物细胞中miRNA主要是抑制其靶基因mRNA的翻译，并不会诱导mRNA的降解[42]，但也有相关研究检测到了miRNA诱导的mRNA降解[43]。或许Djuranovic等的研究可以更好地解释这一问题，他们利用果蝇S2细胞建立了一个可控的表达系统，并以此研究了miRNA对其靶基因调控的动力学特征，发现miRNA不论跟天然的或者人工合成的3′ UTR区靶位点结合，首先是抑制其翻译，进而引起mRNA的脱腺苷化和降解，所以miRNA可能是先抑制其靶基因mRNA的翻译，然后再诱导其降解来实现其对靶基因的调控，若在抑制阶段检测细胞中mRNA的水平，则会得出miRNA不会诱导mRNA降解的结论[44- 45]。本研究中上调巴什拜羊骨骼肌卫星细胞中miR-486的水平，至24 h其靶基因和的mRNA水平未见显著降低（＞0.05），但是至48和72 h时和的mRNA水平显著低于对照组（＜0.05），与此不同的是和的mRNA水平至24 h时即可检测到显著降低（＜0.05）。同为miR-486的靶基因，其调控方式却存在不同，其机理尚有待进一步研究。通过靶位点分析发现，miR-486在的mRNA上有3个靶位点，其中一个位于CDS区，在的mRNA上有2个靶位点，而在和的mRNA上仅检测到1个靶位点，这种调控方式的差异是否与靶位点的位置和数量有关尚有待后续实验验证。另外，当下调巴什拜羊骨骼肌卫星细胞中miR-486时，仅检测到的mRNA水平较对照组显著升高，其他3个靶基因mRNA水平未见显著变化（＞0.05），其机理亦有待进一步阐明。

4 结论

本研究采用构建cDNA文库结合高通量测序的方法首次对巴什拜羊不同发育阶段骨骼肌中miR-486的靶基因进行了深入研究，在胎儿期和出生后骨骼肌中分别检测到123个和118个miR-486的候选靶基因，功能分析结果表明候选靶基因显著富集于PI3k-Akt、Wnt和MAPK等已证明与肌肉发育调控密切相关的信号通路。其中、、、、、和7个候选靶基因在胎儿期巴什拜羊骨骼肌中呈高表达，出生后其表达量显著下调，而候选靶基因呈相反的表达趋势，提示这些基因在巴什拜羊胎儿期和出生后骨骼肌发育状态转换中可能发挥着重要作用。选择、、和4个候选靶基因进行实验验证，确认miR-486与所选4个候选靶基因存在靶向调控关系，说明本研究获得的miR-486的靶基因数据可靠性高。

[1] 努尔拉提汗·木哈买提别克, 赛力克·木勒达汗, 努尔古丽·热孜别克, 决肯·阿尼瓦什. 新疆巴什拜羊青色和白色毛绒品质比较分析. 新疆畜牧业, 2015, 11:25-27.

MUHAMAITIBIEKE N, MULEDAHAN S, REZIBIEKE N, ANIWASHI J. Analysis of white and cyan wool quality of Bashbay sheep. Xinjiang Animal Husbandry, 2015, 11:25-27. (in Chinese)

[2] 决肯·阿尼瓦什, 克木尼斯汗, 海拉提, 季文凭, 杜曼.巴什拜羊瘦肉型新品系的培育及其生产性能分析. 新疆农业科学, 2010, 47(2): 406-409.

ANIWASHI J, KEMUNISIHAN, HAILATI, JI W P, DU M. Analysis of breeding and production performance of new strain lean meat type Bashbay sheep. Xinjiang Agricultural Sciences, 2010, 47(2): 406-409. (in Chinese)

[3] 韩业东, 决肯·阿尼瓦什, 帕娜尔·依都拉, 李齐发, 谢庄. 巴什拜羊群体遗传多样性与遗传分化的研究. 南京农业大学学报, 2011, 34(1): 123-127.

HAN Y D, ANIWASHI J, YIDULA P, LI Q F, XIE Z. Study on genetic diversity and genetic differentiation among Bashbay sheep populations. Journal of Nanjing Agricultural University, 2011, 34(1): 123-127. (in Chinese)

[4] 刘真, 王慧华, 刘瑞凿, 吴明明, 张淑珍, 张莉, 赵福平, 杜立新, 魏彩虹．不同尾型绵羊全基因组选择信号检测．畜牧兽医学报，2015, 46(10): 1721-1732.

LIU Z, WANG H H, LIU R Z, WU M M, ZHANG S Z, ZHANG L, ZHAO F P, DU L X, WEI C H. Genome-wide detection of selection signatures of distinct tail types in sheep populations. Chinese Journal of Animal and Veterinary Science, 2015, 46(10): 1721-1732. (in Chinese)

[5] MIRETTI S, MARTIGNANI E, ACCORNERO P, BARATTA M. Functional effect of miR-27b on myostatin expression: a relationship in piedmontese cattle with double-muscled phenotype. BMC Genomics, 2013, 14(1):1-8.

[6] CLOP A, MARCQ F, TAKEDA H, PIROTTIN D, TORDOIR X, BIBÉ B, BOUIX J, CAIMENT F, ELSEN J M, EYCHENNE F, LARZUL C, LAVILLE E MEISH F, MILENKOVIC D, TOBIN J, CHARLIER C, GEORGES M. A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep. Nature Genetics, 2006, 38(7): 813-818.

[7] ALEXANDER M S, CASAR J C, MOTOHASHI N, MYERS J A, EISENBERG I, GONZALEZ R T, ESTRELLA E A, KANG P B, KAWAHARA G, KUNKEL L M. Regulation of DMD pathology by an ankyrin-encoded miRNA. Skeletal Muscle, 2011, 1(1): 27. DOI: 10.1186/2044-5040-1-27.

[8] HITACHI K, NAKATANI M, TSUCHIDA K. Myostatin signaling regulates Akt activity via the regulation of miR-486 expression. The International Journal of Biochemistry＆Cell Biology, 2014, 47: 93-103.

[9] 张伟, 王世银, 邓双义, 杨力伟, 石国庆, 甘尚权. 巴什拜羊骨骼肌不同发育阶段差异表达miRNA研究. 农业生物技术学报, 2018, 26(1): 104-112.

ZHANG W, WANG S Y, DENG S Y, YANG L W, SHI G Q, GAN S Q. Research of differentially expressed miRNA of skeletal muscle during different development stages in Bashbay sheep (). Journal of Agircultural Biotechnology, 2018, 26(1): 104-112. (in Chinese)

[10] 张伟, 王世银, 石国庆, 邓双义, 刘晓娜, 杨力伟. 巴什拜羊miR-486多态性及其表达规律研究. 农业生物技术学报, 2020, 28(1): 92-100.

ZHANG W, WANG S Y, SHI G Q, DENG S Y, LIU X N, YANG L W. Polymorphism of miR-486 and its expression pattern in Bashbay sheep (). Journal of Agircultural Biotechnology, 2020, 28(1): 92-100. (in Chinese)

[11] 甘尚权, 高蕊, 王立民, 沈敏, 王新华, 刘守仁. 获取miRNA候选靶基因的方法及其专用反转录引物[P]. 中国, 发明专利, CN103194441A, 2013-07-10.

GAN S Q, GAO R, WANG L M, SHEN M, WANG X H, LIU S R. The method and special primers of obtaining target genes of miRNA [P]. China, Patent for Invention, CN103194441A, 2013-07-10. (in Chinese)

[12] CAMON E, MAGRANE M, BARRELL D, LEE V, DIMMER E, MASLEN J, BINNS D, HARTE N, LOPEZ R, APWEILER R. The gene ontology annotation (GOA) database: sharing knowledge in uniprot with gene ontology. Nucleic Acids Research, 2004, 32 (Database issue): D262.

[13] HATTORI M, ITOH M, ARAKI M, HIRAKAWA M, KAWASHIMA S, OKUDA S, GOTO S, KATAYAMA T, TOKIMATSU T, YAMANISHI Y, KANEHISA M. KEGG for linking genomes to life and the environment. Nucleic Acids Research, 2007, 36(Suppl.1): D480-D484.

[14] PICARD B, LEFAUCHEUR L, BERRI C, DUCLOS M J. Muscle fiber ontogenesis in farm animal species. Reproduction Nutrition Development, 2002, 42(5): 415-431.

[15] STICKLAND N C, DEMIRTAS B, CLELLAND A K, ASHTON C. Genetic and nutritional influence on muscle growth in farm animals. Comparative Biochemistry and Physiology Part A: Molecular＆Integrative Physiology, 2000, 126, 141. DOI: 10.1016/S1095- 6433(00)80279-2.

[16] ABMAYR S M, PAVLATH G K. Myoblast fusion: Lessons from flies and mice. Development, 2012, 139(4): 641-656.

[17] WILSON S J, MCEWAN J C, SHEARD P W, HARRIS A J. Early stages of myogenesis in a large mammal: formation of successive generations of myotubes in sheep tibialis muscle. Journal of Muscle Research ＆ Cell Motility, 1992, 13(9): 534-550.

[18] DRAEGER A, WEEDS A G, FITZSIMONS R B. Primary, secondary and tertiary myotubes in the developing skeletal muscle: a new approach to analysis of human myogenesis. Journal of the Neurological Sciences, 1987, 81(1): 19-43.

[19] LEFAUCHEUR L, EDOM F, ECOLAN P, BUTLER-BROWNE G S. Pattern of muscle fiber type formation in the pig. Development Dynamics, 1995, 203(1): 27-41.

[20] MASCARELLO J, STECHINI M L, ROWLERSON A, BALLOCHI E. Tertiary myotubes in postnatal growing pig muscle detected by their myosin isoform composition. Journal of Animal Science, 1992, 70(6): 1806-1813.

[21] 李雪娇, 刘晨曦, 杨开伦, 刘明军. 德美羊与中美羊胎儿期骨骼肌组织学结构发育特征差异性研究. 草食家畜, 2017(4): 1-6.

LI X J, LIU C X, YANG K L, LIU M J. Study on differentiation of fetal skeletal muscle development characteristics between German and Chinese Merino Sheep. Grass-feeding Livestock, 2017(4): 1-6. (in Chinese)

[22] 李雪娇, 刘晨曦, 孙亚伟, 杨开伦, 刘明军. 德国美利奴羊胎儿期骨骼肌组织学结构发育特征研究. 西北农林科技大学学报(自然科学版), 2018, 46(5): 7-13.

LI X J, LIU C X, SUN Y W, YANG K L, LIU M J. Study on structure development characteristics of German Merino sheep fetal skeletal muscle tissue. Journal of Northwest A ＆ F University (Natural Science Edition), 2018, 46(5): 7-13. (in Chinese)

[23] 石田培, 王欣悦, 侯浩宾, 赵志达, 尚明玉, 张莉. 基于全转录组测序的绵羊胚胎不同发育阶段骨骼肌circRNA的分析与鉴定. 中国农业科学, 2020, 53(6): 642-657.

SHI T P, WANG X Y, HOU H B, ZHAO Z D, SHANG M Y, ZHANG L. Analysis and identification of circRNAs of skeletal muscle at different stages of sheep embryos based on whole transcriptome sequencing. Scientia Agricultura Sinica, 2020, 53(6): 642-657. (in Chinese)

[24] SMALL E M, O′ROURKE J R, MORESI V, SUTHERLAND L B, MCANALLY J, GERARD R D, RICHARDSON J A, OLSON E N. Regulation of PI3-kinase/Akt signaling by muscle-enriched microRNA- 486. PNAS, 2010, 107(9): 4218-4223.

[25] PORTER N C, RESNECK W G, O’NEILL A, VAN ROSSUM D B, STONE M R, BLOCH R J. Association of small ankyrin 1 with the sarcoplasmic reticulum. Molecular Membrane Biology, 2005, 22(5): 421-432.

[26] SONG M S, SALMENA L, PANDOLFI P P. The functions and regulation of the PTEN tumour suppressor. Nature Reviews Molecular Cell Biology, 2012, 13: 283-296.

[27] YUE F, BI P P, WANG C, SHAN T Z, NIE Y H, RATLIFF T L, GAVIN T P, KUANG S H. Pten is necessary for the quiescence and maintenance of adult muscle stem cells. Nature Communications, 2017, 8, 14328. https://doi.org/10.1038/ncomms14328.

[28] HAKUNO F, YAMAUCHI Y, KANEKO G, YONEYAMA Y, NAKAE J, CHIDA K, KADOWAKI T, YAMANOUCHI K, NISHIHARA M, TAKAHASHI S. Constitutive expression of insulin receptor substrate (IRS)-1 inhibits myogenic differentiation through nuclear exclusion of Foxo1 in L6 myoblasts. PLoS One, 2011, 6: e25655.

[29] XU M, CHEN X L, CHEN D W, YU B, HUANG Z Q. FoxO1: a novel insight into its molecular mechanisms in the regulation of skeletal muscle differentiation and fiber type specification. Oncotarget, 2017, 8(6): 10662-10674.

[30] ALEXANDER M S, CASAR J C, MOTOHASHI N, VIEIRA N M, EISENBERG I, MARSHALL J L, GASPERINI M J, LEK A, MYERS J A, ESTRELLA E A, KANG P B, SHAPIRO F, RAHIMOV F, KAWAHARA G, WIDRICK J J, KUNKEL L M. MicroRNA-486- dependent modulation of DOCK3/PTEN/AKT signaling pathways improves muscular dystrophy-associated symptoms. The Journal of Clinical Investigation, 2014, 124(6): 2651-2667.

[31] XU J, LI R S, WORKENEH B, DONG Y L, WANG X N, HU Z Y. Transcription factor FoxO1, the dominant mediator of muscle wasting in chronic kidney disease, is inhibited by microRNA-486. Kidney International, 2012, 82: 401-411.

[32] WANG X H, MITCH W E. Mechanisms of muscle wasting in chronic kidney disease. Nature Reviews Nephrology, 2014, 10(9): 504-516.

[33] PONNUSAMY A, SINHA S, HYDE G D, BORLAND S J, TAYLOR R F, POND E, EYRE H J, INKSON C A, GILMORE A, ASHTON N, KALRA P A, CANFIELD A. FTI-277 inhibits smooth muscle cell

calcification by up-regulating PI3K/Akt signaling and inhibiting apoptosis. PLoS ONE, 2018, 13(4): e0196232.

[34] GLASS D. Molecular mechanisms modulating muscle mass. TRENDS in Molecular Medicine, 2003, 9(8): 344-350.

[35] DUTT V, GUPTA S, DABUR R, INJETI E, MITTAL A. Skeletal muscle atrophy: Potential therapeutic agents and their mechanisms of action. Pharmacological Research, 2015, 99: 86-100.

[36] GLASS D L. Skeletal muscle hypertrophy and atrophy signaling pathways. The International Journal of Biochemistry＆Cell Biology, 2005, 37(10): 1974-1984.

[37] STITT T N, DRUJAN D, CLARKE B A, PANARO F, TIMOFEYVA Y, KLINE W O, GONZALEZ M, YANCOPOULOS G D, GLASS D J. The IGF-1/PI3K/Akt Pathway Prevents Expression of Muscle Atrophy-Induced Ubiquitin Ligases by Inhibiting FOXO Transcription Factors. Molecular Cell, 2004, 14(3): 395-403.

[38] MARGOLIA L M, BERRYMAN C E, MURPHY N E, CARRIGAN C T, YOUNG A J, CARBONE J W, PASIAKOS S M. PI3K-AKT- FOXO1 pathway targeted by skeletal muscle microRNA to suppress proteolytic gene expression in response to carbohydrate intake during aerobic exercise. Physiological Reports, 2018, 6(23): e13931.

[39] YU Y H, CHU W L, CHAI J K, LI X, LIU L Y, MA L. Critical role of miRNAs in mediating skeletal muscle atrophy. Molecular Medicine Reports, 2016, 13: 1470-1474.

[40] GEBERT L F R, MACRAE I J. Regulation of microRNA function in animals. Nature Reviews Molecular Cell Biology, 2019, 20(1): 21-37.

[41] SEOK H, HAM J, JANG E S, CHI S W. MicroRNA target recognition: insights from transcriptome-wide non-canonical interactions. Molecules and Cells, 2016, 39(5): 375-381.

[42] ANDERSON C, CATOE H, WERNER R. miR-206 regulates connexin43 expression during skeletal muscle development. Nucleic Acids Research. 2006, 34(20): 5863-5871.

[43] GUO H, INGOLIA N T, WEISSMAN J S, BARTEL D P. Mammalian microRNAs predominantly act to decrease target mRNA levels. Nature, 2010, 466(7308): 835-841.

[44] DJURANOVIC S, NAHVI A, GREEN R. miRNA-mediated gene silencing by translational repression followed by mRNA deadenylation and decay. Science, 2012, 336(6078): 237-240.

[45] HAUSSER J, ZAVOLAN M. Identification and consequences of miRNA-target interactions-beyond repression of gene expression. Nature Reviews Genetics, 2014, 15(9): 599-612.

Investigation of miR-486 Target Genes in Skeletal Muscle of Bashbay Sheep in Different Development Periods

ZHANG Wei1, 2, WANG ShiYin1, GAO Li1, YANG LiWei1, DENG ShuangYi1, LIU XiaoNa1, SHI GuoQing2, GAN ShangQuan2

1Xinjiang Agricultural Professional Technological College, Changji 831100, Xinjiang;2State Key Laboratory of Sheep Genetic Improvement and Healthy Production/Institute of Animal Husbandry and Veterinary, Xinjiang Academy of Agricultural and Reclamation Sciences, Shihezi 832000, Xinjiang

【】The aim of this study was to deeply reveal the target genes of miR-486 in skeletal muscle of Bashbay sheep during different development periods, and to get basic data for finally uncover the molecular regulation mechanism of excellent meat traits of this sheep breed, then finally to support to further breeding. 【】The skeletal muscle of Bashbay sheep were collected during 40 d, 50 d, 60 d, 80 d, 100 d and 120 d of fetal period, and 1, 2, and 3 months old of new born sheep. The total RNA was extracted, and two cDNA libraries of miR-486 target genes were constructed using mRNA of fetal and postnatal period, respectively, then the library was sequenced applied high-throughput sequencing technology. Based on the function analysis, 10 candidate target genes were selected, and their express patterns in Bashbay sheep skeletal muscle of 10 development periods mentioned above were detected by qRT-PCR. By analyzing the expression pattern of these 10 candidate target genes and miR-486 in skeletal muscle of Bashbay sheep, their target regulation relationship and functions were primarily verified. Finally, 4 candidate target genes were chosen to confirm their regulation relationship using the double luciferase reporter assay and target regulation assay in satellite cell of skeletal muscle of Bashbay sheep. 【】Total 123 and 118 target genes of fetal and postnatal period were obtained, respectively. Due to these target genes were not confirmed by further function assay, so they were called candidate target genes temporarily. Among of them, 96 genes were expressed in skeletal muscle of fetal and postnatal period, 27 and 22 genes were specially expressed in these two periods respectively. The result of GO and KEGG analysis showed that these candidate target genes regulated mass pathways related to muscle cell differentiation and development, just like PI3k-Akt, MAPK, Wnt, Adherens junction and Regulation of actin cytoskeleton, etc. All 10 candidate target genes were expressed in skeletal muscle of Bashbay sheep, but their expression patterns were different. Among of them,,,,,,andwere expressed at a relatively high level in the fetal period skeletal muscle of Bashbay sheep, but were significantly down regulated during postnatal period, andshowed an opposite expression pattern compared with above 7 genes. The expression ofandgene did not found significantly change in all 10 development periods of Bashbay sheep. The result of double luciferase reporter assay showed that miR-486 could bind the target sites of,,andefficiently, and significantly suppressed the activity of firefly luciferase. In skeletal muscle satellite cell, mir-486 also could down-regulate mRNA of these four genes significantly, and finally suppressed their biology functions. So it was ultimately confirmed that miR-486 could regulate these 4 target genes indeed in skeletal muscle of Bashbay sheep. 【】The investigation deeply revealed the target genes of miR-486 in skeletal muscle of Bashbay sheep during different development periods for the first time, and comprehensively analyzed the biology process and signaling pathway they participated. The further function assay proved that the data of target genes were reliable. These data would help to uncover the molecular regulation mechanisms related to excellent meat traits of Bashbay sheep.

Bashbay sheep; miR-486; target gene; skeletal muscle; cDNA library