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

Whole-Transcriptome Analysis Sheds Light on the Biological Contexts of Intramuscular Fat Deposition in Ningxiang Pigs

The quality of pork is significantly impacted by intramuscular fat (IMF). However, the regulatory mechanism of IMF depositions remains unclear. We performed whole-transcriptome sequencing of the longissimus dorsi muscle (IMF) from the high (5.1 ± 0.08) and low (2.9 ± 0.51) IMF groups (%) to elucidate potential mechanisms. In summary, 285 differentially expressed genes (DEGs), 14 differentially expressed miRNAs (DEMIs), 83 differentially expressed lncRNAs (DELs), and 79 differentially expressed circRNAs (DECs) were identified. DEGs were widely associated with IMF deposition and liposome differentiation. Furthermore, competing endogenous RNA (ceRNA) regulatory networks were constructed through co-differential expression analyses, which included circRNA-miRNA-mRNA (containing 6 DEMIs, 6 DEGs, 47 DECs) and lncRNA-miRNA-mRNA (containing 6 DEMIs, 6 DEGs, 36 DELs) regulatory networks. The circRNAs sus-TRPM7_0005, sus-MTUS1_0004, the lncRNAs SMSTRG.4269.1, and MSTRG.7983.2 regulate the expression of six lipid metabolism-related target genes, including PLCB1, BAD, and GADD45G, through the binding sites of 2-4068, miR-7134-3p, and miR-190a. For instance, MSTRG.4269.1 regulates its targets PLCB1 and BAD via miRNA 2_4068. Meanwhile, sus-TRPM7_0005 controls its target LRP5 through ssc-miR-7134-3P. These findings indicate molecular regulatory networks that could potentially be applied for the marker-assisted selection of IMF to enhance pork quality.

Transcriptome Analysis Reveals the Age-Related Developmental Dynamics Pattern of the Longissimus Dorsi Muscle in Ningxiang Pigs

The growth and development of the Longissimus Dorsi muscle are complex, playing an important role in the determination of pork quality. The study of the Longissimus Dorsi muscle at the mRNA level is particularly crucial for finding molecular approaches to improving meat quality in pig breeding. The current study utilized transcriptome technology to explore the regulatory mechanisms of muscle growth and intramuscular fat (IMF) deposition in the Longissimus Dorsi muscle at three core developmental stages (natal stage on day 1, growing stage on day 60, and finishing stage on day 210) in Ningxiang pigs. Our results revealed 441 differentially expressed genes (DEGs) in common for day 1 vs. day 60 and day 60 vs. day 210, and GO (Gene Ontology) analysis showed that candidate genes RIPOR2, MEGF10, KLHL40, PLEC, TBX3, FBP2, and HOMER1 may be closely related to muscle growth and development, while KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis showed that DEGs (UBC, SLC27A5, RXRG, PRKCQ, PRKAG2, PPARGC1A, PLIN5, PLIN4, IRS2, and CPT1B) involved the PPAR (Peroxisome Proliferator-Activated Receptor) signaling pathway and adipocytokine signaling pathway, which might play a pivotal role in the regulation of IMF deposition. PPI (Protein-Protein Interaction Networks) analysis found that the STAT1 gene was the top hub gene. Taken together, our results provide evidence for the molecular mechanisms of growth and development and IMF deposition in Longissimus Dorsi muscle to optimize carcass mass.

Identification of transcriptional regulatory variants in pig duodenum, liver, and muscle tissues

BACKGROUND

In humans and livestock species, genome-wide association studies (GWAS) have been applied to study the association between variants distributed across the genome and a phenotype of interest. To discover genetic polymorphisms affecting the duodenum, liver, and muscle transcriptomes of 300 pigs from 3 different breeds (Duroc, Landrace, and Large White), we performed expression GWAS between 25,315,878 polymorphisms and the expression of 13,891 genes in duodenum, 12,748 genes in liver, and 11,617 genes in muscle.

RESULTS

More than 9.68 × 1011 association tests were performed, yielding 14,096,080 significantly associated variants, which were grouped in 26,414 expression quantitative trait locus (eQTL) regions. Over 56% of the variants were within 1 Mb of their associated gene. In addition to the 100-kb region upstream of the transcription start site, we identified the importance of the 100-kb region downstream of the 3'UTR for gene regulation, as most of the cis-regulatory variants were located within these 2 regions. We also observed 39,874 hotspot regulatory polymorphisms associated with the expression of 10 or more genes that could modify the protein structure or the expression of a regulator gene. In addition, 2 motifs (5'-GATCCNGYGTTGCYG-3' and a poly(A) sequence) were enriched across the 3 tissues within the neighboring sequences of the most significant single-nucleotide polymorphisms in each cis-eQTL region.

CONCLUSIONS

The 14 million significant associations obtained in this study are publicly available and have enabled the identification of expression-associated cis-, trans-, and hotspot regulatory variants within and across tissues, thus shedding light on the molecular mechanisms of regulatory variations that shape end-trait phenotypes.

