Genetic Effects of Polymorphisms in Myogenic Regulatory Factors on Chicken Muscle Fiber Traits

Article information

Asian-Australas J Anim Sci.. 2015;28(6):782-787
Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu Campus, Chengdu, Sichuan 611130, China
1College of Animal Science and Technology, Sichuan Agricultural University, Ya’an, Sichuan, 625014, China.
*Corresponding Author: Hua-Dong Yin. Tel: +86-835-2882006, Fax: +86-835-2886080, E-mail: yinhuadong@aliyun.com
Received 2014 September 25; Revised 2014 December 22; Accepted 2015 January 12.

Abstract

The myogenic regulatory factors is a family of transcription factors that play a key role in the development of skeletal muscle fibers, which are the main factors to affect the meat taste and texture. In the present study, we performed candidate gene analysis to identify single-nucleotide polymorphisms in the MyoD, Myf5, MyoG, and Mrf4 genes using polymerase chain reaction-single strand conformation polymorphism in 360 Erlang Mountain Chickens from three different housing systems (cage, pen, and free-range). The general linear model procedure was used to estimate the statistical significance of association between combined genotypes and muscle fiber traits of chickens. Two polymorphisms (g.39928301T>G and g.11579368C>T) were detected in the Mrf4 and MyoD gene, respectively. The diameters of thigh and pectoralis muscle fibers were higher in birds with the combined genotypes of GG-TT and TT-CT (p<0.05). Moreover, the interaction between housing system and combined genotypes has no significant effect on the traits of muscle fiber (p>0.05). Our findings suggest that the combined genotypes of TT-CT and GG-TT might be advantageous for muscle fiber traits, and could be the potential genetic markers for breeding program in Erlang Mountain Chickens.

INTRODUCTION

Because of relatively lower content of fat and cholesterol, and higher level of long-chain polyunsaturated fatty acids, many nutritionists believe that the white meat like poultry, seafood and fish is judged healthier than red meat, such as pork, beef and mutton. Combined with the lower price and rare religious restrictions, the consumption demand of chicken is increasing gradually (Yin et al., 2013). Nowadays, with the improvement of life level and the development of social economy, people begin to lay stress on food quality (Wang et al., 2009). The requirement of free-range chicken is much popular, because it is considered to be more healthier and more suitable for traditional eating habits that require the meat has higher chewiness, especially in southwest China (Dou, 2009). Previous studies have shown that the muscle fiber traits are the main factors to affect the chicken taste and texture (Branciari et al., 2009). The muscle fiber traits are regulated by several genes, the traditional phenotypic selection is very difficult to achieve the genetic improvement in a short time. However, the candidate gene approach is a powerful method to find the quantitative trait loci (QTL) for improving the rate of genetic progress in the meat quality traits (Yin et al., 2011).

In 1987, MyoD was first discovered by subtractive hybridization using myoblast cDNA libraries, which is a basic helix-loop-helix transcription factor that has ability to transform non-muscle cells, such as fibroblasts, into cells that are able to fuse into myotubes (Davis et al., 1987). Subsequently, three related skeletal muscle specific transcriptional factors—Myf5, myogenin (MyoG), and Mrf4—which are also capable of inducing myoblast features in non-muscle cell lines, were identified (Braun et al., 1989; Edmondson and Olson, 1989; Braun et al., 1990). Collectively, MyoD, Myf5, MyoG, and Mrf4 are specifically expressed in skeletal muscle linage and are therefore referred to as myogenic regulatory factors (MRFs) (Weintraub et al., 1991). Proteins of MRFs share a highly conserved basic helix-loop-helix domain that mediates DNA binding. Helix-loop-helix motif is essential for heterodimerization with a ubiquitously expressed E-protein (Massari and Murre, 2000). In the last few decades, studies have made significant advances in distinguishing the function of individual MRFs in myogenesis. It is shown that MyoD and Myf5 are factors involved in determination of myogenic cells, whereas MyoG and Mrf4 are more likely associated with terminal differentiation and homeostasis of myofibers (Singh and Dilworth, 2013).

Due to the very important roles of MRFs in control of the myogenic process, the genetic difference of the MRFs could affect the muscle fibers. To verify this hypothesis, we investigated the genetic polymorphisms of the MyoD, Myf5, MyoG, and Mrf4 genes in chickens and analyzed their association with muscle fibers traits.

