Skip to main content

Integrative analysis of transcriptomics and proteomics of skeletal muscles of the Chinese indigenous Shaziling pig compared with the Yorkshire breed



The Shaziling pig (Sus scrofa) is a well-known indigenous breed in China. One of its main advantages over European breeds is its high meat quality. However, little genetic information is available for the Shaziling pig. To screen for differentially expressed genes and proteins that might be responsible for the meat quality, the longissimus dorsi muscles from Shaziling and Yorkshire pig breeds were investigated using an integrative analysis of transcriptomics and proteomics, involving high-throughput sequencing, the two-dimensional gel electrophoresis, and mass spectrometry.


Sequencing produced 79,320 unigenes by de novo assembly, and 488 differentially expressed genes in the longissimus dorsi muscle of Shaziling pig compared with the Yorkshire breed were identified. Gene Ontology term enrichment of biological functions and Kyoto Encyclopedia of Genes and Genomes analysis showed that the gene products were mainly involved in metabolism, protein binding, and regulation of skeletal muscle development. At the protein level, 23 differentially expressed proteins were identified, which were potentially associated with fatty acid metabolism, the glycolytic pathway, and skeletal muscle growth. Eight differentially expressed genes were confirmed by real-time PCR. These results give an insight into the mechanisms underlying the formation of skeletal muscle in the Shaziling pig.


Certain differentially expressed genes and proteins are involved in fatty acid metabolism, intramuscular fat deposition, and skeletal muscle growth in the Shaziling pig. These results provide candidate genes for improving meat quality and will promote further transcriptomic research in Shaziling pigs.


Pigs are important domestic animals used for meat production worldwide. Pork quality is influenced by many factors, including breed, nutrition, and post-slaughter handling [1]. Among these factors, the breed is the most important. There are more pig breeds in China than in any other country in the word [2]. In China, 118 indigenous pig breeds are listed as the Domestic Animal Diversity in the World index [3]. Over a long period, breeders have attempted to increase muscle yield and decrease carcass fatness, and great progress for these traits has been made in swine breeding. For example, Landrace pigs and Yorkshire pigs grow more quickly and have higher lean meat than other pig breeds. However, some studies suggest that such intensive selection for increased muscle growth and decreased carcass fatness has led to a deterioration in meat quality [46]. Compared with European breeds, Chinese native breeds have higher intramuscular fat (IMF), and increased tenderness and meat quality [710]. The Shaziling pig is a well-known breed, which has evolved for centuries in Hunan Province, China, where it feeds on grains, tubers, and wild herbs. It is a Chinese fat-type line with high IMF, superior meat quality, and strong resistance to general diseases. Compared with the Shaziling breed, the Yorkshire has a faster growth rate and a leaner meat ratio.

Skeletal muscle contains several fiber types [11]. Muscle fiber types and the proportion of fiber types affect meat quality directly [12]. Therefore, studies on development and growth are beneficial to improve meat quality. Skeletal muscle development is very complicated, and comprises several stages: determination of myoblasts, proliferation of myoblasts, differentiation and fusion of myoblasts into myotubes and myofibers, and growth and maturation [13, 14]. Analyzing these stages would provide a good basis for understanding muscle fiber development. Previous reports showed that muscle fibers could be classified into red and white fibers [15]. These compositional differences between fibers determine their distinct metabolic type and physiological functions and affect meat quality [16, 17]. Skeletal muscles have been explored extensively using molecular biology [18, 19], and proteomics and transcriptomics techniques have been applied to study the porcine skeletal muscle of different breeds [20, 21].

Over the past decade, proteomic technologies have been used successfully to study skeletal muscle [11, 22]. Proteomic analysis based on two-dimensional gel electrophoresis (2-DE) and mass spectrometry (MS) is a classical method in quantitative proteomics to separate mixtures of proteins into two dimensions and has been a powerful tool in meat science [23, 24]. Many reports concerning differential proteomics among different pig breeds have been published [2527]. Another technology for characterizing molecular changes in skeletal muscle is analysis of the transcriptome. Recently, transcriptome studies have been applied to analyze differentially expressed genes (DEGs), identify novel genes, describe metabolic pathways, and forecast the relationship between genotypes and phenotypes [2830]. Next-generation sequencing technology has provided a new tool to quantify transcriptomes and analyze gene expression on a global scale. Transcriptional and proteomic methods could be used to analyze the changes from the mRNA expression to the protein abundance. In addition, post-transcriptional regulation is very important for mRNA stability, translation initiation, and protein stability [31]. Thus, it is necessary to combine proteomic and transcriptional methods simultaneously to analyze skeletal muscle growth and development. Currently, several reports have presented results of both proteomic and transcriptional analyses. For instance, longissimus dorsi muscle (LM) proteome and transcriptome profiles of Yorkshire pig and Casertana pig breeds were compared using 2-DE and a microarray. As a result, a large number of genes were identified that are involved in glycolytic metabolism and skeletal muscle growth [32]. In addition, Kim et al. [33] also compared the LM proteome and transcriptome profiles of different pig breeds using 2-DE and a microarray.

In the present study, we performed transcriptomic and proteomic analyses, along with functional enrichment of Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, to characterize the expression profiles in the LM of Shaziling and Yorkshire pigs. The aim of this study was to reveal the differences in breed-related protein and transcript expression patterns between the two breeds. These results could provide an insight into the mechanisms of growth and development of porcine skeletal muscle.


Sample preparation

