Dynamic transcriptome profiles of skeletal muscle tissue across 11 developmental stages for both Tongcheng and Yorkshire pigs

Temporal Transcriptome Dynamics of Longissimus dorsi Reveals the Mechanism of the Differences in Muscle Development and IMF Deposition between Fuqing Goats and Nubian Goats

In this study, we measured the growth performance and intramuscular fat (IMF) content of the Longissimus dorsi (LD) of Fuqing goats (FQs) and Nubian goats (NBYs), which exhibit extreme phenotypic differences in terms of their production and meat quality traits. RNA-Seq analysis was performed, and transcriptome data were obtained from the LD tissue of 3-month fetuses (E3), 0-month lambs (0M), 3-month lambs (3M), and 12-month lambs (12M) to reveal the differences in the molecular mechanisms regulating the muscle development and IMF deposition between FQs and NBYs. The results showed that a higher body weight and average daily gain were observed in the NBYs at three developmental stages after birth, whereas a higher IMF content was registered in the FQs at 12M. Additionally, transcriptome profiles during the embryonic period and after birth were completely different for both FQs and NBYs. Moreover, DEGs (KIF23, CCDC69, CCNA2, MKI67, KIF11, RACGAP1, NUSAP1, SKP2, ZBTB18, NES, LOC102180034, CAPN6, TUBA1A, LOC102178700, and PEG10) significantly enriched in the cell cycle (ko04110) at E3 (FQs vs. NBYs), and DEGs (MRPS7, RPS8, RPL6, RPL4, RPS11, RPS10, RPL5, RPS6, RPL8, RPS13, RPS24, RPS15, RPL23) significantly enriched in ribosomes (ko03010) at 0M (FQs vs. NBYs) related to myogenic differentiation and fusion were identified. Meanwhile, the differences in glucose and lipid metabolism began at the E3 timepoint and continued to strengthen as growth proceeded in FQs vs. NBYs. DEGs (CD36, ADIROQR2, ACACA, ACACB, CPT1A, IGF1R, IRS2, LDH-A, PKM, HK2, PFKP, PCK1, GPI, FASN, FADS1, ELOVL6, HADHB, ACOK1, ACAA2, and ACSL4) at 3M (FQs vs. NBYs) and 12M (FQs vs. NBYs) significantly enriched in the AMPK signaling pathway (ko04152), insulin resistance (ko04931), the insulin signaling pathway (ko04910), fatty acid metabolism (ko01212), and glycolysis/gluconeogenesis (ko00010) related to IMF deposition were identified. Further, the results from this study provide the basis for future studies on the mechanisms regulating muscle development and IMF deposition in different breeds of goats, and the candidate genes identified could be used in the selection process.

Dynamic transcriptome profiles of skeletal muscle growth and development in Shaziling and Yorkshire pigs using RNA-sequencing

BACKGROUND

We previously demonstrated that Shaziling and Yorkshire pigs differ in growth rate and meat quality. However, the molecular mechanisms responsible for such phenotypic differences remain unclear. In the present study, we performed a transcriptomic analysis of 36 longissimus dorsi (LM) and 36 soleus (SM) muscle samples from Shaziling and Yorkshire pigs at six postnatal stages (30, 60, 90, 150, 210 and 300 days) to explore the differences in postnatal skeletal muscle of Shaziling and Yorkshire pigs.

RESULTS

Muscle morphological changes and the number of differentially expressed genes indicated the two stages of 60-90 days and 150-210 days were critical for the muscle growth and development in Shaziling pigs. Genes such as FLNC, COL1A1, NRAP, SMYD1, TNNI3, CRYAB and PDLIM3 played vital roles in the muscle growth, and genes such as CCDC71L, LPIN1, CPT1A, UCP3, NR4A3 and PDK4 played dominant roles in the lipid metabolism. Additionally, in contrast to the LM, the percentage of slow-twitch muscle fibers in the SM of both breeds consistently decreased from 30 to 150 days of age, but there was a significant rebound at 210 days of age. However, the percentage of slow-twitch muscle fibers in the SM of Shaziling pigs was higher than that in Yorkshire pigs, which may be associated with the calcium signaling pathway and the PPARβ/δ signaling pathway.

CONCLUSION

The present study detected two critical periods and many functional genes for the muscle growth and development of Shaziling pigs, and showed differences in muscle fiber characteristics between Shaziling and Yorkshire pigs. © 2024 Society of Chemical Industry.

Transcriptome profiling of longissimus dorsi during different prenatal stages to identify genes involved in intramuscular fat deposition in lean and obese pig breeds

BACKGROUND

There was significant difference in muscle development between fat-type and lean-type pig breeds.

METHODS AND RESULTS

In current study, transcriptome analysis and bioinformatics analysis were used to compare the difference in longissimus dorsi (LD) muscle at three time-points (38 days post coitus (dpc), 58 dpc, and 78 dpc ) between Huainan (HN) and Large white (LW) pig breeds. A total of 24500 transcripts were obtained in 18 samples, and 2319, 2799, and 3713 differently expressed genes (DEGs) were identified between these two breeds at 38 dpc, 58 dpc, and 78 dpc, respectively. And the number and foldchange of DEGs were increased, the alternative splice also increased. The cluster analysis of DEGs indicated the embryonic development progress of LD muscle between these two breeds was different. There were 539 shared DEGs between HN and LW at three stages, and the top-shared DEGs were associated with muscle development and lipid deposition, such as KLF4, NR4A1, HSP70, ZBTB16 and so on.

CONCLUSIONS

The results showed DEGs between Huainan (HN) and Large white (LW) pig breeds, and contributed to the understanding the muscle development difference between HN and LW, and provided basic materials for improvement of meat quality.

Transcriptomic Profiling of Meat Quality Traits of Skeletal Muscles of the Chinese Indigenous Huai Pig and Duroc Pig

The Huai pig is a well-known indigenous pig breed in China. The main advantages of Huai pigs over Western commercial pig breeds include a high intramuscular fat (IMF) content and good meat quality. There are significant differences in the meat quality traits of the same muscle part or different muscle parts of the same variety. To investigate the potential genetic mechanism underlying the meat quality differences in different pig breeds or muscle groups, longissimus dorsi (LD), psoas major (PM), and biceps femoris (BF) muscle tissues were collected from two pig breeds (Huai and Duroc). There were significant differences in meat quality traits and amino acid content. We assessed the muscle transcriptomic profiles using high-throughput RNA sequencing. The IMF content in the LD, PM, and BF muscles of Huai pigs was significantly higher than that in Duroc pigs (p < 0.05). Similarly, the content of flavor amino acids in the three muscle groups was significantly higher in Huai pigs than that in Duroc pigs (p < 0.05). We identified 175, 110, and 86 differentially expressed genes (DEGs) between the LD, PM, and BF muscles of the Huai and Duroc pigs, respectively. The DEGs of the different pig breeds and muscle regions were significantly enriched in the biological processes and signaling pathways related to muscle fiber type, IMF deposition, lipid metabolism, PPAR signaling, cAMP signaling, amino acid metabolism, and ECM-receptor interaction. Our findings might help improve pork yield by using the obtained DEGs for marker-assisted selection and providing a theoretical reference for evaluating and improving pork quality.

Integrated analysis of circRNA, lncRNA, miRNA and mRNA to reveal the ceRNA regulatory network of postnatal skeletal muscle development in Ningxiang pig

Introduction: The development of skeletal muscle is regulated by regulatory factors of genes and non-coding RNAs (ncRNAs). Methods: The objective of this study was to understand the transformation of muscle fiber type in the longissimus dorsi muscle of male Ningxiang pigs at four different growth stages (30, 90, 150, and 210 days after birth, n = 3) by histological analysis and whole transcriptome sequencing. Additionally, the study investigated the expression patterns of various RNAs involved in muscle fiber transformation and constructed a regulatory network for competing endogenous RNA (ceRNA) that includes circular RNA (circRNA)/long non-coding RNA (lncRNA)-microRNA (miRNA)-messenger RNA (mRNA). Results: Histomorphology analysis showed that the diameter of muscle fiber reached its maximum at 150 days after birth. The slow muscle fiber transformation showed a pattern of initial decrease followed by an increase. 29,963 circRNAs, 2,683 lncRNAs, 986 miRNAs and 22,411 mRNAs with expression level ≥0 were identified by whole transcriptome sequencing. Furthermore, 642 differentially expressed circRNAs (DEc), 505 differentially expressed lncRNAs (DEl), 316 differentially expressed miRNAs (DEmi) and 6,090 differentially expressed mRNAs (DEm) were identified by differential expression analysis. Functions of differentially expressed mRNA were identified by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). GO enrichment analysis indicates that 40 known genes and 6 new genes are associated with skeletal muscle development. Additionally, KEGG analysis shows that these genes regulate skeletal muscle development via MAPK, FoxO, Hedgehog, PI3K-Akt, Notch, VEGF and other signaling pathways. Through protein-protein interaction (PPI) and transcription factor prediction (TFP), the action mode of skeletal muscle-related genes was explored. PPI analysis showed that there were stable interactions among 19 proteins, meanwhile, TFP analysis predicted 22 transcription factors such as HMG20B, MYF6, MYOD1 and MYOG, and 12 of the 19 interacting proteins were transcription factors. The regulatory network of ceRNA related to skeletal muscle development was constructed based on the correlation of various RNA expression levels and the targeted binding characteristics with miRNA. The regulatory network included 31 DEms, 59 miRNAs, 667 circRNAs and 224 lncRNAs. conclusion: Overall, the study revealed the role of ceRNA regulatory network in the transformation of skeletal muscle fiber types in Ningxiang pigs, which contributes to the understanding of ceRNA regulatory network in Ningxiang pigs during the skeletal muscle development period.