MATERIALS AND METHODS

Chicken populations

The line SD02×SD03 of Erlang mountainous chicken was employed in the current study, which is a breeding population by Sichuan Agricultural University. From 1 to 42 d, 2,000 chicks with sexes intermingled were reared in concrete floor pens with shavings as litter. At d43, 1,200 chickens with similar body weight (314.6±24.1 g) were then randomly assigned to 3 different housing systems (400 birds/per group), including conventional cages that contained plastic hollowed pads, indoor floor pens with paddock, and outdoor free-range. Ten chickens were reared in each of 20 cages (1.9×1.0×0.7 m; 5.3 birds/m2). Indoor concrete floor pens were covered with wood shavings (3×5×2.7 m, 3.3 birds/m2), and chickens in this housing system had access to the paddock. The area ratio of indoor to paddock was 1:1. Chickens in the free-range treatment were raised under pines (6,670 m2) with a simple house (6.7 birds/m2) for shelter, and chickens were herded to pine grove for exercise in the daytime. Chickens in different systems were fed the same diet. Feed and tap water were freely available, and all diets were formulated to the feeding standards for yellow-feather chickens in China (NY/T33, 2004).

Sample collection and trait measurement

Blood samples were collected from the brachial vein in another 360 chickens (120 each randomly selected from three housing systems at d84, and the ratio of male and female was 1:1). Genomic DNA was extracted by the standard phenol-chloroform method.

The diameter and density of the muscle fibers were measured on the left pectoralis major muscle and thigh muscle of all selected chickens. The muscle slices (3 cm×2 cm) were cut to follow the direction of muscle fibers in the center of the muscle, then fixed overnight with 10% formalin. Fixed tissues embedded in paraffin were sectioned at 6 μm mounted on slides and stained with hematoxylin-eosin. Fiber diameter measurements were performed on each section under light-microscopy (40×) and the number of myofibers per square millimeter was computed.

Single-nucleotide polymorphism identification

Genetic variants in the MRFs genomic DNA sequences were analyzed using the polymerase chain reaction (PCR)-single strand conformation polymorphism method (SSCP). Twenty primer pairs were designed according to the Gallus gallus MRF genes sequences (Table 1) to amplify the eleven exons fully containing the coding sequences of the four genes. The PCR reactions included 0.8 μL (100 ng/μL) of template DNA, 5 μL of 2×Taq PCR MasterMix (TianGen Biochemical Technology, Shanghai, China), 0.4 μL of each primer (10 pmol/μL), and ddH2O to 10 μL reaction volume using the following condition: 95°C for 10 min, 35 cycles of 94°C for 30 s, annealing from 57°C to 62°C for 45 s, 72°C for 1 min, and a final extension at 72°C for 10 min. After PCR, single-nucleotide polymorphism (SNP) analysis was performed according to Orita et al. (1989). Six DNA samples showing identical patterns on the SSCP gel were further amplified in a 50 μL reaction volume, electrophoresed and purified with a gel extraction kit (Takara, Dalian, China), and sequenced by Shanghai Yingjun Biology Technique Corporation. The primer pairs used to amplify the fragments were also used for sequencing.

Primers for SNP detection

Data analysis

Combined genotypes were constructed based on the detected SNPs in all 360 experimental chickens. Genetic effect of the combined genotype was evaluated with the following model:

Yijk=μ+Ci+Hj+Sk+(C×H)ij+(H×S)jk+(C×S)ik+(C×H×S)ijk+eijk,

where Yijk is the trait measured in chicken; μ is the population means of that trait; Ci is the fixed effect associated with the combined genotypes i; Hj is the fixed effect associated with the housing system j, Sk is the fixed effect of sex k; eijk is the residual random effect, and the C×H, H×S, C×S, and C×H×S is the interaction effect of combined genotype i by housing system j, housing system j by sex k, combined genotype i by sex k, and combined genotype i by housing system j by sex k, respectively. Multiple comparisons were carried out using Tukey’s test. Pearson’s chi-square test was used to check for Hardy-Weinberg equilibrium of the SNP detected in the sample population. Polymorphism information content (PIC) was evaluated by Botstein’s method (Orita et al., 1989).

