基于全基因组测序的圈养阿拉伯瞪羚三个谱系近交程度与基因组独特性比较分析

《Global Ecology and Conservation》：Whole-genome sequencing revealed inbreeding and genomic distinctiveness: A comparative analysis of three captive Gazella arabica lineages

　　

圈养于沙特阿拉伯King Khalid野生动物研究中心（KKWRC）的瞪羚包含阿拉伯瞪羚（Gazella arabica）的三个谱系，即大陆谱系、法拉桑谱系（G. a. farasani）和深灰谱系。研究人员对7只圈养个体进行了高覆盖度全基因组测序（WGS）并生成基于SNP的数据，以表征个体水平的基因组结构、杂合度和纯合片段（ROH），从而为上述谱系的遗传变异提供初步基线。经质量过滤后，保留2,180万个双等位基因SNP。主成分分析（PCA）和基于SNP的系统发育树将采样个体按管理群组区分。以Gazella cuvieri为外群构建的有根系统发育树将深灰瞪羚置于支持度极高的单系分支中，而大陆瞪羚与法拉桑瞪羚聚类关系更近。深灰瞪羚表现出最高的基于ROH的平均近交系数（FROH = 0.344），高于大陆瞪羚（0.103）和法拉桑瞪羚（0.045），且其ROH的数量和总长度也最大。在过滤后的多态性SNP集合中，深灰瞪羚的观测杂合度最低。描述性加权成对FST值在深灰瞪羚与大陆瞪羚之间最大（0.586），其次为深灰瞪羚与法拉桑瞪羚之间（0.513），而大陆瞪羚与法拉桑瞪羚之间的分化较低（0.157）。由于样本量小且非随机、管理历史不同以及使用了异源参考基因组，这些结果不能解决分类学、种群历史或物种边界问题。本研究数据为扩大采样、长期遗传监测及KKWRC未来育种管理提供了初步基因组基线和第一步基础。

论文解读

一、研究背景与科学问题

阿拉伯瞪羚（Gazella arabica）是阿拉伯半岛最具代表性的瞪羚物种之一，被国际自然保护联盟（IUCN）列为易危（Vulnerable）物种。长期以来，瞪羚属（Gazella）的分类学问题极为复杂，不同学者基于不同物种概念提出的物种数量从10种到23种不等，形态学和表型特征因个体差异大而难以作为可靠的分类依据。在沙特阿拉伯King Khalid野生动物研究中心（KKWRC），圈养着阿拉伯瞪羚的三个管理谱系：大陆谱系、法拉桑谱系（G. a. farasani）和深灰谱系。其中深灰谱系仅存27只个体（3只雌性、24只雄性），其地理来源不明，且此前基于线粒体和微卫星标记的研究未能明确其分类归属。法拉桑谱系虽然长期隔离于法拉桑群岛，但其遗传多样性状况尚不清楚。从圈养繁殖管理角度出发，阐明这些谱系的进化历史和遗传结构对指导保护工作至关重要。然而，针对这些圈养谱系的全基因组数据此前完全缺失。本研究旨在回答三个问题：这三个谱系是否具有遗传独特性；基于SNP的全基因组遗传多样性在各谱系间如何变化；各谱系在近交水平上是否存在差异。

二、研究内容与主要结论

研究人员对7只代表性个体（3只深灰瞪羚G1-G3、2只大陆瞪羚G4-G5、2只法拉桑瞪羚G6-G7）进行了高覆盖度全基因组测序，以达马瞪羚（Nanger dama）基因组为参考进行比对和变异检测，经严格质控后保留了21,852,137个高质量双等位基因SNP。主要研究结论如下：基于SNP的系统发育分析和主成分分析均清晰地将三个谱系区分为独立的遗传簇，深灰瞪羚形成支持度极高的单系分支，且与大陆和法拉桑谱系的分化程度远大于后两者之间的分化；深灰瞪羚表现出极高的近交水平（F_ROH = 0.344），远超大陆瞪羚（0.104）和法拉桑瞪羚（0.045）；深灰瞪羚的全基因组杂合度极低（约0.04-0.045），而法拉桑瞪羚最高（0.14-0.15），大陆瞪羚居中（0.12）；等位基因频率谱分析显示深灰瞪羚存在大量固定位点（74.6%），法拉桑瞪羚则以单拷贝变异为主（57.7%）；成对F_ST值表明深灰瞪羚与大陆瞪羚之间分化最大（加权值0.586），深灰瞪羚与法拉桑瞪羚之间次之（0.513），而大陆与法拉桑之间分化较低（0.157）。这些结果首次为KKWRC三个圈养谱系提供了个体水平的基因组基线，揭示了深灰瞪羚因长期封闭繁殖导致的严重遗传贫乏，为后续遗传管理和保护决策提供了重要科学依据。

三、关键技术与方法

