NC重磅|豇豆泛基因组解析农艺性状
文章导读：
非洲正面临严峻的粮食短缺与营养不良问题，适应本土环境的乡土作物对保障粮食安全至关重要。豇豆耐热、耐旱且富含蛋白质，是非洲重要豆科作物。然而，其丰富的遗传多样性仍未被充分挖掘，限制了高产、优质品种培育。
中国农业科学院研究团队于2026年5月20日在国际著名期刊Nature Communications发表标题为“Pan-genome analysis uncovers extensive structural variations in Africa’s indigenous legume crop cowpea”的研究性论文，本研究通过构建豇豆泛基因组，系统揭示了结构变异对其遗传多样性与农艺性状形成的驱动作用，为分子育种提供了理论基础。
主要研究成果
1.豇豆高质量基因组组装与注释
研究整合了全球五大基因库的20份代表性豇豆种质资源与6个已发表基因组，构建了由26个基因组组成的豇豆泛基因组。采用PacBio测序与Hifiasm组装技术，获得了多个高质量基因组（557.03–606.59 Mb），BUSCO完整性评估达98.70%，组装指标整体表现优异。以FC6的T2T基因组为参考，scaffold平均锚定率达到95.75%，Hi-C图谱分析进一步验证了染色体结构的高度一致性。基因注释覆盖90%以上功能数据库，转座元件在基因组中平均占比为42.26%。
图1：豇豆种质资源的地理分布与基因组特征分析
2.豇豆泛基因组构建与多样性分析
本研究构建的豇豆泛基因组共鉴定出29,557个直系同源基因簇，表现出典型的封闭型特征，表明豇豆基因组整体保守性较强。其中，核心基因占比73.2%，在所有26个基因组中普遍存在，其编码序列更长、遗传多样性更低、表达水平更高，主要参与基础生命活动。相比之下，可变基因与私有基因更多关联于生物与环境胁迫响应、次生代谢等途径，表明这些基因可能在豇豆适应环境和快速进化过程中起重要作用。
图2：26份豇豆种质资源的泛基因组分析
3.豇豆结构变异的基因组景观
本研究共鉴定出豇豆基因组中62,591个非冗余结构变异（SVs），包括插入、缺失、倒位等类型，主要分布于异染色质区。平均每个基因组包含约9,358个SVs，且泛SVs数量未达饱和。部分大片段SV可包含多个基因，涉及抗病、发育等重要功能。大多数SV为低频变异，并显著影响基因编码序列和表达。转座元件（TEs）特别是LTR与TIR类型与SV分布高度相关，表明其是SV形成的关键驱动因素。SV数量与重组率呈负相关，高重组区有助于清除有害SV。总体表明，SV广泛影响豇豆基因功能与适应性，其形成受转座活动与重组共同调控。
图3：豇豆结构变异图谱与泛结构变异分析
4.豇豆亚种分化相关结构变异
本研究基于26个基因组构建的图形化泛基因组，系统解析了结构变异（SV）在豇豆亚种分化与适应中的作用。通过分析619份样本，鉴定出超过6万个SVs和420万个SNPs。SV-GWAS结合选择性清除分析筛选出上千个候选基因，涉及胁迫响应、代谢调控等过程，部分基因与豆荚发育、开花时间等农艺性状相关。例如，VuELF4启动子缺失可能通过降低表达延迟开花，且该变异在亚种间频率差异显著。研究进一步发现，多个在亚种间固定的SVs与花发育、环境适应等功能相关，表明SV在豇豆适应性进化与人工选择中具有重要作用。
图4：鉴定导致豇豆亚种分化的结构变异
5.结构变异影响豇豆重要农艺性状
研究基于结构变异的全基因组关联分析鉴定出多个与荚长和单荚种子数相关的QTL位点。其中，在PL9.1区域的VuWAK基因启动子中存在的插入/缺失变异，能够调控基因表达水平，不同单倍型间荚长差异显著，且优良单倍型主要分布于Vs亚种，说明这些SV可能是导致亚种间荚长差异的重要原因。另一关键基因VuGA2ox2中的51 bp插入会改变蛋白质结构，降低赤霉素活性，进而抑制豆荚伸长。针对单荚种子数，研究在VuPPR基因中发现一处83 bp缺失，该缺失可形成提前终止密码子，显著降低种子数，且该变异主要存在于Vu亚种。综上所述，结构变异可通过影响基因表达与蛋白质功能，调控豇豆产量及形态性状。
图5：与豇豆荚长相关的结构变异
图6：豇豆单荚种子数相关结构变异的全基因组关联分析
文章小结
本研究成功构建了包含26个基因组的豇豆图结构泛基因组，通过基于结构变异的全基因组关联分析，研究鉴定出多个与亚种分化、荚长和单荚种子数显著相关的候选基因。这些结构变异可能通过调控基因表达或改变蛋白质功能，进而影响开花时间、豆荚伸长和种子形成等生物学过程。研究表明，结构变异是驱动豇豆遗传多样性和环境适应性的重要因素，为豇豆的遗传改良以及高产优质新品种的选育提供了重要的理论依据和分子资源。