Data were presented as least square means±standard error of the mean, and differences were considered significant at p<0.05.

RESULTS

Single-nucleotide polymorphism identification and genetic diversity analysis

Twenty expected fragments (five fragments per gene) covering the full length of the coding sequences of Myf5, MyoD, Mrf4, and MyoG were amplified by the primers. Among them, two polymorphisms were detected by directly sequencing the polymorphic fragment based on the SSCP banding patterns, one of which was g.39928301T>G in the Mrf4 gene (Figure 1), and another is the mutation of g.11579368C>T in the MyoD gene (Figure 2). The variant of g.39928301T>G was a non-synonymous mutation and caused a GTG (Histidine) to GGG (Proline) change, while g.11579368C>T was a synonymous mutation that was not associated with an amino acid substitution.

Figure 1

The partial sequencing result of Mrf4 gene, and g. 39928301T>G was marked by rectangle.

Figure 2

The partial sequencing result of MyoD gene, and g.11579368C>T was marked by rectangle.

Two SNPs, their corresponding allele and genotype frequencies, PIC, and p-values for the Hardy-Weinberg equilibrium test are shown in Table 2. Three genotypes each from two SNPs were identified. Genotypes TT, TG, and GG of g.39928301T>G had an average frequency of 42.22%, 35.56%, and 22.22%; while genotypes CC, CT, and TT of g.11579368C>T with frequencies of 46.11%, 32.78%, and 21.11%, respectively. The allele T was the dominant allele due to the higher frequency (60.00%) compared to the allele G (40.00%) in the SNP of g.39928301T>G. For SNP-g.11579368C>T the allele C represented the higher allele frequency (62.50%). Both g.39928301T>G and g.11579368C>T were out of Hardy-Weinberg equilibrium (p<0.05), and the PIC of them was 0.3648, and 0.3593, respectively.

Frequency of alleles and genotypes, PIC, and the Hardy-Weinberg (H-W) equilibrium test at polymorphic sites

Effect of combined genotypes on the muscle fiber traits

We constructed the combined genotypes from each genotypes of g.39928301T>G and g.11579368C>T (Table 3). In the total 360 samples, there are 9 combined genotypes, those combined genotypes with low frequencies (TT-TT [n = 13, 3.61%], TG-TT [n = 13, 3.61%], GG-CC [n = 13, 3.61%], and GG-CT [n = 17, 4.72%]) were not included for further analysis. Association of the combined genotypes of Myf4 and MyoD with muscle fiber traits is shown in Table 3. The diameters of thigh and pectoralis muscle fibers were higher in birds with the combined genotypes of GG-TT and TT-CT. The density of pectoralis muscle fibers in individuals with GG-TT was lower than TG-CC birds, but the five evaluated combined genotypes had similar densities of thigh muscle fibers. Additionally, there was no interaction of housing system and combined genotype on the traits of muscle fiber (p>0.05).

Association of the combined genotypes with muscle fiber traits in chickens

DISCUSSION

The diameter and density of muscle fiber directly determines the taste of meat, and thereby affects market acceptance. Although more chewy taste in chicken can be achieved through changing the housing system, the feed conversion ratio and growth performances are concurrently decrease (Massari and Murre, 2000; Fanatico et al., 2008; Dou, 2009). Fortunately, selecting individuals with additional information on their genotype for makers associated with QTLs and target traits can effectively balance muscle fiber trait and growth performance in different housing systems. Therefore, we investigated the polymorphisms in the coding regions of Myf5, Mrf4, MyoD, and MyoG, and explored their association with muscle fiber traits. In the present study, two SNPs were detected, locating in exon 3 of the Mrf4 gene, and in exon 1 of the MyoD gene. The PIC used to assess the genetic variation, is divided into three classes, including slight variation (PIC< 0.25), medium polymorphic (0.25<PIC<0.5), and highly informative (PIC>0.5) (Yin et al., 2013). The PIC of 0.3648 and 0.3593 in two SNPs indicated that this population has higher selective space for breeding.

