# Genomic evaluation by including dominance effects and inbreeding depression for purebred and crossbred performance with an application in pigs

Tao Xiang, Ole Fredslund Christensen, Zulma Gladis Vitezica, Andres Legarra

► To cite this version:

Tao Xiang, Ole Fredslund Christensen, Zulma Gladis Vitezica, Andres Legarra. Genomic evaluation by including dominance effects and inbreeding depression for purebred and crossbred performance with an application in pigs. Genetics Selection Evolution, 2016, 48 (1), pp.92. 10.1186/s12711-016-0271-4. hal-01479227

HAL Id: hal-01479227

https://hal.science/hal-01479227

Submitted on 28 Feb 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Genomic evaluation by including dominance effects and inbreeding depression for purebred and crossbred performance with an application in pigs

Tao Xiang<fnref n="1" />,²* , Ole Fredslund Christensen<fnref n="1" />, Zulma Gladis Vitezica³ and Andres Legarra²

## Abstract

**Background:** Improved performance of crossbred animals is partly due to heterosis. One of the major genetic bases of heterosis is dominance, but it is seldom used in pedigree-based genetic evaluation of livestock. Recently, a trivariate genomic best linear unbiased prediction (GBLUP) model including dominance was developed, which can distinguish purebreds from crossbred animals explicitly. The objectives of this study were: (1) methodological, to show that inclusion of marker-based inbreeding accounts for directional dominance and inbreeding depression in purebred and crossbred animals, to revisit variance components of additive and dominance genetic effects using this model, and to develop marker-based estimators of genetic correlations between purebred and crossbred animals and of correlations of allele substitution effects between breeds; (2) to evaluate the impact of accounting for dominance effects and inbreeding depression on predictive ability for total number of piglets born (TNB) in a pig dataset composed of two purebred populations and their crossbreds. We also developed an equivalent model that makes the estimation of variance components tractable.

**Results:** For TNB in Danish Landrace and Yorkshire populations and their reciprocal crosses, the estimated proportions of dominance genetic variance to additive genetic variance ranged from 5 to 11%. Genetic correlations between breeding values for purebred and crossbred performances for TNB ranged from 0.79 to 0.95 for Landrace and from 0.43 to 0.54 for Yorkshire across models. The estimated correlation of allele substitution effects between Landrace and Yorkshire was low for purebred performances, but high for crossbred performances. Predictive ability for crossbred animals was similar with or without dominance. The inbreeding depression effect increased predictive ability and the estimated inbreeding depression parameter was more negative for Landrace than for Yorkshire animals and was in between for crossbred animals.

**Conclusions:** Methodological developments led to closed-form estimators of inbreeding depression, variance components and correlations that can be easily interpreted in a quantitative genetics context. Our results confirm that genetic correlations of breeding values between purebred and crossbred performances within breed are positive and moderate. Inclusion of dominance in the GBLUP model does not improve predictive ability for crossbred animals, whereas inclusion of inbreeding depression does.

<footnote>*Correspondence: Tao.Xiang@mbg.au.dk</footnote>

<footnote n="1">¹Department of Molecular Biology and Genetics, Center for Quantitative Genetics and Genomics, Aarhus University, 8830 Tjele, Denmark
Full list of author information is available at the end of the article</footnote>
## Background

Crossbreeding is primarily and intensively applied in meat production systems [1], especially for swine and poultry. Crossbreeding capitalizes on heterosis effects and complementarity between breeds, and results in an increased performance of crossbred animals compared to purebred animals [1]. In terminal crossbreeding systems, selection on purebred animals to maximize their crossbred performance is the ultimate goal [2, 3]. Due to the existence of genotype-by-environment interaction effects and non-additive genetic effects in combination with different allele frequencies in different breeds [3, 4], the genetic correlation of breeding values between purebred and crossbred performances ($r_{PC}$) is usually lower than 1 [1, 5], and therefore, purebred performance under nucleus conditions may not be an optimal predictor for crossbred performance in commercial animals [4, 6].

One of the major genetic bases of heterosis is dominance [7, 8]. At the level of gene action, dominance is due to interactions between alleles at the same locus [9]. In pedigree-based genetic evaluation, dominance is rarely included because large-scale datasets that comprise a high proportion of full sibs are required to obtain accurate estimates and because the computational complexity is high [10]. With the recent availability of single nucleotide polymorphism (SNP) information and the development of genomic selection, estimation of the dominance effects of SNPs has become more feasible [11, 12].

Genomic evaluation has been successfully used in purebred [13, 14] and crossbred populations [15–17]. However, these studies generally ignore the dominance effects. A number of studies have been carried out on genomic evaluation including dominance effects using either simulated [18] or real purebred data [9, 12].

Recently, several studies [19, 20] have tried to extend genomic evaluation including dominance effects from purebred performance to crossbred performance. However, they either used genomic information on purebred animals only [19] or applied a genomic model that assumed that all animals belong to a single population, and thus the variance components were estimated based only on the genotyped crossbred animals [20]. Nevertheless, combining purebred and crossbred information is essential to implement genetic evaluation for crossbred performance [1, 19]. Furthermore, because of genotype-by-environment interaction effects and different patterns of linkage disequilibrium (LD) between SNPs and quantitative trait loci (QTL), the effects of SNPs may be breed-specific [21]. To overcome these issues, a trivariate genomic best linear unbiased predictor (GBLUP) model that explicitly distinguishes between purebred and crossbred data and includes dominance was recently developed by Vitezica et al. [22]. This model allowed the

estimation of different, yet correlated, additive and dominance marker effects in crossbred and purebred individuals. However, the empirical predictive ability of the trivariate GBLUP model has not been evaluated yet.

Thus, the current study had the following objectives: (1) to show how genomic inbreeding can be meaningfully included in GBLUP, even for crossbred animals; (2) to estimate the variance components of additive and dominance genetic effects by using data on total number of piglets born (TNB) in two Danish purebred and one crossbred pig populations using the trivariate GBLUP model; (3) to show how to derive, from variance component estimates, estimated genetic correlations of breeding values between purebred and crossbred performances in each pure breed, and also correlations of allele substitution effects between the two pure breeds; and (4) to evaluate the impact of dominance effects from genomic information on genomic evaluation by comparing accuracies of estimated genomic values in different cross-validation scenarios.

## Methods

### Animals and genotypes

We begin this section with a short presentation of the data used in the study, with the aim of defining the notation for the methodological developments that follow. For this study, all datasets were provided by the Danish Pig Research Centre. Data from three Danish pig populations were analyzed simultaneously: Landrace (L), Yorkshire (Y) and their reciprocal crosses (LY). Only data on TNB data for the first parity of sows in the three populations were used. In total, there were 2126, 2218 and 5143 genotyped sows with own records on TNB for L, Y and LY, respectively. Instead of using original records, corrected phenotypic values of TNB were used as dependent variables for the trivariate GBLUP model, because the pre-correction for non-genetic effects, such as herd-year-season, month at farrowing, and service sire was more accurately achieved on a larger dataset (293,339 L, 180,112 Y, and 10,974 LY). Among the crossbred animals, 7407 LY had a Landrace sire and a Yorkshire dam, while 3567 LY had a Yorkshire dam and a Landrace sire; L and Y populations were from nucleus farms and LY from a commercial farm. The litters of purebred sows were both purebred and crossbred litters. The relationship between LY-L and LY-Y are comparable since, in both cases, parents of the F1 animals are in the purebred datasets; further details about the model used for the pre-correction are in [17]. All the purebred sows had first farrowing dates between 2003 and 2013, while the crossbred sows first farrowed between 2010 and 2013. Only five of these purebred L and Y sows were dams of the LY.

The pedigrees for both purebred and crossbred sows were available and all crossbred animals were traced
back to their purebred ancestors until 1994 by the DMU Trace program [23], as was done for the larger dataset used for pre-correction. Consequently, 8227 L, 9851 Y and 5143 LY individuals were in the pedigree. The dataset of pre-corrected TNB records for genotyped individuals is termed “full genomic dataset” throughout the whole paper, and it should not be confused with the larger dataset used to do the pre-correction.

For the “full genomic dataset”, purebred sows were genotyped with the Illumina PorcineSNP60 Genotyping BeadChip [24], while the crossbred sows were genotyped with a 8.5 K GGP-Porcine Low Density Illumina Bead SNP chip [25]. SNP quality controls (such as: call rate for individuals ≥80%; call rate for SNPs ≥90%; minor allele frequencies ≥0.01; etc.) were applied on the same dataset in a previous study [26], which provides more details. Then, for the crossbred individuals, imputation from low density to moderate density was done by using a joint reference panel of the two pure breeds [26] using the software Beagle version 3.3.2 [27] (imputation accuracies ≥95% in terms of correlation coefficients and ≥99% in terms of correct rates between imputed and true genotypes). Finally, 41,009 SNPs were available for all the recorded purebred and crossbred sows.

### Considering genomic inbreeding and heterosis

Inbreeding can be defined as the proportion of homozygous SNPs across all loci for each animal, as suggested by several authors (e.g., [28]). If there is directional dominance causing inbreeding depression [29], then inbreeding should be considered in the genetic evaluation models [30]. Otherwise, using pedigree or marker data, estimates of genetic parameters are inflated [30, 31]. In Vitezica et al. [22], genomic inbreeding was fitted as a covariate and, in the current study, we prove this reasoning by using a parametric genomic model, such as a GBLUP.

Theory and evidence of directional dominance (equivalently, inbreeding depression) suggest that dominance effects of genes (here associated to markers) should have a priori a positive value for traits that exhibit inbreeding depression or heterosis. If we call **d** the vector of dominance marker effects, the following prior distribution is plausible:

$$ (\mathbf{d}) \sim N(\mathbf{1}\mu_d, \mathbf{1}\sigma_d^2), $$

where $\mu_d$ is the overall mean of dominance effects, which should be positive if there is heterosis due to dominance. A typical model for genomic prediction is that in Toro and Varona [11]:

$$ \mathbf{y} = \mathbf{X}\beta + \mathbf{Z}\mathbf{a} + \mathbf{W}\mathbf{d} + \mathbf{e}, \quad (1) $$

where $\mathbf{y}$ contains phenotypic values; $\mathbf{X}\beta$ stands for fixed effects and random effects other than additive and dominance effects; $\mathbf{a}$ is the vector of “biological” additive SNP effects, $\mathbf{d}$ is the vector of “biological” dominance SNP effects for each of the markers; matrix $\mathbf{Z}$ has entries 1, 0, -1, for SNP genotypes AA, Aa and aa, respectively, while matrix $\mathbf{W}$ has entries 0, 1, 0 for SNP genotypes AA, Aa and aa, respectively. $\mathbf{e}$ is the vector of overall random residual effects.

