# High-density marker imputation accuracy in sixteen French cattle breeds

Chris Hoze Hozé, Marie-Noëlle M.-N. Fouilloux, Eric Venot, François F.

Guillaume, Romain R. Dassonneville, Sébastien S. Fritz, Vincent Ducrocq, Florence Phocas, Didier Boichard, Pascal Croiseau

► To cite this version:

Chris Hoze Hozé, Marie-Noëlle M.-N. Fouilloux, Eric Venot, François F. Guillaume, Romain R. Dassonneville, et al.. High-density marker imputation accuracy in sixteen French cattle breeds. Genetics Selection Evolution, 2013, 45 (1), pp.1-11. 10.1186/1297-9686-45-33. hal-01001056

HAL Id: hal-01001056

https://hal.science/hal-01001056

Submitted on 28 May 2020

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
High-density marker imputation accuracy
in sixteen French cattle breeds

Chris Hozé¹,<fnref n="2" />,³*, Marie-Noëlle Fouilloux⁴, Eric Venot¹,<fnref n="2" />, François Guillaume¹,<fnref n="2" />,⁴, Romain Dassonneville¹,<fnref n="2" />,⁴, Sébastien Fritz³, Vincent Ducrocq¹,<fnref n="2" />, Florence Phocas¹,<fnref n="2" />, Didier Boichard¹,<fnref n="2" /> and Pascal Croiseau¹,<fnref n="2" />

**Abstract**

**Background:** Genotyping with the medium-density Bovine SNP50 BeadChip® (50K) is now standard in cattle. The high-density BovineHD BeadChip®, which contains 777 609 single nucleotide polymorphisms (SNPs), was developed in 2010. Increasing marker density increases the level of linkage disequilibrium between quantitative trait loci (QTL) and SNPs and the accuracy of QTL localization and genomic selection. However, re-genotyping all animals with the high-density chip is not economically feasible. An alternative strategy is to genotype part of the animals with the high-density chip and to impute high-density genotypes for animals already genotyped with the 50K chip. Thus, it is necessary to investigate the error rate when imputing from the 50K to the high-density chip.

**Methods:** Five thousand one hundred and fifty three animals from 16 breeds (89 to 788 per breed) were genotyped with the high-density chip. Imputation error rates from the 50K to the high-density chip were computed for each breed with a validation set that included the 20% youngest animals. Marker genotypes were masked for animals in the validation population in order to mimic 50K genotypes. Imputation was carried out using the Beagle 3.3.0 software.

**Results:** Mean allele imputation error rates ranged from 0.31% to 2.41% depending on the breed. In total, 1980 SNPs had high imputation error rates in several breeds, which is probably due to genome assembly errors, and we recommend to discard these in future studies. Differences in imputation accuracy between breeds were related to the high-density-genotyped sample size and to the genetic relationship between reference and validation populations, whereas differences in effective population size and level of linkage disequilibrium showed limited effects. Accordingly, imputation accuracy was higher in breeds with large populations and in dairy breeds than in beef breeds. More than 99% of the alleles were correctly imputed if more than 300 animals were genotyped at high-density. No improvement was observed when multi-breed imputation was performed.

**Conclusion:** In all breeds, imputation accuracy was higher than 97%, which indicates that imputation to the high-density chip was accurate. Imputation accuracy depends mainly on the size of the reference population and the relationship between reference and target populations.

**Background**

The development of a high throughput chip of 54 001 single nucleotide polymorphisms (SNPs) for cattle, the Bovine SNP50 BeadChip® (50K), has drastically reduced genotyping costs and strongly contributed to the current implementation of genomic selection (GS). This approach,

which was first proposed by Meuwissen et al. [1], uses a
reference population (usually consisting of progeny-tested
bulls) with both genotypes and phenotypes to estimate
marker effects and then uses these estimates to predict
breeding values for animals without phenotypes. The ac-
curacy of prediction depends mainly on the size of the
reference population, heritability of the phenotypes, and
level of linkage disequilibrium (LD) between markers and
quantitative trait loci (QTL) [2,3]. In some bovine breeds,
a very large number of animals have been genotyped with
the 50K chip and it is now possible to predict breeding
values of animals at birth with high accuracy. In breeds

<footnote>* Correspondence: chris.hoze@jouy.inra.fr

¹INRA, UMR 1313 Génétique Animale et Biologie Intégrative, 78350 Jouy-en-Josas, France</footnote>

<footnote n="2">²AgroParisTech, UMR1313 Génétique Animale et Biologie Intégrative, 75231 Paris 05, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
with a limited number of progeny-tested bulls, assembling a large enough reference population is a real challenge. Under the assumption that LD is conserved across breeds, reference populations from different breeds can be combined to increase the size of the reference population. Conservation of LD, however, requires more than 300 000 informative SNPs and therefore, is not fulfilled with the classically used 50K chip [4]. The BovineHD BeadChip® (HD) that was developed in 2010 and contains 777K SNPs is expected to be sufficiently dense to detect conserved LD across breeds and allow multi-breed GS. Re-genotyping all animals on this HD chip is, however, not economically feasible but prediction (imputation) of HD genotypes from 50K genotypes is possible.

