《Ecology and Evolution》：Unraveling the Evolutionary History of Snow Sheep (Ovis nivicola): Genome-Wide Single Nucleotide Polymorphism Analysis Reveals Genetic Diversity and Clarifies Taxonomy
编辑推荐：
雪羊（Ovis nivicola）是一种栖息于西伯利亚东北部和俄罗斯远东地区高地的野生山地偶蹄目动物。尽管具有生态和保护上的重要性，但由于其分布广泛、形态变异大以及可用材料有限，该物种的种内分类学长期以来一直存在争议。在本研究中，研究人员进行了全面的全基因组单
雪羊（Ovis nivicola）是一种栖息于西伯利亚东北部和俄罗斯远东地区高地的野生山地偶蹄目动物。尽管具有生态和保护上的重要性，但由于其分布广泛、形态变异大以及可用材料有限，该物种的种内分类学长期以来一直存在争议。在本研究中，研究人员进行了全面的全基因组单核苷酸多态性（SNP）分析，涵盖所有已描述的主要种群，以阐明雪羊的进化关系和分类结构。共有来自15个种群的337个个体使用Illumina OvineSNP50和OvineHD 600 K BeadChips进行了基因分型，经过质量控制后保留了1032个常染色体SNP。种群结构分析揭示了八个主要的遗传聚类，它们在很大程度上对应于地理区域：堪察加、科里亚克、鄂霍次克/楚科奇、雅库特（Moma Range）、雅库特（Verkhoyansk Range）、雅布洛诺夫、哈劳拉赫和普托拉纳。成对FST值证实了大多数种群之间存在显著的分化。利用f3统计量，研究人员鉴定出两次显著的基因渗入事件，一次涉及Orulgan种群，另一次涉及堪察加北部种群。基于综合的基因组和历史证据，研究结果支持对堪察加（O. n. nivicola）、普托拉纳（O. n. borealis）、鄂霍次克（O. n. alleni）、科里亚克（O. n. koriakorum）、雅库特（O. n. lydekkeri）和斯坦诺夫/雅布洛诺夫（O. n. potanini）亚种的分类学确认，并提出将哈劳拉赫雪羊确立为一个新亚种，即O. n. ernsti。Kodar种群的特征是杂合度极低，需要进一步采样以解决其分类地位。这些发现为理解雪羊的进化史提供了首个全基因组框架，并为未来的保护管理确立了遗传学基础。
雪羊（Ovis nivicola）是一种分布于西伯利亚东北部和俄罗斯远东地区高山系统的野生偶蹄目动物，在维持高海拔生态系统平衡中发挥着重要的生态作用，并具有极高的保护价值。在系统发育学上，雪羊与北美洲的野生羊（即大角羊Ovis canadensis和白大角羊Ovis dalli）亲缘关系密切，共同被归类于Pachyceros亚属。尽管Ovis属内的种间关系已经相对明确，但雪羊的种内分类学却一直存在争议和模糊之处，特别是关于已描述亚种的数量和有效性长期以来备受辩论。历史上，形态学特征和栖息地是雪羊分类的主要标准，但由于雪羊分布偏远且可用样本极其稀缺，早期的分类学描述常常仅基于单一标本，并且传统形态学分类易受气候、疾病及营养等环境因素的影响，导致遗传相关的个体可能表现出截然不同的表型，从而使可靠的分类学划分变得复杂。
为了解决雪羊种内分类的不确定性并在基因组水平评估其遗传多样性，研究人员开展了一项基于全基因组单核苷酸多态性（SNP）分析的研究。该研究系统地评估了历史上提出的所有雪羊亚种，结合基因组数据为长期的分类学争论提供了明确的诊断结论。研究得出的结论支持将堪察加、普托拉纳、鄂霍次克、科里亚克、雅库特以及斯坦诺夫/雅布洛诺夫种群确认为有效亚种，并提议将哈劳拉赫种群确立为一个新亚种。该研究的重要意义在于为理解雪羊的进化历史建立了首个全基因组框架，并为未来的保护管理规划奠定了坚实的遗传学基础。该论文已发表在学术期刊《Ecology and Evolution》上。
在技术方法上，研究人员从15个主要雪羊种群中共采集了337个样本。基因组DNA提取后，样本使用了两种Illumina绵羊基因分型芯片（OvineSNP50与OvineHD 600 K）进行全基因组SNP基因分型。经过严格的质量过滤，去除了偏离哈代-温伯格平衡及连锁不平衡的位点后，保留了1032个高质量的常染色体SNPs用于后续分析。研究采用了主成分分析（PCA）、ADMIXTURE ancestry分析程序、成对FST遗传距离计算、邻接网算法（Neighbor-Net）最大似然系统发育树重构以及f3统计量分析等多种群体遗传学统计算法，对种群结构、遗传分化、历史基因流及基因渗入事件进行了全面深入地推断。
在研究结果部分，研究揭示了以下核心发现：
1. 遗传种群结构与基因渗入
通过主成分分析（PCA）和祖先成分分析，研究证实所有个体均按其来源种群聚类，并发现了8个对应于地理区域的遗传聚类。进一步分析表明，Orulgan种群和Kamchatka-N种群存在显著的基因渗入，f3统计量证实了这两个混杂事件，分别来源于Kharaulakh与Verkhoyansk谱系、以及Kamchatka-S与Koryak谱系的遗传贡献。
2. 遗传分化与进化史
通过成对FST遗传距离分析，研究发现大多数种群具有显著的遗传分化水平。构建的邻接网络和最大似然系统发育树显示，种群的分支顺序基本符合地理分布特征，并通过TreeMix重构检测到两次主要的历史基因流事件，揭示了Koryak群体的复杂起源和Kamchatka群体最长的独立进化历史。
3. 种群内遗传多样性
通过评估各群体的观察杂合度（HO）、稀有等位基因丰富度（AR）及多位点杂合度（MLH）等指标，研究发现雅库特中心区域的种群遗传多样性最高，反映了历史上连通性较好带来的基因交流。相反，被超过1000公里距离孤立的Putorana群体和样本极度稀缺的Kodar群体表现出极低的杂合度。
在讨论与结论部分，研究人员结合历史分类学和基因组证据，对各个亚种的有效性进行了全面澄清。研究支持了Kamchatkan（O. n. nivicola）等六个传统亚种的有效分类地位，提出Chukotkan雪羊在遗传上与Okhotsk亚种高度相似，应归为O. n. alleni。鉴于极高的遗传分化，研究人员正式提议将Kharaulakh雪羊确立为新亚种O. n. ernsti。对于Kodar雪羊，因仅有的单个样本呈现极低杂合度且受限于样本稀缺，其确切的分类地位仍需后续研究确证。最后，研究指出本次当代遗传多样性分布模式主要受更新世冰期避难所历史及高山系统长期地理隔离的联合效应所塑造，而非单一地理距离所能解释，为羊亚科物种的保护与进化研究提供了新视角。
