Proteomic Assessment of the Relevant Factors Affecting Pork Meat Quality Associated with Longissimus dorsi Muscles in Duroc Pigs

Article information

Asian-Australas J Anim Sci. 2016;29(11):1653-1663
Publication date (electronic) : 2016 March 30
doi : https://doi.org/10.5713/ajas.16.0050
1Department of Dental Pharmacology, School of Dentistry and Institute of Dental Bioscience, BK21 plus, Chonbuk National University, Jeonju 651-756, Korea
2National Marine Biodiversity Institute of Korea, Seocheon 33662, Korea
3Pohang Center for Evaluation of Biomaterials, Pohang 37668, Korea
4Aging Research Institute, Korea Research Institute of Bioscience & BioTechnology, Daejeon 34141, Korea
5Division of Animal Genomics and Bioinformatics, National Institute of Animal science, Rural Development Administration, Suwon, 441-706, Korea
6Department of Animal Science and Technology, Sunchon National University, Suncheon 540-742, Korea
7Department of Animal Science, College of Agricultural and Life Science, Chonbuk National University, Jeonju 651-756, Korea
8Department of Pharmacy, Keimyung University, Daegu 704-701, Korea
9Department of Biomedical Engineering, Konyang University, Daejeon 35365, Korea
10Department of Applied Bioscience, CHA University, Seongnam 463-836, Korea
11Department of Prosthodontics, School of Dentistry and Institute of Oral Bio-Science and Research Institute of Clinical Medicine, Chonbuk National University, Jeonju 561-756, Korea
12Cluster for Craniofacial Development and Regeneration Research, Institute of Oral Biosciences and School of Dentistry, Chonbuk National University, Jeonju 561-756, Korea
13Department of Pharmacy, College of Pharmacy and Natural Medicine Research Institute, Mokpo National University, Mokpo 534-729, Korea
14The China-US (Henan) Hormel Cancer Institute, Zhengzhou, Henan 127, China
*Corresponding Authors: Jung-Hyun Shim. Tel: +82-61-450-2684, Fax: +82-61-450-2684, E-mail: s1004jh@gmail.com/Jung-Il Chae. Tel: +82-63-270-4024, Fax: +82-63-270-4037, E-mail: jichae@jbnu.ac.kr
aThese authors contributed equally to this work.
Received 2016 January 19; Revised 2016 February 21; Accepted 2016 March 29.

Abstract

Meat quality is a complex trait influenced by many factors, including genetics, nutrition, feeding environment, animal handling, and their interactions. To elucidate relevant factors affecting pork quality associated with oxidative stress and muscle development, we analyzed protein expression in high quality longissimus dorsi muscles (HQLD) and low quality longissimus dorsi muscles (LQLD) from Duroc pigs by liquid chromatographytandem mass spectrometry (LC-MS/MS)–based proteomic analysis. Between HQLD (n = 20) and LQLD (n = 20) Duroc pigs, 24 differentially expressed proteins were identified by LC-MS/MS. A total of 10 and 14 proteins were highly expressed in HQLD and LQLD, respectively. The 24 proteins have putative functions in the following seven categories: catalytic activity (31%), ATPase activity (19%), oxidoreductase activity (13%), cytoskeletal protein binding (13%), actin binding (12%), calcium ion binding (6%), and structural constituent of muscle (6%). Silver-stained image analysis revealed significant differential expression of lactate dehydrogenase A (LDHA) between HQLD and LQLD Duroc pigs. LDHA was subjected to in vitro study of myogenesis under oxidative stress conditions and LDH activity assay to verification its role in oxidative stress. No significant difference of mRNA expression level of LDHA was found between normal and oxidative stress condition. However, LDH activity was significantly higher under oxidative stress condition than at normal condition using in vitro model of myogenesis. The highly expressed LDHA was positively correlated with LQLD. Moreover, LDHA activity increased by oxidative stress was reduced by antioxidant resveratrol. This paper emphasizes the importance of differential expression patterns of proteins and their interaction for the development of meat quality traits. Our proteome data provides valuable information on important factors which might aid in the regulation of muscle development and the improvement of meat quality in longissimus dorsi muscles of Duroc pigs under oxidative stress conditions.

INTRODUCTION

Variation in meat quality traits is a well-known problem. Meat quality traits are closely related to biological traits of live animal. Hence, biological sciences including genetics, physiology, cell biology, and biochemistry have been widely employed for decades to characterize the biological mechanisms behind major variability of meat quality traits (Bendixen, 2005). Basic knowledge of these mechanisms is essential to reduce the variation in meat quality traits such as tenderness, water-holding capacity, and color. They are also important to understand the physiology of meat animals, especially on muscle growth and development (Lametsch et al., 2002; Hwang et al., 2005). Understanding and changes related to physiochemical factors, genotypes, and many other factors influence postmortem metabolism (Monin et al., 1995; Brocks et al., 1998; Wheeler et al., 2005). Some previous studies have indicated that meat quality is determined by postmortem muscle metabolism (Pette, 2002; Spangenburg and Booth, 2003). At slaughter, muscles become deprived of oxygen as the circulatory system shuts down. This lack of oxygen results in a shift to glycolytic (anaerobic) metabolism and a buildup of lactic acid, causing a drop in muscle pH (Frisby et al., 2005). Accelerated postmortem glycolysis reduces pH and increases temperature within muscle, resulting in excessive protein denaturation and inferior meat quality (Julve et al., 2000). Although extensively researched, the underlying mechanisms of many different meat quality traits are far from well understood due to many factors affecting the quality of meat (Mullen et al., 2006; Hollung et al., 2007). The proteome expressed from the genome is influenced by environmental conditions. Proteome is the molecular link between the genome and the functional quality characteristics of the meat. Therefore, proteomics is a promising and powerful tool in meat science (Lametsch and Bendixen, 2001; Morzel et al., 2004; Jia et al., 2006; Sayd et al., 2006). However proteomics has been, and still are, used in numerous studies on skeletal muscle (Picard et al., 2010).

