极危状态下的东澳大利亚南古尔棘石龙子（*Nangura spinosa*）的极低遗传多样性与小有效种群大小（*Ne*）

《Conservation Genetics》：Extremely low genetic diversity and small Ne in the Critically Endangered Nangur Spiny Skink (Nangura spinosa) of eastern Australia

　　

在全球范围内，许多物种拥有小而片断化的种群，因此遗传因素成为其持续生存的主要威胁。极危的南古尔棘石龙子（*Nangura spinosa*）栖息于澳大利亚昆士兰州东南部南古尔国家公园与奥克维尤国家公园的藤蔓灌丛中。这两个种群自然隔离，且目前规模极小，南古尔种群的数量估计<50只，奥克维尤种群约为1000只。研究人员利用全基因组单核苷酸多态性（SNP）数据集进行了遗传健康评估，以提供物种恢复的管理行动。这些种群形成两个截然不同的遗传群体，具有高遗传分化。相比之下，每个种群内部存在高连通性，可能由短距离的规律扩散和较长距离的偶发移动驱动。奥克维尤种群具有低杂合度（*H_obs* = 0.139；*H_exp* = 0.142；平均个体*F* = 0.116），并包含72%的SNP。南古尔种群具有极低的杂合度（*H_obs* = 0.114；*H_exp* = 0.112；平均个体*F* = 0.289）。这些种群的*N_e*估计值极小：南古尔约为16–43，奥克维尤约为66–179。奥克维尤种群的遗传健康指标与其他极危澳大利亚石龙子相当，但南古尔种群的情况更糟。南古尔种群的结果表明近交水平高、亲缘关系强且个体杂合度极低。不良的遗传健康可能已经导致适合度下降和种群衰退，尤其是南古尔种群。研究人员建议进行种群内迁地和/或人工管理繁殖以减缓遗传侵蚀。此外，迫切需要通过在种群间进行杂交来实现遗传拯救。这项研究是遗传数据如何辅助物种恢复的明确例证。

**论文解读：东澳大利亚极危南古尔棘石龙子（*Nangura spinosa*）的遗传健康评估与保护管理启示**

**1. 研究背景与问题**

在物种保护中，环境因素（如栖息地丧失、入侵物种、气候变化和疾病）通常被视为主要威胁。然而，对于已受这些因素威胁的物种，遗传过程可能变得同等重要。受威胁物种的种群通常小而孤立，导致遗传多样性快速丧失和近交增加，这被称为遗传健康下降。不良遗传健康可直接降低个体繁殖适合度，削弱种群适应环境变化的韧性，从而增加灭绝风险。若忽视遗传问题，即使应对其他生态和人口威胁（如林火或入侵物种），其效果也可能受限。因此，有效的物种恢复计划必须同时解决遗传与非遗传因素。

南古尔棘石龙子（*Nangura spinosa*）是一种中型石龙子，栖息于澳大利亚东部高度局限区域的半常绿藤蔓灌丛中，仅存两个种群：奥克维尤国家公园（Oakview）和南古尔国家公园（Nangur）。两种群相距约36公里，已因不适宜干燥栖息地隔离至少数十万年。该物种在州和国家层面被列为极危（Critically Endangered），国际层面为濒危（Endangered）。南古尔种群2020-2021年估计仅44只成体，奥克维尤种群约1000只。主要历史威胁是栖息地清除（用于木材采伐、松树种植和放牧），目前仍面临入侵动植物、林火和气候变化等持续威胁。种群规模已接近或低于短期遗传健康推荐的1000只成体（约相当于*N_e* = 100）阈值，遗传侵蚀可能在几代内累积至有害水平，因此亟需开展遗传健康评估。

**2. 主要关键技术方法**

研究人员采用了以下关键技术方法：首先，通过野外采样（2020-2021年从南古尔、奥克维尤及人工圈养种群采集组织样本，另包含2006-2007年历史样本）获取基因组DNA。其次，利用双消化限制性位点相关DNA测序（ddRAD-seq）进行简化基因组测序，生成约100 bp单端读段。随后，使用Stacks v2.62进行单核苷酸多态性（SNP）检测与分型，构建三个数据集：数据集1（每座位单SNP，用于种群结构及有效种群大小[*N_e*]分析）、数据集2和3（含全部位点，包括SNP与不变位点，用于遗传多样性评估）。最后，运用*K*-均值聚类、netviewr网络分析、分子方差分析（AMOVA）评估种群结构与连通性；通过Colony v2.0.7.1推断亲缘关系；利用NeEstimator v2.1基于连锁不平衡（LD）方法估算*N_e*；计算观察杂合度（*H_obs*）、期望杂合度（*H_exp*）、近交系数（*F*）及多位点杂合度（MLH）等遗传健康指标。

**3. 研究结果**

