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

doi:10.16431/j.cnki.1671-7236.2021.06.002

摘要/Abstract

摘要： 试验旨在分析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的表达差异，研究结果为阐明德州驴骨骼肌生长发育的分子调控机制提供参考。

关键词: 德州驴; 骨骼肌发育; 转录组; miRNA

Abstract: To investigate the expression of microRNAs(miRNAs) in skeletal muscles of Dezhou donkeys and further improve the regulatory mechanism of small RNA in the development of skeletal muscle of donkey.In this study,the longissimus dorsi muscle of two Dezhou donkey strains (Wutou donkey and Sanfen donkey) was selected as the research objects,by constructing small RNA sequencing library of Wutou donkey longissimus dorsi muscle (WB) and Sanfen donkey longissimus dorsi muscle (SB),bioinformatics techniques were used to carry out differential expression analysis of miRNAs between WB and SB,and predict and analyze miRNAs that might affect the meat production traits of donkey.In addition,six miRNAs with high expression were selected,and Real-time quantitative PCR was used to verify the accuracy of sequencing results.The results showed that 982 and 1 149 conserved miRNAs and 106 and 143 new miRNAs were identified in WB and SB libraries.In WB and SB comparison group (WB vs SB) based on conditions|log₂(fold_change)|≥2,P≤0.05 get 17 differentially expressed miRNAs,of which 5 were up-regulated and 12 were down-regulated.Through GO annotation,it found that differentially expressed miRNAs were significantly enriched in biological processes such as regulating cell migration and differentiation,cell division and skeletal muscle cell contraction (P<0.05).bta-miR-378d_L+1R-1_1ss22CA and bta-miR-378d_R-2 which related to skeletal muscle development in donkey was found,they target multiple genes such as NPY1R,GPR152,BEST1,KCNH8 and SPAG9.KEGG enrichment results showed that miRNAs were enriched in Wnt,MAPK,ubiquitin proteolysis system and other signaling pathways related to skeletal muscle development.Among them,a total of 207 target genes 17 differentially expressed miRNAs were involved in the Wnt signaling pathway,especially PC-5p-84483_31,hsa-miR-186-3p_R-3 and cfa-miR-329b_1ss19CT had the largest number of target genes.The results of Real-time quantitative PCR were consistent with the results of miRNAs sequencing,and it found that the expression levels of the other five miRNAs were all higher in the SB than in the WB except for eca-miR-181b.In this study,miRNAs of skeletal muscle of two strains of Dezhou donkey were analyzed for the first time,and the differences of RNA-seq between Wutou donkey and Sanfen donkey were revealed,which provided reference for elucidating the molecular regulation mechanism of skeletal muscle growth and development of Dezhou donkey.

Key words: Dezhou donkey; skeletal muscle development; transcriptome; miRNA

中图分类号:

S822

马程, 胡文萍, 石田培, 侯浩宾, 李海静, 张新浩, 杨莉, 张莉. 乌头驴与三粉驴骨骼肌中miRNAs的表达鉴定与分析[J]. 中国畜牧兽医, 2021, 48(6): 1894-1905.

MA Cheng, HU Wenping, SHI Tianpei, HOU Haobin, LI Haijing, ZHANG Xinhao, YANG Li, ZHANG Li. Identification and Analysis of miRNAs Expression in Skeletal Muscle of Wutou Donkey and Sanfen Donkey[J]. China Animal Husbandry and Veterinary Medicine, 2021, 48(6): 1894-1905.

参考文献

