Identification of Candidate Circular RNAs Underlying Intramuscular Fat Content in the Donkey

Intramuscular fat (IMF) content is a crucial indicator of meat quality. Circular RNAs (circRNAs) are a large class of endogenous RNAs that are involved in many physiological processes. However, the expression and function of circRNA in IMF in the donkey remains unresolved. Here we performed an expression profiling of circRNAs in the donkey longissimus dorsi muscle and identified 12,727 candidate circRNAs. Among these, 70% were derived from the exons of protein genes. Furthermore, a total of 127 differentially expressed (DE) circRNAs were identified in high (H) and low (L) IMF content groups, including 63 upregulated and 64 downregulated circRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of the host genes of the DE circRNAs showed that the host genes were enriched in lipid metabolism related GO terms (e.g., fatty acid beta-oxidation using acyl-CoA dehydrogenase and MLL3/4 complex), and signaling pathways (e.g., TGF-beta and lysine degradation signaling pathway). Further analyses indicated that 127 DE circRNAs were predicted to potentially interact with miRNAs, leading to the construction of circRNA-miRNA regulatory network. Multiple circRNAs can potentially function as sponges of miRNAs that regulate the differentiation of adipocytes. Our results provide valuable expression profile information for circRNA in the donkey and new insight into the regulatory mechanisms of circRNAs in the regulation of IMF content.