**3.1 全物种结构与连通性**
*K*-均值聚类和netviewr分析均支持将*N. spinosa*分为两个截然不同的遗传群体（南古尔与奥克维尤），两者间*F_ST*高达0.479，等位基因频率差异（AFD）为0.288，平均基因组距离为0.301。种间AMOVA显示约47%的遗传变异由种群间差异解释，且86%的SNP在其中一个种群中具有私有等位基因。这表明两群体间存在显著遗传分化，且已持续长时间隔离。

**3.2 种群内结构与连通性**
南古尔和奥克维尤内部均未检测到显著亚群结构，暗示各自内部近期近乎随机交配。netviewr网络显示，南古尔内存在细微亚网络（与地理采样区对应），但低*k*-NN值下即相互连接。奥克维尤北部和东部采样区形成明显亚网络，但同样在*k*-NN ≥ 4时连通。AMOVA表明种群内亚结构仅解释少量遗传变异（南古尔≤10%，奥克维尤约5%）。个体间基因组距离与地理距离的关系呈对数曲线，在约500米处趋于平稳，表明扩散主要发生在短距离内。遗传邻域半径估计为约1000米，对应亲本-后代扩散距离标准差约500米，提示个体扩散能力有限。

**3.3 种群间个体关系与扩散**
Colony分析在南古尔鉴定出11对一阶亲缘关系（全同胞或亲子对）和27对二阶关系，形成多个家族群；奥克维尤有18对一阶关系和10对二阶关系。一阶关系个体间地理距离通常<40米，部分可达200-300米。所有一阶和二阶关系个体均相距<500米，进一步支持短距离扩散模式。

**3.4 遗传健康指标**
**种群遗传多样性**：两群体遗传多样性均极低，南古尔所有指标均更差。南古尔观察杂合度（*H_obs* = 0.114）比奥克维尤（0.139）低约20%，多态性SNP少43%，私有等位基因少53%。近交系数*F_IS*南古尔（0.054）远低于奥克维尤（-0.002），但需注意*F_IS*在高度近交群体中可能因等位基因频率偏倚而失真。全基因组杂合度（含不变位点）南古尔为0.00046，奥克维尤为0.00056，远低于其他受威胁石龙子。

**个体遗传多样性与近交**：个体多位点杂合度（MLH）在南古尔为0.081–0.141（均值0.112），奥克维尤为0.089–0.157（均值0.137）。南古尔个体近交系数*F*均值0.289，近四分之三个体*F* > 0.25（相当于全同胞子代近交水平），部分个体接近0.5（克隆子代水平）。奥克维尤*F*均值为0.116。2006-2007年与2020-2021年样本比较未发现杂合度明显下降，表明高近交已持续至少十余年。南古尔个体间最大基因组距离仅相当于奥克维尤二阶关系水平，提示“谱系崩溃”。

**有效种群大小**：南古尔LD *N_e*估计均值16–43，奥克维尤均值67–180。南古尔所有*N_e*估计值均低于短期遗传健康阈值（*N_e* = 100），奥克维尤接近该阈值。去除一阶亲缘关系后*N_e*估计值翻倍，但整体仍显示极小有效种群规模，将加速遗传多样性丧失。

**3.5 圈养种群遗传健康**
基于奥克维尤14只野生个体建立的圈养种群保留了71%的多态性位点，杂合度略低于野生群体（*H_obs* = 0.135），但含一只极近交个体（MLH = 0.089）。为作为保险种群保留95%遗传多样性，需至少再引入26只野生个体。

**4. 讨论与结论**

**讨论总结**：研究结果表明，*N. spinosa*的遗传健康极差，表现为低遗传多样性、高近交和小有效种群大小。南古尔种群问题尤为严重，已处于临界点，可能已导致适合度下降和种群衰退；奥克维尤相对较好但仍在恶化。两者间遗传分化显著，但内部近乎随机交配，扩散距离有限。圈养种群需进一步扩大规模并优化配对。基于发现，研究人员提出三条管理建议：第一，认识到不良遗传健康使物种更易受其他威胁影响，需继续实施传统保护措施（如栖息地管理、入侵防控）；第二，在种群内定期迁移个体（约1000米距离）至不同地理簇以降低近交；第三，遗传拯救（种群间杂交）是短期和长期改善遗传健康的必要措施，应通过小规模试验评估杂交后代适合度，权衡近交衰退与远交衰退风险。最终结论为：若不采取干预，南古尔种群几乎必然灭绝，甚至可能导致整个物种灭绝；遗传数据为物种恢复提供了明确指导。

**研究结论翻译**：这项研究明确展示了遗传数据如何辅助物种恢复。基于对南古尔棘石龙子的遗传健康评估，研究人员发现极低的遗传多样性、高近交水平以及极小的有效种群大小，尤其在Nangur种群已达到临界点，可能已导致适合度损失和种群下降。研究人员建议立即采取种群内迁地管理以及种群间遗传拯救措施，以减缓遗传侵蚀并避免灭绝。