Three 25-day-old, male full-sibs for each of Shaziling pigs and Yorkshire pigs were slaughtered following the Animal Care Guidelines of the Ethics committee of Hunan Agricultural University. Samples of LM were collected from the same area of the ribs and placed in liquid nitrogen immediately. All samples were kept at −80 °C after collection until used.

Protein extraction

The frozen LM tissues (100 mg) from Shaziling pigs and Yorkshire pigs were crushed in a mortar containing liquid nitrogen and were then sonicated for 10 s using a Sonoplus (Bandelin Electronic, Berlin Germany). The crushed tissue was homogenized in 1 mL of cold dissolution buffer containing 7 M urea, 2 M thiourea, 1 % dithiothreitol (DTT), 4 % (w/v) CHAPS, 20 μL protease inhibitor cocktail (BBI, Canada), and 2 % (v/v) pharmalyte (pH 3–10; BioRad, Hercules, CA, USA). The homogenate was centrifuged at 15,000 × g, for 20 min at 4 °C. The supernatant fraction was filtered and kept at −80 °C for subsequent analysis. The total protein content was determined using a Bradford assay kit (Bio-Rad).

2-DE and images analysis of gel

Approximately 100 μg of extracted protein was diluted with rehydration buffer [8 M urea, 2 M thiourea, 50 mM DTT, 4 % (w/v) CHAPS, and 0.5 % carrier ampholytes (pH 3–10, Bio-Rad)]. The mix was loaded onto 13-cm, nonlinear, immobilized pH gradient strips (IPG, pH 3–10; BioRad), which were rehydrated overnight. After rehydration, first-dimension isoelectric focusing was carried out at 20 °C at 30 V for 12 h, 500 V for 1 h, 1000 V for 1 h, and then gradually increased to 8,000 V. Isoelectric focusing was performed using an Ettan IPGphor III system (GE Healthcare, USA) gel apparatus. The focused IPG strips were incubated for 15 min in equilibration buffer containing 6 M urea, 50 mM Tris–HCl, pH 8.8, 30 % glycerol, 2 % sodium dodecyl sulfate (SDS), and 1 % DTT. Then, strips were then incubated again for 15 min in a buffer containing 6 M urea, 50 mM Tris–HCl, pH 8.8, 2 % SDS, 30 % glycerol, and 4 % iodoacetamide. In the second dimension, the equilibrated IPG strips were placed onto SDS polyacrylamide gels (12.5 % T, 2.6 % C) for protein separation using the Ettan DALT six electrophoresis system (GE Healthcare). When the bromophenol blue dye front reached the bottom of the gel, electrophoresis was stopped and analytical gels were subjected to silver staining and Coomassie blue was used to stain preparative gels before identification by MS. Each sample was analyzed in triplicate.

Spot choosing and tryptic digestion

Gel images were scanned using an Image Scanner UMax Powerlook 2110XL (GE Amersham) and analyzed using Image Master 2D Platinum software Hofer SE 600 (Version 5.0; GE Amersham). The protein spots were compared automatically using the software, including matching and quality. Spots whose average density was different by more than 1.5-fold between the two pig breeds were analyzed and marked. The differentially abundant spots were cut out from the preparative gel carefully and washed twice with 200 μL of 50 mM ammonium bicarbonate with 50 % acetonitrile and incubated for 15 min at room temperature. The gel pieces were swollen in a digestion solution containing 5 μL of 25 mM ammonium bicarbonate and 10 ng of trypsin at 4 °C for 30 min; in-gel tryptic digestion was run overnight at 37 °C. Subsequently, the supernatant was extracted twice with solvent A (80 % acetonitrile, 0.1 % trifluoroacetic acid) for 15 min at 3 °C. Finally, the digested tryptic peptides were passed through a Zip-Tip to remove salts, according to the manufacturer’s protocol.

Mass spectrometry

The samples were stored at −70 °C before analysis by MALDI-TOF/TOF. Protein identification was performed using a 5800 MALDI-TOF/TOF mass spectrometer (AB SCIEX) according to the manufacturer’s instructions. Mass spectra were acquired in reflector mode, and recorded in the range of 800–4,500 Da. Eight of the most intense ion signals were selected as precursors for the acquisition of MS/MS. The resulting peptide masses were submitted into the database of the National Center for Biotechnology Information non-redundant (NCBI nr) and the Swiss-Prot database using the Mascot server (Matrix Science, London, UK) to identify proteins. The search parameters were set at ± 100 ppm for peptide-mass mapping (PMF), peptide tolerance and ± 0.4 Da for the MS/MS tolerance.

RNA isolation and transcriptome analysis

Total RNA was extracted from the LM of the two breeds of the pig using the total RNA extraction kit (Qiagen, Valencia, CA, USA), in accordance with the manufacturer’s instructions. The RNA quantity and integrity were checked using a NanoDrop 2000 spectrophotometer and bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). A TruSeq RNA Sample Preparation kit v2 (Illumina, San Diego, CA, USA) was used to construct cDNA libraries. Subsequently, the libraries were sequenced using an Illumina HiSeq 2500 instrument (Illumina, San Diego, CA, USA) that generated paired-end reads of 100 bp.

De novo assembling and functional annotation of reads