Although the function of synonymous mutations remain unknown, recent studies revealed that silent SNPs can lead to the synthesis of protein products with the same amino acid sequence but different structural and functional properties (Komar, 2007), and alter the structure of substrate and inhibitor interaction sites (Kimchi-Sarfaty et al., 2007). Therefore, we presumed that the synonymous mutation of g.11579368C>T in the MyoD gene changes the style of protein folding or it may link to another mutation of g.39928301T>G in the coding regions of the Mrf4 gene that influence the muscle fiber traits. In the present study, the results showed that the combined genotypes of mutations in MyoD and Mrf4 are associated with variations in the muscle fiber traits. The combined genotypes of TT-CT and GG-TT with greater pectoralis and thigh muscle fiber diameters, and could be selected in a breeding program to improvement of the chicken muscle fiber characteristics.

Environmental sensitivity is the ability of a genotype to alter its phenotypic expression in response to environmental influences. The genetic variation in environmental sensitivity is used as a definition of genotype by environment interaction (Kolmodin et al., 2003). There is a reranking of individuals for traits expression because of the varying magnitude of phenotypic change in the different environments. Consequently, the selection for target performance in one environment may not result in enhanced performance under other environmental conditions (Sing, 1990). In the present study, we screened the polymorphisms in the coding sequences of MRFs of chickens from three housing environments, and analyzed their relationship with muscle fiber traits. Fortunately, we found that there was no interaction between housing system and combined genotypes, which showed that the advantageous combined genotypes of TT-CT and GG-TT may promote the diameter of muscle fiber under different rearing conditions.

CONCLUSION

In conclusion, the present study showed that the combined genotypes of TT-CT and GG-TT may serve as molecular makers for breeding program. Although we can use marker-assisted selection to increase muscle fiber diameter, a greater diameter does not necessarily equate to superior meat quality. Thus, we should find the balance between the diameter of the muscle fiber and meat quality in the further studies.

ACKNOWLEDGMENTS

This work was supported by the China Agriculture Research System (CARS-41), the Twelfth Five Year Plan for breeding programs in Sichuan—Selective breeding of new breeds and the synthetic strains in laying hens (2011NZ0099-7).

References

Branciari R, Mugnai C, Mammoli R, Miraglia D, Ranucci D, Dal Bosco A, Castellini C. 2009;Effect of genotype and rearing system on chicken behavior and muscle fiber characteristics. J Anim Sci 87:4109–4117.
Braun T, Bober E, Winter B, Rosenthal N, Arnold HH. 1990;Myf-6, a new member of the human gene family of myogenic determination factors: Evidence for a gene cluster on chromosome 12. EMBO J 9:821–831.
Braun T, Buschhausen-Denker G, Bober E, Tannich E, Arnold HH. 1989;A novel human muscle factor related to but distinct from myod1 induces myogenic conversion in 10t1/2 fibroblasts. EMBO J 8:701–709.
Davis RL, Weintraub H, Lassar AB. 1987;Expression of a single transfected cdna converts fibroblasts to myoblasts. Cell 51:987–1000.
Dou TC, Shi S, Sun R, Wang HJ, Sun HJ, Wang KH. 2009;Growth rate, carcass traits and meat quality of slow-growing chicken grown according to three raising systems. Anim Sci Papers Rep 27:361–369.
Edmondson DG, Olson EN. 1989;A gene with homology to the myc similarity region of myod1 is expressed during myogenesis and is sufficient to activate the muscle differentiation program. Genes Dev 3:628–640.
Fanatico AC, Pillai PB, Hester PY, Falcone C, Mench JA, Owens CM, Emmert JL. 2008;Performance, livability, and carcass yield of slow- and fast-growing chicken genotypes fed low-nutrient or standard diets and raised indoors or with outdoor access. Poult Sci 87:1012–1021.
Kimchi-Sarfaty C, Oh JM, Kim I-W, Sauna ZE, Calcagno AM, Ambudkar SV, Gottesman MM. 2007;A “silent” polymorphism in the mdr1 gene changes substrate specificity. Science 315(5811):525–528.
Kolmodin R, Strandberg E, Jorjani H, Danell B. 2003;Selection in the presence of a genotype by environment interaction: Response in environmental sensitivity. Anim Sci 76:375–386.
Komar AA. 2007;Silent snps: Impact on gene function and phenotype. Pharmacogenomics 8:1075–1080.
Massari ME, Murre C. 2000;Helix-loop-helix proteins: Regulators of transcription in eucaryotic organisms. Mol Cell Biol 20:429–440.
Orita M, Suzuki Y, Sekiya T, Hayashi K. 1989;Rapid and sensitive detection of point mutations and DNA polymorphisms using the polymerase chain reaction. Genomics 5:874–879.
Sing CF. 1990;Introduction to quantitative genetics. Am J Hum Gent 46(6):1231.
Singh K, Dilworth FJ. 2013;Differential modulation of cell cycle progression distinguishes members of the myogenic regulatory factor family of transcription factors. FEBS J 280:3991–4003.
Wang KH, Shi SR, Dou TC, Sun HJ. 2009;Effect of a free-range raising system on growth performance, carcass yield, and meat quality of slow-growing chicken. Poult Sci 88:2219–2223.
Weintraub H, Davis R, Tapscott S, Thayer M, Krause M, Benezra R, Blackwell TK, Turner D, Rupp R, Hollenberg S. 1991;The myod gene family: Nodal point during specification of the muscle cell lineage. Science 251(4995):761–766.
Yin HD, Wang Y, Zhang ZC, Liu YP, Chen SY, Zhu Q. 2013;Characterization of the expression profile and genetic polymorphism of the cellular retinol-binding protein (crbp iv) gene in erlang mountainous chickens. Int J Mol Sci 14:4432–4443.
Yin HD, Gilbert ER, Chen SY, Wang Y, Zhang ZC, Zhao XL, Zhang Y, Zhu Q. 2013;Effect of hybridization on carcass traits and meat quality of erlang mountainous chickens. Asian Australas J Anim Sci 26:1504–1510.
Yin H, Zhang Z, Lan X, Zhao X, Wang Y, Zhu Q. 2011;Association of myf5, myf6 and myoG gene polymorphisms with carcass traits in chinese meat type quality chicken populations. J Anim Vet Adv 10:704–708.

