The identification of novel regions for reproduction trait in Landrace and Large White pigs using a single step genome-wide association study

  • Received : 2018.01.18
  • Accepted : 2018.05.15
  • Published : 2018.12.01


Objective: The purpose of this study was to investigate a single step genome-wide association study (ssGWAS) for identifying genomic regions affecting reproductive traits in Landrace and Large White pigs. Methods: The traits included the number of pigs weaned per sow per year (PWSY), the number of litters per sow per year (LSY), pigs weaned per litters (PWL), born alive per litters (BAL), non-productive day (NPD) and wean to conception interval per litters (W2CL). A total of 321 animals (140 Landrace and 181 Large White pigs) were genotyped with the Illumina Porcine SNP 60k BeadChip, containing 61,177 single nucleotide polymorphisms (SNPs), while multiple traits single-step genomic BLUP method was used to calculate variances of 5 SNP windows for 11,048 Landrace and 13,985 Large White data records. Results: The outcome of ssGWAS on the reproductive traits identified twenty-five and twenty-two SNPs associated with reproductive traits in Landrace and Large White, respectively. Three known genes were identified to be candidate genes in Landrace pigs including retinol binding protein 7, and ubiquitination factor E4B genes for PWL, BAL, W2CL, and PWSY and one gene, solute carrier organic anion transporter family member 6A1, for LSY and NPD. Meanwhile, five genes were identified to be candidate genes in Large White, two of which, aldehyde dehydrogenase 1 family member A3 and leucine rich repeat kinase 1, associated with all of six reproduction traits and three genes; retrotransposon Gag like 4, transient receptor potential cation channel subfamily C member 5, and LHFPL tetraspan subfamily member 1 for five traits except W2CL. Conclusion: The genomic regions identified in this study provided a start-up point for marker assisted selection and estimating genomic breeding values for improving reproductive traits in commercial pig populations.


  1. Cleveland MA, Hickey JM, Forni S. A common data set for genomic analysis of livestock populations. G3: Genes Genomes Genetics 2012;2:429-35.
  2. Wang H, Misztal I, Aguilar I, Legarra A, Muir WM. Genomewide association mapping including phenotypes from relatives without genotypes. Genet Res Camb 2012;94:73-83.
  3. Onteru SK, Fan B, Du ZQ, et al. A whole-genome association study for pig reproductive traits. Anim Genet 2012;43:18-26.
  4. Wu P, Yang Q, Wang K, et al. Single step genome-wide association studies based on genotyping by sequence data reveals novel loci for the litter traits of domestic pigs. Genomics 2018;110:171-9.
  5. Misztal I, Tsuruta S, Strabel T, et al. BLUPF90 and related programs (BGF90). In: Proceedings of the 7th World Congress on Genetics Applied to Livestock Production; 2002 Aug 19-23: Montpellier, France.
  6. Aguilar I, Misztal I, Johnson DL, et al. Hot topic: a unified approach to utilize phenotypic, full pedigree, and genomic information for genetic evaluation of Holstein final score. J Dairy Sci 2010;93:743-52.
  7. VanRaden PM. Efficient methods to compute genomic predictions. J Dairy Sci 2008;91:4414-23.
  8. Vitezica ZG, Aguilar I, Misztal I, Legarra A. Bias in genomic predictions for populations under selection. Genet Res (Camb) 2011;93:357-66.
  9. Zhang Z, Liu J, Ding X, et al. Best linear unbiased prediction of genomic breeding values using a trait-specific marker-derived relationship matrix. PLoS ONE 2010;5:e12648.
  10. Beissinger TM, Rosa GJM, Kaeppler SM, Gianola D, de Leon, N. Defining window-boundaries for genomic analyses using smoothing spline techniques. Genet Sel Evol 2015;47:30.
  11. Diniz DB, Lopes MS, Broekhuijse MLW, et al. A genome-wide association study reveals a novel candidate gene for sperm motility in pigs. Anim Reprod Sci 2014;151:201-7.
  12. Hu ZL, Park CA, Reecy JM. Developmental progress and current status of the Animal QTLdb. Nucleic Acids Res 2016;44:D827-33.
  13. Howard JT, Jiao S, Tiezzi F, et al. Genome-wide association study on legendre random regression coefficients for the growth and feed intake trajectory on Duroc Boars. BMC Genet 2015;16:59.
  14. Suzuki T, Toyohara T, Akiyama Y, et al. Transcriptional regulation of organic anion transporting polypeptide SLCO4C1 as a New Therapeutic Modality to prevent chronic kidney disease. J Pharm Sci 2011;100: 3696-707.
  15. Lee SY, Williamson B, Cabllero OL, et al. Identification of the gonad-specific anion transporter SLCO6A1 as a cancer/testis (CT) antigen expressed in human lung cancer. Cancer Immun 2004;4:13.
  16. Gong WH, Tang ZL, Han JL, et al. Mapping, tissue distribution and polymorphism of porcine retinol binding protein genes (RBP5 and RBP7). Asian-Australas J Anim Sci 2008;21:1544-50.
  17. Kim M, Seo H, Choi Y, et al. Microarray analysis of gene expression in the uterine endometrium during the implantation period in pigs. Asian-Australas J Anim Sci 2012;25:1102-16.
  18. Hu C, Keen HL, Lu KT, et al. Retinal binding protein 7 mediates an anti-oxidant response to cardiovascular stressors by regulating $PPAR{\gamma}$ activity and adiponectin in endothelium. FASEB J 2017;31:1015.14.
  19. Kaneko-Oshikawa C, Nakagawa T, Yamada M, et al. Mammalian E4 is required for cardiac development and maintenance of the nervous system. Mol Cell Biol 2005;25:10953-64.
  20. Zage PE, Sirisaengtaksin N, Liu Y, et al. UBE4B levels are correlated with clinical outcomes in neuroblastoma patients and with altered neuroblastoma cell proliferation and sensitivity to EGFR inhibitors. Cancer 2013;119:915-23.
  21. Duan JJ, Cai J, Guo YF, Bian XW, Yu SC. ALDH1A3, a metabolic target for cancer diagnosis and therapy. Int J Cancer 2016;139:965-75.
  22. Canestro C, Catchen JM, Rodriguez-Mari A, Yokoi H, Postlethwait JH. Consequences of lineage-specific gene loss on functional evolution of surviving paralogs: ALDH1A and retinoic acid signaling in vertebrate genomes. PLoS Genet 2009;5:e1000496.
  23. Dupe V, Matt N, Garnier JM, et al. A newborn lethal defect due to inactivation of retinaldehyde dehydrogenase type 3 is prevented by maternal retinoic acid treatment. Proc Natl Acad Sci USA 2003;100:14036-41.
  24. Bottje W, Kong BW. Cell biology symposium: feed efficiency: mitochondrial function to global gene expression. J Anim Sci 2013;91:1582-93.
  25. Iida A, Xing W, Docx MKF, et al. Identification of biallelic LRRK1 mutations in osteosclerotic metaphyseal dysplasia and evidence for locus heterogeneity. J Med Genet 2016;53:568-74.
  26. Irie M, Yoshikawa M, Ono R, et al. Cognitive function related to the Sirh11/Zcchc16 gene acquired from an LTR retrotransposon in Eutherians. PLoS Genet 2015;11:e1005521.
  27. Henke C, Strissel PL, Schubert MT, et al. Selective expression of sense and antisense transcripts of the sushi-ichi-related retrotransposon - derived family during mouse placentogenesis. Retrovirology 2015;12:9.
  28. Lau OC, Shen B, Wong CO, et al. TRPC5 channels participate in pressure-sensing in aortic baroreceptors. Nat Commun 2016;7:11947.
  29. Gao Y, Yao T, Deng Z, et al. TrpC5 mediates acute leptin and serotonin effects via Pomc neurons. Cell Rep 2017;18:583-92.
  30. Guo YM, Zhang XF, Ren J, et al. A joint genome wide association analysis of pig leg weakness and its related traits in an F2 population and a Sutai population. J Anim Sci 2013;91:4060-8.
  31. Huang C, Guo J, Liu S, et al. Isolation, tissue distribution and prokaryotic expression of a novel human X-linked gene LHFPL1. DNA Seq 2004;15:299-302.