Department of Equine Sciences, Sorabol College, Gyeongju, 780-711 Republic of Korea.
The Animal Genomics and Breeding Center, Hankyong National University, Anseong, 456-749 Republic of Korea.
J Anim Sci Technol. 2014 Dec 5;56:28. doi: 10.1186/2055-0391-56-28. eCollection 2014.
This study was conducted to estimate the effective population size using SNPs data of 240 Jeju horses that had raced at the Jeju racing park. Of the total 61,746 genotyped autosomal SNPs, 17,320 (28.1%) SNPs (missing genotype rate of >10%, minor allele frequency of <0.05 and Hardy-Weinberg equilibrium test P-value of <10(-6)) were excluded after quality control processes. SNPs on the X and Y chromosomes and genotyped individuals with missing genotype rate over 10% were also excluded, and finally, 44,426 (71.9%) SNPs were selected and used for the analysis. The measures of the LD, square of correlation coefficient (r(2)) between SNP pairs, were calculated for each allele and the effective population size was determined based on r(2) measures. The polymorphism information contents (PIC) and expected heterozygosity (HE) were 0.27 and 0.34, respectively. In LD, the most rapid decline was observed over the first 1 Mb. But r(2) decreased more slowly with increasing distance and was constant after 2 Mb of distance and the decline was almost linear with log-transformed distance. The average r(2) between adjacent SNP pairs ranged from 0.20 to 0.31 in each chromosome and whole average was 0.26, while the whole average r(2) between all SNP pairs was 0.02. We observed an initial pattern of decreasing Ne and estimated values were closer to 41 at 1 ~ 5 generations ago. The effective population size (41 heads) estimated in this study seems to be large considering Jeju horse's population size (about 2,000 heads), but it should be interpreted with caution because of the technical limitations of the methods and sample size.
本研究旨在利用在济州赛马公园参赛的 240 匹济州马的 SNP 数据估算有效种群大小。在总共 61746 个全基因组 SNP 中,经过质量控制后,有 17320 个 SNP(缺失基因型率>10%、次要等位基因频率<0.05 和 Hardy-Weinberg 平衡检验 P 值<10(-6))被排除在外。X 和 Y 染色体上的 SNP 以及缺失基因型率超过 10%的个体也被排除在外,最终选择了 44426 个(71.9%)SNP 用于分析。计算了每个等位基因之间 SNP 对的 LD 度量,即相关系数平方(r(2)),并根据 r(2)度量确定了有效种群大小。多态信息含量(PIC)和期望杂合度(HE)分别为 0.27 和 0.34。在 LD 中,观察到第一个 1 Mb 内最快的下降。但是 r(2)随距离的增加下降较慢,在 2 Mb 距离后保持恒定,下降几乎呈对数变换距离的线性关系。每个染色体和整个平均的相邻 SNP 对之间的平均 r(2)范围为 0.20 到 0.31,而所有 SNP 对之间的整个平均 r(2)为 0.02。我们观察到 Ne 呈下降趋势,估计值在 1 到 5 代前更接近 41。考虑到济州马的种群数量(约 2000 头),本研究中估计的有效种群大小(41 头)似乎很大,但由于方法和样本量的技术限制,应谨慎解释。