Article information Continued

Figure 1

The partial sequencing result of Mrf4 gene, and g. 39928301T>G was marked by rectangle.

Figure 2

The partial sequencing result of MyoD gene, and g.11579368C>T was marked by rectangle.

Table 1

Primers for SNP detection

Gene Primer Primer sequence (5′ → 3′) Annealing temperature (°C) Product length (bp) Amplified region a Location of exonb
Myf5 P1-F ACTTTGTCTCGGGTCGCTAAT 57.4 240 g.39931832-39932071 exon 1
P1-R TCAGGGCAGGCAGGCTCA (g.39931917-39932423)
P2-F CCCACTGAGCCTGCCTGC 57.2 294 g.39932049-39932342 exon 1
P2-R TTCCCTCAAGAGCTCCTGG (g.39931917-39932423)
P3-F GGAACGCCATCAGGTACATCG 56.9 141 g.39932296-39932436 exon 1
P3-R GCCCTGCTCTTACCATCACATC (g.39931917-39932423)
P4-F CCCCTGCTCGCCATTTGT 61.0 240 g.39932612-39932851 exon 2
P4-R CCTTACCGTGGGGCATCT (g.39932770-39932845)
P5-F TAAGCGGCTCGCCCCGGTGC 61.7 221 g.39932915-39933135 exon 3
P5-R CCCCTCATAGCGCCTGGT (g.39932938-3993131)
Mrf4 P1-F TCTGGGTTGCTCCTCGGGGTGT 62.2 253 g.39927179-39927431 exon 1
P1-R TGCTCCTCGCCGCTGCTGTC (g.39927250-39927768)
P2-F CCCCGTGTCAGGACCAAC 59.5 240 g.39927374-39927613 exon 1
P2-R CCACAGTCCGCCTTTTTCAG (g.39927250-39927768)
P3-F TCTCAAAAGGCGGACTGTGG 58.4 253 g.39927594-39927846 exon 1
P3-R GGAGGAGGAGGGGAAGGAG (g.39927250-39927768)
P4-F CGCTGCTCGCCTCGCTAA 63.0 286 g.39927933-39928218 exon 2
P4-R GGCCGACGACTCCACCTA (g.39927984-39928074)
P5-F CCGCTCCCTAACCCCTCT 60.1 260 g.39928161-39928420 exon 3
P5-R AAAACGGGAGATGCGGAG (g.39928193-39928311)
MyoD P1-F CAGTCGCCCCCATGGACTTACTGG 57.9 306 g.11578803-11579108 exon 1
P1-R TGGTGGTCTTCCTCTTGCA (g.11578814-11579374)
P2-F AGAGGAAGACCACCAACG 55.8 246 g.11579094-11579339 exon 1
P2-R CATCTGACTCCCCGCTGT (g.11578814-11579374)
P3-F TGCGTGAGCAGGAGGATG 54.8 148 g.11579280-11579422 exon 1
P3-R GCTGGACTACAAGGAGGAA (g.11578814-11579374)
P4-F ATTCTGGCTGCTTCATTC 55.4 280 g.11580360-11580639 exon 2
P4-R CCCATCCTCCGTGCTTCA (g.11580512-11580590)
P5-F CAGTGGATAAATGGATCG 56.2 300 g.11581171-11581470 exon 3
P5-R CCTTTATAGCACTTGGTAG (g.11581208-11581467)
MyoG P1-F ATGGCACCGAGCAGTTGG 57.9 242 g.960256-960497 exon 1
P1-R GGGCAGCGTCGAGTCCTT (g.960321-960800)
P2-F TGACCCTGTGCCCTGAA 57.9 275 g.960442-960716 exon 1
P2-R GCGCTCGATGTACTGGATGG (g.960321-960800)
P3-F TTCGAGGCTCTGAAACGC 58.6 274 g.960624-960897 exon 1
P3-R CGCCCACAGGGACAACAA (g.960321-960800)
P4-F GTTTGGCTGCTGAAGGTG 55.6 168 g.962628-962795 exon 2
P4-R GTCCCTGAGATGGATGCT (g.962676-962757)
P5-F AGTAGGGTCGTGGGGTCA 56.4 274 g.963161-963434 exon 3
P5-R CACAGTGTCGGAGGGGTA (g.963210-963331)

