《Genome Biology》：Dissecting the genetic basis of yield stability in faba bean by multi-environment analysis
编辑推荐：
背景：蚕豆是一种全球适应的豆科蛋白作物，具有高产潜力。目前，跨环境的产量变异限制了更广泛种植，其潜在遗传学仍知之甚少。结果：研究人员鉴定了蚕豆产量和产量稳定性的主要QTL（数量性状位点）。研究人员对ProFaba多样性群体进行高分辨率基因分型，并在欧洲开展了协
背景：蚕豆是一种全球适应的豆科蛋白作物，具有高产潜力。目前，跨环境的产量变异限制了更广泛种植，其潜在遗传学仍知之甚少。结果：研究人员鉴定了蚕豆产量和产量稳定性的主要QTL（数量性状位点）。研究人员对ProFaba多样性群体进行高分辨率基因分型，并在欧洲开展了协调的多年度/地点试验。基于这些数据，研究人员鉴定了超过一百个与14个复杂性状（包括产量）的平均表现和稳定性相关的位点。实验验证支持候选基因Vfaba.Hedin2.R2.1g002122参与植株结构，其基因表达与第一荚位和株高显著相关。此外，研究人员引入了一种基于随机回归混合模型（random regression mixed model）整合环境数据以分析性状稳定性的方法，该方法能够预测在未测试环境中的表现。结论：研究提供了对蚕豆产量、产量稳定性和基因型-环境互作（GxE）遗传结构的见解。基因组资源、候选位点和基于天气的分析框架为预测跨环境表现和加速培育抗逆高产蛋白作物提供了实用工具。
蚕豆（Vicia faba L.）是一种全球适应的豆科蛋白作物，具有高产潜力，但产量稳定性受环境变异限制，阻碍了其更广泛种植和可持续农业系统的应用。此前，产量变异的潜在遗传学仍知之甚少，尤其是基因型-环境互作（GxE）对复杂性状的影响缺乏系统解析。针对这一空白，研究人员在《Genome Biology》发表研究，通过多环境分析和高分辨率基因分型，系统解析了蚕豆产量稳定性的遗传基础。研究利用ProFaba多样性群体（包含234个来自欧洲、北非、拉丁美洲、中东、远东和澳大利亚的品系，覆盖春、冬类型及三种种子大小类），结合基因分型测序（GBS）和单引物富集技术（SPET）获得54万个高质量SNP标记，并在欧洲5个国家（丹麦、爱尔兰、法国、德国、西班牙）的12个环境中开展协调田间试验，评估了14个复杂性状（包括种子产量、形态、生物胁迫等）。通过引入天气衍生的合成协变量（synCov）构建反应范模型，研究揭示了产量稳定性的遗传结构，鉴定了多个候选基因，并实现了对未测试环境中基因型表现的有效预测。该研究为作物育种提供了重要基因组资源和分析框架。
本研究采用的主要关键技术方法包括：1）基因分型：利用GBS和SPET技术对ProFaba多样性群体进行高分辨率基因分型，结合10倍全基因组重测序（WGS）数据验证，获得54万个高质量SNP标记；2）群体结构分析：通过ADMIXTURE和主成分分析（PCA）鉴定4个亚群，并评估连锁不平衡（LD）衰减；3）反应范模型：分别使用未知协变量（unkCov）和天气衍生的合成协变量（synCov，基于主成分分析构建，代表冷湿到暖干环境梯度）拟合随机回归混合模型，提取截距和斜率；4）全基因组关联研究（GWAS）：采用单阶段G×E交互模型（线性混合模型）、贝叶斯稀疏线性混合模型（BSLMM）、BLINK方法以及基于基因和1 Mb窗口的集合分析，鉴定与平均表现和稳定性相关的位点；5）候选基因验证：对8个候选基因进行RT-qPCR表达分析，检测与表型关联；6）基因组预测：使用GBLUP模型进行5倍交叉验证和留一环境交叉验证，评估预测能力。样本队列来源为ProFaba多样性群体，田间试验在欧洲多国协作完成。
研究结果部分主要发现如下：
（1）The ProFaba panel：通过ADMIXTURE分析鉴定出4个亚群（sp1-sp4），分别对应地中海、中欧西欧、亚洲和法国德国品系。PCA和基因组关系网络显示亚群间分离，LD衰减至半最大值约44.08 kb，全基因组核苷酸多样性（π）为0.25，sp3略低。该群体遗传多样性适宜遗传参数估计和基因发现。
（2）Unsupervised learning of environmental relationships：通过主成分分析9个天气变量（温度、土壤湿度、降水、辐射等）构建合成协变量synCov，第一主成分解释43.4%方差，代表从冷湿到暖干的环境梯度。将实际试验地点投影到该空间，提取synCov作为反应范模型的环境指标。
（3）Variance explained by genetic factors：通过BSLMM估计SNP标记解释的表型方差比例（PVE）和稀疏效应解释的方差比例（PGE）。截距PVE范围为0.70-0.99，PGE为0.21-0.83。产量稳定性的稀疏效应明显，约50个大效应SNP解释超过一半的遗传方差，表明其遗传结构相对简单。利用未知协变量和合成协变量均检测到显著的GxE，尤其对于锈病、活力、种子产量等性状，且遗传方差随环境梯度变化（尺度型GxE）。通过遗传回归校正后，发现大多数性状存在秩次型GxE（交叉互作），例如sp1（南欧）品系在暖干环境中产量提高，而sp2和sp3相反。
（4）Single-stage genome-wide scan of QTLxE interaction：通过单阶段交互模型（LMM）鉴定出43个显著SNP，多数显示条件中性效应（仅部分环境显著），14个显示拮抗多效性（效应方向随环境改变）。关键信号包括：位于FR-1位点的基因Vfaba.Hedin2.R2.1g002127（编码DREB转录因子），其冬季类型等位基因在冷湿环境（爱尔兰、德国）提高产量，但在暖干环境（西班牙）降低产量，与FR-1位点耐冻性功能一致；另一个信号位于编码果胶甲酯酶抑制剂的基因，与锈病抗性相关；此外，位于FLOWERING LOCUS T（FT）基因附近的SNP也显示拮抗多效性，影响花后株高。
（5）Association analysis of reaction norm coefficients：对反应范模型的截距和斜率（来自unkCov和synCov）进行GWAS，共鉴定75个显著关联（69个独特SNP）。截距相关关联最多，包括59个候选基因。产量平均表现方面，发现3个显著SNP，分别解释16.42%、20.68%和19.80%的方差。稳定性方面，利用unkCov斜率检测到22个关联，其中最显著峰位于染色体1a，单一基因Vfaba.Hedin2.R2.1g004519（编码MAIN-LIKE 2蛋白）解释43.77%的稳定性方差；利用synCov斜率检测到12个关联，其中位于TUBBY-like F-box蛋白基因的SNP解释近50%方差，携带该等位基因的品系（主要来自sp1）在暖干环境中产量高于群体平均。集合分析进一步验证了FT基因（Vfaba.Hedin2.R2.1g004008）对开花时间的显著贡献，该基因位于染色体1a，与已知染色体5的FT基因相似，系PEBP基因家族重复。表达分析显示，候选基因Vfaba.Hedin2.R2.1g002122（编码ABC转运蛋白ABCB28）的表达与第一荚位和株高显著负相关，支持其调节植株结构的作用。
（6）Genome-wide prediction of mean performance and stability：通过GBLUP模型进行5倍交叉验证，截距预测能力为0.50-0.86，斜率预测能力较低但仍有较高值（如产量未知协变量斜率0.76，合成协变量斜率0.55）。留一环境交叉验证显示，基于合成协变量的反应范模型与标准乘法反应范模型预测能力相当，但能直接预测未测试环境表现。以2020年为例，利用synCov截距和斜率绘制四个代表性品系的生产潜力地图，直观展示了适应区域。
讨论部分总结：本研究系统解析了蚕豆产量和稳定性遗传结构，鉴定了多个与GxE相关的QTL，其中FR-1位点等关键基因在环境特异性表现中起重要作用。通过引入天气衍生的合成协变量，克服了传统反应范模型依赖环境均值的局限，实现了对未测试环境的预测。该框架为加速培育适应气候变化的抗逆高产蛋白作物提供了实用工具。
结论（翻译）：该研究推进了对蚕豆性状和性状稳定性遗传基础的理解。它还提供了一个将复杂气候数据纳入反应范模型的通用框架，取代环境均值，从而允许预测未测试环境中的基因型表现。这现在可以启动对蚕豆和其他产量受环境变异（包括预测的未来气候变化）挑战的作物进行产量稳定性改良的育种。
