RNA editing generates genetic diversity in mammals by altering amino acid sequences, miRNA targeting site sequences, influencing the stability of targeted RNAs, and causing changes in gene expression. However, the extent to which RNA editing affect gene expression via modifying miRNA binding site remains unexplored. Here, we first profiled the dynamic A-to-I RNA editome across tissues of Duroc and Luchuan pigs. The RNA editing events at the miRNA binding sites were generated. The biological function of the differentially edited gene in skeletal muscle was further characterized in pig muscle-derived satellite cells. RNA editome analysis revealed a total of 171,909 A-to-I RNA editing sites (RESs), and examination of its features showed that these A-to-I editing sites were mainly located in SINE retrotransposons PRE-1/Pre0_SS element. Analysis of differentially edited sites (DESs) revealed a total of 4,552 DESs across tissues between Duroc and Luchuan pigs, and functional category enrichment analysis of differentially edited gene (DEG) sets highlighted a significant association and enrichment of tissue-developmental pathways including TGF-beta, PI3K-Akt, AMPK, and Wnt signaling pathways. Moreover, we found that RNA editing events at the miRNA binding sites in the 3′-UTR of HSPA12B mRNA could prevent the miRNA-mediated mRNA downregulation of HSPA12B in the muscle-derived satellite (MDS) cell, consistent with the results obtained from the Luchuan skeletal muscle. This study represents the most systematic attempt to characterize the significance of RNA editing in regulating gene expression, particularly in skeletal muscle, constituting a new layer of regulation to understand the genetic mechanisms behind phenotype variance in animals. Abbreviations: A-to-I: Adenosine-to-inosine; ADAR: Adenosine deaminase acting on RNA; RES: RNA editing site; DEG: Differentially edited gene; DES: Differentially edited site; FDR: False discovery rate; GO: Gene Ontology; KEGG: Kyoto Encyclopaedia of Genes and Genomes; MDS cell: musclederived satellite cell; RPKM: Reads per kilobase of exon model in a gene per million mapped reads; UTR: Untranslated coding regions

The microRNAs (miRNAs) play an important role in regulating myogenesis by targeting mRNA. However, the understanding of miRNAs in skeletal muscle development and diseases is unclear. In this study, we firstly performed the transcriptome profiling in differentiating C2C12 myoblast cells. Totally, we identified 187 miRNAs and 4260 mRNAs significantly differentially expressed that were involved in myoblast differentiation. We carried out validation of microarray data based on 5 mRNAs and 5 miRNAs differentially expressed and got a consistent result. Then we constructed and validated the significantly up- and down-regulated mRNA-miRNA interaction networks. Four interaction pairs (miR-145a-5p-Fscn1, miR-200c-5p-Tmigd1, miR-27a-5p-Sln and miR-743a-5p-Mob1b) with targeted relationships in differentiated myoblast cells were demonstrated. They are all closely related to myoblast development. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated cell cycle signals important for exploring skeletal muscle development and disease. Functionally, we discovered that miR-743a targeting gene Mps One Binder Kinase Activator-Like 1B (Mob1b) gene in differentiated C2C12. The up-regulated miR-743a can promote the differentiation of C2C12 myoblast. While the down-regulated Mob1b plays a negative role in differentiation. In addition, the expression profile of miR-743a and Mob1b are consistent with skeletal muscle recovery after Cardiotoxin (CTX) injury. Our study revealed that miR-743a-5p regulates myoblast differentiation by targeting Mob1b involved in skeletal muscle development and regeneration. Our findings made a further exploration for mechanisms in myogenesis and might provide potential possible miRNA-based target therapies for skeletal muscle regeneration and disease in the near future.