Several imputation methods have been implemented and are widely used to impute genotypes from one chip to another. They are based either on population LD (Fastphase [5], Beagle [6], MaCH [7], IMPUTE2 [8]), or on a combination of LD and family information (Fimpute [9], Dagphase [10], AlphaImpute [11], FindHap [12]). Many studies have compared these methods for imputation from low-density panels to the 50K chip. Beagle 3.3.0 has been shown to be an accurate software package [13-15] and is commonly used in bovine datasets.

In this study, Beagle 3.3.0 was used to study the accuracy of imputation from the Illumina 50K to the HD chip in 16 French cattle breeds and to investigate the main factors that affect this accuracy.

## Methods

### Genotypes

The dataset comprised 5153 animals genotyped with the HD chip. The cryopreserved semen, or blood samples of the animals included in our study, which were used for genotyping, were procured from various commercial AI organizations and breeder organizations through their routine practice in the framework of breeding programs. Therefore, no ethical approval was required for sampling of biological material. The HD chip contains 777 964 markers with an average probe spacing of 3.45 kilobases (kb) [16]. Animals belonged to 16 breeds (seven dairy breeds and nine beef breeds). The number of genotypes was not equally distributed across breeds (Table 1) but depended on population size. Animals genotyped with the HD chip were chosen based on their marginal contribution to the population, as defined by Boichard et al. [17], and computed using the PEDIG software [18]. If possible, two to three progeny of each founder were also genotyped in order to increase phasing accuracy, and the most influential progeny were chosen for HD-genotyping. In some situations, and in particular in breeds with small populations, the number of HD-genotyped animals was constrained by the availability of DNA from ancestors. The family structure of each breed is detailed in Table 1. The mean number of HD-genotyped male progeny per bull was higher in dairy breeds (2.5) than in beef breeds (2). In the Holstein breed, a large proportion of genotypes

<table><caption>Table 1 Number of high-density genotyped animals and population structure per breed</caption><thead><tr><th></th><th>Nb of genotyped animals</th><th>Nb of genotyped families (sire + progeny)</th><th>Mean nb of genotyped progeny per sire</th><th>Nb of effective ancestors</th></tr></thead><tbody><tr><td colspan="5"><strong>Dairy breeds</strong></td></tr><tr><td>Abondance (ABO)</td><td>209</td><td>54</td><td>3.69</td><td>15</td></tr><tr><td>Brown Swiss (BSW)</td><td>99</td><td>52</td><td>1.90</td><td>28</td></tr><tr><td>Holstein (HOL)</td><td>788</td><td>204</td><td>2.30</td><td>21</td></tr><tr><td>Montbéliarde (MON)</td><td>530</td><td>139</td><td>3.77</td><td>18</td></tr><tr><td>Normande (HOR)</td><td>551</td><td>138</td><td>3.82</td><td>23</td></tr><tr><td>Simmental (SIM)</td><td>125</td><td>55</td><td>2.24</td><td>39</td></tr><tr><td>Tarentaise (TAR)</td><td>185</td><td>65</td><td>2.77</td><td>15</td></tr><tr><td colspan="5"><strong>Beef breeds</strong></td></tr><tr><td>Aubrac (AUB)</td><td>254</td><td>116</td><td>2.17</td><td>112</td></tr><tr><td>Bazadaise (BAZ)</td><td>89</td><td>60</td><td>1.45</td><td>46</td></tr><tr><td>Blonde d'Aquitaine (BLA)</td><td>327</td><td>187</td><td>1.74</td><td>78</td></tr><tr><td>Charolais (CHA)</td><td>672</td><td>310</td><td>2.14</td><td>249</td></tr><tr><td>Gasconne (GAS)</td><td>163</td><td>76</td><td>2.12</td><td>197</td></tr><tr><td>Limousine (LIM)</td><td>462</td><td>235</td><td>1.96</td><td>185</td></tr><tr><td>Parthenaise (PAR)</td><td>304</td><td>97</td><td>3.02</td><td>89</td></tr><tr><td>Rouge des Prés (RDP)</td><td>149</td><td>80</td><td>1.83</td><td>99</td></tr><tr><td>Salers (SAL)</td><td>246</td><td>186</td><td>1.31</td><td>99</td></tr></tbody></table>
originated from the Eurogenomics consortium [19,20]. An-
other dataset, containing 33 746 animals genotyped with
the Bovine SNP50 BeadChip® [21], was also available. The
genotypes originated either from the national genomic se-
lection program or from complementary research pro-
grams, except for most of the Holstein genotypes, which
were from the Eurogenomics consortium.

HD and 50k genotypes were used for parentage testing
and to check the clustering quality of the HD genotypes
by concordance analysis of the genotypes of the 1838 in-
dividuals genotyped on both chips.

Data editing

