- Research article
- Open Access
iTRAQ-based quantitative proteomic analysis reveals the distinct early embryo myofiber type characteristics involved in landrace and miniature pig
BMC Genomics volume 17, Article number: 137 (2016)
Pig (Sus scrofa) is a major source of dietary proteins for human consumption and is becoming a valuable model in agricultural and biomedical research. The recently developed isobaric tag for relative and absolute quantitation (iTRAQ) method allows sensitive and accurate protein quantification. Here, we performed the first iTRAQ-based quantitative proteomic analyses of Landrace (LR) and Wuzhishan (WZS) pig longissimus dorsi muscle tissues during early embryonic development.
The iTRAQ-based early embryonic longissimus dorsi muscle study of LR and WZS ranging from 21 to 42 days post coitus (dpc) identified a total of 4431 proteins from17,214 peptides, which were matched with 36,4025 spectra at a false discovery rate of 5 %. In both WZS and LR, the largest amount of differentially expressed proteins (DEPs) were found between 28 and 35 dpc. 252 breed-DEPs were selected by GO analysis, including 8 myofibrillar proteins. Only MYHCI/IIA mRNA were detected due to early embryonic stages, and significantly higher expression of them were found in WZS during these 4 stages. MYHCI was first found in WZS at 28 dpc and expressed in both breeds at later stages, while MYHCII protein was not detected until 35 dpc in both breeds. Thus, 33 myogenic breed-DEPs selected from last two stages were analyzed by STRING, which showed that some myofibrillar proteins (MYH1, TPM4, MYH10, etc.) and functional proteins (CSRP2, CASQ2, OTC, etc.), together with candidate myogenic proteins (H3F3A, HDGFRP2, etc.), probably participate in the regulatory network of myofiber formation.
Our iTRAQ-based early embryonic longissimus dorsi muscle study of LR and WZS provides new data on the in vivo muscle development distinctions during early embryonic development, which contributes to the improved understanding in the regulation mechanism of early myogenesis in agricultural animals.
Pig is a major dietary protein source for human consumption and valuable model in agricultural and biomedical research [1–3], however, few studies concerning pig proteome were reported recently. With high-throughput technologies for protein profiling developing, great deal of molecular data will contribute to the understanding of pig early embryonic myogenesis. Pig myogenesis is clearly a biphasic phenomenon with primary myofibers forming from 35 to 55 dpc, followed by secondary myofibers forming around each primary myofiber between 50 and 90 dpc . The muscle growth is predominantly determined during prenatal skeletal muscle development [5–7].
Skeletal muscle of postnatal pig consists of four major fiber types with type-specific expression of corresponding myosin heavy chain (MYHC) isoforms and specific metabolic pattern, namely slow-oxidative type I, fast-oxidative type IIa, fast oxidative-glycolytic type IIx and fast-glycolytic IIb [8, 9]. Slow oxidative type I myofibers are hypothesized to be beneficial for pork quality, while PSE (pale, soft, exudative) pork is often associated with higher percentage of MYHC IIb fibers [10–12]. Although myofiber types are influenced by neuromuscular activity, mechanical loading, and hypothyroidism, myofiber composition is of high heritability in pigs [13, 14].
High-throughput and shotgun-based proteomic techniques are increasingly utilized to improve proteome coverage and produce protein catalogues . iTRAQ-mediated shotgun proteomics allows the relative quantification of peptides from eight samples simultaneously, and has proven to be a successful tool in protein biomarker discovery in biomedical field [16, 17]. Recent studies have reported 542 proteins in porcine muscle concerning gender and diet dephytinization by iTRAQ . In other livestock species, the analysis of bovine skeletal muscle by iTRAQ identified changes in protein abundance between tender and tough meat from bovine longissimus thoracic muscle .
With the goal to find out whether early embryo proteomic distinctions exist in pig breeds differing in meat quality and gain further insights into the mechanism. 21, 28, 35, 42 dpc embryos/fetuses of two pig breeds were employed: LR, an improved large pig breed, is characterized by high lean meat percentage, fast growing muscle and high body weight [7, 20, 21] and WZS pig, originates from Hainan Island, China, is known for its superior meat quality, and is considered useful for medical research due to its small size (adult weights < 25 kg) . The longissimus dorsi muscle proteome was analyzed by iTRAQ. Myofiber type characteristics were studied by real-time quantitative PCR and western blot (WB). Through these assays, we wanted to uncover the early embryo proteomic distinctions in two pig breeds and elucidate the mechanism underlying muscle development, which may contribute to the understanding of early embryonic myofiber formation in agricultural animals.
DEPs of adjacent stages in WZS and LR
In total, 4431 proteins were identified from 17,214 peptides, which were matched with 36,4025 spectra at a false discovery rate of 5 % (Fig. 1b, Additional file 1: Figure S1.).
The largest amount of up or down-regulated proteins within WZS (168 proteins up and 388 proteins down-regulated) and LR (185 proteins up and 278 proteins down-regulated) were found in the comparison between 28 and 35 dpc (Fig. 1c), and morphologic observations also showed that myofibers didn’t appear until 35 dpc in both breeds (Additional file 2: Figure S2). Taken together, we deduced that 28–35 dpc is important for embryonic muscle development for pig breeds, during which a number of proteins participated in and the myofibers began to form. In the comparisons between LR and WZS at same stages, the largest number of up-regulated DEPs (133 up-regulated DEPs) were found at 28 dpc and the largest number of down-regulated DEPs (124 down-regulated DEPs) were at 35 dpc (Fig. 1d), which indicated that early development of WZS and LR embryos were distinctly regulated by more DEPs during this period than the other two stages. Thus we speculated that 28–35 dpc is also important for embryonic development between breeds. Significantly more DEPs were found in WZS at adjacent stages studied (Fig. 1d, e). Over 9 myogenic terms (muscle cell differentiation, skeletal muscle myosin thick filament assembly, etc.) for biological process (BP) were found significantly enriched in WZS and only two (muscle system process and muscle contraction) were found significantly enriched in LR at adjacent stages (Additional file 3: Table S1), suggesting that distinct skeletal muscle development pattern has arisen at early embryonic stages between WZS and LR.
DEPs between WZS and LR at same stages
Protein quantification based on the relative amounts of the different iTRAQ labels was obtained for all the identified proteins with two or more unique peptides. The datasets for differences in protein levels between WZS and LR at same stages were analyzed then. The largest amount of DEPs were found in 28 dpc (252 DEPs) (Fig. 2a, Additional file 4: Figure S3, Additional file 5: Table S2). While at 21 and 35 dpc, smaller but still great amount of DEPs were present (219 and 215 DEPs, respectively). The smallest amount of DEPs were found at 42 dpc (175 DEPs).
After analyzing common DEPs between WZS and LR at same stages, 16 significant DEPs were selected with 4 trends based on ratio of DEPs’ expression level (WZS/LR) and GO terms for BP were also analyzed in each trend (Fig. 2b, Additional file 6: Figure S4, Additional file 3: Table S1, Additional file 5: Table S2): 1) proteins down-regulated at 21–28 dpc and up-regulated at 35–42 dpc (MDH1, PYCR2, PF4); involved GO terms for BP: energy production/conversion, amino acid transport and metabolism, 2) proteins up-regulated at 21–35 dpc and down regulated at 35–42 dpc (CC2D1B, DEK, HDGFRP2, etc.); involved GO terms for BP: signal transduction and DNA-dependent transcription, 3) proteins down-regulated at 21–42 dpc (MARCKSL1, BANF1, HSD17B14, etc.); involved GO terms for BP: positive regulation of cell proliferation and DNA integration, 4) proteins up-regulated at 21–42 dpc (FDPS, NDUFV3, CFB, etc.); involved GO terms for BP: coenzyme transport and metabolism and metabolic process (Additional file 3: Table S1).
Myogenic DEPs between WZS and LR at same stages
In term of COG database, the function classification of identified protein sequences were listed (Additional file 7: Figure S5, Additional file 8: Table S3, Additional file 9: Table S4). Top myogenic GO terms for BP with the minimal p values were chosen between WZS and LR at same stages, such as muscle system process, muscle cell differentiation and muscle fiber development, etc. The highest enrichment of myogenic BP were found at 42 dpc (Fig. 3a), six of which with p value less than 0.05 (Additional file 10: Table S5). While at 35 dpc, lower enrichment of myogenic BP was present and only one term with p value less than 0.05. At 21 and 28 dpc, no terms with p value less than 0.05 were found.
From the proteins involved in the top 5 myogenic GO terms for BP, 17 myogenic DEPs between WZS and LR were chosen (1.5-fold change) (Additional file 11: Table S6). Taken together, most myogenic DEPs (WZS/LR) were down-regulated at 21 or 28 dpc and up-regulated at 35 dpc, among which 8 myofibrillar proteins showed distinct expression pattern (Fig. 3b). Only 3–5 myogenic DEPs were found at 21 and 28 dpc, while at stage 35 and 42 dpc, 9–12 myogenic DEPs were present (Additional file 11: Table S6, Additional file 12: Table S7).
Myofiber type distinction at early embryonic stages
There are 8 myofibrillar proteins showing distinct expression between WZS and LR. In order to validate whether myofiber type distinction existed, 4 types of MYHC isoforms were examined by qPCR and only MYHCI/IIA were detected in WZS and LR. Both of them were found with significantly higher expression in WZS covering 4 stages.
To further validate the distinction of myofiber development, MYHCI/II proteins were detected by WB. The expression of MYHCII protein was not detected until 35 dpc in both breeds, while MYHCI was first found in WZS at 28 dpc and expressed in both breeds at later stages. It’s interesting that both MYHCI and MYHCII proteins were significantly higher expressed in WZS at 35 dpc, while until 42 dpc, MYHCI was higher expressed in LR (Fig. 4b). Morphologic observations showed that myofibers didn’t appear until 35 dpc in both breeds and myofibers in WZS at 35–42 dpc tended to be larger in size and more in number than those in LR (Additional file 2: Figure S2).
Taken together, at early embryonic stages, myofiber type distinction has been existed between two breeds and WZS showed stronger trend of myofiber differentiation during 35–42 dpc.
Protein interaction networks of significantly DEPs
Protein interactions were predicted in the website of STRING and the interaction network was illustrated by Cytoscape software. There are two cycles in the figure (Fig. 5). Myofibrillar proteins are in the inner cycle; other proteins that were predicted to have interactions with inner cycle proteins are in the outer cycle. As predicted, several strong interactions were found between outer (Ornithine Carbamoyltransferase, Angiotensinogen, Sorcin, Calsequestrin 2) and inner cycle myofibrillar proteins (Troponin C, Myosin regulatory light chain 7, Myosin 4, Myosin 1, Tropomyosin alpha-4) (Fig. 5). Since most of the outer cycle proteins were up or down regulated at both 35 and 42 dpc, and most of the inner cycle proteins were up-regulated at 35 dpc and down-regulated at 42 dpc, it can be reasonably speculated that the outer cycle proteins regulate the expression profile of inner cycle proteins and the distinct expression of outer cycle proteins in WZS and LR may result in distinct characteristics of myofiber formation.
In this study, a comprehensive evaluation of proteomic profile in 21, 28, 35, and 42 dpc of LR and WZS pig were performed, providing new data about distinctions in vivo during early embryonic muscle development. The results showed higher and earlier expression of type I/II MYHC proteins in WZS, which indicate distinct regulation network of early embryonic myofiber development between pig breeds differing in meat quality.
Distinct muscle development pattern in WZS and LR
In this study, porcine embryonic longissimus dorsi muscle proteomic analysis covered the process of primary myofiber formation. In comparison of WZS and LR at adjacent stages, much more DEPs and significantly enriched myogenesis BP terms were found in WZS. In comparison between WZS and LR at same stages, although large amount of DEPs were found in 21 and 28 dpc, no significant myogenic GO terms were found. While at 35 and 42 dpc, though the number of DEPs dropped significantly, higher enrichment of myogenic GO terms were found. These findings indicated that 21–28 dpc may be important stages for the preparation of early muscle development and the DEPs involved in these stages probably participate in the distinct muscle development process between WZS and LR since 35 dpc.
Distinct myofiber characteristics between WZS and LR
To further investigate the myogenic DEPs between WZS and LR, we analyzed the proteins involved in top myogenic GO terms for BP. Many myofibrillar proteins (myosin, troponin, etc.) are known to exist as multiple isoforms . It has been demonstrated that some myofibrillar isoforms are differentially expressed in various muscles and fiber types. Myosin makes up about one-third of the total muscle proteins and is known to produce multiple isoforms. Different myosins involved in differentiate progress of distinct muscle fiber types: type I (slow-twitch, red muscle, oxidative) (MYH7), type IIa (fast-twitch, red muscle, oxidative) (MYH2), and type IIb/IIx (fast-twitch, white muscle, glycolytic) (MYH1) [8, 24].
In this study, multiple differentially expressed myosin proteins were found, among which MYH11, MYL9, MYH6, TNNI1, MYL7, TPM3, TNNI3 were up-regulated in WZS at least 1 stage we studied. In addition, mRNA and protein level of MYHCI/II isoforms also showed higher expression in WZS at all 4 stages. Longissimus dorsi muscle is a fast-twitch glycolytic muscle involved in voluntary movements of the back. It was reported that faster growing commercial pig breeds with more muscle contain high proportions of type IIb/IIx glycolytic fibers (75–86 %) . No studies concerning WZS myofiber composition were reported, while previous studies indicated that skeletal muscles of Meishan pig, another Chinese indigenous pig breed, weighed only 36.1 kg in average before mating , comprises of higher proportion of slow fibers at slaughter, which was hypothesized to be the major attribute to the superior meat quality . As we mentioned before, myofiber composition is of high heritability in pigs. Previous study concerning later embryonic stages to adult pig shows that myofiber types at different stages and different part of the longissimus dorsi muscle have distinct composition of myofiber types , while this study focused on the early embryo stages including primary myofiber formation. Thus, we speculated that the myofiber type distinction began to show up between LR and WZS in early embryonic stages and changing through later developmental stages, which determined the distinct adult myofiber composition between LR and WZS.
DEPs involved in the regulatory network of early muscle development
Several DEPs showed significant changes during these stages. CSRP2 associates with meat quality traits with CSRP1 and CSRP3 in pigs . This protein was down-regulated at 28 dpc and up-regulated at 42 dpc in WZS. H3F3A (Histone H3.3) dropped 74-96 % from 28–35 dpc but increased 2.5 folds at 42 dpc in WZS. Hepatoma-derived growth factor-related protein 2 (HDGFRP2) was up-regulated 2.4 fold at 35 dpc and dropped 60 % at 42 dpc in WZS, and it includes an N-terminal PWWP domain that binds to methyl-lysine-containing histones, with specific binding of this protein to tri-methylated lysines 36 and 79 of histone H3, and di- and tri-methylated lysine 20 of histone H4 . It is reported that myoblast determination protein (MyoD) determines cell fate and facilitates differentiation-dependent gene expression through chromodomain helicase DNA binding protein 2 (CHD2) -dependent deposition of H3.3 at myogenic loci prior to differentiation . Recent study found that selective H3.3 incorporation is essential for establishing specific modifications in myogenic genes, suggesting that the incorporation of specific histone H3 variants determines the lineage potential of myogenesis . Thus the distinct expression of H3F3A and HDGFRP were very likely to regulate the distinct early muscle development between WZS and LR.
Several myofibrillar proteins and functional proteins involving in strong interaction networks were found by STRING analysis (Fig. 5). For example, SRI plays a role in the translation of mechanical cues applied to myofibers into paracrine signals, which in turn modulate satellite cell functions and support muscle growth . It was down-regulated at 28 dpc and up-regulated at 35 dpc in WZS. CASQ2, a calcium binding protein that stores calcium for muscle function , was up-regulated since 21 dpc in WZS. The gene of OTC was mapped closely to Duchenne muscular dystrophy, which may play a role in myogenesis . Taken together, most myogenic DEPs were down-regulated at 21 or 28 dpc and up-regulated at 35 or 42 dpc in WZS, therefore, it’s reasonable to believe that these functional proteins participated in the regulatory network of distinct myofiber development process. In adult mice, a selective accumulation of MyoD transcript in the fast-twitch muscles and myogenin transcripts in the slow-twitch skeletal muscles suggests that fibre characteristics may be affected by the ratio of the different muscle regulatory factors [35, 36], but few studies about regulation mechanism of myofiber formation in agricultural animals are reported.
Based on the above findings, further studies could focus on when and how the myogenic related factors differentially express between the two pig breeds in embryos. As in our study, the expression of H3F3A in LR was significantly higher than that in WZS (about 20 folds) at 21 dpc. Thus we speculate that this histone modification may play a key role in determining muscle development in embryos. As recent study shows that mice embryonic early skeletal muscle lineage potential is established through a selective incorporation of specific H3 variants that governs the balance of histone modifications .
In summary, this iTRAQ-based early embryonic muscle study of WZS and LR showed that the distinction of myofiber characteristics began in the early embryonic stages, which probably result in distinct myofiber type composition after birth. Our study contributes to the understanding of early embryonic myofiber formation and differentiation.
All experiments were carried out according to China Council on Animal Care and the protocols used were approved by the Animal Care and Use Committee of Guangdong Province, China. Approval ID or permit numbers are SCXK (Guangdong) 2011–0029 and SYXK (Guangdong) 2011–0112 (Additional file 13).
Twenty-four purebred sows with the same genetic background of LR and WZS were artificially inseminated with semen from the same purebred boars. For each breed in prenatal stages, three sows per time point were slaughtered at 21, 28, 35 and 42 dpc after insemination, and embryos/fetuses were collected as described before . Whole embryos of 21 dpc (Additional file 14: Figure S6) and longissimus dorsi muscle tissues from fetuses of left stages were used. These samples were snap-frozen in liquid nitrogen and stored until further use.
Protein isolation and iTRAQ labeling
For protein analyzing, tissues were processed as described before . As the amount of protein in one embryo was small (especially at 18 & 21 dpc), we used a single pooling strategy in our study as previous similar studies reported [39–41]. Previous study also pointed that the pooling strategy is appropriate in experiments with low sample yield, and could reduce the biological variance, increasing the power to detect changes . Three embryos/fetuses of different sows were mixed as a pool for each breeds at every stages. 100 μg protein for each sample was digested with Trypsin Gold (Promega, USA) at a protein: trypsin ratio of 20:1 at 37 °C for 4 h. Then peptides were dried via vacuum centrifugation. The samples were labeled with 113 (LR 21 dpc), 114 (LR 28 dpc), 115 (LR 35 dpc), 116 (LR 42 dpc), 117 (WZS 21 dpc), 118 (WZS 28 dpc), 119 (WZS 35 dpc) and 121 (WZS 42 dpc) iTRAQ reagent, respectively, by incubation at room temperature for 2 h. To reduce sample complexity, strong cation exchange (SCX) chromatography, using an off-line SCX cartridge and buffers supplied by ABI, was carried out for cell lysate-derived peptides.
LC-ESI-MS/MS analysis by LT Q-Orbitrap HCD
The peptides were subjected to nanoelectrosprayionization followed by tandem mass spectrometry (MS/MS) in a LTQ Orbitrap Velos (Thermo Fisher, USA) coupled online to the HPLC (Shimadzu, Kyoto, Japan) and processed as described before . For MS scans, the m/z scan range was 350 to 2000 Da. The experiment was repeated three times, and the results were categorized as 113–121 groups.
Database search, protein identification and quantification
The MS/MS data were searched against the Ensembl: Sus_scrofa.Sscrofa1 0.2.69.pep.all.fa (25152 sequences) database for peptide identification and quantification using the Mascot 2.3.02.
For protein identification, a mass tolerance of 6 Da (ppm) was permitted for intact peptide masses and 0.5 Da for fragmented ions, with allowance for one missed cleavages in the trypsin digests. Gln- > pyro-Glu (N-term Q), Oxidation (M), Deamidated (NQ) as the potential variable modifications, and Carbamidomethyl (C), iTRAQ8plex (N-term), iTRAQ8plex (K) as fixed modifications. The charge states of peptides were set to 2+, 3+, 4+. Specifically, an automatic decoy database search was performed in Mascot by choosing the decoy checkbox in which a random sequence of database is generated and tested for raw spectra as well as the real database. To reduce the probability of false peptide identification, only peptides at the 95 % confidence interval by a Mascot probability analysis greater than “identity” were counted as identified. After this filter, the FDR was very low, less than 1.3 %. And each confident protein identification involve at least one unique peptide.
For protein quantitation, it was required that a protein contains at least two unique spectra. The quantitative protein ratios were weighted and normalized by the median ratio in Mascot. We set a 1.5-fold change as the threshold and p-value (t-test) < 0.05 to identify significant changes.
RNA Extraction, and real-time quantitative PCR
To validate the iTRAQ data that myofibrillar proteins distincted from two breeds, we use real-time quantitative PCR to detect the mRNA level of type I/II MYHC isoforms. RNA was extracted from the same samples of iTRAQ analysis using NucleoSpin RNA II Extract kit (Macherey-Nagel, USA), and the quality assessed with a NanoDrop ND-1000 (Thermo Scientific, USA). cDNA was synthesized using the Transcriptor First Strand cDNA Synthesis Kit (Roche, USA).
Real-time quantitative PCR reactions were carried out in triplicate using the LightCycler 480 System (Roche, USA). The expression of each gene was normalized to that of GAPDH transcripts. Results are given as mean ± SD. The single (*), double (**) and triple (***) asterisks represent the p values, p < 0.05, p < 0.005 and p < 0.001, respectively for Student’s unpaired t tests. The following oligonucleotides were used:
MYHCI, (forward) CACTTGCTAAGAGGGACCTCTGAGTTCA; (reverse) ATCCAGGCTGCGTAACGCTCTTTGAGGTTGTAMYHCIIA, (forward) AGCCTCTTTCTTCTCCCAGGGACATTC; (reverse) ATCCAGGCTGCGTAACGCTCTTTGAGGTTGTAGAPDH, (forward) TACATGGTCTACATGTTCCAGTATG; (reverse) CAGGAGGCATTGCTGACAATCTTG
For analyzing type I/II MYHC protein isoforms, proteins were extracted and run on a 12 % SDS-PAGE gel and transferred onto a PVDF membrane. Anti-fast skeletal myosin antibody (Abcam, China) and anti-Slow Skeletal Myosin Heavy chain antibody (Abcam, China) were used for MYHC and anti-GAPDH-71.1 antibody (Sigma, USA) for GAPDH. The anti-mouse and anti-rabbit secondary antibody (CST, USA) was used for each experiment. Blots were visualized using a commercial enhanced chemiluminescene (ECL, USA) detection Kit (Thermo Scientific, USA).
Bioinformatics and statistical analysis
The cellular component (CC), molecular function (MF) and biological process (BP) of the proteins identified by iTRAQ were annotated by the Blast2GO software (https://www.blast2go.com/); GO functional classification and enrichment analysis were also performed to identify GO terms that were significantly enriched in DEPs using DAVID analysis (http://david.abcc.ncifcrf.gov/). The myogenic DEPs and predicted protein interaction network was illustrated by STRING (http://string-db.org/) and Cytoscape 3.1.0 (http://www.cytoscape.org/). Relative expression of proteins was based on the ratio of the reporter ions of peptides using Mascot software. T-test was used to analyze the p-value. Moreover, the protein with the 1.5-fold change as well as the p value less than 0.05 was designated the DEPs. Statistical analyses were conducted with SPSS 10.0.
Availability of supporting data
The datasets supporting the results of this article are included within the article (and its additional file(s)). The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD002481.
Cysteine and glycine-rich protein 2
Differentially expressed proteins
Days post coitus
High mobility group protein B1
Strong cation exchange
Schook L, Beattie C, Beever J, Donovan S, Jamison R, Zuckermann F, Niemi S, Rothschild M, Rutherford M, Smith D. Swine in biomedical research: creating the building blocks of animal models. Anim Biotechnol. 2005;16(2):183–90.
Patterson JK, Lei XG, Miller DD. The pig as an experimental model for elucidating the mechanisms governing dietary influence on mineral absorption. Exp Biol Med. 2008;233(6):651–64.
Lunney JK. Advances in swine biomedical model genomics. Int J Biol Sci. 2007;3(3):179–84.
Picard B, Berri C, Lefaucheur L, Molette C, Sayd T, Terlouw C. Skeletal muscle proteomics in livestock production. Brief Funct Genomics. 2010;9(3):259–78.
Picard B, Lefaucheur L, Berri C, Duclos MJ. Muscle fibre ontogenesis in farm animal species. Reprod Nutr Dev. 2002;42(5):415–31.
Ashmore CR, Addis PB, Doerr L. Development of muscle fibers in the fetal pig. J Anim Sci. 1973;36(6):1088–93.
Tang Z, Li Y, Wan P, Li X, Zhao S, Liu B, Fan B, Zhu M, Yu M, Li K. LongSAGE analysis of skeletal muscle at three prenatal stages in Tongcheng and Landrace pigs. Genome Biol. 2007;8(6):R115.
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.
Lefaucheur L, Ecolan P, Lossec G, Gabillard JC, Butler-Browne GS, Herpin P. Influence of early postnatal cold exposure on myofiber maturation in pig skeletal muscle. J Muscle Res Cell Motil. 2001;22(5):439–52.
Hambrecht E, Eissen JJ, Newman DJ, Smits CHM, Verstegen MWA, den Hartog LA. Preslaughter handling effects on pork quality and glycolytic potential in two muscles differing in fiber type composition. J Anim Sci. 2005;83(4):900–7.
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.
Scheffler TL, Gerrard DE. Mechanisms controlling pork quality development: the biochemistry controlling postmortem energy metabolism. Meat Sci. 2007;77(1):7–16.
Pette D, Staron RS. Myosin isoforms, muscle fiber types, and transitions. Microsc Res Tech. 2000;50(6):500–9.
Larzul C, Lefaucheur L, Ecolan P, Gogue J, Talmant A, Sellier P, Le Roy P, Monin G. Phenotypic and genetic parameters for longissimus muscle fiber characteristics in relation to growth, carcass, and meat quality traits in Large White pigs. J Anim Sci. 1997;75(12):3126–37.
Rehfeldt C, Te Pas MF, Wimmers K, Brameld JM, Nissen PM, Berri C, Valente LM, Power DM, Picard B, Stickland NC, et al. Advances in research on the prenatal development of skeletal muscle in animals in relation to the quality of muscle-based food. I. Regulation of myogenesis and environmental impact. Animal. 2011;5(5):703–17.
Shi X, Chen S, Zheng H, Wang L, Wu Y. Abnormal DNA methylation of imprinted loci in human preimplantation embryos. Reprod Sci. 2014;21(8):978–83.
Carvalho RN, Lettieri T. Proteomic analysis of the marine diatom Thalassiosira pseudonana upon exposure to benzo(a)pyrene. BMC Genomics. 2011;12:11.
Hakimov HA, Walters S, Wright TC, Meidinger RG, Verschoor CP, Gadish M, Chiu DKY, Stromvik MV, Forsberg CW, Golovan SP. Application of iTRAQ to catalogue the skeletal muscle proteome in pigs and assessment of effects of gender and diet dephytinization. Proteomics. 2009;9(16):4000–16.
Bjarnadottir SG, Hollung K, Hoy M, Bendixen E, Codrea MC, Veiseth-Kent E. Changes in protein abundance between tender and tough meat from bovine Longissimus thoracis muscle assessed by isobaric Tag for Relative and Absolute Quantitation (iTRAQ) and 2-dimensional gel electrophoresis analysis. J Anim Sci. 2012;90(6):2035–43.
J-q L, Chen Z-M, Liu D-W, X-h L, Sun B-L, Ling F, Zhang H, Chen Y-S. Genetic effects of IGF-1 gene on the performance in LandraceXLantang pig resource population. Acta Genet Sin. 2003;30(9):835–9.
Newcom DW, Stalder KJ, Baas TJ, Goodwin RN, Parrish FC, Wiegand BR. Breed differences and genetic parameters of myoglobin concentration in porcine longissimus muscle. J Anim Sci. 2004;82(8):2264–8.
Liu C, Guo Y, Lu T, Li X, Guan W, Ma Y. Establishment and genetic characteristics analysis of in vitro culture a fibroblast cell line derived from Wuzhishan miniature pig. Cryobiology. 2006;68(2):281–7.
Schiaffino S, Reggiani C. Molecular diversity of myofibrillar proteins: Gene regulation and functional significance. Physiol Rev. 1996;76(2):371–423.
Ashmore CR, Tompkins G, Doerr L. Postnatal development of muscle fiber types in domestic animals. J Anim Sci. 1972;34(1):37–41.
Wang JQ, Li X, Yang XJ, Sun QW, Huang RH, Xing J, Zhao RQ. Maternal dietary protein induces opposite myofiber type transition in Meishan pigs at weaning and finishing stages. Meat Sci. 2011;89(2):221–7.
Lefaucheur L, Milan D, Ecolan P, Le Callennec C. Myosin heavy chain composition of different skeletal muscles in Large White and Meishan pigs. J Anim Sci. 2004;82(7):1931–41.
Lefaucheur L, Edom F, Ecolan P, Butler-Browne GS. Pattern of muscle fiber type formation in the pig. Dev Dyn. 1995;203(1):27–41.
Xu XW, Qiu HF, Du ZQ, Fan B, Rothschild MF, Yuan F, Liu B. Porcine CSRP3: polymorphism and association analyses with meat quality traits and comparative analyses with CSRP1 and CSRP2. Mol Biol Rep. 2010;37(1):451–9.
Wang H, Jurado KA, Wu XL, Shun MC, Li X, Ferris AL, Smith SJ, Patel PA, Fuchs JR, Cherepanov P, et al. HRP2 determines the efficiency and specificity of HIV-1 integration in LEDGF/p75 knockout cells but does not contribute to the antiviral activity of a potent LEDGF/p75-binding site integrase inhibitor. Nucleic Acids Res. 2012;40(22):11518–30.
Harada A, Okada S, Konno D, Odawara J, Yoshimi T, Yoshimura S, Kumamaru H, Saiwai H, Tsubota T, Kurumizaka H, et al. Chd2 interacts with H3.3 to determine myogenic cell fate. Embo J. 2012;31(13):2994–3007.
Harada A, Maehara K, Sato Y, Konno D, Tachibana T, Kimura H, Ohkawa Y. Incorporation of histone H3.1 suppresses the lineage potential of skeletal muscle. Nucleic Acids Res. 2015;43(2):775–86.
Guerci A, Lahoute C, Hebrard S, Collard L, Graindorge D, Favier M, Cagnard N, Batonnet-Pichon S, Precigout G, Garcia L, et al. Sri-dependent paracrine signals produced by myofibers control satellite cell-mediated skeletal muscle hypertrophy. Cell Metab. 2012;15(1):25–37.
Torrado M, Nespereira B, Lopez E, Centeno A, Castro-Beiras A, Mikhailov AT. ANKRD1 specifically binds CASQ2 in heart extracts and both proteins are co-enriched in piglet cardiac Purkinje cells. J Mol Cell Cardiol. 2005;38(2):353–65.
Davies KE, Briand P, Ionasescu V, Ionasescu G, Williamson R, Brown C, Cavard C, Cathelineau L. Gene for OTC: characterisation and linkage to Duchenne muscular dystrophy. Nucleic Acids Res. 1985;13(1):155–65.
Voytik SL, Przyborski M, Badylak SF, Konieczny SF. Differential expression of muscle regulatory factor genes in normal and denervated adult rat hindlimb muscles. Dev Dyn. 1993;198(3):214–24.
Hughes SM, Taylor JM, Tapscott SJ, Gurley CM, Carter WJ, Peterson CA. Selective accumulation of MyoD and myogenin mRNAs in fast and slow adult skeletal muscle is controlled by innervation and hormones. Development. 1993;118(4):1137–47.
Cagnazzo M, te Pas MFW, Priem J, de Wit AAC, Pool MH, Davoli R, Russo V. Comparison of prenatal muscle tissue expression profiles of two pig breeds differing in muscle characteristics. J Anim Sci. 2006;84(1):1–10.
Nicholson LV, Davison K, Falkous G, Harwood C, O’Donnell E, Slater CR, Harris JB. Dystrophin in skeletal muscle. I. Western blot analysis using a monoclonal antibody. J Neurol Sci. 1989;94(1–3):125–36.
Huang W, Khatib H. Comparison of transcriptomic landscapes of bovine embryos using RNA-Seq. BMC Genomics. 2010;11:10.
Gupta MK, Jang JM, Jung JW, Uhm SJ, Kim KP, Lee HT. Proteomic analysis of parthenogenetic and in vitro fertilized porcine embryos. Proteomics. 2009;9(10):284–2860.
Nie JZ, An L, Miao K, Hou ZC, Yu Y, Tan K, Sui LL, He SZ, Liu Q, Lei X, Wu ZH, Tian JH. Comparative analysis of dynamic proteomic profiles between in vivo and in vitro produced mouse embryos during postimplantation period. J Proteome Res. 2013;12(9):3843–56.
Karp NA, Spencer M, Lindsay H, O’Dell K, Lilley KS. Impact of replicate types on proteomic expression analysis. J Proteome Res. 2005;4(5):1867–71.
Ma CY, Zhou JW, Chen GX, Bian YW, Lv DW, Li XH, Wang ZM, Yan YM. iTRAQ-based quantitative proteome and phosphoprotein characterization reveals the central metabolism changes involved in wheat grain development. BMC Genomics. 2014;15:1029.
This research was supported by the National Natural Science Foundation of China (31272417), Joint Funds of NSFC-Guangdong (U1201213) and China Agriculture Research System (CASR-36). We thank the Beijing Genomics Institute (BGI) Shenzhen for providing us with technical assistance in iTRAQ and bioinformatics analysis.
The authors declare that they have no competing interests.
DM and YC conceived and designed the experiments. XZ performed the experiments and wrote the paper. JP and XW contributed to experimental animals’ artificial insemination and sample collection. XL helped with experiments and sample collection. HC, XZ and ZY helped interpret the results. All authors reviewed the manuscript and agreed with its publication.
Additional file 1: Figure S1.
Basic information of iTRAQ. (A) Basic information statistics; (B) Distribution of protein’s sequences coverage; (C) Peptide length distribution; (D) Peptide number distribution; (E) Protein mass distribution; (F) Length of peptide (in amino acids); (G) Spectrogram matching mass error distribution. (TIF 501 kb)
Additional file 2: Figure S2.
(A) Hematoxylin-eosin staining and (B) immunohistochemical staining of LR and WZS in different stages. MYHC-fast (Abcam, ab7784) antibody was used in (B) which was marked by red. Scale bar = 50 μm. (TIF 7013 kb)
Additional file 3: Table S1.
Complete list of up and down-regulated proteins identified by iTRAQ. (XLS 1326 kb)
Additional file 4: Figure S3.
Protein ratio distribution between LR and WZS at 4 stages. (TIF 143 kb)
Additional file 5: Table S2.
List of significantly changed DEPs identified by iTRAQ between LR and WZS at 4 stages. (XLSX 20 kb)
Additional file 6: Figure S4.
GO categories of identified proteins. (TIF 1081 kb)
Additional file 7: Figure S5.
COG function classification of identified proteins. (TIF 416 kb)
Additional file 8: Table S3.
Detail of GO, GO2proteins and protein2GO identified by iTRAQ. (XLS 3474 kb)
Additional file 9: Table S4.
Terms from the muscle process ontology with p-value as good or better than 1 between LR and WZS at 4 stages. (XLSX 11 kb)
Additional file 10: Table S5.
Detail of COG identified by iTRAQ. (XLS 7958 kb)
Additional file 11: Table S6.
Complete list of significantly changed myogenic DEPs identified by iTRAQ between LR and WZS at 4 stages. (XLSX 45 kb)
Additional file 12: Table S7.
Detail of some important myogenic proteins’ interaction network identified by STRING (http://string-db.org/). (XLSX 14 kb)
Additional file 13:
ARRIVE checklist. (PDF 1104 kb)
Additional file 14: Figure S6.
Embryo picture of 21 dpc (LR&WZS). (TIF 164 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Zhang, X., Chen, Y., Pan, J. et al. iTRAQ-based quantitative proteomic analysis reveals the distinct early embryo myofiber type characteristics involved in landrace and miniature pig. BMC Genomics 17, 137 (2016). https://doi.org/10.1186/s12864-016-2464-1
- Embryonic myogenesis
- Muscle proteome
- Quantitative proteomics