扬子鳄人工圈养繁殖降低近交水平但增加遗传负荷

《Integrative Zoology》：Captive Breeding of the Chinese Alligator Reduces Inbreeding but Increases Genetic Load

　　

近交衰退（inbreeding depression）是小种群长期生存能力的主要威胁。尽管圈养繁殖（captive breeding）与遗传拯救（genetic rescue）是保护生物学中的主要替代方案，但其基因组层面的影响在很大程度上仍属未知，且这些手段可能通过引入新的有害突变而加剧衰退。研究人员基于247只极度濒危扬子鳄（Chinese alligator）的全基因组测序数据，构建了高分辨率遗传变异图谱，鉴定出三个不同的遗传种群，厘清了其种群历史，并进一步支持了美国种群具有独立起源。圈养与半野生种群间的比较表明，圈养管理在降低近交水平与提高遗传多样性方面优于半野生管理，但与此同时增加了有害突变的累积，这可能源于净化选择（purifying selection）的放松。尽管该物种的近交系数显著高于其他濒危物种，遗传净化（purging）可能在减少隐性有害突变数量方面发挥了重要作用。前向时间模拟（forward-in-time simulations）表明，种群混合可降低种群衰退风险，然而当前种群对未来环境变化仍较为脆弱。这些发现为制定扬子鳄平衡且有效的保护策略提供了宝贵的基因组学证据。

扬子鳄保护基因组学研究：圈养管理对近交、遗传负荷与种群生存力的综合影响

一、研究背景与科学问题

扬子鳄（Alligator sinensis）作为中国特有的珍稀爬行动物，属于自中生代延续至今的古老爬行类谱系，历史上曾广泛分布于长江中下游地区的湿地、沼泽和池塘。然而，受气候变化和人类活动的影响，其自然栖息地已严重破碎化，目前仅残存于安徽、江苏和浙江三省交界区域。为阻止该物种灭绝，20世纪80年代起，安徽和浙江两省分别以212只和11只奠基个体启动了圈养繁殖项目，并取得了显著成效——目前安徽和浙江圈养种群规模已分别达到约15000只和6600只，且已有超过2500只圈养个体被放归野外。此外，美国也建立了一个约一个世纪前从中国引入的扬子鳄种群，但其奠基者来源尚不明确。

尽管圈养种群规模庞大且持续增长，但有限的奠基个体不可避免地导致了进一步近交，对扬子鳄的长期生存构成潜在威胁。安徽和浙江圈养种群已表现出明显的近交衰退迹象，包括繁殖成功率下降和后代畸形率增加。当前中国圈养扬子鳄种群由核心圈养繁殖种群（cbAH）和半野生自然繁殖种群（swAH、swZJ、swAM）组成，两类种群均依据系谱或遗传数据进行亲缘评估以尽量降低近交水平。然而，圈养管理在保护成效方面的遗传学评估仍存在大量空白，这限制了圈养管理措施的进一步改进。此外，爬行动物对气候变化尤为敏感，全球气温上升将对扬子鳄施加新的选择压力，而其在遗传上可能难以快速适应。因此，亟需全面了解扬子鳄各群体的遗传特征，为制定科学的保护策略提供依据。

二、研究设计与主要技术方法

研究人员收集了来自半野生种群（浙江swZJ、安徽swAH、美国swAM）和核心圈养繁殖种群（安徽cbAH）共247只扬子鳄个体的全基因组测序数据，其中包含52只高深度（>10×）个体。数据来源于国家基因组科学数据中心（NGDC）和中国国家基因库数据库（CNGBdb）。研究采用了多项关键技术方法：利用GATK流程进行变异检测与质控；通过主成分分析（PCA）、ADMIXTURE和最大似然系统发育树解析种群遗传结构；基于线粒体基因组组装进行单倍型网络分析；运用fastsimcoal2推断种群历史动态；通过D统计量、qpGraph和TreeMix分析基因流；计算遗传多样性指标（π、杂合度）、近交系数（F_ROH）和突变负荷（LOF、dnsSNP、GERP评分）；使用梯度森林（GF）模型预测基因组偏移（genomic offset）；并运用SLiM进行前向时间遗传模拟。

三、研究结果

3.1 大规模种群基因组分析确认三个独立种群

基于247个个体的全基因组数据，研究人员鉴定出101,386,409个变异位点，经严格质控后保留675,402个SNP用于遗传结构分析。PCA、系统发育树和ADMIXTURE分析一致表明，扬子鳄目前分为三个明显不同的种群，分别对应安徽、浙江和美国种群，其中cbAH与swAH在PCA中良好混合。