Quality control was performed within-breed on the HD and 50K genotypes, using the same criteria for both chips. Genotyped animals with a call rate lower than 0.95 were removed from the analysis. Only markers mapped on the UMD3.1 assembly covering the 29 bovine autosomes were used for analysis. SNPs showing departure from Hardy-Weinberg equilibrium (p-value < 0.001) or with more than 10% missing genotypes were removed. In addition, genotype consistency was checked using 1838 animals that were genotyped on both chips and 352 markers that were discordant for more than 1% of these animals were excluded. After editing, 708 771 and 44 580 SNPs were retained on the HD and 50K sets, respectively. The 37 634 SNPs present on both HD and 50K sets were used to mimic 50K genotypes. Parentage was tested following the French routine genomic selection procedure, using both 50K and HD datasets (S. Fritz, personal communication). This procedure uses 500 informative markers and a parentage error was concluded if more than 10 incompatibilities were detected. In case of inconsistency, the progeny was removed, except when at least two progeny from the same sire were found incompatible with their sire. In this case, the sire was removed. Genotypes were checked for Mendelian inconsistencies between compatible parents and offspring. The genotype of the male parent was deleted if more than 20% of its progeny showed contradiction, in other cases the genotype of the progeny was set to missing.

Assessing within-breed imputation accuracy

The accuracy of imputation from 50k to HD genotypes was computed within each breed. For this purpose, the HD-dataset was divided into two parts. The oldest animals were used as a training population to mimic a reference population genotyped with the HD chip. The 20% youngest animals formed the validation population. For animals in the validation population, markers that were only present on the HD chip were masked to mimic a target population genotyped with the 50K.

Beagle 3.3.0 was used to impute mimicked 50K geno-
types up to high-density. This software uses a population-
based method called “localized haplotype-cluster method”,

as described by Browning and Browning [6]. The method
builds and clusters haplotypes along the whole chromo-
some and then uses an underlying variable length Markov
chain based on haplotypes counts (and consequently on
local LD patterns) to determine transition probabilities
from one marker to the next. The scale and shift parame-
ters were set to 2 and 0.1, respectively, and no pedigree in-
formation was taken into account.

Imputation accuracy was estimated based on the com-
parison between imputed and known HD genotypes and
was defined as the allelic imputation error rate, computed
as the ratio between the number of falsely imputed alleles
and the total number of imputed alleles [22].

