Search

Although several studies have investigated the factors affecting imputation accuracy, most of these studies involved a large number of genotyped animals. Thus, results from these studies cannot be directly applied to small populations, since the population structure affects imputation accuracy. In addition, factors affecting imputation accuracy may also be intensified in small populations. Therefore, we aimed to compare different imputation strategies for the Portuguese Holstein cattle population considering several commercially available single nucleotide polymorphism (SNP) panels in a relatively small number of genotyped animals. Data from 1359 genotyped animals were used to evaluate imputation in 7 different scenarios. In the S1 to S6 scenarios, imputations were performed from LDv1, 50Kv1, 57K, 77K, HDv3 and Ax58K panels to 50Kv2 panel. In these scenarios, the bulls in 50Kv2 were divided into reference (352) and validation (101) populations based on the year of birth. In the S7 scenario, the validation population consisted of 566 cows genotyped with the Ax58K panel with their genotypes masked to LDv1. In general, all sample imputation accuracies were high with correlations ranging from 0.94 to 0.99 and concordance rate ranging from 92.59 to 98.18%. SNP-specific accuracy was consistent with that of sample imputation. S4 (40.32% of SNPs imputed) had higher accuracy than S2 and S3, both with less than 7.59% of SNPs imputed. Most probably, this was due to the high number of imputed SNPs with minor allele frequency (MAF) < 0.05 in S2 and S3 (by 18.43% and 16.06% higher than in S4, respectively). Therefore, for these two scenarios, MAF was more relevant than the panel density. These results suggest that genotype imputation using several commercially available SNP panels is feasible for the Portuguese national genomic evaluation.

The objectives of this study were to estimate linkage disequilibrium (LD), describe and scan a haplotype block for the presence of genes that may affect milk production traits in Portuguese Holstein cattle. Totally 526 animals were genotyped using the Illumina BovineSNP50 BeadChip, which contained a total of 52 890 single nucleotide polymorphisms (SNPs). The final set of markers remaining after considering quality control standards consisted of 37 031 SNPs located on 29 autosomes. The LD parameters historical recombinations through allelic association (D') and squared correlation coefficient between locus alleles frequencies ( r2) were estimated and haplotype block analyses were performed using the Haploview software. The averages of D' and r2 values were 0.628 and 0.122, respectively. The LD value decreased with increasing physical distance. The D' and r2 values decreased respectively from 0.815 and 0.283 at the distance of 0-30 kb to 0.578 and 0.090 at the distance of 401-500 kb. The identified total number of blocks was 969 and consisted of 4259 SNPs that covered 159.06 Mb (6.24% of the total genome) on 29 autosomes. Several genes inside the haplotype blocks were detected; CSN1S2 gene in haplotype block 51 on BTA 6, IL6 and B4GALT1 genes in haplotype blocks 6 and 33 on BTA 8, IL1B and ID2 genes in haplotype blocks 19 and 29 on BTA 11, and DGAT1 gene in haplotype block 1 on BTA 14. The extension of LD using BovineSNP50 BeadChip did not exceed 500 kb and its parameters r2 and D' were less than 0.2 and 0.70, respectively, after 70-100 kb. Consequently, the 50K BeadChip would have a poor power in genome wide association studies at distances between adjacent markers lower than 70 kb.