RESEARCH ARTICLE

Comparison of growth performance and related gene expression of muscle and fat from Landrace, Yorkshire, and Duroc and Woori black pigs

Bosung Kim1,#https://orcid.org/0000-0002-5417-0238, Yejin Min2,#https://orcid.org/0000-0002-3083-1513, Yongdae Jeong2https://orcid.org/0000-0002-1985-583X, Sivasubramanian Ramani1https://orcid.org/0000-0001-9370-6552, Hyewon Lim1https://orcid.org/0000-0001-9458-9941, Yeonsu Jo1https://orcid.org/0000-0002-1763-4419, Woosang Kim1https://orcid.org/0000-0003-3742-5603, Yohan Choi2,*https://orcid.org/0000-0003-4710-4731, Sungkwon Park1,*https://orcid.org/0000-0002-7684-9719
Author Information & Copyright
1Department of Food Science and Biotechnology, Sejong University, Seoul 05006, Korea
2Rural Development Administration, National Institute of Animal Science, Cheonan 31000, Korea
*Corresponding author Yohan Choi, Rural Development Administration, National Institute of Animal Science, Cheonan 31000, Korea. Tel: +82-41-580-3451, E-mail: cyh6150@korea.kr
*Corresponding author Sungkwon Park, Department of Food Science and Biotechnology, Sejong University, Seoul 05006, Korea. Tel: +82-2-3408-2906, E-mail: sungkwonpark@sejong.ac.kr

#These authors contributed equally to this work.

© Copyright 2023 Korean Society of Animal Science and Technology. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Jul 19, 2022; Revised: Sep 26, 2022; Accepted: Oct 25, 2022

Published Online: Jan 31, 2023

Abstract

The purpose of this study was to compare marbling score, meat quality, juiciness, sarcomere length, and skeletal muscle satellite cell (SMSC) growth and related gene expression between Woori black pig (WB) and the Landrace, Yorkshire, and Duroc (LYD) crossbreed at different body weights (b.w.). WB was developed to improve meat quality and growth efficiency by crossbreeding Duroc with Korean native black pig. A total of 24 pigs were sacrificed when their b.w. reached about 50, 75, 100, and 120 kg. SMSC were isolated from the femoris muscles, and muscle and adipose tissues were sampled from the middle and the subcutaneous part of the femoris of hind legs, respectively. Expression levels of genes including Myoblast determination protein 1 (MyoD), Paired box gene 3 (Pax3), Myosin heavy chain (MyHC), and Myogenin, which are responsible for the growth and development of SMSC, were higher in LYD than the WB. Muscle growth inhibitor myostatin (MSTN), however, was expressed more in WB compared to LYD (p < 0.01). Numbers of SMSC extracted from femoris muscle of LYD at 50, 75, 100, and 120 kg b.w. were 8.5 ± 0.223, 8.6 ± 0.245, 7.2 ± 0.249, and 10.9 ± 0.795, and those from WB were 6.2 ± 0.32, 6.2 ± 0.374, 5.3 ± 0.423, and 17.1 ± 0.315, respectively. Expression of adipogenic genes in adipose tissue including CCAAT/enhancer-binding protein (CEBP)-β, peroxisome proliferator activated receptor (PPAR)-γ, and fatty acid synthase (FASN), were greater in WB when compared with LYD (p < 0.01). Results from the current study suggest that different muscle cell numbers between 2 different breeds might be affected by related gene expression and this warrants further investigation on other growth factors regulating animal growth and development.

Keywords: Pigs; Genetic analysis; Cell growth; Skeletal muscle; Fat

INTRODUCTION