Identification of SNPs with high error rates is important because they may induce errors in QTL detection and gen- omic selection. To identify these SNPs, only the breeds with the largest populations (Blonde d'Aquitaine, Charolais, Holstein, Limousine, Montbéliarde and Normande) were considered in order to avoid erroneous identification of er- rors due to the use of a low number of genotypes when computing the error rate for a SNP.

Factors affecting imputation accuracy

Given the large number of breeds included here, it was
possible to analyze the factors that may have an effect
on imputation accuracy. Two types of variables that dif-
fer between breeds were studied: genetic diversity indica-
tors such as the level of linkage disequilibrium and the
number of effective ancestors, and variables specific to
our datasets, such as the number of genotypes in the ref-
erence population and the relationship between refer-
ence and validation populations.

To avoid bias due to sampling, we used the indicators computed by Danchin-Burge [23] using the complete pedigree file of each breed. The number of effective ancestors was preferred to effective population size since the former is less sensitive to the depth of the pedigree [23].

Linkage disequilibrium was computed within-breed as a
squared correlation coefficient based on phased genotypes
for each marker pair [24], as defined in equation (1):

$$
r^2 = \frac{(P_{A1B1} - P_{A1}P_{B1})^2}{P_{A1}P_{A2}P_{B1}P_{B2}} \quad (1)
$$

where A1, A2, B1, B2 are the alleles of SNP A and B, pA1,
pA2, pB1, and pB2 are the corresponding allelic frequencies,
and pA1B1 is the frequency of the A1B1 haplotype.

For the Abondance breed, LD was computed for SNPs
on *Bos taurus* (BTA) chromosomes 5, 10, 15, 20 and 25
and no differences between chromosomes were observed
(data not shown). Therefore, for other breeds, LD com-
putation was only performed on BTA5. Pairs of SNPs
for which one SNP has a minor allele frequency (MAF)
lower than 5% in the considered breed were discarded
<figure><img src="image_2.png" /><figcaption>Figure 1 Linkage disequilibrium decay in 16 beef (dashed lines) and dairy (solid lines) cattle breeds.</figcaption></figure>

because it has been shown that LD tends to be very low in such cases [25]. Furthermore, removing SNPs with low MAF facilitates comparisons between populations [25,26]. Then, LD for all marker pairs were averaged by the distance between SNPs. We focused particularly on the level of LD between SNPs 70 kb apart, which corresponds to the average distance between two informative SNPs on the 50K chip.

Average relationship levels between training and validation populations were measured based on pedigree information using the PEDIG [18] software. A relationship coefficient was computed for each pair of individuals in each population and then averaged for each breed.

The effect of the different factors on imputation accuracy was assessed by multiple regression across all evaluated datasets (including both dairy and beef breeds).

<table><caption>Table 2 Within-breed imputation error rate and others parameters affecting imputation error rate</caption><thead><tr><th></th><th>Training population size</th><th>Validation population size</th><th>Allelic imputation error rate (%)</th><th>LD level at 70 kb</th><th>Average R<sub>T/V</sub></th></tr></thead><tbody><tr><td><strong>Dairy breeds</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Abondance (ABO)</td><td>169</td><td>40</td><td>0.75</td><td>0.217</td><td>0.146</td></tr><tr><td>Brown Swiss (BSW)</td><td>79</td><td>20</td><td>1.92</td><td>0.255</td><td>0.074</td></tr><tr><td>Holstein (HOL)</td><td>634</td><td>154</td><td>0.73</td><td>0.255</td><td>0.078</td></tr><tr><td>Montbéliarde (MON)</td><td>424</td><td>106</td><td>0.51</td><td>0.196</td><td>0.116</td></tr><tr><td>Normande (HOR)</td><td>444</td><td>107</td><td>0.33</td><td>0.233</td><td>0.104</td></tr><tr><td>Simmental (SIM)</td><td>100</td><td>25</td><td>2.55</td><td>0.209</td><td>0.050</td></tr><tr><td><strong>Beef breeds</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Aubrac (AUB)</td><td>204</td><td>50</td><td>2.03</td><td>0.177</td><td>0.028</td></tr><tr><td>Bazadaise (BAZ)</td><td>72</td><td>17</td><td>2.07</td><td>0.239</td><td>0.038</td></tr><tr><td>Blonde d'Aquitaine (BLA)</td><td>262</td><td>65</td><td>1.80</td><td>0.175</td><td>0.038</td></tr><tr><td>Charolais (CHA)</td><td>539</td><td>133</td><td>0.68</td><td>0.176</td><td>0.018</td></tr><tr><td>Gasconne (GAS)</td><td>131</td><td>32</td><td>2.26</td><td>0.174</td><td>0.026</td></tr><tr><td>Limousine (LIM)</td><td>370</td><td>92</td><td>1.09</td><td>0.164</td><td>0.014</td></tr><tr><td>Parthenaise (PAR)</td><td>245</td><td>59</td><td>1.88</td><td>0.161</td><td>0.024</td></tr><tr><td>Rouge des Prés (RDP)</td><td>119</td><td>30</td><td>2.39</td><td>0.206</td><td>0.028</td></tr><tr><td>Salers (SAL)</td><td>197</td><td>49</td><td>1.27</td><td>0.213</td><td>0.024</td></tr></tbody></table>

Size of the training and validation populations, within-breed imputation error rate, level of linkage disequilibrium (LD, r²)) at 70 kb and average relationship between training and validation populations (R<sub>T/V</sub>).
Tested factors were the reference population size, the number of effective ancestors, the linkage disequilibrium level at 70 kb and the relationship between training and validation populations. Factors were considered significant if the regression coefficient was different from zero (p-value < 0.05).

## Computation of multi-breed imputation accuracy

Since our aim was to improve imputation accuracy in breeds with a small reference population size, multi-breed imputation was tested. Training and validation populations from different breeds were combined and then imputation accuracy was computed as described above. Considering the diversity of the breeds involved in this project, it did not seem relevant to combine all breeds. Instead, breeds were grouped based on Reynolds genetic distances between breeds, which were computed by Gautier et al. [27] using allele frequencies on the 50K chip and represented in a Neighbor-Joining tree. One branch, which was composed of the Abondance, Montbéliarde and Tarentaise breeds, was taken as an example of closely related breeds and these breeds were merged for the multi-breed study. Imputation accuracy

was then computed per animal and averaged for each breed.

## Results

### Linkage disequilibrium decay

The results of LD decay are presented in Figure 1. Within the first 10 kb, LD strongly decreased in all breeds, thereafter it decreased at a lower rate and differences between breeds became noticeable. LD levels at 70 kb (Table 2) varied from 0.16 in the Parthenaise breed to 0.26 in the Holstein breed; on average LD levels at 70 kb were higher by 0.04 in dairy breeds than in beef breeds, in agreement with the difference in number of effective ancestors between these breeds.

### Imputation accuracy

Error rates were low in most breeds, with an overall mean error rate of 1.36% (Table 2). However, large differences were observed between breeds, with a minimum error rate of 0.31% in the Normande breed and a maximum of 2.41% in the Simmental breed. The error rates were less than 0.7% for breeds with more than 500 animals in the reference population.

<figure><img src="image_3.png" /></figure>
Error rates were lower in dairy breeds (mean error rate = 1.02%) than in beef breeds (mean error rate = 1.62%). In beef breeds, the error rate ranged from 0.64% in the Charolais breed to 2.26% in the Rouge des Prés breed. In dairy breeds, error rates were higher in the Simmental and Brown Swiss breeds, which are regional breeds in France with many imported ancestors and limited numbers of genotyped animals.

