《Plant Molecular Biology Reporter》：Genetic Diversity and Population Structure of Advanced Lima Bean Breeding Lines Revealed by SNP Markers
编辑推荐：
利马豆（Phaseolus lunatus L.）是具有较高农艺与营养价值的豆科作物，尤其对巴西东北部小农生产意义重大。巴西皮奥伊联邦大学（UFPI）育种计划通过双亲杂交选育具有决定型生长习性、早熟性及产量等农艺性状的品系。研究人员利用基因型分型测序（geno
利马豆（Phaseolus lunatus L.）是具有较高农艺与营养价值的豆科作物，尤其对巴西东北部小农生产意义重大。巴西皮奥伊联邦大学（UFPI）育种计划通过双亲杂交选育具有决定型生长习性、早熟性及产量等农艺性状的品系。研究人员利用基因型分型测序（genotyping-by-sequencing, GBS）产生的3002个单核苷酸多态性（single-nucleotide polymorphism, SNP）标记，评估34个第九世代（F9）品系的遗传多样性与群体结构。观测杂合度（observed heterozygosity, HO）低于期望杂合度（expected heterozygosity, HE），均值分别为0.01与0.19。品系H25-66拥有最多私有等位基因且遗传分化显著；H46群体品系高度同质，其余群体呈中等分化。主成分分析（principal component analysis, PCA）、邻接法（Neighbor-joining, NJ）聚类、热图及稀疏非负矩阵分解（sparse non-negative matrix factorization, sNMF）祖源图谱等多变量分析一致揭示明确遗传类群。高度同质组（如H46）与强分化材料（尤指H25-66）并存，表明该计划同时保有适于品种开发的稳定遗传背景与可用于亲本选配的差异化变异源，通过分化基因型杂交可创制新等位组合。上述发现为战略育种决策及加速适应巴西生产系统的改良利马豆品种培育提供支撑。
研究背景与立题依据
利马豆（Phaseolus lunatus L.）为自花授粉为主的多年生豆科作物，在巴西东北部由小农主导种植，2023年该国产量约10372吨，但全球范围内改良品种极少，巴西尚未正式释放任何利马豆栽培种。UFPI联合加州大学戴维斯分校（UC Davis）自2015年起以阿根廷、巴西、美国亲本进行双亲杂交，经改良混合选择至F9世代，前期仅完成形态农艺评价与遗传参数估计，缺乏分子水平遗传背景解析。研究人员假设连续自交与选择会降低高级品系多样性并形成可检出的群体结构，因此亟需利用基因组工具明确亲缘关系以指导亲本选配与拓宽遗传基础。
主要技术方法概述
研究人员以34个F9决定型生长习性品系为材料，这些品系源自5个双亲杂交群体（含全同胞、半同胞及非亲缘关系），种子由UFPI Phaseolus种质库提供，于ESALQ/USP田间单株取嫩叶，CTAB法提取DNA并经山梨醇缓冲液预洗去除多糖酚类；采用PstI-MspI双酶切GBS建库，Illumina NextSeq 2000单端测序；Stacks v2.62进行去复用与SNP发掘，比对至利马豆参考基因组（Garcia等2021），VCFtools过滤（最小深度5、缺失率控制）后保留3002个SNP；群体遗传参数由VCFtools与poppr计算，结构分析依次采用pheatmap热图、adegenet PCA、poppr NJ树（Rogers距离，1000次自举）及LEA包sNMF（K=1~5，交叉熵定K）。
研究结果
遗传多样性与SNP分布
GBS产生超4.19亿条reads，平均深度38.81×，Ts/Tv比1.23。3002个SNP在11条染色体呈非均匀分布，200 kb窗口密度热图显示基因丰富区标记密集、近着丝粒区稀疏。所有品系HO（均约0.016，范围0~0.08）远低于HE（均约0.19），符合F9近交纯合预期；H25-66私有等位基因数最高，H81-34与H94-29次之。
群体结构多变量分析
基于欧氏距离热图显示H46群内紧凑（除H46-50略偏中），H25-66整体距离最大，H81部分品系内部距离较高。PCA前两轴解释41.6%变异，H25-66与H81-34在散点中位置孤立，H46占据特定区间。NJ树中H46紧密成簇，H25亚簇中H25-66依旧最远，H81/H50-86/H94呈中等分化。sNMF交叉熵最低为K=3（0.395），K=3时三类群分别主要由H46、H25+H81（半同胞）、H94+H81+H50-86+H25混合组成；K=2与K=4下H46均保持均匀、H25-66始终为独立组分。
讨论部分总结
研究人员指出SNP染色体异质分布源于GBS偏好取样高重组基因区，与Garcia等观察一致；HO?HE印证连续自交致纯合化，且低于野生与栽培资源报道的HE=0.25，反映育种群体基础收窄。私有等位基因提示品系间分化由重组固定所致，H25-66、H81-34等可作广基亲本。多方法聚合显示同群体全同胞相似度高、H46因亲本表型对比仅涉少数主效位点而基因组均一，H25-66分化来自BGP UFPI 628×G25236组合特异等位固定。K=3结构对应初始亲本贡献与选育史，H46窄基利于品系稳定但需引入分化亲本重组。
结论部分翻译
UFPI育种计划的F9利马豆品系表现出低遗传多样性与由初始杂交、定向选择及连续自交塑造的明晰群体结构。GBS基因组范围SNP标记高效揭示了基因组变异，区分相似与分化品系并厘清育种群体关系。H46等高度同质组与H25-66等显著分化材料并存，证明该计划同时具备供品种开发用的稳定背景与供战略育种用的变异源。本研究为亲本选配与育种策略设计提供基因组依据，品系间遗传分化可辅助未来杂交组合拓宽遗传基础；但需配合配合力研究验证分化亲本农艺价值，并与形态、生理、营养及加工数据整合，以加速适应巴西条件的改良利马豆品种释放，强化这一被忽视作物的公共育种工作。论文发表于《Plant Molecular Biology Reporter》。
