Skip to main content Accessibility help

Genome-wide association study of bone mineral density trait among three pig breeds

  • B. Jiang (a1), M. Wang (a1) (a2), Z. Tang (a1), X. Du (a3), S. Feng (a1), G. Ma (a2), D. Ye (a1), H. Cheng (a1), H. Wang (a3) and X. Liu (a1) (a2)...


Leg weakness (LW) issues are a great concern for pig breeding industry. And it also has a serious impact on animal welfare. To dissect the genetic architecture of limb-and-hoof firmness in commercial pigs, a genome-wide association study was conducted on bone mineral density (BMD) in three sow populations, including Duroc, Landrace and Yorkshire. The BMD data were obtained by ultrasound technology from 812 pigs (including Duroc 115, Landrace 243 and Yorkshire 454). In addition, all pigs were genotyped using genome-by-sequencing and a total of 224 162 single-nucleotide polymorphisms (SNPs) were obtained. After quality control, 218 141 SNPs were used for subsequent genome-wide association analysis. Nine significant associations were identified on chromosomes 3, 5, 6, 7, 9, 10, 12 and 18 that passed Bonferroni correction threshold of 0.05/(total SNP numbers). The most significant locus that associated with BMD (P value = 1.92e−14) was detected at approximately 41.7 Mb on SSC6 (SSC stands for Sus scrofa chromosome). CUL7, PTK7, SRF, VEGFA, RHEB, PRKAR1A and TPO that are located near the lead SNP of significant loci were highlighted as functionally plausible candidate genes for sow limb-and-hoof firmness. Moreover, we also applied a new method to measure the BMD data of pigs by ultrasound technology. The results provide an insight into the genetic architecture of LW and can also help to improve animal welfare in pigs.


Corresponding author


