The SWine IMputation (SWIM) haplotype reference panel enables nucleotide resolution genetic mapping in pigs

Communications Biology volume 6, Article number: 577 (2023) Cite this article

192 Accesses

2 Altmetric

Metrics details

Genetic mapping to identify genes and alleles associated with or causing economically important quantitative trait variation in livestock animals such as pigs is a major goal in animal genetic improvement. Despite recent advances in high-throughput genotyping technologies, the resolution of genetic mapping in pigs remains poor due in part to the low density of genotyped variant sites. In this study, we overcame this limitation by developing a reference haplotype panel for pigs based on 2259 whole genome-sequenced animals representing 44 pig breeds. We evaluated software combinations and breed composition to optimize the imputation procedure and achieved an average concordance rate in excess of 96%, a non-reference concordance rate of 88%, and an r2 of 0.85. We demonstrated in two case studies that genotype imputation using this resource can dramatically improve the resolution of genetic mapping. A public web server has been developed to allow the pig genetics community to fully utilize this resource. We expect this resource to facilitate genetic mapping and accelerate genetic improvement in pigs.

The domestic pig (Sus scrofa) is an important livestock species and a model organism for biomedical research1. Historically, domestication and intense artificial selection have created many pig breeds that are genetically and phenotypically distinct from each other and from their wild relatives2,3,4. More recently, high-throughput DNA sequencing and genotyping technologies5 have facilitated the genetic improvement of pigs. For example, hundreds of genome-wide association and quantitative trait locus (QTL) mapping studies have identified numerous genomic regions associated with various production, physiological, and behavioral phenotypes6. These studies are important for understanding the genetic and biological basis of economically and biomedically important traits such as growth7, fertility8, and disease resistance9.

The resolution of genetic mapping in pigs remains poor due in part to the low density of single nucleotide polymorphism (SNP) genotyping arrays. One proven, cost-effective approach to overcome the limitation in resolution is through genotype imputation, leveraging linkage disequilibrium to infer genotypes at unobserved polymorphic loci10. With large haplotype reference panels created by whole-genome sequencing, imputation has the potential to provide sequence-level genotypes11. In livestock animals, where QTL identification and genetic prediction are two major goals, and linkage disequilibrium is extensive, sequence-level genotype imputation has been successfully applied with a relatively small number of reference haplotypes but decent accuracy12, 13. In pigs, in particular, at least two public imputation servers are available14, 15. However, they either contained a very limited number of animals in the reference panel14 or lacked good representation from major commercial breeds15, limiting their applications. In addition, although many studies have demonstrated improvement in mapping resolution16 and genomic prediction accuracy17, none of these can be publicly accessed.

In this study, we produced whole-genome sequence data from 1530 newly sequenced pigs and combined them with 729 additional animals from public databases to call variants and develop by far the largest and most diverse reference panel of haplotypes in pigs to date. This substantial increase in the number of available genomes allowed us to impute SNP array genotypes to whole genome sequences rapidly and accurately. We evaluated the accuracy of imputation and demonstrated the utility of this haplotype reference panel in genome-wide association mapping. We introduce a new public web server (swimgeno.org) where users may submit array genotypes and retrieve imputed whole-genome sequence-level genotypes. This resource will greatly improve access to high-accuracy genotype imputation, facilitating potentially nucleotide resolution genetic mapping in pigs.

We consolidated whole-genome sequence data from newly sequenced animals (n = 1530) and publicly available data (Supplementary Data 1 and 2) for a total of 2259 pigs, representing 44 different breeds (Supplementary Data 1). The majority of animals were Landrace (n = 651), Yorkshire (n = 543), and Duroc (n = 485), three major commercial breeds. The uniquely aligned sequence depth was approximately 12.86 X averaged across all animals (Supplementary Data 1). We called variants using the GATK pipeline and calibrated variant quality scores with known variant sets compiled from commercial SNP arrays. After filtering out variants of low quality and excessive heterozygosity and missingness, 47.86 M autosomal variants remained. Sub-sampling of animals indicated that the increase in the number of discovered variants quickly diminished (Fig. 1a). More than 95% of all variants could be recovered using only 1000 randomly selected animals.

a Number of variants discovered as a function of the number of animals in the discovery cohort. The curve is generated by randomly subsetting the population and counting DNA variants that remain polymorphic. The number of variants discovered using 1000 and the whole set of 2259 animals are marked. b Pairwise linkage disequilibrium in four domestic breeds and wild boars from three regions. Average r2 is plotted against the distance between variants. LD was calculated after removing low-frequency variants (MAF < 0.05) and close relatives (GRM > 0.5) in 435 Durocs, 522 Landraces, 493 Yorkshires, 36 Meishans, 24 European wild boars, and 27 Asian wild boars. c Scatter plot of first two principal components of genotype matrix for common (MAF > 0.05) and LD-pruned variants. Points are color-coded according to their reported breed information. A preliminary principal component analysis was performed to visually inspect and remove clear outliers from clusters, which indicated errors in breed information. d Ancestries of pigs were estimated with variable (K = 2, 4, 6) numbers of postulated ancestral populations using the ADMIXTURE software. Estimated ancestries were plotted as stacked bar charts with breeds annotated on the top. In addition to annotations above the bar chart, broad geographical locations are also annotated below the bar chart for K = 6.

