全基因组分析揭示濒危哥伦比亚盆地侏兔种群的遗传独特性

为评估遗传拯救对濒危物种的影响，研究人员利用RADseq技术对哥伦比亚盆地(Columbia Basin, CB)侏兔的祖先种群开展了全基因组分析。研究明确了该物种范围内的四个主要遗传类群，证实CB种群具有最高的遗传分化水平，并鉴定了可用于监测混合种群祖源组成的SNP位点，为濒危种群的管理和遗传监测提供了关键基因组学工具。

在美国西部辽阔的山艾树草原上，生活着世界上最小的兔子——侏兔。这种小生灵的命运与其赖以为生的山艾树栖息地紧密相连，然而，由于土地利用变化、入侵物种和火灾等因素，超过50%的栖息地已经丧失。栖息地的丧失和破碎化，对位于物种分布区边缘的种群打击尤为沉重，华盛顿州哥伦比亚盆地的侏兔种群就是其中最岌岌可危的代表。这个种群在2001年野外仅存16只个体，并于2003年被列为濒危物种。为了挽救这个种群，保护工作者们启动了一项雄心勃勃的“遗传拯救”计划，从爱达荷州等其他种群引入个体，与本地个体进行杂交，以期通过增加遗传多样性来扭转其衰退的命运。然而，一个重要的问题随之浮现：在进行遗传拯救和种群混合之前，哥伦比亚盆地的原始祖先种群，其遗传特征究竟是什么？它在整个侏兔物种的遗传谱系中占据着怎样的独特地位？为了回答这些问题，也为未来的保护监测提供精准的分子工具，一项覆盖整个物种分布范围的基因组学研究应运而生。

为了深入探究侏兔的种群遗传历史，Nerkowski, S. A., Waits, L. P., Warheit, K. I., Bacchiocchi, I. 和 Hohenlohe, P. A. 所领导的研究团队采用了限制性位点相关DNA测序技术。他们收集了来自美国西部多个州（包括华盛顿、爱达荷、蒙大拿、怀俄明、犹他、内华达、俄勒冈和加利福尼亚）的侏兔样本，其中最关键的是在遗传拯救和混合发生前采集的、具有纯正哥伦比亚盆地血统的个体样本。通过提取基因组DNA、构建RADseq文库并在Illumina平台上进行测序，研究人员获得了海量的遗传数据。随后，他们利用生物信息学流程（包括Stacks软件进行基因分型，Bowtie2将序列比对到欧洲兔参考基因组）对数据进行处理，最终从123个个体中鉴定出9794个高质量的单核苷酸多态性位点，构建了用于后续分析的数据集。这套数据集为全面评估种群遗传结构、分化程度和潜在的适应性变异奠定了基础。

遗传结构与多样性

研究人员采用了多种分析方法来解析侏兔的种群遗传结构。主成分分析结果显示，哥伦比亚盆地种群沿着解释25%变异的第一主成分轴与其他所有种群明显分开。基于最大似然法的群体结构分析和稀疏非负矩阵分解算法均一致支持将整个物种范围内的侏兔划分为四个主要的遗传类群：1) 哥伦比亚盆地；2) 犹他州北部/怀俄明州；3) 犹他州南部；4) 物种分布区的中心区域（包括内华达、俄勒冈、爱达荷、蒙大拿和加利福尼亚）。在中心区域内部，斯内克河似乎构成了一道基因流屏障，将爱达荷州北部/蒙大拿州的种群与南部的大盆地种群区分开来。遗传多样性分析显示，哥伦比亚盆地种群的核苷酸多样性（π = 0.07）和观测杂合度（H_o= 0.10）均为最低，而中心区域种群的遗传多样性最高。成对遗传分化分析进一步证实，哥伦比亚盆地种群与其他三个类群（特别是犹他州南部、犹他州北部/怀俄明州类群）之间的遗传分化系数最高（F_ST约为0.52-0.53），凸显了其极端的遗传独特性。

适应性位点与基因本体

为了探究不同种群间是否存在与局部适应相关的遗传变异，研究人员使用了pcadapt、OutFLANK和BayeScan三种方法搜寻受到歧化选择的“离群”位点。结果在pcadapt分析中鉴定出827个推定具有适应性的SNP位点，其中约46.7%位于欧洲兔参考基因组已注释的基因内。这些位点广泛分布在各个染色体上，并未聚集在某个特定区域−8, while the blue line corresponds to the suggestive threshold of p = 1.00 × 10⁻⁵. Alternating blue and orange colors represent chromosomes, using alignment to the European rabbit reference genome assembly.">，这表明种群间的适应性分化可能是多基因控制的。然而，针对这些候选适应性位点进行的基因本体富集分析并未发现任何显著富集的功能类别，意味着未能识别出与特定生物学过程相关的强烈信号。

综合研究结果，该项研究得出了几个关键结论。首先，全基因组数据清晰地勾勒出侏兔物种范围内存在四个主要的进化谱系，其中濒危的哥伦比亚盆地种群显示出最高程度的遗传独特性和分化水平，这从基因组层面支持了其当初被列为濒危物种的遗传学依据。其次，哥伦比亚盆地种群的遗传多样性极低，这与其长期的遗传隔离和近代严重的种群瓶颈效应密切相关。第三，尽管鉴定出数百个潜在的适应性位点，但未发现显著富集的基因功能类别，暗示局部适应可能涉及复杂的多基因机制，或许与应对不同地区山艾树中次生代谢物的化学差异有关。研究人员推测，更新世冰期的冰川、湖泊（如邦纳维尔湖、拉洪坦湖）以及历史栖息地的分布，可能是塑造当前种群遗传格局的重要历史因素。

这项发表于《Genes》的研究具有重要的保护生物学意义。它首次利用全基因组范围的单核苷酸多态性标记，描绘了侏兔的种群遗传图谱，并成功刻画了其祖先哥伦比亚盆地种群的遗传“基线”。研究所鉴定出的特异性单核苷酸多态性位点，为未来设计高通量基因分型监测面板（如GT-Seq或SNP芯片）提供了宝贵的资源，使得保护管理者能够长期、精准地追踪混合种群中不同祖源成分的比例、评估遗传拯救的效果、甚至探索在哥伦比亚盆地环境下可能受选择的适应性基因渗入。此外，研究结果也提示，不仅哥伦比亚盆地，其他如犹他州南部、北部/怀俄明以及加利福尼亚莫诺盆地等边缘种群，也同样表现出相当的遗传独特性，这为整个物种的保育单元划分和优先保护策略的制定提供了关键的基因组学依据。这项研究不仅深化了我们对一个濒危物种演化历史的理解，更将基因组学工具实实在在地应用于濒危物种的保护实践之中。