In this study, we focus on its use in the study of livestock muscle development and meat quality with a focus on the differential expression patterns of proteins and their interactions for the development of meat quality traits.

MATERIALS AND METHODS

Animals and sample collection

The meat quality characteristics were assessed from 200 randomly selected great grandparent Duroc pigs raised from October 2011 to March 2012 for one production cycle. The live weight ranged from 100 to 120 kg. The carcasses were kept in a freezer (0°C) for 24 h after slaughtering. The frozen carcasses were thawed, deboned, and trimmed. The left side loin was transferred to the laboratory and placed in a deep-freezer (−45°C) for analysis

Gel electrophoresis and silver staining

High quality longissimus dorsi muscles (HQLD) and low quality longissimus dorsi muscles (LQLD) tissues were collected from Duroc pigs. Total protein isolation was performed using PRO-PREP protein extraction solution (iNtRON biotechnology, Sungnam, Korea) according to the manufacturer’s instructions. Concentrations of eluted proteins were measured using Pierce BCA Protein Assay Kit (Thermo scientific, Rockford, IL, USA). Equal amounts of protein samples were precipitated with cold acetone. Protein pellets dissolved in 1× sodium dodecyl sulfate (SDS) sample buffer were separated by 8% and 12% SDS-polyacrylamide gel electrophoresis (PAGE). Following SDS-PAGE, protein spots were visualized using protocols described in PlusOne Silver staining kit (GE Healthcare Bio-Sciences, Uppsala, Sweden). The complete protocol was followed to analyze gels. To prepare gels, the protocol was modified so that glutaraldehyde was omitted from the sensitization step and formaldehyde was omitted from the silver reaction step (Yan et al., 2000). Silver-stained gels were scanned (UMAX PowerLook 2100KL Imaging system, UMAX, Taiwan) and protein profiles were compared.

Liquid chromatography-tandem mass spectrometry