SNP, single-nucleotide polymorphism.

a

The location of the Gallus gallus genomic sequence.

b

The exon location was determined from the published DNA sequence in GenBank, including NC_006088.3 (Myf5), NC_006088.3 (Mrf4), NC_006113.3 (MyoG), and NC_006092.3 (MyoD1).

Table 2

Frequency of alleles and genotypes, PIC, and the Hardy-Weinberg (H-W) equilibrium test at polymorphic sites

Gene SNP Genotype N Frequency Allele Frequency PIC H-W test, P
Mrf4 g.39928301 TT 152 0.4222 T 0.6000 0.3648 0.015
T>G TG 128 0.3556 G 0.4000
GG 80 0.2222
MyoD g.11579368 CC 166 0.4611 C 0.6250 0.3593 0.006
C>T CT 118 0.3278 T 0.3250
TT 76 0.2111

PIC, polymorphism information content; SNP, single-nucleotide polymorphism.

Table 3

Association of the combined genotypes with muscle fiber traits in chickens

Combined genotypes N Diameter of pectoralis muscle fiber1 (μm) Density of pectoralis muscle fiber1 (N/mm2) Diameter of thigh muscle fiber1 (μm) Density of thigh muscle fiber1 (N/mm2)
TT-CC 114 24.55±0.31 c 1,120.35±12.03 ab 26.88±0.41 c 938.16±28.64 a
TT-CT 26 26.28±0.43 a 1,121.71±34.99 ab 28.62±0.64 a 893.48±80.41 a
TG-CC 38 24.86±0.45 bc 1,185.01±19.12 a 27.41±0.66 b 896.33±40.72 a
TG-CT 76 25.27±0.39 b 1,117.23±18.74 ab 27.48±0.54 b 949.05±43.21 a
GG-TT 50 26.33±0.43 a 1,092.03±26.41 b 28.91±0.56 a 901.95±33.31 a
p-value
H×CG2 0.626 0.496 0.391 0.431
1

Means in a column and source of variation without a common letter differ, p<0.05.

2

H×CG represent interaction effects of housing system and combined genotypes.