miR-381-3p Inhibits Intramuscular Fat Deposition through Targeting FABP3 by ceRNA Regulatory Network

Intramuscular fat (IMF) deposition is an important determinant of pork quality and a complex process facilitated by non-coding ceRNAs. In this study, 52 Berkshire × Anqing Sixwhite crossbred pigs were slaughtered to measure eight carcass and pork quality traits. Whole-transcriptome sequencing analysis was performed using longissimus dorsi samples of six low- and high-IMF samples; 34 ceRNA networks, based on 881, 394, 158 differentially expressed (DE) lncRNAs, miRNAs, and mRNAs, were constructed. Following weighted gene co-expression network analysis between the low and high IMF, only one ceRNA, lncRNA4789/miR-381-3p/FABP3, that showed similar DE trend in longissimus dorsi tissue was retained. Dual-luciferase reporter assays further indicated that FABP3 was a direct, functional target of miR-381-3p, where miR-381-3p overexpression inhibited the mRNA and protein expression of FABP3. In addition, overexpressed lncRNA4789 attenuated the effect of miR-381-3p on FABP3 by sponging miR-381-3p. Cell function verification experiment demonstrated that miR-381-3p suppressed IMF deposition by inhibiting preadipocyte cell differentiation and lipid droplet deposition via the suppression of FABP3 expression in the peroxisome proliferator-activated receptor signalling pathway, whereas lncRNA4789 rescued FABP3 expression by sponging miR-381-3p. Our study may aid in identifying novel molecular markers for its optimization in IMF which is of importance in breeding for improving pork quality.

Genome-wide characterization of lncRNAs and mRNAs in muscles with differential intramuscular fat contents

Meat quality is one of the most important traits in pig production. Long non-coding RNAs (lncRNAs) have been involved in diverse biological processes such as muscle development through regulating gene expression. However, studies on lncRNAs lag behind and a comparatively small number of lncRNAs have been identified in pigs. Also, the effects of lncRNAs on meat quality remain to be characterized. Here, we analyzed lncRNAs in longissimus thoracis (LT) and semitendinosus (ST) muscles, being different in meat quality, with RNA-sequencing technology. A total of 500 differentially expressed lncRNAs (DELs) and 2,094 protein-coding genes (DEGs) were identified. Through KEGG analysis on DELs, we first made clear that fat deposition might be the main reason resulting in the differential phenotype of LT and ST, for which cGMP–PKG and VEGF signaling pathways were the most important ones. In total, forty-one key DELs and 50 DEGs involved in the differential fat deposition were then characterized. One of the key genes, cAMP-response element binding protein 1, was selected to confirm its role in porcine adipogenesis with molecular biology methods and found that it promotes the differentiation of porcine preadipocytes, consistent with its higher expression level and intramuscular fat contents in LT than that in ST muscle. Furthermore, through integrated analysis of DELs and DEGs, transcription factors important for differential fat deposition were characterized among which BCL6 has the most target DEGs while MEF2A was targeted by the most DELs. The results provide candidate genes crucial for meat quality, which will contribute to improving meat quality with molecular-breeding strategies.

Circular RNA Profiling Identifies Novel circPPARA that Promotes Intramuscular Fat Deposition in Pigs

Intramuscular fat (IMF) content plays an important role in pork quality. Circular RNAs (circRNAs) implicate various biological processes; however, the regulatory mechanisms and functions of circRNAs in porcine IMF remains elusive. Hence, the study assessed the circRNA expression profiling in the longissimus dorsi muscle of pigs with high (H) and low (L) IMF content to unravel their regulatory functions in improving meat quality. The RNA sequencing analysis identified 29,732 circRNAs from six sampled pigs, most of which were exon-derived. In the muscle, 336 were differentially expressed (DE) between the H and L IMF groups; 196 circRNAs were upregulated, and 140 were downregulated. Subsequent qRT-PCR validation of 10 DE circRNAs revealed expression patterns consistent with the RNA-seq data. Gene ontology and KEGG enrichment analysis revealed that most significantly enriched DE circRNAs' host genes were linked to lipid metabolism and adipogenesis processes. The circRNA-miRNA regulatory network analysis found several circRNAs targeting miRNAs associated with adipogenesis. Finally, a novel circRNA, circPPARA, was identified with the expression positively correlated with the IMF content. Detailed analysis revealed that circPPARA was formed via head-to-tail splicing and was more stable than the linear PPARA, predominantly located in the cytoplasm. Functional studies using overexpression and siRNA constructs demonstrated that circPPARA promotes differentiation and hinders the proliferation of porcine intramuscular preadipocytes. Moreover, the dual-luciferase assay revealed that circPPARA adsorbed miR-429 and miR-200b, thereby promoting intramuscular adipogenesis in pigs. Our results identified a candidate circRNA, circPPARA, that affects porcine IMF content. The study provides knowledge of the regulatory functions of circRNAs in intramuscular adipogenesis and abundant resource for future research on circRNAs in pigs.