哈萨克羊三个品种的全基因组SNP特征分析：哈萨克脂尾粗毛羊、杰格里斯羊与埃蒂美利奴羊

《Animals》：Genome-Wide SNP Characterisation of Three Kazakh Sheep Breeds: Kazakh Fat-Tailed Coarse-Wool, Degeres, and Etti Merino

　　

哈萨克斯坦的羊产业在国家羊肉与羊毛生产中占据核心地位，但若干主要品种此前从未在DNA层面被表征。本研究旨在对三个具有经济重要性的哈萨克品种——杰格里斯羊（DE）、埃蒂美利奴羊（EM）与哈萨克脂尾粗毛羊（KKG）——进行基因组多样性、群体结构及全球系统发育定位的表征，并判断它们是否构成需要独立管理的不同基因池。研究人员呈现了首项同时纳入DE、EM与KKG的全基因组SNP表征。来自七个生产畜群的1497只动物（DE=354，EM=642，KKG=501）经基因分型，质控后于42279个SNP（其中22766个LD剪枝标记用于主成分分析与AMOVA）上分析。应用主成分分析（PCA）、两两FST、分子方差分析（AMOVA）、邻接系统发育、基于模型的祖源估计（ADMIXTURE）及Hill数多样性分析，并将这些品种投影至全球Ovine SNP50 HapMap面板（74个参考品种，2819只动物，37685个共享SNP）。三个品种均维持均匀的高品种内多样性（期望杂合度0.413–0.417），固定指数近零。AMOVA将94.03%方差分配于品种内（ΦST=0.060，p<0.001）。PCA、系统发育与ADMIXTURE在K=3时一致解析出三个品种特异性簇，数据集内品种间最大FST为0.038。在全球面板中，EM在遗传上最接近美利奴及美利奴衍生参考品种（合并FST=0.017），而距西南亚羊群较远（FST=0.045）；DE与KKG呈相反格局（距西南亚FST=0.027和0.020，距美利奴组0.052）。DE额外表现出杂合子过剩与部分混合，符合其未完全固定的合成品种身份。结果勾勒出三个不同基因池，对哈萨克羊品种管理与基因组多样性保护具直接意义。

研究背景与立项依据

哈萨克斯坦地处欧亚草原与丝绸之路要冲，羊业支撑全国大部分羊肉与羊毛产出。三个经济品种各具代表性：哈萨克脂尾粗毛羊（KKG）为古老土著脂尾粗毛型，适应半荒漠极端气候；杰格里斯羊（DE）为近期由哈萨克脂尾羊与引进品种（如Shropshire、Precoce）杂交并“自交巩固”的肉毛兼用合成种；埃蒂美利奴羊（EM）为2011年用德国肉用美利奴（Merinofleischschaf）公羊与哈萨克细毛母羊育成的肉用细毛种。尽管三者农业地位关键，此前微卫星研究仅触及部分品种，SNP芯片研究未同时纳入DE与EM，更未将美利奴衍生种与土著脂尾种置于同一框架比较，故“三者是否为可区分基因池、EM是否仍带美利奴基因组印记、DE是否仍未巩固”均属空白。研究人员因此设计全基因组SNP研究，为品种管理提供基线。

主要技术方法概要

研究人员从哈国七家牧场采集1497只12–13月龄未产母羊与2–5岁公羊的毛囊DNA，用OvineSNP50v3_HTS芯片分型；质控后保留42279个SNP（LD剪枝至22766个用于PCA/AMOVA）。核心算法包括：SNPRelate做PCA与两两F_ST（Weir-Cockerham估计）；pegas包做AMOVA（1000次置换）；ape与ggtree构邻接（NJ）树；ADMIXTURE无监督（K=2–8）与有监督（K=6，固定全球参考群）祖源分析；dartR算Hill数（q=0,1,2）。全球定位阶段将数据与ISGC OvineSNP50 HapMap面板（74品种2819只）合并，保留37685共享SNP重做PCA/NJ/有监督ADMIXTURE。

