《Molecular Ecology Resources》：Design and Application of a Genome-Wide SNP Array to Improve Conservation Outcomes in the Critically Endangered Southern Corroboree Frog
编辑推荐：
本文介绍了一项开创性研究：为极危物种南部科罗波里蛙(Pseudophryne corroboree)成功开发了首个两栖类定制高密度SNP（单核苷酸多态性）基因芯片。该芯片包含48,386个位点，成功应用于910只圈养个体的基因分型，揭示了种群的历史瓶颈、近交历史及高连锁不平衡(LD)水平。研究证明该芯片具有跨物种（3个近缘物种）应用潜力，并为通过基因组选择(GS)等手段探究壶菌病(Bd)抗性的遗传基础、制定精准繁殖策略以改善重引入效果，提供了关键工具，标志着两栖类保护基因组学迈入新阶段。
1. 引言
两栖动物是全球受威胁最严重的脊椎动物类群，其中两栖壶菌(Batrachochytrium dendrobatidis, Bd)是导致种群崩溃的主要威胁之一。极危的南部科罗波里蛙（Pseudophryne corroboree）因对Bd高度敏感，野生种群近乎完全崩溃，现主要依靠圈养保种群体维持。为应对此类威胁，保护管理策略正从传统的种群监测转向更主动的、基于基因组的干预。单核苷酸多态性(SNP)芯片作为一种精确、可重复的高通量基因分型工具，在农业和水产养殖中已成熟用于性状遗传解析和选择育种，但在非模型野生动物，尤其是两栖类中的应用仍属空白。开发物种特异的基因组工具，对于理解濒危物种的遗传多样性、适应力以及疾病抗性等性状的遗传基础至关重要。此外，SNP芯片具备跨近缘物种应用的潜力，能为资源有限的类群提供经济高效的基因组学研究方案。对科罗波里蛙而言，其密集管理的圈养种群为研究疾病易感性的遗传基础、并最终实施旨在提高保护成效的选择性育种计划提供了独特机会。
2. 材料与方法
2.1 全基因组测序与SNP挖掘
研究利用23只科罗波里蛙（包括20只野生捕获奠基个体）的全基因组测序(WGS)数据进行变异检测，初始获得约5300万个潜在SNP。经过严格的质量过滤（如移除低质量位点、低复杂度区域、次要等位基因频率(MAF)<10%的位点等），并确保位点在基因组中均匀分布（每Kb保留1个SNP），最终得到537,368个高质量SNP用于芯片设计。
2.2 SNP芯片设计与制备
将过滤后的SNP提交给Thermo Fisher Scientific进行“in silico”评分，从中筛选出设计评分最高的48,386个SNP，成功放置在定制的Axiom myDesign 384HT芯片上，制成了首个两栖类SNP芯片（AXIOM PSCO SNP 384S384）。该芯片平均密度为4.6个SNP/兆碱基(Mb)。
2.3 样本与基因分型
共对1143份样本进行了基因分型，主要包括995只圈养繁殖的F1代科罗波里蛙幼体。此外，还纳入了118份历史项目样本以测试芯片对保存状态不佳样本的分型能力，以及5份颊拭子和皮肤拭子以评估非致命取样效果。为测试跨物种扩增能力，还对三种近缘两栖类（P. pengilleyi, P. bibroni, Crinia signifera）各6个个体进行了分型。基因分型在Ramaciotti基因组学中心完成，使用Axiom Analysis Suite (v5.3.0.45)进行数据分析与质量控制。
2.4 种群遗传学分析
对通过质控的910只F1代个体（来自54个谱系家庭）进行种群遗传学分析。内容包括：主成分分析(PCA)和网络分析(NETVIEW P)评估种群/家庭结构；计算全基因组连锁不平衡(LD)，并以r2衡量；基于LD模式，使用SNP1101、NeEstimator和SNeP等多种软件估算历史及当代有效种群大小(Ne)；利用PLINK检测全基因组范围内的纯合子片段(ROH)，并计算基于ROH的近交系数(FROH)和基于杂合度的近交系数(F)。
3. 结果
3.1 SNP芯片性能
芯片性能优异。在1115份科罗波里蛙组织样本中，1087份（97.5%）通过了≥95%检出率的质控阈值。在成功分型的42,058个位点中，39,701个（占芯片总位点的82.1%）为多态性SNP，其余为单态性或“无次要等位基因纯合子”位点。重复样本的基因分型重复性高达99.6%。
3.2 样本类型、数量与质量的影响
组织类型显著影响分型成功率，趾甲样本的检出率显著高于头部和肝脏组织。DNA浓度对样本能否通过检出率阈值无显著影响，但DNA质量（260/280吸光度比）有显著影响，质量越高的样本通过率反而越低，这可能与测定干扰有关。颊拭子能产生基因型数据，但检出率（80.65%-93.97%）未达到95%的质控阈值，主要受限于极低的DNA浓度(<1 ng/μL)。皮肤拭子则未能成功。
3.3 跨物种扩增
尽管三个近缘物种的样本均未达到标准质控的DQC阈值，但在不过滤的情况下，仍有36,372个位点（占75%）成功跨物种扩增。其中21,077个位点（43.6%）为多态性，展示了该芯片在近缘物种中应用的广泛潜力。
3.4 种群多样性、结构与历史
PCA和网络分析揭示了基于家庭的精细遗传结构，但未发现强烈的亚群分化，这与圈养繁殖计划中采用的跨地理种群交配策略一致。F1代个体的平均近交系数F为-0.018，表明群体平均杂合度略高于预期，符合远交策略的结果。
基于ROH的分析显示，个体基因组中ROH的平均比例(FROH)为0.04（范围0.02-0.08），且绝大多数为短片段（<4 Mb），表明存在历史久远的（而非近期的）瓶颈和近交事件。
连锁不平衡分析显示，F1群体中LD水平非常高，相邻SNP间的平均r2高达0.67。LD衰减缓慢，中度LD (r2≥ 0.2) 延伸至约900 Kb。这种高LD与小的有效种群规模一致。2衡量。">
有效种群大小分析揭示该种群经历了两次历史瓶颈。一次发生在约100代以前（300-500年前），另一次更近的瓶颈发生在约10代以前（30-50年前），Ne从约100急剧下降至当前的29-33。这次近期下降与Bd在1984年传入该物种栖息地的时间点基本吻合。
4. 讨论
本研究成功开发并验证了首个用于两栖类（南部科罗波里蛙）的定制高密度SNP芯片。该芯片表现出高转换率、高重复性和对低质量/低浓度DNA的稳健性，为非致命取样（颊拭子）和历史样本的基因分型提供了可能。跨物种测试证明了其在近缘类群中的实用价值，为其他两栖类保护基因组学研究提供了经济高效的工具。
对圈养F1代群体的基因组分析揭示了该物种面临的严峻遗传状况：高的连锁不平衡、小的有效种群规模以及反映历史瓶颈的ROH模式。这些发现证实了种群长期衰退的遗传印记，而近期的Ne急剧下降很可能与Bd的入侵直接相关。
高密度的全基因组SNP数据为在科罗波里蛙中实施基因组学驱动的保护策略奠定了坚实基础。特别是高的LD水平意味着，该芯片所包含的位点数量（近4万个多态性SNP）已远超进行精确的基因组选择(GS)所需的理论阈值（约1万个）。这使其成为一个强大工具，可用于：（1）持续监测圈养种群的遗传健康，优化谱系管理；（2）探究Bd抗性等适应性状的遗传架构和数量性状位点(QTL)；（3）最终实施针对疾病抗性的基因组选择育种计划，旨在提高圈养个体的适应性，并最终提升野化重引入的成功率。
总之，科罗波里蛙SNP芯片的开发不仅标志着两栖类保护基因组学资源的重要进展，其应用将直接推动针对该极危物种的、基于实证的精准保护管理，并为应对全球两栖类危机提供了可借鉴的基因组学解决方案范式。