Linkage disequilibrium (LD) between variants in this population was extensive but differed by breed (Fig. 1b). LD in wild boars declined more rapidly as the distance between variants increased than in domestic breeds, consistent with the high level of inbreeding among intensively selected domestic breeds (Fig. 1b). Genetic variation present in the pig genome separated breeds into distinct clusters that represented geographic differentiation (Fig. 1c, d). The first principal component of the genotypes separated Asian breeds and wild boars from their European counterparts, while the second separated Durocs from other breeds (Fig. 1c). Estimated ancestries of the breeds also indicated clearly separated clusters according to their geographical locations (Fig. 1d). Taken together, the diverse and rich genetic variation in the 2259 pig genomes included in this study provides a strong foundation for whole-genome imputation.

We focused on the ~34 M autosomal variants (30,489,782 SNPs and 4,125,579 indels) segregating at a minor allele frequency (MAF) > 0.005 to construct the haplotype reference panel. To investigate factors that influence imputation accuracy, we considered different combinations of commonly used phasing and imputation software, including SHAPEIT4/IMPUTE5, Beagle5.2/Beagle5.2, and Eagle2.4/Minimac4. We defined imputation accuracy using three metrics, the overall concordance rate between imputed and observed genotypes, non-reference concordance rate summarizing accuracy for non-reference genotypes only, and squared correlation (r2) between imputed and observed genotypes. We focused on Landrace as the target set because it has the largest number of animals in the dataset. We held out 100 Landrace pigs sequenced at high coverage (>15X) and compared observed genotypes with imputed genotypes starting from sequencing-based genotypes at sites on a 50 K SNP array (GeneSeek GGP). Regardless of breed composition in the haplotype reference panel of fixed size, SHAPEIT4/IMPUTE5 outperformed Beagle5.2/Beagle5.2 and Eagle2.4/Minimac4 in all three metrics (Fig. 2a–c). SHAPEIT4/IMPUTE5 was therefore chosen for all subsequent analyses.

a Concordance rate, non-reference concordance rate, and r2 of imputed versus observed genotypes using different software combinations with 550 Landraces as the reference panel. b Same analysis but in a reference panel consisting of 250 Landraces, 150 Durocs, and 150 Yorkshires. c Same analysis but in a reference panel consisting of 250 Landraces and 300 other breeds (not Duroc or Yorkshire).

In cattle, imputation using multi-breed reference panels appeared to be more accurate than using a single-breed panel12,18. However, multi-breed panels are confounded by larger sample sizes. We asked whether imputation using reference panels of the same size from a single breed and from a mixture of multiple breeds made a difference (Fig. 3a, compare L, DLY, and LO). This question was important as it informs whether to use a multi-breed or breed-specific reference panel to achieve optimal accuracy. We again considered 100 Landrace animals as the target set because of its relatively larger sample size. We found imputation accuracy measured by all three metrics to be remarkably similar (Fig. 3b–d) when the reference panel size was equal. Reference panel derived from the same breed as the target set had a very slight advantage (Fig. 3b–d). However, multi-breed panels are useful because reference from the same breed alone (but smaller sample size) was not able to achieve the same accuracy (Fig. 3, compare L-250 with others). Because the vast majority of Landrace pigs were from a single population, the imputation accuracy may not reflect a realistic scenario when new target sets are derived from other populations. We evaluated imputation accuracy using 550 animals as the reference set but 41 Landrace pigs from the SRA as the target set, thus representing a situation where the target sets are distant from the reference. Imputation accuracies were lower, and the multi-breed panel appeared to hold a small advantage (Supplementary Fig. 1). Expanding the reference panel to 2218 animals increased the accuracy substantially (Supplementary Fig. 2). The lower accuracies may be due to a combination of the small number of target animals as well as further genetic distance from the reference panel. Taken together, although the comparison between multi-breed and breed-specific panels of the same size depends on specific situations, a multi-breed reference panel is desired as opposed to a breed-specific reference panel in most cases as it maximizes reference panel size.

a Experimental design to investigate the effect of breed composition of haplotype reference panel on imputation accuracy. Three reference panels were tested, including ‘L’: 550 Landrace animals; ‘DLY’: 550 pigs from the Duroc, Landrace, and Yorkshire breeds; ‘LO’: 550 pigs from Landraces and other non-Duroc or Yorkshire breeds; ‘L-250’: 250 Landrace animals only. One hundred Landraces were used as the target set. Concordance rate (b), non-reference concordance rate (c), and r2 (d) of imputed versus observed genotypes using different breed compositions of the haplotype reference panel.

We compared our SWine IMputation (SWIM) resource using the multi-breed reference panel with an imputation server for pigs (PHARP) that utilized 1006 animals publicly available in the SRA15. We evaluated imputation accuracy among variants that were present in both reference panels. PHARP contained relatively few major commercial breeds, including 115 Yorkshires, 85 Durocs, and 48 Landraces. We considered target sets from Landrace, Duroc, and Yorkshire, in which the vast majority of GWAS are conducted (Fig. 4a). When evaluating imputation accuracy, we held out 100 animals as the target set and used the remainder (n = 2159) as the haplotype reference panel. While the overall concordance rate was uniformly high (>94.24%), imputation using the SWIM panel developed in the present study was consistently higher than PHARP within each breed (Fig. 4b). The improvement was much more pronounced when considering the non-reference concordance rate and r2, two metrics that more faithfully reflect the accuracy, especially at low frequency (Fig. 4c, d). The difference between SWIM and PHARP could simply be a sample size difference, especially for the breeds evaluated. The final reference haplotype panel consisting of all 2259 animals is expected to achieve a concordance rate in excess of 95.84%, a non-reference concordance rate of 88.26%, and an r2 of 0.85.