[1] CALLIS T E,CHEN J F,WANG D Z.MicroRNAs in skeletal and cardiac muscle development[J].DNA and Cell Biology,2007,26(4):219-225.
[2] JIA B,LIU Y,LI Q,et al.Altered miRNA and mRNA expression in sika deer skeletal muscle with age[J].Genes,2020,11(2):172.
[3] LEWIS B P,BURGE C B,BARTEL D P.Conserved seed pairing,often flanked by adenosines,indicates that thousands of human genes are microRNA targets[J].Cell,2005,120(1):15-20.
[4] HOU X H,YANG Y L,ZHU S Y,et al.Comparison of skeletal muscle miRNA and mRNA profiles among three pig breeds[J].Molecular Genetics and Genomics,2016,291(2):559-573.
[5] MIAO Z G,WANG S,WANG Y M,et al.Comparison of microRNAs in the intramuscular adipose tissue from Jinhua and Landrace pigs[J].Journal of Cellular Biochemistry,2019,120(1):192-200.
[6] HORAK M,NOVAK J,BIENERTOVA-VASKU J.Muscle-specific microRNAs in skeletal muscle deve-lopment[J].Developmental Biology,2015,410(1):1-13.
[7] GUI Y,LI D,CHEN J,et al.Soluble epoxide hydrolase inhibitors,t-AUCB,downregulated miR-133 in a mouse model of myocardial infarction[J].Lipids in Health & Disease,2018,17(1):129.
[8] CHEN X,ZHAO C,DOU M,et al.Deciphering the miRNA transcriptome of Rongchang pig longissimus dorsi at weaning and slaughter time points[J].Journal of Animal Physiology and Animal Nutrition,2020,104(3):1-11.
[9] ZHAO Q,KANG Y,WANG H Y,et al.Expression profiling and functional characterization of miR-192 throughout sheep skeletal muscle development[J].Scientific Reports,2016,6(1):1-12.
[10] MCGIVNEY B A,GRIFFIN M E,GOUGH K F,et al.Evaluation of microRNA expression in plasma and skeletal muscle of thoroughbred racehorses in training[J].BMC Veterinary Research,2017,13(1):347.
[11] 侯浩宾.德州驴NCAPG,DCAF16,FAM184B,TBX3基因多态性与早期生长性状的关联分析[D].北京:中国农业科学院,2019. HOU H B.Association analysis of NCAPG,DCAF16,FAM184B,TBX3 gene polymorphism and early growth traits in Dezhou donkey[D].Beijing:Chinese Academy of Agricultural Sciences,2019.(in Chinese)
[12] ZHOU X L,XIAO G L.Supplementation of high protein and high fat concentrates affects milk yield and milk composition of donkeys(Equus asinus)[J].Chinese Journal of Animal Nutrition,2012,24(3):577-582.
[13] DENG L,DUAN H,ZHANG X,et al.Advances in the research and application of artificial insemination to equids in China:1935-2012[J].Journal of Equine Veterinary Science,2014,34(3):351-359.
[14] SUN Y,JIANG Q,YANG C,et al.Characterization of complete mitochondrial genome of Dezhou donkey (Equus asinus) and evolutionary analysis[J].Current Genetics,2016,62(2):383-390.
[15] WANG Z,BIAN L,MO C,et al.Targeted quanti-fication of lipid mediators in skeletal muscles using restricted access media-based trap-and-elute liquid chromatography-mass spectrometry[J].Analytica Chimica Acta,2017,984:151-161.
[16] BO J L,PING H L,WANG J W,et al.Progresses in research of the mechanisms of skeletal muscle fiber formation[J].Scientia Agricultura Sinica,2014,47(6):1200-1207.
[17] GALIMOV A,MERRY B T L,LUCA C E,et al.microRNA-29a in adult muscle stem cells controls skeletal muscle regeneration during injury and exercise downstream of fibroblast growth factor-2[J].Stem Cells,2016,34(3):768-780.
[18] SHI T P,HU W P,HOU H B,et al.Identification and comparative analysis of long non-coding RNA in the skeletal muscle of two Dezhou donkey strains[J].Genes,2020,11(5):508.
[19] 孙伟,唐中林,谭林,等.HMGCS1基因对猪肌肉生长发育的影响及miRNA-18a/b对HMGCS1基因的调控[J].中国农业科学,2013,46(12):2543-2549. SUN W,TANG Z L,TAN L,et al.Effect of HMGCS1 gene on muscle growth and development of pigs and regulation of miRNA-18a/b on HMGCS1 gene[J].Chinese Journal of Agricultural Sciences,2013,46(12):2543-2549.(in Chinese)
[20] FILOMENA N D.Epigenetic regulators:Polycomb-miRNA circuits in cancer[J].Biochimica et Biophysica Acta,2016,1859(5):697-704.
[21] CONNOLLY M,GARFIELD B E,CROSBY A,et al.miR-1-5p targets TGF-βR1 and is suppressed in the hypertrophying hearts of rats with pulmonary arterial hypertension[J].PLoS One,2020,15(2):1-10.
[22] ZHOU B,LIU H L,SHI F X,et al.microRNA expression profiles of porcine skeletal muscle[J].Animal Genetics,2010,41(5):499-508.
[23] JI H,WANG H,JI Q M,et al.Differential expression profile of microRNA in yak skeletal muscle and adipose tissue during development[J].Genes & Genomics,2020,42(11):1347-1359.
[24] WUST S,DROSE S,HEIDLER J,et al.Metabolic maturation during muscle stem cell differentiation is achieved by miR-1/133a-mediated inhibition of the Dlk1-Dio3 mega gene cluster[J].Cell Metabolism,2018,27(5):1026-1039.
[25] CHEN X,HE L,ZHAO Y,et al.Malat1 regulates myogenic differentiation and muscle regeneration through modulating MyoD transcriptional activity[J].Cell Discovery,2017,3:1-23.
[26] NINSUWON J,WAIYAMITRA P,ROONGSITTHICHAI A,et al.Expressions of miR-155 and miR-181 and predictions of their structures and targets in pigs (Sus scrofa)[J].Veterinary World,2020,13(8):1667-1673.
[27] WEI X,LI H,ZHANG B,et al.miR-378a-3p promotes differentiation and inhibits proliferation of myoblasts by targeting HDAC4 in skeletal muscle development[J].RNA Biology,2016,13(12):1300-1309.
[28] WANG B,ZHANG A,WANG H,et al.miR-26a limits muscle wasting and cardiac fibrosis through exosome-mediated microRNA transfer in chronic kidney disease[J].Theranostics,2019,9(7):1864-1877.
[29] CHEN M H,WANG J Y,WANG Y P,et al.Genome-wide detection of selection signatures in Chinese indigenous Laiwu pigs revealed candidate genes regulating fat deposition in muscle[J].BMC Genetics,2018,19(1):31.
[30] SCHWARTZ A L.Ubiquitin-mediated proteolysis:Biological regulation via destruction[J].Bioassays,2015,22(5):442-451.
[31] LAWAN A,MIN K,ZHANG L,et al.Keletal muscle-specific deletion of MKP-1 reveals a p38MAPK/JNK/Akt signaling node that regulates obesity-induced insulin resistance[J].Diabetes,2018,67(4):624-635.
[32] HAN S,CUI C,HE H,et al.Myoferlin regulates Wnt/β-catenin signaling-mediated skeletal muscle development by stabilizing dishevelled-2 against autophagy[J].International Journal of Molecular Sciences,2019,20(20):5130.
[33] LOGAN C Y,NUSSE R.The Wnt signaling pathway in development and disease[J].Annual Review of Cell and Developmental Biology,2004,20(1):781-810.
[34] LACOUR F,VEZIN E,BENTZINGER F C,et al.R-spondin1 regulates muscle progenitor cell differentiation and fusion through dual regulation of canonical and non-canonical Wnt signaling pathways[J].Cell Reports,2017,18(10):2320-2330.
[35] 王建清.小尾寒羊miR-221和miR-329b-3p功能研究及其与IRS1基因的靶向关系验证[D].兰州:甘肃农业大学,2020. WANG J Q.Functional studies of miR-221 and miR-329b-3P in Small-tailed Han sheep and validation of their targeted relationship with IRS1 gene[D].Lanzhou:Gansu Agricultural University,2020.(in Chinese)
[36] HE M,JIN Q,CHEN C,et al.The miR-186-3p/EREG axis orchestrates tamoxifen resistance and aerobic glycolysis in breast cancer cells[J].Oncogene,2019,38(28):5551-5565.
[37] YANG R,WEI M J,YANG F,et al.Diosbulbin B induced G2/M cell cycle arrest in hepatocytes by miRNA-186-3p and miRNA-378a-5p-mediated the decreased expression of CDK1[J].Toxicology & Applied Pharmacology,2018,357:1-9.