研究结果

3.1 群体结构（Population Structure）

PCA显示PC1（3.2%）完全分离EM与另两品种，PC2（0.9%）进一步分DE/KKG。两两F_ST在-0.0008至0.0407，跨品种家庭对最高。AMOVA表明94.03%方差在品种内（Φ_ST=0.060），94.42%在家庭内。NJ树与PCA一致成三簇。Hill数各家庭q=0≈2.0、q=1≈1.80–1.84、q=2≈1.73–1.75，无显著差异。结论：三品种内部多样性高、间分化浅但可辨。

3.2 全球群体结构（Global Population Structure）

投影至全球面板后，EM在PCA中紧靠美洲/澳洲美利奴组，DE与KKG邻西南/东亚/南亚组。F_ST量化：EM对美利奴组0.017、对西南亚0.045；DE/KKG对西南亚0.027/0.020、对美利奴组0.052。NJ种群树（补充）显示EM姊妹群为中国美利奴与Rambouillet，DE/KKG互为姊妹嵌于亚洲羊群。结论：EM保留强美利奴基因组签名，DE/KKG归亚洲脂尾-合成集群。

3.3 祖源估计（Ancestry Estimation）

无监督ADMIXTURE在K=3时ΔK峰值，EM与KKG各占主导成分，DE整体携显著KKG关联成分（尤Akterek家庭），印证其未完全巩固。有监督K=6分析显示EM 77.6%美洲美利奴成分、亚洲仅14.8%；DE 95.8%亚洲、残留4.1%美利奴（来自Precoce）；KKG 99.3%亚洲、美利奴0.6%。结论：三方法互证EM美利奴导向、DE部分混合、KKG纯亚洲土著。

讨论部分总结

研究人员指出，三品种同框全基因组表征为首报。EM虽在哈国土培育超十年，仍持美利奴签名（F_ST对美利奴组0.017），与全球美利奴衍生种“祖先持久性”规律吻合；DE杂合子过剩（F=-0.0123）、跨家庭分化、携KKG成分与残美利奴成分，符合“未巩固合成种”典型信号；KKG与DE归亚洲群，呼应青铜时代西草原祖先流入后的中亚共享背景。三品种期望杂合度0.413–0.417、F近零、AMOVA品种间仅5.97%，并非品种界限模糊，而是绵羊基因组普遍“品种内方差主导”常态，不影响三池可分。管理上应：分池记录选育；EM借用美利奴参考群与选种指数，DE/KKG需建中亚脂尾参考；94%方差在品种内说明育种增益靠品种内选种而非继续杂交；DE需核心群固定。

结论部分翻译

本研究给出DE、EM、KKG首次同框全基因组表征，得出三结论直接服务品种管理：EM保留足够强的美利奴基因组签名，使其在合并全球面板中比对哈萨克同场品种（F_ST=0.041）或西南亚地理邻群（F_ST=0.045）更接近美利奴及美利奴衍生参考种（合并F_ST=0.017），尽管本地育种已逾十年；DE携带未完全巩固合成种的多重签名（杂合子过剩、K=3时非DE成分、家庭间分化），与其文献记载多品种起源一致；三品种均维持均匀高品种内多样性（H_e=0.413–0.417），绝大多数方差在品种内（Φ_ST=0.060）或家庭内（Φ_ST=0.056），均无近交或瓶颈迹象。对哈萨克羊群结构有四重实操含义：三品种作独立种群记录选育；EM用美利奴型选种体系、DE/KKG用中亚脂尾参考基；约94%遗传方差在品种内，育种进步靠品种内选种且可提强度而近交风险低；DE最需巩固，设核心群比现行“自交”更快固定品种。高多样性基线亦使三者均可直接做选择信号作图，尤其围绕KKG亲属已刻画的脂尾脂肪生成位点（BMP2、PDGFD、VRTN）。论文发表于《Animals》。