揭秘埃塞俄比亚高地特有物种狮尾狒的遗传多样性：景观地形如何塑造其种群结构与保护格局

为了揭示埃塞俄比亚高原特有灵长类——狮尾狒（Theropithecus gelada）的遗传变异格局以及高原关键地理特征如何影响其种群连通性，研究者扩大了此前有限的中央狮尾狒基因组数据，对来自四个地点的十个中央狮尾狒个体和两个南方狮尾狒进行了全基因组测序。研究发现，结合海拔偏好的景观模型比单纯基于直线地理距离的模型能更好地解释观测到的基因组距离，表明海拔结构强烈制约了其扩散。该研究明确了高原地形在塑造狮尾狒种群结构中的关键作用，为埃塞俄比亚灵长类的研究与保护规划提供了重要信息。

在非洲大陆的屋脊——埃塞俄比亚高原上，生活着一种独特而迷人的灵长类动物：狮尾狒（Theropithecus gelada）。它们是目前仅存的狮尾狒属物种，凭借对高海拔稀薄空气的卓越适应力，在1800至4500米的崎岖高原上开辟了自己的生存空间。白天，它们在高原草甸上采食禾草，夜晚则撤退到陡峭的崖壁缝隙中躲避捕食者，这种对特殊“悬崖-草原”复合生境的依赖，构成了它们生存的基本生态背景。长期以来，研究者们根据地理分布和形态差异，将狮尾狒划分为三个亚种：北方狮尾狒（T. g. gelada）、中央狮尾狒（T. g. obscurus）和南方狮尾狒（T. g. arsi）。然而，在这些大致分类之下，尤其是对于分布最广的中央狮尾狒亚种，其内部真实的基因组多样性、种群结构以及塑造这些格局的生态和地理力量，仍然笼罩在重重迷雾之中。这团迷雾不仅阻碍了我们理解这种独特灵长类的演化历史，也为其科学保护带来了挑战。由于此前中央狮尾狒的基因组数据仅来自单一地点的三个个体，这种关键信息的缺失，使得我们无法精确评估不同地理屏障（如河流、深谷）和生境梯度如何影响了种群间的基因流动，从而制约了制定更有效、更具针对性的保护策略。为了拨开迷雾，一支国际研究团队开展了一项深入的基因组学研究，其成果发表于《International Journal of Primatology》。

为了解答上述问题，研究者们综合运用了一系列基因组学和生物信息学技术。他们从先前采集的血液样本中，选取了来自四个中埃塞俄比亚地点的十只中央狮尾狒和两只南方狮尾狒，利用Illumina NovaSeq 6000平台进行了高覆盖度的全基因组测序。随后，他们将新产生的数据与已发表的来自北方狮尾狒、其他中央狮尾狒及动物园个体的基因组数据整合，构建了包含84个个体的数据集。通过BWA-MEM2将测序数据比对到狮尾狒参考基因组（Tgel_1.0），并利用GATK和ANGSD（针对低深度样本）进行变异检测。基于高质量的单核苷酸多态性位点，研究者进行了主成分分析、群体结构分析、线粒体基因组和Y染色体系统发育重建，以揭示种群分化和历史。他们还通过运行纯合性分析和成对连续马尔可夫链算法估算了近交水平和历史有效种群大小变化。最后，核心的创新在于将遗传数据与景观特征结合，研究者利用埃塞俄比亚高地的数字高程模型，构建了结合海拔偏好、距悬崖距离和距永久水体距离的成本表面模型，计算了种群间的最小成本路径，并将其与观测到的遗传距离进行比较，以检验景观地形对基因流的实际影响。

研究结果

• •
  中央狮尾狒的基因组多样性图谱
  通过对新测序的十二个中央和南方狮尾狒个体进行分析，并结合已有数据，研究首次相对全面地描绘了中央狮尾狒的遗传变异格局。主成分分析清晰地将所有个体划分为与三个亚种对应的三大集群，证实了宏观分类的遗传学基础。有趣的是，在中央狮尾狒群体内部，来自不同地理位置的个体（如Guassa、Goshmeda、Magdera、Debra Libanos）在主成分空间中没有形成明显分离的亚群，这暗示了在亚种水平之下，地理隔离可能尚未导致强烈的遗传分化。然而，线粒体DNA系统发育分析显示，之前研究中识别出的中央狮尾狒两个母系分支（C1和C2）在现有基因组数据中主要由C1代表，这表明当前的取样可能仍未完全捕捉到中央群体的全部线粒体多样性。
