Makina Sithembile O, Taylor Jeremy F, van Marle-Köster Este, Muchadeyi Farai C, Makgahlela Mahlako L, MacNeil Michael D, Maiwashe Azwihangwisi
Agricultural Research Council-Animal Production Institute Pretoria, South Africa ; Department of Animal and Wildlife Sciences, University of Pretoria Pretoria, South Africa.
Division of Animal Sciences, University of Missouri Columbia, MO, USA.
Front Genet. 2015 Dec 1;6:337. doi: 10.3389/fgene.2015.00337. eCollection 2015.
Knowledge on the extent of linkage disequilibrium (LD) in livestock populations is essential to determine the minimum distance between markers required for effective coverage when conducting genome-wide association studies (GWAS). This study evaluated the extent of LD, persistence of allelic phase and effective population size (Ne) for four Sanga cattle breeds in South Africa including the Afrikaner (n = 44), Nguni (n = 54), Drakensberger (n = 47), and Bonsmara breeds (n = 46), using Angus (n = 31) and Holstein (n = 29) as reference populations. We found that moderate LD extends up to inter-marker distances of 40-60 kb in Angus (0.21) and Holstein (0.21) and up to 100 kb in Afrikaner (0.20). This suggests that genomic selection and association studies performed within these breeds using an average inter-marker r (2)≥ 0.20 would require about 30,000-50,000 SNPs. However, r (2)≥ 0.20 extended only up to 10-20 kb in the Nguni and Drakensberger and 20-40 kb in the Bonsmara indicating that 75,000 to 150,000 SNPs would be necessary for GWAS in these breeds. Correlation between alleles at contiguous loci indicated that phase was not strongly preserved between breeds. This suggests the need for breed-specific reference populations in which a much greater density of markers should be scored to identify breed specific haplotypes which may then be imputed into multi-breed commercial populations. Analysis of effective population size based on the extent of LD, revealed Ne = 95 (Nguni), Ne = 87 (Drakensberger), Ne = 77 (Bonsmara), and Ne = 41 (Afrikaner). Results of this study form the basis for implementation of genomic selection programs in the Sanga breeds of South Africa.
了解家畜群体中的连锁不平衡(LD)程度对于在进行全基因组关联研究(GWAS)时确定有效覆盖所需标记之间的最小距离至关重要。本研究评估了南非四个桑加牛品种(包括阿非利卡牛(n = 44)、恩古尼牛(n = 54)、德拉肯斯伯格牛(n = 47)和邦斯玛拉牛(n = 46))的LD程度、等位基因相位的持续性和有效种群大小(Ne),以安格斯牛(n = 31)和荷斯坦牛(n = 29)作为参考群体。我们发现,在安格斯牛(r² = 0.21)和荷斯坦牛(r² = 0.21)中,中等程度的LD延伸至标记间距离40 - 60 kb,在阿非利卡牛中(r² = 0.20)延伸至100 kb。这表明在这些品种中使用平均标记间r²≥0.20进行基因组选择和关联研究将需要约30,000 - 50,000个单核苷酸多态性(SNP)。然而,在恩古尼牛和德拉肯斯伯格牛中,r²≥0.20仅延伸至10 - 20 kb,在邦斯玛拉牛中延伸至20 - 40 kb,这表明在这些品种中进行GWAS需要75,000至150,000个SNP。相邻位点等位基因之间的相关性表明,各品种之间的相位保存不强。这表明需要特定品种的参考群体,其中应标记更高密度的标记以识别特定品种的单倍型,然后可将其推算到多品种商业群体中。基于LD程度对有效种群大小的分析显示,恩古尼牛Ne = 95,德拉肯斯伯格牛Ne = 87,邦斯玛拉牛Ne = 77,阿非利卡牛Ne = 41。本研究结果为南非桑加牛品种实施基因组选择计划奠定了基础。