Global meat production is being increased, estimated to 337.2 million tons in 2020 [1]. Furthermore, following the COVID-19 pandemic, meat production is expected to rise to 373 million tons by 2030 [2]. According to the global trend in 2018, 127.31 million tons of chicken, 4.46 million tons of duck, 120.88 million tons of pork, and 71.61 million tons of beef were consumed among meat groups [3]. In addition, beef and buffalo meat accounted for about 22% of meat production from 1961 to 2018, which was reduced by half, but pork production remained constant at approximately 35%–40% [3]. Landrace is a large white pig crossbred from Denmark. This breed requires a slower feeding cycle but has a faster growth rate [4]. The Yorkshire breed, which originates from England, has good muscles not only in the pork belly, but also in the thigh due to its good growth rate and long body. Hence, many countries have used this type as a crossbreed to produce meat [5]. The Duroc breed originated from the United States shows a greater capacity of accumulating intramuscular fat (IMF), which shows better eating quality, flavor and consumer preference [6,7]. The crossbreed of Landrace, Yorkshire, and Duroc (LYD) is raised as a meat-producing pig because of its excellent fertility and high-quality meat [8]. The Korean traditional pig, a black and short-haired breed, was bred since the 1970s [9,10]. Compared to the LYD crossbreed, it has a lower performance, but a stronger disease prevention capacity [10], but has a better texture and harder fat than the LYD crossbreed [11]. Woori black pig (WB) was developed to improve meat quality and growth efficiency by crossbreeding Duroc with Korean native black pig [12]. Among factors that affect the meat quality including juiciness, color, pH, IMF, and sarcomere length, sarcomere length is positively related to the meat quality or marbling score [13]. Another important factor muscle fiber type, which is categorized into several different types depending on myosin heavy chains (type I, II, IIa, IIx, or IIb), energy metabolism (oxidative or glycolytic), or speed of contraction (fast or slow), affects muscle growth and development as well as meat quality [14]. Skeletal muscle satellite cells (SMSC), also known as muscle stem cells located between the basal lamina of muscle fibers and the fascia, play crucial roles in skeletal muscle hypertrophy and regeneration [15]. Characterization of SMSC, therefore, can indirectly provide the capacity of muscle growth, development, and metabolism ultimately affecting the animal growth efficiency and meat quality [16,17]. There are several myogenic genes including Myosin heavy chain (MyHC), paired box gene 3 (Pax3), myoblast determination protein 1 (MyoD), and myogenin [18]. MyHC defines the fiber type and contraction speed [19]. Myostatin, encoded by the MSTN gene inhibits the growth and differentiation of muscle cells [20]. Pleomorphic adenoma gene 1 (PLAG1) induces cellular growth and in some cases it is related to specific cancers [21]. Adipogenic genes include fatty acid synthase (FASN) a key enzyme related to the synthesis of fatty acids. Stearoyl-CoA desaturase-1 (SCD), which inhibits the growth of adipocytes, and Peroxisome proliferator activated receptor-γ (PPAR-γ), the major regulator involved in differentiation of adipocytes [22]. CCAAT/enhancer-binding protein (CEBP) exerts its functions during the initial stage of adipogenesis, and adiponectin plays a role in body fat reduction [23]. Although growth rate and meat quality of WB has been investigated [24], there has not been many studies which compare and analyze the growth and differentiation of muscle cells and relevant gene expression. The current study therefore compares the number of SMSC and related gene expression patterns of muscle and fat tissues of LYD to those of WB.

MATERIALS AND METHODS

Materials

Dulbecco’s Modified Eagle Medium-F12 (DMEM-F12), fetal bovine serum (FBS), horse serum (HS), and antibiotic-antimycotic (AA) were purchased from Gibco (Thermo Fisher Scientific, Waltham, MA, USA). Trizol reagent (AccuzolTM Total RNA Extraction Reagent, Bioneer, Seoul, Korea), diethylpyrocarbonate (DEPC) water (Bioneer), cDNA transcription kit (AccuPower CycleScript RT PreMix, Bioneer), qPCR MasterMix (Bioneer), and nuclease-free-water (Ambion®, Austin, TX, USA) were purchased for RNA extraction, cDNA synthesis, and quantitative real-time polymerase chain reaction (PCR). The experimental protocols for this research were reviewed and approved by the Institutional Animal Care and Use Committee at the National Institute of Animal Science (NIAS-2020-437).

Porcine muscle resection

The WB and LYD sire pigs were provided by the National Institute of Animal Science, Rural Development Administration (RDA), Korea. Animal sampling methods were followed upon ethical clearance. The pigs were slaughtered for muscle sample collection at RDA.

A total of 24 pigs (12 LYD and 12 WB) were used in this study. Three pigs in each breed were sacrificed when their body weights (b.w.) reached 50 ± 0.20, 75 ± 2.89, 100 ± 1.78, and 120 ± 2.37 kg of LYD, and 50 ± 1.31, 75 ± 1.30, 100 ± 2.20, and 120 ± 0.95 kg of WB. Muscle and adipose tissues were obtained from the hind leg. Skin of hind leg was sterilized with 70% ethanol, subcutaneous fat sampled. Tendons were transected at both proximal and distal sides while holding one tendon with forceps and transecting the other tendon with the scalpel and femoris muscle was harvested.

Culture procedures of muscle cells

SMSC were isolated from the femoris muscle using Pronase enzyme digestion followed by centrifugation at 1,200×g and cells were then seeded at 2×104 cells/mL in a T-25 flask, cultured in growth media consisting of DMEM-F12 containing 1% of AA and 10% of FBS. After 24 hours, the media was replaced with fresh growth media, and the number of cells was counted on day 1. Then, the media was replaced with fresh media every 48 hours. Trypsinization for cell counting was performed using 0.05% of Trypsin-EDTA (Gibco, Gaithersburg, TN, USA) on days 1, 3, 5, and 7. The number of cells were counted using a hemocytometer (Counting Chamber, Paul Marienfeld GmbH & Co., Am Wöllerspfad, Germany). This was performed in triplicate for each sample and averaged each day. Population doubling time (PDT) was calculated by PDT Calculating software.

Muscle and adipose tissue RNA isolation

