The Identification of Novel Regions for Reproduction Trait in Landrace and Large White Pigs Using a Single Step Genome-Wide Association Study
Total Page:16
File Type:pdf, Size:1020Kb
Open Access Asian-Australas J Anim Sci Vol. 31, No. 12:1852-1862 December 2018 https://doi.org/10.5713/ajas.18.0072 pISSN 1011-2367 eISSN 1976-5517 The identification of novel regions for reproduction trait in Landrace and Large White pigs using a single step genome-wide association study Rattikan Suwannasing1, Monchai Duangjinda1,*, Wuttigrai Boonkum1, Rutjawate Taharnklaew2, and Komson Tuangsithtanon3 * Corresponding Author: Monchai Duangjinda Objective: The purpose of this study was to investigate a single step genome-wide association Tel: +66-43-202362, Fax: +66-43-202361, E-mail: [email protected] study (ssGWAS) for identifying genomic regions affecting reproductive traits in Landrace and Large White pigs. 1 Department of Animal Science, Faculty of Agriculture, Methods: The traits included the number of pigs weaned per sow per year (PWSY), the Khon Kaen University, Khon Kaen 40002, Thailand 2 Research and Development Center Betagro Group, number of litters per sow per year (LSY), pigs weaned per litters (PWL), born alive per litters Pathumthani 12120, Thailand (BAL), non-productive day (NPD) and wean to conception interval per litters (W2CL). A 3 Betagro Hybrid International Company Limited, total of 321 animals (140 Landrace and 181 Large White pigs) were genotyped with the Illumina Bangkok 10210, Thailand Porcine SNP 60k BeadChip, containing 61,177 single nucleotide polymorphisms (SNPs), ORCID while multiple traits single-step genomic BLUP method was used to calculate variances of Rattikan Suwannasing 5 SNP windows for 11,048 Landrace and 13,985 Large White data records. https://orcid.org/0000-0002-6950-4384 Monchai Duangjinda Results: The outcome of ssGWAS on the reproductive traits identified twenty-five and twenty- https://orcid.org/0000-0001-7044-8271 two SNPs associated with reproductive traits in Landrace and Large White, respectively. Three Wuttigrai Boonkum known genes were identified to be candidate genes in Landrace pigs including retinol binding https://orcid.org/0000-0003-2239-2300 Rutjawate Taharnklaew protein 7, and ubiquitination factor E4B genes for PWL, BAL, W2CL, and PWSY and one https://orcid.org/0000-0001-5526-3324 gene, solute carrier organic anion transporter family member 6A1, for LSY and NPD. Mean- while, five genes were identified to be candidate genes in Large White, two of which, aldehyde Submitted Jan 18, 2018; Revised Mar 20, 2018; dehydrogenase 1 family member A3 and leucine rich repeat kinase 1, associated with all of Accepted May 15, 2018 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. Keywords: Genomics; Single Step Genome-wide Association Study (ssGWAS); Candidate Genes; Single Nucleotide Polymorphisms (SNPs); Swine INTRODUCTION The genetic improvement is one approach to improve reproductive performance. However, the reproduction traits are characteristic of low heritability and difficult using conventional selection method to improve. The conventional selection method may provide a lower accu- racy, in comparison to the whole genome selection [1]. In the current development of molecular technique, such as single nucleotide polymorphism (SNP) chip has been widely used in for genome-wide association study (GWAS) to be a powerful tool in the identification of genomic regions or quantitative trait loci (QTL) related to an economically important trait. Single-step GWAS (ssGWAS) is the new GWAS approach which utilized all data (genotypes, phenotypes, Copyright © 2018 by Asian-Australasian Journal of Animal Sciences This is an open-access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and repro- 1852 duction in any medium, provided the original work is properly cited. www.ajas.info 4 4 79 because79 of becausethe lack ofof thedata lack for ofsome data animals. for some A animals.total of 11,048 A total Landraceof 11,048 and Landrace 13,985 and Large 13,985 White Large pigs White were collectedpigs were collected 80 data80 records data with records 13,351 with and 13,351 16,731 and pedigree 16,731 recordspedigree respectively records respectively. These sows. These were sows raised were in sixraised farms. in sixThe farms. The 81 reproduction81 reproduction traits included traits PWSY, included LSY, PWSY, NPD, LSY, W2C NPD,L, BAL W,2C andL, PBALWL, wereand P recordedWL were to recorded a maximum to a maximumof 10 parities of 10 parities 82 between82 200between6 through 200 2016 through5. PWSY 201 was5. PWSY calculated was ascalculated litter per as sow litter per per year sow multiplied per year multipliedby pigs weaned by pigs per weaned litter; LSY per litter; LSY 83 was83 calculated was calculatedas (the number as (the days number of gestating days of divided gestating by divided115 days) by divided115 days) by divided(the number by (t heof numberdays in theof days breeding in the breeding 84 herd84 divided herd by divided365 days); by 365W2CL, days); BAL, W2CL, and PWLBAL, were and PWLcollected were on collected every litter on every and divided litter and by dividedtotal litters; by total the NPDlitters; the NPD 85 was85 calculated was ascalculated 365 minuses as 365 productive minuses productivedays multiplied days multipliedby LSY, productive by LSY, productivedays was the day totals was number the total of numberdays that of days that 86 all 86gilts andall sowsgilts wereand sows either were gestating either orgestating lactating. or Thenlactating. checks Then the checks normal the distribution normal distribution before basic before statistical basic statistical 87 analyses87 suchanalyses as mean, such standard as mean, deviation standard and deviation coefficient and coefficientof variation. of variation. 88 88 89 Genotype89 Genotypedata and quality data and control quality control Suwannasing et al (2018) Asian-Australas J Anim Sci 31:1852-1862 90 Used90 321 Usedanimals; 321 140 animals; Landrace 140 andLan drace181 Large and 181White Large pigs White were pigsgenotyped were genotypedwith the Illumina with the Porcine Illumina SNP Porcine 60k SNP 60k 91 BeadChip, contained 61,177 SNPs. The quality control for genotypes SNPs of each breed following criteria was: SNP and pedigree information) jointly in one step, proposed91 BeadChip, by containof eached 61,177 breed SNPs. following The quality criteria control was: for genotypes SNP call SNPs rates of <0.90, each breed following criteria was: SNP Wang et al [2]. This approach can use for many92 models call and92 rates <call0.90, genotyperates genotype <0.90, call call genotype rates rates <0.90,< call0.90 rates, minorminor <0.90 allele allele, m freqinor frequencies uenciesallele freq < 0.05,<0.05,uencies Monomorphic Mono <0.05,- Monomorphic and checks and parent checks-progeny parent -progeny computing is fast and simplicity [2]. In pigs have been GWAS morphic and checks parent-progeny Mendelian conflicts were study using SNP chip in reproduction traits especially93 Mendelian litter93 Mendelianconflictsselected were conflicts for selected further were for selected analysis.further foranalysis A further total. A analysis oftotal 129 of .Landrace129 A total Landrace of 129and and Landrace175 175 Large and 175White Large pigs White with a pigs total with a total trait such as the number of born alive (NBA), total94 numberof 47,59094 andof 47,590 47,865Large and WhiteSNPs 47,865 respectively pigs SNPs with respectively a were total available of 47,590 were for available and the 47,865 genome for theSNPs-wide genome respecassociation-wide- association analyses in analyses this study in. this study. born, mummy (MUM), stillborn (SB) and total litter birth tively were available for the genome-wide association analyses weight [3,4]. While, no previous literature has 95studied the95 in this study. GWAS of pig weaned per sow per year (PWSY), litter96 perGenome sow96 -wideGenome association-wide association analysis analysis per year (LSY), pigs weaned per litters (PWL), born alive per Genome-wide association analysis litters (BAL), non-productive day (NPD), and wean97 to concepSingle97- -stepSingle genomeSingle-step-step-wide genome association genome-wide-wide association study :association The studygenome: The -studywide genome association: The-wide genome-wide association analysis was analysis estimated was by estimated using single by using-step single -step tion interval per litters (W2CL). Therefore, finding98 the genomicgenomic98 BLUPgenomicassociation (ssGBLUP) BLUP (ssGBLUP)analysis [2]. GWAS was [2] by estimated. GWASssGBLUP by by ssGBLUPcan using be called single-step can ssGWAS. be called genomic InssGWAS. this methodology, In this methodology, multivariate multivariate and and regions and candidate genes in the regions of signifi cant SNPs BLUP (ssGBLUP) [2]. GWAS by ssGBLUP can be called ss- from the whole genome that related to the reproduction99 separate trait99 breedseparate GWAS.analyzed breed Inthe analyzed this data methodology,. The the statistical data. The modelmultivariatestatistical was modelused: and was separate used: breed can be used as a powerful tool for selection to obtain100 the high100 analyzed the data. The statistical