在物种大灭绝的背景下，仍有少数幸运儿成功地从灭绝边缘“反弹”了回来。然而，对于这些物种而言，除了外部保护措施的加强，其自身内在的遗传和进化机制在其中扮演了何种角色，人们知之甚少。海南长臂猿(Nomascus hainanus)就是这样一个神秘反弹的案例：作为全世界最稀有的灵长类动物，其种群数量从上世纪50年代的大约2000只锐减至2003年左右的仅存13只，而后却奇迹般地稳步增长，目前已恢复到约42只。这种戏剧性的变化为我们探究濒危物种恢复的内在驱动力提供了一个绝佳的研究窗口。一个核心的科学问题是：在经历了如此剧烈的种群瓶颈后，是什么内在的遗传机制帮助海南长臂猿成功“触底反弹”？解答这个问题，不仅关乎这一物种的未来，也为其他濒危物种的保护策略提供了关键见解。
为了回答上述问题，研究人员在《SCIENCE ADVANCES》上发表了一项开创性的研究。他们面对从野外极度濒危、难以接近的灵长类动物身上获取高质量DNA样本（如血液、组织）的巨大挑战，转而开发了一套系统性的非侵入性研究方法。首先，他们从海南长臂猿当前分布区内采集了42份新鲜粪便样本，通过微卫星系统识别出5个现存家庭群中的18个个体（约占整个种群的40%）。其次，利用基于胞嘧啶-磷酸-鸟嘌呤(CpG)甲基化的富集方法，显著提升了粪便样本中宿主DNA的比例，建立了可靠的下一代测序(NGS)基因组数据生成、质量评估和基因分型错误估计（如等位基因丢失率ADO）流程。这项研究创造性地结合了当前的粪便基因组数据和历史博物馆样本（4个）的基因组数据，获得了跨越长臂猿种群两个不同人口统计阶段（衰退期和恢复期）的群体基因组数据。研究人员运用这一研究系统，重点探讨了两个核心问题：(i) 长期的人口统计历史在海南长臂猿近期种群反弹中扮演了何种角色？(ii) 自然基因流是否以及如何促进了这一反弹？
关键研究方法包括：
- 1.
样本收集与个体鉴定：通过野外跟踪和声音定位，系统采集海南长臂猿粪便样本，并利用微卫星标记进行个体识别和性别分子鉴定。
- 2.
粪便基因组学技术：建立包括宿主DNA富集、文库构建、高通量测序及严格数据质控（评估基因分型错误、GC含量偏差等）在内的完整分析流程，从低比例宿主DNA的粪便样本中生成可靠的群体基因组单核苷酸多态性(SNP)数据集。
- 3.
比较基因组学分析：将海南长臂猿的遗传参数（如有效种群大小Ne、杂合度、近交水平、遗传负荷）与从公共数据库获取的其他23种灵长类动物（18种受威胁，5种未受威胁）的数据进行比较。
- 4.
人口历史重建：使用多种方法（如PSMC、Stairway Plot、PopSizeABC、GONE）重建海南长臂猿从百万年前至今，特别是近期的有效种群大小变化历史。
- 5.
进化模拟与功能分析：通过正向时间模拟(SLiM)评估历史人口波动对近交和遗传负荷积累的影响；利用GC4（四倍简并位点的GC含量）作为局部重组率的代理，分析重组与遗传负荷及平衡选择的关系；利用群体结构分析(ADMIXTURE)、系统发育分析和亲缘关系重建(RELPAIR)揭示种群内的隐秘谱系及其杂交事件。
- 6.
种群生存力分析(PVA)：结合基因组结果和物种生活史特征，使用SLiM进行基于个体的随机模拟，预测不同交配模式（谱系间杂交 vs. 随机交配）对未来100代种群灭绝率、增长率和近交水平的影响。
研究结果部分：
1. 海南长臂猿粪便基因组数据的生成
研究建立了一套系统的粪便基因组学流程。从42份粪便样本中成功识别出18个个体，并利用CpG甲基化富集法有效提高了宿主DNA比例。通过评估等位基因丢失(ADO)率，确定在10×测序深度下，基因分型错误率极低，数据可靠。与博物馆样本的GC含量对比显示，富集方法未引入显著GC偏好性，测序读数在基因组窗口上呈正态分布，表明测序片段具有代表性。
2. 海南长臂猿当前的遗传状况
基于SNP数据重建了覆盖所有5个家庭群、跨越三代的较完整谱系。分析显示，在种群衰退期（1950-2003年），海南长臂猿的有效种群大小(Ne)估计为2，是所有已知拥有群体基因组信息的灵长类物种中最小的。其基因组水平的杂合度也较低。然而，出乎意料的是，海南长臂猿的近交水平(FE)和遗传负荷（包括功能丧失LoF突变和有害突变）在受威胁灵长类中属于最低之列，甚至与大多数未受威胁的灵长类相当。
3. 千年尺度的种群扩张缓解了长期末次盛冰期(LGM)相关瓶颈的负面遗传效应
人口历史重建表明，海南长臂猿的Ne从约100万年前开始下降，在末次盛冰期(LGM)前后达到最低点。随后，种群开始扩张，并在约5000年前达到峰值，之后再次下降至今。不同于大多数此前报道的持续衰退的濒危物种，海南长臂猿在约5千年前经历了一次独特的近期扩张。正向时间模拟显示，伴随着LGM后Ne的增加，海南长臂猿的近交水平和遗传负荷显著降低，并在5千年前达到最低值。尽管此后随Ne再次下降，遗传负荷略有回升，但其当前水平仍远低于假设的、没有这次扩张的经典岛屿种群历史模型下的预期。这表明，约5千年前的这次近期Ne大幅增长，帮助海南长臂猿缓解了长期LGM相关瓶颈的负面遗传后果。
4. 高局部重组有助于降低遗传负荷并维持功能变异
使用GC4作为局部重组率的代理，研究发现海南长臂猿在高重组区域的遗传负荷显著低于其他受威胁灵长类，甚至低于未受威胁的灵长类。通过计算ß统计量，在全基因组范围内鉴定出503个受平衡选择的位点，其中28%位于基因区（22个基因）。这些受平衡选择基因内部的重组率显著高于其侧翼序列或其余基因。这表明高重组可能通过提高整体选择效率，帮助这个小型种群维持了适应性相关基因的高度变异。在鉴定出的22个基因中，包括与生殖相关的TEKT3和CCT6B基因，以及与碳水化合物代谢相关的βGLU基因，后者可能帮助长臂猿适应不同季节的饮食变化。
5. 两个古老谱系的自然杂交促进了海南长臂猿种群的反弹
尽管生态调查认为海南长臂猿仅存一个种群，但核基因组SNP分析意外地发现了两个隐秘但不同的谱系（谱系A和谱系B）。这两个谱系对应着2003年仅存的两个最古老家庭群（A群和B群）。CoalHMM模拟估计这两个谱系大约在450年前（~157至576年前）发生分化。谱系结构分析、主成分分析(PCA)和系统发育分析均支持这一划分。研究发现，2011年形成的新家庭群C的奠基者确实来自这两个不同的遗传谱系，且模拟证实从谱系分离到C群形成期间，两个谱系间没有基因流。因此，C群是由谱系A和B的迁出个体杂交建立的。对杂交后代C21的分析显示，其杂合度显著高于其父母，并与1960年（即最近种群衰退开始前）采集的博物馆样本杂合度相当。这种谱系间的杂交为种群注入了更高的遗传多样性。此外，在2015年和2019年新形成的D群和E群中也观察到了谱系间的杂交事件。种群生存力分析(PVA)模拟预测，在未来100代中，如果所有或一半的交配发生在不同谱系个体间，种群的灭绝率（分别为5%和10%）将远低于不考虑谱系杂交的随机交配情景（40%），且种群增长潜力更大，近交水平显著更低。
结论与讨论：
本研究通过整合粪便与历史基因组数据，揭示了海南长臂猿这一极危灵长类物种在神秘反弹过程中多种进化力量之间的相互作用。研究结论指出，该物种的种群恢复得益于其遗传遗产，主要包括一个千年尺度的种群扩张缓解了长期LGM相关瓶颈的负面遗传后果，以及基因组局部高重组率降低了遗传负荷并帮助维持了功能变异。此外，近期两个长期分离的遗传谱系间的自然基因流和杂交促进了这一最濒危灵长类物种的种群增长。
这项研究的意义在于：首先，它系统性地建立并验证了粪便基因组学在保护遗传学中的应用流程，克服了传统方法的限制，为难以获取样本的濒危物种研究开辟了新途径。其次，研究挑战了传统观念，即极度濒危的小种群必然伴随高近交和高遗传负荷。海南长臂猿的案例表明，特定的人口历史（如末次冰盛期后的扩张）和基因组特性（如高重组率）可以共同作用，缓解小种群的遗传困境。再者，研究首次在野生濒危灵长类中证实，长期分离的亚群或谱系间的自然杂交（而非人为遗传拯救）可以成为种群恢复的重要驱动力，提高了种群遗传多样性并降低了灭绝风险。最后，该研究强调了在保护生物学中，除了关注外部威胁和人为保护措施，深入理解物种自身的人口统计历史、基因组架构和行为调节机制同样至关重要。这些发现为评估其他濒危物种的恢复潜力及制定更有针对性的保护策略（如促进适宜栖息地间的个体迁移以鼓励自然基因流）提供了重要的理论依据和实践视角。