## Factors affecting imputation accuracy

Sizes of training and validation populations are presented in Table 2. Error rates decreased linearly when the number of animals in the reference population increased (Figure 2). Again, we observed a higher accuracy in dairy breeds: an error rate lower than 1% was achieved with only 100 to 200 genotyped animals in dairy breeds, whereas approximately 400 animals were required for beef breeds. With more than 400 genotyped animals, the impact of reference population size on accuracy leveled off, which suggests that the effect of size of the reference population is non-linear when a critical size is reached.

The relationship between the number of effective ancestors and imputation error rate is illustrated in Figure 3.

The number of effective ancestors appeared to have an opposite effect on dairy versus beef breeds. In dairy breeds, the number of effective ancestors was low (< 50) and its relationship with imputation error rate was positive. However, the two breeds with the highest number of effective ancestors (Simmental and Brown Swiss) among the dairy breeds had the highest imputation error rate but also the smallest reference populations and many foreign ancestors were not genotyped. In beef breeds, a negative relationship between the number of effective ancestors and imputation error rate was observed, but the breeds with the highest number of effective ancestors (Charolais and Limousine) also had the largest reference populations. This suggests that the effect of the number of effective ancestors was masked by the effect of the reference population size. However, when considering populations with equal reference population sizes (e.g. Limousine and Montbeliarde), the error rates increased with number of effective ancestors.

Figure 4 presents the relationship of the imputation error rate with the degree of genetic relationship between training and validation populations ($R_{T/V}$). The first observation is that $R_{T/V}$ was clearly lower in beef than in dairy breeds,

<figure><img src="image_4.png" /></figure>
<figure><img src="image_5.png" /></figure>

which may be explained by a higher number of effective ancestors in beef breeds. The second observation was that when breeds with the largest populations were discarded, the relationship between R<sub>T/V</sub> and imputation error rate appeared to be linear. Among breeds with large populations, the relationship between training and validation populations had only a limited effect on the imputation error rate. The Charolais and Limousine breeds had low imputation error rates despite their lower R<sub>T/V</sub> but these breeds have large reference populations. However, the lower R<sub>T/V</sub> in the Holstein breed probably explains why the error rate was higher in this breed than in the Normande or Montbeliarde breeds.

Finally, multiple linear regression was performed to better quantify the impact of each factor on imputation accuracy. Results of the multiple linear regressions are in Table 3. Reference population size and R<sub>T/V</sub> had a significant effect (p < 0.05) on imputation error rate, in contrast to the number of effective ancestors and the level of LD. About half of the variation in imputation error rates was explained by reference population size and one quarter by R<sub>T/V</sub>. Using the regression coefficients, we can predict that increasing R<sub>T/V</sub> by 1% reduces imputation error rate by 0.12% and that adding 100 animals to the reference population decreases the error rate by 0.26% (Table 3). This suggests that the size of the reference population is the major

<table><caption>Table 3 Results of the multiple linear regression model</caption><thead><tr><td></td><td>Part of variance explained</td><td>Regression coefficient</td><td>p-value</td></tr></thead><tbody><tr><td>Reference population size</td><td>52%</td><td>-0.0026 ± 0.0006</td><td>0.002</td></tr><tr><td>Relationship between reference and validation populations</td><td>25%</td><td>-12.6042 ± 3.8937</td><td>0.008</td></tr><tr><td>Level of linkage disequilibrium at 70 kb</td><td>2.5%</td><td>-0.0028 ± 0.0026</td><td>0.305</td></tr><tr><td>Number of effective ancestors</td><td>0.03%</td><td>-0.0576 ± 4.4641</td><td>0.899</td></tr></tbody></table>

The model tests the effect of reference population size, relationship between reference and validation populations, level of linkage disequilibrium at 70 kb, and number of effective ancestors on imputation error rate.
factor affecting imputation error rate but the relationship
between training and validation populations must also be
taken into consideration. Despite the difference between
breeds, our results suggest that the number of effective an-
cestors and the level of LD are not major factors affecting
imputation accuracy.

SNP by SNP analysis of imputation accuracy

Imputation error rate was computed for each SNP in order
to detect SNPs that were mismapped in the UMD3.1 as-
sembly. This analysis was performed for the six major
breeds (three beef breeds and three dairy breeds). Results
in the Montbéliarde breed are presented in Figure 5. Des-
pite an overall mean error rate of 0.51%, 13 104 and 6030
SNPs had error rates greater than 5 and 20%, respectively.
Consequently, the error rate dropped from 0.51 to 0.40%
after removing potentially mismapped SNPs with error
rates greater than 5%.

The relationship between MAF and imputation error
rate is illustrated in Figure 6 for the Montbéliarde breed.
SNPs were divided in two groups based on their im-
putation error rates. No relationship with MAF was
detected for SNPs with error rates less than 0.1, whereas
error rates increased with MAF for SNPs with high error
rates. Assuming that the LD between a mismapped SNP
and its direct neighbors is low, a mismapped SNP is im-
puted by chance and therefore its error rate is related to
its MAF.

