- Research article
- Open Access
Imputation to whole-genome sequence using multiple pig populations and its use in genome-wide association studies
© The Author(s) 2019
- Received: 12 October 2018
- Accepted: 10 January 2019
- Published: 24 January 2019
Use of whole-genome sequence data (WGS) is expected to improve identification of quantitative trait loci (QTL). However, this requires imputation to WGS, often with a limited number of sequenced animals for the target population. The objective of this study was to investigate imputation to WGS in two pig lines using a multi-line reference population and, subsequently, to investigate the effect of using these imputed WGS (iWGS) for GWAS.
Phenotypes and genotypes were available on 12,184 Large White pigs (LW-line) and 4943 Dutch Landrace pigs (DL-line). Imputed 660 K and 80 K genotypes for the LW-line and DL-line, respectively, were imputed to iWGS using Beagle v.4.1. Since only 32 LW-line and 12 DL-line boars were sequenced, 142 animals from eight commercial lines were added. GWAS were performed for each line using the 80 K and 660 K SNPs, the genotype scores of iWGS SNPs that had an imputation accuracy (Beagle R2) higher than 0.6, and the dosage scores of all iWGS SNPs.
For the DL-line (LW-line), imputation of 80 K genotypes to iWGS resulted in an average Beagle R2 of 0.39 (0.49). After quality control, 2.5 × 106 (3.5 × 106) SNPs had a Beagle R2 higher than 0.6, resulting in an average Beagle R2 of 0.83 (0.93). Compared to the 80 K and 660 K genotypes, using iWGS led to the identification of 48.9 and 64.4% more QTL regions, for the DL-line and LW-line, respectively, and the most significant SNPs in the QTL regions explained a higher proportion of phenotypic variance. Using dosage instead of genotype scores improved the identification of QTL, because the model accounted for uncertainty of imputation, and all SNPs were used in the analysis.
Imputation to WGS using the multi-line reference population resulted in relatively poor imputation, especially when imputing from 80 K (DL-line). In spite of the poor imputation accuracies, using iWGS instead of a lower density SNP chip increased the number of detected QTL and the estimated proportion of phenotypic variance explained by these QTL, especially when dosage scores were used instead of genotype scores. Thus, iWGS, even with poor imputation accuracy, can be used to identify possible interesting regions for fine mapping.
Use of whole-genome sequence (WGS) data is expected to improve the detection of quantitative trait loci (QTL) because such data are expected to contain most of the causal single nucleotide polymorphisms (SNPs), as was shown in dairy cattle populations by using WGS data of 234 bulls . Improved QTL detection is even more important in pig breeding populations, since the QTL can be used to improve the accuracy of prediction in across- or multi-population scenarios , which is especially relevant for pig breeding programs where cross-breeding is practised.
To benefit from WGS data, a large population of animals with such data is needed. In spite of the decreasing costs of WGS , it is still relatively expensive to sequence a large number of animals. A less expensive approach to increasing the number of animals with WGS data is to impute from lower density SNP chips to WGS. With imputation, a smaller group of sequenced animals is required and the majority of the population can be genotyped with a lower density and cheaper SNP panel. Then, low-density SNP genotypes are imputed to WGS using the sequenced animals as reference population. In dairy cattle, several studies have shown that imputation to WGS is reliable even with a limited number of sequenced animals in the reference population, e.g. [4–7]. For example, for imputation from 777 K SNP genotypes to WGS, van Binsbergen et al.  obtained an imputation accuracy evaluated by cross-validation within the reference population of 0.83 with 90 sequenced Holstein bulls. Bouwman and Veerkamp  obtained an imputation accuracy, based on the correlation between true and imputed WGS (iWGS) of 0.83 for imputation using a multi-breed reference population consisting of 20 Holstein and 60 animals of three different breeds, and demonstrated that other breeds can improve imputation accuracy when the number of sequenced individuals from the target breeds is small. For most breeding companies, only a small number of animals is sequenced per line because, often, sequencing expenses must be divided across lines. In those cases, it might be beneficial to combine the WGS data across lines into one reference population for imputation. In addition to the size and composition of the reference population, accuracy of imputation depends on the size of the genotyping array, the extent of linkage disequilibrium (LD) in the reference and target population and the minor allele frequency (MAF) of the SNPs on the genotyping array [4–9]. As a result, combining populations for imputation may not provide sufficient imputation accuracy in all populations.
Inaccurate imputation can influence the results of follow-up analyses such as genome-wide association studies (GWAS), especially when the accuracy of imputation is ignored in those analyses. Two approaches to account for imputation errors are to filter SNPs based on imputation accuracy prior to analysis or to use dosage scores in the analyses. Using dosage scores means that all imputed SNPs are included in the analysis, although the power to detect associations with poorly imputed SNPs will be low compared to using accurately imputed SNPs. So far, there is little information on the accuracy of imputation to iWGS using data of a commercial pig breeding population and on the effect of using iWGS genotype or dosage scores in a GWAS. Therefore, the objectives of this study were (1) to investigate the accuracy of imputation to WGS in two pig lines using a multi-line reference population and a limited number of sequenced animals available for the target lines; and (2) to investigate the effect of using imputed WGS genotypes versus dosage scores in a GWAS on QTL detection.
The dataset used in this study was provided by Topigs Norsvin. Phenotypes for the number of teats recorded after birth on 12,184 Large White (LW-line) and 4943 Dutch Landrace (DL-line) pigs were available. We investigated this trait because records were available for both sexes and its heritability is relatively high, i.e. 0.4 . The phenotypic records were pre-corrected for fixed effects, i.e. herd-year-of-birth, sex, and the random effect of litter, which were estimated with a pedigree-based linear model by Lopes et al. . After correction the average numbers of teats (± SD) were 15.68 (± 0.98) and 15.71 (± 1.04) for the LW-line and DL-line, respectively.
For both the LW-line and the DL-line, Geneseek-Neogen GPPHD 80 K SNP genotypes were available for all animals with phenotypic records. In addition, genotypes from the Affymetrix Axiom porcine 660 K SNP chip were available for the 120 sires with the largest number of offspring in the set of genotyped animals of the LW-line. Using the 120 sires as reference population, within-line imputation from 80 K to 660 K SNP genotypes was performed for all phenotyped animals from the LW-line using FImpute v2.2 software  with the pedigree option. This resulted in an average imputation accuracy of ~ 0.99. Imputation from 80 K to 660 K SNP genotypes was not possible for the DL-line because none of the DL-line animals were genotyped with the 660 K SNP chip. To avoid confounding of results, imputation of the DL-line to 660 K using the LW-line 660 K reference population was not considered but for the subsequent GWAS, 660 K genotypes for the DL-line were generated from the iWGS dataset.
Quality control of the within-line genotypes (80 K and 660 K) consisted of excluding (1) insertions and deletions, (2) SNPs with a MAF lower than 0.01, (3) SNPs with a frequency lower than 0.1 for either one of the three genotypes, and (4) SNPs with missing map information (based on Sscrofa10.2).
Reference population for whole-genome sequence data
Whole-genome sequence data was available for 168 of the most influential boars with many offspring for 10 commercial Dutch and Norwegian (Topigs Norsvin) lines, including 36 individuals originating from the Landrace breed, 39 from the Large White breed, 60 from the Duroc breed, 16 from a synthetic breed, 13 from the Pietrain breed, one Large White Dutch Landrace crossbred animal, and three animals of unknown origin. The reference population included 12 DL-line animals and 32 LW animals, corresponding to the target populations. One DL-line animal and nine LW-line animals had both whole-genome sequence data and high-density genotypes.
Raw sequence data were mapped to the pig genome build Sscrofa10.2 (Ensembl72)  using the Burrows-Wheeler Aligner (BWA)-mem algorithm . The average sequence coverage across the complete reference population was 11.6 fold. SNPs, short insertions and deletions were called with the GATK unified genotype-caller  for the complete reference population using default settings, but, in addition, the standard minimum confidence threshold was set to 30.0, the standard emittance confidence threshold was set to 20.0, and the target coverage threshold for down-sampling to coverage was set to 200. Subsequently, all detected variants were filtered using VCFtools , retaining variants with read depth values (per individual) ranging from 4 to 35 and variants with an overall Phred Quality score higher than 20, excluding variants with more than 20% missing data, and removing insertions and deletions. In addition, the variants were thinned such that the distance between variant sites was not shorter than 3 bp. The final step included phasing the reference population and imputing missing genotypes in the sequence data using Beagle 4.1 with 10 phasing iterations .
Imputation to iWGS
Imputation of the 80 K genotypes of the DL-line and the 660 K genotypes of the LW-line to the iWGS was performed with Beagle 4.1  using the multi-line reference population. The default parameter settings for Beagle were used, except for setting the effective population size to 300 instead of the default of 1 million  because it is much smaller in livestock than in humans [19, 20]. The accuracy of imputation at the SNP level was assessed by the Beagle R2, which is the estimated squared correlation of the genotype score with the true genotype .
For further analyses, the same quality controls as for the 80 K and 660 K genotypes were applied to the iWGS data but one additional quality filter was applied to the iWGS scores to account for possible imputation errors by removing SNPs with a Beagle R2 ≤ 0.6. This threshold was chosen to maintain a balance between the average imputation accuracy and the number of SNPs removed. This filter was not applied to allele dosage scores, which were coded as any real value between 0 and 2, because dosage scores account for imputation uncertainty.
Genome-wide association study
To account for population structure, the GWAS p-values for each SNP density were adjusted by their corresponding genomic inflation factors , which were calculated for each SNP density as the median of the observed Chi squared test statistics divided by the expected median of the corresponding Chi squared distribution assuming 1 degree of freedom. The Chi square test statistics were calculated from the p-values. Significance thresholds were then established by applying a Bonferroni correction by dividing the expected probability of a type-1 error (α = 0.05) by the number of independent tests. Following Duggal et al.  and Ricard et al. , we assumed that the number of independent tests was equal to the number of independent chromosome segments, which was calculated using the formula proposed by Goddard et al. . For the LW-line, the number of independent chromosomal segments was 776.4, 648.1 and 782.3 for the 80 K and 660 K SNP panels and the iWGS genotype scores, respectively. For the DL-line, the number of independent chromosomal segments was 249.5, 277.4, and 280.6 for the 80 K and 660 K SNP panels and the iWGS genotype scores, respectively. Since, within a line, the numbers of independent chromosomal segments were similar between SNP densities, the same significance threshold was used for all densities within a line. As a result, a SNP was considered significant when it was associated with a − log10 (p value) higher than 4.2 and 3.7 for the LW-line and the DL-line, respectively.
To identify QTL regions, SNPs on each chromosome were ranked based on their p values and, starting with the SNP with the highest − log10 (p value), all significant SNPs within a 0.5-Mb region to the left and right of the SNP were assigned to that QTL region. This procedure was repeated until all significant SNPs were assigned to a QTL region. We chose this definition for a QTL region and assumed that significant SNPs that are more than 0.5 Mb apart belong to independent QTL regions because the average LD of commercial pig lines decreases to less than 0.3 when the SNPs are more than 0.5 Mb apart [28–30].
Number of SNPs used for GWAS after quality control for different SNP densities and use of imputed whole-genome sequence dosage or genotype scores
80 K SNP genotypes
660 K SNP genotypes
Descriptive statistics of results of the GWAS for the two lines using different SNP densities and imputed whole-genome sequence dosage or genotype scores
Number of QTL regions
Genomic inflation factors
80 K SNP genotypes
660 K SNP genotypes
The Manhattan plots for each SNP density and for both lines showed a clear peak on chromosome 7 (Figs. 3, 4), reaching a -log10(adjusted p-value) of at least 10 for each SNP density for both lines. The peak was located between 103 and 105 Mb but, within this window, the position of the most significant SNP differed between SNP densities. In addition, for the LW-line, strong significant QTL regions were identified on chromosomes 10 and 16 for all densities, and on chromosomes 2 and 12 for the higher SNP density. For the DL-line, strong significant QTL regions were identified on chromosome 12 and 16 for all densities, and on chromosomes 2 and 10 for the higher SNP density.
Very high genomic inflation factors were observed, especially for the LW-line (> 3). For both lines, genomic inflation factors slightly increased as SNP density increased from 80 K to 660 K and then dropped to the same level as for the 80 K chip when iWGS genotype scores were used (Table 2). With iWGS dosage scores, genomic inflation factors decreased even more, i.e., from 1.80 for 80 K SNP genotypes to 1.51 for iWGS dosage scores for the DL-line, and from 3.50 for 80 K to 3.18 for iWGS dosage scores for the LW-line. Thus, inflation of the test statistic was lowest when dosage scores were used.
The objective of this study was to investigate the accuracy of imputation to WGS in two pig lines when using a multi-line reference population and the numbers of sequenced animals that belonged to the target lines were small, and subsequently, to investigate the effect of using the resulting imputed WGS genotype or dosage scores in a GWAS. Imputation from 80 K or 660 K to WGS using a multi-line reference population resulted in only 40 to 50% of the SNPs having a Beagle R2 higher than 0.60 in the iWGS data. Nonetheless, when accounting for imputation inaccuracy by filtering iWGS genotypes or by using dosage scores, the number of QTL detected and the estimated proportion of phenotypic variance explained by these were larger compared to when conducting the GWAS using lower density SNP chip genotypes, especially when using iWGS dosage scores. In the following section, the factors that influenced imputation accuracy are discussed, followed by a discussion on the effect of using iWGS genotype or dosage scores on the results of GWAS.
Factors affecting accuracy of imputation
The average Beagle R2 found in this study was 0.39 for the DL-line and 0.49 for the LW-line, which were relatively low compared to other studies that have used multi-line reference populations but with a larger size (more than 242 individuals) [5–7], but similar to the accuracy of 0.46 obtained with 90 sequenced Holstein bulls by van Binsbergen et al. . In our study, 168 individuals were included in the reference population and not all haplotypes present in our target populations were represented in the reference population, which increased imputation errors [4, 7]. Van Binsbergen et al.  investigated imputation accuracy in three scenarios that differed in the number of animals in a single breed reference population. They showed that for imputation from 50 K to iWGS, the accuracy was 0.37 with a reference population of 45 Holstein cows, and that it increased to 0.46 when the reference population increased to 90 Holstein cows. A similar increase of imputation accuracy was observed when imputing from 777 K to iWGS. These imputation accuracies were comparable or even higher than those found in our study, although we used a larger reference population. However, our reference population consisted of animals from 10 lines and the two target populations were only represented by 12 DL-line animals and 32 LW-line animals. The latter is probably the main reason for the disappointing imputation accuracy obtained in this study.
Our hypothesis was that adding animals from other lines and using a multi-line reference population would improve imputation accuracy, which was demonstrated previously especially for SNPs with a low MAF in the target population but that are segregating in other breeds [4, 6, 31–33]. For example, Bouwman and Veerkamp  showed that adding 60 individuals from the Jersey, Brown Swiss and Nordic Red Dairy cattle breeds to a reference population of 20 Holstein individuals increased the imputation accuracy of 777 K genotypes to iWGS for Holsteins from 0.71 to 0.83. In our study, the original small reference population was augmented by adding animals from other lines but imputation accuracy was still low for many chromosome regions. The genetic distances of the target populations with the other lines in the reference population were maybe too large because of, for example, different breeding goals. Given the small number of reference animals from the target lines (12 DL-line and 32 LW-line), accuracy of imputation would probably benefit more from additional sequenced DL-line and LW-line animals than from additional animals of other lines.
The imputation accuracy was lower for the DL-line than for the LW-line. Apart from the fact that there were only 12 DL-line individuals in the reference population compared to 32 LW-line individuals, the difference between starting density (i.e. 80 K genotypes for the DL-line and 660 K (imputed) genotypes for the LW-line) and target density (i.e. WGS) also contributed to the lower accuracy for the DL-line versus the LW-line. With a lower density SNP panel, the LD between the SNPs on the genotyping panel and WGS is lower and there is less information to identify shared haplotypes precisely, which increases the uncertainty of imputed genotypes [4, 33]. Stepwise imputation from 80 to 660 K to iWGS, instead of from 80 K to iWGS directly, as suggested by van Binsbergen et al. , and as performed for the LW-line, could improve imputation accuracy for the DL-line. However, this was not possible here because the number of available 660 K genotypes for the DL-line was not sufficient.
Other possible reasons for the difference in imputation accuracies between the two lines are differences in population structure or genetic architecture between the lines. For example, selection or different effective population sizes could have resulted in different LD decays and numbers of independent chromosome segments between the lines . Populations with a smaller number of independent chromosome segments are expected to have a higher imputation accuracy because they are expected to have less LD decay across the genome and a greater number of shared haplotypes. In this study, the DL-line had a smaller number of independent chromosomal segments (280.6) than the LW-line (782.3), which was expected to increase imputation accuracy for the DL-Line. However, due to other factors such as the starting density and smaller representation in the reference population, we did not observe a higher imputation accuracy for the DL-Line.
In addition to the above-mentioned factors, imputation accuracy could also be affected by the unequal distribution of SNPs on the genotyping array along the genome and by mapping errors. The latter complicate imputation because incorrect positions of SNPs lead to errors in haplotypes and LD structure in the region they are incorrectly mapped to. Thus, an improved reference genome should increase imputation accuracy [1, 35].
GWAS with iWGS genotype and dosage scores
Iso-Touru et al.  and Daetwyler et al.  identified new QTL when using iWGS genotype and dosage scores compared to SNP panel genotypes. Similarly, we also found new QTL regions when iWGS genotype scores were used. For the DL-line and the LW-line, 48.9 and 64.4%, respectively, of the QTL detected with iWGS genotype scores were not identified with the 80 K SNP genotypes. The QTL that were identified with all genotype densities were also reported in other GWAS studies that used DL-line or LW-line SNP genotypes [37–41]. In all these studies, the most significant QTL region was located on chromosome 7 and can be linked to the Vert[n]in gene. This gene is important for vertebral development and is positively correlated with number of teats in pigs [37, 39, 41]. Several new QTL regions identified with iWGS genotype and dosage scores include candidate genes for number of teats. For example, at approximately 125 Mb on chromosome 2, a QTL region was identified for the LW-line based on iWGS that includes the casein kinase 1 gamma 3 gene, which plays a regulatory role in the Wnt signalling pathway, which is essential for early mammary gland formation [42–44]. For the DL-line, we found a QTL region on chromosome 2 at about 76.5 Mb that was close to another casein kinase 1 gamma gene (casein kinase 1 gamma 2).
Compared to the DL-line, more QTL regions were identified for the LW-line when 660 K genotypes and genotype scores were used but fewer QTL regions were identified when 80 K genotypes or dosage scores were used for GWAS. This was unexpected because the DL-line has a smaller population size and more DL-line iWGS genotype scores were removed because of low imputation accuracy. Therefore, the power to detect associations was expected to be lower in the DL-line than in the LW-line . The larger number of QTL regions identified for the DL-line with 80 K genotypes and dosage scores might be because it has a smaller number of independent chromosomal segments (280.6) compared to the LW-line (782.3). A smaller number of independent chromosomal segments is expected to increase imputation accuracy because of less LD decay across the genome. As a result, power to detect associations might have been higher in the DL-line when using the 80 K genotypes and dosage scores. For 660 K and genotype scores, the power was lower in the DL-line compared to the LW-line because many SNPs with low imputation accuracy were removed. However, another possibility is that the smaller number of independent chromosomal segments for the DL-line, nay have resulted in the identification of a larger number of false positive QTL.
Of the 94 QTL regions found for the DL-line with iWGS genotype scores, 24% overlapped with QTL regions identified for the LW-line with iWGS genotype scores. Although power was limited in each line and we did not expect to detect each QTL in each line, differences might also be caused by differences in genetic architecture between the lines, such as MAF and LD patterns . Comparing the MAF of the most significant SNPs that were identified in the DL-line with the same SNPs in the LW-line (see Additional file 2: Figure S3), clearly showed that the QTL regions had different MAF between the two lines. For example, a SNP on chromosome 8 located at 104.1 Mb had a MAF of 0.005 and 0.36 in the DL-line and LW-line, respectively. Assuming that the SNP is associated with the phenotype in both lines, its low MAF reduces the power to detect it in the DL-line.
Although a large number of new QTL were found using iWGS genotype, whether they are indeed new associations or artefacts of the definition of the QTL region can be questioned. Here, a QTL region was defined as the 0.5-Mb region to the left and right of the most significant SNP in a region because multiple studies have found that the average LD in commercial pig lines decreases below 0.3 when the SNPs are more than 0.5 Mb apart [28–30]. This definition may increase the number of QTL regions for iWGS because iWGS consists of many SNPs that are in very high LD. In addition, LD decay can vary across and within chromosomes  and, therefore, this definition can be too strict for some QTL regions and neighbouring regions could be merged into one region or the other way around. To test this, we increased the size of the QTL region to 1 Mb to the right and left of the most significant SNP in that region. This resulted in a reduction in the number of QTL detected with iWGS genotype scores from 104 to 64 for the DL-line and from 94 to 56 for the LW-line. Although this is a significant reduction, it is still larger than the number of QTL regions based on the original QTL region definition found with 80 K and 660 K.
In addition to the definition of the QTL region, the newly identified QTL regions could also be an artefact due to the significance threshold that we applied in this study. Here, the significance threshold was set by applying a Bonferroni correction using the number of independent chromosomal segments instead of the commonly used total number of SNPs, which does not take LD between SNPs into account. The number of independent chromosomal segments has often been used for Bonferroni correction for GWAS in human [26, 46, 47], plant [48–50], and animal [27, 51] genetics. However, so far there has been little consensus about the most appropriate approach for testing significance for GWAS. This should be a topic for future research.
Finally, even with low imputation accuracy, GWAS using iWGS can be beneficial in several cases. For example, the QTL regions identified in this study could help to pre-select SNPs to improve the accuracy of genomic predictions, as shown by [52–54]. In addition, QTL regions identified could be used as indicators for fine-mapping of possible interesting regions, even if they are based on poorly imputed WGS data. For example, SNPs that are in the detected QTL regions could be added to high-density SNP chips, or animals with extreme phenotypes could be sequenced for the detected QTL regions.
Genomic inflation factor
In this study, the genomic inflation factors were lower when using dosage scores compared to using 80 K and 660 K genotypes and iWGS genotype scores. The scenario using dosage scores included SNPs with a low imputation accuracy, and hence SNPs with a low MAF. In the other scenarios, there is ascertainment bias, which is caused by the preference of SNPs on a chip that are more likely to be common . In addition, SNPs selected for iWGS genotype scores have ascertainment bias, because accurately imputed SNPs in general have a higher MAF. By (indirectly) selecting SNPs with a high MAF, it is easier to detect effects for these SNPs, and their surrounding SNPs, and therefore it is likely that more significant SNPs are found than expected based on the theoretical distribution of the test statistic distribution under the null hypothesis. Linkage disequilibrium between the SNPs will also result in more significant SNPs within the region that surrounds a causal variant. Both the higher frequency of SNPs with high MAF and long-range LD increase genomic inflation factors [56–58] and the rate of false positives. Use of dosage scores leads to less inflation of the test statistic because the ascertainment bias is partly removed, leading to less biased SNP effects compared to the use of lower density SNP chips and iWGS genotype scores.
Use of a multi-line reference population resulted in relatively poor imputation accuracy to iWGS, with 52.7% of the SNPs on the 660 K array in the LW-line and 39.1% of the SNPs on the 80 K array in the DL-line having Beagle R2 lower than 0.6. Imputation from 660 K to iWGS was more accurate than imputation from 80 K to iWGS, which suggests that step-wise imputation, i.e. first imputing from 80 to 660 K and then from 660 K to iWGS, could increase the accuracy of imputation. Although imputation accuracy was poor, using iWGS instead of genotypes from a lower density SNP chip increased the number of detected QTL regions and the expected proportion of phenotypic variance that they explained. When using iWGS dosage scores instead of genotype scores, even more QTL regions were detected because all SNPs were used in the analysis and the uncertainty of imputation was taken into account by the model.
SvdB, RFV, JV and FAvE contributed to the study design. SvdB did the data analyses, JV adapted GCTA to perform GWAS with dosage scores and ACB conducted the variant calling within the WGS data and assisted with the imputation. MSL provided the genotype and phenotype data. SvdB wrote the first draft of the manuscript and RFV, JV, ACB, MSL and FAvE provided feedback and helped to revise the manuscript. All authors read and approved the final manuscript.
The authors want to acknowledge Topigs Norsvin for providing the data.
The authors declare that they have no competing interests.
Consent for publication
The data used for this study were collected as part of routine data recording for a commercial breeding program. Samples collected for DNA extraction were used only for the breeding program. Data recording and sample collection were conducted strictly in line with the Dutch law on the protection of animals (Gezondheids- en welzijnswet voor dieren).
This study was financially supported by NWO-TTW, the Dutch Ministry of Economic Affairs (TKI Agri & Food project 16022, code BO-22.04-011-001-ASG-LR), and the Breed4Food partners Cobb Europe, CRV, Hendrix Genetics and Topigs Norsvin. The use of the HPC cluster was made possible by CAT-AgroFood (Shared Research Facilities, Wageningen University and Research).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Daetwyler HD, Capitan A, Pausch H, Stothard P, van Binsbergen R, Brøndum RF, et al. Whole-genome sequencing of 234 bulls facilitates mapping of monogenic and complex traits in cattle. Nat Genet. 2014;46:858–65.PubMedView ArticleGoogle Scholar
- Kizilkaya K, Fernando RL, Garrick DJ. Genomic prediction of simulated multibreed and purebred performance using observed fifty thousand single nucleotide polymorphism genotypes. J Anim Sci. 2014;88:544–51.View ArticleGoogle Scholar
- Wetterstrand KA. DNA Sequencing costs: data from the NHGRI Genome Sequencing Program (GSP). www.genome.gov/sequencingcostsdata. Accessed 26 June 2017.
- van Binsbergen R, Bink MCAM, Calus MPL, van Eeuwijk FA, Hayes BJ, Hulsegge I, et al. Accuracy of imputation to whole-genome sequence data in Holstein Friesian cattle. Genet Sel Evol. 2014;46:41.PubMedPubMed CentralView ArticleGoogle Scholar
- van Binsbergen R, Calus MPL, Bink MCAM, van Eeuwijk FA, Schrooten C, Veerkamp RF. Genomic prediction using imputed whole-genome sequence data in Holstein Friesian cattle. Genet Sel Evol. 2015;47:71.PubMedPubMed CentralView ArticleGoogle Scholar
- Bouwman AC, Veerkamp RF. Consequences of splitting whole-genome sequencing effort over multiple breeds on imputation accuracy. BMC Genet. 2014;15:105.PubMedPubMed CentralView ArticleGoogle Scholar
- Brøndum RF, Guldbrandtsen B, Sahana G, Lund MS, Su G. Strategies for imputation to whole genome sequence using a single or multi-breed reference population in cattle. BMC Genomics. 2014;15:728.PubMedPubMed CentralView ArticleGoogle Scholar
- Hayes BJ, Bowman PJ, Daetwyler HD, Kijas JW, Van der Werf JH. Accuracy of genotype imputation in sheep breeds. Anim Genet. 2012;43:72–80.PubMedView ArticleGoogle Scholar
- Zhang Z, Druet T. Marker imputation with low-density marker panels in Dutch Holstein cattle. J Dairy Sci. 2010;93:5487–94.PubMedView ArticleGoogle Scholar
- McKay RM, Rahnefeld GW. Heritability of teat number in swine. Can J Anim Sci. 1990;70:425–30.View ArticleGoogle Scholar
- Lopes MS, Bovenhuis H, van Son M, Nordbø Ø, Grindflek EH, Knol EF, et al. Using markers with large effect in genetic and genomic predictions. J Anim Sci. 2017;95:59–71.PubMedGoogle Scholar
- Sargolzaei M, Chesnais JP, Schenkel FS. FImpute-An efficient imputation algorithm for dairy cattle populations. J Dairy Sci. 2011;94:421.Google Scholar
- Groenen MA, Archibald AL, Uenishi H, Tuggle CK, Takeuchi Y, Rothschild MF, et al. Analyses of pig genomes provide insight into porcine demography and evolution. Nature. 2012;491:393–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv:1303.3997v1 [q-bio.GN]. 2013.
- McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297–303.PubMedPubMed CentralView ArticleGoogle Scholar
- Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, et al. The variant call format and VCFtools. Bioinformatics. 2011;27:2156–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Browning SR, Browning BL. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet. 2007;81:1084–97.PubMedPubMed CentralView ArticleGoogle Scholar
- Browning BL, Browning SR. Genotype imputation with millions of reference samples. Am J Hum Gen. 2016;98:116–26.View ArticleGoogle Scholar
- Hall SJG. Effective population sizes in cattle, sheep, horses, pigs and goats estimated from census and herdbook data. Animal. 2016;10:1778–85.PubMedView ArticleGoogle Scholar
- Pocrnic I, Lourenco DAL, Masuda Y, Legarra A, Misztal I. The dimensionality of genomic information and its effect on genomic prediction. Genetics. 2016;203:573–81.PubMedPubMed CentralView ArticleGoogle Scholar
- Browning BL, Browning SR. A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. Am J Hum Genet. 2009;84:210–23.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang J, Lee SH, Goddard ME, Visscher PM. GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet. 2011;88:76–82.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang J, Zaitlen NA, Goddard ME, Visscher PM, Price AL. Advantages and pitfalls in the application of mixed-model association methods. Nat Genet. 2014;46:100–6.PubMedPubMed CentralView ArticleGoogle Scholar
- Goddard ME, Hayes BJ, Meuwissen THE. Using the genomic relationship matrix to predict the accuracy of genomic selection. J Anim Breed Genet. 2011;128:409–21.PubMedView ArticleGoogle Scholar
- Yang J, Weedon MN, Purcell S, Lettre G, Estrada K, Willer CJ, et al. Genomic inflation factors under polygenic inheritance. Eur J Hum Genet. 2011;19:807–12.PubMedPubMed CentralView ArticleGoogle Scholar
- Duggal P, Gillanders EM, Holmes TN, Bailey-Wilson JE. Establishing an adjusted p-value threshold to control the family-wide type 1 error in genome wide association studies. BMC Genomics. 2008;9:516.PubMedPubMed CentralView ArticleGoogle Scholar
- Ricard A, Robert C, Blouin C, Baste F, Torquet G, Morgenthaler C, et al. Endurance exercise ability in the horse: a trait with complex polygenic determinism. Front Genet. 2017;8:89.PubMedPubMed CentralView ArticleGoogle Scholar
- Grossi DA, Jafarikia M, Brito LF, Buzanskas ME, Sargolzaei M, Schenkel FS. Genetic diversity, extent of linkage disequilibrium and persistence of gametic phase in Canadian pigs. BMC Genet. 2017;18:6.PubMedPubMed CentralView ArticleGoogle Scholar
- Veroneze R, Lopes PS, Guimarães SEF, Silva FF, Lopes MS, Harlizius B, et al. Linkage disequilibrium and haplotype block structure in six commercial pig lines. J Anim Sci. 2013;91:3493–501.PubMedView ArticleGoogle Scholar
- Badke YM, Bates RO, Ernst CW, Schwab C, Steibel JP. Estimation of linkage disequilibrium in four US pig breeds. BMC Genomics. 2012;13:24.PubMedPubMed CentralView ArticleGoogle Scholar
- Ni G, Strom TM, Pausch H, Reimer C, Preisinger R, Simianer H, et al. Comparison among three variant callers and assessment of the accuracy of imputation from SNP array data to whole-genome sequence level in chicken. BMC Genomics. 2015;16:824.PubMedPubMed CentralView ArticleGoogle Scholar
- VanRaden PM, Null DJ, Sargolzaei M, Wiggans GR, Tooker ME, Cole JB, et al. Genomic imputation and evaluation using high-density Holstein genotypes. J Dairy Sci. 2013;96:668–78.PubMedView ArticleGoogle Scholar
- Pausch H, MacLeod IM, Fries R, Emmerling R, Bowman PJ, Daetwyler HD, et al. Evaluation of the accuracy of imputed sequence variants and their utility for causal variant detection in cattle. Genet Sel Evol. 2017;49:24.PubMedPubMed CentralView ArticleGoogle Scholar
- Goddard ME. Genomic selection: prediction of accuracy and maximisation of long term response. Genetica. 2009;136:245–57.PubMedView ArticleGoogle Scholar
- Erbe M, Hayes BJ, Matukumalli LK, Goswami S, Bowman PJ, Reich CM, et al. Improving accuracy of genomic predictions within and between dairy cattle breeds with imputed high-density single nucleotide polymorphism panels. J Dairy Sci. 2012;95:4114–29.PubMedView ArticleGoogle Scholar
- Iso-Touru T, Sahana G, Guldbrandtsen B, Lund M, Vilkki J. Genome-wide association analysis of milk yield traits in Nordic Red Cattle using imputed whole genome sequence variants. BMC Genet. 2016;17:55.PubMedPubMed CentralView ArticleGoogle Scholar
- Duijvesteijn N, Veltmaat JM, Knol EF, Harlizius B. High-resolution association mapping of number of teats in pigs reveals regions controlling vertebral development. BMC Genomics. 2014;15:542.PubMedPubMed CentralView ArticleGoogle Scholar
- Bidanel JP, Rosendo A, Iannuccelli N, Riquet J, Gilbert H, Caritez JC, et al. Detection of quantitative trait loci for teat number and female reproductive traits in Meishan × Large White F2 pigs. Animal. 2008;2:813–20.PubMedView ArticleGoogle Scholar
- Ding N, Guo Y, Knorr C, Ma J, Mao H, Lan L, et al. Genome-wide QTL mapping for three traits related to teat number in a White Duroc × Erhualian pig resource population. BMC Genet. 2009;10:6.PubMedPubMed CentralView ArticleGoogle Scholar
- Guo YM, Lee GJ, Archibald AL, Haley CS. Quantitative trait loci for production traits in pigs: a combined analysis of two Meishan × Large White populations. Anim Genet. 2008;39:486–95.PubMedView ArticleGoogle Scholar
- Lopes MS, Bastiaansen JW, Harlizius B, Knol EF, Bovenhuis H. A genome-wide association study reveals dominance effects on number of teats in pigs. PLoS One. 2014;9:e105867.PubMedPubMed CentralView ArticleGoogle Scholar
- Chu EY, Hens J, Andl T, Kairo A, Yamaguchi TP, Brisken C, et al. Canonical WNT signaling promotes mammary placode development and is essential for initiation of mammary gland morphogenesis. Development. 2004;131:4819–29.PubMedView ArticleGoogle Scholar
- McKay RM, Peters JM, Graff JM. The casein kinase I family in Wnt signaling. Dev Biol. 2001;235:388–96.PubMedView ArticleGoogle Scholar
- Lindvall C, Zylstra CR, Evans N, West RA, Dykema K, Furge KA, et al. The Wnt co-receptor Lrp6 is required for normal mouse mammary gland development. PLoS One. 2009;4:e5813.PubMedPubMed CentralView ArticleGoogle Scholar
- Klein RJ. Power analysis for genome-wide association studies. BMC Genet. 2007;8:58.PubMedPubMed CentralView ArticleGoogle Scholar
- Demerath EW, Guan W, Grove ML, Aslibekyan S, Mendelson M, Zhou YH, et al. Epigenome-wide association study (EWAS) of BMI, BMI change and waist circumference in African American adults identifies multiple replicated loci. Hum Mol Genet. 2015;24:4464–79.PubMedPubMed CentralView ArticleGoogle Scholar
- Luciano M, Huffman JE, Arias-Vásquez A, Vinkhuyzen AA, Middeldorp CM, Giegling I, et al. Genome-wide association uncovers shared genetic effects among personality traits and mood states. Am J Med Genet B Neuropsychiatr Genet. 2012;159:684–95.View ArticleGoogle Scholar
- Bouchet S, Servin B, Bertin P, Madur D, Combes V, Dumas F, et al. Adaptation of maize to temperate climates: mid-density genome-wide association genetics and diversity patterns reveal key genomic regions, with a major contribution of the Vgt2 (ZCN8) locus. PLoS One. 2013;8:e71377.PubMedPubMed CentralView ArticleGoogle Scholar
- Tello J, Torres-Pérez R, Grimplet J, Carbonell-Bejerano P, Martínez-Zapater JM, Ibáñez J. Polymorphisms and minihaplotypes in the VvNAC26 gene associate with berry size variation in grapevine. BMC Plant Biol. 2015;15:253.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang W, Guo Z, Huang C, Duan L, Chen G, Jiang N, et al. Combining high-throughput phenotyping and genome-wide association studies to reveal natural genetic variation in rice. Nat Commun. 2014;5:5087.PubMedPubMed CentralView ArticleGoogle Scholar
- Xie L, Luo C, Zhang C, Zhang R, Tang J, Nie Q, et al. Genome-wide association study identified a narrow chromosome 1 region associated with chicken growth traits. PLoS One. 2012;7:e30910.PubMedPubMed CentralView ArticleGoogle Scholar
- Raymond B, Bouwman AC, Schrooten C, Houwing-Duistermaat J, Veerkamp RF. Utility of whole-genome sequence data for across-breed genomic prediction. Genet Sel Evol. 2018;50:27.PubMedPubMed CentralView ArticleGoogle Scholar
- Brøndum RF, Su G, Janss L, Sahana G, Guldbrandtsen B, Boichard D, et al. Quantitative trait loci markers derived from whole genome sequence data increases the reliability of genomic prediction. J Dairy Sci. 2015;98:4107–16.PubMedView ArticleGoogle Scholar
- Shi J, Park JH, Duan J, Berndt ST, Moy W, Yu K, et al. Winners curse correction and variable thresholding improve performance of polygenic risk modeling based on summary-level data from genome-wide association studies. PLoS Genet. 2016;12:e1006493.PubMedPubMed CentralView ArticleGoogle Scholar
- Howie B, Fuchsberger C, Stephens M, Marchini J, Abecasis GR. Fast and accurate genotype imputation in genome-wide association studies through pre-phasing. Nat Genet. 2012;44:955–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Devlin B, Roeder K. Genomic control for association studies. Biometrics. 1999;55:997–1004.PubMedPubMed CentralView ArticleGoogle Scholar
- Hinrichs AL, Larkin EK, Suarez BK. Population stratification and patterns of linkage disequilibrium. Genet Epidemiol. 2009;33:S88–92.PubMedPubMed CentralView ArticleGoogle Scholar
- Reich DE, Goldstein DB. Detecting association in a case-control study while correcting for population stratification. Genet Epidemiol. 2001;20:4–16.PubMedView ArticleGoogle Scholar