The total RNA was isolated using Trizol reagent (ACcuzolTM Total RNA Extraction Reagent, Bioneer). Specifically, the Trizol reagent was used at a concentration of 1 mL per 1 g of muscle or adipose tissue sample. Next, 200 µL of chloroform per 1 mL of Trizol reagent was added and the sample was shaken vigorously for 15 seconds. The mixture was then incubated on an ice block rack for 5 minutes and centrifuged at 12,000×g for 15 minutes at 4°C. Then, take off the upper aqueous phase and transferred to a new tube where an equal volume of isopropyl alcohol was added. The sample was then mixed and incubated at −20°C for 10 minutes. Next, the sample was centrifuged at 12,000×g for 10 minutes at 4°C. The supernatant was then removed upon which 1 mL of 80% ethanol was added. The sample was mixed well and then centrifuged at 12,000×g for 5 minutes at 4°C. The supernatant was then removed, and the pellet was dried. Finally, the pellet was resuspended in 20 µL of DEPC water (Bioneer) and the purity and RNA concentration were checked for use in cDNA synthesis using a Microplate Spectrophotometer (Multiskan Sky, Thermo Fisher Scientific) and µDrop plate (µDropTM, Thermo Fisher Scientific). The concentration of the obtained RNA samples was then adjusted to 1 µg. The RNA samples were then reverse transcribed into cDNA using the cDNA reverse transcription kit (AccuPower CycleScript RT PreMix, Bioneer) with a GeneAmp PCR System 9700 machine (Applied Biosystem, Singapore). RNA samples were added to the PreMix tub, which was filled with 20 µL of DEPC water. The machine was then operated based on the cDNA synthesis condition, which was composed of a synthesis step at 45°C for 60 minutes and a heat inactivation step at 95°C for 5 minutes.

Quantitative real-time polymerase chain reaction

The cDNA was used in the comparative cycle threshold experiment for relative gene expression. A quantitative real-time PCR was performed using the AccuPower® 2X Greenstar qPCR MasterMix (Bioneer) and StepOnePlus Real-Time PCR system (Applied Biosystems). First, the cDNA concentration was adjusted to 100 ng/µL with nuclease-free-water (Ambion®). Then, the 20 µL total volume was made by mixing 1 µL of diluted cDNA, 10 µL of the master mix (containing fluorescent material), 0.5 µL of the forward primer, 0.5 µL of the reverse primer, and 8 µL of nuclease-free-water. The mixtures were then stored in reaction tubes (MicroAmp® Fast Reaction Tubes, Applied Biosystems). qPCR reaction conditions were implemented based on the following protocol. During the first stage, the initial denaturation was completed at 95°C for 5 minutes. Denaturation was then performed at 95°C for 15 seconds. The next step involved 40 repeated cycles (the annealing temperatures and times are shown in Table 1), and an extension at 72°C for 45 seconds. The last step was the melt curve, which was included in the final extension and hold phase. The final extension was performed at 95°C for 15 seconds, 60°C for 1 minute, and 95°Cfor 15 seconds. During the last hold stage, the temperature remained at 4°C for infinity. The primer sequences of the housekeeping and target genes are listed in Table 1.

Table 1. The primer sequences used for quantitative real-time PCR of muscle and adipose tissue
Gene Primer sequence Temp (°C) Time (s)
GAPDH F: CTCAACGGGAAGCTCACTGG 55.9 30
R: TGTCGTACGAGGAAATGAGC
Pax 3 F: GCAGCACCGTTCACAGACCT 57.9 30
R: CGGGGTTCATGGGGTTGGAG
MyoD F: TGCGTATTCTCAACCCCTTC 53.8 30
R: AGTATGCAAGGGTGGAGTGG
Myogenin F: GTGAATGAGGCCTTTGAGGC 53.8 30
R: TGTGGGAACTGCATTCACTG
MHC F: TCAAGGGGAGATCACTGTCC 53.8 30
R: TCAGCAACTTCTGTGCCATC
MSTN F: TCTCGATGCTGTCGTTACCCT 59.9 30
R: GCACCAAGCAAACCCCAGA
PLAG1 F: ACCCGTTCGGTTCTACCTCA 53.8 30
R: TTAGACGACGACGCTGGAGA
FASN F: GTCCTGCTGAAGCCTAACTC 53.8 30
R: TCCTTGGAACCGTCTGTG
CEBP-β F: TGTGTACAGATGAATGATAAACTCTGC 55.2 30
R: GGTTTCGAAGTTGATGCAATC
SCD F: CTACACAACCACCACTACCATCAC 57.4 30
R: GCAAACGCCCAGAGCAAGG
PPAR-γ F: ATGGGTGAAACTCTGGGAGA 53.8 30
R: TCAAAGGAGTGGGAGTGGTC
Adiponectin F: TCCACGTCACGGTCTACTTG 53.8 30
R: TTCTCTTCATCCCCGTATGC

PCR, polymerase chain reaction.

Download Excel Table
Determination of sarcomere length and muscle fiber cross section area

Each sample was 1 × 1 × 1 cm cut in the orientation of the muscle fiber on a flat surface and stored at -85°C in a deep freezer (TSE320GPD, Thermo Fisher Scientific). The frozen sample was cut into 10 μm sections at -25°C with a cryostat cryocut micro-tome (CM3050 S, LEICA®, Germany), then the sarcomere length was observed under a high- resolution field emission scanning microscope (MIRA3-LM, Tesan, Czech Republic). The muscle fiber cross section area was measured using a laser confocal scanning microscope (ZEISS LSM 800, ZEISS, Oberkochen, Germany).