We looked for SNPs with high error rates in each
breed based on a threshold that was defined as the mean

error rate plus 3 times its standard deviation for that
breed. Three thousand and eighty-three, 1980, and 1146
SNPs had high error rates in at least two, three or six
breeds, respectively. SNPs with high error rates in at
least three breeds are likely mismapped and are listed in
Additional file 1: Table S1.

Multi-breed imputation accuracy

Since imputation accuracy depends highly on the number
of HD genotypes, we increased the size of the reference
population by combining breeds. The results for single-
and multi-breed situations are in Table 4. No difference
in accuracy was found between single- and multi-breed
imputation for the Abondance and Montbéliarde breeds.
For the Tarentaise breed, imputation error rate was
slightly higher in the multi-breed analysis than in the
single-breed analysis.

Discussion

In this study, we evaluated accuracy of imputation from
50K to HD-genotypes for 16 cattle breeds and we inves-
tigated the corresponding causes of variation. We ob-
served large differences in imputation accuracy between
breeds. Several factors may explain these differences, i.e.
size of the reference population and the relationship be-
tween training and validation populations (closely related
to population structure). Results from the multiple linear
regression performed in this study, combined with other
published results, lead us to propose several hypotheses on
the impact of each factor.

<figure><img src="image_6.png" /><figcaption>Figure 5 Allelic imputation error rate along the genome in Montbéliarde breed.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>Figure 6 Relationship between allelic imputation error rate and minor allele frequency in Montbéliarde breed.</figcaption></figure>

The number of HD genotypes in the reference population ranged from 72 in the Bazadaise breed to 634 in the Holstein breed and is the major factor that explains differences in imputation accuracy between breeds. The imputation error rate decreased by 0.26% when 100 animals were added to the reference population. However, other studies have shown a non-linear effect of reference population size on imputation error rate. In Holstein cattle, Schrooten et al. (unpublished data) reported that imputation error rate decreased by between 0.17 and 0.04% when moving from 200 to 500 HD genotypes by steps of 100, while in sheep, the decrease was 5% when moving from 50 to 150 individuals in the reference population but only 2% when moving from 204 to 2512 individuals [28]. In our study, the effect of the size of the reference population on imputation error rate was linear. This effect would probably have been found non-linear

<table><caption>**Table 4 Imputation error rate using single-breed populations compared to a multi-breed reference population for three breeds**</caption><thead><tr><th></th><th>Abondance</th><th>Montbéliarde</th><th>Tarentaise</th></tr></thead><tbody><tr><td>Size of training / validation population</td><td>159 / 40</td><td>422 / 106</td><td>146 / 36</td></tr><tr><td>Single-breed imputation error rate (%)</td><td>0.755</td><td>0.487</td><td>0.763</td></tr><tr><td>Multi-breed imputation error rate (%)</td><td>0.753</td><td>0.485</td><td>0.824</td></tr></tbody></table>

if more breeds with a large reference population had been included. In fact, the size of the reference population had a limited effect when the number of HD genotypes was greater than a minimum threshold that was estimated at 200–400 animals.

Hayes et al. [28] reported that most of the differences in imputation accuracy are due to differences in the relationship between the reference and target populations and the genetic diversity of the breed. Schrooten et al. (unpublished data) used traceability, defined as “the expected contribution of HD-genotyped ancestors to the genotype of an animal”, as a measure of the relationship between the reference population and one animal from the target population. Imputation error rates were lower for animals with higher traceability, meaning that a higher average traceability, i.e. a higher relationship between the reference and target populations, will result in lower error rates. We reached the same conclusion, i.e. the imputation error rate was decreased by 0.12% when the average relationship increased by 0.01.

Increasing the size of the reference population decreases the probability to miss a haplotype in the reference population. For a fixed reference population size, an increase in the number of effective ancestors, i.e. an increase in the number of haplotypes in the total population, increases the probability to miss a haplotype and thus increases the error rate. This explains why reference population size and number of effective ancestors had opposite effects on imputation accuracy and compensate for
each other. In our study, differences in within-breed diversity explain why the error rate was higher in beef breeds than in dairy breeds and why more than 99% accuracy was achieved with only 200 animals for dairy breeds, while 400 animals were necessary for beef breeds. The poor results obtained with the Simmental and Brown Swiss breeds were also due to lower relationships between the reference and target populations; in these breeds, key ancestors mainly originate from abroad and were not included in our data. Combining reference populations from different breeds did not improve imputation accuracy, which confirms the results on multi-breed imputation of Hayes et al. [28] and Erbe et al. [29]. In fact, multi-breed imputation is expected to improve imputation accuracy only when 50K haplotypes are conserved across breeds, which is quite unlikely given the history of the breeds and their estimated divergence time, even for the most closely related breeds.

Although observed imputation error rates were low in all breeds, 1980 SNPs had particularly high error rates (Figure 5 and Additional file 1: Table S1), which suggests that errors exist in the marker map. Erbe et al. [29] identified 1231 SNPs with genotype error rates greater than 20%. When some of these SNPs were remapped using LD, error rates dropped. In our study, no remapping was performed but removing SNPs with high error rates resulted in a 0.1% drop in error rate, which suggests that lower error rates can be achieved with a more accurate map. However, Erbe et al. [29] still found 630 poorly imputed SNPs after remapping, which means that other reasons, such as recombination hot spots or regions on the 50K panel with lower SNP density, explain the high imputation error rates for some SNPs. SNPs with high imputation error rates likely also impact the quality of genomic selection and QTL detection. This has not been specifically investigated, but some studies have focused on the impact of imputation from low-density panels on reliability of genomic selection [14,30] and concluded that an imputation error rate between 2 and 3% leads to a mean loss of reliability of 2%. For imputation from 50K to the HD chip, the mean error rate is close to 1%, which suggests that the impact on reliability of genomic selection is even lower. However, because of the large number of markers available after imputation, it is preferable to discard markers with high error rates.

## Conclusions

The mean error rate for imputation from the Illumina Bovine50K® to the BovineHD® was around 1%. Differences in error rates between breeds were large and ranged from 2.41% in the Simmental breed to 0.31% in the Normande breed. These differences were mainly due to the size of the reference population and the relationship between the reference and target populations. This means that imputation accuracy could be increased by

increasing the number of HD genotypes and by improving the reference population to maximize its relationship with the population to impute. Using 50K genotypes to impute HD genotypes is possible, which implies that a large HD imputed reference population can be available for genomic selection at low cost. However, new HD genotypes are likely required in the future generations in order to maintain relationship links between the reference and target populations and limit the imputation error rate.

## Additional file

**Additional file 1: High_Error_Rate_SNP.** Csv file with header and separator ' ' containing the list of markers with high imputation error rate. Columns are as follows: chromosome number, SNP name, number of breeds for which a high error rate was detected and average error rates across breeds.

## Competing interests

The authors declare that they have no competing interests.

## Authors' contributions

CH and MNF performed the analysis. DB, VD, FP and PC designed the study. EV and SF managed the databases. FG and RD provided computer programmes. CH, DB, VD, PC wrote the manuscript. All authors read and approved the final manuscript.

## Acknowledgements

GEMBAL project is funded by the Agence Nationale de la Recherche (ANR-10-GENM-0014), APISGENE, Races de France and INRA "AIP Bioressources". The EUROGENOMICS consortium provided most of the Holstein genotypes. Most 50K genotypes originated from the Cartofine-ANR-05-GENANIMAL-007 project funded by ANR (French National Research Agency) and ApisGene, from the Qualvigene project funded by Apisgene, from the EUROGENOMICS consortium, and from genomic selection activity generated by the French cattle breeding companies.

## Author details

¹INRA, UMR 1313 Génétique Animale et Biologie Intégrative, 78350 Jouy-en-Josas, France. ²AgroParisTech, UMR1313 Génétique Animale et Biologie Intégrative, 75231 Paris 05, France. ³Union Nationale des Coopératives agricoles d'Elevage et d'Insémination Animale, 149 rue de Bercy, 75595 Paris Cedex 12, France. ⁴Institut de l'Elevage, 149 rue de Bercy, 75595 Paris Cedex 12, France.

Received: 8 February 2013 Accepted: 19 July 2013

Published: 3 September 2013

## References

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

2. Hayes BJ, Goddard ME: Technical note: prediction of breeding values using marker-derived relationship matrices. *J Anim Sci* 2008, 86:2089-2092.

3. Goddard M: Genomic selection: prediction of accuracy and maximisation of long term response. *Genetica* 2009, 136:245-257.

4. 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-1512.

5. Scheet P, Stephens M: A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. *Am J Hum Genet* 2006, 78:629-644.
6. Browning SR, Browning BL: Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. *Am J Hum Genet* 2007, **81**:1084-1097.

7. Li Y, Willer CJ, Ding J, Scheet P, Abecasis GR: MaCH: using sequence and genotype data to estimate haplotypes and unobserved genotypes. *Genet Epidemiol* 2010, **34**:816-834.

8. Howie BN, Donnelly P, Marchini J: A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. *PLoS Genet* 2009, 5:e1000529.

9. Sargolzaei M, Chesnais JP, Schenkel F: Flmpute - An efficient imputation algorithm for dairy cattle populations. *J Dairy Sci* 2011, **94**:421.

10. Druet T, Georges M: A hidden Markov model combining linkage and linkage disequilibrium information for haplotype reconstruction and quantitative trait locus fine mapping. *Genetics* 2010, **184**:789-798.

11. Hickey JM, Kinghorn BP, Tier B, Wilson JF, Dunstan N, van der Werf JH: A combined long-range phasing and long haplotype imputation method to impute phase for SNP genotypes. *Genet Sel Evol* 2011, **43**:12.