The resulting tryptic peptides were separated and analyzed using reversed-phase capillary high-performance liquid chromatography directly coupled to a Thermo LTQ Orbitrap mass spectrometer using published procedure described by Zuo et al. (2001) with slight modifications. Briefly, a 0.075×20 mm trapping column and a 0.075×120 mm resolving column were packed with C18AQ 218MS low formic acid C18 beads (5 μm in size, 200Å pore size; C18AQ, Michrom BioResources, Auburn, CA, USA) and placed in-line. Peptides were bound to the trapping column for 10 min with 2% (vol/vol) aqueous acetonitrile containing 0.1% (vol/vol) formic acid. The bound peptides were then eluted with a 67 min gradient of 2% to 90% (vol/vol) acetonitrile containing 0.1% (vol/vol) formic acid at a flow rate of 0.2 μL/min. For tandem mass spectrometry, the full mass scan range mode was set at m/z = 50 to 2,000 Da. After determining the charge states of the ion zoom scans, product ion spectra were acquired in MS/MS mode with relative collision energy of 55%. The individual spectra from MS/MS were processed using Protein discoverer 2.1 software (Thermo scientific, USA). The generated peak list files were used to query either the MSDB or the NCBI database using the MASCOT program (http://www.matrixscience.com). We considered modifications of methionine and cysteine, peptide mass tolerance at 2 Da, MS/MS ion mass tolerance at 0.8 Da, allowance of missed cleavage at 2, and charge states (namely, +1, +2, and +3). Only significant hits as defined by MASCOT probability analysis were initially considered.

Cell culture

Mitotic C2C12 mouse myoblasts were obtained from Chonbuk University (Jeonju, Korea). C2C12 were passaged as subconfluent monolayers in growth medium (GM) using Dulbecco’s modified Eagle’s medium (Invitrogen, Grand Island, NY, USA) supplemented with 20% fetal bovine serum, 200 mM L-glutamine, 10 units/mL penicillin, and 10 μg/mL streptomycin. Confluent (90%) myoblasts were differentiated into myotubes by culturing the cells in differentiation medium (DM) with Dulbecco’s modified Eagle’s medium supplemented with 2% horse serum.

3-(4, 5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium (MTS) assay

The effects of H2O2 on cell viability were estimated using MTS Assay Kit (Promega, Madison, WI, USA). C2C12 cells were seeded into a 96-well plate for 24 h and treated with H2O2 (12.5 μM to 1 mM) for 24 h or 48 h. MTS solution was added to the plates and incubated at 37°C with 5% CO2 for 2 h. Absorbance at 490 nm was recorded using a GloMax-Multi Microplate Multimode Reader (Promega, USA).

Immunocytochemistry

To detect myosin heavy chain (MYH), myogenic differentiation (MyoD), and myogenin (Myog), cells were blocked with 1% bovine serum albumin and incubated with monoclonal anti-MYH (B-5), anti-MyoD (E-1), or anti-Myog (M-225) antibody at 4°C overnight (Santa Cruz Biotechnology, Santa Cruz, CA, USA). Probed cells were reacted with a 488-conjugated anti-mouse or 594-conjugated anti-rabbit secondary antibody. For nucleus staining, cells were treated with mounting medium with 4′-6-diamidino-2-phenylindole (DAPI) (Vector Laboratories, Inc. Burlingame, CA, USA). The cells were visualized using a FluoView confocal laser microscope (Fluoview FV10i, Olympus Corporation, Tokyo, Japan).

Reverse transcription polymerase chain reaction and real-time polymerase chain reaction analysis

Total RNA isolation was performed using TRIzol reagent (Invitrogen, USA) according to the manufacturer's instructions. Briefly, total RNA levels were quantified by absorbance at 260 nm. RNA integrity was evaluated by 1.2% (w/v) agarose gel. Total RNA (2 μg amounts) was reverse-transcribed into cDNA using QuantiTect Reverse Transcription Kit (Qiagen, Chatsworth, CA, USA) according to the manufacturer’s instructions. Real-time polymerase chain reaction analysis (PCR) was performed with SYBR green Premix Ex Taq II (Takara, Dalian, China) using Applied Biosystems StepOne Plus Real-time PCR System (Applied Biosystems, Carlsbad, CA, USA). Relative quantification analysis was performed using the comparative Ct (2 (−ΔΔCT)) method. The expression of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and β-actin was used as endogenous control for the detection of mRNA expression levels. Primers used in the study were listed in Table 4.

Primer sequences used to generate templates for reverse transcription polymerase chain reaction

Kinetic determination of LDH activity

Commercially available kits for lactate dehydrogenase (DLDH-100, QuantiChrom Lactate Dehydrogenase Kit, Gentaur Molecular, Hayward, CA, USA) were used according to the manufacturers’ instructions (Stentz et al., 2010).

RESULTS AND DISCUSSION

Animals and phenotypes

The average weight before and after slaughtering of the GGP pigs were 98.23 kg and 89.98 kg, respectively. Difference of 10 kg before and after slaughtering was recorded. Important economic traits such as lean percent and eye muscle area showed an average value of 54.74% and 26.03 cm2, respectively (Table 1).

Phenotypic record of meat quality traits in Duroc pigs

Protein profiles in HQLD and LQLD from Duroc pigs

Meat qualities were evaluated by Korea Institute for Animal Products Quality Evaluation (KAPE) authorized by South Korean government to perform animal products grading service. The normality test was applied to show normal distribution of the traits. The highest and lowest meat grades for pH and water holding capacity were identified from the sample which accounted for 10% (20 heads) of the total population (Table 2). To obtain a comprehensive overview of protein components in HQLD and LQLD from 12 individuals, protein profiles of whole lysate of HQLD and LQLD separated by 8% and 12% SDS-PAGE were assessed by silver-stained image analysis (Figure 1A). Patterns of total protein components in whole lysate of HQLD and LQLD were similar among the individual six groups. However, significantly fewer proteins (i.e. band a and b) were expressed in HQLD compared to in LQLD. Two different protein spots were identified by a mass spectrometric analysis. Myosin binding protein C (MYBPC2) expressed during skeletal muscle development (Gurnett et al., 2010) had higher expression in HQLD than in LQLD. However, lactate dehydrogenase A (LDHA) catalyzing the conversion of pyruvate to lactate during glycolysis (Fan et al., 2011) had higher expression in LQLD than in HQLD (Figure 1B).

High- or low-meat quality traits of longissimus dorsi muscles in Duroc pigs

Figure 1

Protein profiles of high quality longissimus dorsi muscles (HQLD) and low quality longissimus dorsi muscles (LQLD) from Duroc pigs by image analysis. (A) The overall patterns of total protein bands from individuals. All gels were visualized by sliver staining. (B) Two different protein spots were identified by a mass spectrometric analysis. a*: MYBPC2, myosin binding protein C; b#: LDHA, lactate dehydrogenase A.

Protein identification and gene ontological classification by LC-MS/MS-based proteomic analysis

Next, we performed liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based proteomic analysis to elucidate the proteins involved in longissimus dorsi muscle (i.e. HQLD and LQLD) properties involved in meat quality. Among the total 24 proteins identified, 10 and 14 proteins were confirmed to be highly expressed in HQLD and LQLD, respectively. All identified proteins were clustered into the following seven categories (Figure 2A) based on information obtained from DAVID gene ontology (GO) database (http://david.abcc.ncifcrf.gov) and UniProt (http://www.uniprot.org): Catalytic activity (31%), ATPase activity (19%), oxidoreductase activity (13%), cytoskeletal protein binding (13%), actin binding (12%), calcium ion binding (6%) and structural constituent of muscle (6%) (Supplementary Table S1). The GO analysis was performed using DAVID Bioinformatics Resources 6.7 categories for both molecular function (MF) and biological process (BP). Depending on the MF in which the proteins were involved, they were categorized into the following three groups (Figure 2B): Cytoskeletal protein binding (40%), actin binding (40%), and structural constituent of muscle (20%) (Supplementary Table S2). Depending on the BP in which the proteins were involved, they were categorized into the following five groups (Figure 2C): Primary metabolic process (26%), cellular metabolic process (26%), catabolic process (18%), nitrogen compound metabolic process (19%), and oxidation reduction (11%) (Supplementary Table S3). The expression changes of the up- and down-regulated proteins in HQLD and LQLD of Duroc pigs were summarized in Table 3. LDHA was selected and subjected to further analysis by LDH activity assay and in vitro study of myogenesis under oxidative stress conditions.

Figure 2

Ontological classifications of differentially regulated proteins in high quality longissimus dorsi muscles (HQLD) and low quality longissimus dorsi muscles (LQLD) from Duroc pigs. Of the total 24 identified proteins, 10 and 14 proteins were highly expressed in HQLD and LQLD, respectively. (A) The identified proteins were clustered into 7 categories based on information obtained from DAVID gene ontology (GO) database; (B) The identified proteins were clustered into 3 categories based on their molecular function; (C) The identified proteins were clustered into 5 categories based on their biological processes.

List of differentially expressed genes on high- and low quality of longissimus dorsi muscles in duroc pigs

Gross changes in C2C12 myoblasts in response to myogenic differentiation

C2C12 myoblasts serve as an experimentally tractable model system for investigating the molecular basis of skeletal muscle cell specification and development (Kislinger et al., 2005). A temporally well-defined myogenic differentiation program can be selectively triggered in cultured C2C12 myoblasts upon withdrawal of GM and mitogens (Gramolini and Jasmin, 1999). When switched to DM, mitotic C2C12 myoblasts rapidly cease proliferation and initiate a synchronously terminal differentiation program (Figure 3A). To investigate the patterns of protein expression and efficiency of myotube formation during myogenesis under oxidative stress condition, undifferentiation C2C12 cells were treated with various concentrations of H2O2 (12.5 μM to 1 mM). Cytotoxicity was negligible with 200 μM H2O2. However, up to 1 mM H2O2 did reduce viability which was confirmed by 3-(4, 5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium) (MTS) assay (Figure 3B). Control cells exhibited striking morphological changes over the course of 3 to 7 days, eventually fusing into mature multinucleated myotubes (i.e. by day 5). However, H2O2-treated cells exhibited thread-like shape without fusing into mature multinucleated myotubes (Figure 3C). Myotubes were identified by immunocytochemistry with anti-myosin heavy chain (MYH) antibody (Figure 3D).

Figure 3

Gross changes and cell viability in myogenic differentiation under oxidative stress condition. (A) Morphological changes of myoblast during myogenesis. Arrow head and dashed lines indicating myoblast with mono nucleus and myotube with three or more nuclei, respectively; (B) Cytotoxic effects of H2O2 in C2C12 treated with the various concentrations of H2O2 (12.5 μM to 1 mM). Data represented mean value of three independent experiments; bars, standard deviation. * Significantly different compared to untreated conditions (24 h); #, significantly different compared to untreated conditions (48 h) (n = 3, p<0.0001); (C) Light microscopy-based images of undifferentiated (day 0) proliferating C2C12 myoblasts and differentiating cells at various time points (up to day 7) under normal and oxidative stress conditions (H2O2 [200 μM]); (D) Immunocytochemistry presented the MYH expression levels of undifferentiated (day 0) and differentiated (day 5). GM, growth medium; DM, differentiation medium. (E) Oxidative stress induced down-regulation of Myog gene expression and continuous up-regulation of MyoD gene expression with glyceraldehyde-3-phosphate dehydrogenase (GAPDH) as a control; (F) Oxidative stress induced down-regulation of Myog protein expression and continuous up-regulation of MyoD protein expression with β-actin as a control. Student’s t test was performed to evaluate statistical significance (ns, non-significant; *** p<0.0001; mean±standard error of the mean; n = 3).

In vitro model of myogenesis under oxidative stress condition

The expression of myogenic regulatory factors consisting of MyoD and Myog characterizes various phases of skeletal muscle development, including myoblast proliferation, cell-cycle exit, cell fusion, and the maturation of myotubes to form myofibers (Lee et al., 2014). MyoD, the chief regulatory molecule of myogenic differentiation (Langen et al., 2004), plays an important role in cell cycle exit of differentiating myoblasts (Guo et al., 1995; Halevy et al., 1995) Terminal differentiation of myoblast, driven by expression of Myog, is essential for the formation of functional multinucleated myofibers (Liu et al., 2012). In order to study the patterns of MyoD and Myog expression (gene and protein levels) depending on oxidative stress condition during myogenesis, we used western blot and RT-PCR analyses. Our data revealed that, among various genes subjected to comparison, MyoD had significant (p<0.0001) higher expression under oxidative stress condition than under normal condition on day 3. The expression of MyoD on day 3 was not significantly different from that on day 2 under oxidative stress condition. Myog had significant (p<0.0001) higher expression under normal condition than under oxidative stress condition at day 6. Under oxidative stress condition, the expression of Myog was significantly (p<0.05) decreased at day 6 compared to that at day 5 (Figure 3E). The mRNA expression levels for selected genes were analyzed by quantitative real-time PCR with specific primer (Table 4). Western blot analysis data revealed that, of different proteins subjected to comparison, MyoD had significantly (p<0.0001) higher expression under oxidative stress condition than under normal condition at day 2. Under oxidative stress condition, the expression of MyoD at 12 h was not significantly different from that at 6 h. Moreover, Myog had significantly (p<0.0001) higher expression under normal condition than under oxidative stress condition at day 4. Under oxidative stress condition, the expression of Myog was significantly more decreased (p = 0.0016) at day 4 than at day 3 (Figure 3F). These results indicate that H2O2-induced oxidative stress inhibits myogenesis through the down-regulation of Myog and the continuous up-regulation of MyoD.

Relationships between LDHA gene expression and myogenesis

Previous studies reported that porcine myogenic differentiation 1 (MyoD1) gene has been mapped to swine chromosome 2p14-p17 which is involved in the regulation of the proliferation and differentiation of skeletal muscle cells. LDHA genes mapped close to MyoD are involved in energy metabolism and protein transport processes. LDHA genes might play important roles in muscle development (Qiu et al., 2010). However, little is known about porcine LDHA genes. Therefore, we determined the relationships between LDHA gene expression and myogenesis under normal and oxidative stress condition. To investigate whether oxidative stress regulated LDHA expression at genetic levels, we used quantitative real-time PCR. The mRNA expression levels of selected genes were subjected to quantitative real-time PCR with specific primers (Table 4). Figure 4A showed that LDHA genes were increased up to day 5 during myogenesis under normal and oxidative stress condition. LDHA expression was significantly (p = 0.0088) higher under oxidative stress condition. These results indicate that up-regulated LDHA genes induced by oxidative stress might play dysfunctional roles in myogenesis.

Figure 4

Gene and protein expression levels on high quality longissimus dorsi muscles (HQLD) and low quality longissimus dorsi muscles (LQLD) depending on oxidative stress. Myoblast and myotube cDNA and protein were determined by real-time PCR analysis and western blot analysis, respectively. (A) The quantitative differences of lactate dehydrogenase A (LDHA) at the transcriptional level were measured by real-time PCR during myogenesis with β-actin as control; (B) Cytotoxic effects of resveratrol (RSV) in C2C12 treated with various concentrations of RSV (6.25 μM to 100 μM). Data represented mean value of three independent experiments; bars, standard deviation. * Significantly different compared with untreated conditions (24 h); #, significantly different compared to untreated conditions (48 h) (n = 3, p<0.0001); (C) Quantitative differences of LDHA at the transcriptional level were measured by real-time PCR in RSV-treated C2C12 during myogenesis.

Antioxidant properties of resveratrol

Resveratrol (RSV), a well-known phytocompound and food component, has antioxidative and multifunctional bioactivities (Wu et al., 2013). Previous studies have reported that RSV in skeletal muscle acts on protein catabolism and muscle function and confers resistance against oxidative stress, injury, and cell death. However, its action mechanisms and protein targets in myogenesis process are not completely understood (Montesano et al., 2013). Therefore, we determined the effect of RSV on LDHA gene expression in myogenesis under oxidative stress condition. C2C12 cells were treated with various concentrations of RSV (6.25 μM to 100 μM). Cytotoxicity was negligible with RSV (25 μM) under normal and oxidative stress conditions. However, up to 100 μM RSV did reduce viability which was confirmed by MTS assay (Figure 4B). However, after RSV treatment, there was no significant difference in mRNA expression of LDHA between normal condition and oxidative stress condition (Figure 4C).

Confirmation of lactate dehydrogenase activity

Elevation of plasmatic LDH levels are characteristic responses to strenuous exercise which are often used as indicators of muscle damage (Kanter et al., 1988; Bouzid et al., 2014). However, difference of LDH activity between HQLD and LQLD of Duroc pigs is not well determined. Our data showed that LDH activity was significantly (p = 0.0003) higher in LQLD than in HQLD of Duroc pigs (Figure 5A). Moreover, higher LDH activity was positively correlated with in vitro model of myogenesis under oxidative stress condition. In addition, LDH activity was significantly reduced by RSV treatment (Figure 5B). We also confirmed the patterns of MyoD and Myog expression under oxidative stress condition and RSV treatment during myogenesis by immunocytochemistry. Oxidative stress induced down-regulation of Myog and continuous up-regulation of MyoD. The down-regulation induced by oxidative stress was recovered by RSV treatment. The up-regulation of MyoD induced by the oxidative stress was reduced by the treatment of RSV in myogenesis (Figure 5C). In addition, there was a significant correlation between MyoD and Myog expression (gene and protein levels) under oxidative stress condition during myogenesis (Figure 5D). These results indicate that high activity of LDH by oxidative stress will result in dysfunction of myogenesis and that RSV treatment will result in its functional recovery.

Figure 5

Antioxidant effect of resveratrol in lactate dehydrogenase activity and myogenesis. (A) Difference of lactate dehydrogenase (LDH) activity between high quality longissimus dorsi muscles (HQLD) and low quality longissimus dorsi muscles (LQLD) from individual six groups of Duroc pigs. (B) In vitro model of myogenesis under oxidative stress condition and changes of lactate dehydrogenase activity by resveratrol (RSV) treatment. Student’s t test was performed to evaluate statistical significance (*** p<0.0001; mean±standard error of the mean; n = 3); (C) Patterns of myogenic differentiation (MyoD) and myogenin (Myog) expression depending on oxidative stress condition and RSV treatment during myogenesis by immunocytochemistry; (D) Patterns of MyoD and Myog expression (gene and protein levels) depending on oxidative stress condition and RSV treatment during myogenesis by reverse transcription polymerase chain reaction and western blot analysis.

In conclusion, our data demonstrated that up- or down-regulation of genes and proteins are involved in muscle development, muscle function, actin organization, oxidative stress, cell proliferation, cell differentiation, and cell growth. In this paper, differential expression patterns of genes and their interaction are found to be important for the development of meat quality traits. Our proteome data provided valuable information on differentially expressed genes (LDHA) and activity of LDH in HQLD and LQLD from Duroc pigs, which may aid in the regulation of muscle development. Our study provided experimental evidence for RSV as an important regulator to improve meat quality grades in porcine. However, further study is required to determine the relationship between differential expression of genes or proteins and their direct effects on meat quality.

ACKNOWLEDGMENTS

This work was carried out with the support of “Cooperative Research Program for Agriculture Science & Technology Development (Project No. PJ011682)” Rural Development Administration, Republic of Korea.

Notes

CONFLICT OF INTEREST

We certify that there is no conflict of interest with any financial organization regarding the material discussed in the manuscript.

References

Bendixen E. 2005;The use of proteomics in meat science. Meat Sci 71:138–149.
Bouzid MA, Hammouda O, Matran R, Robin S, Fabre C. 2014;Changes in oxidative stress markers and biological markers of muscle injury with aging at rest and in response to an exhaustive exercise. PLoS One 9:e90420.
Brocks L, Hulsegge B, Merkus G. 1998;Histochemical characteristics in relation to meat quality properties in the Longissimus Lumborum of fast and lean growing lines of Large White pigs. Meat Sci 50:411–420.
Fan J, Hitosugi T, Chung TW, Xie J, Ge Q, Gu TL, Polakiewicz RD, Chen GZ, Boggon TJ, Lonial S, Khuri FR, Kang S, Chen J. 2011;Tyrosine phosphorylation of lactate dehydrogenase A is important for NADH/NAD(+) redox homeostasis in cancer cells. Mol Cell Biol 31:4938–4950.
Frisby J, Raftery D, Kerry JP, Diamond D. 2005;Development of an autonomous, wireless pH and temperature sensing system for monitoring pig meat quality. Meat Sci 70:329–336.
Gramolini AO, Jasmin BJ. 1999;Expression of the utrophin gene during myogenic differentiation. Nucleic Acids Res 27:3603–3609.
Guo K, Wang J, Andres V, Smith RC, Walsh K. 1995;MyoD-induced expression of p21 inhibits cyclin-dependent kinase activity upon myocyte terminal differentiation. Mol Cell Biol 15:3823–3829.
Gurnett CA, Desruisseau DM, McCall K, Choi R, Meyer ZI, Talerico M, Miller SE, Ju JS, Pestronk A, Connolly AM, Druley TE, Weihl CC, Dobbs MB. 2010;Myosin binding protein C1: a novel gene for autosomal dominant distal arthrogryposis type 1. Hum Mol Genet 19:1165–1173.
Halevy O, Novitch BG, Spicer DB, Skapek SX, Rhee J, Hannon GJ, Beach D, Lassar AB. 1995;Correlation of terminal cell cycle arrest of skeletal muscle with induction of p21 by MyoD. Science 267:1018–1021.
Hollung K, Veiseth E, Jia X, Faergestad EM, Hildrum KI. 2007;Application of proteomics to understand the molecular mechanisms behind meat quality. Meat Sci 77:97–104.
Hwang IH, Park BY, Kim JH, Cho SH, Lee JM. 2005;Assessment of postmortem proteolysis by gel-based proteome analysis and its relationship to meat quality traits in pig longissimus. Meat Sci 69:79–91.
Jia X, Hollung K, Therkildsen M, Hildrum KI, Bendixen E. 2006;Proteome analysis of early post-mortem changes in two bovine muscle types: M. longissimus dorsi and M. semitendinosis. Proteomics 6:936–944.
Julve J, Escola-Gil JC, Marzal-Casacuberta A, Ordonez-Llanos J, Gonzalez-Sastre F, Blanco-Vaca F. 2000;Increased production of very-low-density lipoproteins in transgenic mice overexpressing human apolipoprotein A-II and fed with a high-fat diet. Biochim Biophys Acta 1488:233–244.
Kanter MM, Lesmes GR, Kaminsky LA, La Ham-Saeger J, Nequin ND. 1988;Serum creatine kinase and lactate dehydrogenase changes following an eighty kilometer race. Relationship to lipid peroxidation. Eur J Appl Physiol Occup Physiol 57:60–63.
Kislinger T, Gramolini AO, Pan Y, Rahman K, MacLennan DH, Emili A. 2005;Proteome dynamics during C2C12 myoblast differentiation. Mol Cell Proteomics 4:887–901.
Lametsch R, Bendixen E. 2001;Proteome analysis applied to meat science: characterizing postmortem changes in porcine muscle. J Agric Food Chem 49:4531–4537.
Lametsch R, Roepstorff P, Bendixen E. 2002;Identification of protein degradation during post-mortem storage of pig meat. J Agric Food Chem 50:5508–5512.
Langen RC, Van Der Velden JL, Schols AM, Kelders MC, Wouters EF, Janssen-Heininger YM. 2004;Tumor necrosis factor-alpha inhibits myogenic differentiation through MyoD protein destabilization. FASEB J 18:227–237.
Lee EJ, Malik A, Pokharel S, Ahmad S, Mir BA, Cho KH, Kim J, Kong JC, Lee DM, Chung KY, Kim SH, Choi I. 2014;Identification of genes differentially expressed in myogenin knock-down bovine muscle satellite cells during differentiation through RNA sequencing analysis. PLoS ONE 9:e92447.
Liu QC, Zha XH, Faralli H, Yin H, Louis-Jeune C, Perdiguero E, Pranckeviciene E, Munoz-Canoves P, Rudnicki MA, Brand M, Perez-Iratxeta C, Dilworth FJ. 2012;Comparative expression profiling identifies differential roles for Myogenin and p38alpha MAPK signaling in myogenesis. J Mol Cell Biol 4:386–397.
Monin G, Lambooy E, Klont R. 1995;Influence of temperature variation on the metabolism of pig muscle in situ and after excision. Meat Sci 40:149–158.
Montesano A, Luzi L, Senesi P, Mazzocchi N, Terruzzi I. 2013;Resveratrol promotes myogenesis and hypertrophy in murine myoblasts. J Transl Med 11:310.
Morzel M, Chambon C, Hamelin M, Sante-Lhoutellier V, Sayd T, Monin G. 2004;Proteome changes during pork meat ageing following use of two different pre-slaughter handling procedures. Meat Sci 67:689–696.
Mullen AM, Stapleton PC, Corcoran D, Hamill RM, White A. 2006;Understanding meat quality through the application of genomic and proteomic approaches. Meat Sci 74:3–16.
Pette D. 2002;The adaptive potential of skeletal muscle fibers. Can J Appl Physiol 27:423–448.
Picard B, Berri C, Lefaucheur L, Molette C, Sayd T, Terlouw C. 2010;Skeletal muscle proteomics in livestock production. Brief Funct Genomics 9:259–278.
Qiu H, Xu X, Fan B, Rothschild MF, Martin Y, Liu B. 2010;Investigation of LDHA and COPB1 as candidate genes for muscle development in the MYOD1 region of pig chromosome 2. Mol Biol Rep 37:629–636.
Sayd T, Morzel M, Chambon C, Franck M, Figwer P, Larzul C, Le Roy P, Monin G, Cherel P, Laville E. 2006;Proteome analysis of the sarcoplasmic fraction of pig semimembranosus muscle: implications on meat color development. J Agric Food Chem 54:2732–2737.
Spangenburg EE, Booth FW. 2003;Molecular regulation of individual skeletal muscle fibre types. Acta Physiol Scand 178:413–424.
Stentz R, Bongaerts RJ, Gunning AP, Gasson M, Shearman C. 2010;Controlled release of protein from viable Lactococcus lactis cells. Appl Environ Microbiol 76:3026–3031.
Wheeler TL, Cundiff LV, Shackelford SD, Koohmaraie M. 2005;Characterization of biological types of cattle (Cycle VII): Carcass, yield, and longissimus palatability traits. J Anim Sci 83:196–207.
Wu RE, Huang WC, Liao CC, Chang YK, Kan NW, Huang CC. 2013;Resveratrol protects against physical fatigue and improves exercise performance in mice. Molecules 18:4689–4702.
Yan JX, Wait R, Berkelman T, Harry RA, Westbrook JA, Wheeler CH, Dunn MJ. 2000;A modified silver staining protocol for visualization of proteins compatible with matrix-assisted laser desorption/ionization and electrospray ionization-mass spectrometry. Electrophoresis 21:3666–3672.
Zuo X, Echan L, Hembach P, Tang HY, Speicher KD, Santoli D, Speicher DW. 2001;Towards global analysis of mammalian proteomes using sample prefractionation prior to narrow pH range two-dimensional gels and using one-dimensional gels for insoluble and large proteins. Electrophoresis 22:1603–1615.

Article information Continued

Figure 1

Protein profiles of high quality longissimus dorsi muscles (HQLD) and low quality longissimus dorsi muscles (LQLD) from Duroc pigs by image analysis. (A) The overall patterns of total protein bands from individuals. All gels were visualized by sliver staining. (B) Two different protein spots were identified by a mass spectrometric analysis. a*: MYBPC2, myosin binding protein C; b#: LDHA, lactate dehydrogenase A.

Figure 2

Ontological classifications of differentially regulated proteins in high quality longissimus dorsi muscles (HQLD) and low quality longissimus dorsi muscles (LQLD) from Duroc pigs. Of the total 24 identified proteins, 10 and 14 proteins were highly expressed in HQLD and LQLD, respectively. (A) The identified proteins were clustered into 7 categories based on information obtained from DAVID gene ontology (GO) database; (B) The identified proteins were clustered into 3 categories based on their molecular function; (C) The identified proteins were clustered into 5 categories based on their biological processes.

Figure 3

Gross changes and cell viability in myogenic differentiation under oxidative stress condition. (A) Morphological changes of myoblast during myogenesis. Arrow head and dashed lines indicating myoblast with mono nucleus and myotube with three or more nuclei, respectively; (B) Cytotoxic effects of H2O2 in C2C12 treated with the various concentrations of H2O2 (12.5 μM to 1 mM). Data represented mean value of three independent experiments; bars, standard deviation. * Significantly different compared to untreated conditions (24 h); #, significantly different compared to untreated conditions (48 h) (n = 3, p<0.0001); (C) Light microscopy-based images of undifferentiated (day 0) proliferating C2C12 myoblasts and differentiating cells at various time points (up to day 7) under normal and oxidative stress conditions (H2O2 [200 μM]); (D) Immunocytochemistry presented the MYH expression levels of undifferentiated (day 0) and differentiated (day 5). GM, growth medium; DM, differentiation medium. (E) Oxidative stress induced down-regulation of Myog gene expression and continuous up-regulation of MyoD gene expression with glyceraldehyde-3-phosphate dehydrogenase (GAPDH) as a control; (F) Oxidative stress induced down-regulation of Myog protein expression and continuous up-regulation of MyoD protein expression with β-actin as a control. Student’s t test was performed to evaluate statistical significance (ns, non-significant; *** p<0.0001; mean±standard error of the mean; n = 3).

Figure 4

Gene and protein expression levels on high quality longissimus dorsi muscles (HQLD) and low quality longissimus dorsi muscles (LQLD) depending on oxidative stress. Myoblast and myotube cDNA and protein were determined by real-time PCR analysis and western blot analysis, respectively. (A) The quantitative differences of lactate dehydrogenase A (LDHA) at the transcriptional level were measured by real-time PCR during myogenesis with β-actin as control; (B) Cytotoxic effects of resveratrol (RSV) in C2C12 treated with various concentrations of RSV (6.25 μM to 100 μM). Data represented mean value of three independent experiments; bars, standard deviation. * Significantly different compared with untreated conditions (24 h); #, significantly different compared to untreated conditions (48 h) (n = 3, p<0.0001); (C) Quantitative differences of LDHA at the transcriptional level were measured by real-time PCR in RSV-treated C2C12 during myogenesis.

Figure 5

Antioxidant effect of resveratrol in lactate dehydrogenase activity and myogenesis. (A) Difference of lactate dehydrogenase (LDH) activity between high quality longissimus dorsi muscles (HQLD) and low quality longissimus dorsi muscles (LQLD) from individual six groups of Duroc pigs. (B) In vitro model of myogenesis under oxidative stress condition and changes of lactate dehydrogenase activity by resveratrol (RSV) treatment. Student’s t test was performed to evaluate statistical significance (*** p<0.0001; mean±standard error of the mean; n = 3); (C) Patterns of myogenic differentiation (MyoD) and myogenin (Myog) expression depending on oxidative stress condition and RSV treatment during myogenesis by immunocytochemistry; (D) Patterns of MyoD and Myog expression (gene and protein levels) depending on oxidative stress condition and RSV treatment during myogenesis by reverse transcription polymerase chain reaction and western blot analysis.

Table 1

Phenotypic record of meat quality traits in Duroc pigs

Traits N Mean Min Max Median Std
Eye muscle area (cm2) 199 28.13 20.99 33.00 28.15 0.218
Age at 90 kg (d) 199 142.97 123 171 142 0.960
Average daily gain (g) 199 643.34 493.30 780.00 641.75 5.806
Backfat thickness (mm) 199 14.78 10.26 24.23 14.48 0.184
Loin depth (mm) 199 55.61 50.40 60.40 55.80 0.189
Weight at end of test day (kg) 199 99.18 75.00 120.00 100.00 0.861
Meat quality characteristics
 pH 24h 199 5.71 5.43 6.03 5.71 0.109
 Brightness 199 54.50 46.71 62.25 54.40 2.981
 Redness 199 16.98 9.75 21.10 17.27 1.984
 Yellowness 199 11.31 5.45 16.29 11.52 2.076
 Water holding capacity 199 73.55 56.99 89.04 72.80 6.874
 Cooking loss 199 14.53 6.65 22.95 14.64 2.998
 Moisture 199 72.76 63.67 79.33 73.00 1.922
 Color, sensory test 199 5.37 2.50 7.50 5.42 0.873
 Flavor, sensory test 199 5.28 1.25 8.00 5.33 1.187
 Tenderness, sensory test 199 5.28 1.25 8.50 5.33 1.450
 Juiciness, sensory test 199 5.10 1.25 9.75 5.25 1.506
 Palatability, sensory test 199 5.13 1.25 8.50 5.25 1.491
 Shear force 199 6.24 2.46 12.94 5.63 2.215
Biochemical measures
 Palmitic acid 199 25.58 22.74 28.08 25.68 1.044
 Oleic acid 199 35.89 10.89 44.91 40.00 10.118
 Stearic acid 199 17.19 9.72 43.71 13.17 9.799
 Linolenic acid 199 8.42 5.31 13.06 8.06 1.743

Table 2

High- or low-meat quality traits of longissimus dorsi muscles in Duroc pigs

Traits N Mean Min Max Median Std
High quality of longissimus dorsi muscles
 pH 24h 20 5.82 5.73 6.03 5.80 0.083
 Water holding capacity 20 73.96 65.60 85.75 73.42 5.699
 Cooking loss 20 14.40 7.30 21.45 13.82 3.457
 Moisture 20 72.95 70.17 75.33 72.67 1.369
 Shear force 20 6.44 3.83 10.55 6.02 2.229
Low quality of longissimus dorsi muscles
 pH 24h 20 5.60 5.55 5.64 5.60 0.031
 Water holding capacity 20 73.76 62.33 87.35 72.60 7.135
 Cooking loss 20 15.39 9.27 22.95 15.41 2.850
 Moisture 20 72.74 68.00 75.00 73.00 1.713
 Shear force 20 6.14 2.46 10.29 5.60 2.036

Table 3

List of differentially expressed genes on high- and low quality of longissimus dorsi muscles in duroc pigs

No UniProt1 UniGene2 Protein identified Gene name pI MW (kDa) Seq. Cov (%) Individual ion score3

HQLD LQLD
Highly expressed genes in HQLD
 1 Q5S1S4 Ssc.10960 Carbonic anhydrase 3 CA3 7.85 29.4 53.46 150.33 82.96
 2 F1RKI3 Ssc:100518898 histidine triad nucleotide-binding protein 1 LOC100518898 6.87 13.7 11.11 3.48 0
 3 Q9TSX9 Ssc.2979 Peroxiredoxin-6 PRDX6 6.01 25 7.59 2.67 0
 4 D3GGC9 Ssc.3835 Actinin-associated LIM protein 3 PDLIM3 8.12 30.5 5.42 5.96 2.24
 5 F1RGK5 Ssc.51787 Tropomyosin alpha-3 chain TPM3 4.75 28.9 5.24 7.13 0
 6 P34930 Ssc.5145 Heat shock 70 kDa protein 1A HSPA1A 5.73 70 3.12 12.17 4.39
 7 Q04967 Ssc.114 Heat shock 70 kDa protein 6 HSPA6 6.06 71.1 2.95 9.69 3.98
 8 F1SEN8 Ssc.97236 LIM domain-binding protein 3 LDB3 7.78 75.8 1.68 2.83 2.13
 9 F1RH20 Ssc.83876 Myosin-binding protein C, fast-type MYBPC2 6.55 127.6 1.58 2.57 0
 10 F1SMN5 Ssc.46794 Filamin-C FLNC 5.96 290.2 0.44 4.59 0
Highly expressed genes in LQLD
 11 P00355 Ssc.16135 Glyceraldehyde-3-phosphate dehydrogenase GAPDH 8.35 35.8 46.85 81.05 447.46
 12 P00339 Ssc.50275 L-lactate dehydrogenase A LDHA 8.07 36.6 39.76 129.85 230.24
 13 F1SKJ8 Ssc.26154 Parvalbumin 1 PVALB1 5.07 12.1 20.91 15.89 22.65
 14 F1SLA0 Ssc.279 ATP synthase subunit beta ATP5B 5.27 56.3 8.71 0 4.89
 15 Q9GJT2 Ssc.217 S-formylglutathione hydrolase ESD 7.02 31.5 7.45 0 2.32
 16 F1RFH9 Ssc.55270 Sarcoplasmic/endoplasmic reticulum calcium ATPase 1 ATP2A1 5.29 109.1 6.45 1.96 21.11
 17 Q06AA5 Ssc.26272 Tetraspanin-9 TSPAN9 7.44 26.8 6.28 0 2.32
 18 F1S6Q7 Ssc:100523423 ATP synthase subunit delta, mitochondrial -like LOC100523423 5.25 17.5 5.36 0 4.15
 19 A1X898 Ssc.97027 Procollagen-proline 2-oxoglutarate-4-dioxygenase P4HA1 6.01 60.9 4.49 0 2.47
 20 F1RFU5 Ssc.3588 Aspartate aminotransferase GOT2 8.73 24.1 3.72 0 4.7
 21 E7EBY5 Ssc.26469 MACRO domain containing protein 1 MACROD1 9.22 35.1 3.69 2.39 4.3
 22 Q2XQV4 Ssc.11147 Aldehyde dehydrogenase, mitochondrial ALDH2 6.87 56.9 3.26 3.59 4.95
 23 I3LL15 Ssc.16302 Uricase UOX 8.32 34.9 2.63 1.78 4.37
 24 F1SLR6 Ssc.5041 Putative ribosomal RNA methyltransferase NOP2 NOP2 9.06 90.2 1.95 0 6.07

pI, isoelectric point of the protein; MW, molecular weight of the protein; Seq. Cov, percentage of sequence coverage; HQLD, high quality of longissimus dorsi muscles; LQLD, low quality of longissimus dorsi muscles.

1

UniProt, Accession number in the UniProt database.

2

UniGene, UniGene number from NCBI (National Center for Biotechnology Information) database.

3

Individual ion score, TurboSEQUEST or gMASCOT score.

Table 4

Primer sequences used to generate templates for reverse transcription polymerase chain reaction

Gene name Symbol GenBank ID Primer sequence (5′→3′) Size (bp)
Myogenic differentiation MyoD NM_010866 F: GAT GGC ATG ATG GAT TAC AGC
R: GAC TAT GTC CTT TCT TTG GGG
528
Myogenin Myog NM_031189 F: GCT CAG CTC CCT CAA CCA G
R: ATG TGA ATG GGG AGT GGG GA
424
Glyceraldehyde-3-phosphate dehydrogenase GAPDH NM_008084 F: ACC ACA GTC CAT GCC ATC AC
R: TAC AGC AAC AGG GTG GTG GA
452