Systematic Identification and Comparison of the Expressed Profiles of lncRNAs, miRNAs, circRNAs, and mRNAs with Associated Co-Expression Networks in Pigs with Low and High Intramuscular Fat

Abstract

Intramuscular fat (IMF) content is a complex trait that affects meat quality and determines pork quality. In order to explore the potential mechanisms that affect the intramuscular fat content of pigs, a Large white × Min pigs F2 resource populations were constructed, then whole-transcriptome profile analysis was carried out for five low-IMF and five high-IMF F2 individuals. In total, 218 messenger RNA (mRNAs), 213 long non-coding RNAs (lncRNAs), 18 microRNAs (miRNAs), and 59 circular RNAs (circRNAs) were found to be differentially expressed in the longissimus dorsi muscle. Gene ontology analysis and Kyoto Encyclopedia of Genes and Genomes annotations revealed that these differentially expressed (DE) genes or potential target genes (PTGs) of DE regulatory RNAs (lncRNAs, miRNAs, and circRNAs) are mainly involved in cell differentiation, fatty acid synthesis, system development, muscle fiber development, and regulating lipid metabolism. In total, 274 PTGs were found to be differentially expressed between low- and high-IMF pigs, which indicated that some DE regulatory RNAs may contribute to the deposition/metabolism of IMF by regulating their PTGs. In addition, we analyzed the quantitative trait loci (QTLs) of DE RNAs co-located in high- and low-IMF groups. A total of 97 DE regulatory RNAs could be found located in the QTLs related to IMF. Co-expression networks among different types of RNA and competing endogenous RNA (ceRNA) regulatory networks were also constructed, and some genes involved in type I diabetes mellitus were found to play an important role in the complex molecular process of intramuscular fat deposition. This study identified and analyzed some differential RNAs, regulatory RNAs, and PTGs related to IMF, and provided new insights into the study of IMF formation at the level of the genome-wide landscape.