12. VanRaden PM, O'Connell JR, Wiggans GR, Weigel KA: Genomic evaluations with many more genotypes. *Genet Sel Evol* 2011, **43**:10.

13. Sun C, Wu XL, Weigel KA, Rosa GJM, Bauck S, Woodward BW, Schnabel RD, Taylor JF, Gianola D: An ensemble-based approach to imputation of moderate-density genotypes for genomic selection with application to Angus cattle. *Genet Res (Camb)* 2012, **94**:133-150.

14. Mulder HA, Calus MPL, Druet T, Schrooten C: Imputation of genotypes with low-density chips and its effect on reliability of direct genomic values in Dutch Holstein cattle. *J Dairy Sci* 2012, **95**:876-889.

15. Calus MPL, Veerkamp RF, Mulder HA: Imputation of missing single nucleotide polymorphism genotypes using a multivariate mixed model framework. *J Anim Sci* 2011, **89**:2042-2049.

16. Rincon G, Weber KL, Van Eenennaam AL, Golden BL, Medrano JF: Hot topic: performance of bovine high-density genotyping platforms in Holsteins and Jerseys. *J Dairy Sci* 2011, **94**:6116-6121.

17. Boichard D, Maignel L, Verrier E: The value of using probabilities of gene origin to measure genetic variability in a population. *Genet Sel Evol* 1997, **29**:5-23.

18. Boichard D: Pedig: a fortran package for pedigree analysis suited for large populations. In *Proceedings of the 7th World Congress on Genetics applied to Livestock Production (WCGALP):19-23 August 2002; Montpellier*. CD-ROM communication No. 28-13; 2002.

19. Lund MS, de Roos APW, de Vries AG, Druet T, Ducrocq V, Fritz S, Guillaume F, Guldbrandtsen B, Liu Z, Reents R, Schrooten C, Seefried F, Su G: A common reference population from four European Holstein populations increases reliability of genomic predictions. *Genet Sel Evol* 2011, **43**:43.

20. Su G, Brøndum RF, Ma P, Guldbrandtsen B, Aamand GP, Lund MS: Comparison of genomic predictions using medium-density (approximately 54,000) and high-density (approximately 777,000) single nucleotide polymorphism marker panels in Nordic Holstein and Red Dairy Cattle populations. *J Dairy Sci* 2012, **95**:4657-4665.

21. Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, Heaton MP, O'Connell J, Moore SS, Smith TPL, Sonstegard TS, Van Tassell CP: Development and characterization of a high density SNP genotyping assay for cattle. *PLoS One* 2009, 4:e5350.

22. Zhang Z, Druet T: Marker imputation with low-density marker panels in Dutch Holstein cattle. *J Dairy Sci* 2010, **93**:5487-5494.

23. Danchin-Burge C: Estimation de la variabilité génétique de 19 races bovines à partir de leurs généalogies. Paris: Institut de l'Elevage; 2009.

24. Hill WG, Robertson A: Linkage disequilibrium in finite populations. *Theor Appl Genet* 1968, **38**:226-231.

25. Goddard KAB, Hopkins PJ, Hall JM, Witte JS: Linkage disequilibrium and allele-frequency distributions for 114 single-nucleotide polymorphisms in five populations. *Am J Hum Genet* 2000, **66**:216-234.

26. Reich DE, Cargill M, Bolk S, Ireland J, Sabeti PC, Richter DJ, Lavery T, Kouyoumjian R, Farhadian SF, Ward R, Lander ES: Linkage disequilibrium in the human genome. *Nature* 2001, **411**:199-204.

27. Gautier M, Laloë D, Moazami-Goudarzi K: Insights into the genetic history of French cattle from dense SNP data on 47 worldwide breeds. *PLoS One* 2010, **5**:e0013038.

28. Hayes BJ, Bowman PJ, Daetwyler HD, Kijas JW, Van der Werf JHJ: Accuracy of genotype imputation in sheep breeds. *Anim Genet* 2012, **43**:72-80.

29. Erbe M, Hayes BJ, Matukumalli LK, Goswami S, Bowman PJ, Reich CM, Mason BA, Goddard ME: 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-4129.

30. Dassonneville R, Brøndum RF, Druet T, Fritz S, Guillaume F, Guldbrandtsen B, Lund MS, Ducrocq V, Su G: Effect of imputing markers from a low-density chip on the reliability of genomic breeding values in Holstein populations. *J Dairy Sci* 2011, **94**:3679-3686.

doi:10.1186/1297-9686-45-33

Cite this article as: Hozé et al.: High-density marker imputation accuracy in sixteen French cattle breeds. *Genetics Selection Evolution* 2013 **45**:33.

Submit your next manuscript to BioMed Central and take full advantage of:

• Convenient online submission

• Thorough peer review

• No space constraints or color figure charges

• Immediate publication on acceptance

• Inclusion in PubMed, CAS, Scopus and Google Scholar

• Research which is freely available for redistribution

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central