Transcriptome Profiling of Different Developmental Stages on Longissimus Dorsi to Identify Genes Underlying Intramuscular Fat Content in Wannanhua Pigs

Intramuscular fat (IMF) is a key index to measure the tenderness and flavor of pork. Wannanhua pig, a famous indigenous pig breed in Anhui Province, is renowned for its high lipid deposition and high genetic divergence, making it an ideal model for investigating the lipid position trait mechanisms in pigs. However, the regulatory mechanisms of lipid deposition and development in pigs remain unclear. Furthermore, the temporal differences in gene regulation are based on muscle growth and IMF deposition. The purpose of this study was to study the expression changes of longissimus dorsi (LD) at different growth stages of WH pigs at the molecular level, to screen the candidate genes and signaling pathways related to IMF during development by transcriptome sequencing technology, and to explore the transcriptional regulation mechanism of IMF deposition-related genes at different development stages. In total, 616, 485, and 1487 genes were differentially expressed between LD60 and LD120, LD120 and LD240, and LD60 and LD240, respectively. Numerous differentially expressed genes (DEGs) associated with lipid metabolism and muscle development were identified, and most of them were involved in IMF deposition and were significantly up-regulated in LD120 and LD240 compared to LD60. STEM (Short Time-series Expression Miner) analysis indicated significant variations in the mRNA expression across distinct muscle development stages. The differential expression of 12 selected DEGs was confirmed by RT-qPCR. The results of this study contribute to our understanding of the molecular mechanism of IMF deposition and provide a new way to accelerate the genetic improvement of pork quality.

Identification and characterization of circRNAs related to meat quality during embryonic development of the longissimus dorsi muscle in two pig breeds

Meat quality, an important economic trait, is regulated by many factors, especially by genetic factors, including coding genes, miRNAs, and lncRNAs. Recent studies have elucidated that circRNAs also play a key role in muscle development and lipid deposition. However, the functions and regulatory mechanisms of circRNAs in meat quality remain mostly unknown. The circRNA expression profiles between Huainan pigs (Chinese indigenous pigs, fat-type, Huainan HN) and Large White pigs (Western commercial pigs, lean-type, LW) in the longissimus dorsi (LD) muscle at 38, 58, and 78 days post conception (dpc) were compared by sequencing. In total, 39,887 circRNAs were identified in 18 samples, and 60, 78, and 86 differentially expressed circRNAs (DECs) were found at the three stages mentioned above between these two breeds. The parent genes of DECs were enriched in myogenesis, proliferation, adipogenesis and muscle fiber-type transition. The circRNA-miRNA interaction networks included 38 DECs and 47 miRNAs, and these miRNAs were involved in muscle development and lipid metabolism. Two shared DECs (circ_0030593 and circ_0032760) of these three stages were selected, their head-to-tail junction sites were validated by Sanger sequencing, and RT‒qPCR results suggested that these two DECs might be involved in intramuscular fat deposition. These findings provide a basis for understanding the role of circRNAs in meat quality.

Long-range interaction within the chromatin domain determines regulatory patterns in porcine skeletal muscle

Spatial chromatin structure is crucial for understanding the early growth and development of porcine skeletal muscle. However, its characteristic of 3D architecture and elaborate regulation of gene transcription remains unclear. In this study, ChIA-PET method is used to study the changes of early chromatin three-dimensional structure in skeletal muscle of lean type Yorkshire pig and fat type Meishan pig. Integrating the in situ Hi-C data revealed the 3D architecture and long-range interaction of the porcine muscle. The results showed the CTCF/RNAPII mediated long-range interaction shapes the different chromatin architecture and dominates the unique regulation of enhancers. In addition, the results revealed that key myogenic genes like ssc-mir-1 had a unique enhancer regulation function in myogenesis. Interestingly, the FGF6 gene is of breed-specific regulation, implying the difference between two breeds in skeletal muscle development. Our research thus may provide a clue for the porcine genetic improvement of skeletal muscle.

Genome-Wide Association Study of Growth Traits in a Four-Way Crossbred Pig Population

Growth traits are crucial economic traits in the commercial pig industry and have a substantial impact on pig production. However, the genetic mechanism of growth traits is not very clear. In this study, we performed a genome-wide association study (GWAS) based on the specific-locus amplified fragment sequencing (SLAF-seq) to analyze ten growth traits on 223 four-way intercross pigs. A total of 227,921 highly consistent single nucleotide polymorphisms (SNPs) uniformly dispersed throughout the entire genome were used to conduct GWAS. A total of 53 SNPs were identified for ten growth traits using the mixed linear model (MLM), of which 18 SNPs were located in previously reported quantitative trait loci (QTL) regions. Two novel QTLs on SSC4 and SSC7 were related to average daily gain from 30 to 60 kg (ADG30–60) and body length (BL), respectively. Furthermore, 13 candidate genes (ATP5O, GHRHR, TRIM55, EIF2AK1, PLEKHA1, BRAP, COL11A2, HMGA1, NHLRC1, SGSM1, NFATC2, MAML1, and PSD3) were found to be associated with growth traits in pigs. The GWAS findings will enhance our comprehension of the genetic architecture of growth traits. We suggested that these detected SNPs and corresponding candidate genes might provide a biological foundation for improving the growth and production performance of pigs in swine breeding.

Notch signaling leads to a slower progression of embryonic myogenic differentiation in Landrace than in Langtang pigs

Delving into porcine embryonic myogenesis is the key to elucidate the complex regulation of breed-specific differences in growth performance and meat production. Increasing evidence proves that pigs with less meat production show earlier embryonic myogenesis, but little is known about the underlying mechanisms. In this study, we examine the longissimus dorsi muscle (LDM) by immunohistochemistry and confirm that the differentiation of myogenic progenitors is increased ( P<0.05) in Lantang (LT, fatty) pigs compared with that in Landrace (LR, lean) pigs, which results in more ( P<0.001) differentiated myoblasts (Pax7 ^-/MyoD ⁺) and less ( P<0.001) myogenic progenitors (Pax7 ⁺/MyoD ^-) in LT pigs at 35 days post-conception (35dpc). Additionally, embryonic myogenic progenitors isolated from LT pigs show greater ( P<0.001) differentiation capacity with earlier expression of MyoD compared with those from LR pigs. Moreover, Notch signaling is more active ( P<0.05) in LR pig myogenic progenitors than in LT pig myogenic progenitors. Inhibition of Notch signaling in LR myogenic progenitors suppresses Pax7 expression and increases MyoD expression, thus promoting myogenic differentiation. Consistently, the process of myogenic progenitors differentiating into myoblasts in ex vivo embryo limbs is accelerated when Notch signaling is inhibited. These results indicate that Notch signaling facilitates the maintenance of myogenic progenitors and antagonizes myogenic differentiation by promoting Pax7 expression and preventing MyoD expression in LR pigs.

Genetic regulation and variation of expression of miRNA and mRNA transcripts in fetal muscle tissue in the context of sex, dam and variable fetal weight

BACKGROUND

Impaired skeletal muscle growth in utero can result in reduced birth weight and pathogenesis of intrauterine growth restriction. Fetal and placental growth is influenced by many factors including genetic, epigenetic and environmental factors. In fact, the sex and genotype of the fetus itself, as well as the mother providing it with a suitable environment, influence the growth of the fetus. Hence, our goal was to decipher and elucidate the molecular pathways of developmental processes mediated by miRNAs and mRNAs in fetal muscle tissue in the context of sex, dam, and fetal weight. Therefore, we analyse the variation of miRNA and mRNA expression in relation to these factors. In addition, the coincidence of genetic regulation of these mRNAs and miRNAs, as revealed by expression quantitative trait loci (eQTL) analyses, with sex-, mother- and weight-associated expression was investigated.

METHODS

A three-generation pig F2 population (n = 118) based on reciprocal crossing of German Landrace (DL) and Pietrain (Pi) was used. Genotype information and transcriptomic data (mRNA and miRNA) from longissimus dorsi muscle (LDM) of pig fetuses sampled at 63 days post-conception (dpc) were used for eQTL analyses.

RESULTS

The transcript abundances of 13, 853, and 275 probe-sets were influenced by sex, dam and fetal weight at 63 dpc, respectively (FDR < 5%). Most of significant transcripts affected by sex were located on the sex chromosomes including KDM6A and ANOS1 or autosomes including ANKS1B, LOC100155138 and miR-153. The fetal muscle transcripts associated with fetal weight indicated clearer metabolic directions than maternally influenced fetal muscle transcripts. Moreover, coincidence of genetic regulation (eQTL) and variation in transcript abundance due to sex, dam and fetal weight were identified.

CONCLUSIONS

Integrating information on eQTL, sex-, dam- and weight-associated differential expression and QTL for fetal weight allowed us to identify molecular pathways and shed light on the basic biological processes associated with differential muscle development in males and females, with implications for adaptive fetal programming.

Profiling and Functional Analysis of mRNAs during Skeletal Muscle Differentiation in Goats

Skeletal myogenesis is a complicated biological event that involves a succession of tightly controlled gene expressions. In order to identify novel regulators of this process, we performed mRNA-Seq studies of goat skeletal muscle satellite cells (MuSCs) cultured under proliferation (GM) and differentiation (DM1/DM5) conditions. A total of 19,871 goat genes were expressed during these stages, 198 of which represented novel transcripts. Notably, in pairwise comparisons at the different stages, 2551 differentially expressed genes (DEGs) were identified (p < 0.05), including 1560 in GM vs. DM1, 1597 in GM vs. DM5, and 959 in DM1 vs. DM5 DEGs. The time-series expression profile analysis clustered the DEGs into eight gene groups, three of which had significantly upregulated and downregulated patterns (p < 0.05). Functional enrichment analysis showed that DEGs were enriched for essential biological processes such as muscle structure development, muscle contraction, muscle cell development, striated muscle cell differentiation, and myofibril assembly, and were involved in pathways such as the MAPK, Wnt and PPAR signaling pathways. Moreover, the expression of eight DEGs (MYL2, DES, MYOG, FAP, PLK2, ADAM, WWC1, and PRDX1) was validated. These findings offer novel insights into the transcriptional regulation of skeletal myogenesis in goats.

Profiling of open chromatin in developing pig (Sus scrofa) muscle to identify regulatory regions

There is very little information about how the genome is regulated in domestic pigs (Sus scrofa). This lack of knowledge hinders efforts to define and predict the effects of genetic variants in pig breeding programs. To address this knowledge gap, we need to identify regulatory sequences in the pig genome starting with regions of open chromatin. We used the "Improved Protocol for the Assay for Transposase-Accessible Chromatin (Omni-ATAC-Seq)" to identify putative regulatory regions in flash-frozen semitendinosus muscle from 24 male piglets. We collected samples from the smallest-, average-, and largest-sized male piglets from each litter through five developmental time points. Of the 4661 ATAC-Seq peaks identified that represent regions of open chromatin, >50% were within 1 kb of known transcription start sites. Differential read count analysis revealed 377 ATAC-Seq defined genomic regions where chromatin accessibility differed significantly across developmental time points. We found regions of open chromatin associated with downregulation of genes involved in muscle development that were present in small-sized fetal piglets but absent in large-sized fetal piglets at day 90 of gestation. The dataset that we have generated provides a resource for studies of genome regulation in pigs and contributes valuable functional annotation information to filter genetic variants for use in genomic selection in pig breeding programs.

Changes in Biceps femoris Transcriptome along Growth in Iberian Pigs Fed Different Energy Sources and Comparative Analysis with Duroc Breed

Simple Summary

The genetic mechanisms that regulate biological processes, such as skeletal muscle development and growth, or intramuscular fat deposition, have attracted great interest, given their impact on production traits and meat quality. In this sense, a comparison of the transcriptome of skeletal muscle between phenotypically different pig breeds, or along growth, could be useful to improve the understanding of the molecular processes underlying the differences in muscle metabolism and phenotypic traits, potentially driving the identification of causal genes, regulators and metabolic pathways involved in their variability.

Abstract

This experiment was conducted to investigate the effects of developmental stage, breed, and diet energy source on the genome-wide expression, meat quality traits, and tissue composition of biceps femoris muscle in growing pure Iberian and Duroc pigs. The study comprised 59 Iberian (IB) and 19 Duroc (DU) animals, who started the treatment at an average live weight (LW) of 19.9 kg. The animals were kept under identical management conditions and fed two diets with different energy sources (6% high oleic sunflower oil or carbohydrates). Twenty-nine IB animals were slaughtered after seven days of treatment at an average LW of 24.1 kg, and 30 IB animals plus all the DU animals were slaughtered after 47 days at an average LW of 50.7 kg. The main factors affecting the muscle transcriptome were age, with 1832 differentially expressed genes (DEGs), and breed (1055 DEGs), while the effect of diet on the transcriptome was very small. The results indicated transcriptome changes along time in Iberian animals, being especially related to growth and tissue development, extracellular matrix (ECM) composition, and cytoskeleton organization, with DEGs affecting relevant functions and biological pathways, such as myogenesis. The breed also affected functions related to muscle development and cytoskeleton organization, as well as functions related to solute transport and lipid and carbohydrate metabolism. Taking into account the results of the two main comparisons (age and breed effects), we can postulate that the Iberian breed is more precocious than the Duroc breed, regarding myogenesis and muscle development, in the studied growing stage.

Muscle transcriptome provides the first insight into the dynamics of gene expression with progression of age in sheep

The dynamic synergy of genes and pathways in muscles in relation to age affects the muscle characteristics. Investigating the temporal changes in gene expression will help illustrate the molecular mechanisms underlying muscle development. Here we report the gene expression changes in skeletal muscles through successive age groups in Bandur, a meat type sheep of India. RNA sequencing data was generated from the longissimus thoracis muscles from four age groups, ranging from lamb to adult. Analysis of 20 highest expressed genes common across the groups revealed muscle protein, phosphorylation, acetylation, metal binding and transport as significant functions. Maximum differentiation was observed after 2.5–3 years on transition from lambs to adult. Transcriptional regulation by the TFAP2 transcription factors, IL-6 signaling and PI3K/AKT signaling pathways were enriched in younger animals. The gene-protein network demarcated key interactive genes involved in muscle development and proliferation that can be used as candidates for future research on improvement of muscle characteristics.

The genome variation and developmental transcriptome maps reveal genetic differentiation of skeletal muscle in pigs

Natural and artificial directional selections have resulted in significantly genetic and phenotypic differences across breeds in domestic animals. However, the molecular regulation of skeletal muscle diversity remains largely unknown. Here, we conducted transcriptome profiling of skeletal muscle across 27 time points, and performed whole-genome re-sequencing in Landrace (lean-type) and Tongcheng (obese-type) pigs. The transcription activity decreased with development, and the high-resolution transcriptome precisely captured the characterizations of skeletal muscle with distinct biological events in four developmental phases: Embryonic, Fetal, Neonatal, and Adult. A divergence in the developmental timing and asynchronous development between the two breeds was observed; Landrace showed a developmental lag and stronger abilities of myoblast proliferation and cell migration, whereas Tongcheng had higher ATP synthase activity in postnatal periods. The miR-24-3p driven network targeting insulin signaling pathway regulated glucose metabolism. Notably, integrated analysis suggested SATB2 and XLOC_036765 contributed to skeletal muscle diversity via regulating the myoblast migration and proliferation, respectively. Overall, our results provide insights into the molecular regulation of skeletal muscle development and diversity in mammals.

Dynamic transcriptome profiles of postnatal porcine skeletal muscle growth and development

BACKGROUND

Skeletal muscle growth and development are closely associated with the quantity and quality of pork production. We performed a transcriptomic analysis of 12 Longissimus dorsi muscle samples from Tibetan piglets at four postnatal stages of 0, 14, 30, and 60 days using RNA sequencing.

RESULTS

According to the pairwise comparisons between the libraries of the muscle samples at the four postnatal stages, a total of 4115 differentially expressed genes (DEGs) were identified in terms of |log2(fold change)| ≥ 1 and an adjusted P value < 0.01. Short-time series expression miner (STEM) analysis of the DEGs identified eight significantly different expression profiles, which were divided into two clusters based on the expression pattern. DEGs in cluster I displayed a pattern of decreasing to a nadir, and then a rise, and the significantly enriched gene ontology (GO) terms detected using them were involved in multiple processes, of which the cell cycle, immunocyte activation and proliferation, as well as actin cytoskeleton organization, were the top three overrepresented processes based on the GO terms functional classification. DEGs in cluster II displayed a pattern of increasing to a peak, then declining, which mainly contributed to protein metabolism. Furthermore, besides the pathways related to immune system, a few diseases, and protein metabolism, the DEGs in clusters I and II were significantly enriched in pathways related to muscle growth and development, such as the Rap1, PI3K-Akt, AMPK, and mTOR signaling pathways.

CONCLUSIONS

This study revealed GO terms and pathways that could affect the postnatal muscle growth and development in piglets. In addition, this study provides crucial information concerning the molecular mechanisms of muscle growth and development as well as an overview of the piglet transcriptome dynamics throughout the postnatal period in terms of growth and development.

In Utero Fetal Weight in Pigs Is Regulated by microRNAs and Their Target Genes

Impaired skeletal muscle growth in utero can result in reduced birth weight and poor carcass quality in pigs. Recently, we showed the role of microRNAs (miRNAs) and their target genes in prenatal skeletal muscle development and pathogenesis of intrauterine growth restriction (IUGR). In this study, we performed an integrative miRNA-mRNA transcriptomic analysis in longissimus dorsi muscle (LDM) of pig fetuses at 63 days post conception (dpc) to identify miRNAs and genes correlated to fetal weight. We found 13 miRNAs in LDM significantly correlated to fetal weight, including miR-140, miR-186, miR-101, miR-15, miR-24, miR-29, miR-449, miR-27, miR-142, miR-99, miR-181, miR-199, and miR-210. The expression of these miRNAs decreased with an increase in fetal weight. We also identified 1315 genes significantly correlated to fetal weight at 63 dpc, of which 135 genes were negatively correlated as well as identified as potential targets of the above-listed 13 miRNAs. These miRNAs and their target genes enriched pathways and biological processes important for fetal growth, development, and metabolism. These results indicate that the transcriptomic profile of skeletal muscle can be used to predict fetal weight, and miRNAs correlated to fetal weight can serve as potential biomarkers of prenatal fetal health and growth.

Prenatal Skeletal Muscle Transcriptome Analysis Reveals Novel MicroRNA-mRNA Networks Associated with Intrauterine Growth Restriction in Pigs

Intrauterine growth restriction (IUGR) occurs in 15–20% of pig neonates and poses huge economic losses to the pig industry. IUGR piglets have reduced skeletal muscle growth, which may persist after birth. Prenatal muscle growth is regulated by complex molecular pathways that are not well understood. MicroRNAs (miRNAs) have emerged as the main regulators of vital pathways and biological processes in the body. This study was designed to identify miRNA–mRNA networks regulating prenatal skeletal muscle development in pigs. We performed an integrative miRNA–mRNA transcriptomic analysis in longissimus dorsi muscle from IUGR fetuses and appropriate for gestational age (AGA) fetuses at 63 days post conception. Our data showed that 47 miRNAs and 3257 mRNAs were significantly upregulated, and six miRNAs and 477 mRNAs were significantly downregulated in IUGR compared to AGA fetuses. Moreover, 47 upregulated miRNAs were negatively correlated and can potentially target 326 downregulated genes, whereas six downregulated miRNAs were negatively correlated and can potentially target 1291 upregulated genes. These miRNA–mRNA networks showed enrichment in biological processes and pathways critical for fetal growth, development, and metabolism. The miRNA–mRNA networks identified in this study can potentially serve as indicators of prenatal fetal growth and development as well as postnatal carcass quality.

Key Genes Regulating Skeletal Muscle Development and Growth in Farm Animals

Simple Summary

Skeletal muscle mass is an important economic trait, and muscle development and growth is a crucial factor to supply enough meat for human consumption. Thus, understanding (candidate) genes regulating skeletal muscle development is crucial for understanding molecular genetic regulation of muscle growth and can be benefit the meat industry toward the goal of increasing meat yields. During the past years, significant progress has been made for understanding these mechanisms, and thus, we decided to write a comprehensive review covering regulators and (candidate) genes crucial for muscle development and growth in farm animals. Detection of these genes and factors increases our understanding of muscle growth and development and is a great help for breeders to satisfy demands for meat production on a global scale.

Abstract

Farm-animal species play crucial roles in satisfying demands for meat on a global scale, and they are genetically being developed to enhance the efficiency of meat production. In particular, one of the important breeders’ aims is to increase skeletal muscle growth in farm animals. The enhancement of muscle development and growth is crucial to meet consumers’ demands regarding meat quality. Fetal skeletal muscle development involves myogenesis (with myoblast proliferation, differentiation, and fusion), fibrogenesis, and adipogenesis. Typically, myogenesis is regulated by a convoluted network of intrinsic and extrinsic factors monitored by myogenic regulatory factor genes in two or three phases, as well as genes that code for kinases. Marker-assisted selection relies on candidate genes related positively or negatively to muscle development and can be a strong supplement to classical selection strategies in farm animals. This comprehensive review covers important (candidate) genes that regulate muscle development and growth in farm animals (cattle, sheep, chicken, and pig). The identification of these genes is an important step toward the goal of increasing meat yields and improves meat quality.

Skeletal Muscle Transcriptome Analysis of Hanzhong Ma Duck at Different Growth Stages Using RNA-Seq

As one of the most important poultry worldwide, ducks (Anas platyrhynchos) are raised mainly for meat and egg products, and muscle development in ducks is important for meat production. Therefore, an investigation of gene expression in duck skeletal muscle would significantly contribute to our understanding of muscle development. In this study, twenty-four cDNA libraries were constructed from breast and leg muscles of Hanzhong Ma ducks at day 17, 21, 27 of the embryo and postnatal at 6-month-old. High-throughput sequencing and bioinformatics were used to determine the abundances and characteristics of transcripts. A total of 632,172,628 (average 52,681,052) and 637,213,938 (average 53,101,162) reads were obtained from the sequencing data of breast and leg muscles, respectively. Over 71.63% and 77.36% of the reads could be mapped to the Anas platyrhynchos genome. In the skeletal muscle of Hanzhong duck, intron variant (INTRON), synonymous variant (SYNONYMOUS_CODING), and prime 3' UTR variant (UTR_3_PRIME) were the main single nucleotide polymorphisms (SNP) annotation information, and "INTRON", "UTR_3_PRIME", and downstream-gene variant (DOWNSTREAM) were the main insertion-deletion (InDel) annotation information. The predicted number of alternative splicing (AS) in all samples were mainly alternative 5' first exon (transcription start site)-the first exon splicing (TSS) and alternative 3' last exon (transcription terminal site)-the last exon splicing (TTS). Besides, there were 292 to 2801 annotated differentially expressed genes (DEGs) in breast muscle and 304 to 1950 annotated DEGs in leg muscle from different databases. It is worth noting that 75 DEGs in breast muscle and 49 DEGs in leg muscle were co-expressed at all developmental points of comparison, respectively. The RNA-Seq data were confirmed to be reliable by qPCR. The identified DEGs, such as CREBL2, RHEB, GDF6, SHISA2, MYLK2, ACTN3, RYR3, and STMN1, were specially highlighted, indicating their strong associations with muscle development in the Hanzhong Ma duck. KEGG pathway analysis suggested that regulation of actin cytoskeleton, oxidative phosphorylation, and focal adhesion were involved in the development of skeletal muscle. The findings from this study can contribute to future investigations of the growth and development mechanism in duck skeletal muscle.

RNA-seq Analysis Reveals Mitochondrial and Cuticular Protein Genes Are Associated with Phosphine Resistance in the Rusty Grain Beetle (Coleoptera:Laemophloeidae)

The rusty grain beetle, Cryptolestes ferrugineus (Stephens), is a serious pest of stored grain, which has developed high levels of resistance to phosphine. In this study, five geographically distant populations of C. ferrugineus had been collected in China, specifically in granaries where phosphine fumigant is used for pest control, and they showed a high resistance ratio up to 1,907 (LC50 = 21.0 mg/liter). Then, a reference transcriptome was constructed to use as a basis for investigating the molecular mechanisms of phosphine resistance in this species, which consisted of 47,006 unigenes with a mean length of 1,090. Subsequently, the RNA-Seq analysis of individuals from the most susceptible and resistant populations led to the identification of 54 genes that are differentially expressed. GO and KEGG analysis demonstrated that genes associated with mitochondrial and respiration functions were significantly enriched. Also, the 'structural constituent of cuticle' term was annotated in the GO enrichment analysis and further qRT-PCR confirmed that the expression levels of nine cuticular protein genes were significantly increased in the resistant population. In conclusion, we present here a transcriptome-wide overview of gene expression changes between resistant and susceptible populations of C. ferrugineus, and this in turn documents that mitochondria and cuticular protein genes may play together a crucial role in phosphine resistance. Further gene function analysis should enable the provision of advice to expedite resistance management decisions.

LncRNA IMFlnc1 promotes porcine intramuscular adipocyte adipogenesis by sponging miR-199a-5p to up-regulate CAV-1

Background

Local Chinese local pig breeds have thinner muscle fiber and higher intramuscular-fat (IMF) content. But its regulation mechanism has not been discussed in-depth. Studies indicated that long non coding RNAs (lncRNAs) play important role in muscle and fat development.

Results

The lncRNAs expressional differences in the longissimus dorsi (LD) muscle were identified between Huainan pigs (local Chinese pigs, fat-type, HN) and Large White pigs (lean-type, LW) at 38, 58, and 78 days post conception (dpc). In total, 2131 novel lncRNAs were identified in 18 samples, and 291, 305, and 683 differentially expressed lncRNAs (DELs) were found between these two breeds at three stages, respectively. The mRNAs that co-expressed with these DELs were used for GO and KEGG analysis, and the results showed that muscle development and energy metabolism were more active at 58 dpc in HN, but at 78 dpc in LW pigs. Muscle cell differentiation and myofibril assembly might associated with earlier myogenesis and primary-muscle-fiber assembly in HN, and cell proliferation, insulin, and the MAPK pathway might be contribute to longer proliferation and elevated energy metabolism in LW pigs at 78 dpc. The PI3K/Akt and cAMP pathways were associated with higher IMF deposition in HN. Intramuscular fat deposition-associated long noncoding RNA 1 (IMFlnc1) was selected for functional verification, and results indicated that it regulated the expressional level of caveolin-1 (CAV-1) by acting as competing endogenous RNA (ceRNA) to sponge miR-199a-5p.

Conclusions

Our data contributed to understanding the role of lncRNAs in porcine-muscle development and IMF deposition, and provided valuable information for improving pig-meat quality.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-020-00324-8.

The potential role of PSMA6 in modulating fat deposition in pigs by promoting preadipocyte proliferation and differentiation

To investigate whether the proteasome subunit alpha 6 (PSMA6) gene has an effect on fat deposition, the gene expression profile was first detected in Berkshire pigs and Jinhua pigs (JHP). The results demonstrated that significantly higher levels of mRNA expression were identified in adipose tissues and the liver. Interestingly, when compared to the longissimus dorsi muscle (LDM) in each breed, it was discovered that the expression levels of the PSMA6 gene in these tissues of JHP were considerably higher than those in Berkshire pigs. Meantime, some significant correlations of PSMA6 mRNA expression in lipid metabolism-related tissues such as the liver and LDM with the marbling score, as well as the content of intramuscular fat (IMF), in pigs were found by correlation coefficient analysis. To further explore the effects of PSMA6 expression on fat deposition, we performed PSMA6 overexpression in 3T3-L1 cells via Lentivirus infection. Our results indicated that PSMA6 could promote cell proliferation and accelerate cell division. It was also found that the transcription factors CCAAT/enhancer-binding protein alpha (CEBPA) and peroxisome proliferator-activated receptor gamma (PPARG), as well as the key genes related to adipogenesis, were upregulated, while the genes related to fat oxidation were significantly downregulated, which indicated that the PSMA6 gene could stimulate the differentiation of preadipocytes.

Functional contexts of adipose and gluteal muscle tissue gene co-expression networks in the domestic horse

