《New Crops》：Meiotic recombination assembles complementary parental alleles in an elite soybean line
编辑推荐：
本研究旨在阐明大豆（Glycine max）优良性状形成的分子基础。研究人员采用了一个亲本-子代 trio 策略，包含两个亲本品系及其优良后代品系 HX3。通过对 trio 进行全基因组重测序，重建了 HX3 的减数分裂重组图谱和高分辨率基因组结构；同时结合转
本研究旨在阐明大豆（Glycine max）优良性状形成的分子基础。研究人员采用了一个亲本-子代 trio 策略，包含两个亲本品系及其优良后代品系 HX3。通过对 trio 进行全基因组重测序，重建了 HX3 的减数分裂重组图谱和高分辨率基因组结构；同时结合转录组分析，整合基因型、表型和基因表达数据，追溯了与关键性状相关的候选基因的 parental origin，并分析了非加性基因表达模式与优良表型的关联。
测序结果显示，亲本品系桂早1号（GZ1）和巴西13（BX13）间存在 1,260,418 个高质量单核苷酸多态性（SNP）位点，平均间距 752 bp。GZ1 携带的 SNP、小片段插入缺失（indel）和结构变异（SV）数量均多于 BX13。在蛋白编码序列（CDS）中，GZ1 有 52,137 个 SNP（3.4%），BX13 有 41,597 个（3.5%），其中约 57% 为非同义 SNP。此外，GZ1 和 BX13 分别检测到 225,786 个和 178,943 个 small indel，以及 108,457 个和 87,398 个 SV，其中以重复序列为主（分别占 84.5% 和 95.4%）。 Direct comparison 发现 25,052 个非同义 SNP 影响 10,878 个基因，其中 243 个基因含 10 个以上非同义 SNP，被定义为高分化基因。Gene Ontology（GO）富集分析显示，这些高分化基因显著富集于端粒维持（telomere maintenance）、DNA 修复（DNA repair）、信号转导（signal transduction）、DNA 解旋酶活性（DNA helicase activity）和 ADP 结合（ADP binding）等功能。通过计算 Ka/Ks 比值，184 个高分化基因中 70 个表现出强纯化选择（Ka/Ks < 0.3），其余 114 个显示出松弛的负选择（Ka/Ks 0.3-1）。
基于 GZ1 与 BX13 间的 SNP 标记，HX3 基因组被重建为 55.8% 母源（GZ1）和 41.0% 父源（BX13）的镶嵌结构（mosaic structure），并检测到 84 个重组交换（crossover, CO）事件，涉及 308 个基因。CO 断点区间大小不一（100 bp – 316 kb），大多位于基因丰富、转座子贫乏的区域。此外，2.9% 的基因组为杂合区域，包含 402 个基因，可能代表选择压力下保留的残余杂合性。
针对 89 个已知功能的重要农艺性状标记基因进行基因型溯源发现：调控株高、主茎节数和开花时间的基因（如 E2、GmGID1-2、GmGBP1、GmLUX2、GmEID1、GmZTL3、GmFT2a/2b 等）主要从父本 BX13 遗传而来；而与光合作用、花色、籽粒品质和抗逆性相关的基因（如 GmCCD4、GmRCA11、atpG、Glyma.18G107000、GmSW14、GmAPX6、GmCOL1a、GmLMM2-1、Rsv3 等）则主要来源于母本 GZ1。
转录组分析显示，HX3 中 59.2% 的注释基因表达，鉴定出 1,944 个差异表达基因（differentially expressed genes, DEGs），占注释基因的 3.5%。通过 k-means 聚类分析，DEGs 被分为 9 个共表达簇（cluster）。其中 383 个基因（cluster C1）在 HX3 中呈超显性表达（overdominant expression），显著高于双亲，富集于细胞壁组织、发酵、糖酵解、碳水化合物代谢、氮代谢、光合作用和四吡咯生物合成等功能，代表基因包括调控株高和开花时间的 E2。147 个基因（clusters C5 和 C8）呈低显性表达（underdominant expression），在 HX3 中表达低于双亲，富集于 RNA 代谢、次生代谢、信号转导和应激响应等功能，代表基因包括盐胁迫响应基因 GmSST1 和结瘤共生相关基因 GmCHI1B。546 个基因（clusters C2、C4、C7）呈现母源样表达模式，在 GZ1 和 HX3 中表达相似，富集于发育、脂质代谢、次生代谢、应激响应、光合作用、重金属处理等功能，代表基因包括花色调控基因 GmCCD4、光合作用相关基因 psb27 抗性基因 Rsv3 等。493 个基因（clusters C3、C6）呈现父源样表达模式，在 BX13 和 HX3 中表达相似，富集于细胞过程、激素代谢、脂质代谢、次生代谢和应激响应等功能，代表基因包括开花调控转录因子 GmFUL2a/GmTof5 和 GmMDE6，以及低磷胁迫响应基因 GmPAP12、GmPAP17 和 GmWRKY6 等。315 个基因（cluster C9）在双亲间差异表达但在 HX3 中呈现亲本中间型或偏向某一亲本的表达模式，富集于信号转导、应激响应、转运和金属处理等功能，代表基因包括磷效率相关基因 GmVPE1、GmVPE2 和 GmPHT1。
整合基因组和转录组证据发现：E2 基因及其表达模式与 BX13 更为接近，这与 HX3 在株高和开花时间上表现 BX13 样表型一致；光合作用相关基因（psb27、atpG、GmRCA11）及其表达模式与 GZ1 更为接近，这与 HX3 叶绿素 a 含量与 GZ1 相似的表型一致； Grain weight 和 oil content 相关基因 Glyma.18G107000 的母源遗传和高表达，与 HX3 在这些性状上的母源样表现一致；GmTHI1（GmPGL1）的母源样表达、GmFUL2a 和 GmMDE6 的父源样表达、以及 GmBEN1 的低表达模式，均与相应表型特征吻合。
这些发现深化了对减数分裂重组如何组装互补亲本等位基因以形成优良表型的理解，并为大豆及其他作物的功能基因组学和分子育种提供了有价值的基因组资源和分析框架。
本研究围绕大豆优良品系 HX3 的遗传构成展开，采用亲本-子代 trio 策略，系统解析了减数分裂重组在组装互补亲本等位基因、形成优良农艺性状中的作用机制，为理解作物育种中的遗传重组规律提供了重要见解。
研究背景方面，作物杂交育种依赖于减数分裂重组（meiotic recombination）产生新的等位基因组合，从而聚合优良农艺性状。然而，重组在基因组和转录组层面塑造优良表型的具体机制仍知之甚少。大豆作为世界上最重要的作物之一，其约 1 Gb 的基因组经历了至少两轮全基因组复制事件，导致大量序列冗余，增加了复杂性状遗传基础解析的难度。尽管基于群体的数量性状位点（quantitative trait locus, QTL）鉴定和分子标记开发已取得进展，但群体方法在识别衍生后代中具体等位基因及其组合方面存在局限。因此，亲本-子代 trio 策略为追溯优良等位基因起源、表征重组图谱和阐明有利等位基因组合的形成机制提供了有效途径。
研究人员以母本桂早1号（GZ1）和父本巴西13（BX13）及其优良后代 HX3 构成的 trio 为研究对象，系统开展了以下研究：首先，对三个品系进行全基因组重测序，重建 HX3 的减数分裂重组图谱和高分辨率基因组结构；其次，通过转录组测序比较三者的基因表达谱；最后，整合基因组、表型和基因表达数据，追溯关键性状相关位点的亲本起源，解析优良性状的遗传基础。
研究结果表明，HX3 在株高、主茎节数和开花时间等性状上超亲分离（transgressive segregation），表现出优于双亲的表型，而在百粒重、叶绿素 a 含量、可育荚比例、籽粒蛋白质和油分含量等性状上则接近或介于双亲之间。基因组分析揭示，HX3 基因组中 55.8% 来源于母本 GZ1，41.0% 来源于父本 BX13，另有 2.9% 的杂合区域。共检测到 84 个重组交换事件，断点区间大小从 100 bp 到 316 kb 不等，大多位于基因丰富、转座子贫乏的常染色质区域。
基于 89 个重要农艺性状标记基因的基因型溯源分析，研究揭示了功能分区化的亲本贡献模式：调控株高、主茎节数和开花时间的基因主要从父本 BX13 遗传，包括 E2、GmGID1-2、GmGBP1、GmLUX2、GmEID1、GmZTL3 和 GmFT2a/2b 等；而与光合作用、籽粒品质（粒重和油分）、花色、抗逆性（盐旱胁迫、疫霉根腐病抗性、大豆花叶病毒抗性）相关的基因则主要从母本 GZ1 遗传，包括 GmCCD4、GmRCA11、atpG、Glyma.18G107000 和 GmSW14、GmAPX6、GmCOL1a、GmLMM2-1 和 Rsv3 等。此外，部分低磷胁迫响应基因如 GmACP1 和 GmPHR25 则来源于父本 BX13。
转录组分析进一步揭示了广泛的非加性基因表达现象。在 1,944 个差异表达基因中，383 个基因在 HX3 中呈超显性表达，显著高于双亲，富集于细胞壁组织、糖酵解、光合作用等功能，代表基因 E2 的表达模式与父本 BX13 一致；147 个基因呈低显性表达，富集于 RNA 代谢和应激响应等功能；546 个基因呈现母源样表达模式，493 个基因呈现父源样表达模式，315 个基因呈亲本中间型表达模式。整合分析证实，基因型、表达模式和表型特征高度一致：光合作用相关基因的母源遗传和高表达与 HX3 叶绿素 a 含量的母源样表现一致；籽粒品质和油分相关基因的母源贡献与其表型一致；而株高和开花时间相关基因的父源遗传则与其相应的父源样表型吻合。
讨论部分指出，现代作物育种依赖于减数分裂重组产生有利等位基因的新组合，但重组塑造优良表型的遗传机制尚不完全清楚。本研究通过 trio 策略重建了 HX3 的高分辨率基因型图谱，描绘了重组产生的基因组镶嵌结构，并追溯了与优良性状相关的候选基因的 parental origin。研究发现，父源贡献主要体现在株高和开花时间等植物结构相关性状，而母源贡献则主要体现在光合作用效率和籽粒品质等性状，暗示了亲本等位基因的功能分区化贡献。同时，多个抗逆相关基因也主要来源于母本。研究人员也指出了方法的局限性：转录组分析仅基于叶片组织，而许多重要农艺性状（如花序转变、产量形成、种子品质）由分生组织、花器官、发育中的荚果和种子等组织特异性基因调控，因此有限的组织取样可能遗漏了部分重要的性状相关位点。此外，由于转录丰度低、取样发育时期或技术阈值等原因，某些功能基因可能未被检测为差异表达，但其功能重要性不可忽视。
研究结论表明，减数分裂重组通过重组交换重塑了亲本单倍型（haplotype），将来自不同亲本的互补优势等位基因聚合于后代基因组中，从而组装出具有优良农艺性状的基因型。父本来源的等位基因主要贡献于 HX3 的株高和开花时间调控，而母本来源的等位基因则主要贡献于光合作用效率、籽粒品质（粒重和油分）以及多种生物和非生物胁迫抗性。这些发现为理解作物育种中的遗传重组规律提供了新见解，并为大豆及其他作物的分子设计育种提供了有价值的基因组资源和可广泛应用的分析框架。
