Genome-Wide Expression Profiling of mRNAs, lncRNAs and circRNAs in Skeletal Muscle of Two Different Pig Breeds

Construction of LncRNA-Related ceRNA Networks in Longissimus Dorsi Muscle of Jinfen White Pigs at Different Developmental Stages

The development of skeletal muscle in pigs might determine the quality of pork. In recent years, long non-coding RNAs (lncRNAs) have been found to play an important role in skeletal muscle growth and development. In this study, we investigated the whole transcriptome of the longissimus dorsi muscle (LDM) of Jinfen White pigs at three developmental stages (1, 90, and 180 days) and performed a comprehensive analysis of lncRNAs, mRNAs, and micro-RNAs (miRNAs), aiming to find the key regulators and interaction networks in Jinfen White pigs. A total of 2638 differentially expressed mRNAs (DE mRNAs) and 982 differentially expressed lncRNAs (DE lncRNAs) were identified. Compared with JFW_1d, there were 497 up-regulated and 698 down-regulated DE mRNAs and 212 up-regulated and 286 down-regulated DE lncRNAs in JFW_90d, respectively. In JFW_180d, there were 613 up-regulated and 895 down-regulated DE mRNAs and 184 up-regulated and 131 down-regulated DE lncRNAs compared with JFW_1d. There were 615 up-regulated and 477 down-regulated DE mRNAs and 254 up-regulated and 355 down-regulated DE lncRNAs in JFW_180d compared with JFW_90d. Compared with mRNA, lncRNA has fewer exons, fewer ORFs, and a shorter length. We performed GO and KEGG pathway functional enrichment analysis for DE mRNAs and the potential target genes of DE lncRNAs. As a result, several pathways are involved in muscle growth and development, such as the PI3K-Akt, MAPK, hedgehog, and hippo signaling pathways. These are among the pathways through which mRNA and lncRNAs function. As part of this study, bioinformatic screening was used to identify miRNAs and DE lncRNAs that could act as ceRNAs. Finally, we constructed an lncRNA-miRNA-mRNA regulation network containing 26 mRNAs, 7 miRNAs, and 17 lncRNAs; qRT-PCR was used to verify the key genes in these networks. Among these, XLOC_022984/miR-127/ENAH and XLOC_016847/miR-486/NRF1 may function as key ceRNA networks. In this study, we obtained transcriptomic profiles from the LDM of Jinfen White pigs at three developmental stages and screened out lncRNA-miRNA-mRNA regulatory networks that may provide crucial information for the further exploration of the molecular mechanisms during skeletal muscle development.

The emerging role of circRNAs on skeletal muscle development in economical animals

CircRNAs are a novel type of closed circular molecules formed through a covalent bond lacking a 5'cap and 3' end tail, which mainly arise from mRNA precursor. They are widely distributed in plants and animals and are characterized by stable structure, high conservativeness in cells or tissues, and showed the expression specificity at different stages of development in different tissues. CircRNAs have been gradually attracted wide attention with the development of RNA sequencing, which become a new research hotspot in the field of RNA. CircRNAs play an important role in gene expression regulation. Presently, the related circRNAs research in the regulation of animal muscle development is still at the initial stage. In this review, the formation, properties, biological functions of circRNAs were summarized. The recent research progresses of circRNAs in skeletal muscle growth and development from economic animals including livestock, poultry and fishes were introduced. Finally, we proposed a prospective for further studies of circRNAs in muscle development, and we hope our research could provide new ideas, some theoretical supports and helps for new molecular genetic markers exploitation and animal genetic breeding in future.

Integrated transcriptome analysis reveals roles of long non-coding RNAs (lncRNAs) in caprine skeletal muscle mass and meat quality

Long non-coding RNAs (lncRNAs) play important roles in the growth and development of skeletal muscle. However, there is limited information on goats. In this study, expression profiles of lncRNAs in Longissimus dorsi muscle from Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats with divergent meat yield and meat quality were compared using RNA-sequencing. Based on our previous microRNA (miRNA) and mRNA profiles obtained from the same tissues, the target genes and binding miRNAs of differentially expressed lncRNAs were obtained. Subsequently, lncRNA-mRNA interaction networks and a ceRNA network of lncRNA-miRNA-mRNA were constructed. A total of 136 differentially expressed lncRNAs were identified between the two breeds. Fifteen cis target genes and 143 trans target genes were found for differentially expressed lncRNAs, and they were enriched in muscle contraction, muscle system process, muscle cell differentiation, and p53 signaling pathway. A total of 69 lncRNA-trans target gene pairs were constructed, with close relationship with muscle development, intramuscular fat deposition, and meat tenderness. A total of 16 lncRNA-miRNA-mRNA ceRNA pairs were identified, of which some reportedly associated with skeletal muscle development and fat deposition were found. The study will provide an improved understanding of the roles of lncRNAs in caprine meat yield and meat quality.

Integrated Analysis of Long Non-Coding RNA and mRNA to Reveal Putative Candidate Genes Associated with Backfat Quality in Beijing Black Pig

