《Plants》：Genome-Wide Identification and Abiotic Stress Response Analysis of the Isopentenyl Transferase (IPT) Gene Family in Soybean (Glycine max L.)
Zhihao Zhang,
Haorang Wang,
Mujeeb Ur Rehman,
Chunling Pei,
Yongzhe Gu,
Yingpeng Han and
Lijuan Qiu
编辑推荐：
本文通过对中国和乌拉圭大豆种质进行比较基因组学分析，系统揭示了导致两地大豆含油量差异的遗传学基础。研究鉴定出GmDGAT1和GmbZIP123等关键基因，其有利单倍型在乌拉圭高油种质中显著富集，这为通过分子模块精准导入来提升我国大豆的油脂产量提供了明确靶点和育种路径。
本研究聚焦于大豆（Glycine max L.）异戊烯基转移酶（Isopentenyltransferase, IPT）基因家族的系统鉴定与功能分析。IPT是细胞分裂素生物合成途径中的限速酶，在植物适应非生物胁迫过程中扮演关键角色。研究首先通过全基因组水平鉴定，在大豆中识别出IPT家族成员，并对其进行全面的生物信息学分析。结合多环境表型数据，评估了该基因家族成员在大豆应对干旱、盐碱等非生物胁迫条件下的表达模式与功能关联，旨在为解析大豆抗逆性的分子机制及分子设计育种提供基因资源。
1. 引言
大豆是全球最重要的经济作物之一，提供约28%的世界植物油和70%的膳食蛋白质。大豆原产于中国，由野生大豆在大约5000年前驯化而来。随着全球生物能源产业的快速发展和植物油需求的持续增长，提高种子含油量已成为现代大豆育种的核心目标。尽管中国拥有全球最丰富的野生大豆和地方品种资源，但近半个世纪以来，大豆的主要产区和育种中心已转移至美洲。长期的地理隔离和异质的育种目标导致了中国与国外种质在农艺性状上的显著分化。历史上，受饮食文化影响，中国育种优先考虑高蛋白质含量和适合豆制品加工的品质。相比之下，美洲的育种体系主要服务于工业榨油和饲料生产，推动了对高产量和高出油率的强烈选择。这种差异导致了含油量的显著分化，中国品种的含油量普遍比美洲商品品种低2-3个百分点。本研究旨在阐明中乌种质间含油量分化的遗传基础，解析人工选择的基因组印记，并为分子设计育种提供理论依据。
2. 结果
2.1. 中乌大豆群体种子蛋白与油分含量变异
研究分析了三个不同大豆群体：中国种质、乌拉圭常规种质和乌拉圭转基因种质的种子蛋白和油分含量。统计结果显示，中国群体平均蛋白含量最高，为43.61%（±1.18），显著高于两个乌拉圭群体。乌拉圭常规种质平均蛋白含量为40.24%（±1.18），而转基因种质最低，为36.30%（±2.14）。与此相反，油分含量则呈现出相反的趋势，这与大豆种子中油分与蛋白质积累普遍存在的负相关性一致。两个乌拉圭群体的平均油分含量显著高于中国种质，分别为21.56%（±0.55）和21.41%（±1.03），而中国种质的平均油分含量为19.42%（±1.05）。相关性分析证实，中乌群体间的种子蛋白与油分含量存在极显著的负相关。
2.2. 油分含量的广义遗传力分析
对三个群体的油分含量进行广义遗传力分析，结果显示，中国群体、乌拉圭常规群体和转基因群体的广义遗传力估计值分别为0.87、0.82和0.98，表明该性状在不同环境下具有高度的遗传稳定性。当将72份常规种质（合并48份中国种质和24份乌拉圭常规种质）作为一个“合并”群体分析时，广义遗传力增加到0.92。这反映了两个基因库在油分性状上存在显著的遗传分化，表型差异主要由可固定的遗传因素驱动，而非环境可塑性。
2.3. 中乌大豆种质间油分与蛋白质含量的表型分化
主成分分析进一步可视化了三个大豆群体的整体表型结构及其分化。PC1解释了绝大部分变异（90.23%），主要反映了蛋白质与油分含量之间的强负相关。沿PC1轴，中国种质与乌拉圭群体表现出明显的分离。中国种质集中在PC1负向一侧（高蛋白/低油），而乌拉圭群体则集中在正向一侧（高油/低蛋白）。PC2解释了剩余的9.77%变异。乌拉圭转基因与常规种质群体高度重叠且紧密聚类，这表明乌拉圭品系共享的遗传背景是决定种子品质性状的主要因素，转基因事件的渗入并未导致其在蛋白和油分组成上与常规轮回亲本产生显著的表型分化。
2.4. 中乌种质遗传多样性的全基因组分布
基于158,327个SNP标记，绘制了三个群体全基因组期望杂合度的分布。中国种质展现出最高的遗传多样性水平，平均He为0.198，其分布曲线相对宽而平，表明在起源中心的基础亲本中保持了丰富的等位基因变异库。相比之下，两个乌拉圭群体的遗传多样性均显著降低，乌拉圭转基因群体表现出最低的遗传多样性（平均He = 0.135），低于乌拉圭常规群体（平均He = 0.151）。这种模式反映了大豆引入南美过程中的“奠基者效应”和遗传瓶颈。转基因群体多样性的进一步降低，表明转基因事件（如草甘膦耐受性）的渗入及后续对商业性状的严格选择施加了额外的遗传瓶颈，导致更大比例的基因组位点被固定。
2.5. 全基因组选择信号与关键油分QTL GmDGAT1分析
为探究与育种历史和适应性相关的基因组分化与选择信号，研究计算了中国与乌拉圭大豆种质间的群体分化指数。全基因组扫描将排名前1%的FST窗口定义为候选选择清除区域。在17号染色体上鉴定到一个与油分含量相关的极显著选择信号（FST= 0.51）。这个高峰区域包含了GmDGAT1，这是大豆油脂合成途径中的一个关键基因。对该区域进行的局部连锁不平衡分析显示，乌拉圭群体在GmDGAT1及其侧翼区域存在广泛的高r2LD区块，表明该区域经历了强烈的人工选择，导致优良单倍型快速固定，遗传多样性显著降低。这些基因组层面的发现支持了GmDGAT1是中乌大豆种质油分差异主要贡献者的假说。
2.6. 关键油脂与适应性基因的单倍型分化与有利等位基因聚合分析
研究筛选出两个高置信度的候选基因进行详细的单倍型分析：GmDGAT1和GmbZIP123。在GmDGAT1位点，鉴定出三个主要单倍型。在中国种质中，Hap2是优势单倍型（频率0.650），而Hap1频率较低（0.257）。相比之下，Hap1在乌拉圭常规种质中频率显著升高（0.783），并在乌拉圭转基因种质中达到完全固定（1.000）。在GmbZIP123位点，中国种质以Hap5为主（频率0.750），Hap4频率较低（0.170）。与此形成鲜明对比的是，Hap4在乌拉圭常规种质中完全固定（1.000），在转基因种质中也占主导（0.880）。这些结果表明，乌拉圭种质在多基因位点上存在高度一致的等位基因富集趋势。在长期的集约化育种过程中，乌拉圭种质不仅被直接选择高油基因GmDGAT1，同时也可能固定了对增强代谢适应性有益的GmbZIP123等位基因。这种聚合“高油+高适应性”优良等位基因的模式，构成了其优良综合农艺性状的重要遗传基础。
2.7. 中乌群体间关键油分基因的等位基因分化
对文献汇总的42个已知油脂相关基因座进行统计，比较了中国与乌拉圭大豆种质间参与油脂生物合成与调控的候选基因等位基因频率。结果显示，不同位点间的分化程度存在显著变异。其中，GmDGAT1基因关联位点的频率变化最为剧烈。在该基因的一个核心SNP位点上，有利等位基因频率在中国基础种质中相对较低（约0.26），但在引进的乌拉圭种质中接近固定（约0.90-1.00）。这种鲜明对比表明，GmDGAT1在大豆引入乌拉圭或当地高强度高油育种过程中经历了强烈的正向选择，这很可能是导致这些优良种质油分含量更高的一个关键遗传因素。
3. 讨论
3.1. 驯化选择下油分性状的演化与遗传分化
本研究量化了中乌大豆种质之间显著的油分差异（约2.1%），证实了在南美新兴大豆产区已经发生了显著的遗传分化。这种表型分化深刻地反映了两地不同的驯化和育种历史。由于大豆种子中蛋白质和油分含量存在显著的负相关，由能量代谢约束驱动的这种“进化权衡”意味着对高蛋白的选择往往伴随着低油基因的保留。乌拉圭育种体系深受北美模式影响，主要受工业榨油效率驱动。对高油率的追求打破了原有的蛋白-油脂平衡，导致了高油基因的快速积累。本研究中的FST数据强有力地支持了这一观点：乌拉圭种质在油分相关位点上的遗传分化。这表明这种分化是高强度定向人工选择的结果，而非随机遗传漂变。
3.2. GmDGAT1的选择印记与功能演化
本研究检测到多个染色体上存在显著的分化信号。尽管全基因组范围内最高的分化峰位于其他区域，但本研究所关注的油分性状在17号染色体的GmDGAT1区域显示出特异且显著的选择信号。GmDGAT1是催化三酰甘油生物合成最后一步的关键限速酶。该位点的显著分化为中乌种质间油分含量的显著表型差异提供了直接的功能基因组学证据。本研究发现GmDGAT1基因在乌拉圭种质中接近固定。乌拉圭种质很可能通过引入北美高油亲本，成功固定了GmDGAT1的高油单倍型。由于长期高油选择压力，GmDGAT1的高油单倍型在乌拉圭群体中被系统地保留和富集。相比之下，该有利等位基因在中国种质中频率相对较低，这解释了中国品种普遍较低的油分含量，并意味着该位点在中国育种历史上尚未受到足够的正向选择压力。
3.3. 遗传瓶颈效应下的多样性丧失与遗传增益
尽管乌拉圭种质表现出极高的油分含量，但其遗传多样性显著降低。这种源自少数优良亲本的“奠基者效应”在美洲大豆中普遍存在。遗传多样性的丧失，即“遗传阻力”，在固定优良性状的同时，也可能增加作物对极端气候的脆弱性。相比之下，中国种质保留了更丰富的变异库，为应对未来气候变化提供了宝贵的遗传缓冲。本研究还发现乌拉圭常规与转基因种质之间遗传差异极小，这表明当地的转基因育种主要是将抗性基因导入到少数优良高油的遗传背景中。
3.4. 基于基因组分化的分子育种与种质改良策略
基于中乌大豆种质间鉴定的显著遗传分化区域和关键固定等位基因，研究提出了整合的分子育种策略以提升中国大豆的油分含量。首先，应优先发展精准的标记辅助选择。针对GmDGAT1等区域开发高通量KASP或SNP单倍型标记，将有助于高效选择携带乌拉圭优良等位基因的个体。鉴于油分与蛋白质含量之间的强负相关，应利用基因组选择模型来打破高油与低蛋白基因之间的连锁累赘。此外，基因编辑技术为改良提供了直接途径。鉴于GmDGAT1经历了强烈的人工选择，利用CRISPR/Cas9技术修饰中国品种中的顺式调控元件，以模拟高油单倍型的调控模式，可望快速提升油分含量。总之，乌拉圭种质资源是拓宽中国大豆遗传背景、改良油分性状的重要基因库。反之，中国种质在丰富乌拉圭大豆的高蛋白特性方面扮演着关键角色。
4. 材料与方法
4.1. 植物材料与表型调查
本研究构建了一个包含98份大豆种质材料的代表性核心集合。该集合包括48份代表出口级中国大豆的骨干亲本，以及50份代表乌拉圭当前农业生产主力的核心商业品种。田间试验采用随机完全区组设计，在海南、浙江和江西三地进行，设两次重复。收获后，使用近红外光谱分析仪测定种子的粗脂肪和粗蛋白含量。通过线性混合模型估计广义遗传力，并计算方差分析。
4.2. 基因分型与群体遗传分析
使用“中豆芯1号”SNP芯片对乌拉圭新引进种质和出口种质进行基因分型，平均获得158,327个变异位点。使用PLINK 1.9过滤低质量标记，利用主成分分析进行群体结构分析，并使用VCFtools计算群体间的遗传分化指数和核苷酸多样性，以100-kb窗口和10-kb步长扫描全基因组选择信号。
4.3. 单倍型与有利等位基因分析
基于芯片数据，提取关键候选基因区域内的SNP。结合文献整理的42个油脂相关基因列表，统计每份种质材料携带的“高油”有利等位基因总数，并使用t检验比较群体间差异。
5. 结论
本研究证实，乌拉圭大豆种质代表了一个由高强度定向选择形成的、区别于中国起源中心的专业化高油基因库。研究鉴定出GmDGAT1和GmbZIP123为关键靶点，其有利单倍型在乌拉圭品系中已固定，但在中国种质中仍存在变异。这些发现意味着，两地群体间的油分差距是由特定、可操作的遗传因素驱动的。因此，我们建议将乌拉圭种质作为“预组装”高油遗传模块的资源库加以利用。未来的育种必须超越简单的杂交，聚焦于将这些模块精准导入适应性广泛的中国遗传背景中，利用标记辅助选择或基因编辑等策略。这一策略为打破蛋白质与油分之间的负相关，最终保障更高的植物油生产效率提供了可行的路线图。
