Genetic effects of polymorphisms in myogenic regulatory factors on chicken muscle fiber traits.
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
The line SD02xSD03 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 [+ or -] 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.9x1.0x0.7 m; 5.3 birds/[m.sup.2]). Indoor concrete floor pens were covered with wood shavings (3x5x2.7 m, 3.3 birds/[m.sup.2]), 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 [m.sup.2]) with a simple house (6.7 birds/[m.sup.2]) 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 cmx2 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 pm mounted on slides and stained with hematoxylineosin. Fiber diameter measurements were performed on each section under light-microscopy (40x) 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 [micro]L (100 ng/[micro]L) of template DNA, 5 [micro]L of 2xTaq PCR MasterMix (TianGen Biochemical Technology, Shanghai, China), 0.4 [micro]L of each primer (10 pmol/[micro]L), and ddH2O to 10 [micro]L reaction volume using the following condition: 95[degrees]C for 10 min, 35 cycles of 94[degrees]C for 30 s, annealing from 57[degrees]C to 62[degrees]C for 45 s, 72[degrees]C for 1 min, and a final extension at 72[degrees]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 [micro]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.
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:
[Y.sub.ijk] = [mu] + [C.sub.i] + [H.sub.j] + [S.sub.k] + [(C x H).sub.ij] + [(H x S).sub.jk] + [(C x S).sub.ik] + [(C x H x S).sub.ijk] + [e.sub.ijk],
where [Y.sub.ijk] is the trait measured in chicken; [mu] is the population means of that trait; [C.sub.i] is the fixed effect associated with the combined genotypes i; [H.sub.j] is the fixed effect associated with the housing system j, [S.sub.k] is the fixed effect of sex k; [e.sub.ijk] is the residual random effect, and the C x H, H x S, C x S, and C x H x 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 [+ or -] standard error of the mean, and differences were considered significant at p<0.05.
[FIGURE 1 OMITTED]
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.
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.
[FIGURE 2 OMITTED]
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).
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.
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.
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).
Branciari, R., C. Mugnai, R. Mammoli, D. Miraglia, D. Ranucci, A. Dal Bosco, and C. Castellini. 2009. Effect of genotype and rearing system on chicken behavior and muscle fiber characteristics. J. Anim. Sci. 87:4109-4117.
Braun, T., E. Bober, B. Winter, N. Rosenthal, and H. H. Arnold. 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., G. Buschhausen-Denker, E. Bober, E. Tannich, and H. H. Arnold. 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, R. L., H. Weintraub, and A. B. Lassar. 1987. Expression of a single transfected cdna converts fibroblasts to myoblasts. Cell 51:987-1000.
Dou, T. C., S. Shi, R. Sun, H. J. Wang, H. J. Sun, and K. H. Wang. 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, D. G. and E. N. Olson. 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, A. C., P. B. Pillai, P. Y. Hester, C. Falcone, J. A. Mench, C. M. Owens, and J. L. Emmert. 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., J. M. Oh, I.-W. Kim, Z. E. Sauna, A. M. Calcagno, S. V. Ambudkar, and M. M. Gottesman. 2007. A "silent" polymorphism in the mdr1 gene changes substrate specificity. Science 315(5811):525-528.
Kolmodin, R., E. Strandberg, H. Jorjani, and B. Danell. 2003. Selection in the presence of a genotype by environment interaction: Response in environmental sensitivity. Anim. Sci. 76:375-386.
Komar, A. A. 2007. Silent snps: Impact on gene function and phenotype. Pharmacogenomics 8:1075-1080.
Massari, M. E. and C. Murre. 2000. Helix-loop-helix proteins: Regulators of transcription in eucaryotic organisms. Mol. Cell. Biol. 20:429-440.
Orita, M., Y. Suzuki, T. Sekiya, and K. Hayashi. 1989. Rapid and sensitive detection of point mutations and DNA polymorphisms using the polymerase chain reaction. Genomics 5:874-879.
Sing, C. F. 1990. Introduction to quantitative genetics. Am. J. Hum. Gent. 46(6): 1231.
Singh, K. and F. J. Dilworth. 2013. Differential modulation of cell cycle progression distinguishes members of the myogenic regulatory factor family of transcription factors. FEBS. J. 280: 3991-4003.
Wang, K. H., S. R. Shi, T. C. Dou, and H. J. Sun. 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., R. Davis, S. Tapscott, M. Thayer, M. Krause, R. Benezra, T. K. Blackwell, D. Turner, R. Rupp, and S. Hollenberg. 1991. The myod gene family: Nodal point during specification of the muscle cell lineage. Science 251(4995): 761-766.
Yin, H. D., Y. Wang, Z. C. Zhang, Y. P. Liu, S. Y. Chen, and Q. Zhu. 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, H. D., E. R. Gilbert, S. Y. Chen, Y. Wang, Z. C. Zhang, X. L. Zhao, Y. Zhang, and Q. Zhu. 2013. Effect of hybridization on carcass traits and meat quality of erlang mountainous chickens. Asian Australas. J. Anim. Sci. 26:1504-1510.
Yin, H., Z. Zhang, X. Lan, X. Zhao, Y. Wang, and Q. Zhu. 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.
Submitted Sept. 25, 2014; Revised Dec. 22, 2014; Accepted Jan. 12, 2015
Zhi-Qin Yang, Ying Qing, Qing Zhu, Xiao-Ling Zhao (1), Yan Wang, Di-Yan Li, Yi-Ping Liu, and Hua-Dong Yin *
Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu Campus, Chengdu, Sichuan 611130, China
* Corresponding Author: Hua-Dong Yin. Tel: +86-835-2882006, Fax: +86-835-2886080, E-mail: email@example.com
(1) College of Animal Science and Technology, Sichuan Agricultural University, Ya'an, Sichuan, 625014, China.
Table 1. Primers for SNP detection Gene Primer Primer sequence Annealing Product (5' [right arrow] 3') temperature length ([degrees]C) (bp) Myf5 P1-F ACTTTGTCTCGGGTCGCTAAT 57.4 240 P1-R TCAGGGCAGGCAGGCTCA P2-F CCCACTGAGCCTGCCTGC 57.2 294 P2-R TTCCCTCAAGAGCTCCTGG P3-F GGAACGCCATCAGGTACATCG 56.9 141 P3-R GCCCTGCTCTTACCATCACATC P4-F CCCCTGCTCGCCATTTGT 61.0 240 P4-R CCTTACCGTGGGGCATCT P5-F TAAGCGGCTCGCCCCGGTGC 61.7 221 P5-R CCCCTCATAGCGCCTGGT Mrf4 P1-F TCTGGGTTGCTCCTCGGGGTGT 62.2 253 P1-R TGCTCCTCGCCGCTGCTGTC P2-F CCCCGTGTCAGGACCAAC 59.5 240 P2-R CCACAGTCCGCCTTTTTCAG P3-F TCTCAAAAGGCGGACTGTGG 58.4 253 P3-R GGAGGAGGAGGGGAAGGAG P4-F CGCTGCTCGCCTCGCTAA 63.0 286 P4-R GGCCGACGACTCCACCTA P5-F CCGCTCCCTAACCCCTCT 60.1 260 P5-R AAAACGGGAGATGCGGAG MyoD P1-F CAGTCGCCCCCATGGACTTACTGG 57.9 306 P1-R TGGTGGTCTTCCTCTTGCA P2-F AGAGGAAGACCACCAACG 55.8 246 P2-R CATCTGACTCCCCGCTGT P3-F TGCGTGAGCAGGAGGATG 54.8 148 P3-R GCTGGACTACAAGGAGGAA P4-F ATTCTGGCTGCTTCATTC 55.4 280 P4-R CCCATCCTCCGTGCTTCA P5-F CAGTGGATAAATGGATCG 56.2 300 P5-R CCTTTATAGCACTTGGTAG MyoG P1-F ATGGCACCGAGCAGTTGG 57.9 242 P1-R GGGCAGCGTCGAGTCCTT P2-F TGACCCTGTGCCCTGAA 57.9 275 P2-R GCGCTCGATGTACTGGATGG P3-F TTCGAGGCTCTGAAACGC 58.6 274 P3-R CGCCCACAGGGACAACAA P4-F GTTTGGCTGCTGAAGGTG 55.6 168 P4-R GTCCCTGAGATGGATGCT P5-F AGTAGGGTCGTGGGGTCA 56.4 274 P5-R CACAGTGTCGGAGGGGTA Gene Primer Amplified Location of exon (b) region (a) Myf5 P1-F g.39931832-39932071 exon 1 P1-R (g.39931917-39932423) P2-F g.39932049-39932342 exon 1 P2-R (g.39931917-39932423) P3-F g.39932296-39932436 exon 1 P3-R (g.39931917-39932423) P4-F g.39932612-39932851 exon 2 P4-R (g.39932770-39932845) P5-F g.39932915-39933135 exon 3 P5-R (g.39932938-3993131) Mrf4 P1-F g.39927179-39927431 exon 1 P1-R (g.39927250-39927768) P2-F g.39927374-39927613 exon 1 P2-R (g.39927250-39927768) P3-F g.39927594-39927846 exon 1 P3-R (g.39927250-39927768) P4-F g.39927933-39928218 exon 2 P4-R (g.39927984-39928074) P5-F g.39928161-39928420 exon 3 P5-R (g.39928193-39928311) MyoD P1-F g.11578803-11579108 exon 1 P1-R (g.11578814-11579374) P2-F g.11579094-11579339 exon 1 P2-R (g.11578814-11579374) P3-F g.11579280-11579422 exon 1 P3-R (g.11578814-11579374) P4-F g.11580360-11580639 exon 2 P4-R (g.11580512-11580590) P5-F g.11581171-11581470 exon 3 P5-R (g.11581208-11581467) MyoG P1-F g.960256-960497 exon 1 P1-R (g.960321-960800) P2-F g.960442-960716 exon 1 P2-R (g.960321-960800) P3-F g.960624-960897 exon 1 P3-R (g.960321-960800) P4-F g.962628-962795 exon 2 P4-R (g.962676-962757) P5-F g.963161-963434 exon 3 P5-R (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 Mrf4 g.39928301 TT 152 0.4222 T T > G TG 128 0.3556 G GG 80 0.2222 MyoD g. 11579368 CC 166 0.4611 C C > T CT 118 0.3278 T TT 76 0.2111 Gene SNP Frequency PIC H-W test, P Mrf4 g.39928301 0.6000 0.3648 0.015 T > G 0.4000 MyoD g. 11579368 0.6250 0.3593 0.006 C > T 0.3250 PIC, polymorphism information content; SNP, single-nucleotide polymorphism. Table 3. Association of the combined genotypes with muscle fiber traits in chickens Combined N Diameter of Density of genotypes pectoralis pectoralis muscle fiber muscle fiber (1) ([micro]m) (1) (N/ [mm.sup.2]) TT-CC 114 24.55 [+ or -] 1,120.35 [+ or -] 0.31 (c) 12.03 (a)(b) TT-CT 26 26.28 [+ or -] 1,121.71 [+ or -] 0.43 (a) 34.99 (a)(b) TG-CC 38 24.86 [+ or -] 1,185.01 [+ or -] 0.45 (bc) 19.12 (a) TG-CT 76 25.27 [+ or -] 1,117.23 [+ or -] 0.39 (b) 18.74 (a)(b) GG-TT 50 26.33 [+ or -] 1,092.03 [+ or -] 0.43 (a) 26.41 (b) p-value HxCG (2) 0.626 0.496 Combined Diameter of Density of genotypes thigh muscle thigh muscle fiber (1) fiber1 (N/ ([micro]m) [mm.sup.2]) TT-CC 26.88 [+ or -] 938.16 [+ or -] 0.41 (c) 28.64 (a) TT-CT 28.62 [+ or -] 893.48 [+ or -] 0.64 (a) 80.41 (a) TG-CC 27.41 [+ or -] 896.33 [+ or -] 0.66 (b) 40.72 (a) TG-CT 27.48 [+ or -] 949.05 [+ or -] 0.54 (b) 43.21 (a) GG-TT 28.91 [+ or -] 901.95 [+ or -] 0.56 (a) 33.31 (a) HxCG (2) 0.391 0.431 (1) Means in a column and source of variation without a common letter differ, p < 0.05. (2) HxCG represent interaction effects of housing system and combined genotypes.
|Printer friendly Cite/link Email Feedback|
|Author:||Yang, Zhi-Qin; Qing, Ying; Zhu, Qing; Zhao, Xiao-Ling; Wang, Yan; Li, Di-Yan; Liu, Yi-Ping; Yin, Hua|
|Publication:||Asian - Australasian Journal of Animal Sciences|
|Date:||May 27, 2015|
|Previous Article:||Multiple genes related to muscle identified through a joint analysis of a two-stage genome-wide association study for racing performance of 1,156...|
|Next Article:||Proteomic analysis of bovine pregnancy-specific serum proteins by 2D fluorescence difference gel electrophoresis.|