Statistical analysis

Statistical analysis was performed using Prism 9.4.0 (GraphPad). The data are presented as mean ± standard deviation. Two-way analyses of variance (ANOVA) followed by Tukey significant difference test were performed. Significant differences were considered by p< 0.05 (*p< 0.05, **p< 0.01, and ***p< 0.001) in all Figures.

RESULTS AND DISCUSSION

Comparison of growth performance

Result of growth performance of the LYD and WB is shown in Table 2. The initial and final b.w. of the LYD crossbreed were 26.20 ± 0.60 kg and 124.37 ± 2.37 kg, those of WB were 24.86 ± 2.40 kg and 120.63 ± 0.95 kg, respectively, with no significant difference. Average daily gain, however, was greater in LYD crossbreed than WB (975.70 ± 105.15 g vs. 768.49 ± 69.17 g; p< 0.001). Daily feed intake was not significantly different between 2 groups (2,736.04 ± 120.35 g vs. 2,736.52 ± 131.72 g; p=0.974).

Table 2. Comparison of growth performance of LYD and the WB
Category LYD WB p-value
Entire period (1–90 days)
Initial weight (average 65 days, kg) 26.20 ± 0.60 24.86 ± 2.40 0.097
Final weight (average 155 days, kg) 124.37 ± 2.37 120.63 ± 0.95 0.126
Daily gain (g) 975.70 ± 105.15 768.49 ± 69.17 0.001
Daily feed intake (g) 2,736.04 ± 120.35 2,736.52 ± 131.72 0.974

LYD, Landrace, Yorkshire, and Duroc; WB, Woori black pig.

Download Excel Table
Cell proliferation and doubling time

Number of cells harvested, growth rate and related gene expression is crucial to understand the mechanisms by which genetic traits regulate the muscle growth and development [16]. Average muscle cell yields were tended to be higher in LYD than WB at 50, 75, or 100 kg b.w. groups, but WB at 120 kg showed more cells than LYD (p < 0.01) (Fig. 1). Although there was a difference in SMSC number between LYD and WB, cell morphology and doubling time was not different (Figs. 2 and 3). Cell doubling time refers to the time taken to double the number of cells, and the difference can stem from various reasons including breeds, age, genetic traits, gender of animal, or cell culture condition. To this end, genes related to myogenesis were analyzed.

jast-65-1-160-g1
Fig. 1. Cell proliferation between LYD and WB. Results were expressed as Mean ± SD, **p < 0.01. LYD, Landrace, Yorkshire, and Duroc; WB, Woori black pig.
Download Original Figure
jast-65-1-160-g2
Fig. 2. Cell doubling time between LYD and WB. Results were expressed as Mean ± SD, and there was no significant difference between both breeds. LYD, Landrace, Yorkshire, and Duroc; WB, Woori black pig.
Download Original Figure
jast-65-1-160-g3
Fig. 3. Cell morphology taken by 0, 1, 3, 5 and 7 days within passage 2, LYD (left), WB (right). LYD, Landrace, Yorkshire, and Duroc; WB, Woori black pig.
Download Original Figure
Gene expression analysis

Cell number is increased through the proliferation process. With external stimuli and growth factors, muscle forms its structure from multiple elongated muscle fibers via cell fusion or differentiation [25]. This developmental process involves many factors and genes [26]. During maturation, orchestration of these gene expressions render the myoblasts fused into myotubes and matured [27]. MyoD acts as a transcriptional activator and thus is engaged in muscle hyperplasia and hypertrophy. Myogenic factor 5 (Myf5) and myogenin (MyoG) are involved in promoting gene expression and muscle development [28]. MyoG accelerates the transcription of muscle-specific genes. The gene expression pattern in skeletal muscle tissue shows that MyoD is expressed throughout all stages [29]. As shown in Fig. 4A, muscle from LYD crossbreed had a greater expression levels of MyoD, MyHC, and Pax3 than the WB. In particular, MyoD expression was higher in LYD than in WB at 50, 75, and 100 kg b.w. groups (p < 0.001).

jast-65-1-160-g4
Fig. 4. Relative gene expression of LYD and WB muscle tissues were analyzed by qPCR. The results were expressed as Mean ± SD (n=3), (*p < 0.05, **p < 0.01, ***p < 0.001). (A) Gene expression of LYD and WB compared with MyoD. (B) Gene expression of LYD and WB compared with MyHC. (C) Gene expression of LYD and WB compared with myogenin. (D) Gene expression of LYD and WB compared with Pax3. (E) Gene expression of LYD and WB compared with MSTN. LYD, Landrace, Yorkshire, and Duroc; WB, Woori black pig; qPCR, quantitative polymerase chain reaction.
Download Original Figure

