- Research
- Open Access
Empirical and deterministic accuracies of across-population genomic prediction
- Yvonne CJ Wientjes^{1, 2}Email author,
- Roel F Veerkamp^{1, 2},
- Piter Bijma^{2},
- Henk Bovenhuis^{2},
- Chris Schrooten^{3} and
- Mario PL Calus^{1}
https://doi.org/10.1186/s12711-014-0086-0
© Wientjes et al.; licensee BioMed Central. 2015
- Received: 14 April 2014
- Accepted: 17 December 2014
- Published: 6 February 2015
Abstract
Background
Differences in linkage disequilibrium and in allele substitution effects of QTL (quantitative trait loci) may hinder genomic prediction across populations. Our objective was to develop a deterministic formula to estimate the accuracy of across-population genomic prediction, for which reference individuals and selection candidates are from different populations, and to investigate the impact of differences in allele substitution effects across populations and of the number of QTL underlying a trait on the accuracy.
Methods
A deterministic formula to estimate the accuracy of across-population genomic prediction was derived based on selection index theory. Moreover, accuracies were deterministically predicted using a formula based on population parameters and empirically calculated using simulated phenotypes and a GBLUP (genomic best linear unbiased prediction) model. Phenotypes of 1033 Holstein-Friesian, 105 Groninger White Headed and 147 Meuse-Rhine-Yssel cows were simulated by sampling 3000, 300, 30 or 3 QTL from the available high-density SNP (single nucleotide polymorphism) information of three chromosomes, assuming a correlation of 1.0, 0.8, 0.6, 0.4, or 0.2 between allele substitution effects across breeds. The simulated heritability was set to 0.95 to resemble the heritability of deregressed proofs of bulls.
Results
Accuracies estimated with the deterministic formula based on selection index theory were similar to empirical accuracies for all scenarios, while accuracies predicted with the formula based on population parameters overestimated empirical accuracies by ~25 to 30%. When the between-breed genetic correlation differed from 1, i.e. allele substitution effects differed across breeds, empirical and deterministic accuracies decreased in proportion to the genetic correlation. Using a multi-trait model, it was possible to accurately estimate the genetic correlation between the breeds based on phenotypes and high-density genotypes. The number of QTL underlying the simulated trait did not affect the accuracy.
Conclusions
The deterministic formula based on selection index theory estimated the accuracy of across-population genomic predictions well. The deterministic formula using population parameters overestimated the across-population genomic accuracy, but may still be useful because of its simplicity. Both formulas could accommodate for genetic correlations between populations lower than 1. The number of QTL underlying a trait did not affect the accuracy of across-population genomic prediction using a GBLUP method.
Keywords
- Quantitative Trait Locus
- Genetic Correlation
- Genomic Prediction
- Selection Candidate
- Estimate Breeding Value
Background
For genomic prediction, a reference population that consists of individuals with phenotypes and marker genotypes is used to estimate marker effects and to predict breeding values for another group of genotyped individuals, called selection candidates. The accuracy of predicting breeding values for selection candidates within one population is influenced by the level of linkage disequilibrium (LD) between markers, i.e. single-nucleotide polymorphisms (SNPs) and quantitative trait loci (QTL) that influence the trait, and by the level of family relationships [1-4]. Across populations, there are differences in LD, allele frequencies [4-6], and allele substitution effects of QTL [7,8], and close family relationships between individuals of different populations are absent. Therefore, the potential accuracy of predicting breeding values when the predicted population is not included in the reference population is likely to be limited. Indeed, in dairy cattle breeding, several empirical studies showed that the potential of using information across breeds was limited, e.g. [9-11]. The concept of combining individuals of different breeds in cattle is essentially similar to combining individuals from different lines in other animal and plant species, e.g. [4,12,13] or from different subpopulations in humans, e.g. [3,6] because close family relationships are absent and the extent of LD is limited across breeds, lines, and subpopulations.
A higher marker density may increase the consistency in LD phase across populations, since at short distances (5 to 30 kb) LD phases are conserved across populations [5]. However, several empirical studies showed that an increase in marker density resulted only in a small increase in accuracy using multiple populations in the reference population [14,15]. This small effect of marker density on accuracy indicates that other factors are also important, such as differences in segregating QTL or in the effect of QTL across populations due to differences in genetic background between populations [7,8]. DGAT1 (diacylglycerol O-acyltransferase 1) is one example of a gene with different effects across populations in dairy cattle. Allele substitution effects of a QTL in the DGAT1 locus on milk yield and fat yield have been found to be on average 0.8 and 0.5 times, respectively, as large in Jersey than in Holstein-Friesian populations in New Zealand [7] and 0.7 and 1.2 times, respectively, as large in Fleckvieh than in Holstein-Friesian populations in Germany [8]. Since the SNP that was analysed is considered to be the causal polymorphism, which rules out incomplete LD, these results demonstrate that large differences in allele substitution effects can exist across populations.
Another factor that may affect accuracy of genomic prediction across populations is the number of QTL underlying the trait. For genomic prediction based on one population, accuracy is shown to be independent of the number of QTL underlying the trait when a genomic best linear unbiased prediction method (GBLUP) is used [16,17], at least in situations for which there are no QTL that explain an extremely large part of the genetic variance. However, those studies only looked at the effect of the number of QTL on accuracy of genomic prediction within one population and not across populations.
For genomic prediction within one population, different deterministic formulas have been proposed to calculate the accuracy [1,2]. The formula of Daetwyler et al. [1] uses population and trait parameters, i.e. size of the reference population, heritability and number of effective chromosome segments. If the number of effective chromosome segments is calculated from the variation of genomic relationships around their expectations based on pedigree information, the formula of Daetwyler et al. [1] can also be applied for populations with a complex family structure [18]. The formula of VanRaden [2] can be derived both from selection index theory and prediction error variance of the mixed model equation and it estimates the accuracy using the relationships within the reference population and between selection candidates and the reference population. Hayes et al. [9] showed that applying the formula based on prediction error variance in multi-population situations without rescaling the genomic relationships across populations resulted in overestimation of the accuracy. This indicates that formulas for estimating the accuracy of genomic prediction using multiple populations need further investigation to define the best way to calculate genomic relationships across populations.
The first objective of this study was to develop a deterministic formula to estimate the accuracy of across-population genomic prediction. The second objective was to investigate the effect of differences in allele substitution effects of QTL across populations on accuracy of across-population genomic prediction. The last objective was to investigate the effect of the number of QTL underlying a trait on accuracy of across-population genomic prediction. Two deterministic formulas were evaluated and empirical accuracies were calculated using simulated phenotypes based on real genotypes from three cattle breeds representing different populations. Phenotypes were simulated using different correlations between allele substitution effects across breeds and different numbers of QTL underlying the trait. The reason for simulating the phenotypes of the individuals was to be able to investigate the actual effects of differences in allele substitution effects of QTL across populations and of the number of QTL by changing one factor at a time without changing the other factors, which would not be possible with real data.
Methods
Across-population genomic prediction
For genomic prediction based on one population, breeding values are predicted for individuals using a reference population of individuals from the same population. In most genomic prediction models, the QTL effects that underlie the traits of interest are assumed to be additive, e.g. [19]. For across-population genomic prediction, breeding values are predicted for individuals using a reference population of individuals from one or more different populations. Due to differences in allele frequencies across populations, the presence of non-additive effects can result in differences in allele substitution effects of QTL [20]. Therefore, the models used for across-population genomic prediction should include non-additive effects or allow for differences in allele substitution effects across populations. Since it is difficult to accurately estimate non-additive effects, e.g. [21,22], assuming additive gene action and, at the same time, allowing for differences in allele substitution effects may be a good first step and is the focus of this study. The correlation between allele substitution effects across populations can be considered as the genetic correlation between the populations [20,23].
Under the assumption that SNPs are representative of QTL, i.e. that characteristics such as allele frequency are the same for SNPs and QTL, resulting in usable LD between SNPs and QTL, a genomic relationship matrix based on SNPs can be used to represent the relationships between breeding values of the individuals. To calculate the genomic relationships, covariances between the individuals of both populations need to be calculated. The mathematical definition of a covariance, \( Cov\left(x,y\right)=E\left[\left(x-\overline{x}\right)\left(y-\overline{y}\right)\right] \), indicates that both components are corrected for their own mean. For the genomic relationships, this can be achieved by correcting the SNP genotypes of the individuals using the allele frequencies of their own population. Thus, the genotype of individual i from population j at locus k, g _{ ijk }, is standardized as \( {x}_{ijk}=\frac{g_{ijk}-2{p}_{jk}}{\sqrt{2{p}_{jk}\left(1-{p}_{jk}\right)}} \), where p _{ jk } is the allele frequency of population j at locus k, and the standardized genotypes are used to calculate the genomic relationship matrices using the method of Yang et al. [24], which will be described later.
This expression for the estimated breeding value (EBV) will subsequently be used in the next section to derive the accuracy.
Deterministic accuracy of across-population genomic prediction based on selection index theory
Deterministic accuracy of across-population genomic prediction using multiple populations in the reference population based on selection index theory
Although Equation 19 is derived for across-population genomic prediction, this formula can also be applied to estimate the accuracy of multi population genomic prediction for which one of the reference populations is the population of the selection candidates. Moreover, it is interesting to note that when one population is included in the reference population and selection candidates are from the same population as the reference individuals, Equation 19 becomes equivalent to the expression derived by VanRaden [2].
Deterministic accuracy of across-population genomic prediction based on population parameters
Simulations
Genotypes
Genotypes were available for 1285 dairy cows from the Netherlands that originated from three breeds (1033 Holstein-Friesian (HF), 105 Groninger White Headed (GWH), and 147 Meuse-Rhine-Yssel (MRY)). All individuals were pure-bred animals since at least 87.5% of their genes originated from one of the three breeds.
Individuals from the breeds GWH and MRY were genotyped with the Illumina BovineHD Beadchip (777 k, Illumina, San Diego, CA). Quality controls consisted in removing genotypes with a GenCall (GC) score lower than 0.2, SNPs with a call rate smaller than 95% in one of the breeds and SNPs with an unknown map position or located on the sex chromosomes. The HF individuals were genotyped with the Illumina BovineSNP50 Beadchip (50 k, Illumina, San Diego, CA), and imputed to high-density (777 k) using a reference population of 3150 HF individuals as described by Pryce et al. [26]. Quality control consisted in removing SNPs with a call rate smaller than 95% or with an unknown map position or located on the sex chromosomes. After editing the imputed genotypes, the mean Beagle R ^{2} value, which reflects the accuracy of imputation, was equal to 0.96 across imputed loci, which indicates that imputation was highly accurate.
Loci for which the genotypes passed the quality control of both the HF dataset and the combined GWH and MRY dataset were retained in the entire dataset. From this entire dataset, SNPs with a minor allele frequency equal to or lower than 0.5%, SNPs for which only two genotypes were observed, and SNPs in complete LD (r^{2} = 1) with an adjacent SNP were removed. To increase the power of accurately estimating genomic breeding values, arbitrarily, we took only three chromosomes, namely chromosomes 13, 23 and 28 that contained about 10% of the remaining high-density SNPs into account. According to the literature, the LD pattern of those chromosomes is comparable to the LD pattern of the entire cattle genome [27,28]. After editing, a total of 31 503 SNPs remained across the three chromosomes.
Simulation of phenotypes
Phenotypes of the individuals were simulated using different scenarios with two variables i.e. (1) the number of QTL underlying the simulated trait and (2) the correlation between allele substitution effects of the QTL underlying the simulated trait in the different populations, i.e. the genetic correlation between populations [20,23]. From the 31 503 SNPs available after editing, 5000 were randomly selected to become candidate QTL, regardless of the chromosome. In each replicate, the actual QTL with an effect on the trait were randomly sampled from those candidate QTL. The remaining (31 503 – 5000 =) 26 503 SNPs composed the group of markers used in all analyses. Using this approach allowed us to keep the set of markers constant across all replicates but still made it possible to randomly select the QTL from the group of candidate QTL within each replicate. The numbers of QTL underlying the simulated trait were equal to 3000 (~10% of all SNPs), 300 (~1%), 30 (~0.1%) or 3 (~0.01%).
The allele substitution effects of QTL were sampled from a multinormal distribution with mean 0 and standard deviation 1, assuming a correlation of 1, 0.8, 0.6, 0.4, or 0.2 between the allele substitution effects across all three pairs of breeds. This was simulated by sampling random numbers from a normal distribution with mean 0 and standard deviation 1 and multiplying those numbers with the Cholesky decomposition of the covariance matrix between the allele substitution effects of the breeds.
For each of the individuals, the TBV was calculated by multiplying the simulated allele substitution effects with the genotypes of the 3000, 300, 30, or 3 QTL coded as 0, 1, and 2. Only additive effects and no dominance effects or epistatic interactions were simulated, therefore, the effects were summed over all QTL. Finally, TBV of all individuals of the three breeds were rescaled to a mean of 0 and variance of 1 across breeds. By rescaling the TBV in this way, their mean and variance were the same for each replicate and for the different numbers of QTL, which indicates that when the number of QTL was higher, each QTL explained a smaller part of the variance.
Allele frequencies for simulated QTL (sampled from the SNPs) differed for each of the three breeds, resulting in differences in average TBV between the breeds. To simulate environmental effects for each individual assuming equal heritability for the three breeds, TBV were first adjusted by subtracting the average TBV of the individual’s breed before the genetic variance across TBV was calculated. Thereafter, the environmental effect per individual was sampled for the three breeds from a normal distribution with mean 0 and variance \( \left(\frac{1}{h^2}-1\right) \)*(variance of TBV corrected for mean TBV within breed). For each individual, the phenotype was calculated as the sum of its TBV and the randomly sampled environmental effect. Note that the within-breed TBV means were only subtracted from the TBV to calculate the environmental variance, the TBV itself, and therefore the phenotypes as well, still included the within-breed TBV mean.
Scenarios to evaluate accuracy of genomic prediction
Overview of the breeds used in the different reference populations and as selection candidates
Reference population | Predicted individuals | |||
---|---|---|---|---|
Scenario | Breed(s) | Nb of individuals | Breed | Nb of individuals |
Base | HF | 1032/981-982^{1} | HF | 1/51-52^{1} |
1 | HF | 1033 | GWH | 105 |
2 | HF + MRY | 1180 | GWH | 105 |
3 | HF | 1033 | MRY | 147 |
4 | HF + GWH | 1138 | MRY | 147 |
Empirical accuracy based on simulated phenotypes
For the empirical estimation of the accuracy, a GBLUP–model type, called GREML, was run in ASReml [30]. This GREML model used a genomic relationship matrix (G) and simulated phenotypes based on 3000, 300, 30 or 3 QTL underlying the simulated trait. In this model, breed was included as a fixed effect. This model is termed GREML, because it has the same features as the commonly known GBLUP model, however variances were not assumed to be known but were estimated simultaneously with the breeding values using REML. Accuracy was calculated for each population as the correlation between EBV from this model and TBV. Since simulated phenotypes were different per replicate, averages and standard errors of empirical accuracies were calculated across replicates.
The G matrix used in GREML contained all reference individuals and selection candidates and was calculated based on the method of Yang et al. [24]; \( {\mathbf{G}}_{SNPs}=\frac{\mathbf{XX}\hbox{'}}{n} \). In this equation, n represents the number of SNPs (26 503) and the X matrix contains standardized genotypes (one locus per column) of each individual (one individual per row). For the empirical estimation of the accuracy, standardized genotypes were calculated as \( {x}_{ij}=\frac{g_{ij}-2{p}_j}{\sqrt{2{p}_j\left(1-{p}_j\right)}} \), where g _{ ij } codes the genotype for individual i at marker locus j as 0, 1 and 2, and p _{ j } is the allele frequency at marker locus j for the second allele averaged over all breeds. To calculate the average allele frequency per locus, the allele frequency per locus was calculated per breed and thereafter averaged over the three breeds, with an equal weight for each of the breeds. In that way, average allele frequency is not dominated by the breed with the largest number of genotyped individuals. Note that for each scenario, the G _{ SNPs } matrix contained only the reference individuals and selection candidates (and the SNPs segregating in that group), so four different G _{ SNPs } matrices were calculated that contained (1) all HF individuals (26 486 SNPs), (2) all HF and GWH individuals (26 500 SNPs), (3) all HF and MRY individuals (26 498 SNPs), and (4) all HF, GWH and MRY individuals (26 503 SNPs).
Since the level of family relationships influences the sampling error on the elements in G, the regression coefficient b was calculated separately for bins of family relationships in A (0–0.10, > 0.10-0.25, > 0.25-0.50 and > 0.5) within each breed and for each combination of breeds. Across-breed relationships were indeed 0 in A, so in that case Var(Ĝ − A) approximately reduced to Var(Ĝ). Parent-offspring relationships and self-relationships were not or hardly affected by sampling error and therefore excluded from the regression. The regression coefficient b was always above 0.95, and, in most cases, even above 0.99. Therefore, the effect of regressing the G matrix back to the A matrix was limited.
The inbreeding level in A depends on the depth of the pedigree, which indicates that different pedigree depths across populations can cause differences in inbreeding levels across the populations. To remove these differences in pedigree depth, the pedigree was cut off at seven generations for all individuals. Based on the pedigree, small relationships between some animals of the different breeds occurred, with a maximum relationship of 0.035 between HF and GWH, 0.034 between HF and MRY, and 0.029 between GWH and MRY. These relationships resemble more or less the relationship between an individual and one of its ancestors five generations back.
Deterministic accuracies of genomic prediction
For each scenario, accuracies of genomic prediction were deterministically derived using the two methods explained before; one method based on selection index theory (Equation 19) and one method based on population parameters (Equation 24). It is interesting to note that the formula based on selection index theory provides a single accuracy for each selection candidate, while the formula using population parameters provides an accuracy that applies to all selection candidates of the same population. Both deterministic methods calculate the accuracy based on genomic relationships and do not use phenotypes. Since the subset of SNPs was constant across all replicates and scenarios with different numbers of QTL, only one accuracy was calculated that applied to all replicates and numbers of QTL. Therefore, it was not possible to calculate standard errors across replicates for the deterministic accuracies.
Estimating genetic correlations between populations
In this simulation study, the genetic correlation between populations was known. In studies using real data, this is usually not the case and the genetic correlation needs to be estimated from the data. We investigated how accurate the genetic correlations between HF and GWH, and between HF and MRY were when using a multi-trait model in ASReml [30] in which the same trait in different breeds was treated as different traits. Within the multi-trait model, the same G matrix was used as in the GBLUP model, the environmental correlation was set to 0 and genetic and environmental variances of GWH and MRY animals were fixed at the simulated values, because the small number of animals in those breeds made it difficult to estimate variance components reliably.
Results
Differences between populations
In this study, accuracy of genomic prediction was evaluated by using genotypes of three cattle breeds. In cases where allele substitution effects were equal across breeds, differences in accuracy between single- and across-breed genomic predictions were due to differences in allele frequencies, relationships and LD pattern across breeds. The correlation between allele frequencies of all 26 503 SNPs was 0.67 for HF and GWH, 0.73 for HF and MRY, and 0.65 for GWH and MRY. Correlations of allele frequencies of SNPs and candidate QTL across breeds were similar.
Based on pedigree information, there were few differences in average relationships between breeds with average relationships of 0.0004 between HF and GWH (ranging from 0 to 0.035, 0.0004 between HF and MRY (ranging from 0 to 0.034), and 0.0005 between GWH and MRY (ranging from 0 to 0.029). Based on genotype data, differences in average relationships across breeds became more pronounced, with average relationships of −0.084 between HF and GWH (ranging from −0.194 to +0.115), −0.050 between HF and MRY (ranging from −0.151 to +0.125), and −0.098 between GWH and MRY (ranging from −0.184 to +0.088).
Equal allele substitution effects across populations
Accuracies for the base scenario, for which breeding values of HF individuals were predicted using a reference population of HF individuals, were very high (>0.9). Empirically derived accuracies were the same for the different numbers of QTL underlying the trait, which indicates that the number of QTL did not affect empirical accuracy in single-breed genomic prediction. With both deterministic methods, accuracies were in good agreement with the empirically-derived accuracies.
Accuracies with the other four scenarios, for which across-breed genomic prediction was applied, were much lower than those with the base scenario, but still ranged from 0.4 to 0.5. In each scenario, empirical accuracies using different numbers of QTL underlying the trait were very similar, which indicates that there is no effect of number of QTL on empirical accuracy. As with single-breed genomic prediction, estimated accuracies based on selection index theory were in good agreement with empirical accuracies for all four scenarios of across-breed genomic prediction. The deterministic prediction formula using population parameters overestimated empirical accuracies by about 25%.
Empirical accuracies as well as deterministic accuracies were slightly higher for selection candidates from breed GWH than for those from breed MRY. For both breeds, empirical and deterministic accuracies slightly increased when the other breed was added to the HF reference population, thus maintaining a near constant difference in accuracy between GWH and MRY individuals.
Different allele substitution effects across populations
As expected, deterministic and empirical accuracies were about equal to the accuracies obtained with equal allele substitution effects across breeds multiplied by the correlation between allele substitution effects. Empirical accuracies across the different numbers of QTL underlying the trait were again very similar, although those obtained with the 3-QTL scenario seemed to differ slightly from the other scenarios. This is in agreement with the much higher standard error across the replicates obtained with the 3–QTL scenario than with the 3000-, 300- or 30-QTL scenarios.
As in scenarios with equal allele substitution effects across breeds, accuracies obtained with the formula based on selection index theory were in good agreement with empirical accuracies. This indicates that this formula can be used to estimate the accuracy even when the genetic correlation between populations differs from 1. The formula using population parameters overestimated empirical accuracies by about 25% to 30%, regardless of the genetic correlation between breeds.
Estimated genetic correlations between populations
Simulated and estimated genetic correlations (standard errors across replicates) between the populations
Simulated genetic correlation | Estimated genetic correlation (s.e.) | ||||||||
---|---|---|---|---|---|---|---|---|---|
Populations | 3000 QTL | 300 QTL | 30 QTL | 3 QTL | |||||
HF - GWH | 1.0 | 0.913 | (0.012) | 0.915 | (0.012) | 0.889 | (0.014) | 0.860 | (0.022) |
HF - GWH | 0.8 | 0.791 | (0.016) | 0.785 | (0.014) | 0.769 | (0.021) | 0.555 | (0.049) |
HF - GWH | 0.6 | 0.605 | (0.022) | 0.601 | (0.023) | 0.568 | (0.027) | 0.527 | (0.053) |
HF - GWH | 0.4 | 0.472 | (0.024) | 0.507 | (0.026) | 0.440 | (0.031) | 0.309 | (0.061) |
HF - GWH | 0.2 | 0.194 | (0.028) | 0.218 | (0.027) | 0.203 | (0.039) | 0.156 | (0.067) |
HF - MRY | 1.0 | 0.888 | (0.013) | 0.893 | (0.012) | 0.911 | (0.012) | 0.851 | (0.024) |
HF - MRY | 0.8 | 0.806 | (0.017) | 0.778 | (0.016) | 0.806 | (0.016) | 0.690 | (0.039) |
HF - MRY | 0.6 | 0.614 | (0.022) | 0.688 | (0.018) | 0.616 | (0.024) | 0.456 | (0.054) |
HF - MRY | 0.4 | 0.440 | (0.022) | 0.450 | (0.025) | 0.444 | (0.028) | 0.281 | (0.062) |
HF - MRY | 0.2 | 0.244 | (0.024) | 0.253 | (0.025) | 0.239 | (0.038) | 0.235 | (0.061) |
Discussion
Deterministic accuracy of across-population genomic prediction
The first objective of this study was to develop a deterministic formula to investigate the accuracy of across-population genomic prediction. Our study as other previous studies [2,18,32] shows that the formula based on selection index theory (Equation 19) and the formula using population parameters (Equation 24) can accurately estimate the accuracy of genomic prediction within one population using relationship matrices. By setting up across-population genomic relationship matrices based on population-specific allele frequencies, it was also possible to accurately estimate the accuracy of across-population genomic prediction based on selection index theory. The application of the prediction formula using population parameters, as described in our study, overestimated the empirical accuracy for across-population genomic prediction in all scenarios by about 25 to 30%.
The genetic correlation in the deterministic formulas accounts for differences in allele substitution effects across populations. These differences may also lead to differences in genetic variances across populations, i.e. heterogeneous variances. For example, among populations, the genetic variance tends to be larger for the population with the highest mean for a given trait [33,34]. In addition, differences in allele frequencies across populations may also lead to heterogeneous variances; for example, a QTL may only segregate in one of the populations, which results in differences in the genetic variance explained by that QTL across populations although the actual allele substitution effects could be the same. Moreover, environmental variances may be different across populations when deregressed proofs of bulls are used as phenotypes, since the heritability of those proofs depends on the number of daughters of the bull, which can differ across populations. Heterogeneous variances across populations, which are not properly accounted for, may affect bias and accuracy of EBV. The deterministic formula based on selection index theory can take those heterogeneous variances into account as well, in contrast to the application of the formula based on populations parameters described here. Makgahlela et al. [35] empirically showed that accuracies of multi-breed genomic prediction can be increased by accounting for those heterogeneous variances across breeds in a multi-trait random regression model [35,36].
The genomic relationship matrix used in the deterministic formulas was calculated based on population-specific allele frequencies. Harris and Johnson [37] already mentioned that differences in allele frequencies should be taken into account to calculate genomic covariances and relationships between individuals of different populations. Not using population-specific allele frequencies results in average genomic relationships across populations different from 0 [38], large differences in average diagonal elements across populations [12,37] and overestimation of the accuracies [9]. In our study, using population-specific allele frequencies resulted in average genomic relationship close to 0, i.e. equal to 0.00003 with a standard deviation of 0.023 between HF and GWH, and 0.00003 with a standard deviation of 0.020 between HF and MRY.
The deterministic formula based on selection index theory (Equation 19) estimated the accuracy of across-population genomic prediction accurately for all scenarios. With a genetic correlation of 0.8, 0.6, 0.4, or 0.2, empirical and deterministic accuracies were respectively 80%, 60%, 40%, or 20% of the accuracies achieved with a genetic correlation of 1. This indicates that the deterministic formula can be used to estimate genetic correlations between populations (but does not provide information about the mechanism underlying this correlation); for example when the empirical accuracy is only 60% of the accuracy estimated assuming a genetic correlation of 1, the actual genetic correlation between populations is expected to be 0.6. Using this deterministic formula to estimate the genetic correlation between populations can be especially attractive when only one of the populations has a small number of genotyped individuals.
Overestimation of accuracies with the formula using population parameters for the across-population scenarios is probably due to the inability of the SNPs to capture all the genetic variance in the selection candidates [39,40], which is an underlying assumption of this formula. The empirical accuracy was about 80% of the predicted accuracy, both when GWH individuals or MRY individuals were used as selection candidates. This indicates that only 80% of the genetic variance in the selection candidates was captured by the markers in the reference population, due to differences in LD and allele frequencies of QTL between the reference population and the selection candidates. This proportion of the genetic variance in the selection candidates captured by SNPs in the reference population is the maximum accuracy of genomic prediction for those populations based on the used SNP chip [39].
By using an estimation of the genetic variance in the validation population that can be captured by SNPs in the reference population, the formula based on population parameters becomes a useful formula to predict the accuracy of across-population genomic prediction. This formula is very simple to use and can assess expected accuracies before individuals are genotyped. However, an important question remains regarding which values to use for M _{ e } and the genetic correlation. In this study, M _{ e } were estimated based on the variation in genomic relationships between reference and selection individuals around their expectations based on pedigree information. Similarly to the single-population scenario, M _{ e } of the across-population scenarios were estimated based on the relationships across population. Using this approach, an M _{ e } of about 1800 was estimated when GWH individuals were used as selection candidates, and 2400 when MRY individuals were used as selection candidates, both when HF individuals were used as reference population. Since only 10% of the genome was taken into account, this M _{ e } should be multiplied by 10 to get the actual M _{ e } across those populations. In a previous study, an M _{ e } of 11 500 was obtained when reference individuals and selection candidates shared allele frequencies and LD patterns and of 122 000 when reference individuals and selection candidates shared only allele frequencies [18]. Across breeds, allele frequencies are different, but LD patterns may be partly the same, therefore, M _{ e } across breeds was indeed expected to fall within the values of those groups. This suggests that perhaps an M _{ e } of about 20 000 could be used to predict the accuracy of across-population genomic prediction for closely related cattle breeds and an M _{ e } of about 40 000 or more for more distantly related cattle breeds.
The actual genetic correlation between populations, which is needed in the prediction formula, is in practice not known and depends on the traits and populations of interest. However, we showed that this genetic correlation can be estimated quite accurately using a multi-trait model and high-density genotypes. Thus, it may be possible to estimate this genetic correlation in a limited number of animals and to use it to predict the accuracies of genomic selection for different scenarios.
Empirical accuracies of genomic prediction
The second objective of this study was to investigate the effect of differences in allele substitution effects of QTL between populations, i.e. genetic correlations that differ from 1, on accuracy of across-population genomic prediction. Our results showed that genetic correlations between populations that are smaller than 1 resulted in a reduced accuracy of across-population genomic prediction that is proportional to the genetic correlation.
In this study, it was assumed that SNPs are representative of QTL, i.e. that SNPs and QTL have the same characteristics. Regarding this assumption, we know that for most complex traits, QTL minor allele frequencies are expected to be low [24,41,42]. However, the SNPs on the chip were selected to have an intermediate allele frequency [43], resulting in ascertainment bias of these SNPs. These differences in allele frequencies indicate that, in practice, QTL and SNPs have other characteristics, thereby reducing LD between QTL and SNPs in empirical studies. In our study, QTL were selected from the SNPs on the chip, which did not completely cover the range of expected allele frequencies of the actual QTL. Therefore, LD between QTL and SNPs may be overestimated, which results in higher accuracies of genomic prediction. In a future study, we will investigate the effect of different QTL allele frequencies on the accuracy of multi-population genomic prediction using loci with different allele frequencies and representative of the whole genome.
Another assumption used in this study was that the trait of interest was only influenced by additive effects. Due to the existence of non-additive effects, the average effects of allele substitution depend on the QTL allele frequencies [20], and might therefore be different across populations. In this study, different effects were considered by simulating genetic correlations between populations that differed from 1. In general, empirical studies use additive models for across-population genomic prediction and provide much lower accuracies than those obtained in this study for a genetic correlation of 1, e.g. [9,11]. This suggests that either SNPs do not represent QTL or that non-additive effects are important for the traits of interest in empirical studies, or a combination of both, which is important biological information.
In this study, genetic correlations between populations of 1, 0.8, 0.6, 0.4, and 0.2 were used to simulate phenotypes. Our results showed that genetic correlations between populations can be estimated quite accurately from the data using a multi-trait model. To date, this was done only in a few empirical studies [38,44]. Karoui et al. [38] reported estimated genetic correlations between French dairy cattle breeds that ranged from 0 (fertility; Montbéliarde – Normande) to 0.79 (milk; Montbéliarde – Holstein), with only two out of nine estimated genetic correlations above 0.6. These empirical results show that genetic correlation between populations can differ from 1 and depends on the trait of interest.
Results of this study clearly show that genetic correlation between populations is an important parameter for across-population genomic prediction. The true genetic correlation between populations is not influenced by differences in LD between QTL and SNPs. It is worth noting that apart from differences in allele substitution effects, the genetic correlation can also differ from 1 because of different QTL for the same trait. In terms of accuracy, the value of the genetic correlation is important and not the underlying cause of this genetic correlation. In fact, the genetic correlation specifies the maximum accuracy that can be obtained with across-population genomic prediction, provided that the reference population is very large and the number of SNPs is large enough to find a consistent linkage phase across populations.
Effect of number of QTL
The third objective of this study was to investigate the effect of the number of QTL underlying a trait on accuracy of across-population genomic prediction, which was studied using a GBLUP method. The results showed that changing the number of QTL without changing any other parameter had no effect on the accuracy.
In the case of genomic prediction within one population, different studies have already shown that accuracies of genomic prediction using GBLUP do not depend on number of QTL underlying the trait [16,17]. If variable selection models were used for genomic prediction, higher numbers of QTL resulted in lower accuracies [16,17,45]. One of these studies also showed that variable selection models have an advantage over GBLUP when the number of QTL is below M _{ e } in genomic prediction within one population [16]. In across-population situations, M _{ e } is much larger than within one population [18], which suggests that, in those situations, it will be easier to have a number of QTL smaller than M _{ e } and, thus it is expected that the use of variable selection models will be beneficial.
Conclusions
The deterministic formula based on selection index theory, that was derived in this study, can accurately estimate the accuracy of across-population genomic prediction by using population-specific allele frequencies to set-up genomic relationship matrices. Another deterministic formula using population parameters overestimates the accuracy of across-population genomic prediction, because the SNPs in the reference population cannot capture all of the genetic variance in the selection candidates. However, this formula may still be useful because of its simplicity, and is expected to be much more accurate when the proportion of genetic variance in the selection candidates is known with reasonable accuracy and included in the formula. Moreover, the results of this study show that differences in allele substitution effects across populations reduce the accuracy of across-population genomic prediction, with a proportion equal to the correlation between allele substitution effects across populations. The number of QTL underlying a trait does not affect the accuracy of across-population genomic prediction when a GBLUP method is used.
Declarations
Acknowledgements
This study was financially supported by Breed4Food (KB-12-006.03-005-ASG-LR), a public-private partnership in the domain of animal breeding and genomics, and CRV (Arnhem, The Netherlands). The RobustMilk project and the National Institute of Food and Agriculture (NIFA) are acknowledged for providing the 50 k genotypes of the HF cows, and the gDMI consortium is acknowledged for imputing those to 777 k genotypes. The Dutch Milk Genomics Initiative and the project 'Melk op Maat', funded by Wageningen University (the Netherlands), the Dutch Dairy Association (NZO, Zoetermeer, the Netherlands), the cooperative cattle improvement organization CRV BV (Arnhem, the Netherlands), the Dutch Technology Foundation (STW, Utrecht, the Netherlands), the Dutch Ministry of Economic Affairs (The Hague, the Netherlands) and the Provinces of Gelderland and Overijssel (Arnhem, the Netherlands), are thanked for providing the 777 k genotypes of the GWH and MRY cows. The authors acknowledge Myrthe Maurice-van Eijndhoven for collecting the data of the GWH and MRY cows, and the herd owners for their help in collecting the data.
Authors’ Affiliations
References
- Daetwyler HD, Villanueva B, Woolliams JA. Accuracy of predicting the genetic risk of disease using a genome-wide approach. PLoS ONE. 2008;3:e3395.PubMed CentralPubMedView ArticleGoogle Scholar
- VanRaden PM. Efficient methods to compute genomic predictions. J Dairy Sci. 2008;91:4414–23.PubMedView ArticleGoogle Scholar
- De los Campos G, Vazquez AI, Fernando R, Klimentidis YC, Sorensen D. Prediction of complex human traits using the genomic best linear unbiased predictor. PLoS Genet. 2013;9:e1003608.PubMed CentralView ArticleGoogle Scholar
- Zhong S, Dekkers JCM, Fernando RL, Jannink JL. Factors affecting accuracy from genomic selection in populations derived from multiple inbred lines: A barley case study. Genetics. 2009;182:355–64.PubMed CentralPubMedView ArticleGoogle Scholar
- De Roos APW, Hayes BJ, Spelman RJ, Goddard ME. Linkage disequilibrium and persistence of phase in Holstein-Friesian, Jersey and Angus cattle. Genetics. 2008;179:1503–12.PubMed CentralPubMedView ArticleGoogle Scholar
- De los Campos G, Klimentidis YC, Vazquez AI, Allison DB. Prediction of expected years of life using whole-genome markers. PLoS ONE. 2012;7:e40964.PubMed CentralPubMedView ArticleGoogle Scholar
- Spelman RJ, Ford CA, McElhinney P, Gregory GC, Snell RG. Characterization of the DGAT1 gene in the New Zealand dairy population. J Dairy Sci. 2002;85:3514–7.PubMedView ArticleGoogle Scholar
- Thaller G, Krämer W, Winter A, Kaupe B, Erhardt G, Fries R. Effects of DGAT1 variants on milk production traits in German cattle breeds. J Anim Sci. 2003;81:1911–8.PubMedGoogle Scholar
- Hayes BJ, Bowman PJ, Chamberlain AJ, Verbyla K, Goddard ME. Accuracy of genomic breeding values in multi-breed dairy cattle populations. Genet Sel Evol. 2009;41:51.PubMed CentralPubMedView ArticleGoogle Scholar
- Schrooten C, Schopen GCB, Beatson P: Across-breed genomic evaluation based on BovineHD genotypes, and phenotypes of bulls and cows. In Proceedings of the 64th EAAP annual meeting 26–30 August 2013; Nantes. Wageningen: Wageningen Academic Publishers; 2013:229.Google Scholar
- Pryce JE, Gredler B, Bolormaa S, Bowman PJ, Egger-Danner C, Fuerst C, et al. Short communication: genomic selection using a multi-breed, across-country reference population. J Dairy Sci. 2011;94:2625–30.PubMedView ArticleGoogle Scholar
- Simeone R, Misztal I, Aguilar I, Vitezica ZG. Evaluation of a multi-line broiler chicken population using a single-step genomic evaluation procedure. J Anim Breed Genet. 2012;129:3–10.PubMedView ArticleGoogle Scholar
- Ibánẽz-Escriche N, Fernando RL, Toosi A, Dekkers JCM. Genomic selection of purebreds for crossbred performance. Genet Sel Evol. 2009;41:12.PubMed CentralPubMedView ArticleGoogle Scholar
- Harris BL, Creagh FE, Winkelman AM, Johnson DL. Experiences with the Illumina high density Bovine BeadChip. Interbull Bull. 2011;44:3–7.Google Scholar
- Erbe M, Hayes BJ, Matukumalli LK, Goswami S, Bowman PJ, Reich CM, et al. Improving accuracy of genomic predictions within and between dairy cattle breeds with imputed high-density single nucleotide polymorphism panels. J Dairy Sci. 2012;95:4114–29.PubMedView ArticleGoogle Scholar
- Daetwyler HD, Pong-Wong R, Villanueva B, Woolliams JA. The impact of genetic architecture on genome-wide evaluation methods. Genetics. 2010;185:1021–31.PubMed CentralPubMedView ArticleGoogle Scholar
- Clark SA, Hickey JM, Van Der Werf JHJ. Different models of genetic variation and their effect on genomic evaluation. Genet Sel Evol. 2011;43:18.PubMed CentralPubMedView ArticleGoogle Scholar
- Wientjes YCJ, Veerkamp RF, Calus MPL. The effect of linkage disequilibrium and family relationships on the reliability of genomic prediction. Genetics. 2013;193:621–31.PubMed CentralPubMedView ArticleGoogle Scholar
- Meuwissen THE, Hayes BJ, Goddard ME. Prediction of total genetic value using genome-wide dense marker maps. Genetics. 2001;157:1819–29.PubMed CentralPubMedGoogle Scholar
- Falconer DS, Mackay TFC. Introduction to quantitative genetics. Harlow: Pearson Education Limited; 1996.Google Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Wittenburg D, Melzer N, Reinsch N. Including non-additive genetic effects in Bayesian methods for the prediction of genetic values based on genome-wide markers. BMC Genet. 2011;12:74.PubMed CentralPubMedView ArticleGoogle Scholar
- Bohren BB, Hill WG, Robertson A. Some observations on asymmetrical correlated responses to selection. Genet Res. 1966;7:44–57.PubMedView ArticleGoogle Scholar
- Yang J, Benyamin B, McEvoy BP, Gordon S, Henders AK, Nyholt DR, et al. Common SNPs explain a large proportion of the heritability for human height. Nat Genet. 2010;42:565–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Goddard ME, Hayes BJ, Meuwissen THE. Using the genomic relationship matrix to predict the accuracy of genomic selection. J Anim Breed Genet. 2011;128:409–21.PubMedView ArticleGoogle Scholar
- Pryce JE, Johnston J, Hayes BJ, Sahana G, Weigel KA, McParland S, et al. Imputation of genotypes from low density (50,000 markers) to high density (700,000 markers) of cows from research herds in Europe, North America, and Australasia using 2 reference populations. J Dairy Sci. 2014;97:1799–811.PubMedView ArticleGoogle Scholar
- McKay SD, Schnabel RD, Murdoch BM, Matukumalli LK, Aerts J, Coppieters W, et al. Whole genome linkage disequilibrium maps in cattle. BMC Genet. 2007;8:74.PubMed CentralPubMedView ArticleGoogle Scholar
- Khatkar MS, Nicholas FW, Collins AR, Zenger KR, Cavanagh JAL, Barris W, et al. Extent of genome-wide linkage disequilibrium in Australian Holstein-Friesian cattle based on a high-density SNP panel. BMC Genom. 2008;9:187.View ArticleGoogle Scholar
- Mrode RA, Thompson R. Linear models for the prediction of animal breeding values. Wallingford: CABI Publishing; 2005.View ArticleGoogle Scholar
- Gilmour AR, Gogel B, Cullis B, Thompson R, Butler D, Cherry M, et al. ASReml User Guide Release 3.0. Hemel Hempstead: VSN International Ltd; 2009.Google Scholar
- Powell JE, Visscher PM, Goddard ME. Reconciling the analysis of IBD and IBS in complex trait studies. Nat Rev Genet. 2010;11:800–5.PubMedView ArticleGoogle Scholar
- Clark SA, Hickey JM, Daetwyler HD, van der Werf JHJ. The importance of information on relatives for the prediction of genomic breeding values and the implications for the makeup of reference data sets in livestock breeding schemes. Genet Sel Evol. 2012;44:4.PubMed CentralPubMedView ArticleGoogle Scholar
- Legates JE. Heritability of fat yields in herds with different production levels. J Dairy Sci. 1962;45:990–3.View ArticleGoogle Scholar
- Boldman KG, Freeman AE. Adjustment for heterogeneity of variances by herd production level in dairy cow and sire evaluation. J Dairy Sci. 1990;73:503–12.View ArticleGoogle Scholar
- Makgahlela ML, Mäntysaari EA, Strandén I, Koivula M, Nielsen US, Sillanpää MJ, et al. Across breed multi-trait random regression genomic predictions in the Nordic Red dairy cattle. J Anim Breed Genet. 2013;130:10–9.PubMedView ArticleGoogle Scholar
- Strandén I, Mäntysaari EA. Use of random regression model as an alternative for multibreed relationship matrix. J Anim Breed Genet. 2013;130:4–9.PubMedView ArticleGoogle Scholar
- Harris BL, Johnson DL. Genomic predictions for New Zealand dairy bulls and integration with national genetic evaluation. J Dairy Sci. 2010;93:1243–52.PubMedView ArticleGoogle Scholar
- Karoui S, Carabaño M, Díaz C, Legarra A. Joint genomic evaluation of French dairy cattle breeds using multiple-trait models. Genet Sel Evol. 2012;44:39.PubMed CentralPubMedView ArticleGoogle Scholar
- Daetwyler HD: Genome-wide evaluation of populations. PhD thesis. Wageningen: Wageningen University, Animal Breeding and Genomics Centre; 2009.Google Scholar
- Erbe M, Gredler B, Seefried FR, Bapst B, Simianer H. A function accounting for training set size and marker density to model the average accuracy of genomic prediction. PLoS ONE. 2013;8:e81046.PubMed CentralPubMedView ArticleGoogle Scholar
- Goddard ME, Hayes BJ. Mapping genes for complex traits in domestic animals and their use in breeding programmes. Nat Rev Genet. 2009;10:381–91.PubMedView ArticleGoogle Scholar
- Kemper KE, Goddard ME. Understanding and predicting complex traits: knowledge from cattle. Hum Mol Genet. 2012;21:R45–51.PubMedView ArticleGoogle Scholar
- Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, Heaton MP, et al. Development and characterization of a high density SNP genotyping assay for cattle. PLoS ONE. 2009;4:e5350.PubMed CentralPubMedView ArticleGoogle Scholar
- Legarra A, Baloche G, Barillet F, Astruc J, 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.PubMedView ArticleGoogle Scholar
- Coster A, Bastiaansen JWM, Calus MPL, Van Arendonk JAM, Bovenhuis H. Sensitivity of methods for estimating breeding values using genetic markers to the number of QTL and distribution of QTL variance. Genet Sel Evol. 2010;42:9.PubMed CentralPubMedView ArticleGoogle Scholar
Copyright
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.