Pigs' backfat quality has an important impact on the quality of pork and pork products and has a strong relationship with nutrition and sensory characteristics. This study aimed to identify the related candidate genes of backfat quality and to preliminary clarify the molecular regulatory mechanism underlying pig backfat quality phenotypes. Expression assessments of long non-coding RNA (lncRNA) and mRNA profiling in backfat from high-quality (firm and white) and low-quality (soft and yellow) Beijing Black pigs were performed by RNA sequencing. Significantly different expressions were observed in 610 protein-coding genes and 290 lncRNAs between the two groups. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway annotation showed that some candidate differentially expressed genes that participate in lipid-related pathways and pigmentation terms may play a role in backfat quality in pigs. The cis-target and trans-target genes were predicted to explore the regulatory function of lncRNAs, and integrative analyses of different expression lncRNAs targets and different expression genes were performed. The results showed the regulatory networks of lncRNA-mRNA related to backfat quality, and our study obtained strong candidate genes for backfat quality: ELOVL5, SCD, DGAT2, SLC24A5, and TYRP1, which were involved in fat metabolism, adipogenesis regulation, and pigmentation. To our knowledge, this study is the first to demonstrate the molecular genetic mechanisms of backfat quality in pigs, and these findings improve the current understanding of backfat quality mechanisms and provide a foundation for further studies.

Comparative Transcriptomic Analysis of mRNAs, miRNAs and lncRNAs in the Longissimus dorsi Muscles between Fat-Type and Lean-Type Pigs

In pigs, meat quality and production are two important traits affecting the pig industry and human health. Compared to lean-type pigs, fat-type pigs contain higher intramuscular fat (IMF) contents, better taste and nutritional value. To uncover genetic factors controlling differences related to IMF in pig muscle, we performed RNA-seq analysis on the transcriptomes of the Longissimus dorsi (LD) muscle of Laiwu pigs (LW, fat-type pigs) and commercial Duroc × Landrace × Yorkshire pigs (DLY, lean-type pigs) at 150 d to compare the expression profiles of mRNA, miRNA and lncRNA. A total of 225 mRNAs, 12 miRNAs and 57 lncRNAs were found to be differentially expressed at the criteria of |log2(foldchange)| > 1 and q < 0.05. The mRNA expression of LDHB was significantly higher in the LD muscle of LW compared to DLY pigs with log2(foldchange) being 9.66. Using protein interaction prediction method, we identified more interactions of estrogen-related receptor alpha (ESRRA) associated with upregulated mRNAs, whereas versican (VCAN) and proenkephalin (PENK) were associated with downregulated mRNAs in LW pigs. Integrated analysis on differentially expressed (DE) mRNAs and miRNAs in the LD muscle between LW and DLY pigs revealed two network modules: between five upregulated mRNA genes (GALNT15, FKBP5, PPARGC1A, LOC110258214 and LOC110258215) and six downregulated miRNA genes (ssc-let-7a, ssc-miR190-3p, ssc-miR356-5p, ssc-miR573-5p, ssc-miR204-5p and ssc-miR-10383), and between three downregulated DE mRNA genes (IFRD1, LOC110258600 and LOC102158401) and six upregulated DE miRNA genes (ssc-miR1379-3p, ssc-miR1379-5p, ssc-miR397-5p, ssc-miR1358-5p, ssc-miR299-5p and ssc-miR1156-5p) in LW pigs. Based on the mRNA and ncRNA binding site targeting database, we constructed a regulatory network with miRNA as the center and mRNA and lncRNA as the target genes, including GALNT15/ssc-let-7a/LOC100523888, IFRD1/ssc-miR1379-5p/CD99, etc., forming a ceRNA network in the LD muscles that are differentially expressed between LW and DLY pigs. Collectively, these data may provide resources for further investigation of molecular mechanisms underlying differences in meat traits between lean- and fat-type pigs.

Population Genetic Structure and Selection Signature Analysis of Beijing Black Pig

Beijing Black pig is an excellent cultivated black pig breed in China, with desirable body shape, tender meat quality and robust disease resistance. To explore the level of admixture and selection signatures of Beijing Black pigs, a total number of 90 individuals covering nine pig breeds were used in our study, including Beijing Black pig, Large White, Landrace, Duroc, Lantang pig, Luchuan pig, Mashen pig, Huainan pig and Min pig. These animals were resequenced with 18.19 folds mapped read depth on average. Generally, we found that Beijing Black pig was genetically closer to commercial pig breeds by population genetic structure and genetic diversity analysis, and was also affected by Chinese domestic breeds Huainan pig and Min pig. These results are consistent with the cross-breeding history of Beijing Black pig. Selection signal detections were performed on three pig breeds, Beijing Black pig, Duroc and Large White, using three complementary methods (F_ST, θπ, and XP-EHH). In total, 1,167 significant selected regions and 392 candidate genes were identified. Functional annotations were enriched to pathways related to immune processes and meat and lipid metabolism. Finally, potential candidate genes, influencing meat quality (GPHA2, EHD1, HNF1A, C12orf43, GLTP, TRPV4, MVK, and MMAB), reproduction (PPP2R5B and MAP9), and disease resistance (OASL, ANKRD13A, and GIT2), were further detected by gene annotation analysis. Our results advanced the understanding of the genetic mechanism behind artificial selection of Beijing Black pigs, and provided theoretical basis for the subsequent breeding and genetic research of this breed.