MyHC is a major structural protein in muscle that converts chemical energy into mechanical energy through the hydrolysis of ATP and expressed in proportion to the amount of muscle [30]. In Fig. 4B, muscle from LYD showed a higher MyHC expression.

Pax3, a paired box transcription factor that regulates proliferation, migration, and cellular apoptosis, plays a key role in controlling myoblast fusion and skeletal muscle fiber development and differentiation [31]. Pax3 expression was higher in femoris muscles from LYD at b.w. of 50, 75, 100, and 120 kg compare to those from WB (Fig. 4C; p < 0.001). Based on its role in muscle development, relatively lower b.w. of WB pigs might be related to this lower Pax3 expression in the current study.

Myogenin is engaged in the functions of MyoD, Myf5, and myogenic regulatory factor 4 (MRF4) as a MyoD family [32]. Myogenin expressed more in muscle tissue from 50 kg LYD (Fig. 4D) compare to other groups (p < 0.001).

MSTN binds various transforming growth factor (TGF)-beta receptors that induce activation of SMAD sequence transcription factors down-regulating skeletal muscle cell proliferation and differentiation [33]. In Fig. 4E, MSTN expression was higher in WB group at 50, 100, 120 kg b.w. groups than in LYD (p < 0.001) which can partially explain the relatively lower daily gain in BW in the current study.

Adipogenesis is a process involving the development preadipocytes derived from multipotent mesenchymal stem cells (MSCs) [34]. It is a multistep process with the sequential activation of numerous transcription factors containing CEBP and PPAR-γ [35]. To reach maturity, these cells must go through three major well-defined steps: commitment of MSCs to the adipocyte lineage; somatic cell division and expansion involving DNA and cell duplication as well as terminal differentiation, including gene expression and transcription factors (such as CEBP and PPAR-γ); and an increase in lipid formation and the introduction of lipogenic genes, including acetyl CoA, carboxylase, FASN, and adipocyte fatty acid binding protein [36]. PLAG1 is an RNA-specific polymerase which regulates the activity of DNA-binding transcriptional factors and proximal promoter DNA-binding transcription, which induces the up-regulation of target genes, a higher proliferation rate, and transformation. Some stimulators include PPAR-γ, macrophage colony stimulating factor, prostaglandins, insulin-like growth factor I (IGF-l), glucocorticoids, and fatty acids [37]. CEBP-β is considered the most important factor and is induced rapidly after the induction of adipogenic stimuli [38].

As shown in Fig. 5A, LYD muscle from 75, 100, 120 kg showed greater CEBP-β expression (p < 0.001).

jast-65-1-160-g5
Fig. 5. Relative gene expression of LYD and WB adipose tissues were analyzed by qPCR. The results were expressed as Mean ± SD (n=3), (*p < 0.05, **p < 0.01, ***p < 0.001) (A)Gene expression of LYD and WB compared with CEBP-β. (B) Gene expression of LYD and WB compared with PPAR-γ. (C) Gene expression of LYD and WB compared with PLAG1. (D) Gene expression of LYD and WB compared with Adiponectin. (E) Gene expression of LYD and WB compared with FASN. FASN, fatty acid synthase; LYD, Landrace, Yorkshire, and Duroc; WB, Woori black pig; qPCR, quantitative polymerase chain reaction.
Download Original Figure

Adipose tissue from LYD at 75, 100, and 120 kg b.w. showed a higher expression than the WB (Fig. 5A; p<0.001). As a master regulator of fat cell differentiation, PPAR-γ is working with CEBP [39]. Gene expression pattern of PPAR-γ and CEBP was similar in Figs. 5A and 5B, but these patterns are not always mirroring the actual adipogenesis. Since their expression is precisely regulated in a time-dependent manner, expression levels of these genes can be up- or down-regulated very rapidly.

PLAG1 (Fig. 5C), known as the zinc finger transcription factor, has an association with muscle growth showing that growth of PLAG1 knock-out animals decreased by 50% compared to the normal group [40]. Expression of PLAG1 in the current study, however, is somewhat inconsistent and tended to decrease as b.w. increased in both breeds. Since gene expression is a rapid and complex process to cope with the environmental cues, integration of gene expression will be necessary to precisely understand this expression plasticity.

Adiponectin expression was greatest in muscle from 50 kg LYD among groups (Fig. 5D; p<0.001). Adiponectin plays a role in improving insulin resistance linked to body fat reduction [41], so, its expression pattern might mirror the speed of muscle growth. Further investigation is needed to elucidate this relationship.

Based on the expression pattern of FASN in Fig. 5E, more fat accumulation could be expected in WB compare to LYD since it is involved in fat synthesis [42].

Determination of sarcomere length and muscle fiber cross section area