A gene's response to an environment is tightly bound to the underlying genetic variation present in an individual's genome and varies greatly depending on the tissue it is being expressed in. Gene co-expression networks provide a mechanism to understand and interpret the collective transcriptional responses of genes. Here, we use the Camoco co-expression network framework to characterize the transcriptional landscape of adipose and gluteal muscle tissue in 83 domestic horses (Equus caballus) representing 5 different breeds. In each tissue, gene expression profiles, capturing transcriptional response due to variation across individuals, were used to build two separate, tissue-focused, genotypically-diverse gene co-expression networks. The aim of our study was to identify significantly co-expressed clusters of genes in each tissue, then compare the clusters across networks to quantify the extent that clusters were found in both networks as well as to identify clusters found in a single network. The known and unknown functions for each network were quantified using complementary, supervised and unsupervised approaches. First, supervised ontological enrichment was utilized to quantify biological functions represented by each network. Curated ontologies (GO and KEGG) were used to measure the known biological functions present in each tissue. Overall, a large percentage of terms (40.3% of GO and 41% of KEGG) were co-expressed in at least one tissue. Many terms were co-expressed in both tissues, however a small proportion of terms exhibited single tissue co-expression suggesting functional differentiation based on curated, functional annotation. To complement this, an unsupervised approach not relying on ontologies was employed. Strongly co-expressed sets of genes defined by Markov clustering identified sets of unannotated genes showing similar patterns of co-expression within a tissue. We compared gene sets across tissues and identified clusters of genes the either segregate in co-expression by tissue or exhibit high levels of co-expression in both tissues. Clusters were also integrated with GO and KEGG ontologies to identify gene sets containing previously curated annotations versus unannotated gene sets indicating potentially novel biological function. Coupling together these transcriptional datasets, we mapped the transcriptional landscape of muscle and adipose setting up a generalizable framework for interpreting gene function for additional tissues in the horse and other species.

Impact of genotype, body weight and sex on the prenatal muscle transcriptome of Iberian pigs

Growth is dependent on genotype and diet, even at early developmental stages. In this study, we investigated the effects of genotype, sex, and body weight on the fetal muscle transcriptome of purebred Iberian and crossbred Iberian x Large White pigs sharing the same uterine environment. RNA sequencing was performed on 16 purebred and crossbred fetuses with high body weight (340±14g and 415±14g, respectively) and 16 with low body weight (246±14g and 311±14g, respectively), on gestational day 77. Genotype had the greatest effect on gene expression, with 645 genes identified as differentially expressed (DE) between purebred and crossbred animals. Functional analysis showed differential regulation of pathways involved in energy and lipid metabolism, muscle development, and tissue disorders. In purebred animals, fetal body weight was associated with 35 DE genes involved in development, lipid metabolism and adipogenesis. In crossbred animals, fetal body weight was associated with 60 DE genes involved in muscle development, viability, and immunity. Interestingly, the results suggested an interaction genotype*weight for some DE genes. Fetal sex had only a modest effect on gene expression. This study allowed the identification of genes, metabolic pathways, biological functions and regulators related to fetal genotype, weight and sex, in animals sharing the same uterine environment. Our findings contribute to a better understanding of the molecular events that influence prenatal muscle development and highlight the complex interactions affecting transcriptional regulation during development.

Genome-Wide Analysis of Circular RNAs Mediated ceRNA Regulation in Porcine Embryonic Muscle Development

Many circular RNAs (circRNAs) have been discovered in various tissues and cell types in pig. However, the temporal expression pattern of circRNAs during porcine embryonic muscle development remains unclear. Here, we present a panorama view of circRNA expression in embryonic muscle development at 33-, 65-, and 90-days post-coitus (dpc) from Duroc pigs. An unbiased analysis reveals that more than 5,000 circRNAs specifically express in embryonic muscle development. The amount and complexity of circRNA expression is most pronounced in skeletal muscle at day 33 of gestation. Our circRNAs annotation analyses show that “hot-spot” genes produce multiple circRNA isoforms and RNA binding protein (RBPs) may regulate the biogenesis of circRNAs. Furthermore, we observed that host genes of differentially expressed circRNA across porcine muscle development are enriched in skeletal muscle function. A competing endogenous RNA (ceRNA) network analysis of circRNAs reveals that circRNAs regulate muscle gene expression by functioning as miRNA sponges. Finally, our experimental validation demonstrated that circTUT7 regulate the expression of HMG20B in a ceRNA mechanism. Our analyses show that circRNAs are dynamically expressed and interacting with muscle genes through ceRNA manner, suggesting their critical functions in embryonic skeletal muscle development.

Integrated miRNA-mRNA transcriptomic analysis reveals epigenetic-mediated embryonic muscle growth differences between Wuzhishan and Landrace pigs1

Pig is one of the major dietary protein sources for human consumption, from which muscle is the largest protein origin. However, molecular mechanisms concerning early porcine embryonic muscle development distinctions between pig breeds are still unclear. In this study, an integrated analysis of transcriptome and miRNAome was conducted using longissimus dorsi muscle of 4 early embryonic stages around the primary myofiber formation time (18-, 21-, 28-, and 35-d post coitus) from 2 pig breeds (Landrace [LR] and Wuzhishan [WZS]) differing in meat mass. The global miRNA/mRNA expression profile showed that WZS prepared for myogenic developmental processes earlier than LR. After identifying and analyzing the interaction network of top 100 up-/down-regulated miRNA and their target genes, we were able to find 3 gene clusters: chromatin modification-related (Chd2, H3f3a, Chd6, and Mll1), myogenesis-related (Pax3, Pbx1, Mef2a, and Znf423), and myosin component-related (Mylk, Myo5a, Mylk4, Myh9, and Mylk2) gene clusters. These genes may involve in miRNA-gene myogenic regulatory network that plays vital role in regulating distinct early porcine embryonic myogenic processes between LR and WZS. In summary, our study reveals an epigenetic-mediated myogenic regulatory axial that will help us to decipher molecular mechanisms concerning early porcine embryonic muscle development distinctions between pig breeds.

Earlier demethylation of myogenic genes contributes to embryonic precocious terminal differentiation of myoblasts in miniature pigs

Here, we performed whole-genome bisulfite sequencing of longissimus dorsi muscle from Landrace and Wuzhishan (WZS) miniature pigs during 18, 21, and 28 d postcoitum. It was uncovered that in regulatory regions only around transcription start sites (TSSs), gene expression and methylation showed negative correlation, whereas in gene bodies, positive correlation occurred. Furthermore, earlier myogenic gene demethylation around TSSs and earlier hypermethylation of myogenic genes in gene bodies were considered to trigger their earlier expression in miniature pigs. Furthermore, by analyzing the methylation pattern of the myogenic differentiation 1(MyoD) promoter and distal enhancer, we found that earlier demethylation of the MyoD distal enhancer in WZSs contributes to its earlier expression. Moreover, DNA demethylase Tet1 was found to be involved in the demethylation of the myogenin promoter and promoted immortalized mouse myoblast cell line (C2C12) and porcine embryonic myogenic cell differentiation. This study reveals that earlier demethylation of myogenic genes contributes to precocious terminal differentiation of myoblasts in miniature pigs.-Zhang, X., Nie, Y., Cai, S., Ding, S., Fu, B., Wei, H., Chen, L., Liu, X., Liu, M., Yuan, R., Qiu, B., He, Z., Cong, P., Chen, Y., Mo, D. Earlier demethylation of myogenic genes contributes to embryonic precocious terminal differentiation of myoblasts in miniature pigs.

Comparative analyses of longissimus muscle miRNAomes reveal microRNAs associated with differential regulation of muscle fiber development between Tongcheng and Yorkshire pigs

Tongcheng (TC) and Yorkshire (YK) are two pig breeds with distinctive muscle morphology. Porcine microRNAome (miRNAome) of the longissimus muscle during five developmental stages (40, 55, 63, 70, and 90 days post coitum (dpc)) was explored by Solexa sequencing in the present study to find miRNAs involved in the different regulation of skeletal muscle development between the two breeds. A total of 320 known porcine miRNAs, 64 miRNAs corresponding to other mammals, and 224 potentially novel miRNAs were identified. Principal component analysis (PCA) and hierarchical cluster analysis (HCA) suggested that the factor “pig breed” affected the miRNA expression profiles to a lesser extent than the factor “developmental stage”. Fifty-seven miRNAs were differentially expressed (DE) between the neighbor developmental stages in TC and 45 such miRNAs were found in YK, 34 in common; there were more down-regulated stage-DE miRNAs than up-regulated. And a total of 23, 30, 12, 6, and 30 breed-DE miRNAs between TC and YK were identified at 40, 55, 63, 70, and 90 dpc, respectively, which were mainly involved in cellular protein modification process, protein transport, and metabolic process. As the only highly expressed breed-DE miRNA found in no less than four developmental stages, and also a stage-DE miRNA found both in TC and YK, miR-499-5p could bind the 3’-UTR of a myofibrillogenesis regulator, destrin/actin depolymerizing factor (DSTN), as validated in dual luciferase reporter assay. The results suggested that miR-499-5p possibly play a noteworthy role in the breed-distinctive porcine muscle fiber development associated with the regulation of DSTN.

A new approach of gene co-expression network inference reveals significant biological processes involved in porcine muscle development in late gestation