a Composition of the haplotype reference panels, including different scenarios of SWIM and PHARP, and target set. SWIM-100Y, SWIM-100D, and SWIM-100L are holding out 100 Yorkshires, 100 Durocs, and 100 Landraces as the target set, respectively. For the PHARP reference panel, the same 100 Yorkshires, 100 Durocs, and 100 Landraces are used to evaluate imputation accuracy. b Concordance rate of imputed versus observed genotypes using different haplotype reference panels. The mean concordance rate across all variants is also indicated on the plot for each reference panel. c Non-reference concordance rate of imputed versus observed genotypes using different haplotype reference panels. The mean non-reference concordance rate across all variants is also indicated on the plot for each reference panel. d r2 of imputed versus observed genotypes using different haplotype reference panels. Mean r2 across all variants is also indicated on the plot for each reference panel.

We also assessed the performance of different starting SNP chips, including the GeneSeek GGP 50K, Affymetrix Wens 55K, and Affymetrix Axiom PigHD 660K. These chips were chosen because the Wens 55K and GGP 50K have a similar number of SNPs but share fewer SNPs, and the Axiom PigHD represents a higher density. The imputation accuracies were evaluated in 100 Durocs and using 2159 animals as the reference (Supplementary Fig. 3a). After removal of SNPs whose probes did not map uniquely to the reference genome or were monomorphic, 39,491, 48,337, and 561,111 SNPs overlapped with the haplotype reference panel for the GeneSeek GGP, Wens, and Axiom PigHD, respectively (Supplementary Fig. 3b). As expected, higher density of SNPs led to higher imputation accuracy (Supplementary Fig. 3c–e) in all three metrics, with the Affymetrix PigHD 660K SNP chip achieving remarkably high accuracy at 99.50% overall concordance rate (Supplementary Fig. 3c), 98.63% non-reference concordance rate (Supplementary Fig. 3d), and 0.98 r2 (Supplementary Fig. 3e).

To demonstrate the usefulness of sequence-level genotype imputation in genetic mapping, we performed genome-wide association studies (GWAS) for two important growth traits in pigs, using both SNP arrays and imputed genotypes. The two traits, backfat thickness and body length, were chosen because putative causal genes and mutations have been previously well characterized. Our objective was to see if imputation-based GWAS was able to find previously validated functional genes and variants.

Backfat thickness (BF) is one of the most important economic traits in pigs and has been intensively interrogated for its genetic basis. Genomic heritabilities estimated using either array SNPs or imputed SNPs were similar and indicated a moderately heritable trait (Fig. 5a). Alleles in several genes, including IGF219,20, MC4R21, and LEPR22, have been consistently associated with BF variation in pigs. In particular, a missense mutation in the MC4R gene (chr1:160773437:G>A) has been suggested as the causative mutation21 and extensively replicated in multiple genetic backgrounds23. Furthermore, mutations in MC4R are strongly associated with early onset obesity in humans24, and its role in the regulation of energy homeostasis is well established25. Importantly, the putative causal mutation in MC4R has been included in one of the commercially available SNP genotyping arrays, the Geneseek GGP Porcine 50K SNP Chip (Neogen, Lincoln, NE). However, the same SNP is not present in the more widely used Illumina PorcineSNP60 chip. To see if genotype imputation was able to correctly impute the genotypes of this SNP, we excluded the MC4R SNP and imputed whole-genome genotypes from a population of 3769 Duroc pigs genotyped using the GGP Porcine 50K SNP arrays. Remarkably, the concordance rate and r2 between the imputed and array MC4R SNP genotypes were 99.71% and 0.9916, respectively. We performed GWAS using array and imputed genotypes; both showed a major peak on chromosome 1 (Fig. 5a, Supplementary Data 3 and 4) and a clear deviation of P-value distribution from the null (Supplementary Fig. 4a). Using imputed genotypes, the highest hit from imputed SNPs (chr1:161511936:T > C, P = 2.98 × 10−13) explained 2.85% of the total phenotypic variance (Fig. 5a). Under this peak in a 4-Mb region (158.5–162.5 Mb), there were 7138 variants within 22 genes. Linkage disequilibrium in this region was extensive, with 1050 variants in strong LD (r2 > 0.8) with the top hit, including the MC4R SNP (Fig. 5b). The highest hit was an intronic SNP in the gene CCBE1 (Fig. 5b). However, the extensive LD in this region makes it difficult to pinpoint a causative mutation by genetic data alone. Additional functional information and genetic data that break the LD are necessary to further fine-map causative genes and mutations. Nevertheless, the ability to identify the putative MC4R causative SNP as one of the top associated variants in a long stretch of high LD region clearly demonstrated the improvement of resolution using imputed genotypes. In our analysis, the MC4R SNP was initially removed and would otherwise be invisible without the imputation, as would be the case if the Illumina PorcineSNP60 chips were used.

a Manhattan plot of genome-wide association studies (GWAS) for backfat thickness. The gray (dark and light) points on the background are from GWAS using imputed genotypes, while the blue (light and dark) points are from GWAS using SNP chips. Genomic heritabilities calculated using array and imputed genotypes are indicated. The most significant SNPs from GWAS using imputed and array genotypes are indicated by circles and arrows. b Association within the 158.5–162.5 Mb region of chromosome 1, where the top hits in GWAS are located. Points indicate −log10(P-value) along the chromosome using imputed genotypes and SNPs where arrays also have genotypes are marked by crosses. The top SNPs from GWAS using imputed and array genotypes are marked by circles and arrows. r2 between the SNPs and the top SNP (chr1:161511936:T > C) is indicated by a gradient of blue color. Locations of genes are indicated in the box below the plot, where blue boxes and gene names with a left arrowhead (<) indicate genes transcribed on the reverse strand, and red boxes and gene names with a right arrowhead (>) indicate genes transcribed from the forward strand. Genes that are not marked do not have gene symbols. Gene locations are based on the Ensembl Release 98 annotation.

We next considered body length. We imputed genotypes from an Affymetrix 55K SNP chip (Wens55K) to a whole genome sequence using our imputation platform and performed GWAS in a population of 1694 Yorkshire boars (Fig. 6a). The trait has a moderately high heritability, as estimated using both array (h2 ~ 0.32) and imputed (h2 ~ 0.34) genotypes (Fig. 6a). Using GWAS (Supplementary Fig. S4b), we found a highly significant peak on chromosome 17 (Fig. 6a, Supplementary Data 5 and 6) where the lead variant was an intergenic SNP upstream of the BMP2 gene (chr17:15643342:C>T, P = 3.45 × 10−39). Remarkably, this variant explained 13.65% of the total phenotypic variance, and the homozygous C/C animals were, on average, 4.01 cm longer than the T/T homozygotes (Fig. 6b, c). BMP2 has been repeatedly shown to be associated with growth traits in pigs. A recent study implicated a regulatory variant upstream of the BMP2 gene and validated its functional impact using reporter genes26. This regulatory variant was the third most significant SNP under this peak in our analysis. Whether one or both of these potentially regulatory variants are the causative mutations remains to be determined. Given the strong association, high MAF of these SNPs, and less extensive LD in this region, it is unlikely that these regulatory variants were tagging protein-coding and less common variants in the BMP2 gene. In addition to the genetic support from this Yorkshire population, the body length increasing C allele was much more prevalent in Landrace than in other breeds. A hallmark of the Landrace breed is its long body size; thus, regulatory variation of the BMP2 gene may be a major contributor to the phenotypic differentiation between pig breeds. In contrast, although the SNP chip was able to broadly identify this region, the most significant SNP (chr17:15827832:T>G, P = 1.58 × 10−25) in an SNP chip-based GWAS was about 184 kb away from the lead SNP and explained a substantially smaller variance (8.22% versus 13.65%).

a Manhattan plot of genome-wide association studies (GWAS) for body length. The gray (dark and light) points on the background are from GWAS using imputed genotypes, while the blue (light and dark) points are from GWAS using SNP chips. Genomic heritabilities calculated using array and imputed genotypes are indicated. b Association within the 15.3–16.3 Mb region of chromosome 17, where the top hits in GWAS are located. Points indicate −log10(P-value) along the chromosome using imputed genotypes and SNPs where arrays also have genotypes are marked by crosses. The top SNPs from GWAS using imputed and array genotypes are marked by circles and arrows. r2 between the SNPs and the top SNP (chr17:15643342:C>T) are indicated by a gradient of blue color. Locations of genes are indicated in the box below the plot and according to the Ensembl Release 98 annotation. All three genes are colored in red and transcribed from the forward strand. The only gene with a symbol in this region is BMP2. c Scatter and box plots of body length (in cm) for the three genotypes of the chr17:15643342:C>T SNP. The lower and upper boundaries of the box are, respectively, 25% and 75% quantiles of the data, the midline median, and the whiskers minimum and maximum. d Allele frequencies of the chr17:15643342:C>T SNP in different breeds.