Typically, genetic models require $\mathbf{a}$ and $\mathbf{d}$ to have zero means, which is not true for $\mathbf{d}$ when directional dominance exist. Defining $\mathbf{d}^* = \mathbf{d} - E(\mathbf{d})$, then $E(\mathbf{d}^*) = \mathbf{0}$, and Eq. (1) can be written as:

$$ \begin{aligned} \mathbf{y} &= \mathbf{X}\beta + \mathbf{Z}\mathbf{a} + \mathbf{W}(\mathbf{d}^* + \mathbf{E}(\mathbf{d})) + \mathbf{e} \\ &= \mathbf{X}\beta + \mathbf{Z}\mathbf{a} + \mathbf{W}\mathbf{d}^* + \mathbf{W}\mathbf{1}\mu_d + \mathbf{e}. \end{aligned} $$

The term $\mathbf{W}\mathbf{1}\mu_d$ is actually an average of dominance effects for each individual and is equal to $\mathbf{h}\mu_d$, where $\mathbf{h} = \mathbf{W}\mathbf{1}$ contains the row-sums of $\mathbf{W}$, i.e. individual heterozygosities (it should be noted that $\mathbf{W}$ has a value of 1 at heterozygous loci for an individual). Inbreeding coefficients $\mathbf{f}$ can be calculated as:

$$ \mathbf{f} = \mathbf{1} - \mathbf{h}/N, $$

where $N$ is the number of SNPs. Then, the prior means $\mathbf{h}\mu_d$ can be rewritten as:

$$ \mathbf{h}\mu_d = (\mathbf{1} - \mathbf{f})N\mu_d = \mathbf{1}N\mu_d + \mathbf{f}(-N\mu_d). $$

The term $\mathbf{1}N\mu_d$ is confounded with the overall mean of the model ($\mu$), while the term $\mathbf{f}(-N\mu_d)$ models the inbreeding depression and $b = (-N\mu_d)$ is the inbreeding depression parameter summed over the SNPs, which has to be estimated. Thus, the linear model including genomic inbreeding is, finally:

$$ \mathbf{y} = \mathbf{X}\beta + \mathbf{f}b + \mathbf{Z}\mathbf{a} + \mathbf{W}\mathbf{d}^* + \mathbf{e}. $$

Thus, we have proven why fitting overall homozygosity for the individual as a measure of inbreeding depression accounts for directional dominance.

### Estimating genetic (co)variances of markers with additive and dominance effects

A trivariate model based on “biological” (genotypic) additive and dominance effects of SNPs [22, 32], and including genomic inbreeding as above, was applied considering TNB as a different trait in each population:

$$ \begin{aligned} \mathbf{y}_L &= \mathbf{1}\mu_L + \mathbf{f}_L b_L + \mathbf{Z}_L \mathbf{a}_L + \mathbf{W}_L \mathbf{d}_L + \mathbf{e}_L, \\ \mathbf{y}_Y &= \mathbf{1}\mu_Y + \mathbf{f}_Y b_Y + \mathbf{Z}_Y \mathbf{a}_Y + \mathbf{W}_Y \mathbf{d}_Y + \mathbf{e}_Y, \\ \mathbf{y}_{LY} &= \mathbf{1}\mu_{LY} + \mathbf{f}_Y b_Y + \mathbf{Z}_{LY} \mathbf{a}_{LY} + \mathbf{W}_{LY} \mathbf{d}_{LY} + \mathbf{e}_{LY}, \end{aligned} \quad (2) $$

where $\mathbf{y}_L$, $\mathbf{y}_Y$ and $\mathbf{y}_{LY}$ contain corrected phenotypic values for purebred L, purebred Y and crossbred LY sows,
respectively; $\mu_L, \mu_Y$ and $\mu_{LY}$ are the respective means; $\mathbf{a}_L, \mathbf{a}_Y$ and $\mathbf{a}_{LY}$ are the “biological” additive SNP effects and $\mathbf{d}_L, \mathbf{d}_Y$ and $\mathbf{d}_{LY}$ are the “biological” dominance SNP effects for each of the SNPs for L, Y and LY, respectively; matrices $\mathbf{Z}$ and $\mathbf{W}$ are as above; $\mathbf{f}_L b_L, \mathbf{f}_Y b_Y$ and $\mathbf{f}_{LY} b_{LY}$ model the inbreeding depression for L, Y and LY populations; $\mathbf{e}_L, \mathbf{e}_Y$ and $\mathbf{e}_{LY}$ are the overall random residual effects.

Note that “biological” is used here to refer to the genotypic additive and dominance values of the SNPs, to distinguish them from the traditional treatment of quantitative genetics in terms of “statistical” effects (breeding values and dominance deviations) [32].

The above equations can be reformulated to genotypic values of individuals instead of SNPs, in order to be compatible with the classical GBLUP model and animal breeding software, such as BLUPF90 [33] and DMU [34]:

$$ 
\begin{aligned}
\mathbf{y}_L &= \mathbf{1}\mu_L + \mathbf{f}_L b_L + \mathbf{u}_L + \mathbf{v}_L + \mathbf{e}_L, \\
\mathbf{y}_Y &= \mathbf{1}\mu_Y + \mathbf{f}_Y b_Y + \mathbf{u}_Y + \mathbf{v}_Y + \mathbf{e}_Y, \\
\mathbf{y}_{LY} &= \mathbf{1}\mu_{LY} + \mathbf{f}_{LY} b_{LY} + \mathbf{u}_{LY} + \mathbf{v}_{LY} + \mathbf{e}_{LY}.
\end{aligned} 
\quad (3) $$

Note that $\mathbf{u}$ and $\mathbf{v}$ are vectors of genotypic additive and dominance effects and therefore cannot be directly compared to breeding values and dominance deviations in the pedigree-based genetic evaluation. In addition, $\mathbf{f}$ is a vector of genomic inbreeding coefficients and $b$ is a population-specific inbreeding depression parameter per unit of genomic inbreeding, respectively. Note that there is potentially inbreeding depression at the level of the crossbred animals, although, first, the numeric values of the vector $\mathbf{f}$ should be smaller since crossbred animals have a higher level of heterozygosity, and second, the estimates of the inbreeding depression parameters (b) do not need to be identical across the three populations, which thus gives considerable flexibility.

In terms of the genotypic additive effects $\mathbf{u}$, the variances within each breed are:

$$ 
\begin{aligned}
\mathrm{Var}(\mathbf{u}_L) &= \mathrm{var}(\mathbf{Z}_L \mathbf{a}_L) = \mathbf{Z}_L \mathbf{Z}'_L \sigma_{a_L}^2, \\
\mathrm{Var}(\mathbf{u}_Y) &= \mathrm{var}(\mathbf{Z}_Y \mathbf{a}_Y) = \mathbf{Z}_Y \mathbf{Z}'_Y \sigma_{a_Y}^2, \\
\mathrm{Var}(\mathbf{u}_{LY}) &= \mathrm{var}(\mathbf{Z}_{LY} \mathbf{a}_{LY}) = \mathbf{Z}_{LY} \mathbf{Z}'_{LY} \sigma_{a_{LY}}^2,
\end{aligned} 
$$

where $\sigma_{a_L}^2, \sigma_{a_Y}^2$ and $\sigma_{a_{LY}}^2$ are the additive variances of SNP effects in breeds L, Y and LY, respectively. The covariances between the genotypic additive effects $\mathbf{u}$ are:

$$ 
\mathrm{Cov} \begin{pmatrix} \mathbf{u}_L \\ \mathbf{u}_Y \\ \mathbf{u}_{LY} \end{pmatrix} = \begin{pmatrix} \mathbf{Z}_L \mathbf{Z}'_L \sigma_{a_L}^2 & \mathbf{Z}_L \mathbf{Z}'_Y \sigma_{a_{LY}}^2 & \mathbf{Z}_L \mathbf{Z}'_{LY} \sigma_{a_{LY,LY}} \\ \mathbf{Z}_Y \mathbf{Z}'_L \sigma_{a_{LY,LY}} & \mathbf{Z}_Y \mathbf{Z}'_Y \sigma_{a_Y}^2 & \mathbf{Z}_Y \mathbf{Z}'_{LY} \sigma_{a_{LY,LY}} \\ \mathbf{Z}_{LY} \mathbf{Z}'_L \sigma_{a_{LY,LY}} & \mathbf{Z}_{LY} \mathbf{Z}'_Y \sigma_{a_{LY,Y}} & \mathbf{Z}_{LY} \mathbf{Z}'_{LY} \sigma_{a_{LY}}^2 \end{pmatrix}, \quad (4) $$

where $\sigma_{a_{LY},Y}, \sigma_{a_{LY,LY}}$ and $\sigma_{a_{LY,Y,LY}}$ are the additive covariances of SNP effects between populations L and Y, populations

L and LY, and populations Y and LY, respectively. Analogous structures exist for dominance genotypic effects:

$$ 
\mathrm{Cov} \begin{pmatrix} \mathbf{v}_L \\ \mathbf{v}_Y \\ \mathbf{v}_{LY} \end{pmatrix} = \begin{pmatrix} \mathbf{W}_L \mathbf{W}'_L \sigma_{d_L}^2 & \mathbf{W}_L \mathbf{W}'_Y \sigma_{d_{LY}} & \mathbf{W}_L \mathbf{W}'_{LY} \sigma_{d_{LY,LY}} \\ \mathbf{W}_Y \mathbf{W}'_L \sigma_{d_{LY}} & \mathbf{W}_Y \mathbf{W}'_Y \sigma_{d_Y}^2 & \mathbf{W}_Y \mathbf{W}'_{LY} \sigma_{d_{Y,LY}} \\ \mathbf{W}_{LY} \mathbf{W}'_L \sigma_{d_{LY,LY}} & \mathbf{W}_{LY} \mathbf{W}'_Y \sigma_{d_{Y,LY}} & \mathbf{W}_{LY} \mathbf{W}'_{LY} \sigma_{d_Y}^2 \end{pmatrix}. $$

### Estimation of marker-based variance components using an equivalent model

The variance components $\sigma_{a_L}^2, \sigma_{a_Y}^2, \sigma_{a_{LY}}^2$ and $\sigma_{a_{LY,LY}}, \sigma_{a_{LY,Y}}$ cannot be estimated by regular methods or software (i.e. REML or Gibbs sampling) because they cannot be factorized out from Eq. (4). To fit such a multivariate structure, we used an equivalent model. Additional effects need to be defined, even if they are of no interest per se. For instance, the vectors of hypothetical genotypic additive effects of the genotypes of the L breed on the scale of breed Y ($\mathbf{u}_{L,Y}$) and LY ($\mathbf{u}_{LY,LY}$) have variance-covariance matrices $\mathbf{Z}_L \mathbf{Z}'_L \sigma_{a_L}^2$ and $\mathbf{Z}_L \mathbf{Z}'_L \sigma_{a_{LY}}^2$, respectively. Thus, as a whole, the genetic variance and covariance structure for the genotypic additive effects $\mathbf{u}$ are:

$$ 
\begin{aligned}
\mathrm{Var}(\mathbf{u}) &= \mathrm{var} \begin{bmatrix} \mathbf{u}_L \\ \mathbf{u}_{L,Y} \\ \mathbf{u}_{L,LY} \\ \mathbf{u}_{Y,L} \\ \mathbf{u}_{Y,Y} \\ \mathbf{u}_{LY,L} \\ \mathbf{u}_{LY,Y} \\ \mathbf{u}_{LY} \end{bmatrix} = \mathrm{var} \begin{bmatrix} \mathbf{Z}_L \mathbf{a}_L \\ \mathbf{Z}_L \mathbf{a}_Y \\ \mathbf{Z}_L \mathbf{a}_{LY} \\ \mathbf{Z}_Y \mathbf{a}_L \\ \mathbf{Z}_Y \mathbf{a}_Y \\ \mathbf{Z}_Y \mathbf{a}_{LY} \\ \mathbf{Z}_{LY} \mathbf{a}_L \\ \mathbf{Z}_{LY} \mathbf{a}_Y \\ \mathbf{Z}_{LY} \mathbf{a}_{LY} \end{bmatrix} \\
&= \begin{bmatrix} \mathbf{Z}_L \mathbf{Z}'_L & \mathbf{Z}_L \mathbf{Z}'_Y & \mathbf{Z}_L \mathbf{Z}'_{LY} \\ \mathbf{Z}_Y \mathbf{Z}'_L & \mathbf{Z}_Y \mathbf{Z}'_Y & \mathbf{Z}_Y \mathbf{Z}'_{LY} \\ \mathbf{Z}_{LY} \mathbf{Z}'_L & \mathbf{Z}_{LY} \mathbf{Z}'_Y & \mathbf{Z}_{LY} \mathbf{Z}'_{LY} \end{bmatrix} \otimes \begin{bmatrix} \sigma_{a_L}^2 & \sigma_{a_{LY}} & \sigma_{a_{LY,LY}} \\ \sigma_{a_{LY,L}} & \sigma_{a_Y}^2 & \sigma_{a_{Y,LY}} \\ \sigma_{a_{LY,Y}} & \sigma_{a_{LY,Y}} & \sigma_{a_{LY}}^2 \end{bmatrix} \\
&= \mathbf{Z}\mathbf{Z}' \otimes \begin{bmatrix} \sigma_{a_L}^2 & \sigma_{a_{LY}} & \sigma_{a_{LY,LY}} \\ \sigma_{a_{LY,L}} & \sigma_{a_Y}^2 & \sigma_{a_{Y,LY}} \\ \sigma_{a_{LY,Y}} & \sigma_{a_{LY,Y}} & \sigma_{a_{LY}}^2 \end{bmatrix},
\end{aligned}
$$

where matrix $\mathbf{Z}$ contains elements 1, 0, -1 for the three genotypes, and is defined across the three breeds,

$$ \mathbf{Z} = \begin{bmatrix} \mathbf{Z}_L \\ \mathbf{Z}_Y \\ \mathbf{Z}_{LY} \end{bmatrix}. $$