The integration of genetic information in the cellular and nuclear environments is crucial for deciphering the way in which the genome functions under different physiological conditions. Experimental techniques of 3D nuclear mapping, a high-flow approach such as transcriptomic data analyses, and statistical methods for the development of co-expressed gene networks, can be combined to develop an integrated approach for depicting the regulation of gene expression. Our work focused more specifically on the mechanisms involved in the transcriptional regulation of genes expressed in muscle during late foetal development in pig. The data generated by a transcriptomic analysis carried out on muscle of foetuses from two extreme genetic lines for birth mortality are used to construct networks of differentially expressed and co-regulated genes. We developed an innovative co-expression networking approach coupling, by means of an iterative process, a new statistical method for graph inference with data of gene spatial co-localization (3D DNA FISH) to construct a robust network grouping co-expressed genes. This enabled us to highlight relevant biological processes related to foetal muscle maturity and to discover unexpected gene associations between IGF2, MYH3 and DLK1/MEG3 in the nuclear space, genes that are up-regulated at this stage of muscle development.

A transcriptomics resource reveals a transcriptional transition during ordered sarcomere morphogenesis in flight muscle

Muscles organise pseudo-crystalline arrays of actin, myosin and titin filaments to build force-producing sarcomeres. To study sarcomerogenesis, we have generated a transcriptomics resource of developing Drosophila flight muscles and identified 40 distinct expression profile clusters. Strikingly, most sarcomeric components group in two clusters, which are strongly induced after all myofibrils have been assembled, indicating a transcriptional transition during myofibrillogenesis. Following myofibril assembly, many short sarcomeres are added to each myofibril. Subsequently, all sarcomeres mature, reaching 1.5 µm diameter and 3.2 µm length and acquiring stretch-sensitivity. The efficient induction of the transcriptional transition during myofibrillogenesis, including the transcriptional boost of sarcomeric components, requires in part the transcriptional regulator Spalt major. As a consequence of Spalt knock-down, sarcomere maturation is defective and fibers fail to gain stretch-sensitivity. Together, this defines an ordered sarcomere morphogenesis process under precise transcriptional control - a concept that may also apply to vertebrate muscle or heart development.

Genome-wide profiling of Sus scrofa circular RNAs across nine organs and three developmental stages

The spatio-temporal expression patterns of Circular RNA (circRNA) across organs and developmental stages are critical for its function and evolution analysis. However, they remain largely unclear in mammals. Here, we comprehensively analysed circRNAs in nine organs and three skeletal muscles of Guizhou miniature pig (S. scrofa), a widely used biomedical model animal. We identified 5,934 circRNAs and analysed their molecular properties, sequence conservation, spatio-temporal expression pattern, potential function, and interaction with miRNAs. S. scrofa circRNAs show modest sequence conservation with human and mouse circRNAs, are flanked by long introns, exhibit low abundance, and are expressed dynamically in a spatio-temporally specific manner. S. scrofa circRNAs show the greatest abundance and complexity in the testis. Notably, 31% of circRNAs harbour well-conserved canonical miRNA seed matches, suggesting that some circRNAs act as miRNAs sponges. We identified 149 circRNAs potentially associated with muscle growth and found that their host genes were significantly involved in muscle development, contraction, chromatin modification, cation homeostasis, and ATP hydrolysis-coupled proton transport; moreover, this set of genes was markedly enriched in genes involved in tight junctions and the calcium signalling pathway. Finally, we constructed the first public S. scrofa circRNA database, allowing researchers to query comprehensive annotation, expression, and regulatory networks of circRNAs.

Dynamic transcriptome and DNA methylome analyses on longissimus dorsi to identify genes underlying intramuscular fat content in pigs

Background

The intramuscular fat content (IMF) refers to the amount of fat within muscles, including the sum of phospholipids mainly found in cell membranes, triglycerides and cholesterol, and is determined both by hyperplasia and hypertrophy of adipocyte during the development of pigs. The IMF content is an important economic trait that is genetically controlled by multiple genes. The Laiwu pig is an indigenous fatty pig breed distributed in North China, characterized by excessively higher level of IMF content (9%~12%), therefore, is suitable for the identification of genes controlling IMF variations. To identify genes underlying IMF deposition, we performed genome-wide transcriptome and methylome analyses on longissimus dorsi (LD) muscle in Laiwu pigs across four developmental stages.

Results

A total of 22,524 expressed genes were detected and 1158 differentially expressed genes (DEGs) were hierarchically clustered in the LD muscle over four developmental stages from 60 d to 400 d. These genes were significantly clustered into four temporal expression profiles, and genes participating in fat cell differentiation and lipid biosynthesis processes were identified. From 120 d to 240 d, the period with the maximum IMF deposition rate, the lipid biosynthesis related genes (FOSL1, FAM213B and G0S2), transcription factors (TFs) (EGR1, KLF5, SREBF2, TP53 and TWIST1) and enriched pathways (steroid biosynthesis and fatty acid biosynthesis) were revealed; and fat biosynthesis relevant genes showing differences in DNA methylation in gene body or intergenic region were detected, such as FASN, PVALB, ID2, SH3PXD2B and EGR1.

Conclusions

This study provides a comprehensive landscape of transcriptome of the LD muscle in Laiwu pigs ranging from 60 to 400 days old, and methylome of the LD muscle in 120 d and 240 d Laiwu pigs. A set of candidate genes and TFs involved in fat biosynthesis process were identified, which were probably responsible for IMF deposition. The results from this study would provide a reference for the identification of genes controlling IMF variation, and for exploring molecular mechanisms underlying IMF deposition in pigs.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4201-9) contains supplementary material, which is available to authorized users.

Dynamic transcriptomic analysis in hircine longissimus dorsi muscle from fetal to neonatal development stages

Muscle growth and development from fetal to neonatal stages consist of a series of delicately regulated and orchestrated changes in expression of genes. In this study, we performed whole transcriptome profiling based on RNA-Seq of caprine longissimus dorsi muscle tissue obtained from prenatal stages (days 45, 60, and 105 of gestation) and neonatal stage (the 3-day-old newborn) to identify genes that are differentially expressed and investigate their temporal expression profiles. A total of 3276 differentially expressed genes (DEGs) were identified (Q value < 0.01). Time-series expression profile clustering analysis indicated that DEGs were significantly clustered into eight clusters which can be divided into two classes (Q value < 0.05), class I profiles with downregulated patterns and class II profiles with upregulated patterns. Based on cluster analysis, GO enrichment analysis found that 75, 25, and 8 terms to be significantly enriched in biological process (BP), cellular component (CC), and molecular function (MF) categories in class I profiles, while 35, 21, and 8 terms to be significantly enriched in BP, CC, and MF in class II profiles. KEGG pathway analysis revealed that DEGs from class I profiles were significantly enriched in 22 pathways and the most enriched pathway was Rap1 signaling pathway. DEGs from class II profiles were significantly enriched in 17 pathways and the mainly enriched pathway was AMPK signaling pathway. Finally, six selected DEGs from our sequencing results were confirmed by qPCR. Our study provides a comprehensive understanding of the molecular mechanisms during goat skeletal muscle development from fetal to neonatal stages and valuable information for future studies of muscle development in goats.

Transcriptomic profiling in muscle and adipose tissue identifies genes related to growth and lipid deposition

Growth performance and meat quality are important traits for the pig industry and consumers. Adipose tissue is the main site at which fat storage and fatty acid synthesis occur. Therefore, we combined high-throughput transcriptomic sequencing in adipose and muscle tissues with the quantification of corresponding phenotypic features using seven Chinese indigenous pig breeds and one Western commercial breed (Yorkshire). We obtained data on 101 phenotypic traits, from which principal component analysis distinguished two groups: one associated with the Chinese breeds and one with Yorkshire. The numbers of differentially expressed genes between all Chinese breeds and Yorkshire were shown to be 673 and 1056 in adipose and muscle tissues, respectively. Functional enrichment analysis revealed that these genes are associated with biological functions and canonical pathways related to oxidoreductase activity, immune response, and metabolic process. Weighted gene coexpression network analysis found more coexpression modules significantly correlated with the measured phenotypic traits in adipose than in muscle, indicating that adipose regulates meat and carcass quality. Using the combination of differential expression, QTL information, gene significance, and module hub genes, we identified a large number of candidate genes potentially related to economically important traits in pig, which should help us improve meat production and quality.

Phenotypes, antioxidant responses, and gene expression changes accompanying a sugar-only diet in Bactrocera dorsalis (Hendel) (Diptera: Tephritidae)

BACKGROUND

Diet composition (yeast:carbohydrate ratio) is an important determinant of growth, development, and reproduction. Recent studies have shown that decreased yeast intake elicits numerous transcriptomic changes and enhances somatic maintenance and lifespan, which in turn reduces reproduction in various insects. However, our understanding of the responses leading to a decrease in yeast ratio to 0% is limited.

RESULTS

In the present study, we investigated the effects of a sugar-only diet (SD) on the gene expression patterns of the oriental fruit fly, Bactrocera dorsalis (Hendel), one of the most economically important pests in the family Tephritidae. RNA sequencing analyses showed that flies reared on an SD induced significant changes in the expression levels of genes associated with specific metabolic as well as cell growth and death pathways. Moreover, the observed upregulated genes in energy production and downregulated genes associated with reproduction suggested that SD affects somatic maintenance and reproduction in B. dorsalis. As expected, we observed that SD altered B. dorsalis phenotypes by significantly increasing stress (starvation and desiccation) resistance, decreasing reproduction, but did not extend lifespan compared to those that received a normal diet (ND) regime. In addition, administration of an SD resulted in a reduction in antioxidant enzyme activities and an increase in MDA concentrations, thereby suggesting that antioxidants cannot keep up with the increase in oxidative damage induced by SD regime.

CONCLUSIONS

The application of an SD diet induces changes in phenotypes, antioxidant responses, and gene expressions in B. dorsalis. Previous studies have associated extended lifespan with reduced fecundity. The current study did not observe a prolongation of lifespan in B. dorsalis, which instead incurred oxidative damage. The findings of the present study improve our understanding of the molecular, biochemical, and phenotypic response of B. dorsalis to an SD diet.