Hide All
Bethel, M, Barnes, CLT, Taylor, AF, Cheng, YH, Chitteti, BR, Horowitz, MC, Bruzzaniti, A, Srour, EF and Kacena, MA 2015. A novel role for thrombopoietin in regulating osteoclast development in humans and mice. Journal of Cellular Physiology 230, 21422151.
Blott, SC, Swinburne, JE, Sibbons, C, Fox-Clipsham, LY, Helwegen, M, Hillyer, L, Parkin, TD, Newton, J and Vaudin, M 2014. A genome-wide association study demonstrates significant genetic variation for fracture risk in Thoroughbred racehorses. BMC Genomics 15, 147152.
Browning, SR and Browning, BL 2007. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. American Journal of Human Genetics 81, 10841097.
Chen, H, Wang, C, Conomos, MP, Stilp, AM and Lin, X 2016. Control for population structure and relatedness for binary traits in genetic association studies via logistic mixed models. American Journal of Human Genetics 98, 653666.
Dan, H, O’Sullivan, J, Urquhart, J, Daly, S, Bhaskar, SS, Biesecker, LG, Cole, T, Kirk, J, Chandler, K and Donnai, D 2011. Exome sequencing identifies CCDC8 mutations in 3-M syndrome, suggesting that CCDC8 contributes in a pathway with CUL7 and OBSL1 to control human growth. American Journal of Human Genetics 89, 148153.
de Sevilla, XF, Fabrega, E, Tibau, J and Casellas, J 2009. Genetic background and phenotypic characterization over two farrowings of leg conformation defects in Landrace and Large White sows. Journal of Animal Science 87, 16061612.
Duan, X, Murata, Y, Liu, Y, Nicolae, C, Olsen, BR and Berendsen, AD 2015. Vegfa regulates perichondrial vascularity and osteoblast differentiation in bone development. Development 142, 19841991.
Fan, B, Onteru, SK, Mote, BE, Serenius, T, Stalder, KJ and Rothschild, MF 2009. Large-scale association study for structural soundness and leg locomotion traits in the pig. Genetics Selection Evolution 41, 1418.
Fukawa, K and Kusuhara, S 2001. The genetic and non-genetic aspects of leg weakness and osteochondrosis in pigs. Journal of Animal Science 14, 114122.
Guo, YM, Ai, HS, Ren, J, Wang, GJ, Wen, Y, Mao, HR, Lan, LT, Ma, JW, Brenig, B and Rothschild, MF 2009. A whole genome scan for quantitative trait loci for leg weakness and its related traits in a large F2 intercross population between White Duroc and Erhualian. Journal of Animal Science 87, 15691575.
Gurgul, A, Mikszacybulska, A, Szmatoła, T, Jasielczuk, I, Piestrzyńskakajtoch, A, Fornal, A, Semikgurgul, E and Bugnoponiewierska, M 2018. Genotyping-by-sequencing performance in selected livestock species. Genomics 43, 620628.
Hayes, M, Gao, X, Yu, LX, Paria, N, Henkelman, RM, Wise, CA and Ciruna, B 2014. ptk7 mutant zebrafish models of congenital and idiopathic scoliosis implicate dysregulated Wnt signalling in disease. Nature Communications 5, 47774780.
Hu, K, Olsen, BR, Hu, K and Olsen, BR 2016. Osteoblast-derived VEGF regulates osteoblast differentiation and bone formation during bone repair. Journal of Clinical Investigation 126, 509.
Ibeagha-Awemu, EM, Peters, SO, Akwanji, KA, Imumorin, IG and Zhao, X 2016. High density genome wide genotyping-by-sequencing and association identifies common and low frequency SNPs, and novel candidate genes influencing cow milk traits. Scientific Reports 6, 3110931115.
Jensen, TB, Kristensen, HH and Toft, N 2012. Quantifying the impact of lameness on welfare and profitability of finisher pigs using expert opinions. Livestock Science 149, 209214.
Joshi, NA and Fass, JN 2011. Sickle: A sliding-window, adaptive, quality-based trimming tool for FastQ files (Version 1.33) [Software]. Retrieved from
Komori, T 2018. Runx2, an inducer of osteoblast and chondrocyte differentiation. Histochemistry and Cell Biology 149, 313323.
Laenoi, W, Uddin, MJ, Cinar, MU, Großebrinkhaus, C, Tesfaye, D, Jonas, E, Scholz, AM, Tholen, E, Looft, C and Wimmers, K 2011. Quantitative trait loci analysis for leg weakness-related traits in a Duroc × Pietrain crossbred population. Genetics Selection Evolution 43, 17.
Li, H and Durbin, R 2009. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 17541760.
Li, X, Masinde, G, Gu, W, Wergedal, J, Mohan, S and Baylink, DJ 2002. Genetic dissection of femur breaking strength in a large population (MRL/MpJ×SJL/J) of F2 mice: single QTL Effects, epistasis, and pleiotropy. Genomics 79, 734740.
Li, Z, Pei, XH, Yan, J, Yan, F, Cappell, KM, Whitehurst, AW and Xiong, Y 2014. CUL9 mediates the functions of the 3M complex and ubiquitylates survivin to maintain genome integrity. Molecular Cell 54, 805819.
Liu, X, Huang, M, Fan, B, Buckler, ES and Zhang, Z 2016. Iterative usage of fixed and random effect models for powerful and efficient genome-wide association studies. Plos Genetics 12, e1005767.
Mckenna, A, Hanna, M, Banks, E, Sivachenko, A, Cibulskis, K, Kernytsky, A, Garimella, K, Altshuler, D, Gabriel, S and Daly, M 2010. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Research 20, 12971303.
Meazza, C, Lausch, E and Pagani, S 2013. 3-M syndrome associated with growth hormone deficiency: 18. Italian Journal of Pediatrics 39, 16.
Moris, M, Peretz, A, Tjeka, R, Negaban, N and Bergmann, P 1995. Quantitative ultrasound bone measurements: normal values and comparison with bone mineral density by dual X-ray absorptiometry. Calcified Tissue International 57, 610.
Pavel, E, Nadella, K, Towns, WH and Kirschner, LS 2008. Mutation of Prkar1a causes osteoblast neoplasia driven by dysregulation of protein kinase A. Molecular Endocrinology 22, 430435.
Pértille, F, Guerrero-Bosagna, C, Silva, VHD, Boschiero, C, Nunes, JDRDS, Ledur, MNC, Jensen, P and Coutinho, LL 2016. High-throughput and cost-effective chicken genotyping using next-generation sequencing. Scientific Reports 6, 2692926935.
Peter, P, Charles, YC, Gonzales, J, Anderson, J, Sakhaee, K and Rubin, C 1993. Measurement of intrinsic bone quality in vivo by reflection ultrasound: correction of impaired quality with slow-release sodium fluoride and calcium citrate. Journal of Bone and Mineral Research 8, 234240.
Purcell, S, Neale, B, Todd-Brown, K, Thomas, L, Ferreira, MAR, Bender, D, Maller, J, Sklar, P, Bakker, PIWD and Daly, MJ 2007. PLINK: a tool set for whole-genome association and population-based linkage analyses. American Journal of Human Genetics 81, 559575.
Racine, JS 2012. RStudio: a platform-independent IDE for R and sweave. Journal of Applied Econometrics 27, 167172.
Rangkasenee, N, Murani, E, Brunner, RM, Schellander, K, Cinar, MU, Luther, H, Hofer, A, Stoll, M, Witten, A, Ponsuksili, S and Wimmers, K 2013. Genome-wide association identifies TBX5 as candidate gene for osteochondrosis providing a functional link to cartilage perfusion as initial factor. Frontiers in Genetics 4, 7882.
Rothammer, S, Bernau, M, Kremer-Rücker, PV, Medugorac, I and Scholz, AM 2017. Genome-wide QTL mapping results for regional DXA body composition and bone mineral density traits in pigs. Archives Animal Breeding 60, 5159.
Shaikh, S, Shettigar, SKG, Kumar, S, Kantharia, S, Kurva, J and Cherian, S 2019. Novel mutation in Cul7 gene in a family diagnosed with 3M syndrome. Journal of Genetics 98, 2126.
Shanks 2016. A cross-sectional study of radial and tibial ultrasonography between different ethnic and physically active groups. Master thesis, Essex University, Kochester, Britain.
Tan, C, Wu, Z, Ren, J and Huang, Z 2017. Genome-wide association study and accuracy of genomic prediction for teat number in Duroc pigs using genotyping-by-sequencing. Genetics Selection Evolution 49, 3540.
Wallner, C, Schira, J, Wagner, JM, Schulte, M, Fischer, S, Hirsch, T, Richter, W, Abraham, S, Kneser, U and Lehnhardt, M 2015. Application of VEGFA and FGF-9 enhances angiogenesis, osteogenesis and bone remodeling in type 2 diabetic long bone regeneration. PLoS ONE 10, e0118823.
Wang, Y, Xuemin, C, Yiqiang, Z, Jing, F, Xiaoxiang, H, Ning, L and Peng, X 2017. Optimized double-digest genotyping by sequencing (ddGBS) method with high-density SNP markers and high genotyping accuracy for chickens. PLoS ONE 12, e0179073.
Watts, NB, Lewiecki, EM, Miller, PD and Baim, S 2008. National Osteoporosis Foundation 2008 Clinician’s guide to prevention and treatment of osteoporosis and the World Health Organization Fracture Risk Assessment Tool (FRAX): what they mean to the bone densitometrist and bone technologist. Journal of Clinical Densitometry the Official Journal of the International Society for Clinical Densitometry 11, 473477.
Zhang, M, Manchanda, PK, Wu, D, Wang, Q and Kirschner, LS 2014. Knockdown of PRKAR1A, the gene responsible for Carney complex, interferes with differentiation in osteoblastic cells. Molecular Endocrinology 28, 295300.
Zhou, Z, Liu, G and Su, L 2016. A new approach to detect epistasis utilizing parallel implementation of ant colony optimization by MapReduce framework. International Journal of Computer Mathematics 93, 511523.


Type Description Title
Supplementary materials

Jiang et al. supplementary material
Figures S1-S2 and Tables S1-S7

 Word (595 KB)
595 KB

Genome-wide association study of bone mineral density trait among three pig breeds

  • B. Jiang (a1), M. Wang (a1) (a2), Z. Tang (a1), X. Du (a3), S. Feng (a1), G. Ma (a2), D. Ye (a1), H. Cheng (a1), H. Wang (a3) and X. Liu (a1) (a2)...


Altmetric attention score

Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed.