BackgroundIn many traits, not only individual trait levels are under genetic control, but also the variation around that level. In other words, genotypes do not only differ in mean, but also in (residual) variation around the genotypic mean. New statistical methods facilitate gaining knowledge on the genetic architecture of complex traits such as phenotypic variability. Here we study litter size (total number born) and its variation in a Large White pig population using a Double Hierarchical Generalized Linear model, and perform a genome-wide association study using a Bayesian method.ResultsIn total, 10 significant single nucleotide polymorphisms (SNPs) were detected for total number born (TNB) and 9 SNPs for variability of TNB (varTNB). Those SNPs explained 0.83 % of genetic variance in TNB and 1.44 % in varTNB. The most significant SNP for TNB was detected on Sus scrofa chromosome (SSC) 11. A possible candidate gene for TNB is ENOX1, which is involved in cell growth and survival. On SSC7, two possible candidate genes for varTNB are located. The first gene is coding a swine heat shock protein 90 (HSPCB = Hsp90), which is a well-studied gene stabilizing morphological traits in Drosophila and Arabidopsis. The second gene is VEGFA, which is activated in angiogenesis and vasculogenesis in the fetus. Furthermore, the genetic correlation between additive genetic effects on TNB and on its variation was 0.49. This indicates that the current selection to increase TNB will also increase the varTNB.ConclusionsTo the best of our knowledge, this is the first study reporting SNPs associated with variation of a trait in pigs. Detected genomic regions associated with varTNB can be used in genomic selection to decrease varTNB, which is highly desirable to avoid very small or very large litters in pigs. However, the percentage of variance explained by those regions was small. The SNPs detected in this study can be used as indication for regions in the Sus scrofa genome involved in maintaining low variability of litter size, but further studies are needed to identify the causative loci.

In the era of genome-wide selection (GWS), genotype-by-environment (G×E) interactions can be studied using genomic information, thus enabling the estimation of SNP marker effects and the prediction of genomic estimated breeding values (GEBV) for young candidates for selection in different environments. Although G×E studies in pigs are scarce, the use of artificial insemination has enabled the distribution of genetic material from sires across multiple environments. Given the relevance of reproductive traits, such as the total number born (TNB) and the variation in environmental conditions encountered by commercial dams, understanding G×E interactions can be essential for choosing the best sires for different environments. The present work proposes a two-step reaction norm approach for G×E analysis using genomic information. The first step provided estimates of environmental effects (herd-year-season, HYS), and the second step provided estimates of the intercept and slope for the TNB across different HYS levels, obtained from the first step, using a random regression model. In both steps, pedigree ( A: ) and genomic ( G: ) relationship matrices were considered. The genetic parameters (variance components, h(2) and genetic correlations) were very similar when estimated using the A: and G: relationship matrices. The reaction norm graphs showed considerable differences in environmental sensitivity between sires, indicating a reranking of sires in terms of genetic merit across the HYS levels. Based on the G: matrix analysis, SNP by environment interactions were observed. For some SNP, the effects increased at increasing HYS levels, while for others, the effects decreased at increasing HYS levels or showed no changes between HYS levels. Cross-validation analysis demonstrated better performance of the genomic approach with respect to traditional pedigrees for both the G×E and standard models. The genomic reaction norm model resulted in an accuracy of GEBV for "juvenile" boars varying from 0.14 to 0.44 across different HYS levels, while the accuracy of the standard genomic prediction model, without reaction norms, varied from 0.09 to 0.28. These results show that it is important and feasible to consider G×E interactions in evaluations of sires using genomic prediction models and that genomic information can increase the accuracy of selection across environments.

BackgroundUse 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.MethodsPhenotypes 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.ResultsFor 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.ConclusionsImputation 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.Electronic supplementary materialThe online version of this article (10.1186/s12711-019-0445-y) contains supplementary material, which is available to authorized users.

Dominance has been suggested as one of the genetic mechanisms explaining heterosis. However, using traditional quantitative genetic methods it is difficult to obtain accurate estimates of dominance effects. With the availability of dense SNP (Single Nucleotide Polymorphism) panels, we now have new opportunities for the detection and use of dominance at individual loci. Thus, the aim of this study was to detect additive and dominance effects on number of teats (NT), specifically to investigate the importance of dominance in a Landrace-based population of pigs. In total, 1,550 animals, genotyped for 32,911 SNPs, were used in single SNP analysis. SNPs with a significant genetic effect were tested for their mode of gene action being additive, dominant or a combination. In total, 21 SNPs were associated with NT, located in three regions with additive (SSC6, 7 and 12) and one region with dominant effects (SSC4). Estimates of additive effects ranged from 0.24 to 0.29 teats. The dominance effect of the QTL located on SSC4 was negative (−0.26 teats). The additive variance of the four QTLs together explained 7.37% of the total phenotypic variance. The dominance variance of the four QTLs together explained 1.82% of the total phenotypic variance, which corresponds to one-fourth of the variance explained by additive effects. The results suggest that dominance effects play a relevant role in the genetic architecture of NT. The QTL region on SSC7 contains the most promising candidate gene: VRTN. This gene has been suggested to be related to the number of vertebrae, a trait correlated with NT.