Phenotypic plasticity, trade-offs and gene expression changes accompanying dietary restriction and switches in Bactrocera dorsalis (Hendel) (Diptera: Tephritidae)

In this study, we investigated the effects of dietary restriction (DR) and variable diets on phenotypes and gene expression in oriental fruit fly, Bactrocera dorsalis (Hendel), one of the most economically important pests in the family Tephritidae around the world. As expected, we found that DR altered the B. dorsalis phenotypes by significantly increasing stress resistance and lifespan, but reduced egg production when compared with the control diet. The results suggested a trade-off between reproduction versus somatic maintenance (stress resistance) and lifespan in B. dorsalis. Diet also had a significant effect on hatchability, and DR could increase the egg hatching success of B. dorsalis. Furthermore, DR up-regulated metabolic pathways involved in energy homeostasis and down-regulated pathways in egg production, which might mediate trade-offs between somatic maintenance and reproduction under DR regimes. The gene expression profiles in response to the acute dietary switches indicated that the digestive and metabolic pathways maybe involved in the adaptability of flies to variable dietary resources. In summary, the research facilitates a better understanding of the molecular mechanisms responsible for the B. dorsalis’ phenotypic adjustments to the different qualities of the available diets.

Comparative analysis of DNA methylome and transcriptome of skeletal muscle in lean-, obese-, and mini-type pigs

DNA methylation plays a pivotal role in biological processes by affecting gene expression. However, how DNA methylation mediates phenotype difference of skeletal muscle between lean-, obese-, and mini-type pigs remains unclear. We systematically carried out comparative analysis of skeletal muscle by integrating analysis of genome-wide DNA methylation, mRNA, lncRNA and miRNA profiles in three different pig breeds (obese-type Tongcheng, lean-type Landrace, and mini-type Wuzhishan pigs). We found that the differentially methylated genes (DMGs) were significantly associated with lipid metabolism, oxidative stress and muscle development. Among the identified DMGs, 253 genes were related to body-size and obesity. A set of lncRNAs and mRNAs including UCP3, FHL1, ANK1, HDAC4, and HDAC5 exhibited inversely changed DNA methylation and expression level; these genes were associated with oxidation reduction, fatty acid metabolism and cell proliferation. Gene regulatory networks involved in phenotypic variation of skeletal muscle were related to lipid metabolism, cellular movement, skeletal muscle development, and the p38 MAPK signaling pathway. DNA methylation potentially influences the propensity for obesity and body size by affecting gene expression in skeletal muscle. Our findings provide an abundant information of epigenome and transcriptome that will be useful for animal breeding and biomedical research.

Developmental Stage, Muscle and Genetic Type Modify Muscle Transcriptome in Pigs: Effects on Gene Expression and Regulatory Factors Involved in Growth and Metabolism

Iberian pig production includes purebred (IB) and Duroc-crossbred (IBxDU) pigs, which show important differences in growth, fattening and tissue composition. This experiment was conducted to investigate the effects of genetic type and muscle (Longissimus dorsi (LD) vs Biceps femoris (BF)) on gene expression and transcriptional regulation at two developmental stages. Nine IB and 10 IBxDU piglets were slaughtered at birth, and seven IB and 10 IBxDU at four months of age (growing period). Carcass traits and LD intramuscular fat (IMF) content were measured. Muscle transcriptome was analyzed on LD samples with RNA-Seq technology. Carcasses were smaller in IB than in IBxDU neonates (p < 0.001), while growing IB pigs showed greater IMF content (p < 0.05). Gene expression was affected (p < 0.01 and Fold change > 1.5) by the developmental stage (5,812 genes), muscle type (135 genes), and genetic type (261 genes at birth and 113 at growth). Newborns transcriptome reflected a highly proliferative developmental stage, while older pigs showed upregulation of catabolic and muscle functioning processes. Regarding the genetic type effect, IBxDU newborns showed enrichment of gene pathways involved in muscle growth, in agreement with the higher prenatal growth observed in these pigs. However, IB growing pigs showed enrichment of pathways involved in protein deposition and cellular growth, supporting the compensatory gain experienced by IB pigs during this period. Moreover, newborn and growing IB pigs showed more active glucose and lipid metabolism than IBxDU pigs. Moreover, LD muscle seems to have more active muscular and cell growth, while BF points towards lipid metabolism and fat deposition. Several regulators controlling transcriptome changes in both genotypes were identified across muscles and ages (SIM1, PVALB, MEFs, TCF7L2 or FOXO1), being strong candidate genes to drive expression and thus, phenotypic differences between IB and IBxDU pigs. Many of the identified regulators were known to be involved in muscle and adipose tissues development, but others not previously associated with pig muscle growth were also identified, as PVALB, KLF1 or IRF2. The present study discloses potential molecular mechanisms underlying phenotypic differences observed between IB and IBxDU pigs and highlights candidate genes implicated in these molecular mechanisms.

Gene expression profile indicates involvement of NO in Camellia sinensis pollen tube growth at low temperature

BACKGROUND

Nitric oxide (NO) functions as a critical signaling molecule in the low-temperature stress responses in plants, including polarized pollen tube growth in Camellia sinensis. Despite this, the potential mechanisms underlying the participation of NO in pollen tube responses to low temperature remain unclear. Here, we investigate alterations to gene expression in C. sinensis pollen tubes exposed to low-temperature stress and NO using RNA-Seq technology, in order to find the potential candidate genes related to the regulation of pollen tube elongation by NO under low-temperature stress.

RESULTS

Three libraries were generated from C. sinensis cv. 'Longjingchangye' pollen tubes cultured at 25 °C (CsPT-CK) and 4 °C (CsPT-LT) or with 25 μM DEA NONOate (CsPT-NO). The number of unigenes found for the three biological replications were 39,726, 40,440 and 41,626 for CsPT-CK; 36,993, 39,070 and 39,439 for CsPT-LT; and 39,514, 38,298 and 39,061 for CsPT-NO. A total of 36,097 unique assembled and annotated sequences from C. sinensis pollen tube reads were found in a BLAST search of the following databases: NCBI non-redundant nucleotide, Swiss-prot protein, Kyoto Encyclopedia of Genes and Genomes, Cluster of Orthologous Groups of proteins, and Gene Ontology. The absolute values of log2Ratio > 1 and probability > 0.7 were used as the thresholds for significantly differential gene expression, and 766, 497 and 929 differentially expressed genes (DEGs) were found from the comparison analyses of the CK-VS-LT, CK-VS-NO and LT-VS-NO libraries, respectively. Genes related to metabolism and signaling pathways of plant hormones, transcription factors (TFs), vesicle polarized trafficking, cell wall biosynthesis, the ubiquitination machinery of the ubiquitin system and species-specific secondary metabolite pathways were mainly observed in the CK-VS-LT and CK-VS-NO libraries.

CONCLUSION

Differentially expressed unigenes related to the inhibition of C. sinensis pollen tube growth under low temperature and NO are identified in this study. The transcriptomic gene expression profiles present a valuable genomic tool to improve studying the molecular mechanisms underlying low-temperature tolerance in pollen tube.

Genome-wide identification and characterization of long non-coding RNAs in developmental skeletal muscle of fetal goat

Background

Long non-coding RNAs (lncRNAs) have been studied extensively over the past few years. Large numbers of lncRNAs have been identified in mouse, rat, and human, and some of them have been shown to play important roles in muscle development and myogenesis. However, there are few reports on the characterization of lncRNAs covering all the development stages of skeletal muscle in livestock.

Results

RNA libraries constructed from developing longissimus dorsi muscle of fetal (45, 60, and 105 days of gestation) and postnatal (3 days after birth) goat (Capra hircus) were sequenced. A total of 1,034,049,894 clean reads were generated. Among them, 3981 lncRNA transcripts corresponding to 2739 lncRNA genes were identified, including 3515 intergenic lncRNAs and 466 anti-sense lncRNAs. Notably, in pairwise comparisons between the libraries of skeletal muscle at the different development stages, a total of 577 transcripts were differentially expressed (P < 0.05) which were validated by qPCR using randomly selected six lncRNA genes. The identified goat lncRNAs shared some characteristics, such as fewer exons and shorter length, with the lncRNAs in other mammals. We also found 1153 lncRNAs genes were neighbored 1455 protein-coding genes (<10 kb upstream and downstream) and functionally enriched in transcriptional regulation and development-related processes, indicating they may be in cis-regulatory relationships. Additionally, Pearson’s correlation coefficients of co-expression levels suggested 1737 lncRNAs and 19,422 mRNAs were possibly in trans-regulatory relationships (r > 0.95 or r < −0.95). These co-expressed mRNAs were enriched in development-related biological processes such as muscle system processes, regulation of cell growth, muscle cell development, regulation of transcription, and embryonic morphogenesis.

Conclusions

This study provides a catalog of goat muscle-related lncRNAs, and will contribute to a fuller understanding of the molecular mechanism underpinning muscle development in mammals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3009-3) contains supplementary material, which is available to authorized users.

Transcriptome analysis revealed chimeric RNAs, single nucleotide polymorphisms and allele-specific expression in porcine prenatal skeletal muscle

