Effect of SNP origin on analyses of genetic diversity in cattle

Abstract

The methods of single nucleotide polymorphism (SNP) identification can lead to ascertainment bias, which will affect population genetic analyses based on those data. In livestock species, the methods of SNP identification through genome sequencing are likely to suffer from this ascertainment bias. In the present study, a subset of data from the Bovine HapMap Project was re-analysed to quantify the effects of ascertainment bias on a range of common analyses and statistics. Data from 189 animals of the zebu breeds Brahman, Nelore and Gir, taurine beef Angus, Limousin and Hereford and taurine dairy Holstein, Jersey and Brown Swiss were analysed. There were 141 SNPs each of Angus, Brahman and Holstein origin, giving a total of 423 SNPs organised in 141 triplets. Each triplet consisted of one SNP of each breed, separated on average by 0.75 Mb within each triplet and where triplets were separated by 14.96 Mb to ensure that each triplet was unaffected by linkage disequilibrium. The minor allele frequency distribution, estimates of the F-statistic, FST, the partitioning of variance and population substructure were relatively unaffected by breed of origin of the SNPs. Estimates of heterozygosity were significantly affected by breed of origin of the SNPs. The clustering of animals of closely related breeds varied in the principal component analyses (PCA). However, in the PCA the effect of breed of origin of 141 SNPs was similar to the effect of using different panels of 141 SNPs of all three breeds, so the differences found in the PCA may not be all due to bias by the origin of the SNPs. Based on these results, analyses that depend on FST, including signatures of selection, gene flow and effective population size are unlikely to be strongly affected by SNP origin. Analyses that partition genetic variance and some analyses of population substructure will also be largely unaffected. However, analyses that are dependent on locus heterozygosity, which can be used for studying population bottlenecks, or those that study selection using extended haplotype homozygosity may be significantly affected by breed of origin of the SNPs.