Sarcomere is a major player in muscular contraction and contains 28 or more proteins including myosin, actin, titin, tropomyosin, troponin, and nebulin. Among these, actin and myosin play an important role in muscle contraction together with tropomyosin and troponin. In the relaxed muscle of a living animal, the sarcomere length is estimated to be around 2.5 µm. Meanwhile, in the rigor mortis condition where the energy inside the muscles after death is exhausted, actin and myosin filaments cannot be detached from one another, and the corresponding sarcomere length is reduced by half the typical length [43]. Interestingly, a study conducted by Ward et al. indicated that sarcomere length is positively related to the marbling scores [44]. Average sarcomere length of muscles from WB and LYD were 1.710 μm and 1.655 μm, respectively with no statistical difference (p= 0.618) (Fig. 6).

jast-65-1-160-g6
Fig. 6. Sarcomere length between 120 kg of LYD and WB. There was no significant difference between both breeds. p = 0.618. LYD, Landrace, Yorkshire, and Duroc; WB, Woori black pig.
Download Original Figure
Muscle fiber cross-sectional area

Size of skeletal muscle depends on the size and number of the muscle fiber. In the early stage of growth, the size of the muscle is affected by the number of muscle fibers, and the amount of muscle increases with growing muscle fibers afterward [45]. As 75%–90% of the muscle is composed of muscle fibers [46], total number of fibers and cross-sectional areas play a key role in deciding the muscle weight, juiciness and flavor of the meat [47]. Research conducted by Choi and Oh showed that when there exists little or no difference in types I, IIa, and IIb within the muscle fiber, no significant difference was found in meat tenderness, juiciness, and flavor [48]. Concomitant to these results, there was no difference in the cross-sectional fiber areas in muscles from either WB or LYD (Fig. 7) (p<0.663) indicating that meat quality is not supposed to be different between these breeds.

jast-65-1-160-g7
Fig. 7. Muslce fiber cross-sectional area between LYD and WB. There was no significant difference between both breeds. p=0.663. LYD, Landrace, Yorkshire, and Duroc; WB, Woori black pig.
Download Original Figure

CONCLUSION

WB is a crossbreed for better meat quality as well as growth efficiency. When comparing growth performance between the two breeds, the final weight of the LYD and WB was 124.37 ± 2.37 kg and 120.63 ± 0.95 kg, respectively. Results of SMSC analysis confirmed that the LYD had more number of SMSC, and PDT was faster than WB. Moreover, as no significant difference was found in sarcomere length, muscle fiber number, and cross-sectional areas, a similar muscle growth rate in the LYD and WB is expected. Gene expression results showed that reduction in myogenesis in SMSC of WB could be induced by down regulation of Pax3 via MSTN (Fig. 8). Expression of adipogenic genes including CEBP-β and FASN, which affect the growth, development, and accumulation of fat, was higher in WB than in LYD. Interestingly, PPAR-γ, which regulates the differentiation of adipose tissue, was expressed lower in WB compare to LYD. Altogether, relatively lower daily weight gain might stem from myostatin activation in WB and results from our current study warrant further investigation of comparative analysis for different pig breeds.

jast-65-1-160-g8
Fig. 8. Myostatin (MSTN) inhibition on Pax3, myosin heavy chain (MyHC), Myogenin, and MyoD.
Download Original Figure

Competing interests

No potential conflict of interest relevant to this article was reported.

Funding sources

This work was carried out with the support of “Cooperative Research Program for Agriculture Science and Technology Development (Project No. PJ01491801)” Rural Development Administration, Republic of Korea and Technology Innovation Program (20012411, Alchemist Project) funded by the Ministry of Trade, Industry and Energy (MOTIE).

Acknowledgements

This study was supported by 2023 the RDA Fellowship Program of National Institute of Animal Science, Rural Development Administration, Korea.

Availability of data and material

Upon reasonable request, the datasets of this study can be available from the corresponding author.

Authors’ contributions

Conceptualization: Kim B, Choi Y, Park S.

Formal analysis: Kim B, Ramani S.

Methodology: Jeong Y, Ramani S.

Validation: Kim B, Min Y, Park S.

Investigation: Lim H, Jo Y, Kim W, Choi Y.

Writing - original draft: Kim B, Park S.

Writing - review & editing: Kim B, Min Y,

Jeong Y, Ramani S, Lim H, Jo Y, Kim W,

Choi Y, Park S.

Ethics approval and consent to participate

The experimental protocols for this research were reviewed and approved by the Institutional Animal Care and Use Committee at the National Institute of Animal Science (NIAS-2020-437).

REFERENCES

1.

FAO [Food and Agricultural Organization of the United nations]. Meat market review [Internet]. FAO of the United Nations. 2019 [[cited 2022 June 15]]. http://www.fao.org/3/ca3880en/ca3880en.pdf.

2.

OECD [Organisation for Economic Co-operation Development] / FAO [Food and Agricultural Organization of the United nations]. Meat. OECD-FAO Agric outlook 2021-2030. Paris: OECD. 2021; p. 163-77

3.

Ritchie H, Rosado P, Roser M. Meat and dairy production [Internet]. Our World Data. 2017 [[cited 2022 June 15]]. https://ourworldindata.org/meat-production.

4.

Taylor G, Roese G, Hermesch S. Breeds of pigs — Landrace. Primefact. 2005; 63:6148

5.