• •
  种群遗传结构与多样性
  对野生群体的遗传多样性分析显示，南方狮尾狒群体（T. g. arsi）的平均杂合度最低，而纯合性片段总长最高，这与它们种群规模最小、分布范围最狭窄且面临较高灭绝风险的现状相符。北方和中央狮尾狒群体显示出更高的遗传多样性。基于F_ST（群体间遗传分化系数）构建的系统树支持北方、中央、南方三个群体的划分，其中北方与中央群体的遗传距离相对较近，而南方群体则与前两者分化更为明显。此外，对动物园个体的分析发现，它们大多在遗传上聚类于中央狮尾狒，但拥有独特的Y染色体单倍型，提示圈养种群可能源于特定的、或已在野外稀少的遗传谱系。
• •
  历史种群动态
  通过成对连续马尔可夫链分析对历史有效种群大小进行回溯，研究发现所有狮尾狒群体在距今约10万至1万年前（晚更新世）都经历了显著的种群规模下降。这一下降趋势与已知的该时期非洲气候波动，特别是高原地区环境的干旱化和草地收缩时期相吻合，表明历史气候变化是塑造狮尾狒种群历史的重要驱动力。
• •
  景观特征对基因流的约束作用
  这是本研究的核心发现。研究者检验了经典的“距离隔离”模型，即遗传分化是否与地理直线距离成正比。然而，曼特尔检验显示，在狮尾狒中，单纯的直线地理距离并不能很好地解释观察到的遗传距离模式。于是，他们转向构建“景观遗传学”模型。考虑到狮尾狒对高海拔（偏好2000-3000米）和临近悬崖生境的特殊依赖，研究者创建了一个“成本表面”地图。在这张地图上，移动的“成本”或“阻力”不是距离，而是由海拔和距悬崖距离共同决定——偏离偏好海拔带或远离悬崖的区域，移动成本更高。当比较基于这种“成本距离”（即沿最小成本路径累积的阻力）的模型与基于直线距离的模型时，发现结合海拔偏好的成本距离模型与观测到的遗传距离矩阵（F_ST）拟合得显著更好。这一结果强有力地证明，埃塞俄比亚高原复杂的地形，特别是海拔梯度，是制约狮尾狒个体扩散和基因流动的关键景观过滤器，其影响力超过了单纯的平面地理间隔。
• •
  推断未取样种群的亲缘关系
  利用上述优化的景观连通性模型，研究者能够预测目前缺乏基因组数据的狮尾狒种群的遗传归属。例如，对于提格雷地区（Tigray）的狮尾狒种群，模型推断支持它们与北方狮尾狒（而非中央狮尾狒）具有更紧密的遗传联系。更重要的是，分析为北方与中央狮尾狒之间的分界提供了新见解：模型结果更支持一条以塔卡泽河-青尼罗河峡谷系统为标志的东西向分隔线，而不是简单的南北划分。这为厘清狮尾狒亚种间的进化关系提供了新的地理框架。

研究结论与意义

这项研究通过对狮尾狒，特别是此前数据匮乏的中央亚种，进行大规模基因组测序与深入分析，极大地丰富了我们对该物种遗传多样性的认知。其最关键的发现是，埃塞俄比亚高地的地形结构，尤其是海拔分布，是塑造狮尾狒当前种群遗传格局的首要景观力量。狮尾狒并非在二维平面上自由扩散，它们的移动和基因交流被限制在符合其特定海拔和生境偏好的“走廊”中。这一结论将景观遗传学的分析框架成功应用于高山灵长类，揭示了局地环境因子如何精细地调控物种水平的基因流动。

研究的发现具有多重重要意义。在理论层面，它深化了我们对高山特化物种演化与适应的理解，展示了如何将基因组数据与高分辨率环境变量结合，来解析种群历史与动态。在分类学上，研究为北方与中央狮尾狒之间基于河流峡谷的东西向分化提供了遗传学支持，可能对未来亚种界定产生影响。在保护实践上，这项工作的价值尤为突出。首先，它确认了南方狮尾狒极低的遗传多样性和高近交水平，为其被列为濒危物种（Endangered）提供了紧迫的基因组证据。其次，研究所构建的、结合了海拔和悬崖依赖性的景观连通性模型，可以直接转化为保护工具。保护规划者可以利用此模型，识别出对维持种群间基因流至关重要的生态走廊，或预测未来气候变化导致适宜栖息地收缩时，哪些种群可能最先面临隔离风险。最后，对圈养种群遗传独特性的认识，提示了对其进行专门遗传管理以保存潜在稀有基因库的必要性。总而言之，这项研究不仅填补了狮尾狒遗传学的关键空白，更重要的是，它提供了一套方法论范例，展示了如何利用基因组学和景观模型，为生活在复杂地形中的濒危物种制定前瞻性、基于科学证据的保护策略。