本研究采用了以下核心技术方法：对7只圈养个体（3只深灰瞪羚、2只大陆瞪羚、2只法拉桑瞪羚，均来自KKWRC）进行全基因组测序（WGS），采用Illumina NovaSeq X平台PE150双端测序，平均深度约30×；以达马瞪羚（Nanger dama）参考基因组（GenBank: GCA_917880005.1）进行读段比对（BWA-MEM）和变异检测（GATK最佳实践流程）；以Cuvier's瞪羚（Gazella cuvieri）全基因组数据（ENA项目PRJEB41985）作为外群进行系统发育分析，包括邻接法（NJ）和最大似然法（ML，IQ-TREE v2.2）；主成分分析（PCA）用于可视化群体结构；基于纯合片段（ROH）分析估计近交系数（F_ROH）；计算全基因组杂合度；基于次要等位基因计数（MAC）进行等位基因频率谱分析；采用Weir-Cockerham方法计算成对F_ST。

四、研究结果

4.1 谱系基因组关系

基于全基因组SNP的系统发育分析（NJ和ML两种方法）均显示，三个谱系的个体按其管理群组分别聚类，深灰瞪羚形成支持度极高的单系分支，大陆瞪羚和法拉桑瞪羚聚类关系更近，G. cuvieri作为外群提供了明确的系统发育根。这一结果证实了三个谱系之间存在显著的基因组分化。

4.2 主成分分析

PCA结果显示，PC1解释了42.39%的遗传方差，将深灰瞪羚与大陆和法拉桑瞪羚完全分离；PC2解释了23.46%的方差，将大陆瞪羚与法拉桑瞪羚进一步区分。三个谱系在PCA空间中形成三个独立簇，与其管理分组完全对应，表明深灰瞪羚与另外两个谱系之间存在显著遗传分化，大陆与法拉桑谱系之间也存在分化。

4.3 纯合片段与近交分析

在100 kb阈值下，深灰瞪羚的平均近交系数（F_ROH = 0.344）远高于大陆瞪羚（0.104）和法拉桑瞪羚（0.045）。深灰瞪羚的ROH段数（平均5,698.67）和ROH总长度（平均894 Mb）也远超大陆瞪羚（1,770.5段、270 Mb）和法拉桑瞪羚（750.5段、118 Mb）。由于参考基因组仅为contig水平，所有检测到的ROH片段均短于0.5 Mb。不同阈值（300-1000 kb）下的分析结果一致，证实深灰瞪羚基因组中分布着高密度的纯合片段。

4.4 全基因组杂合度比较

深灰瞪羚的基因组杂合度极低，个体估计值集中在0.04-0.045之间；大陆瞪羚为中等水平（0.12）；法拉桑瞪羚最高（0.14-0.15）。这一结果与ROH分析一致，进一步证实深灰瞪羚经历了严重的遗传多样性丧失。

4.5 基于次要等位基因计数的等位基因频率谱

深灰瞪羚表现出固定位点的强烈富集（74.6%，16,311,208个SNP），单拷贝变异比例极低（8.5%）；法拉桑瞪羚则相反，单拷贝变异比例最高（57.7%，12,603,974个SNP），固定位点比例最低（32.0%）；大陆瞪羚居中（固定位点42.3%，单拷贝变异43.6%）。双拷贝变异仅在深灰瞪羚中检出（6.2%）。这一模式与深灰瞪羚经历严重遗传瓶颈和近交衰退的结论一致。

4.6 成对F_ST分化

深灰瞪羚与大陆瞪羚之间的加权F_ST最高（0.586），深灰瞪羚与法拉桑瞪羚之间次之（0.513），大陆瞪羚与法拉桑瞪羚之间最低（0.157）。由于各群体样本量仅2-3只个体，这些F_ST值仅作为描述性分化指标，而非群体分歧的确定性估计。

五、讨论与结论

讨论部分指出，深灰瞪羚表现出的严重遗传贫乏与其封闭的繁殖历史密切相关。该群体于1992年仅由5只奠基个体（2雄3雌）建立，此后从未引入任何外部个体或基因流，所有原始奠基者均已死亡，当前个体均为封闭群体的后代。奠基者代表性不均、繁殖偏斜、遗传漂变和持续封闭繁殖等因素共同导致了有效种群规模缩小和遗传多样性加速丧失。相比之下，法拉桑瞪羚虽然长期隔离于岛屿环境，但其繁殖计划自2004年起持续从野生法拉桑种群引入个体，且遗传证据表明其与大陆谱系的分化仅发生在约2000年内，加之法拉桑群岛种群历史上避免了大规模狩猎和栖息地破碎化，因而保持了较高的遗传多样性。大陆瞪羚的中间状态则与其更长的圈养历史（自1988年起）和野外种群的小型化、碎片化有关。研究结论强调，鉴于样本量小、圈养来源、异源参考基因组以及缺乏种群历史分析，观察到的基因组差异应仅限于所研究的个体，不能作为物种边界或更广泛进化关系的证据。本研究为KKWRC的个体水平遗传管理提供了初步基础，遗传拯救应作为保护选项之一，但需通过扩大采样、系谱分析和遗传与表型兼容性评估加以验证。