To enable the broad research community to efficiently utilize the resource developed in this study, we developed a public SWine IMputation (SWIM) web server (https://www.swimgeno.org and https://swim.scau.pigselection.com/swim), on which users can upload SNP chip genotypes and retrieve imputed genotypes. The user interface is extremely simple, which only requires users to upload the genotypes in gzipped ped/map format and leave their email addresses. Unlike other servers, such as PHARP, allele matching and flipping are performed on the server end, further simplifying the process on the user end. Imputation status can be monitored and results can be downloaded from a dynamic link without having to register an account. The server is set up to accommodate multiple users at the same time while limiting multiple jobs of the same user. Our tests indicated that a typical job with 2000 individuals and 50K SNP chip genotypes can be completed in approximately 12 h for all chromosomes.

We present here the development of the largest reference haplotype panel in pigs and an accompanying web server for the public to utilize this resource for genotype imputation. The high level of diversity and the large number of animals in the panel enabled us to achieve very high imputation accuracy with concordance rate, non-reference concordance rate, and r2 in excess of 95.84%, 88.26%, and 0.85, respectively, starting from 50K SNP arrays (Fig. 2). The accuracies were comparable to those obtained with medium density SNP arrays within pedigreed populations27. Given the high accuracy and easy access with no requirement for pedigree, we expect this public resource to vastly democratize sequence-level imputation in pigs and accelerate genetic discoveries. The SWIM server only supports SNP chip-based imputation at present. Low-coverage sequencing-based imputation is much more challenging to accommodate on a web server due to its requirement for massive computational resources. Nevertheless, users may implement their low-coverage sequencing-based imputation using the haplotype reference panel we share.

High-throughput genotyping arrays greatly simplified genotyping, and numerous new QTLs have been mapped by association mapping, typically within a breed and with hundreds to thousands of individuals6. However, while the resolution has improved with SNP arrays, causative genes and mutations remain extremely elusive, partly because SNP arrays prioritize assay feasibility, homogeneous spacing, and common SNPs5.

Our evaluations indicated that Shapeit4/Impute5 outperformed other software combinations, higher density of SNP chips led to higher imputation accuracy, and multi-breed haplotype reference panels maximizing sample size were preferred. Importantly, animals that were genetically closer to the haplotype reference panel could be imputed with higher accuracy. This further reinforces the importance of data sharing to increase representation in the haplotype reference panel.

As we have shown with the examples above, imputation is expected to greatly improve the resolution of gene mapping. Given the large number of existing genome-wide association studies in pigs6, we expect this resource to be highly utilized and impactful. Indeed, more than 130,000 genomes were imputed in the first year since the server became public, including a recent study that found SWIM imputed genomes to detect more significant SNPs compared to other platforms28. All existing studies using SNP arrays can be improved by a simple imputation followed by GWAS without additional data. Meta-analysis also becomes possible because a common SNP set can be obtained. Nonetheless, the resolution of genetic mapping depends not only on SNP density but also on experimental design and genetic structure in the mapping population. Sequence-level imputation does not necessarily identify causative mutations in one single step16. The availability of this resource will allow for suitable designs of mapping studies to achieve the highest possible resolution in specific circumstances and potentially nucleotide resolution.

We consolidated WGS data from multiple sources. A total of 1530 animals are first reported in this study using Illumina (n = 863) and BGI (n = 667) platforms with 150 bp paired-end reads. Among them, 610 Landrace, 413 Duroc, 391 Yorkshire, 18 Taiwanhei, and 17 Lichahei were from Wen's Food Group Co., Ltd. (Yunfu, Guangdong, China), 21 Dahuabai, 21 Lantanghei, 20 Guangdong Xiaoerhua, and 19 Yuedonghei from Guangdong Gene Bank of Livestock and Poultry (Guangzhou, Guangdong, China). Additionally, sequences for 729 animals were downloaded from the sequence read archive (SRA). A complete breakdown, including accession numbers, sample sizes, and average sequencing coverage, can be found in Supplementary Data 1 and 2.

We aligned sequence reads to the pig reference genome (Sscrofa11.1, a Duroc pig)29 using BWA-MEM-0.7.1730 and called variants (in GVCF format) using GATK-4.1.8.1 HaplotypeCaller31 after several post-alignment processing steps including duplicate removal using PicardTools-2.23.331, and base quality recalibration using GATK. A population VCF was generated by combining GVCFs across all samples. Variants with excessive heterozygosity ("ExcessHet > 54.69") were removed. Variant quality score recalibration (VQSR) on SNPs was performed with truth SNP sets compiled from commercial SNP arrays, including 50K, 60K, and 80K SNP chips (prior = 15.0) on the Illumina platform and the 660K (prior = 12.0), SowPro90 (prior = 15.0) SNP chips from the Affymetrix platform. SNPs were filtered with a truth sensitivity filter level at 99.0. Without a truth set of indels, we applied hard filtering on them by excluding indels with QD < 2.0, QUAL < 50.0, FS > 100.0, ReadPosRankSum < −20.0, as recommended by GATK's best practices. Additionally, we filtered out animals with a missing rate >0.20, heterozygosity >0.20, and retained bi-allelic sites with a missing rate <0.2 and mean sequencing depth between 5 and 500. Filtering was performed using a combination of VCFtools 0.1.1332 and BCFtools 1.1333 commands.

Linkage disequilibrium was computed using PopLDdecay34 on individuals within the same breed after removing close relatives (GRM > 0.5) and low-frequency variants (MAF < 0.05). To understand the genetic structure in the population, we retained variants with MAF > 0.05 and missing rate <0.1 and pruned SNPs with LD (r2 < 0.3, -indep-pairwise 50 10 0.3) using PLINK 1.935. Principal component analysis (PCA) was performed on the filtered list of 1,223,882 variants using GCTA 1.93.236 for all individuals. Ancestries were estimated using ADMIXTURE 1.337 on 185 individuals randomly selected according to breed representation in the dataset or at least four individuals per breed. The downsampling was necessary to properly visualize population structure.

We further filtered variants prior to phasing haplotypes in the reference population. Variants with missing rate >0.1 and MAF < 0.005 were removed. Additionally, variants with a Hardy–Weinberg equilibrium test P-value < 10−10 implemented separately in PLINK in all three of the Duroc, Landrace, and Yorkshire pigs were removed. Only autosomal variants were retained for imputation.

We extracted 100 Landrace pigs with the highest sequencing depth (17.42 X average sequencing depth, ranging from 14.98 to 63.11 X) and designated these individuals as the target population to evaluate imputation accuracy. To test the effect of breed composition of the reference population, we constructed four reference haplotype panels using different sets of individuals, including All (n = 2159): all individuals except the 100 Landraces; L (n = 550): Landrace pigs only; DLY (n = 550): 250 Landraces + 150 Durocs + 150 Yorkshires; and LO (n = 550): 250 Landraces + 300 randomly selected pigs other than Durocs and Yorkshires. Phasing was independently performed in these reference sets. In addition, we also tested imputation using the PHARP web server (http://alphaindex.zju.edu.cn/PHARP/index.php), which contains reference haplotypes constructed from 1006 individuals in the SRA.

We tested three combinations of software for phasing and imputation, including SHAPEIT 4.238 + IMPUTE5 1.1.539, Beagle 5.240 + Beagle 5.2, and Eagle 2.441 + Minimac 442. All software tools were run with default options and an uninformative linkage map (1 cM per 1 Mb), but the effective population size was set to 100. Imputed genotypes were called by the ones with the highest posterior genotype probability. However, users of the imputation web server also receive genotype probabilities.

We considered three commonly used metrics of imputation accuracy, concordance rate, non-reference concordance rate43, and r2. Concordance rate is defined as the proportion of individuals with imputed genotypes in concordance with observed genotypes. Non-reference concordance rate is similar to the concordance rate but is restricted to only individuals that are not homozygous for the reference allele. r2 is the squared Pearson correlation coefficient between observed and imputed genotypes. We measured concordance rates and r2 on a per SNP basis and averaged them over SNPs in MAF bins or across the whole genome.

To demonstrate the utility of imputation in genetic mapping, we collected phenotypes and genotypes for three populations of pigs, which were managed by three core breeding farms of Wen's Food Group Co., Ltd. (Yunfu, Guangdong, China), all under standard management practices. For backfat thickness, the phenotypes were collected on 3769 Duroc pigs from 2013 to 2018, and SNP genotyping was performed using the Geneseek GGP Porcine 50K SNP chip (Neogen, Lincoln, NE, USA). Backfat thickness was measured between the 10th and 11th ribs using an Aloka 500 V SSD B ultrasound (Corometrics Medical Systems, USA) when live weights of pigs reached about 100 kg (100 ± 5 kg). For body length, phenotypes from a total of 1694 Yorkshire boars were collected from 2012 to 2018, and SNP genotyping was performed using the Affymetrix PorcineWens55K SNP chip (Affymetrix, Santa Clara, CA, United States). Body length was measured from the base of the ear to the base of the tail in pigs at approximately 100 kg (100 ± 5 kg) body weight. All samples were collected according to the guidelines for the care and use of experimental animals approved by the Ministry of Agriculture and Rural Affairs of the People's Republic of China. The ethics committee of South China Agricultural University specifically approved the animal use in this study.

We used GCTA 1.92.1 to perform a mixed linear model (MLM) based association analysis. The following statistical model was used: \(y=\mu +{xb}+g+e\) (Equation 1), where y is the vector of the phenotypic values for all animals, \(\mu\) is the intercept, \(x\) is the design matrix coding genotypes and other incidences of fixed effects, \(b\) is the vector of fixed effects including SNP effect and additional covariates such as sex, pen, year-season effects depending on the traits, and \(g\) is the vector of polygenic random effects with covariance dictated by the genomic relationship matrix, and \(e\) is the vector of random residuals. We used SNPs on the GeneSeek GGP 50 K SNP chip (for backfat thickness) and Affymetrix Wens 55K SNP chip (for body length) to compute the genomic relationship matrix. We used a genome-wide significance threshold of P = 5 × 10−8 to declare significance. Variance explained by a single significant SNP was estimated by fitting a mixed linear model with the genomic relationship matrix determined by a single SNP.

All statistical analyses are performed using either software packages as described or in R 4.2.2. We supply all scripts, including those to generate figures in a GitHub (https://github.com/qgg-lab/swim-public) as well as a Zenodo repository44 (https://doi.org/10.5281/zenodo.7900470). The sample size for the entire SWIM haplotype reference panel is 2259, with subsets selected for the different designs to answer specific questions. Sample sizes for the backfat thickness and body length GWAS were 3769 and 1694, respectively.

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Raw sequence data for 512 animals have been deposited to SRA (PRJNA842867). Additional sequenced animals were proprietary properties of Wen's Food Group Co., Ltd. and Guangdong Gene Bank of Livestock and Poultry. They may be requested by contacting [email protected] and [email protected], respectively. Raw sequence data for a subset of the animals (n = 729) utilized in this study were downloaded from SRA (Supplementary Data 1 and 2). Imputation utilizing the full dataset is delivered as a web service (https://www.swimgeno.org and https://swim.scau.pigselection.com/swim) and is publicly available. Phased haplotypes from all publicly available individuals, including this study (n = 1241), are available as VCF files at https://quantgenet.msu.edu/swim/statistics.php. Source data underlying Figs. 1a, b, 2, 3, 4, and 6c are provided in Supplementary Data 7, 8, 9, 10, 11, and 12, respectively.

All computer codes, including all analyses performed in this study and codes for the SWIM web server, are available at https://github.com/qgg-lab/swim-public and at a Zenodo repository44 (https://doi.org/10.5281/zenodo.7900470).

Lunney, J. K. et al. Importance of the pig as a human biomedical model. Sci. Transl. Med. 13, eabd5758 (2021).

Article CAS PubMed Google Scholar

Groenen, M. A. M. et al. Analyses of pig genomes provide insight into porcine demography and evolution. Nature 491, 393–398 (2012).

Article CAS PubMed PubMed Central Google Scholar

Li, M. et al. Genomic analyses identify distinct patterns of selection in domesticated pigs and Tibetan wild boars. Nat. Genet. 45, 1431–1438 (2013).

Article CAS PubMed Google Scholar

Bosse, M. et al. Genomic analysis reveals selection for Asian genes in European pigs following human-mediated introgression. Nat. Commun. 5, 4392 (2014).

Article CAS PubMed Google Scholar

Ramos, A. M. et al. Design of a high density SNP genotyping assay in the pig using SNPs identified and characterized by next generation sequencing technology. PLoS ONE 4, e6524 (2009).

Article PubMed PubMed Central Google Scholar

Hu, Z.-L., Park, C. A. & Reecy, J. M. Building a livestock genetic and genomic information knowledgebase through integrative developments of Animal QTLdb and CorrDB. Nucleic Acids Res. 47, D701–D710 (2019).

Article CAS PubMed Google Scholar

Onteru, S. K. et al. Whole genome association studies of residual feed intake and related traits in the pig. PLoS ONE 8, e61756 (2013).

Article CAS PubMed PubMed Central Google Scholar

Sell-Kubiak, E. et al. Genome-wide association study reveals novel loci for litter size and its variability in a Large White pig population. BMC Genomics 16, 1049 (2015).

Article CAS PubMed PubMed Central Google Scholar

Boddicker, N. J. et al. Genome-wide association and genomic prediction for host response to porcine reproductive and respiratory syndrome virus infection. Genet. Sel. Evol. 46, 18 (2014).

Article PubMed PubMed Central Google Scholar

Marchini, J. & Howie, B. Genotype imputation for genome-wide association studies. Nat. Rev. Genet. 11, 499–511 (2010).

Article CAS PubMed Google Scholar

Das, S. et al. Next-generation genotype imputation service and methods. Nat. Genet.48, 1284–1287 (2016).

Article CAS PubMed PubMed Central Google Scholar

Daetwyler, H. D. et al. Whole-genome sequencing of 234 bulls facilitates mapping of monogenic and complex traits in cattle. Nat. Genet. 46, 858–865 (2014).

Article CAS PubMed Google Scholar

van den Berg, S. et al. Imputation to whole-genome sequence using multiple pig populations and its use in genome-wide association studies. Genet. Sel. Evol. 51, 2 (2019).

Article PubMed PubMed Central Google Scholar

Yang, W. et al. Animal-ImputeDB: a comprehensive database with multiple animal reference panels for genotype imputation. Nucleic Acids Res. 48, D659–D667 (2020).

Article PubMed Google Scholar

Wang, Z. et al. PHARP: a pig haplotype reference panel for genotype imputation. Sci. Rep. 12, 12645 (2022).

Article CAS PubMed PubMed Central Google Scholar

Yan, G. et al. An imputed whole-genome sequence-based GWAS approach pinpoints causal mutations for complex traits in a specific swine population. Sci. China Life Sci. 65, 781–794 (2022).

Article CAS PubMed Google Scholar

Ros-Freixedes, R. et al. Genomic prediction with whole-genome sequence data in intensely selected pig lines. Genet. Sel. Evol. 54, 65 (2022).

Article CAS PubMed PubMed Central Google Scholar

Rowan, T. N. et al. A multi-breed reference panel and additional rare variants maximize imputation accuracy in cattle. Genet. Sel. Evol. 51, 77 (2019).

Article CAS PubMed PubMed Central Google Scholar

Nezer, C. et al. An imprinted QTL with major effect on muscle mass and fat deposition maps to the IGF2 locus in pigs. Nat. Genet. 21, 155–156 (1999).

Article CAS PubMed Google Scholar

Van Laere, A.-S. et al. A regulatory mutation in IGF2 causes a major QTL effect on muscle growth in the pig. Nature 425, 832–836 (2003).

Article PubMed Google Scholar

Kim, K. S., Larsen, N., Short, T., Plastow, G. & Rothschild, M. F. A missense variant of the porcine melanocortin-4 receptor (MC4R) gene is associated with fatness, growth, and feed intake traits. Mamm. Genome 11, 131–135 (2000).

Article CAS PubMed Google Scholar

OVilo, C. et al. Test for positional candidate genes for body composition on pig chromosome 6. Genet. Sel. Evol. 34, 465–479 (2002).

Article PubMed Google Scholar

Gozalo-Marcilla, M. et al. Genetic architecture and major genes for backfat thickness in pig lines of diverse genetic backgrounds. Genet. Sel. Evol. 53, 76 (2021).

Article CAS PubMed PubMed Central Google Scholar

Farooqi, I. S. et al. Dominant and recessive inheritance of morbid obesity associated with melanocortin 4 receptor deficiency. J. Clin. Invest. 106, 271–279 (2000).

Article CAS PubMed PubMed Central Google Scholar

Krashes, M. J., Lowell, B. B. & Garfield, A. S. Melanocortin-4 receptor-regulated energy homeostasis. Nat. Neurosci. 19, 206–219 (2016).

Article CAS PubMed PubMed Central Google Scholar

Li, J. et al. Identification and validation of a regulatory mutation upstream of the BMP2 gene associated with carcass length in pigs. Genet. Sel. Evol. 53, 94 (2021).

Article CAS PubMed PubMed Central Google Scholar

Whalen, A. & Hickey, J. M. AlphaImpute2: fast and accurate pedigree and population based imputation for hundreds of thousands of individuals in livestock populations. Preprint at bioRxiv https://doi.org/10.1101/2020.09.16.299677 (2020).

Sun, J. et al. Genome-wide association study on reproductive traits using imputation-based whole-genome sequence data in Yorkshire pigs. Genes 14, 861 (2023).

Article CAS PubMed PubMed Central Google Scholar

Warr, A. et al. An improved pig reference genome sequence to enable pig genetics and genomics research. Gigascience 9, giaa051 (2020).

Article PubMed PubMed Central Google Scholar

Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760 (2009).

Article CAS PubMed PubMed Central Google Scholar

DePristo, M. A. et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43, 491–498 (2011).

Article CAS PubMed PubMed Central Google Scholar

Danecek, P. et al. The variant call format and VCFtools. Bioinformatics 27, 2156–2158 (2011).

Article CAS PubMed PubMed Central Google Scholar

Danecek, P. et al. Twelve years of SAMtools and BCFtools. Gigascience 10, giab008 (2021).

Article PubMed PubMed Central Google Scholar

Zhang, C., Dong, S.-S., Xu, J.-Y., He, W.-M. & Yang, T.-L. PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics 35, 1786–1788 (2019).

Article CAS PubMed Google Scholar

Chang, C. C. et al. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 4, 7 (2015).

Article PubMed PubMed Central Google Scholar

Yang, J., Lee, S. H., Goddard, M. E. & Visscher, P. M. GCTA: a tool for genome-wide complex trait analysis. Am. J. Hum. Genet. 88, 76–82 (2011).

Article CAS PubMed PubMed Central Google Scholar

Alexander, D. H., Novembre, J. & Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19, 1655–1664 (2009).

Article CAS PubMed PubMed Central Google Scholar

Delaneau, O., Zagury, J.-F., Robinson, M. R., Marchini, J. L. & Dermitzakis, E. T. Accurate, scalable and integrative haplotype estimation. Nat. Commun. 10, 5436 (2019).

Article PubMed PubMed Central Google Scholar

Rubinacci, S., Delaneau, O. & Marchini, J. Genotype imputation using the Positional Burrows Wheeler Transform. PLoS Genet. 16, e1009049 (2020).

Article CAS PubMed PubMed Central Google Scholar

Browning, B. L., Tian, X., Zhou, Y. & Browning, S. R. Fast two-stage phasing of large-scale sequence data. Am. J. Hum. Genet. 108, 1880–1890 (2021).

Article CAS PubMed PubMed Central Google Scholar

Loh, P.-R., Palamara, P. F. & Price, A. L. Fast and accurate long-range phasing in a UK Biobank cohort. Nat. Genet. 48, 811–816 (2016).

Article CAS PubMed PubMed Central Google Scholar

Howie, B., Fuchsberger, C., Stephens, M., Marchini, J. & Abecasis, G. R. Fast and accurate genotype imputation in genome-wide association studies through pre-phasing. Nat. Genet. 44, 955–959 (2012).

Article CAS PubMed PubMed Central Google Scholar

Li, J. H., Mazur, C. A., Berisa, T. & Pickrell, J. K. Low-pass sequencing increases the power of GWAS and decreases measurement error of polygenic risk scores compared to genotyping arrays. Genome Res. 31, 529–537 (2021).

Article PubMed PubMed Central Google Scholar

qgg-lab. qgg-lab/swim-public: swim-public-v1. Zenodo. https://doi.org/10.5281/zenodo.7900470 (2023).

Download references

This work is supported by a USDA-NIFA project (2021-67021-34149 to W.H., C.G., J.S., and R.Sc.), a USDA-NIFA Hatch project (MICL 02560 to W.H.), a Natural Science Foundation of China project (31972540 to J.Y.), a Natural Science Foundation of Guangdong Province project (2018B030313011 to Z.W.), and a Key Technologies R&D Program of Guangdong Province project (2022B0202090002 to Z.W.). The web server (https://www.swimgeno.org) is supported by the USDA Swine Genome Coordinator Fund (NRSP8).

Rodrigo Savegnago

Present address: Genus IntelliGen Technologies, De Forest, Wisconsin, USA

Suxu Tan

Present address: College of Life Sciences, Qingdao University, Qingdao, Shandong, China

College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, Guangdong, China

Rongrong Ding, Gengyuan Cai, Zhanwei Zhuang, Jie Wu, Ming Yang, Yibin Qiu, Donglin Ruan, Jianping Quan, Enqin Zheng, Huaqiang Yang, Zicong Li, Jie Yang & Zhenfang Wu

Department of Animal Science, Michigan State University, East Lansing, Michigan, USA

Rongrong Ding, Rodrigo Savegnago, Jinding Liu, Jianping Quan, Suxu Tan, Mohammed Bedhane, Juan Steibel, Cedric Gondro & Wen Huang

Yunfu Subcenter of Guangdong Laboratory for Lingnan Modern Agriculture, Yufu, Guandong, China

Rongrong Ding, Cheng Tan & Zhenfang Wu

Academy for Advanced Interdisciplinary Studies, Nanjing Agricultural University, Nanjing, Jiangsu, China

Jinding Liu

Institute for Cyber-Enabled Research, Michigan State University, East Lansing, Michigan, USA

Nanye Long

Guangdong Zhongxin Breeding Technology Co., Ltd, Guangzhou, Guangdong, China

Cheng Tan & Gengyuan Cai

Guangdong Provincial Key Laboratory of Agro-animal Genomics and Molecular Breeding, South China Agricultural University, Guangzhou, Guangdong, China

Zicong Li & Jie Yang

Division of Animal Sciences, University of Missouri, Columbia, Missouri, USA

Robert Schnabel

Department of Fisheries and Wildlife, Michigan State University, East Lansing, Michigan, USA

Juan Steibel

You can also search for this author in PubMed Google Scholar

W.H., Z.W., J.Y., and R.D.: conceptualization and design; R.D., R.Sa., N.L., and W.H.: developed and optimized pipeline; R.D., S.T., and M.B.: analyzed data; J.L. and W.H.: developed web server; R.Sc., C.T., G.C., Z.Z., J.W., M.Y., Y.Q., D.R., J.Q., E.Z., H.Y., Z.L., J.S., and C.G.: contributed tools and data; R.D. and W.H.: wrote the paper, with input from all authors.

Correspondence to Jie Yang, Wen Huang or Zhenfang Wu.

C.T. and G.C. are employees of the Guangdong Zhongxin Breeding Technology Co., Ltd. All other authors declare no competing interests.

Communications Biology thanks the anonymous reviewers for their contribution to the peer review of this work. Primary handling editor: George Inglis. A peer review file is available.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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 images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Ding, R., Savegnago, R., Liu, J. et al. The SWine IMputation (SWIM) haplotype reference panel enables nucleotide resolution genetic mapping in pigs. Commun Biol 6, 577 (2023). https://doi.org/10.1038/s42003-023-04933-9

Download citation

Received: 24 November 2022

Accepted: 12 May 2023

Published: 30 May 2023

DOI: https://doi.org/10.1038/s42003-023-04933-9

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.