土耳其本土绵羊种群遗传背景再探：瓶颈效应与迁移事件的分子证据

这篇综述通过28个高多态性微卫星位点对土耳其四个本土绵羊品种（GKR、KRK、KNG、NRD）进行基因分型，结合Wilcoxon检验（TPM/SMM/IAM模型）和TreeMix算法，揭示了种群未经历近期遗传瓶颈（有效种群规模Ne>采样量），并发现NRD品种因地理隔离呈现遗传独特性（迁移率0.0096至KRK-GKR分支）。研究为保护濒危绵羊遗传资源（如NRD纯种繁育）和可持续育种策略（SNP芯片/NGS技术应用）提供了关键数据支撑。

　　

遗传瓶颈效应分析
通过Bottleneck v1.1.2.02软件对120头土耳其绵羊（GKR/KRK/KNG/NRD各30头）进行突变漂变平衡检测。Wilcoxon符号秩检验在TPM模型下概率值达1.00（GKR/KNG）和0.99（NRD/KRK），配合等位基因频率的L型分布模式（图1），证实四品种均未经历近期种群规模缩减。LDNe软件估算的Ne值显示，GKR群体在最低等位基因频率0.05时Ne达820，而KNG因采样误差出现"无限大"估值，暗示其可能源自更大规模的祖先群体。

种群迁移网络构建
TreeMix算法（OptM包确定最优迁移事件数=3）绘制出三条关键基因流：NRD→KRK-GKR分支（迁移率0.0096）、KNG→GKR（最高迁移率0.0176）。值得注意的是，NRD作为"供体"而非"受体"的迁移模式（图3），与其长期局限于凡省（Van）的地理隔离特性相符。这种单向基因流动解释了早期微卫星研究中NRD的遗传分化现象（F_ST
=0.12）。

可持续育种启示
尽管当前未检测到遗传瓶颈，但NRD的孤立繁育模式可能加剧未来近交衰退风险。建议在保护计划中结合微卫星（如MAF214/DYMS1位点）与高通量技术（Illumina SNP芯片），同时建立迁移走廊以打破NRD的遗传隔离。研究特别强调，KNG品种的"无限大"Ne估值需通过扩大采样量（建议n≥50）验证，避免因连锁不平衡(LD)估算偏差导致误判。

方法论创新
本研究首次将TreeMix算法应用于土耳其绵羊迁移分析，克服了传统F_ST
仅能反映分化程度的局限。通过比较三种突变模型（IAM/TPM/SMM），发现TPM模型下Wilcoxon检验与模式偏移指示器的结果一致性最高，为后续濒危物种瓶颈分析提供了模型选择依据。