武汉植物园在菱属泛基因组和结构变异研究中获进展

　　水稻、玉米和小麦是三大主粮作物，但主粮来源单一成为粮食安全和营养安全的隐患，因此对“孤儿作物”（orphan crop）的重视利用是应对粮食安全危机的有效手段之一。然而，随着城乡现代化进程的持续推进和现代集约化农业的快速发展，“孤儿作物”野生种和地方品种资源正日趋减少，与主粮作物相比，其种质资源的研究及现代育种技术的应用滞后。菱角是典型的“孤儿作物”。新石器时代以来，长江流域的先民开始有意识地大规模采集、驯化菱角。在南宋时期，菱角已成为江南地区的主粮之一，仅太湖地区便形成了包括“乌菱”“南湖菱”等20多个栽培品种。随着基因组学技术的发展，泛基因组研究已成为探讨作物重要性状形成机制的重要手段，被广泛应用于作物育种改良。

　　近日，中国科学院武汉植物园东亚植物演化、保护与利用学科组研究员邱英雄团队，在《园艺研究》（Horticulture Research）上，以Research article形式，在线发表了题为Pangenome of water caltrop reveals structural variations and asymmetric subgenome divergence after allopolyploidization的研究论文。

　　菱属（TrapaL.）包含欧菱（T. natans）和细果野菱（T.incisa）两个种。其中，欧菱包括异源四倍体（AABB）与二倍体（AA）两种倍性，而细果野菱（BB）为二倍体，栽培菱角即驯化于形态较大的二倍体欧菱（图1）。本研究组装了菱属两个二倍体物种的高质量基因组，结合已发表的四个（亚）基因组，构建了基于基因家族的泛基因组（图2）。结果显示：在全部基因家族中，核心基因家族（coregene clusters）、可变基因家族（dispensable gene clusters）与私有基因家族（private gene clusters），分别占所有基因家族的48.05%、28.92%和23.03%。同时，该研究以栽培南湖菱（欧菱）基因组（TnA_NL）为骨架，基于211,598个非冗余的插入缺失变异（PAVs）构建了菱属的图形化基因组。这种图形化基因组为结构变异的精确分型提供了平台。为了探究PAVs是否与二倍体种之间的表型差异有关，研究进一步鉴定到40,453个PAVs在A/B两种基因组类型间完全分化（这些PAVs涉及2,570个基因）。GO富集分析发现，这些基因的功能主要包括器官生长发育、有机物代谢过程、响应外界刺激，且在多个器官中呈现差异表达，表明这些基因可能在二倍体种间表型分化和生殖隔离方面起到关键作用（图3）。此外，研究基于亚基因组和亲本基因组之间的比较和泛转座子（TE）分析发现，异源四倍体欧菱经历了亚基因组不对称演化。其中，B亚基因组为显性亚基因组，而不平衡的PAVs、TE扩增、部分同源互换以及基因表达差异等因素共同驱动异源多倍体的亚基因组不对称演化（图4）。

　　菱属泛基因组构建和结构变异鉴定，为菱角研究提供了重要的资源和平台。上述研究加深了科学家对多倍化过程中基因组演化的认知，并为菱角种质资源的遗传评价和种质创新奠定了基础。研究工作得到浙江省金华市农业科学研究院地合作项目的支持。该研究由武汉植物园和浙江省金华市农业科学研究院合作完成。

　　图1.菱属多样化的果实及演化路线

　　图2.菱属的泛基因组组成及特征