《Plant Biotechnology Journal》：Comparative Genome-Wide Association Studies of Metabolites and Grain-Related Traits in Common Wheat
编辑推荐：
代谢组高度多样且是最接近表型的层面；因此，它通常被视为植物中基因组和表型组之间的桥梁。在本研究中，研究人员利用液相色谱-串联质谱（LC-MS/MS）法对一组多样化自然种质和一组重组自交系（RIL）群体进行了大规模代谢组分析，并测定了33个籽粒相关性状。研究人员
代谢组高度多样且是最接近表型的层面；因此，它通常被视为植物中基因组和表型组之间的桥梁。在本研究中，研究人员利用液相色谱-串联质谱（LC-MS/MS）法对一组多样化自然种质和一组重组自交系（RIL）群体进行了大规模代谢组分析，并测定了33个籽粒相关性状。研究人员确定了947种代谢物与33个籽粒相关性状之间的2286个关联的新网络。通过系统整合代谢全基因组关联分析（mGWAS）和代谢数量性状位点（mQTL）分析，研究鉴定了33 566个显著单核苷酸多态性（SNPs）和3128个mQTL。13个已注释的代谢物在7A染色体上的一个物理区间内共定位。基于代谢物和基于表型的GWAS及QTL整合分析显示，在4A染色体上有一个赤霉素A4（GA4）含量和籽粒圆度重叠区域。对甲基磺酸乙酯（EMS）诱导突变体的表型分析证实了TaSDR在调节GA4含量和籽粒形态中的作用。这些发现为影响关键籽粒相关性状的代谢通路提供了新的见解，并加深了研究人员对调控小麦籽粒代谢物和表型的复杂分子机制的理解。所鉴定的代谢标记和候选基因为旨在提高小麦产量和品质的分子育种计划提供了有价值的目标。
**小麦籽粒代谢物与性状关联的遗传解析——基于多组学整合与功能验证的研究解读**
**1. 研究背景与科学问题**
普通小麦是全球重要的粮食作物，为人类提供约20%的日热量和15%的蛋白质。籽粒形态与品质性状直接决定产量，但遗传机制复杂，且表型鉴定缺乏标准化流程。植物代谢组是连接基因组与表型组的桥梁，能够反映生长发育和生理状态。然而，此前小麦中代谢物与籽粒性状的关联研究多基于单一群体类型，缺乏将代谢组、表型组与全基因组关联分析（GWAS）及数量性状位点（QTL）定位系统整合的比较研究，导致代谢物与籽粒形态间的遗传基础尚不清晰。为此，本研究旨在通过平行开展代谢物和表型的GWAS与QTL分析，在自然群体和重组自交系群体中共同鉴定共定位区间，揭示代谢物与籽粒相关性状的遗传联系，并利用突变体验证关键候选基因。
**2. 研究内容与主要技术方法**
研究人员以349份中国小麦品种构成的关联群体（AP）和138个家系组成的重组自交系（XX-RIL，襄阳市大穗/新麦26）群体为材料，利用液相色谱-串联质谱（LC-MS/MS）测定籽粒代谢组，共获得947种代谢物；同时考察33个籽粒相关性状（12个形态、21个品质）。基因分型使用小麦660K SNP芯片，GWAS采用混合线性模型，QTL定位采用完备区间作图（ICIM），并整合mGWAS、mQTL、pGWAS、pQTL进行共定位分析。对候选基因TaSDR利用甲基磺酸乙酯（EMS）突变体K3716进行功能验证，并开展外源GA4处理试验。
**3. 主要研究结果**
**3.1 籽粒相关性状的变异与相关性**
基于349份AP群体在多年多点的表型分析发现，籽粒形态与面粉品质性状均呈广泛连续变异，变异系数介于1.19%–191.58%。籽粒圆度、粒长、粒宽、粒厚、千粒重变幅分别为0.59–0.86、5.89–7.32 mm、2.98–3.72 mm、2.75–3.45 mm和31.95–56.15 g。品质性状如 setback粘度、损伤淀粉含量、面团稳定时间变幅分别为383.5–1847 cP、3.1%–27.5%和1.65–10.84 min。相关性分析表明，籽粒形态、淀粉糊化特性内部多呈正相关，但淀粉糊化与面团流变学特性之间普遍负相关，揭示了两类性状间的权衡关系。
**3.2 小麦籽粒代谢组特征**
通过LC-MS/MS共鉴定947种代谢物，其中461种获得结构鉴定或推定注释，分为11个生化类别。比较AP和XX-RIL群体发现，AP群体代谢物变异系数范围（11.09%–618.76%）大于XX-RIL群体（9.73%–569.28%），反映其更广泛的遗传背景。461种注释代谢物相关性分析显示，糖代谢物与其他通路多呈负相关或弱相关，而同一生化类别内代谢物强正相关。脂类代谢物构成密切共调控网络，其中亚油酸（mr1054）为节点；同时存在以色氨酸为核心的跨通路共调控网络，其与吲哚类、喹啉生物碱等下游代谢物显著关联。
**3.3 代谢物与籽粒性状的关联**
通过大规模Pearson相关分析，在AP群体中鉴定到2286个代谢物-性状显著相关（|r|>0.3，p<0.05），其中110个代谢物在多环境下与同一性状稳定关联。籽粒圆度与亚油酸、4-羟基苯甲醛及GA4正相关；粒长与多个黄酮类、多酚类代谢物负相关；粒宽、粒厚和千粒重与黄酮类、激素类代谢物负相关。此外，部分代谢物同时关联籽粒品质与形态，如黄酮类mr048/1097/1092与面粉品质正相关而与籽粒性状负相关。这些关联在AP和XX-RIL中均有验证。
**3.4 代谢物mGWAS与mQTL分析及候选基因鉴定**
mGWAS在AP群体中鉴定出33 566个显著SNP，其中A、B亚基因组数量高于D亚基因组，平均表型变异解释率20.4%。XX-RIL群体联合连锁与关联分析获得3128个QTL、2011个显著SNP。整合分析发现167种代谢物共定位于89个基因组区间，其中13种来自不同通路的代谢物显著相关（r>0.5），并一致定位到7A染色体27.19–30.28 Mb区域。该区间内的单倍型Block1与其中9种代谢物显著关联，进一步通过变异重分析确定候选基因TraesCS7A02G057300（编码锌指GRF结构域蛋白，TaZGRF-7A）。单倍型分析显示，TaZGRF-7A-Hap1携带者（占测试材料40.8%）的9种代谢物含量显著低于其他单倍型，提示其负向调控这些代谢物的积累。
**3.5 整合表型GWAS和QTL解析复杂籽粒性状**
pGWAS检测到22 014个与籽粒性状相关的显著SNP，其中84.24%关联籽粒形态，12.73%关联品质。XX-RIL群体共定位到313个pQTL。将代谢物与表型结果整合，发现1792个SNP同时存在于mGWAS和pGWAS中，96个QTL同时存在于mQTL和pQTL中；395个SNP为多效性SNP，29个热点区间。代谢通路与性状间显示出系统关联，例如脂质通路与籽粒颜色和种子面积共定位最强，黄酮类与粒长关联，而面团流变学特性与激素类代谢物共享最多共定位位点。尤其值得注意的是，GA4含量与籽粒圆度在4A染色体621.8–624.7 Mb区间共定位，并鉴定到候选基因TraesCS4A02G343700（短链脱氢酶/还原酶家族蛋白，TaSDR）。TaSDR-Hap1单倍型与较高GA4含量和籽粒圆度相关，在AP中频率为32.9%。
**3.6 EMS突变体验证TaSDR功能**
在EMS诱导的四倍体小麦Kronos突变体K3716中，TaSDR基因携带错义突变（G235A）。与野生型相比，K3716的GA4含量降低3.45%，籽粒圆度降低3.18%，粒宽降低10.94%，粒厚降低9.66%，千粒重降低27.07%，而粒长无显著变化。外源GA4喷施处理可部分恢复K3716的籽粒圆度、粒长、粒宽、粒厚及千粒重。研究表明TaSDR正调控GA4积累，进而影响籽粒形态。
**4. 讨论与结论**
讨论部分指出，代谢物与性状的关联揭示了植物在生长、防御和品质间通过内在调控网络实现平衡。整合mGWAS和mQTL能够跨群体验证遗传位点，并发现TaZGRF-7A可能作为转录因子调节芳香族氨基酸及次级代谢。比较mGWAS和pGWAS进一步将GA4与籽粒圆度关联到TaSDR，结合突变体证据表明TaSDR参与GA4代谢或信号，但其是否直接影响GA4合成、或通过ABA-赤霉素交叉对话发挥作用仍需深入研究。
研究结论表明，这些发现为影响关键籽粒相关性状的代谢途径提供了新见解，加深了调控小麦籽粒代谢物和表型的复杂分子机制的理解。所鉴定的代谢标记和候选基因为旨在改良小麦产量和品质的分子育种提供了有价值的目标。