To construct a relationship matrix similar to the classical G-matrix of GBLUP [35], Vitezica et al. [22] introduced a normalized genomic relationship matrix $\mathbf{G} = \frac{\mathbf{ZZ}'}{\{\mathrm{tr}[\mathbf{ZZ}']\}/n}$, where n is the number of animals across the three populations and the division by $\{\mathrm{tr}[\mathbf{ZZ}']\}/n$ scales the matrix such that the average of the diagonal elements equals 1. This alters the variances across genotypic additive effects $\mathbf{u}$ in the following way:
$$
\begin{align}
\mathrm{Var}(\mathbf{u}) &= \mathrm{var} \left[ \begin{bmatrix} \mathbf{u}_L \\ \mathbf{u}_{L,Y} \\ \mathbf{u}_{L,LY} \\ \mathbf{u}_{Y,L} \\ \mathbf{u}_Y \\ \mathbf{u}_{Y,LY} \\ \mathbf{u}_{LY,L} \\ \mathbf{u}_{LY,Y} \\ \mathbf{u}_{LY} \end{bmatrix} \right] = \mathbf{ZZ}' \otimes \begin{bmatrix} \sigma_{a_L}^2 & \sigma_{a_L,Y} & \sigma_{a_L,LY} \\ \sigma_{a_{Y,L}} & \sigma_{a_Y}^2 & \sigma_{a_Y,LY} \\ \sigma_{a_{LY,L}} & \sigma_{a_{LY,Y}} & \sigma_{a_{LY}}^2 \end{bmatrix} \nonumber \\
&= \left( \mathbf{G} \times \frac{\{\mathrm{tr}[\mathbf{ZZ}']\}}{n} \right) \otimes \begin{bmatrix} \sigma_{a_L}^2 & \sigma_{a_L,Y} & \sigma_{a_L,LY} \\ \sigma_{a_{Y,L}} & \sigma_{a_Y}^2 & \sigma_{a_Y,LY} \\ \sigma_{a_{LY,L}} & \sigma_{a_{LY,Y}} & \sigma_{a_{LY}}^2 \end{bmatrix} \nonumber \\
&= \mathbf{G} \otimes \begin{bmatrix} \sigma_{A_L}^2 & \sigma_{A_L,A_Y} & \sigma_{A_L,A_{LY}} \\ \sigma_{A_Y,A_L} & \sigma_{A_Y}^2 & \sigma_{A_Y,A_{LY}} \\ \sigma_{A_L,A_{LY}} & \sigma_{A_Y,A_{LY}} & \sigma_{A_{LY}}^2 \end{bmatrix} = \mathbf{G} \otimes \mathbf{G}_0, \tag{5}
\end{align}
$$

where $G_0$ are variance components associated to the genotypic additive effects $\mathbf{u}$. This structure (a Kronecker product) is compatible with animal breeding software for BLUP and REML and the variance-covariance component $G_0$ can be estimated in a straightforward manner. Then, the (co)variances of additive genotypic effects of SNPs across populations can be obtained as:

$$
\begin{align}
\begin{bmatrix}
\sigma_{a_L}^2 & \sigma_{a_L,Y} & \sigma_{a_L,LY} \\
\sigma_{a_{Y,L}} & \sigma_{a_Y}^2 & \sigma_{a_Y,LY} \\
\sigma_{a_{LY,L}} & \sigma_{a_{LY,Y}} & \sigma_{a_{LY}}^2
\end{bmatrix}
&= \mathbf{G}_0 / \{\mathrm{tr}[\mathbf{ZZ}']\} / n \nonumber \\
&= \begin{bmatrix}
\sigma_{A_L}^2 & \sigma_{A_L,A_Y} & \sigma_{A_L,A_{LY}} \\
\sigma_{A_Y,A_L} & \sigma_{A_Y}^2 & \sigma_{A_Y,A_{LY}} \\
\sigma_{A_L,A_{LY}} & \sigma_{A_Y,A_{LY}} & \sigma_{A_{LY}}^2
\end{bmatrix} / \{\mathrm{tr}[\mathbf{ZZ}']\} / n. \tag{6}
\end{align}
$$

The variances across genotypic dominance effects $\mathbf{v}$ are altered in a similar way:

$$
\begin{align}
\mathrm{Var}(\mathbf{v}) &= \mathrm{var} \left[ \begin{bmatrix} \mathbf{v}_L \\ \mathbf{v}_{L,Y} \\ \mathbf{v}_{L,LY} \\ \mathbf{v}_{Y,L} \\ \mathbf{v}_Y \\ \mathbf{v}_{Y,LY} \\ \mathbf{v}_{LY,L} \\ \mathbf{v}_{LY,Y} \\ \mathbf{v}_{LY} \end{bmatrix} \right] \nonumber \\
&= \mathbf{D} \otimes \begin{bmatrix} \sigma_{D_L}^2 & \sigma_{D_L,D_Y} & \sigma_{D_L,D_{LY}} \\ \sigma_{D_Y,D_L} & \sigma_{D_Y}^2 & \sigma_{D_Y,D_{LY}} \\ \sigma_{D_L,D_{LY}} & \sigma_{D_Y,D_{LY}} & \sigma_{D_{LY}}^2 \end{bmatrix} = \mathbf{D} \otimes \mathbf{D}_0, \tag{7}
\end{align}
$$

where $\mathbf{D}_0$ contains variances and covariances associated to the genotypic dominance effects $\mathbf{v}$ and $\mathbf{D} = \frac{\mathbf{WW}'}{\{\mathrm{tr}[\mathbf{WW}']\}/n}$, where the matrix $\mathbf{W}$ contains elements 0, 1, 0 for the three genotypes, and is defined

across the three breeds ($\mathbf{W} = \begin{bmatrix} \mathbf{W}_L \\ \mathbf{W}_Y \\ \mathbf{W}_{LY} \end{bmatrix}$) and $\mathbf{W}' = [\mathbf{W}'_L \mathbf{W}'_Y \mathbf{W}'_{LY}]$. Then, the (co)variances of dominance genotypic effects of SNPs are:

$$
\begin{align}
\begin{bmatrix}
\sigma_{d_L}^2 & \sigma_{d_L,Y} & \sigma_{d_L,LY} \\
\sigma_{d_{Y,L}} & \sigma_{d_Y}^2 & \sigma_{d_Y,LY} \\
\sigma_{d_{LY,L}} & \sigma_{d_{LY,Y}} & \sigma_{d_{LY}}^2
\end{bmatrix}
&= \mathbf{D}_0 / \{\mathrm{tr}[\mathbf{WW}']\} / n \tag{8} \\
&= \begin{bmatrix}
\sigma_{D_L}^2 & \sigma_{D_L,D_Y} & \sigma_{D_L,D_{LY}} \\
\sigma_{D_Y,D_L} & \sigma_{D_Y}^2 & \sigma_{D_Y,D_{LY}} \\
\sigma_{D_L,D_{LY}} & \sigma_{D_Y,D_{LY}} & \sigma_{D_{LY}}^2
\end{bmatrix} / \{\mathrm{tr}[\mathbf{WW}']\} / n \nonumber
\end{align}
$$

This approach, which is an extension of Vitezica et al. [22], makes it possible to estimate (co)variances of genotypic effects of SNPs in purebred and crossbred populations under a genomic model with additive and non-additive (dominance) inheritance.

Matrices $\mathbf{Z}$ and $\mathbf{W}$, their crossproducts and the inverses of $\mathbf{G}$ and $\mathbf{D}$ were built using own programs. Genetic parameters were estimated by using average information REML with software iremlf90 [33]. Standard errors on functions of genetic parameters (i.e. standard errors on correlations) were estimated from the average information matrix using the REML-MVN method of Houle and Meyer [36].

### Additive and dominance variances in purebred and crossbred populations

The additive and dominance (co)variances of genotypic effects of SNPs, either within breed or between breeds, were calculated using Eqs. (6) and (8), respectively. Using these calculated additive and dominance (co)variances of SNPs across all the SNPs, the corresponding traditional, individual-based genetic parameters can be obtained as follows. The genetic parameters obtained are directly comparable to pedigree-based estimates [32].

Consider the allele substitution effect $\alpha = a + (q - p)d$. According to [32], the additive genetic variances for purebred performance (mating animals in the same breed) for breed L ($\sigma_{AP_L}^2$) and Y ($\sigma_{AP_Y}^2$) are:

$$
\begin{align*}
\sigma_{AP_L}^2 &= \sum (2p_i^L q_i^L) \sigma_{a_L}^2 + \sum (2p_i^L q_i^L (q_i^L - p_i^L))^2 \sigma_{d_L}^2, \\
\sigma_{AP_Y}^2 &= \sum (2p_i^Y q_i^Y) \sigma_{a_Y}^2 + \sum (2p_i^Y q_i^Y (q_i^Y - p_i^Y))^2 \sigma_{d_Y}^2,
\end{align*}
$$

where $\sigma_a^2$ and $\sigma_d^2$ are the variances of additive and dominance genotypic effects of SNPs in either breed L or Y; $p_i$ and $q_i$ are allele frequencies for SNP $i$; indices $L$ and $Y$ denote the breeds Landrace and Yorkshire, respectively. For crossbred performance of say, Landrace, the allele substitution effect is $\alpha_{AC_L} = a_{AC_L} + (q^Y - p^Y)d_{AC_L}$. Thus, the additive genetic variances within purebred L and Y for crossbred performance (due to gametes from
the L or Y individuals in the crossbred population) are equal to:

$$ 
\begin{aligned}
\sigma_{AC_L}^2 &= \sum (2p_i^L q_i^L) \sigma_{a_{LY}}^2 + \sum (2p_i^L q_i^L (q_i^Y - p_i^Y)^2) \sigma_{d_{LY}}^2, \\
\sigma_{AC_Y}^2 &= \sum (2p_i^Y q_i^Y) \sigma_{a_{LY}}^2 + \sum (2p_i^Y q_i^Y (q_i^L - p_i^L)^2) \sigma_{d_{LY}}^2,
\end{aligned}
$$

where the $\sigma_{AC_L}^2$ represents the additive genetic variance of animals in breed L when mated to animals in breed Y; the $\sigma_{AC_Y}^2$ represents the additive genetic variance of animals in breed Y when mated to animals in breed L; and $\sigma_{a_{LY}}^2$ and $\sigma_{d_{LY}}^2$ are the variances of additive and dominance genotypic effects of SNPs in the crossbred LY population, respectively. The additive genetic variance for animals in the crossbred LY population ($\sigma_{AC_{LY}}^2$) is the sum of the additive genetic variance of Landrace alleles and that of Yorkshire alleles in the crossbred animals [22] as follows:

$$ \sigma_{AC_{LY}}^2 = \frac{1}{2}\sigma_{AC_L}^2 + \frac{1}{2}\sigma_{AC_Y}^2. $$

Note that this variance is not the additive genetic variance of the crossbred animals acting as reproducers (i.e., creating an F2) [37].

The additive genetic covariances between purebred and crossbred performances within breeds L ($\sigma_{AP_L,AC_L}$) and Y ($\sigma_{AP_Y,AC_Y}$) are:

$$ 
\begin{aligned}
\sigma_{AP_L,AC_L} &= \sum (2p_i^L q_i^L) \sigma_{a_{L,LY}} + \sum (2p_i^L q_i^L (q_i^L - p_i^L) (q_i^Y - p_i^Y)) \sigma_{d_{L,LY}}, \\
\sigma_{AP_Y,AC_Y} &= \sum (2p_i^Y q_i^Y) \sigma_{a_{Y,LY}} + \sum (2p_i^Y q_i^Y (q_i^Y - p_i^Y) (q_i^L - p_i^L)) \sigma_{d_{Y,LY}},
\end{aligned}
$$

where $\sigma_{a_{L,LY}}$ and $\sigma_{d_{L,LY}}$ are the covariances of SNP effects between purebred L and crossbred LY populations for additive and dominance, respectively; $\sigma_{a_{Y,LY}}$ and $\sigma_{d_{Y,LY}}$ are the covariances of SNP effects between purebred Y and crossbred LY populations for additive and dominance, respectively.

Therefore, the genetic correlations of breeding values between purebred and crossbred performances within L ($r_{pCL}$) and Y ($r_{pCY}$) are: $r_{pCL} = \frac{\sigma_{AP_L,AC_L}}{\sqrt{\sigma_{AP_L}^2 \sigma_{AC_L}^2}}$ and $r_{pCY} = \frac{\sigma_{AP_Y,AC_Y}}{\sqrt{\sigma_{AP_Y}^2 \sigma_{AC_Y}^2}}$.

According to [22], the dominance genetic variances within purebred populations L and Y are $\sigma_{D_L}^2 = \sum (2p_i^L q_i^L)^2 \sigma_{d_L}^2$ and $\sigma_{D_Y}^2 = \sum (2p_i^Y q_i^Y)^2 \sigma_{d_Y}^2$, respectively. The dominance genetic variance in crossbred LY animals is $\sigma_{D_{LY}}^2 = \sum (4p_i^L q_i^L p_i^Y q_i^Y) \sigma_{d_Y}^2$.

The broad sense heritabilities for purebred performance ($H_p^2$) were calculated as the ratio of total genetic variances for purebred performance ($\sigma_{AP}^2 + \sigma_{D}^2$) to phenotypic variances ($\sigma_{AP}^2 + \sigma_{D}^2 + \sigma_{e}^2$).

## Correlations of allele substitution effects between two breeds

The breeding value of an individual includes the allele substitution effects of all genes and the allele frequencies. For purebred performance, the allele substitution effects of one locus for breed L and Y are:

$$ 
\begin{aligned}
\alpha_L &= a_L + (q_i^L - p_i^L) d_L, \\
\alpha_Y &= a_Y + (q_i^Y - p_i^Y) d_Y,
\end{aligned}
$$

where $a$ is the additive effect and $d$ is the dominance effect for each SNP; $p_i$ and $q_i$ are allele frequencies for SNP $i$, with superscripts denoting breeds L or Y. In the case of purely additive gene action, the covariance between $\alpha_L$ and $\alpha_Y$ is $\sigma_{a_{L,Y}}$, which can be interpreted as a genetic correlation among populations [38–40]. Then, the covariance between the allele substitution effects of one locus is:

$$ 
\begin{aligned}
\operatorname{cov}(\alpha_L, \alpha_Y) &= \operatorname{cov}(a_L + (q_i^L - p_i^L)d_L, a_Y + (q_i^Y - p_i^Y)d_Y) \\
&= \operatorname{cov}(a_L, a_Y) + (q_i^L - p_i^L)(q_i^Y - p_i^Y)\operatorname{cov}(d_L, d_Y) \\
&= \sigma_{a_{L,Y}} + (q_i^L - p_i^L)(q_i^Y - p_i^Y)\sigma_{d_{L,Y}},
\end{aligned}
$$

where $\sigma_{a_{L,Y}}$ and $\sigma_{d_{L,Y}}$ are the additive and dominance covariances of SNP effects between breeds L and Y for additive and dominance, respectively. If we assume that

SNP effects (both additive and dominance) are independent across loci, then the covariance between the allele substitution effects across all $n$ loci is:

$$ \operatorname{cov}(\alpha_L, \alpha_Y) = \sigma_{a_{L,Y}} = \sigma_{a_{L,Y}} + \frac{1}{n} \sum ((q_i^L - p_i^L)(q_i^Y - p_i^Y)) \sigma_{d_{L,Y}} $$

Also, the variances of allele substitution effects across all $n$ loci for breeds L and Y are:

$$ 
\begin{aligned}
\operatorname{var}(\alpha_L) &= \sigma_{\alpha_L}^2 = \sigma_{a_L}^2 + \frac{1}{n} \sum ((q_i^L - p_i^L)^2) \sigma_{d_L}^2, \\
\operatorname{var}(\alpha_Y) &= \sigma_{\alpha_Y}^2 = \sigma_{a_Y}^2 + \frac{1}{n} \sum ((q_i^Y - p_i^Y)^2) \sigma_{d_Y}^2,
\end{aligned}
$$

where $\sigma_a^2$ and $\sigma_d^2$ are the additive and dominance variance of SNPs. Then, the correlation of allele substitution effects for purebred performance between populations L and Y is $r_{\alpha P_L, \alpha P_Y} = \frac{\sigma_{\alpha_{L,Y}}}{\sigma_{\alpha_L} \sigma_{\alpha_Y}}$. If there is no dominance variation, the $r_{\alpha P_L, \alpha P_Y}$ relates to additive genetic variances as $r_{\alpha P_L, \alpha P_Y} = \frac{\sigma_{a_{L,Y}}}{\sigma_{a_L} \sigma_{a_Y}}$.

The correlation of allele substitution effects for crossbred performance between populations L and Y is similar
to that for purebred performance, but the allele frequencies are swapped, as:

$$r_{\alpha C_L, \alpha C_Y} = \frac{\sigma_{\alpha_{L in LY}, Y in LY}}{\sigma_{\alpha_{L in LY}} \sigma_{\alpha_Y in LY}} = \frac{\sigma_{a_{LY}}^2 + \frac{1}{n} \sum ((q_i^L - p_i^L)(q_i^Y - p_i^Y)) \sigma_{d_{LY}}^2}{\sqrt{\sigma_{a_{LY}}^2 + \frac{1}{n} \sum ((q_i^Y - p_i^Y)^2) \sigma_{d_{LY}}^2} \sqrt{\sigma_{a_{LY}}^2 + \frac{1}{n} \sum ((q_i^L - p_i^L)^2) \sigma_{d_{LY}}^2}},$$

where $\sigma_{a_{LY}}^2$ and $\sigma_{d_{LY}}^2$ are the additive and dominance variance of SNPs in the crossbred LY population. If there is no dominance variation, the $r_{\alpha C_L, \alpha C_Y}$ is equal to 1, by assumption in the model.

### Scenarios

Variance components, genetic correlations of breeding values between purebred and crossbred performances ($r_{PC}$) within each pure breed and correlations of allele substitution effects for purebred ($r_{\alpha P_L, \alpha P_Y}$) and crossbred ($r_{\alpha C_L, \alpha C_Y}$) performance between two pure breeds were first investigated using the full genomic dataset. To explore the effects of using genomic information and the inclusion of dominance deviation on the genetic evaluation of crossbred performance in the trivariate model, three different scenarios were compared.

### Nogen

The statistical model was a trivariate BLUP model, similar to Eq. (3), but the dominance deviation was excluded. Instead of using a genomic relationship matrix, a single relationship matrix A was constructed across the three breeds, assuming that they form a single population. Thus, the genetic (co)variances of additive genetic effects u were:

$$Var(\mathbf{u}) = \mathbf{A} \otimes \begin{bmatrix} \sigma_{A_L}^2 & \sigma_{A_L A_Y} & \sigma_{A_L A_{LY}} \\ \sigma_{A_Y A_L} & \sigma_{A_Y}^2 & \sigma_{A_Y A_{LY}} \\ \sigma_{A_L A_{LY}} & \sigma_{A_Y A_{LY}} & \sigma_{A_{LY}}^2 \end{bmatrix} = \mathbf{A} \otimes \mathbf{A}_0,$$

where $\mathbf{A}_0$ were variance components associated to genetic additive effects and not the genotypic additive effects in Eq. (5). Pedigree-based inbreeding depression was also included in the model. The pedigree-based inbreeding coefficients were calculated as in [41] using the software inbupgf90 [33].

### Gen_AM

The statistical model was similar to Eq. (3), but without dominance deviations. Genomic information was used to construct the additive genomic relationship matrix.

### Gen_ADM

The statistical model includes additive and dominance effects as in Eq. (3). Genomic information was used to

construct the additive and dominance genomic relationship matrices.

To explore the impact of genomic information and dominance effects on genomic evaluation for crossbred performance, the full genomic dataset was split into training and validation populations and the predictive ability for crossbred animals in the validation population was investigated in different scenarios. The farrowing date of January 1, 2013 was used as the cut-off date to divide recorded purebred and crossbred sows into training and validation populations. As a result, 6769 sows (1270 L, 1405 Y and 4094 LY) were included in the training population, while the remaining 2716 sows (854 L, 813Y and 1049 LY) were included in the validation population. Predictive ability of crossbreds was measured as the correlations $cor(y_c, \hat{y})$ in the validation population for each scenario, where $y_c$ is the corrected phenotypic records of TNB for crossbred animals; $\hat{y}$ is the predicted corrected observations of TNB for crossbred animals and is equal to the sum of the estimated population mean ($\hat{\mu}$), inbreeding ($f\hat{b}$) and genotypic values ($\hat{g}$); the genotypic value $\hat{g}$ was calculated as the sum of additive and dominance genetic effects in the scenario Gen_ADM. In the other two scenarios, the genotypic value $\hat{g}$ only included the additive genetic effect. Hotelling-Williams t test at a confidence level of 5% was applied to evaluate the significance of the differences in validation correlations in each scenario. Furthermore, to detect the possible biases in the predictions, the regression coefficients of $y_c$ on $\hat{y}$ were explored. Note that no bias implies that a regression coefficient equals 1. In addition, to measure the uncertainty associated with the predictions, 1000 bootstrap samples [42] was applied to estimate the means and standard errors.

For comparison, the predictive ability of crossbred animals was also investigated in a model without inbreeding depression effects, for all three scenarios. The predictive ability was measured as the correlation $cor(y_c, \hat{y})$, where $\hat{y}$ is the sum of the estimated population mean ($\hat{\mu}$) and genotypic value ($\hat{g}$).
<table><caption>Table 1 Variance components of additive and dominance genetic effects for purebred and crossbred animals</caption><thead><tr><th>Scenario</th><th>Breed</th><th>σ<sup>2</sup><sub>AP</sub></th><th>σ<sub>AP,AC</sub></th><th>σ<sup>2</sup><sub>AC</sub></th><th>σ<sup>2</sup><sub>D</sub></th><th>σ<sup>2</sup><sub>e</sub></th><th>σ<sup>2</sup><sub>ACLY</sub></th><th>σ<sup>2</sup><sub>DLY</sub></th><th>σ<sup>2</sup><sub>eLY</sub></th></tr></thead><tbody><tr><td rowspan="2">Nogen</td><td>L</td><td>0.99 (0.31)</td><td>0.17 (0.07)</td><td>0.05 (0.02)</td><td>–</td><td>10.82 (0.43)</td><td>0.05 (0.02)</td><td>–</td><td>7.35 (0.15)</td></tr><tr><td>Y</td><td>1.07 (0.33)</td><td>0.15 (0.07)</td><td>0.05 (0.02)</td><td>–</td><td>8.96 (0.38)</td><td>–</td><td>–</td><td>–</td></tr><tr><td rowspan="2">Gen_AM</td><td>L</td><td>0.87 (0.22)</td><td>0.47 (0.10)</td><td>0.28 (0.07)</td><td>–</td><td>10.89 (0.38)</td><td>0.28 (0.07)</td><td>–</td><td>7.11 (0.15)</td></tr><tr><td>Y</td><td>0.55 (0.20)</td><td>0.17 (0.10)</td><td>0.28 (0.07)</td><td>–</td><td>9.42 (0.33)</td><td>–</td><td>–</td><td>–</td></tr><tr><td rowspan="2">Gen_ADM</td><td>L</td><td>0.86 (0.21)</td><td>0.46 (0.10)</td><td>0.28 (0.06)</td><td>0.04 (0.03)</td><td>10.86 (0.38)</td><td>0.28 (0.06)</td><td>0.02 (0.01)</td><td>7.11 (0.15)</td></tr><tr><td>Y</td><td>0.54 (0.18)</td><td>0.17 (0.09)</td><td>0.28 (0.06)</td><td>0.06 (0.05)</td><td>9.35 (0.33)</td><td>–</td><td>–</td><td>–</td></tr></tbody></table>

Numbers in brackets are the standard errors of the corresponding parameters

$\sigma_{AP}^2$ is the additive genetic variance for purebred performance; $\sigma_{AP,AC}$ is the additive genetic covariance between purebred and crossbred performance; $\sigma_{AC}^2$ is the additive genetic variance for crossbred performance; $\sigma_{D}^2$ is the dominance genetic variance for either purebred animals; $\sigma_{e}^2$ is the residual variance for purebred animals; $\sigma_{ACLY}^2$ is the additive genetic variance for the F1 crossbred animals LY; $\sigma_{DLY}^2$ is the dominance genetic variance for the F1 crossbred animals LY; and $\sigma_{eLY}^2$ is the residual variance for the F1 crossbred animals LY

L Landrace, Y Yorkshire breeds

# Results

## Variance components, heritabilities and correlations

Table 1 shows the estimates of variance components for additive genetic effects for purebred ($\sigma_{AP}^2$) and crossbred ($\sigma_{AC}^2$) performance in different scenarios, and dominance variations ($\sigma_{D}^2$) in the Gen_ADM scenario. For all scenarios, the additive genetic variances for purebred performance ($\sigma_{AP}^2$) were larger than those for their crossbred performance ($\sigma_{AC}^2$). Estimated variance components in the scenarios Gen_AM and Gen_ADM were very close, but different from those obtained in scenarios without using genomic information. In general, estimates had large standard errors in all scenarios, but no obvious differences in standard errors were detected between different scenarios. Residual variance for purebred animals ($\sigma_{e}^2$) was larger than for crossbred animals ($\sigma_{eLY}^2$) in each scenario. For the scenario Gen_ADM, the ratios of dominance genetic variance to additive genetic variance ranged from 5 to 11% for both purebred and crossbred populations.

The broad sense heritabilities for purebred and crossbred animals, genetic correlations between breeding values for purebred and crossbred performances within pure breeds and correlations of allele substitution effects across the two breeds are in Table 2. In different scenarios, the heritabilities of purebred performance ($H_p^2$) ranged from 0.07 (0.03) to 0.08 (0.03) and from 0.06 (0.03) to 0.10 (0.03) for breeds L and Y, respectively. Standard errors of $H_p^2$ were almost consistent across scenarios. Estimated genetic correlations of breeding values between purebred and crossbred performances ($r_{PC}$) increased from 0.76 (0.20) (Nogen) to 0.95 (0.06) (Gen_AM) for breed L and from 0.43 (0.22) (Gen_ADM) to 0.54 (0.30) (Nogen) for breed Y. The $r_{PC}$ was higher for breed L than for breed Y in all scenarios, but the standard errors of $r_{PC}$ were always higher for breed Y than for breed L. With genomic information, the correlations

<table><caption>Table 2 Heritabilities and genetic correlations between breeding values for purebred and crossbred performances</caption><thead><tr><th>Scenario</th><th>Breed</th><th>r<sub>PC</sub></th><th>H<sub>p</sub><sup>2</sup></th></tr></thead><tbody><tr><td rowspan="2">Nogen</td><td>L</td><td>0.76 (0.20)</td><td>0.08 (0.03)</td></tr><tr><td>Y</td><td>0.54 (0.30)</td><td>0.10 (0.03)</td></tr><tr><td rowspan="2">Gen_AM</td><td>L</td><td>0.95 (0.06)</td><td>0.07 (0.03)</td></tr><tr><td>Y</td><td>0.44 (0.20)</td><td>0.06 (0.03)</td></tr><tr><td rowspan="2">Gen_ADM</td><td>L</td><td>0.93 (0.05)</td><td>0.08 (0.03)</td></tr><tr><td>Y</td><td>0.43 (0.22)</td><td>0.06 (0.03)</td></tr></tbody></table>

Numbers between brackets are the standard errors of the corresponding parameters

$r_{PC}$ is the genetic correlation of breeding values between purebred and crossbred performances within the Landrace or Yorkshire breeds; $H_p^2$ is the broad sense heritability for purebred performance for the Landrace and Yorkshire breeds in different scenarios

L Landrace, Y Yorkshire

of allele substitution effects between purebred ($r_{\alpha P_L, \alpha P_Y}$) and crossbred ($r_{\alpha C_L, \alpha C_Y}$) performance between breeds L and Y were estimated, as shown in Table 3. For purebred performance, $r_{\alpha P_L, \alpha P_Y}$ was equal to 0.14 and 0.19 in Gen_AM and Gen_ADM, respectively. However, the standard errors were large, around 0.2 in both scenarios. For crossbred performance, $r_{\alpha C_L, \alpha C_Y}$ was equal to 0.98 in Gen_ADM. This high correlation is a byproduct of assuming that additive biological effects in crossbred animals are the same regardless of the Yorkshire or Landrace origin of the allele. However, the same allele has potentially different effects in the respective Landrace or Yorkshire genetic backgrounds, and the difference is modeled through the correlations, hence the low values of $r_{\alpha P_L, \alpha P_Y}$. Without including the dominance effects in the model Gen_AM, $r_{\alpha C_L, \alpha C_Y}$ was equal to 1 by definition.
<table><caption>**Table 3** Correlations of allele substitution effects for purebred and crossbred performance between Landrace and Yorkshire breeds</caption><thead><tr><th>Scenario</th><th>r<sub>αP<sub>L</sub>αP<sub>Y</sub></sub></th><th>r<sub>αC<sub>L</sub>αC<sub>Y</sub></sub></th></tr></thead><tbody><tr><td>Nogen</td><td>–</td><td>–</td></tr><tr><td>Gen_AM</td><td>0.14 (0.22)</td><td>1</td></tr><tr><td>Gen_ADM</td><td>0.19 (0.24)</td><td>0.98 (0.02)</td></tr></tbody></table>

Numbers between brackets are the standard errors of the corresponding parameters

$r_{αP_L,αP_Y}$ is the correlation of allele substitution effects for purebred performance between the Landrace and Yorkshire breeds; $r_{αC_L,αC_Y}$ is the correlation of allele substitution effects for crossbred performance between the Landrace and Yorkshire breeds. For Gen_AM, $r_{αC_L,αC_Y}$ is equal to 1 by definition

## Predictive abilities

Predictive abilities for crossbred pigs in the validation population are in Table 4. The correlation between the corrected phenotypic values and the predicted observations for TNB ($cor(y_c, \hat{y})$) ranged from 0.010 in the scenario Nogen to 0.056 in scenarios Gen_AM and Gen_ADM. Standard errors of $cor(y_c, \hat{y})$ based on 1000 bootstrap samples were equal to 0.03 across all scenarios. No significant differences in predictive ability between scenarios were detected by the Hotelling-Williams t test at the confidence level of 5%.

The regression coefficients of corrected phenotypic values on the predicted corrected observations for TNB are in the second row of Table 4. Regression coefficients were smaller than 1 for the three scenarios. Among these scenarios, regression coefficients for scenarios with genomic information (Gen_AM and Gen_ADM) were slightly closer to 1 than that for the pedigree-based scenario (Nogen). Except for the Nogen scenario, standard errors of regression coefficients were around 0.39. For the Nogen scenario, the standard error was around 5 times larger than that for other scenarios. Overall, there was no clear trend towards a scenario with less bias.

<table><caption>**Table 4** Predictive ability for crossbred animals in the validation population</caption><thead><tr><th></th><th>Nogen</th><th>Gen_AM</th><th>Gen_ADM</th></tr></thead><tbody><tr><td>cor(y_c, ŷ)<sup>a</sup></td><td>0.010 (0.031)</td><td>0.056 (0.031)</td><td>0.056 (0.031)</td></tr><tr><td>Regression coefficient<sup>b</sup></td><td>0.703 (2.218)</td><td>0.736 (0.386)</td><td>0.730 (0.385)</td></tr></tbody></table>

Numbers between brackets are the standard errors of the corresponding parameters

For comparison, predictive abilities $cor(y_c, ŷ)$ for crossbred pigs in the validation population for the models without the inbreeding depression effect were equal to -0.08 in scenario Nogen, 0.045 in scenario Gen_AM and 0.046 in scenario Gen_ADM. In all cases, these are lower than the predictive abilities in Table 4, and these differences are statistically significant according to the Hotelling-Williams t test.

## Inbreeding depression

Marker-based and pedigree-based inbreeding coefficient (f) for each population and their estimated corresponding inbreeding depression parameters (b) in the different scenarios are in Table 5. Marker-based inbreeding coefficients were almost identical for breeds L and Y, but they were larger than those for LY, which was expected because crossbred animals have a higher level of heterozygosity than purebred animals. However, according to the pedigree-based inbreeding coefficients, the Landrace population was slightly more inbred than the Yorkshire population. In terms of inbreeding depression parameters (b), they were all negative (thus, genomic inbreeding has detrimental effects for TNB even in crossbred animals) but not of the same magnitude across the three populations. Note that for the scenario Nogen, b was estimated based on the pedigree-based inbreeding coefficients. As a whole, breed L had the most negative b, while breed Y had the least negative b, regardless of the scenario. Thus, TNB was more negatively affected by inbreeding in breed L than in breed Y and population LY.

<table><caption>**Table 5** Marker-based and pedigree-based inbreeding coefficients f and estimated inbreeding depression parameter b (piglets per 100% of inbreeding) in different scenarios for each breed</caption><thead><tr><th></th><th>L</th><th>Y</th><th>LY</th></tr></thead><tbody><tr><td>Marker-based inbreeding coefficient <sup>a</sup></td><td>0.695 (0.019)</td><td>0.698 (0.020)</td><td>0.565 (0.012)</td></tr><tr><td>Pedigree-based inbreeding coefficient <sup>b</sup></td><td>0.111 (0.032)</td><td>0.078 (0.031)</td><td>0</td></tr><tr><td>Nogen (b)</td><td>-4.821</td><td>-3.561</td><td>0</td></tr><tr><td>Gen_AM (b)</td><td>-9.656</td><td>-1.924</td><td>-5.122</td></tr><tr><td>Gen_ADM (b)</td><td>-9.731</td><td>-1.878</td><td>-5.055</td></tr></tbody></table>

The inbreeding coefficient is the mean inbreeding coefficient across individuals within each breed

Numbers between brackets are the standard deviations of the mean inbreeding coefficient

For Nogen, the inbreeding depression parameter b is the regression of phenotype on pedigree-based inbreeding. For Gen_AM and Gen_ADM, the inbreeding depression parameter b is the regression of phenotype on marker-based inbreeding

<footnote>a Predictive ability ($cor(y_c, ŷ)$) is given by the correlation coefficient between the corrected phenotypes ($y_c$) and their predictions ($ŷ$) for total number of piglets born (TNB) in crossbred animals</footnote>

<footnote>b Regression coefficient of the corrected phenotypes ($y_c$) on the predicted observations ($ŷ$) in crossbred animals</footnote>

<footnote>a Calculated as the proportion of homozygous loci per individual</footnote>

<footnote>b Calculated as in Meuwissen and Luo [41]</footnote>
Discussion

This study extended the trivariate GBLUP model of Vite-
zica et al. [22] in order to obtain (co)variances of effects
of SNPs, genetic correlations of breeding values between
purebred and crossbred performances and correlations
of allele substitution effects under dominance. We also
evaluated this model using different scenarios for the
genetic evaluation of crossbred performance in Danish
purebred and crossbred pigs. Scenarios that included or
not genomic information were studied to estimate the
genetic correlations of breeding values between purebred
and crossbred performances. To our knowledge, this is
the first study to report correlations of allele substitution
effects between two breeds in the presence of dominance
effects. The results show that the Vitezica model [22] is a
tool that can be used for the genomic evaluation of cross-
bred performance in genotyped animals. In this study, for
TNB, models with dominance deviations did not improve
the genomic evaluation of crossbred performance with
regard to both predictive ability and unbiasedness, but
the inclusion of an inbreeding depression effect in the
models significantly improved predictive ability.

Phenotypic variances were larger for purebred animals (11.76 for breed L and 9.99 for breed Y) than for crossbred animals (7.30 for LY). This could be the reason why the estimated additive genetic variances for purebred performance (σ<sup>2</sup><sub>AP</sub>) were larger than those for crossbred performance (σ<sup>2</sup><sub>AC</sub>). However, compared to results in a previous study that used a much larger Danish purebred and crossbred dataset [17], both estimated additive genetic variances and phenotypic variances in the current study were smaller, which is due to three reasons. (1) The dataset in the current study was a genotyped subset of the population used in the previous study. Purebred genotyped individuals were pre-selected and their performances were more homogeneous than that of the whole population. The preselection process resulted in a loss of about 15% of the purebred phenotypic variation. However, the genotyped crossbred animals were an almost random sample of the whole population and there was only a small loss of about 5% of phenotypic variation for crossbred animals. (2) The phenotypic values for TNB in the current study were pre-corrected for fixed and non-genetic random effects. This pre-correction led to a loss of about 11 and 17% of phenotypic variation for purebreds and crossbreds, respectively. (3) During the pre-correction, some genetic variation may have been allocated to other random effects (e.g. service boar effects), in particular because TNB is a lowly heritable trait.

The estimated heritabilities of TNB for purebred per-
formance ($H_p^2$) were slightly lower than those previously
reported (0.11 and 0.09 for breeds L and Y, respectively)

[17, 22, 43]. Large standard errors of $H_p^2$ implied that the current dataset was not large enough. The consistent standard errors across scenarios indicated that even when genomic information was included, the uncertainty of $H_p^2$ did not decrease. Taking the standard errors into account, the estimated $H_p^2$ across scenarios were not very different. Compared to the results of [17], the lower $H_p^2$ found in the current study was due to the sharp decrease in additive genetic variances ($\sigma_{AP}^2$).

The ratios of estimated dominance genetic variances to
additive genetic variances in the current study (5 to 11%)
were generally a little smaller than in other studies on
TNB. Vitezica et al. [22] reported that this ratio was equal
to about 20% for litter size in both purebred and cross-
bred lines by using the same trivariate GBLUP model.
Esfandyari et al. [19] stated that, by using purebred
genomic information in a univariate Bayesian mixture
model at the SNP level, the ratio between dominance
variance and additive variance for TNB was equal to 15
and 18% for breeds L and Y, respectively. Based on pedi-
gree information, Misztal et al. [10] reported a ratio that
reached about 25% for number of piglets born alive in
a Yorkshire population. However, there are some stud-
ies that did report smaller ratios than those reported
here. For instance, Hidalgo [20] reported that, based on
genotyped crossbred animals, the dominance variance
for TNB accounted for nearly zero of the total genetic
variance and concluded that TNB was not affected by
dominance effects in the Dutch Landrace and Yorkshire
populations. For other traits or species, different ratios of
dominance genetic variance to additive genetic variance
were also reported. For average daily gain in Duroc pigs,
Su et al. [9] estimated a ratio of 15%, but their results
were based on genotypic variance components and can-
not be directly compared to genetic variance components
[32]. For average daily weight gain in Yorkshire and Lan-
drace pigs, Lopes et al. [44] reported ratios of 13.8 and
28%, respectively by including genomic information.
For Fleckvieh cattle, Ertl et al. [12] calculated ratios that
ranged from 3.4% for stature to 69% for protein yield by
using a univariate SNP-BLUP model. Overall, these dif-
ferent ratios of dominance genetic variance to additive
genetic variance may reflect differences in the traits ana-
lyzed and in the type of information used for the estima-
tion [9], and also uncertainty in the estimates.

The genetic correlation of breeding values between
purebred and crossbred performances ($r_{PC}$) is a key
parameter in crossbreeding schemes [2]. In the current
study, the estimated $r_{PC}$ was in line with results reviewed
by Wei et al. [3]. Lutaaya et al. [5] also reported $r_{PC}$
that ranged from 0.32 to 1. Such differences in $r_{PC}$ may
reflect differences in the extent of GxE interactions and
the distance across breeds. In our study, estimated $r_{PC}$
did not vary dramatically across the scenarios, when the standard errors were taken into account. These standard errors were very large, which indicated that the amount of available information was too small to ensure accurate $r_{PC}$ estimates. Across scenarios, standard errors of $r_{PC}$ decreased when genomic information was included, which indicates that including genomic information may reduce the uncertainty of the estimations. $r_{PC}$ was larger for breed L than for breed Y, which was in agreement with a previous study [17] and may be due to the data structure. Among the 5143 crossbred animals, the number of Yorkshire sires (N = 1125) was much smaller than that of Landrace sires (N = 4018). Such a different amount of information affects the accuracy of the estimates, and thus the standard error of $r_{PC}$ was larger for breed Y than for breed L (see Table 2). However, compared to the results reported in [17], the $r_{PC}$ for breed L increased by about 10% while that for breed Y did not change much. Both pre-correction of data and the genotyped subset of original data used may play a role in the differences observed between the current and previous results [17]. In the previous study, a single-step method, which can use pedigree information and genomic information simultaneously, was used. In this study, the use of only phenotypic records on genotyped individuals affected the accuracy of estimates. Our results confirmed the moderate value of the $r_{PC}$ for TNB in breeds L and Y.

To our knowledge, this is the first time that correlations of allele substitution effects for both purebred ($r_{\alpha P_L, \alpha P_Y}$) and crossbred ($r_{\alpha C_L, \alpha C_Y}$) performance between two breeds in the presence of dominance variation are estimated. In genomic selection, SNPs are assumed to be in LD with QTL along the whole genome [45]. The correlation of allele substitution effects between breeds measures the degree of average similarities between SNP effects assuming that the QTL effects are the same in breeds 1 and 2 [38–40]. In practice, the correlation of allele substitution effects between two breeds can be interpreted as indicating “how consistent the SNP substitution effects are across two breeds”. For purebred performance, the estimated SNP substitution effects were based on the within-breed allele frequencies. A high $r_{\alpha P_L, \alpha P_Y}$ correlation means that the estimated SNP substitution effects based on allele frequencies from breed L can be used for breed Y and vice versa. However, $r_{\alpha P_L, \alpha P_Y}$ was not significantly different from 0 in the current study, which demonstrates that SNP effects estimated from a reference population that consists of one pure breed (e.g. Landrace) cannot be readily applied to the other breed (e.g. Yorkshire). This was in agreement with the findings of [46] who reported that prediction based on an across-population reference panel was worse than within-population prediction. In other species, estimated correlations

of allele substitution effects between breeds based on models without dominance, oscillate between 0 and 0.8, and are trait-dependent [38, 47]. For crossbred performance, an $r_{\alpha C_L, \alpha C_Y}$ close to 1 was found in the current study, which indicated that the allele substitution effects based on the allele frequencies from the opposite breeds were very similar for the L and Y breeds. In practice, this suggests that SNP substitution effects that are estimated based on a reference population consisting of crossbred animals can be used to estimate crossbred breeding values for both breeds L and Y.

It was expected that genomic evaluations obtained by including dominance deviations in the model would be improved, especially when records of crossbred animals were included [9]. However, our results showed that inclusion of dominance deviations did not increase the predictive ability for crossbreds. This result was in line with conclusions in [9, 12, 20], but was opposite to those in [18, 19, 48, 49]. Theoretically, estimating dominance genetic effects should be useful because ignoring them will result in less accurate estimates of allele substitution effects and consequently less accurate estimated breeding values in genomic prediction [11]. However, regarding the additive genetic variance, estimates were nearly the same in scenarios *Gen_AM* and *Gen_ADM*, which demonstrated that the additive variances were already well captured by the additive model. Thus, the accuracy of the estimated additive genetic effects was not affected when dominance effects were included in the model [12]. Moreover, a simulation study at the level of the gene action showed that when all gene actions were purely additive, including dominance in addition to the additive effects in the model was not advantageous compared to using an additive model. Hidalgo [20] showed that TNB was not affected by dominance in the Dutch crossbred population. In the current study, we also observed similar results, and dominance variation accounted for a small proportion of the total genetic variation (4 to 10%). The lack of change in predictive ability also indicated the difficulty of distinguishing dominance genetic effects from additive genetic effects [9], but it confirmed a previous simulation study that concluded that the use of a dominance model did not negatively affect genomic evaluation even if the trait was purely additive [18].

Scenarios in which genomic information was included (*Gen_AM* and *Gen_ADM*) showed higher predictive abilities than the pedigree-based scenario (*Nogen*). For the *Nogen* scenario, the relationship matrix was constructed based on a base population that was considered as a mixture of L and Y animals, which was not the case. Therefore, the results of the *Gen_AM* and *Gen_ADM* scenarios were more reliable than those of the *Nogen* scenario. Although predictive abilities were not significantly
different according to the Hotelling-Williams t-test, the results from 1000 bootstrap samples still showed that the predictive abilities of about 90% of the crossbred animals would be higher when genomic information was available (894 of 1000 bootstrap samples showed higher predictive abilities in scenarios that included genomic information than those in the *Nogen* scenario; results not shown). Comparison of the predictive abilities that were estimated in the current study with those from a previous study [17] indicated that the single-step model [16] might be more robust than the Vitezica model [22] used in this paper in terms of both predictive ability and unbiasedness for the crossbred performance. Our results suggested that using a small set of genotyped animals and pre-corrected data to implement genetic evaluation for crossbred performance was less powerful than using the whole dataset, which is similar to the conclusions for purebred performance [43].

The regression coefficients obtained with the Vitezica model were less than 1, which suggests that variations in total genetic effects could be overestimated (inflated). In terms of unbiasedness, there was no clear trend among the scenarios examined, regardless of whether genomic information was included or not. Overall, unbiasedness was not a problem in the current study because the regression coefficients in all scenarios did not significantly differ from 1.

Inbreeding depression for litter size in pigs is a well-known phenomenon [50, 51], and we found that inclusion of inbreeding effects in the model improved predictive abilities of crossbred animals. Estimates of inbreeding depression effects are rarely reported, but our estimates agree with those previously reported for commercial and Iberian pigs [52]. Inbreeding depression was, for the same amount of marker-based inbreeding, more detrimental in the Landrace than in the Yorkshire breed. There are many possible explanations among which the purging of lethal recessive alleles [53]. We also report an estimate of the inbreeding depression parameter for the crossbred animals, which is between the estimates for the parental breeds. To our knowledge, this estimate has never been reported.

The correlation between breeding values and dominance deviations is of theoretical concern [30]. However, this does not apply to the current marker-based analyses for the following reasons. (1) In a pedigree-based analysis, mating in an inbred population produces deviations from the Hardy-Weinberg equilibrium, which generate correlations between breeding values and dominance deviations [30]. However, in our study, SNPs are in Hardy-Weinberg equilibrium if allele frequencies are considered in the current generation. (2) Such a correlation occurs because the pedigree

information forces the genetic model to refer to the base population, since the state of alleles is not known, i.e. only probabilities of IBD are known. In our study, the states of alleles are known and the model can be described as referring to the current generation instead. (3) The equivalent GBLUP models in Eq. (3) used genotypic additive and dominance values, not breeding values and dominance deviations. A reasonable assumption in the model is that additive and dominance effects are unrelated at each SNP. Thus, covariance between additive and dominance genetic effects was ignored in the current study.

## Conclusions

We present for the first time the use of genomic inbreeding in crossbred and purebred genomic evaluation. Estimates are biologically sound and are relevant even for crossbred animals. We also report for the first time, estimated correlations of allele substitution effects in the presence of dominance. For TNB, the dominance genetic variance accounts for only a small proportion of the total genetic variation (4 to 10%). A moderate, positive genetic correlation between breeding values for TNB for purebred and crossbred performances was confirmed. Inclusion of dominance in the GBLUP model did not improve predictive ability for crossbred animals, whereas inclusion of inbreeding depression effects did. An additive GBLUP model is sufficient to capture the additive genetic variances and for genomic evaluation. The GBLUP model [22] was applied successfully for genetic evaluations for crossbred performance in pigs. This model can potentially be a useful tool in genetic evaluation for crossbred performance.

## Authors' contributions

TX performed data analysis and wrote the manuscript. All authors participated in the derivation of the theory. AL and OFC coordinated the project, conceived the study, made substantial contributions for the interpretation of results and revised the manuscript. ZGV improved the manuscript and added valuable comments during the study. All authors read and approved the final manuscript.

## Author details

## Acknowledgements

The work was funded through the Green Development and Demonstration Programme (Grant No. 34009-12-0540) by the Danish Ministry of Food, Agriculture and Fisheries, the Pig Research Centre and Aarhus University. The first author benefits from a joint grant from the European Commission and Aarhus University, within the framework of the Erasmus-Mundus joint doctorate "EGS-ABG". AL and ZGV thank financial support from the INRA SelGen metaprogram projects X-Gen and SelDir. OFC acknowledges funding from the GenSAP project. We are grateful to the Genotoul bioinformatics platform Toulouse Midi-Pyrenees for providing computing resources. Discussions with Luis Varona are gratefully acknowledged.

<footnote n="1">¹Department of Molecular Biology and Genetics, Center for Quantitative Genetics and Genomics, Aarhus University, 8830 Tjele, Denmark. ²UR1388 GenPhySE, INRA, CS-52627, 31326 Castanet-Tolosan, France. ³INP, ENSAT, GenPhySE, Université de Toulouse, 31326 Castanet-Tolosan, France.</footnote>
**Competing interests**

The authors declare that they have no competing interests.

Received: 11 April 2016 Accepted: 15 November 2016

Published online: 25 November 2016

**References**

1. Wei M, Van der Werf JHJ. Maximizing genetic response in cross-breds using both purebred and crossbred information. *Anim Sci.* 1994;59:401–13.

2. Bijma P, Bastiaansen JWM. Standard error of the genetic correlation: how much data do we need to estimate a purebred–crossbred genetic correlation? *Genet Sel Evol.* 2014;46:79.

3. Wei M, Van der Steen HAM. Comparison of reciprocal recurrent selection with pure-line selection systems in animal breeding (a review). *Anim Breed Abstr.* 1991;59:281–98.

4. Dekkers JCM. Marker-assisted selection for commercial crossbred performance. *J Anim Sci.* 2007;85:2104–14.

5. Lutaaya E, Misztal I, Mabry JW, Short T, Timm HH, Holzbauer R. Genetic parameter estimates from joint evaluation of purebreds and crossbreds in swine using the crossbred model. *J Anim Sci.* 2001;79:3002–7.

6. Lo LL, Fernando RL, Grossman M. Genetic evaluation by BLUP in two-breed terminal crossbreeding systems under dominance. *J Anim Sci.* 1997;75:2877–84.

7. Charlesworth D, Willis JH. The genetics of inbreeding depression. *Nat Rev Genet.* 2009;10:783–96.

8. Falconer DS, Mackay TFC. Introduction to quantitative genetics. New York: Longman Group Ltd; 1981.

9. Su G, Christensen OF, Ostersen T, Henryon M, Lund MS. Estimating additive and non-additive genetic variances and predicting genetic merits using genome-wide dense single nucleotide polymorphism markers. *PLoS One.* 2012;7:e45293.

10. Misztal I, Varona L, Culbertson M, Bertrand JK, Mabry J, Lawlor TJ, et al. Studies on the value of incorporating the effect of dominance in genetic evaluations of dairy cattle, beef cattle and swine. *Biotechnol Agron Soc Environ.* 1998;2:227–33.

11. Toro MA, Varona L. A note on mate allocation for dominance handling in genomic selection. *Genet Sel Evol.* 2010;42:33.

12. Ertl J, Legarra A, Vitezica ZG, Varona L, Edel C, Emmerling R, et al. Genomic analysis of dominance effects on milk production and conformation traits in Fleckvieh cattle. *Genet Sel Evol.* 2014;46:40.

13. Fulton JE. Genomic selection for poultry breeding. *Anim Front.* 2012;2:30–6.

14. Loberg A, Dürr JW. Interbull survey on the use of genomic information. *Interbull Bull.* 2009;39:3–14.

15. Christensen OF, Legarra A, Lund MS, Su G. Genetic evaluation for three-way crossbreeding. *Genet Sel Evol.* 2015;47:98.

16. Christensen OF, Madsen P, Nielsen B, Su G. Genomic evaluation of both purebred and crossbred performances. *Genet Sel Evol.* 2014;46:23.

17. Xiang T, Nielsen B, Su G, Legarra A, Christensen OF. Application of single-step genomic evaluation for crossbred performance in pig. *J Anim Sci.* 2016;94:936–48.

18. Zeng J, Toosi A, Fernando RL, Dekkers JCM, Garrick DJ. Genomic selection of purebred animals for crossbred performance in the presence of dominant gene action. *Genet Sel Evol.* 2013;45:11.

19. Esfandyari H, Bijma P, Henryon M, Christensen OF, Sorensen AC. Genomic prediction of crossbred performance based on purebred Landrace and Yorkshire data using a dominance model. *Genet Sel Evol.* 2016;48:40.

20. Hidalgo AM. Exploiting genomic information on purebred and crossbred pigs. PhD thesis. Swedish University of Agricultural Sciences, Uppsala; 2015.

21. Ibánez-Escriche N, Fernando RL, Toosi A, Dekkers JCM. Genomic selection of purebreds for crossbred performance. *Genet Sel Evol.* 2009;41:12.

22. Vitezica ZG, Varona L, Elsen MJ, Misztal I, Herring W, Legarra A. Genomic BLUP including additive and dominant variation in purebreds and F1 crossbreds, with an application in pigs. *Genet Sel Evol.* 2016;48:6.

23. Madsen P. DMU trace, a program to trace the pedigree for a subset of animals from a large pedigree file, version 2. Tjele: Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University; 2012.

24. Ramos AM, Crooijmans RP, Affara NA, Amaral AJ, Archibald AL, Beever JE, 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.* 2009;4:e6524.

25. GeneSeek Company. GGP-for Porcine LD (*GeneSeek Genomic Profiler for Porcine Low Density*). 2012. http://www.neogen.com/Genomics/pdf/Slicks/GGP_PorcineFlyer.pdf.

26. Xiang T, Ma P, Ostersen T, Legarra A, Christensen OF. Imputation of genotypes in Danish purebred and two-way crossbred pigs using low-density panels. *Genet Sel Evol.* 2015;47:54.

27. Browning SR. Missing data imputation and haplotype phase inference for genome-wide association studies. *Hum Genet.* 2008;124:439–50.

28. Silió L, Rodríguez M, Fernández A, Barragán C, Benítez R, Óvilo C, et al. Measuring inbreeding and inbreeding depression on pig growth from pedigree or SNP-derived metrics. *J Anim Breed Genet.* 2013;130:349–60.

29. Lynch M, Walsh B. Genetics and analysis of quantitative traits. 1st ed. Sunderland: Sinauer Assoc; 1998.

30. de Boer UM, Hoeschele I. Genetic evaluation methods for populations with dominance and inbreeding. *Theor Appl Genet.* 1993;86:245–58.

31. Aliloo H, Pryce JE, Gonzalez-Recio O, Cocks BG, Hayes BJ. Accounting for dominance to improve genomic evaluations of dairy cows for fertility and milk production traits. *Genet Sel Evol.* 2016;48:8.

32. Vitezica ZG, Varona L, Legarra A. On the additive and dominant variance and covariance of individuals within the genomic selection scope. *Genetics* 2013;195:1223–30.

33. Misztal I, Tsuruta S, Strabel T, Auvray B, Druet T, Lee DH. BLUPF90 and related programs (BGF90). In *Proceedings of the 7th world congress on genetics applied to livestock production*, Montpellier; 19–23 August, 2002.

34. Madsen P, Jensen J. A user's guide to DMU, version 6, release 5.2. Tjele: Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University; 2013.

35. VanRaden PM. Efficient methods to compute genomic predictions. *J Dairy Sci.* 2008;91:4414–23.

36. Houle D, Meyer K. Estimating sampling error of evolutionary statistics based on genetic covariance matrices using maximum likelihood. *J Evol Biol.* 2015;28:1542–9.

37. Lo LL, Fernando RL, Grossman M. Covariance between relatives in multi-breed populations: additive model. *Theor Appl Genet.* 1993;87:423–30.

38. Karoui S, Caraballo MJ, Díaz C, Legarra A. Joint genomic evaluation of French dairy cattle breeds using multiple-trait models. *Genet Sel Evol.* 2012;44:39.

39. Wientjes YC, Veerkamp RF, Bijma P, Bovenhuis H, Schrooten C, Calus MP. Empirical and deterministic accuracies of across-population genomic prediction. *Genet Sel Evol.* 2015;47:5.

40. Porto-Neto LR, Barendse W, Henshall JM, McWilliam SM, Lehnert SA, Reverter A. Genomic correlation: harnessing the benefit of combining two unrelated populations for genomic selection. *Genet Sel Evol.* 2015;47:84.

41. Meuwissen THE, Luo Z. Computing inbreeding coefficients in large populations. *Genet Sel Evol.* 1992;24:305–13.

42. Mäntysaari EA, Koivula M. GEBV validation test revisited. *Interbull Bull.* 2012;45:11–6.

43. Guo X, Christensen OF, Ostersen T, Wang Y, Lund MS, Su G. Improving genetic evaluation of litter size and piglet mortality for both genotyped and nongenotyped individuals using a single-step method. *J Anim Sci.* 2015;93:503–12.

44. Lopes MS, Bastiaansen JWM, Janss L, Knol EF, Bovenhuis H. Estimation of additive, dominance, and imprinting genetic variance using genomic data. *G3 (Bethesda)*. 2015;5:2629–37.

45. Meuwissen THE, Hayes BJ, Goddard ME. Prediction of total genetic value using genome-wide dense marker maps. *Genetics*. 2001;157:1819–29.

46. Hidalgo AM, Bastiaansen JWM, Lopes MS, Harlizius B, Groenen MAM, de Koning D-J. Accuracy of predicted genomic breeding values in purebred and crossbred pigs. *G3 (Bethesda)*. 2015;5:1575–83.
47. Legarra A, Baloche G, Barillet F, Astruc JM, Soulas C, Aguerre X, et al. Within-and across-breed genomic predictions and genomic relationships for Western Pyrenees dairy sheep breeds Latxa, Manech, and Basco-Béarnaise. J Dairy Sci. 2014;97:3200–12.

48. Moghaddar N, Swan AA, van der Werf JH. Comparing genomic prediction accuracy from purebred, crossbred and combined purebred and cross-bred reference populations in sheep. Genet Sel Evol. 2014;46:58.

49. Sun C, VanRaden PM, Cole JB, O'Connell JR. Improvement of prediction ability for genomic selection of dairy cattle by including dominance effects. PLoS One. 2014;9:e103934.

50. Dickerson GE. Inbreeding and heterosis in animals. In Proceedings of the animal breeding and genetics symposium in honor of Jay L. Lush, 29 July 1972, Blacksburg; 1973. p. 54–77.

51. Leroy G. Inbreeding depression in livestock species: review and meta-analysis. Anim Genet. 2014;45:618–28.

52. Silió L, Barragán C, Fernández AI, García-Casco J, Rodríguez MC. Assessing effective population size, coancestry and inbreeding effects on litter size using the pedigree and SNP data in closed lines of the Iberian pig breed. J Anim Breed Genet. 2016;133:145–54.

53. Hinrichs D, Meuwissen THE, Ødegard J, Holt M, Vangen O, Woolliams JA. Analysis of inbreeding depression in the first litter size of mice in a long-term selection experiment with respect to the age of the inbreeding. Heredity. 2007;99:81–8.

Submit your next manuscript to BioMed Central and we will help you at every step:

* We accept pre-submission inquiries

* Our selector tool helps you to find the most relevant journal

* We provide round the clock customer support

* Convenient online submission

* Thorough peer review

* Inclusion in PubMed and all major indexing services

* Maximum visibility for your research

Submit your manuscript at
www.biomedcentral.com/submit

( ) **BioMed** Central