Landscape of tissue-specific RNA Editome provides insight into co-regulated and altered gene expression in pigs (Sus-scrofa)

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