Raw sequencing reads of each sample were trimmed and assembled de novo using CLC Genomics Workbench (CLC Bio, Aarhus, Denmark). After the adapter sequences, ambiguous bases and sequences less than 20 nucleotides were removed, credible contigs generated by de novo sequence assembly. The assembled contigs were annotated against the UniProt and NCBI non-redundant (nr) database using BlastX alignment with an E-value cut off of <1e-5. Based on the annotation results, GO terms were extracted using Blast2GO ( and the results were classified as biological processes, cellular components, and molecular functions. The EuKaryotic Orthologous Groups (KOG) and KEGG databases were used to predict the functions and define the main metabolic pathways, separately.

Gene expression quantification and differential expression analysis

The expression level of each gene was calculated using the reads per kb per million reads values by the Qualimap v0.5 software [34]. The DEGseq program [35] and R packages were used to identify significantly DEGs between Shaziling pigs and Yorkshire pigs. In this study, the criteria were a fold change (FC) greater than two and cut-off of the false discovery rate of 5 %. For the unigenes that were considered as having differential expression, GO functional enrichment was carried out when the P value was less than 0.05. DEGseq provides statistical routines for determining DEGs.

Validation of DEGs by quantitative real-time PCR (qPCR)

To validate the expression level of the DEGs, eight DEGs from 2-DE results were analyzed by using qPCR. These genes included Enolase 3 (ENO3); ATP synthase subunit beta, mitochondrial (ATP5B); Enolase 1 (ENO1); Triosephosphate Isomerase (TPI); alpha-actin (ACTA1); actin, alpha, cardiac muscle 1 (ACTC1); myosin light chain, phosphorylatable, fast skeletal muscle (MYLPF); and heat shock protein 27 (HSPB1), which mainly belonged to metabolic pathways and structural proteins. Total RNA was extracted from the LM of Shaziling and Yorkshire pigs using the Trizol reagent (Invitrogen, CA, USA), according to the manufacturer’s instructions, and the RNA was treated with RNase-free DNase (Qiagen, Valencia, CA, USA). The quality and concentration of the total RNA was evaluated by gel electrophoresis. The cDNA synthesis of samples was performed using the Quantitect Reverse Transcription kit (Qiagen). The primers (Table 1) were designed by the Roche Universal Probe Library Assay Design Center ( and synthesized by Sangon Biotechnology Corporation (Shanghai, China). qPCR was performed on an Applied Biosystems 7300 Real time PCR system (Applied Biosystems, USA) with the SYBR Premix Ex Taq Kit (Takara Biotechnology, Japan) in a 20-μL PCR mix. GAPDH was used as the reference gene to measure the expression levels of mRNA between the samples and data were calculated by the 2−ΔΔCt comparative CT method [36].

Table 1 Primer sequences for the quantitative real-time PCR amplification of the differential expressed genes in Shaziling and Yorkshire pigs

Results and discussion

Identification of differentially abundant proteins by 2-DE and MS

The Longissimus dorsi proteomes from the two pig breeds were analyzed by the 2-DE technique. Thirty-eight protein spots (Fig. 1) were significantly different (FC >1.5) between the two breeds: 27 protein spots were upregulated in Shaziling pigs and 11 in Yorkshire. The differentially abundant protein spots were identified by MALDI-TOF/TOF-MS: two spots could not be identified successfully; however, the 36 remaining spots were identified by matching peptide data to the UniProt database. The characteristics of the identified proteins and the identification parameters are listed in Table 2. Some of the identified proteins were resolved as multiple spots. For example, spots 2480, 2654, and 2659 represented the same protein and showed similar quantitative trends. Ultimately, 23 different proteins were identified.

Fig. 1
figure 1

2-DE map of longissimus dorsi musle from Yorkshire and Shaziling pig breeds. Scanned 2-DE image of separated using an IPG pH 3–10 strip in the first dimension (12 cm, BioRad, USA), and 12.5 % SDS gel in the second dimension. Tag show 38 spots that were significantly changed between the two breeds [a Yorkshire pigs (up-regulation), b Shaziling pigs (up-regulation)]

Table 2 Protein differentially expressed for Shaziling pig breeds and Yorkshire were identified by 2-DE and MALDI-TOF-MS/MS

Function analysis of differentially expressed proteins

To further determine the functions of the differentially abundant proteins, functional categorization of all the identified proteins was performed using GO annotation. The annotated results of biological process are shown in Fig. 2. The 23 differential proteins were grouped into 11 categories. We focused on the categories of fatty acid metabolism, glycolytic pathway, and skeletal muscle growth and development. The result of enrichment analysis revealed that some of the identified proteins had an important impact on meat quality. For example, TPI plays an important role in ensuring immediate equilibration of the triosephosphates produced by aldolase in glycolysis, which is interconnected with lipid metabolism and to the pentose phosphate pathway [37]. In addition, a previous report showed that a quantitative increase in TPI was positively correlated with more tender meat [38].

Fig. 2
figure 2

GO annotation of different proteinic spots (a: molecular function distribution, b: Biological process distribution)

Enolase 3 was another metabolic enzyme that was differentially abundant between the two breeds. Enolase 3 is a crucial enzyme in glycolysis that catalyzes the interconversion of diphosphoglycerate and phosphoenolpyruvate. Pig muscle Enolase 3 was investigated some time ago [39]. Many isoforms of Enolase 3 have been confirmed to have an influence on IMF in pigs [20, 40].

Muscles are composed primarily of different muscle fibers. Muscle fiber type is an important factor influencing meat quality [41]. For example, increasing the percentage of type IIb fibers could lead to the reduced meat quality because of altered metabolic rates and biochemical processes [42]. ACTA1 is a member of the actin family and is a major constituent of the contractile apparatus in skeletal muscle. Previous reports showed that alpha-actin levels positively correlate with the synthesis of muscle fiber proteins and muscle growth [43]. In this study, ACTA1 was differentially abundant between the two pig breeds.

Heat shock proteins (HSPs) including HSPA1, HSPA8, HSPB1, and other chaperone proteins have been associated with pig meat quality [44, 45]. In the present study, HSPB1 was more abundant in the Shaziling pig compared with the Yorkshire pig. This observation suggested that HSPB1 is correlated with meat quality. In addition to tenderness, further research is needed to confirm whether the IMF content of meat is affected by the expression of HSPB1.

Illumina sequencing and de novo assembly

Sequencing generated 30,751,992 and 35,933,274 100-bp paired-end reads from cDNA libraries of Shaziling and Yorkshire pigs, respectively. After trimming of low-quality reads and the removal of adapters, 29,522,984 and 34,937,710 reads were used for de novo assembly. The reads were assembled into 86,759 contigs (N50 = 939 bp) ranging from 100 to 49,881 bp, with an average length of 672 bp. Thereafter, the contigs were assembled into scaffolds (N50 = 1028 bp) with a mean length of 713 bp. After the final assembly of scaffolds using CAP3 [46], we obtained 79,320 unigenes (N50 = 1112 bp) with an average length of 733 bp and a maximum of 66,767 bp. The length distribution of the assembled unigenes is shown in Fig. 3. The result of the assembly is shown in Table 3. The transcriptome data have been submitted to the NCBI GEO database under accession number GSE70673.

Fig. 3
figure 3

Length distribution of the assembled final unigenes of Illumina sequences

Table 3 Statistical summary of the longissimus dorsi muscle transcriptome for assembling

Functional annotation of longissimus dorsi transcripts

To validate and annotate the protein functions of the 79,320 assembled unigenes, they were searched against the NCBI nr and UniProt protein database using the BLAST program (E-value cutoff <1e−5). Of the 79,320 assembled unigenes, 23,055 sequences (29.07 %) were assigned to the nr database, and 25,784 sequences (32.50 %) were found to have significant similarity with proteins in the UniProt database (Table 4). Approximately 30 % of all assembled unigenes was perfectly matched to the Uniprot database and the nr database, respectively. Thus, most of the assembled unigenes were unknown, indicating that many new genes and non-coding RNA sequences were obtained. Generally, direct sequencing of transcriptomes is considered an effective tool to discover new genes [47, 48] and has been applied to various organisms [49, 50]. These unigenes provided much valuable information for further identification of genes related to development of the longissimus dorsi.

Table 4 Functional annotations using transcript BLAST analyses

Functional classification of unigenes

GO is an extremely useful tool to classify the functions of a large quantity of genes, and is an international standardized gene functional classification system. GO analysis has been widely applied to predict the functions of genes in many organisms. Three ontologies, cellular component, molecular function, and biological process, are shared in the GO database. In this study, of the annotated unigenes, 23,702 could be assigned to the biological process (BP, 111,204 sequences), cellular components (CC, 68,072), and molecular functions (MF, 30,087) (Fig. 4) using the Blast2GO program [51]. In the biological process, the majority of unigenes were involved in catalytic activity and transcription factor activity, indicating that they might play a key role in the regulation of skeletal muscle development. In cellular components, organelle, membrane, cell, and cell part were prominently represented. For molecular functions, “catalytic activity” and “binding” represented the majority of the category.

Fig. 4
figure 4

Column chart presentation of GO classification of unigenes

To further evaluate the function of the unigenes, KOG classification of all the unigenes was performed, and 16,171 unigenes could be assigned to 25 categories (Fig. 5). Among the 25 KOG categories, the highest proportion of unigenes matched to “signal transduction mechanisms” (13.8 %), followed by cluster for “general function prediction only” (12.6 %), “post-translational modification, protein turnover, chaperones” (9.8 %), and “carbohydrate transport and metabolism” (5.8 %). The KEGG pathway database is a powerful tool for the analysis of gene function in terms of gene networks [52]. To identify the biological pathways of transcriptome from the LM, all unigenes were submitted to KEGG pathway analysis. A total of 16,755 unigenes were assigned to 278 KEGG pathways (Additional file 1: Table S1), which fell into five main categories, including metabolism, environmental information processing, genetic information processing, cellular processes, and organism systems. Among the five main categories, the metabolism pathway represented the largest number of the unigenes (1013). These metabolic pathways were mainly carbohydrate metabolism, energy metabolism, amino acid metabolism, and lipid metabolism, which provided a valuable resource to investigate muscle growth and lipid metabolism processes.

Fig. 5
figure 5

KOG functional classification of transcriptome

Identification and analysis of DEGs

To identify the DEGs between the Shaziling and Yorkshire pigs, the relative expression of the genes was estimated using the RPKM value, which is a common method for estimating transcript levels [53]. According to the RNA-seq analysis, DEGs were selected using the criteria of FC ≥2 and false discovery rate ≤0.05; the FC distribution of DEGs is shown in Fig. 6. A total of 488 unigenes were identified as DEGs between the two breeds, of which 297 were upregulated in the Shaziling pig and the 192 genes were upregulated in the Yorkshire pig. To better explore the functions of the DEGs, GO functional analysis was carried out. The results indicated that 488 DEGs were significantly enriched in 208 GO categories (Additional file 2: Table S2). Notably, some of the categories are involved in the metabolic process and the regulation of skeletal muscle development.

Fig. 6
figure 6

The fold change distribution of up- and down-regulated DEGs. Green bars refer to down-regulated DEGs and red bars refer to up-regulated DEGs in Shaziling pigs compared with Yorkshire pigs. The X axis shows fold change of DEGs and the Y axis number of DEGs

qPCR analysis

To evaluate the proteomics and transcriptomics results, the expressions of eight genes were determined in the two pig breeds by qPCR on the basis of their different influence on the formation of meat flavor and skeletal muscle development. Four of the genes encoded proteins involved in glycolysis and fat deposition (ENO1, ENO3, ATP5B, and TPI1), three encoded proteins that correlated with synthesis of muscle fiber proteins (MYLPF, ACTA1, and ACTC1) and the HSPB1 gene, whose encoded protein plays a crucial role in maintaining cellular homeostasis [54] and protects against stress [23]. The fold-changes in expression of the eight genes were compared with the transcriptomics and proteomics analysis results. As shown in Fig. 7, the qPCR results revealed that five genes were in agreement with the RNA-Seq results and six genes showed a similar trend to proteomics results. In the present research, we also observed different expression trends between the proteomics and transcriptome analyses.

Fig. 7
figure 7

qRT-PCR validation of the differentially expressed genes analyzed by RNA-seq and 2-DE. qRT-PCR was performed for eight genes that were identified as differential expressed genes between the Yorkshire and Shaziling pig breeds. The Y axis shows the relative expression levels

Correlation analysis of mRNA and protein expression

RNA-seq analyses identified 488 DEGs (Additional file 3: Table S3), of which 297 were upregulated in the Shaziling pig and 192 in the Yorkshire pig. Proteomics revealed identified 38 differentially abundant proteins, of which 27 protein spots were upregulated in the Shaziling pig and 11 in the Yorkshire pig. Similar to previous reports, the transcriptomics and proteomics data were divergent. In the present study, ENO1 and ACTC1 were overexpressed in Shaziling pigs, and ATP5B was overexpressed in Yorkshire pigs: only for these three genes were the proteomics results consistent with the transcriptomic results. In 2009, Timperio and colleagues performed a comparative analysis of proteomics and transcriptomics from the livers of Chianina and Holstein Friesian cattle. The results indicated that only three of 39 differentially abundant proteins were validated by microarray analyses [55]. Other research also confirmed that proteomics and transcriptomics data seldom overlap [56]. These differences were probably caused by alternative splicing, differential regulation of translation, and annotation errors of databases [57]. Another aspect concerning little overlap between transcriptomics and proteomics data is biological factors [58]. Although proteomics and transcriptomics data have almost no overlap, interaction pathway analyses could indicate shared biological significance [33]. Taking this into consideration, the differentially abundant proteins and DEGs that converged in the same metabolic pathways, especially regulation of skeletal muscle development, were meaningful. Some of the proteins and upregulated gene transcripts in Shaziling pigs were found to be involved in the same metabolic pathways, particularly the glycolytic pathway (ENO1, TPI1, and HSPB1). ENO1, a glycolytic enzyme, is positively correlated with meat tenderness [59]. TPI is also a glycolytic enzyme, and has been shown to correlate with meat tenderness in porcine muscles [60]. Notably, although the TPI1 result of proteomics and transcriptomics data did not match, pathway analyses of either DEG transcripts or proteins for the Shaziling samples were involved in a metabolic network. HSP proteins are related to protein folding and the oxidative stress response. In our research, HSPB1 was overexpressed in Shaziling pigs, and might be positively correlated with meat quality, which agreed with previous studies [61, 62]. Shaziling pigs have excellent meat quality like other Chinese indigenous pig breeds. The IMF content in Shaziling pigs is 3.5 %, in Jinhua pigs it is 3.38 %, and in Lantang pigs it is 2.46 %. By contrast, Yorkshire pigs and Landrace pigs have IMFs of 1.79 % [7, 63]. Increased IMF content can improve meat quality significantly, especially in terms of tenderness [64]. According to our results, genes related to tenderness have a higher level of expression in Shaziling pigs than in Yorkshire pigs, for example TPI1, HSPB1, and ENO1. Further analysis of the DEGs identified a number of novel genes and pathways (Additional file 1: Tables S1, Additional file 2: Table S2 and Additional file 3: Table S3), which have not been reported to affect meat quality previously. Further characterization of these novel genes might reveal the regulatory mechanism underlying meat quality.


The object of this study was to investigate differences in the growth and development of skeletal muscle between Shaziling and Yorkshire pigs. The combined use of proteomic and transcriptomic analyses was effective in detecting DEGs and proteins. As a result, 38 differentially abundant proteins and 488 DEGs were identified by mass spectrometry and RNA-seq analysis, respectively. Some of the proteins and unigenes are associated with lipid metabolism or glycolytic metabolism, according to previously published results. Based on the putative results of GO term enrichment and KEGG pathway analyses, we determine that many of the differential abundant proteins and DEGs are related to lipid mobilization, energy metabolism, the cytoskeleton, and signal transduction. Our study provided valuable information that could contribute to a deeper understanding of the molecular mechanisms regulating the development and formation of skeletal muscle.


2-DE, two-dimensional fluorescence difference gel electrophoresis; DEGs, differentially expressed genes; GO, gene ontology; IMF, intramuscular fat; KEGG, Kyoto Encyclopedia of Genes and Genomes; LM, longissimus dorsi muscle; MS, mass spectrometry; NGS, next-generation sequencing; qPCR, quantitative real-time PCR


  1. Rosenvold K, Andersen HJ. Factors of significance, for pork quality – a review. Meat Sci. 2003;64:219–37.

    Article  PubMed  Google Scholar 

  2. Megens HJ, Crooijmans RP, San Cristobal M, Hui X, Li N, Groenen MA. Biodiversity of pig breeds from China and Europe estimated from pooled DNA samples: differences in microsatellite variation between two areas of domestication. Gen Sel Evol. 2008;40(1):103–28.

    Google Scholar 

  3. Scherf BD. World watch list for animal diversity. 3rd ed. Rome: FAO; 2000.

    Google Scholar 

  4. Lefaucheur L. A second look into fibre typing – relation to meat quality. Meat Sci. 2010;84(2):257–70.

    Article  CAS  PubMed  Google Scholar 

  5. Xu YJ, Jin ML, Wang LJ, Zhang AD, Zuo B, Xu DQ, et al. Differential proteome analysis of porcine skeletal muscles between Meishan and Large White. J Anim Sci. 2009;87(8):2519–27.

    Article  CAS  PubMed  Google Scholar 

  6. Zhao SM, Ren LJ, Chen L, Zhang X, Cheng ML, Li WZ, et al. Differential expression of lipid metabolism related genes in porcine muscle tissue leading to different intramuscular fat deposition. Lipids. 2009;44(11):1029–37.

    Article  CAS  PubMed  Google Scholar 

  7. Wu T, Zhang Z, Yuan Z, Lo LJ, Chen J, Wang Y, et al. Distinctive genes determine different intramuscular fat and muscle fiber ratios of the longissimus dorsi muscles in Jinhua and landrace pigs. PLoS One. 2013;8(1):e53181.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Lu P, Li DF, Yin JD, Zhang LY, Wang ZY. Flavour differences of cooked longissimus muscle from Chinese indigenous pig breeds and hybrid pig breed (Duroc× Landrace× Large White). Food Chem. 2008;107(4):1529–37.

    Article  CAS  Google Scholar 

  9. Liu X, Xiong X, Yang J, Zhou L, Yang B, Ai H, et al. Genome-wide association analyses for meat quality traits in Chinese Erhualian pigs and a Western Duroc×(Landrace× Yorkshire) commercial population. Genet Sel Evol. 2015;47(1):44.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Guo J, Shan T, Wu T, Zhu LN, Ren Y, An S, et al. Comparisons of different muscle metabolic enzymes and muscle fiber types in Jinhua and Landrace pigs. J Anim Sci. 2011;89(1):185–91.

    Article  CAS  PubMed  Google Scholar 

  11. Kim NK, Joh JH, Park HR, Kim OH, Park BY, Lee CS. Differential expression profiling of the proteomes and their mRNAs in porcine white and red skeletal muscles. Proteomics. 2004;4(11):3422–8.

    Article  CAS  PubMed  Google Scholar 

  12. Guo Y, Li J, Shang J, Jin Y. Development of muscle-related genes and their effects on meat quality. Energy Procedia. 2012;16:229–33.

    Article  Google Scholar 

  13. Picard B, Lefaucheur L, Berri C, Duclos MJ. Muscle fibre ontogenesis in farm animal species. Reprod Nutr Dev. 2002;42(5):415–31.

    Article  PubMed  Google Scholar 

  14. Zhao X, Mo D, Li A, Gong W, Xiao S, Zhang Y, et al. Comparative analyses by sequencing of transcriptomes during skeletal muscle development between pig breeds differing in muscle growth rate and fatness. PLoS One. 2011;6(5):e19774.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Karlsson AH, Klont RE, Fernandez X. Skeletal muscle fibres as factors for pork quality. Livest Prod Sci. 1999;60(2):255–69.

    Article  Google Scholar 

  16. Wimmers K, Murani E, Ngu NT, Schellander K, Ponsuksili S. Structural and functional genomics to elucidate the genetic background of microstructural and biophysical muscle properties in the pig. J Anim Breed Genet. 2007;124(1):27–34.

    Article  PubMed  Google Scholar 

  17. Quiroz-Rothe E, Rivero JL. Coordinated expression of myosin heavy chains, metabolic enzymes, and morphological features of porcine skeletal muscle fiber types. Micros Res Techniq. 2004;65(1–2):43–61.

    Article  CAS  Google Scholar 

  18. da Costa N, McGillivray C, Chang KC. Postnatal myosin heavy chain isoforms in prenatal porcine skeletal muscles: insights into temporal regulation. Anat Rec Part A. 2003;273(2):731–40.

    Article  Google Scholar 

  19. Lefaucheur L, Ecolan P, Plantard L, Gueguen N. New insights into muscle fiber types in the pig. J Histochem Cytochem. 2002;50(5):719–30.

    Article  CAS  PubMed  Google Scholar 

  20. Liu J, Damon M, Guitton N, Guisle I, Ecolan P, Vincent A, et al. Differentially-expressed genes in pig Longissimus muscles with contrasting levels of fat, as identified by combined transcriptomics, reverse transcription PCR, and proteomic analyses. J Agri Food Chem. 2009;57(9):3808–17.

    Article  CAS  Google Scholar 

  21. Xu YJ, Qian H, Feng XT, Xiong YZ, Lei MG, Ren ZQ, et al. Differential proteome and transcriptome analysis of porcine skeletal muscle during debelopmemt. J Proteomics. 2012;75(7):2093–108.

    Article  CAS  PubMed  Google Scholar 

  22. Li A, Mo D, Zhao X, Jiang W, Cong P, He Z. Comparison of the longissimus muscle proteome between obese and lean pigs at 180 days. Mamm Genome. 2013;24(1–2):72–9.

    Article  CAS  PubMed  Google Scholar 

  23. Görg A, Weiss W, Dunn MJ. Current two dimensional electrophoresis technology for proteomics. Proteomics. 2004;4(12):3665–85.

    Article  PubMed  Google Scholar 

  24. Timms JF, Cramer R. Difference gel electrophoresis. Proteomics. 2008;8(23–24):4886–97.

    Article  CAS  PubMed  Google Scholar 

  25. Lin CS, Hsu CW. Differentially transcribed genes in skeletal muscle of Duroc and Taoyuan pigs. J Anim Sci. 2005;83(9):2075–86.

    CAS  PubMed  Google Scholar 

  26. Hollung K, Grove H, Færgestad EM, Sidhu MS, Berg P. Comparison of muscle proteome profiles in pure breeds of Norwegian Landrace and Duroc at three different ages. Meat Sci. 2009;81(3):487–92.

    Article  CAS  PubMed  Google Scholar 

  27. Kim NK, Lim JH, Song MJ, Kim QH, Park BY, Kim MJ, et al. Developmental proteomic profiling of porcine skeletal muscle during postnatal development. Asian Austral J Anim. 2007;20(10):1612–7.

    Article  CAS  Google Scholar 

  28. Marguerat S, Bähler J. RNA-seq: from technology to biology. Cell Mol Life Sci. 2010;67(4):569–79.

    Article  CAS  PubMed  Google Scholar 

  29. Ramayo-Caldas Y, Mach N, Esteve-Codina A, Corominas J, Castello A, Ballester M, et al. Liver transcriptome profile in pigs with extreme phenotypes of intramuscular fatty acid composition. BMC Genomics. 2012;13(1):547.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Chen C, Ai H, Ren J, Li W, Li P, Qiao R, et al. A global view of porcine transcriptome in three tissues from a full-sib pair with extreme phenotypes in growth and fat deposition by pairedend RNA sequencing. BMC Genomics. 2011;12(1):448.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Tian Q, Stepaniants SB, Mao M, Weng L, Feetham MC, Doyle MJ. Integrated genomic and proteomic analyses of gene expression in mammalian cells. Mol Cell Proteomics. 2004;3(10):960–9.

    Article  CAS  PubMed  Google Scholar 

  32. Murgiano L, D’Alessandro A, Egidi MG, Crisa A, Prosperini G, Timpero AM, et al. Proteomics and transcriptomics investigation on longissimus muscles in Large White and Casertana pig breeds. J Proteome Res. 2010;9(12):6450–66.

    Article  CAS  PubMed  Google Scholar 

  33. Kim NK, Park HR, Lee HC, Yoon D, Son ES, Kim YS, et al. Comparative studies of skeletal muscle proteome and transcriptome profilings between pig breeds. Mamm Genome. 2010;21(5–6):307–19.

    Article  CAS  PubMed  Google Scholar 

  34. García-Alcalde F, Okonechnikov K, Carbonell J, Cruz LM, Götz S, Tarazona S, et al. Qualimap: evaluating next-generation sequencing alignment data. Bioinformatics. 2012;28(20):2678–9.

    Article  PubMed  Google Scholar 

  35. Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136–8.

    Article  PubMed  Google Scholar 

  36. Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3(6):1101–8.

    Article  CAS  PubMed  Google Scholar 

  37. Orosz F, Oláh J, Ovádi J. Triosephosphate isomerase deficiency: new insights into an enigmatic disease. Biochim Biophys Acta. 2009;1792(12):1168–74.

    Article  CAS  PubMed  Google Scholar 

  38. Laville E, Sayd T, Terlouw C, Chambon C, Damon M, Larzul C, et al. Comparison of sarcoplasmic proteomes between two groups of pig muscles selected for shear force of cooked meat. J Agric Food Chem. 2007;55(14):5834–41.

    Article  CAS  PubMed  Google Scholar 

  39. Farrar WW, Deal Jr WC. Purification and properties of pig liver and muscle enolases. J Protein Chem. 1995;14(6):487–97.

    Article  CAS  PubMed  Google Scholar 

  40. Kwasiborski A, Sayd T, Chambon C, Sante’-Lhoutellier V, Rocha D, Terlouw C. Pig Longissimus lumborum proteome: part II: relationships between protein content and meat quality. Meat Sci. 2008;80(4):982–96.

    Article  CAS  PubMed  Google Scholar 

  41. klont RE, Brocks l, Eikelenboom G. Muscle fibre type and meat quality. Meat Sci. 1998;49:s219–29.

    Article  Google Scholar 

  42. Ryu YC, Kim BC. Comparison of histochemical characteristics in various pork groups categorized by postmortem metabolic rate and pork quality. J Anim Sci. 2006;84(4):894–901.

    CAS  PubMed  Google Scholar 

  43. Helferich W, Jump DB, Anderson DB, Skjaerlund DM, Merkel RA, Bergen WG. Skeletal muscle alpha-actin synthesis is increased pretranslationally in pigs fed the phenethanolamine ractopamine. Endocrinology. 1990;126(6):3096–100.

    Article  CAS  PubMed  Google Scholar 

  44. Guillemin N, Bonnet M, Jurie C, Picard B. Functional analysis of beef tenderness. J Proteomics. 2011;75(2):352–65.

    Article  CAS  PubMed  Google Scholar 

  45. D’Alessandro A, Marrocco C, Zolla V, D’Andrea M, Zolla L. Meat quality of the longissimus lumborum muscle of Casertana and Large White pigs: metabolomics and proteomics intertwined. J Proteomics. 2011;75(2):610–27.

    Article  PubMed  Google Scholar 

  46. Huang X, Madan A. CAP3: a DNA sequence assembly program. Genome Res. 1999;9(9):868–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Wickramasinghe S, Cánovas A, Rincón G, Medrano JF. RNA-sequencing: a tool to explore new frontiers in animal genetics. Livest Sci. 2014;166:206–16.

    Article  Google Scholar 

  48. Qian X, Ba Y, Zhuang Q, Zhong G. RNA-Seq technology and its application in fish transcriptomics. OMICS. 2014;18(2):98–110.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Gui D, Jia K, Xia J, Yang L, Chen J, Wu Y, et al. De novo assembly of the indo-pacific humpback dolphin leucocyte transcriptome to identify putative genes involved in the aquatic adaptation and immune response. PLoS One. 2013;8(8):e72417.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Ruan R, Guo AH, Hao YJ, Zheng JS, De WD. Novo assembly and characterization of narrow-ridged finless porpoise renal transcriptome and identification of candidate genes involved in osmoregulation. Int J Mol Sci. 2015;16(1):2220–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.

    Article  CAS  PubMed  Google Scholar 

  52. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2011;40(Database issue):D109–14.

    PubMed  PubMed Central  Google Scholar 

  53. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.

    Article  CAS  PubMed  Google Scholar 

  54. Liu Y, Steinacker JM. Changes in skeletal muscle heat shock proteins: pathological significance. Front Biosci. 2001;6:D12–25.

    CAS  PubMed  Google Scholar 

  55. Timperio AM, D’Alessandro A, Pariset L, D’Amici GM, Valentini A, Zolla L. Comparative proteomics and transcriptomics analyses of livers from two different Bos taurus breeds: “Chianina and Holstein Friesian”. J Proteomics. 2009;73(7):309–22.

    Article  CAS  PubMed  Google Scholar 

  56. Nie L, Wu G, Culley DE, Schotten JCM, Zhang W. Integrative analysis of transcriptomic and proteomic data: challenges, solutions and applications. Crit Rev Biotechnol. 2007;27(2):63–75.

    Article  CAS  PubMed  Google Scholar 

  57. Hornshøj H, Bendixen E, Conley LN, Andersen PK, Hedegaard J, Panitz F. Transcriptomic and proteomic profiling of two porcine tissues using high-throughput technologies. BMC Genomics. 2009;10(1):30.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Samulin J, Berget I, Lien S, Sundvold H. Differential gene expression of fatty acid binding proteins during porcine adipogenesis. Comp Biochem Physiol B: Biochem Mol Biol. 2008;151(2):147–52.

    Article  Google Scholar 

  59. D’Alessandroa A, Rinalduccia S, Marroccoa C, Zollaa V, Napolitanob F, Zollaa L. Love me tender: an Omics window on the bovine meat tenderness network. J Proteomics. 2012;75(2):4360–80.

    Article  Google Scholar 

  60. Hwang IH, Park BY, Kim JH, Cho SH, Lee JM. Assessment of postmortem proteolysis by gel-based proteome analysis and its relationship to meat quality traits in pig longissimus. Meat Sci. 2005;69(1):79–91.

    Article  CAS  PubMed  Google Scholar 

  61. Bouley J, Meunier B, Chambon C, De Smet S, Hocquette JF, Picard B. Proteomic analysis of bovine skeletal muscle hypertrophy. Proteomics. 2005;2:490–500.

    Article  Google Scholar 

  62. Morzel M, Terlouw C, Chambon C, Micol D, Picard B. Muscle proteome and meat eating qualities of Longissimus thoracis of “Blonde d’ Aquitaine” young bulls: a central role of HSP27 isoforms. Meat Sci. 2008;78:297–304.

    Article  CAS  PubMed  Google Scholar 

  63. Yu K, Shu G, Yuan F, Zhu X, Gao P, Wang S, et al. Fatty acid and transcriptome profiling of longissimus dorsi muscles between pig breeds differing in meat quality. Int J Biol Sci. 2013;9(1):108–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Van Laack R, Stevens S, Stalder K. The influence of ultimate pH and intramuscular fat content on pork tenderness and tenderization. J Anim Sci. 2001;79:392–7.

    PubMed  Google Scholar 

Download references


We would like to thank LY for exceptional assistance in the qPCR and DX in collection of samples of pig tissue.


This study was supported by the National High Technology Research and Development Program of China (2011AA100304), the Provincial Natural Science Foundation of Hunan (13JJ1021) and Science and technology project of hunan province (2014NK4135).

Availability of data and materials

The sequencing data from this study has been archived at the NCBI Sequence Read Archive under GEO database [GSE:70673].

Authors’ contributions

HY performed analysis of the transcriptomics and proteomics and drafted the manuscript. XX and JJ participated in the biological experiments and data analysis. HM managed the whole project. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

The experiment and all its procedures were reviewed and approved by the Institutional Animal Care and Use Committee (IACUC) at Hunan Agricultural University.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Hai-ming Ma.

Additional files

Additional file 1: Table S1.

KEGG Pathway information of unigenes. (XLS 705 kb)

Additional file 2: Table S2.

GO enrichment of DEG. (XLS 77 kb)

Additional file 3: Table S3.

Summary of the differently expressed genes. (XLS 58 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yang, H., Xu, Xl., Ma, Hm. et al. Integrative analysis of transcriptomics and proteomics of skeletal muscles of the Chinese indigenous Shaziling pig compared with the Yorkshire breed. BMC Genet 17, 80 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: