Identification of growth trait related genes in a Yorkshire purebred pig population by genome-wide association studies

  • Meng, Qingli (Beijing Breeding Swine Center) ;
  • Wang, Kejun (Department of Animal Genetics and Breeding, National Engineering Laboratory for Animal Breeding, MOA Laboratory of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University) ;
  • Liu, Xiaolei (Key Laboratory of Agricultural Animal Genetics, Breeding, and Reproduction of Ministry of Education and Key Laboratory of Swine Genetics and Breeding of Ministry of Agriculture, Huazhong Agricultural University) ;
  • Zhou, Haishen (Beijing Breeding Swine Center) ;
  • Xu, Li (Beijing Breeding Swine Center) ;
  • Wang, Zhaojun (Beijing Breeding Swine Center) ;
  • Fang, Meiying (Department of Animal Genetics and Breeding, National Engineering Laboratory for Animal Breeding, MOA Laboratory of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University)
  • Received : 2016.07.19
  • Accepted : 2016.10.24
  • Published : 2017.04.01


Objective: The aim of this study is to identify genomic regions or genes controlling growth traits in pigs. Methods: Using a panel of 54,148 single nucleotide polymorphisms (SNPs), we performed a genome-wide Association (GWA) study in 562 pure Yorshire pigs with four growth traits: average daily gain from 30 kg to 100 kg or 115 kg, and days to 100 kg or 115 kg. Fixed and random model Circulating Probability Unification method was used to identify the associations between 54,148 SNPs and these four traits. SNP annotations were performed through the Sus scrofa data set from Ensembl. Bioinformatics analysis, including gene ontology analysis, pathway analysis and network analysis, was used to identify the candidate genes. Results: We detected 6 significant and 12 suggestive SNPs, and identified 9 candidate genes in close proximity to them (suppressor of glucose by autophagy [SOGA1], R-Spondin 2 [RSPO2], mitogen activated protein kinase kinase 6 [MAP2K6], phospholipase C beta 1 [PLCB1], rho GTPASE activating protein 24 [ARHGAP24], cytoplasmic polyadenylation element binding protein 4 [CPEB4], GLI family zinc finger 2 [GLI2], neuronal tyrosine-phosphorylated phosphoinositide-3-kinase adaptor 2 [NYAP2], and zinc finger protein multitype 2 [ZFPM2]). Gene ontology analysis and literature mining indicated that the candidate genes are involved in bone, muscle, fat, and lung development. Pathway analysis revealed that PLCB1 and MAP2K6 participate in the gonadotropin signaling pathway and suggests that these two genes contribute to growth at the onset of puberty. Conclusion: Our results provide new clues for understanding the genetic mechanisms underlying growth traits, and may help improve these traits in future breeding programs.


Supported by : Beijing San Yuan Breeding Technology CO., LTD, Natural Science Foundation of China(NSFC)


  1. Fontanesi L, Schiavo G, Galimberti G, Calo DG, Russo V. A genomewide association study for average daily gain in Italian Large White pigs. J Anim Sci 2014;92:1385-94.
  2. Wang K, Liu D, Hernandez-Sanchez J, et al. Genome wide association analysis reveals new production trait genes in a male duroc population. PLoS ONE 2015;10:e0139207.
  3. Hu ZL, Park CA, Reecy JM. Developmental progress and current status of the animal QTLdb. Nucleic Acids Res 2016;44:D827-33.
  4. Andersson L. Genome-wide association analysis in domestic animals: a powerful approach for genetic dissection of trait loci. Genetica 2009;136:341-9.
  5. Sanchez MP, Tribout T, Iannuccelli N, et al. A genome-wide association study of production traits in a commercial population of Large White pigs: evidence of haplotypes affecting meat quality. Genet Sel Evol 2014;46:12.
  6. Fan B, Onteru SK, Du ZQ, et al. Genome-wide association study identifies Loci for body composition and structural soundness traits in pigs. PLoS ONE 2011;6:e14726.
  7. Luo W, Cheng D, Chen S, et al. Genome-wide association analysis of meat quality traits in a porcine Large White x Minzhu intercross population. Int J Biol Sci 2012;8:580-95.
  8. Onteru SK, Fan B, Nikkila MT, et al. Whole-genome association analyses for lifetime reproductive traits in the pig. J Anim Sci 2011;89:988-95.
  9. Sahana G, Kadlecova V, Hornshoj H, Nielsen B, Christensen OF. A genome-wide association scan in pig identifies novel regions associated with feed efficiency trait. J Anim Sci 2013;91:1041-50.
  10. Onteru SK, Gorbach DM, Young JM, et al. Whole genome association studies of residual feed intake and related traits in the pig. PLoS ONE 2013;8:e61756.
  11. Bai Y, Zhang JB, Xue Y, et al. Differential expression of CYB5A in Chinese and European pig breeds due to genetic variations in the promoter region. Anim Genet 2015;46:16-22.
  12. Liu X, Huang M, Fan B, Buckler ES, Zhang Z. Iterative usage of fixed and random effect models for powerful and efficient genome-wide association studies. PLoS Genet 2016;12:e1005767.
  13. Pearson TA, Manolio TA. How to interpret a genome-wide association study. JAMA 2008;299:1335-44.
  14. Purcell S, Neale B, Todd-Brown K, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 2007;81:559-75.
  15. Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 2009;4:44-57.
  16. Shan D, Li JL, Wu L, et al. GPAT3 and GPAT4 are regulated by insulin-stimulated phosphorylation and play distinct roles in adipogenesis. J Lipid Res 2010;51:1971-81.
  17. Onteru SK, Fan B, Du ZQ, et al. A whole-genome association study for pig reproductive traits. Anim Genet 2012;43:18-26.
  18. Wang JY, Luo YR, Fu WX, et al. Genome-wide association studies for hematological traits in swine. Anim Genet 2013;44:34-43.
  19. Bouaziz M, Ambroise C, Guedj M. Accounting for population stratification in practice: a comparison of the main strategies dedicated to genome-wide association studies. PLoS ONE 2011;6:e28845.
  20. Cowherd RB, Asmar MM, Alderman JM, et al. Adiponectin lowers glucose production by increasing SOGA. Am J Pathol 2010;177:1936-45.
  21. Do DN, Strathe AB, Ostersen T, Pant SD, Kadarmideen HN. Genome-wide association and pathway analysis of feed efficiency in pigs reveal candidate genes and pathways for residual feed intake. Front Genet 2014;5:307.
  22. Su ZJ, Hahn CN, Goodall GJ, et al. A vascular cell-restricted RhoGAP, p73RhoGAP, is a key regulator of angiogenesis. Proc Natl Acad Sci USA 2004;101:12212-7.
  23. Widmann P, Reverter A, Fortes MR, et al. A systems biology approach using metabolomic data reveals genes and pathways interacting to modulate divergent growth in cattle. BMC Genomics 2013;14:798.
  24. Yingling VR. A delay in pubertal onset affects the covariation of body weight, estradiol, and bone size. Calcif Tissue Int 2009;84:286-96.
  25. Yingling VR, Khaneja A. Short-term delay of puberty causes a transient reduction in bone strength in growing female rats. Bone 2006;38:67-73.
  26. Fujishiro M, Gotoh Y, Katagiri H, et al. Three mitogen-activated protein kinases inhibit insulin signaling by different mechanisms in 3T3-L1 adipocytes. Mol Endocrinol 2003;17:487-97.
  27. Engelman JA, Berg AH, Lewis RY, et al. Constitutively active mitogen-activated protein kinase kinase 6 (MKK6) or salicylate induces spontaneous 3T3-L1 adipogenesis. J Biol Chem 1999;274:35630-8.
  28. Fujishiro M, Gotoh Y, Katagiri H, et al. MKK6/3 and p38 MAPK pathway activation is not necessary for insulin-induced glucose uptake but regulates glucose transporter expression. J Biol Chem 2001;276:19800-6.
  29. O'Carroll SJ, Mitchell MD, Faenza I, Cocco L, Gilmour RS. Nuclear PLCbeta1 is required for 3T3-L1 adipocyte differentiation and regulates expression of the cyclin D3-cdk4 complex. Cell Signal 2009;21:926-35.
  30. Faenza I, Bavelloni A, Fiume R, et al. Expression of phospholipase C beta family isoenzymes in C2C12 myoblasts during terminal differentiation. J Cell Physiol 2004;200:291-6.
  31. Gerber AN, Wilson CW, Li YJ, Chuang PT. The hedgehog regulated oncogenes Gli1 and Gli2 block myoblast differentiation by inhibiting MyoD-mediated transcriptional activation. Oncogene 2007;26:1122-36.
  32. Yokoyama K, Tezuka T, Kotani M, et al. NYAP: a phosphoprotein family that links PI3K to WAVE1 signalling in neurons. EMBO J 2011;30:4739-54.
  33. Jack BH, Crossley M. GATA proteins work together with friend of GATA (FOG) and C-terminal binding protein (CTBP) co-regulators to control adipogenesis. J Biol Chem 2010;285:32405-14.
  34. Zhang Z, Zhang H, Pan RY, et al. Genetic parameters and trends for production and reproduction traits of a Landrace herd in China. J Integr Agric 2016;15:1069-75.
  35. Liu G, Kim JJ, Jonas E, et al. Combined line-cross and half-sib QTL analysis in Duroc-Pietrain population. Mamm Genome 2008;19:429-38.

Cited by

  1. Unravelling the genetic loci for growth and carcass traits in Chinese Bamaxiang pigs based on a 1.4 million SNP array pp.09312668, 2018,
  2. Extent of linkage disequilibrium and effective population size of Korean Yorkshire swine vol.31, pp.12, 2018,
  3. Scans for signatures of selection in Russian cattle breed genomes reveal new candidate genes for environmental adaptation and acclimation vol.8, pp.1, 2018,
  4. Genome-wide association study reveals genetic loci and candidate genes for average daily gain in Duroc pigs vol.31, pp.4, 2018,
  5. Genetic Architecture of Feeding Behavior and Feed Efficiency in a Duroc Pig Population vol.9, pp.1664-8021, 2018,
  6. Multimarker and rare variants genomewide association studies for bone weight in Simmental cattle vol.135, pp.3, 2018,