编辑推荐：
本文通过全基因组测序技术，系统研究了澳大利亚袋鼠岛考拉种群的遗传多样性、近交程度和遗传负荷，揭示了隔离种群在保育方舟（conservation arks）中面临的基因组侵蚀风险。研究发现该种群存在显著的有效种群规模（Ne）下降、长纯合片段（ROH）增加和有害突变负荷升高等问题，强调需要通过遗传拯救（genetic rescue）和长期基因组监测来维持种群长期存续力。
隔离栖息地（如保育公园、围栏保护区和岛屿）已成为许多受威胁物种存续的关键场所。这些"保育方舟"虽能规避当代威胁，但也可能导致遗传多样性丧失和近交衰退，进而加剧灭绝风险。澳大利亚袋鼠岛的考拉种群是典型案例——该种群源于1920年代少数移植个体，现已发展成规模庞大且传染病低发的种群。研究通过高覆盖率全基因组数据，比较该种群与澳大利亚大陆种群的基因组多样性、近交程度和适应性潜力。
岛屿种群不仅特有性较高，还为受入侵物种和土地利用变化威胁的物种提供避难所。然而隔离保护也带来风险：岛屿种群更易受 demographic stochasticity 和环境随机性影响，而基因组侵蚀和近交会进一步增加脆弱性。当自然基因流受限时，封闭种群必须依赖现有遗传变异应对环境变化。但受威胁物种本身多样性较低，若奠基种群规模过小，会引发奠基者效应（founder effect），导致多样性降低和等位基因频率失真。
小种群还面临近交风险。亲缘个体交配会导致纯合度升高，当有害突变以纯合形式存在时，便出现近交衰退（inbreeding depression），降低种群适合度和存活率。虽然近交有时可能通过纯化选择（purging）消除有害等位基因，但遗传漂变（genetic drift）在小种群中可能压倒选择，使轻度有害突变固定（drift load）。这种固定突变只能通过基因流消除，凸显遗传拯救的价值。
考拉作为澳大利亚标志性物种，其保护历史复杂。南部种群因栖息地丧失和皮毛贸易在1900年代初濒临灭绝，南澳大利亚州于1930年代已无野生考拉。早期保护措施包括从维多利亚州向岛屿移植个体，袋鼠岛种群即源于1923-1925年间移植的18-24只个体。与大陆种群不同，袋鼠岛考拉维持着低传染病流行率（如未发现衣原体感染），但疑似存在遗传性疾病（如睾丸发育不全、骨骼异常和肾病）。先前微卫星和简化基因组测序研究均显示该种群遗传多样性降低。
本研究旨在通过全基因组测序：（1）基于基因组数据建立袋鼠岛种群建立过程中的奠基者效应时间和严重程度；（2）通过比较大陆种群，评估瓶颈效应对遗传多样性、近交和适应性的影响；（3）探讨孤立种群管理意义，包括遗传拯救必要性。假设若经历严重奠基者效应和长期隔离，该种群应表现出多样性降低、近交增加（表现为ROH延长）和纯合遗传负荷升高。
研究对103只考拉进行全基因组测序，包括袋鼠岛（n=74）、维多利亚州（n=17）和昆士兰州（n=16）种群。袋鼠岛样本来自两个地点（Parndana和Newland），大陆种群代表南北谱系。采样采用旗套法捕获，麻醉后取皮肤组织活检并进行临床健康评估。全基因组测序平均深度＞30×，使用MagAttract试剂盒提取DNA，Illumina NovaSeq平台测序，比对至考拉参考基因组（phaCin_unsw_v4.1）。另采用DArTseq技术（简化基因组测序）分析袋鼠岛个体亲缘关系。
变异检测获得33,538,007个多态位点，质控后保留16,370,050个常染色体双等位基因SNP。过滤标准包括：次要等位基因计数＞2、映射质量＞50、QUAL/DP＞0.05、平均深度＞8。移除性别连锁区域和长度＜50 kb的contig。连锁不平衡筛选后获得3,483,166个SNP用于分析。
使用VCFTOOLS计算观察杂合度（HO）和期望杂合度（HE），近交系数FIS=1-HO/HE。通过PLINK计算群体间FST，并进行主成分分析（PCA）。
采用GONE软件基于连锁不平衡估计有效种群大小（Ne）变化，重组率设为0.3 cM/Mb，最大c=0.05，运行2000代。假定世代时间6年，将代际时间转换为年份。通过AlphaSimR模拟验证软件性能。
使用PLINK检测＞1 Mb的ROH，参数为：滑动窗口20个SNP，每窗口最多2个杂合SNP，最小片段长度1000 kb。计算FROH（ROH覆盖基因组比例）和近交发生时间（g=100/2rL）。
使用SnpEff注释SNP功能效应，以毛鼻袋熊基因组为外群极化 ancestral/derived 等位基因。计算三类遗传负荷：总负荷、纯合负荷（实现负荷）和杂合负荷（潜在负荷）。通过衍生同义突变标准化计数，采用Rxy比率比较群体间衍生等位基因频率。
基于DArTseq数据，使用KING方法计算 pairwise kinship 系数，识别一度亲缘（kinship＞0.177）等关系。
袋鼠岛种群仅110万SNP多态，显著低于维多利亚（160万）和昆士兰种群（270万）。杂合度指标一致显示袋鼠岛多样性最低（HO=0.0009，HE=0.0010），维多利亚次之（HO=0.0012，HE=0.0013），昆士兰最高（HO=0.0017，HE=0.0018）。PCA显示昆士兰与南方种群明显分化（PC1=35.4%），南方种群内部袋鼠岛与维多利亚形成独立聚类。
GONE分析显示所有种群在约45代前（约275年前）Ne开始下降，与殖民时期吻合。袋鼠岛Ne在22-26代前（约1869-1893年）降至最低（Ne=21），但从13代前（1920年代移植时期）快速恢复。昆士兰种群历史上Ne最大（约7000），但下降幅度也最大（＜100）。模拟测试显示GONE可能高估瓶颈时间约50代。
袋鼠岛个体携带557-695个ROH，总长度1202±84 Mb，FROH=0.387，显著高于维多利亚（625±168 Mb，FROH=0.201）和昆士兰（336±165 Mb，FROH=0.108）。长ROH（≥7 Mb）在袋鼠岛最多，对应近交发生在1877-1995年间。亲缘关系分析显示袋鼠岛整体 kinship 为负值（-0.015±0.051），仅发现8对一度亲缘个体，表明当前近交主要源于历史瓶颈而非当代亲缘交配。
袋鼠岛纯合负荷最高（0.2206 vs 维多利亚0.2017），杂合负荷最低（0.0507 vs 维多利亚0.0673），总负荷无显著差异（0.2713 vs 0.2690）。标准化后，袋鼠岛和维多利亚的错义负荷（missense load）均高于昆士兰（0.4926/0.4905 vs 0.3759），LoF负荷也更高（0.0215/0.0213 vs 0.0195）。Rxy分析显示南方种群衍生有害等位基因频率显著高于昆士兰（LoF:1.26-1.33，错义:1.29-1.89）。
对135只袋鼠岛考拉的临床评估发现：雄性体重10.36±2.52 kg，雌性7.48±1.81 kg。多数个体健康状况良好，但记录到7例睾丸发育不全、1例育儿袋发育不全、4例骨骼异常（包括颌畸形和足部发育不全），提示遗传性疾病负担。
袋鼠岛种群的低多样性、高分化和高ROH表明奠基者效应和遗传漂变持续影响。与维多利亚种群的差异可能低估实际遗传侵蚀，因维多利亚种群本身也可能经历多样性丧失。长ROH表明近亲交配历史悠久，短ROH（1-2 Mb）反映历史Ne较小。昆士兰种群短ROH较少但长ROH较多，提示近期衰退。GONE模拟表明殖民时期前南方种群Ne已较低，1920年代移植后袋鼠岛种群快速扩张。亲缘关系分析支持高近交源于历史瓶颈而非当代行为。
遗传负荷分析显示：袋鼠岛纯合负荷升高，杂合负荷降低，表明潜在负荷转化为实现负荷。总负荷未降低提示纯化选择不足，遗传漂变主导进化过程。错义负荷远高于LoF负荷，说明选择对强有害突变更有效。南方种群衍生有害等位基因频率更高，证明漂变压倒选择。这种模式与历史上小Ne种群已清除强有害突变的理论一致。尽管种群数量大，但遗传负荷可能降低应对气候变化和疾病入侵的韧性。临床发现的遗传性疾病与基因组结果吻合，证实遗传侵蚀的表型代价。
全基因组数据能揭示近交、负荷和群体历史特征，指导 translocation 和遗传拯救策略。类似研究显示近交衰退与 fitness 下降相关（如印度虎、虎鲸），而某些物种（如鼠海豚、伊比利亚猞猁）显示瓶颈后负荷清除。袋鼠岛考拉需遗传拯救提升多样性，但需谨慎选择源种群。模拟软件（如SLiM）可预测遗传救援效果。尽管引入有害突变风险存在，但多样性提升效益通常更大。建议长期监测异质性、近交系数和ROH变化，结合健康评估制定动态管理策略。
保育方舟虽能提供短期避难所，但可能成为遗传陷阱。袋鼠岛考拉种群规模掩盖了基因组风险，需通过遗传拯救和基因组监测维持长期存续力。该案例警示孤立种群管理需平衡 demographic 与 genetic 因素，前瞻性干预避免遗传多样性丧失。
生物通 版权所有
