突破瓶颈：考拉（Phascolarctos cinereus）通过种群扩张实现遗传恢复的途径

《SCIENCE》：Escaping bottlenecks: The demographic path to genetic recovery in koalas (Phascolarctos cinereus)

为应对种群数量锐减（瓶颈）带来的遗传多样性丧失和进化潜力降低等核心问题，研究人员对澳大利亚多地考拉种群开展了全基因组研究。研究发现，历史上经历瓶颈但正在恢复的种群，其有效种群大小（Ne）正在增加，并通过快速的种群扩张重组了遗传变异，减少了有害突变负荷，重新获得了进化潜力。这一研究为受威胁种群的保护策略提供了新的理论依据。

　　

在保护生物学中，有一个令人担忧的经典理论：当一个物种的种群数量骤降到极低的水平，如同通过一个狭窄的“瓶颈”时，随之而来的将是遗传多样性的急剧丧失、近交风险的增加以及有害突变的积累。这些遗传上的困境，很可能将种群拖入“灭绝漩涡”——数量减少和多样性降低相互强化，最终走向灭亡。然而，自然界似乎存在一些“奇迹”，有些种群在经历了近乎毁灭性的瓶颈事件后，不仅种群数量得以恢复，甚至展现出了遗传上的复苏迹象。这究竟是如何发生的？是哪些机制帮助它们逃离了遗传的“死胡同”？为了回答这个关键的科学问题，一支研究团队将目光投向了澳大利亚的“国宝”——考拉。

考拉曾广泛分布于澳大利亚东部和南部，但因栖息地丧失、疾病和历史上的皮毛贸易而数量锐减，其北方种群已被列为濒危物种。然而，维多利亚州的考拉却上演了一段跌宕起伏的生存故事。在19世纪90年代，为保护它们免于皮毛贸易的捕杀，不到10只考拉被转移至维多利亚州的两个岛屿避难所，构成了一个极端的遗传瓶颈。整个20世纪，这些岛屿后裔被重新引入维多利亚州大陆，种群数量从20世纪20年代仅存的约500-1000只，激增至2020年的近50万只，形成了从濒临局部灭绝到种群过盛的戏剧性逆转。这种“极端瓶颈后快速扩张”的模式，为科学家提供了一个检验种群增长如何塑造遗传结果的绝佳案例。

为了深入探究这一过程，研究人员对来自澳大利亚三个州（昆士兰州、新南威尔士州、维多利亚州）27个地理种群的418只考拉个体进行了全基因组测序，并利用这些数据重建了种群历史和遗传变异模式。

在技术方法层面，本研究主要运用了以下关键技术：1) 种群水平的溯祖模型，用于估计种群历史有效大小（N_e）的长期变化（超过100代前）；2) 基于连锁不平衡（Linkage Disequilibrium, LD）的方法，用于分析更近时期（0-100代前）的种群动态；3) 计算全基因组杂合度、近交系数（F_ROH和F_IS）以及分析纯合性片段（Runs of Homozygosity, ROH），以评估遗传多样性和近交水平；4) 按次要等位基因频率（Minor Allele Frequency, MAF）对等位基因进行分组（常见、低频、稀有），分析其分布模式；5) 使用R_X/Y指标和基因组进化速率谱（Genomic Evolutionary Rate Profiling, GERP）分析，评估不同种群间的功能性变异和突变负荷差异；6) 基于观察到的种群历史进行计算机模拟，验证人口统计学过程是否能单独产生观察到的遗传模式。

研究结果

重建考拉种群的种群历史

对历史有效种群大小（N_e）的重建显示，维多利亚州的种群在约91代前开始经历瓶颈，N_e从1162急剧下降到第41代的102（下降约92.2%），随后在35代内扩张至494。相比之下，昆士兰州和新南威尔士州的种群则在过去10-30代内经历了N_e的急剧下降。这些重建结果与历史管理记录相符，例如，从岛屿避难所重新引入后，维多利亚州大部分种群的N_e迅速增加。Cape Otway种群在1981年由75个个体建立，到2013年已超过1万只，其当代N_e估计值是所有采样种群中最高的。

遗传多样性和近交的模式

分析显示，尽管在常染色体杂合度（aH_O）上，维多利亚州种群仍显著低于北部种群，但差距比之前报告的要小。在近交模式上，维多利亚州表现出最高的累积ROH（F_ROH），平均有68.1%的基因组位于大于50 kb的ROH中，远高于新南威尔士州的41.2%和昆士兰州的36.0%。维多利亚州种群的中等长度ROH（>500且<10,000 kb）显著更高，这表明它们可能经历了多次瓶颈事件。