Hu H, Wu C, Ding Y, Zhang X, Yang M, Wen A, et al. Comparative analysis of meat sensory quality, antioxidant status, growth hormone and orexin between Anqingliubai and Yorkshire pigs. J Appl Anim Res. 2019; 47:357-61

6.

Kim JA, Cho ES, Jeong YD, Choi YH, Kim YS, Choi JW, et al. The effects of breed and gender on meat quality of Duroc, Pietrain, and their crossbred. J Anim Sci Technol. 2020; 62:409-19

7.

Smith WC, Pearson G, Purchas RW. A comparison of the Duroc, Hampshire, Landrace, and Large White as terminal sire breeds of crossbred pigs slaughtered at 85 kg liveweight: 1. performance and carcass characteristics. New Zeal J Agric Res. 1990; 33:89-96

8.

Choi JS, Lee HJ, Jin SK, Choi YI, Lee JJ. Comparison of carcass characteristics and meat quality between duroc and crossbred pigs. Korean J Food Sci Anim Resour. 2014; 34:238-44

9.

Chung HY, Ko MS. Characteristics and application of native pigs in Jeju, South Korea. In. Proceedings of the 29th Society for Jeju Stuies National Conference. 2006; Jeju:95-119

10.

Jin SK, Kim CW, Song YM, Jang WH, Kim YB, Yeo JS, et al. Physicochemical characteristics of longissimus muscle between the Korean native pig and Landrace. Korean J Food Sci Anim Resour. 2001; 21:142-8

11.

Choi YS, Park BY, Lee JM, Lee SK. Comparison of carcass and meat quality characteristics between Korean native black pigs and commercial crossbred pigs. Korean J Food Sci Anim Resour. 2005; 25:322-7

12.

Alam M, Chang HK, Lee SS, Choi T. Genetic Analysis of major production and reproduction traits of Korean Duroc, Landrace and Yorkshire pigs. Animals. 2021; 11:1321

13.

Kemp CM, Sensky PL, Bardsley RG, Buttery PJ, Parr T. Tenderness – an enzymatic view. Meat Sci. 2010; 84:248-56

14.

Maltin CA, Warkup CC, Matthews KR, Grant CM, Porter AD, Delday MI. Pig muscle fibre characteristics as a source of variation in eating quality. Meat Sci. 1997; 47:237-48

15.

Park S, Jeong Y. Characterization of muscle satellite cells in vitro. J Agric Life Sci. 2014; 48:59-67

16.

Li BJ, Li PH, Huang RH, Sun WX, Wang H, Li QF, et al. Isolation, culture and identification of porcine skeletal muscle satellite cells. Asian-Australas J Anim Sci. 2015; 28:1171-7

17.

Rhoads RP, Fernyhough ME, Liu X, McFarland DC, Velleman SG, Hausman GJ, et al. Extrinsic regulation of domestic animal-derived myogenic satellite cells II. Domest Anim Endocrinol. 2009; 36:111-26

18.

Calhabeu F, Hayashi S, Morgan JE, Relaix F, Zammit PS. Alveolar rhabdomyosarcoma-associated proteins Pax3/FOXO1A and PAX7/FOXO1A suppress the transcriptional activity of MyoD-target genes in muscle stem cells. Oncogene. 2013; 32:651-62

19.

LaFramboise WA, Guthrie RD, Scalise D, Elborne V, Bombach KL, Armanious CS, et al. Effect of muscle origin and phenotype on satellite cell muscle-specific gene expression. J Mol Cell Cardiol. 2003; 35:1307-18

20.

McCroskery S, Thomas M, Maxwell L, Sharma M, Kambadur R. Myostatin negatively regulates satellite cell activation and self-renewal. J Cell Biol. 2003; 162:1135-47

21.

Wang J, Huang Y, Xu J, Yue B, Wen Y, Wang X, et al. Pleomorphic adenoma gene 1 (PLAG1) promotes proliferation and inhibits apoptosis of bovine primary myoblasts through the PI3K-Akt signaling pathway. J Anim Sci. 2022; 100:skac098

22.

Wang L, Xue K, Wang Y, Niu L, Li L, Zhong T, et al. Molecular and functional characterization of the adiponectin (AdipoQ) gene in goat skeletal muscle satellite cells. Asian-Australas J Anim Sci. 2018; 31:1088-97

23.

Kokta TA, Dodson MV, Gertler A, Hill RA. Intercellular signaling between adipose tissue and muscle tissue. Domest Anim Endocrinol. 2004; 27:303-31

24.

Kim YM, Choi TJ, Cho KH, Cho ES, Lee JJ, Chung HJ, et al. Effects of sex and breed on meat quality and sensory properties in three-way crossbred pigs sired by duroc or by a synthetic breed based on a Korean native breed. Korean J Food Sci Anim Resour. 2018; 38:544-53

25.

Yin H, Price F, Rudnicki MA. Satellite cells and the muscle stem cell niche. Physiol Rev. 2013; 93:23-67

26.