3.2 美国种群的独立起源

模拟研究表明，若美国种群起源于安徽种群并通过遗传漂变分化，至少需要30–60代（约600–1200年）才能达到当前观察到的F_ST值（0.14），这远超美国种群离开中国后经历的不到十代的繁殖时间。因此，美国种群极可能源自一个与现存浙江或安徽种群不同的独立祖先谱系，该谱系的野生种群可能已在野外灭绝。

3.3 种群历史动态

线粒体基因组分析仅鉴定出4个单倍型，其中1个为三个种群共享，表明圈养扬子鳄种群经历了严重的种群瓶颈。滑窗系统发育分析显示近80%的窗口支持安徽种群为祖先类群，浙江与美国种群形成姐妹支。fastsimcoal2最佳模型与滑窗树一致，TreeMix和D统计量分析均检测到美国种群与安徽种群之间存在强烈的基因流，提示美国种群的祖先来源可能与安徽种群在地理上相邻或历史上存在联系。

3.4 保护管理提高遗传多样性并缓解近交

核心圈养种群（cbAH）的平均杂合度和π显著高于三个半野生种群，整体近交水平（F_ROH>100kb）显著更低。短ROH片段（<2 Mb）在cbAH中显著高于半野生种群，表明约25代前的历史近交在安徽野生种群中更为严重；而长ROH片段（>4 Mb）在cbAH中显著更低，体现了近10代内良好繁殖管理对缓解近交的积极作用。然而，即便管理最优的cbAH种群，其最低F_ROH仍超过50%，远高于许多其他濒危物种。

3.5 不同种群的突变负荷与遗传净化

cbAH种群中检测到显著更多的有害突变（LOF和dnsSNP）和更高的突变负荷，但纯合有害突变比例更低。所有扬子鳄种群均表现出高效的遗传净化效应——ROH区域内外的LOF频率相似，且F_ROH与有害突变数量呈负相关。cbAH种群的净化效应相对较弱，与其较低近交水平和较高杂合有害突变保留一致。

3.6 模拟揭示遗传负荷动态

在推断的种群历史下进行全基因组编码区模拟，结果显示随着种群持续衰退，实际遗传负荷（realized genetic load）增加，主要由隐性有害突变纯合度上升驱动。但以总遗传负荷评估时，持续衰退的扬子鳄种群表现出显著的遗传净化效应，主要源于净化选择对最有害和隐性突变的清除。引入美国个体至安徽或浙江种群的模拟显示，两种情景下杂合有害突变均减少、纯合有害突变增加，且浙江种群中的效应更为明显。

3.7 基因组偏移预测种群适应风险

梯度森林模型预测，在高排放情景（SSP5-8.5）下，扬子鳄当前分布区的大部分区域呈现较高的基因组偏移，表明现有种群对未来气候变化较为脆弱。东部浙江和上海地区预计将提供更适宜的气候条件。反向偏移主要出现在西部栖息地，而正向偏移在东部栖息地更为显著。

四、讨论与结论

美国扬子鳄种群保留了现存中国种群所不具备的独特遗传多样性，具有重要的保护价值。前向模拟表明，无论将美国个体引入安徽还是浙江种群，后代种群的有害突变数量均会减少，为未来将美国种群引回中国提供了科学支持。圈养管理有效降低了近交、提高了遗传多样性，但同时导致有害突变在杂合状态下大量保留，这可能源于圈养环境中净化选择的放松。尽管遗传净化已成功消除了大量本应累积的遗传负荷，但扬子鳄种群的遗传威胁仍然严峻——管理最优的种群基因组中仍有超过50%的ROH片段，远高于多数其他濒危物种。辅助基因流（assisted gene flow）被广泛视为遗传拯救的有效策略，但需谨慎规划以避免不当引入导致种群崩溃。模拟结果支持将美国种群个体引入浙江种群可更有效地通过打断长纯合区段来降低突变负荷。自2007年起已有3只美国来源的扬子鳄被引回中国并与崇明岛浙江种群混合，该种群已扩展至约50只，但其基因组层面的混合后果仍有待监测。此外，扬子鳄种群在当前栖息地面临显著的基因组偏移，放归项目应优先选择高适合度个体并仔细甄别适宜栖息地。研究也存在一定局限性，包括半野生与核心圈养种群的采样时间略有差异、美国种群样本量较小（仅6只）以及种群间样本量不平衡可能带来的随机抽样效应，未来需扩大采样以进一步验证主要结论。