蚕豆(Vicia faba L.)群体中已天然存在育种所需的自交系(inbreds)

《Theoretical and Applied Genetics》：Faba bean populations already contain the inbreds needed for breeding

　　

摘要：在蚕豆(Vicia faba L.)育种中，获得自交系(inbred lines)虽至关重要，但因该作物具部分异交(partial allogamy)特性，传统自交系培育过程漫长、昂贵且费力。然而，在任何遗传多样性丰富且开放授粉的蚕豆群体中，此种特定繁殖系统会导致所谓的身份不平衡(identity disequilibrium)，使相当比例的个体成为高度近交(highly inbred)个体。本研究提出一种低成本方法——利用有限数量KASP(Kompetitive Allele Specific PCR)标记从群体中鉴定这些高度近交个体。经模拟研究发现，在标记数较少(10～100个)时估算近交系数(inbreeding coefficient, F)最准确的方法是最大似然估计量FL。在"哥廷根冬蚕豆群体(G?ttingen Winter Bean Population, GWBP)"中，8.4%(728株中的61株)的个体被归类为高近交个体(FL≥ 0.8\overline{3}，即F≥5/6)。研究人员利用模拟数据和经验数据评估了用少量标记分类高近交个体的表现，并提出优化策略。使用50个标记且阈值FL≥0.9时，从群体中鉴定出42个个体，其中36个为真高近交个体(来自已知的61个真阳性)，6个为假阳性。此外，先用少量标记初筛剔除低近交个体，可使基因分型成本降低最高达46%。即便无双单倍体(doubled haploid, DH)技术可用，蚕豆育种者亦可借此方法快速、低成本、低劳力地鉴定高近交个体。

论文解读：蚕豆群体中天然存在可用于育种的高近交个体——基于有限KASP标记鉴定的研究总结

本文发表于《Theoretical and Applied Genetics》，针对蚕豆(Vicia faba L.)因缺乏双单倍体(Doubled Haploid, DH)技术导致自交系(inbred line)选育需多代人工自交、隔离授粉及人工tripping(花瓣刺激自花授粉)，过程漫长费时且成本高昂的问题展开研究。蚕豆具部分异交(partial allogamy, 自交率约0.5)特性，开放授粉群体中因身份不平衡(identity disequilibrium)会自然产生不同近交程度的近交组(inbreeding cohorts)，其中包括相当于F₄及以上(近交系数F理论均值≥0.875)的高度近交(highly inbred)个体，但这些个体无法凭表型识别。本研究旨在开发利用少量低成本KASP(Kompetitive Allele Specific PCR)标记估算个体近交系数(inbreeding coefficient, F)并从中筛选天然高近交个体以替代传统多代自交的新方法，从而加速蚕豆自交系及综合品种(synthetic cultivar)育种进程。

主要关键技术方法：

研究人员以哥廷根冬蚕豆群体(G?ttingen Winter Bean Population, GWBP)——由11个创始自交系混合后经37代开放授粉维持的近平衡部分异交群体——728份个体为经验材料，另设候选春蚕豆综合品种(Syn-4及Syn-8世代)为验证材料。采用Vfaba_v2 Axiom SNP芯片(~60K标记)获取全基因组基因型作基准，模拟三类不同异交率的平衡群体F分布。从芯片数据或模拟数据中抽取5～100个随机或近理想(等位基因频率p≈0.5、连锁不平衡LD低)标记子集，比较六种经典F估计量及最大似然估计量F_L的RMSE。定义高近交阈值为F≥0.8\overline{3}(F₄群体下四分位数)，设置不同近交分类阈值(Inbreeding Classification Threshold, ICT)评估真阳性(TP)与假阳性(FP)，并设计初筛步骤——用极少标记按初步近交分类阈值(Preliminary Inbreeding Classification Threshold, PICT)剔除低近交个体以减少后续基因分型量，计算KASP成本因子评估节约效果。

研究结果：

Inbreeding coefficient estimation（近交系数估算）：

通过模拟部分异交群体(平均自交率S=0.53，最接近GWBP实际情况)发现，随标记数增加各估计量精度上升(5～25标记时提升最显著)；随机等位基因频率下最大似然估计量F_L在≥10标记时RMSE最低、精度最优，Li & Horvitz第一版估计量(F_LH1)次之(需多约10个标记达同等精度)；理想标记(p=0.5且无LD)时所有估计量公式退化为F=1?2H/m(H为观测杂合位点数，m为多态标记总数)。纯合率(Homozygosity rate, HOM)估算精度略低于F_LH1但与常用估计量相当。50个标记时F_L估值与真实F的Pearson相关系数平均达0.94。

Presence and proportion of highly inbred individuals in existing faba bean populations（现有蚕豆群体中高近交个体的存在比例）：

GWBP全芯片数据(F_L估算)：群体平均FF=0.35，推得平均自交率S=0.52；728个体中61株(8.4%)F_L≥0.83，判定为高近交——与模拟部分异交群体预期比例(10.6%)接近。另624份GWBP个体用50个近理想KASP标记估算：34株(5.4%)初判高近交，其中26株(84%)经芯片验证为真高近交。NPZ候选春蚕豆综合品种：Syn-4(F=0.49,S=0.65)高近交率15.4%(70/456)；Syn-8(F=0.37,S=0.54)高近交率13.3%(46/347)，表明开放授粉多代群体均含可检出高近交个体。

Classification of highly inbred individuals（高近交个体分类效能）：

分类效果受群体异交程度、F估算精度及ICT影响。异交率越低(自交率越高)，TP/FP比越高。提高标记数、选用理想标记、提高ICT均可降FP；ICT接近0.83时高精度增TP(减假阴性)，ICT接近1时低精度反而因估值上偏可能增TP。模拟部分异交群体结果与GWBP经验数据吻合良好——如50理想标记ICT=0.90时模拟平均TP=39、FP=4，GWBP实测TP=36、FP=6。零FP需≥80理想标记且ICT=1，但TP极少(≤8)，故应依群体组成权衡标记数与ICT。

The preliminary classification step（初筛步骤）：

以模拟推导的PICT(确保95%概率保留所有真高近交个体)先行初筛，随机标记至少需15个、理想标记至少需10个方可启用。随初筛标记数增加被剔除的低近交个体比例上升(理想标记10个→37%，100个→79%)。成本分析显示两步法(初筛M₁=20标记 + 主分类M₂=100标记)较单步可节省最高46%基因分型费用。初筛后剩余个体进入主分类(更多标记、合适ICT)完成最终高近交判定。

讨论与结论翻译：

研究人员确认蚕豆开放授粉群体中存在显著比例天然高近交个体，可用少量KASP标记结合最大似然估计量F_L识别；通过设置匹配标记集与群体组成的近交分类阈值(ICT)可有效分类；加入初步低标记数初筛步骤在不影响性能前提下进一步降本。育种者可据自身群体预估异交率，参照模拟结果选定标记数与ICT，以低成本从现有群体直接获取自交系原始材料，作为DH技术缺失时的替代方案，亦可与基因组预测联合筛选高产–高近交精英个体用于二阶段育种策略，或偶联重要性状(如vicine-convicine含量)的标记辅助选择。此法理论上亦适用于其他部分异交作物(如甘蓝型油菜Brassica napus、高粱Sorghum bicolor等)。结论：基于文献与本研究发现，用合理数量标记可从温带农业生态系统蚕豆群体中鉴定高近交个体且假发现率低；初筛步骤可再降基因分型成本达46%；本方法使蚕豆育种无需DH技术即可快速、廉价获取自交系，大幅缩短传统自交系选育周期。