Lehka L, Topolewska M, Wojton D, Karatsai O, Alvarez-Suarez P, Pomorski P, et al. Formation of aberrant myotubes by myoblasts lacking myosin VI is associated with alterations in the cytoskeleton organization, myoblast adhesion and fusion. Cells. 2020; 9:1673

27.

Hernández-Hernández O, Ávila-Avilés RD, Hernández-Hernández JM. Chromatin landscape during skeletal muscle differentiation. Front Genet. 2020; 11:578712

28.

Zhu KC, Liu BS, Guo HY, Zhang N, Guo L, Jiang SG, et al. Functional analysis of two MyoDs revealed their role in the activation of myomixer expression in yellowfin seabream (Acanthopagrus latus) (Hottuyn, 1782). Int J Biol Macromol. 2020; 156:1081-90

29.

Hitachi K, Nakatani M, Takasaki A, Ouchi Y, Uezumi A, Ageta H, et al. Myogenin promoter-associated lncRNA Myoparr is essential for myogenic differentiation. EMBO Rep. 2019; 20e47468

30.

Wang L, Liu X, Niu F, Wang H, He H, Gu Y. Single nucleotide polymorphisms, haplotypes and combined genotypes in MYH3 gene and their associations with growth and carcass traits in Qinchuan cattle. Mol Biol Rep. 2013; 40:417-26

31.

Wang Q, Fang WH, Krupinski J, Kumar S, Slevin M, Kumar P. Pax genes in embryogenesis and oncogenesis. J Cell Mol Med. 2008; 12:2281-94

32.

Wilschut KJ, Jaksani S, Van Den Dolder J, Haagsman HP, Roelen BAJ. Isolation and characterization of porcine adult muscle-derived progenitor cells. J Cell Biochem. 2008; 105:1228-39

33.

Fennen M, Pap T, Dankbar B. Smad-dependent mechanisms of inflammatory bone destruction. Arthritis Res Ther. 2016; 18:279

34.

Bahmad HF, Daouk R, Azar J, Sapudom J, Teo JCM, Abou-Kheir W, et al. Modeling adipogenesis: current and future perspective. Cells. 2020; 9:2326

35.

Rosen ED, Hsu CH, Wang X, Sakai S, Freeman MW, Gonzalez FJ, et al. C/EBPα induces adipogenesis through PPARγ : a unified pathway. Genes Dev. 2002; 16:22-6

36.

Zhang K, Yang X, Zhao Q, Li Z, Fu F, Zhang H, et al. Molecular mechanism of stem cell differentiation into adipocytes and adipocyte differentiation of malignant tumor. Stem Cells Int. 2020; 2020:8892300

37.

Kang JJ, Auble DT, Ranish JA, Hahn S. Analysis of the yeast transcription factor TFIIA: distinct functional regions and a polymerase II-specific role in basal and activated transcription. Mol Cell Biol. 1995; 15:1234-43

38.

Reusch JEB, Colton LA, Klemm DJ. CREB activation induces adipogenesis in 3T3-L1 cells. Mol Cell Biol. 2000; 20:1008-20

39.

Wang QA, Zhang F, Jiang L, Ye R, An Y, Shao M, et al. Peroxisome proliferator-activated receptor γ and its role in adipocyte homeostasis and thiazolidinedione-mediated insulin sensitization. Mol Cell Biol. 2018; 38:e00677-17

40.

Juma AR, Damdimopoulou PE, Grommen SVH, Van de Ven WJM, De Groef B. Emerging role of PLAG1 as a regulator of growth and reproduction. J Endocrinol. 2016; 228:R45-56

41.

Achari AE, Jain SK. Adiponectin, a therapeutic target for obesity, diabetes, and endothelial dysfunction. Int J Mol Sci. 2017; 18:1321

42.

Li WZ, Zhao SM, Huang Y, Yang MH, Pan HB, Zhang X, et al. Expression of lipogenic genes during porcine intramuscular preadipocyte differentiation. Res Vet Sci. 2012; 93:1190-4

43.

Ertbjerg P, Puolanne E. Muscle structure, sarcomere length and influences on meat quality: a review. Meat Sci. 2017; 132:139-52

44.

Ward SR, Tomiya A, Regev GJ, Thacker BE, Benzl RC, Kim CW, et al. Passive mechanical properties of the lumbar multifidus muscle support its role as a stabilizer. J Biomech. 2009; 42:1384-9

45.

Rehfeldt C, Fiedler I, Dietl G, Ender K. Myogenesis and postnatal skeletal muscle cell growth as influenced by selection. Livest Prod Sci. 2000; 66:177-88

46.

Ruusunen M, Puolanne E. Comparison of histochemical properties of different pig breeds. Meat Sci. 1997; 45:119-25

47.

Joo ST, Kim GD, Hwang YH, Ryu YC. Control of fresh meat quality through manipulation of muscle fiber characteristics. Meat Sci. 2013; 95:828-36

48.

Choi YM, Oh HK. Carcass performance, muscle fiber, meat quality, and sensory quality characteristics of crossbred pigs with different live weights. Korean J Food Sci Anim Resour. 2016; 36:389-96