Prenatal skeletal muscle development genetically determines postnatal muscle characteristics such as growth and meat quality in pigs. However, the molecular mechanisms underlying prenatal skeletal muscle development remain unclear. Here, we performed the first genome-wide analysis of chimeric RNAs, single nuclear polymorphisms (SNPs) and allele-specific expression (ASE) in prenatal skeletal muscle in pigs. We identified 14,810 protein coding genes and 163 high-confidence chimeric RNAs expressed in prenatal skeletal muscle. More than 94.5% of the chimeric RNAs obeyed the canonical GT/AG splice rule and were trans-splicing events. Ten and two RNAs were aligned to human and mouse chimeric transcripts, respectively. We detected 106,457 high-quality SNPs (6,955 novel), which were mostly (89.09%) located within QTLs for production traits. The high proportion of non-exonic SNPs revealed the incomplete annotation status of the current swine reference genome. ASE analysis revealed that 11,300 heterozygous SNPs showed allelic imbalance, whereas 131 ASE variants were located in the chimeric RNAs. Moreover, 4 ASE variants were associated with various economically relevant traits of pigs. Taken together, our data provide a source for studies of chimeric RNAs and biomarkers for pig breeding, while illuminating the complex transcriptional events underlying prenatal skeletal muscle development in mammals.

Combined Cytological and Transcriptomic Analysis Reveals a Nitric Oxide Signaling Pathway Involved in Cold-Inhibited Camellia sinensis Pollen Tube Growth

Nitric oxide (NO) as a signaling molecule plays crucial roles in many abiotic stresses in plant development processes, including pollen tube growth. Here, the signaling networks dominated by NO during cold stress that inhibited Camellia sinensis pollen tube growth are investigated in vitro. Cytological analysis show that cold-induced NO is involved in the inhibition of pollen tube growth along with disruption of the cytoplasmic Ca(2+) gradient, increase in ROS content, acidification of cytoplasmic pH and abnormalities in organelle ultrastructure and cell wall component distribution in the pollen tube tip. Furthermore, differentially expressed genes (DEGs)-related to signaling pathway, such as NO synthesis, cGMP, Ca(2+), ROS, pH, actin, cell wall, and MAPK cascade signal pathways, are identified and quantified using transcriptomic analyses and qRT-PCR, which indicate a potential molecular mechanism for the above cytological results. Taken together, these findings suggest that a complex signaling network dominated by NO, including Ca(2+), ROS, pH, RACs signaling and the crosstalk among them, is stimulated in the C. sinensis pollen tube in response to cold stress, which further causes secondary and tertiary alterations, such as ultrastructural abnormalities in organelles and cell wall construction, ultimately resulting in perturbed pollen tube extension.

Comparative Analysis of Muscle Transcriptome between Pig Genotypes Identifies Genes and Regulatory Mechanisms Associated to Growth, Fatness and Metabolism

Iberian ham production includes both purebred (IB) and Duroc-crossbred (IBxDU) Iberian pigs, which show important differences in meat quality and production traits, such as muscle growth and fatness. This experiment was conducted to investigate gene expression differences, transcriptional regulation and genetic polymorphisms that could be associated with the observed phenotypic differences between IB and IBxDU pigs. Nine IB and 10 IBxDU pigs were slaughtered at birth. Morphometric measures and blood samples were obtained and samples from Biceps femoris muscle were employed for compositional and transcriptome analysis by RNA-Seq technology. Phenotypic differences were evident at this early age, including greater body size and weight in IBxDU and greater Biceps femoris intramuscular fat and plasma cholesterol content in IB newborns. We detected 149 differentially expressed genes between IB and IBxDU neonates (p < 0.01 and Fold-Change > 1. 5). Several were related to adipose and muscle tissues development (DLK1, FGF21 or UBC). The functional interpretation of the transcriptomic differences revealed enrichment of functions and pathways related to lipid metabolism in IB and to cellular and muscle growth in IBxDU pigs. Protein catabolism, cholesterol biosynthesis and immune system were functions enriched in both genotypes. We identified transcription factors potentially affecting the observed gene expression differences. Some of them have known functions on adipogenesis (CEBPA, EGRs), lipid metabolism (PPARGC1B) and myogenesis (FOXOs, MEF2D, MYOD1), which suggest a key role in the meat quality differences existing between IB and IBxDU hams. We also identified several polymorphisms showing differential segregation between IB and IBxDU pigs. Among them, non-synonymous variants were detected in several transcription factors as PPARGC1B and TRIM63 genes, which could be associated to altered gene function. Taken together, these results provide information about candidate genes, metabolic pathways and genetic polymorphisms potentially involved in phenotypic differences between IB and IBxDU pigs associated to meat quality and production traits.

Polymorphisms in the regulatory region of the porcine MYLPF gene are related to meat quality traits in the Large White breed

The MYLPF gene encodes fast myosin regulatory light chain, and is a positional and functional candidate gene for meat quality. The aim of this study was to identify associations between SNPs in the promoter region of the porcine MYLPF gene and meat quality traits. A total of 22 SNPs were identified in a population of crossbred animals (n=86) and based on minor allele frequency and proximity to the transcription start site, five SNPs were genotyped in purebred; Large White (n=98), Duroc (n=99) and Pietrain (n=98) pigs. No associations were observed in the Pietrain breed, while the Duroc breed was almost monomorphic for all SNPs. In the Large White breed SNP g-1314A>G and linked SNPS g.-871T>G, g.-566T>C, g.-403C>G were associated with ultimate pH and driploss (P<0.05). This study identified associations between MYLPF and meat quality and highlights the importance of considering the genetic background within gene-assisted selection programmes.

Identification of Genes Related to Growth and Lipid Deposition from Transcriptome Profiles of Pig Muscle Tissue

Transcriptome profiles established using high-throughput sequencing can be effectively used for screening genome-wide differentially expressed genes (DEGs). RNA sequences (from RNA-seq) and microRNA sequences (from miRNA-seq) from the tissues of longissimus dorsi muscle of two indigenous Chinese pig breeds (Diannan Small-ear pig [DSP] and Tibetan pig [TP]) and two introduced pig breeds (Landrace [LL] and Yorkshire [YY]) were examined using HiSeq 2000 to identify and compare the differential expression of functional genes related to muscle growth and lipid deposition. We obtained 27.18 G clean data through the RNA-seq and detected that 18,208 genes were positively expressed and 14,633 of them were co-expressed in the muscle tissues of the four samples. In all, 315 DEGs were found between the Chinese pig group and the introduced pig group, 240 of which were enriched with functional annotations from the David database and significantly enriched in 27 Gene Ontology (GO) terms that were mainly associated with muscle fiber contraction, cadmium ion binding, response to organic substance and contractile fiber part. Based on functional annotation, we identified 85 DEGs related to growth traits that were mainly involved in muscle tissue development, muscle system process, regulation of cell development, and growth factor binding, and 27 DEGs related to lipid deposition that were mainly involved in lipid metabolic process and fatty acid biosynthetic process. With miRNA-seq, we obtained 23.78 M reads and 320 positively expressed miRNAs from muscle tissues, including 271 known pig miRNAs and 49 novel miRNAs. In those 271 known miRNAs, 20 were higher and 10 lower expressed in DSP-TP than in LL-YY. The target genes of the 30 miRNAs were mainly participated in MAPK, GnRH, insulin and Calcium signaling pathway and others involved cell development, growth and proliferation, etc. Combining the DEGs and the differentially expressed (DE) miRNAs, we drafted a network of 46 genes and 18 miRNAs for regulating muscle growth and a network of 15 genes and 16 miRNAs for regulating lipid deposition. We identified that CAV2, MYOZ2, FRZB, miR-29b, miR-122, miR-145-5p and miR-let-7c, etc, were key genes or miRNAs regulating muscle growth, and FASN, SCD, ADORA1, miR-4332, miR-182, miR-92b-3p, miR-let-7a and miR-let-7e, etc, were key genes or miRNAs regulating lipid deposition. The quantitative expressions of eight DEGs and seven DE miRNAs measured with real-time PCR certified that the results of differential expression genes or miRNAs were reliable. Thus, 18,208 genes and 320 miRNAs were positively expressed in porcine longissimus dorsi muscle. We obtained 85 genes and 18 miRNAs related to muscle growth and 27 genes and 16 miRNAs related to lipid deposition, which provided new insights into molecular mechanism of the economical traits in pig.

Comparison of skeletal muscle miRNA and mRNA profiles among three pig breeds

The pig is an important source of animal protein, and is also an ideal model for human disease. There are significant differences in growth rate, muscle mass, and meat quality between different breeds. To understand the molecular mechanisms underlying porcine skeletal muscle phenotypes, we performed mRNA and miRNA profiling of muscle from three different breeds of pig, Landrace (lean-type), Tongcheng (obese-type), and Wuzhishan (mini-type) by Solexa sequencing. Forty-three genes and 106 miRNAs were differentially expressed between Landrace and Tongcheng pigs, 92 genes and 151 miRNAs were differentially expressed between Tongcheng and Wuzhishan pigs, and 145 genes and 156 miRNAs were differential expressed between Landrace and Wuzhishan pigs. Gene ontology analysis suggested that genes differentially expressed between Landrace and Tongcheng pigs were mainly involved in the biological processes of oxidative stress and muscle organ development. Meanwhile, for Tongcheng vs Wuzhishan and Landrace vs Wuzhishan pigs, the differentially expressed genes were involved in fatty acid metabolism, oxidative stress, muscle contraction, and muscle organ development, processes that are closely related to meat quality. To investigate the molecular mechanisms underlying meat quality diversity based on differentially expressed genes and miRNAs, interaction networks were constructed, according to target prediction results and integration analysis of up-regulated genes with down-regulated miRNAs or down-regulated genes with up-regulated miRNAs. Our findings identify candidate genes and miRNAs associated with muscle development and indicate their potential roles in muscle phenotype variance between different pig breeds. These results serve as a foundation for further studies on muscle development and molecular breeding.