等位基因的分布

研究人员将杂合度按等位基因频率分组分析发现，在常见等位基因（MAF > 0.05）和低频等位基因（MAF 0.01-0.05）组中，维多利亚州的多样性最低。然而，在稀有等位基因（MAF < 0.01）组中，州级差异减小。分析低频/常见等位基因比率与稀有/低频等位基因比率发现，后者在州级间的差异明显小于前者。这表明，在维多利亚州种群中，稀有等位基因正在快速恢复，这很可能是由种群扩张期间的新突变驱动的。在维多利亚州内部，South Gippsland、Strathbogies和Western Victoria等种群在低频和稀有等位基因箱中表现出更高的杂合度，显示出此前未被认识到的种群内部异质性。

种群模拟

为了验证仅靠人口统计学过程是否能产生观察到的模式，研究人员模拟了维多利亚州种群（50代前发生急性瓶颈）和北部种群（20代前开始渐进性瓶颈）的进化轨迹。模拟结果证实，N_e紧密跟踪种群实际大小，近交在急性瓶颈后急剧增加。杂合度在维多利亚州的严重瓶颈下迅速下降，但在北部渐进性瓶颈下下降更缓慢。等位基因的命运取决于其频率：稀有变异在种群扩张期间恢复最快，而低频等位基因在渐进性衰退中流失更慢。因此，稀有与低频变异之比随扩张而迅速增加，而低频与常见变异之比在瓶颈期间下降，并在恢复后仍保持较低水平。这些动态模拟结果与基因组数据中的等位基因频率模式相符。

功能性遗传多样性

通过对比种群间功能性与基因间变异的R_X/Y指标分析发现，北部考拉显示出更强的净化选择迹象，而维多利亚州种群则携带了更高比例的功能性变异。这与理论一致：快速的种群增长通过放松选择和增加突变输入，使得功能性变异得以积累。在不同等位基因频率组中，人口统计学历史进一步塑造了功能性多样性。新南威尔士州每个样本携带的常见功能性和基因间变异较少，而维多利亚州则拥有最少的稀有和低频变异，这与严重瓶颈相符。维多利亚州丢失的基因间变异比例高于功能性变异，表明基因组不同区域面临着差异化的进化压力。基因组进化速率谱分析也支持这一模式，维多利亚州个体在保守区域携带的有害突变（特别是中到大效应大小的）更少，表明其突变负荷的降低更多是瓶颈期间随机丢失的结果，而非系统性清除。

研究结论与讨论

本研究表明，尽管经历了持久的瓶颈和较低的遗传多样性，维多利亚州的考拉种群显示出恢复的迹象，包括上升的N_e和稀有变异的再现。这一模式与“入侵物种遗传悖论”相似，种群可以通过快速的种群扩张“逃脱”极端瓶颈的遗传限制。重组在这一恢复过程中是关键机制，它通过打乱单倍型，产生新的等位基因组合，类似于新突变，增加了单倍型多样性，并加速了连锁不平衡的衰减，从而提升了N_e。

这些发现对考拉保护具有紧迫的现实意义。长期以来被视为该物种遗传堡垒的北部种群，其N_e正在下降，突变负荷升高，令人担忧。相比之下，维多利亚州种群正在经历种群和遗传上的恢复，积累了经过重组的、可能具有适应性的变异。然而，这种恢复是脆弱的，再次的瓶颈可能抹去新获得的遗传多样性。因此，管理策略必须采用一个前瞻性的框架，考虑稀有等位基因的动态、种群增长和遗传连通性。应考虑通过维多利亚州和北部种群之间的遗传混合来进一步增强物种的恢复力，这可能有助于保护剩余的遗传变异。

总而言之，该研究揭示了种群历史、重组、突变和人为干预以复杂的方式共同塑造了基因组变异的当代模式。经历瓶颈的种群并非一律贫瘠，在扩张条件下，它们可以再生功能性变异，并降低突变负荷。这些发现对保护基因组学有重要意义：首先，低遗传多样性并不总是意味着有限的进化潜力；其次，近交和多样性指标具有背景依赖性，并受到种群历史的塑造；最后，管理策略应纳入近期的种群轨迹以及稀有和功能性等位基因的分布情况。这项研究为受威胁物种的管理提供了新的视角，强调了在保护中考虑动态恢复过程的重要性。相关研究成果发表在《科学》（SCIENCE）杂志上。