水稻粒长与粒宽全基因组关联分析揭示新位点及优异单倍型聚合效应

《Biology》：Genome-Wide Association Study on Grain Length and Grain Width of Rice

本研究通过全基因组关联分析（GWAS）在231份水稻种质中鉴定出调控粒长（GL）和粒宽（GW）的4个稳定数量性状位点（QTL），包括新位点qGL11，并证实qGL3.1（GS3）、qGL3.2（SMG3）和qGW4（OsOFP14）的协同作用。研究通过单倍型分析筛选出具细长粒表型的优异种质（如Newbonnet等），为水稻粒形精准育种提供了关键基因资源和分子标记辅助选择（MAS）路径。

　　

水稻粒长与粒宽全基因组关联研究

引言

水稻（Oryza sativa）粒形是衡量产量和品质的重要性状，通常以粒长、粒宽和粒厚量化。粒形变异不仅影响单株产量，还与稻米加工性能和商品价值紧密相关。此外，粒形差异通过改变淀粉积累和贮藏蛋白分布，影响米饭食味品质。因此，优化粒形兼具增产和提质双重意义。粒形作为多基因控制的复杂数量性状，已报道超过600个相关QTL，近200个基因被克隆验证。然而，多数QTL定位区间仍较宽，且对不同单倍型组合效应的系统性研究尚显不足。为此，本研究利用231份具有广泛遗传多样性的水稻种质，开展粒长和粒宽的全基因组关联分析（GWAS），旨在发现新位点，并重点评估优异单倍型聚合对粒形的调控作用。

材料与方法

研究材料为来自23个国家和地区的231份水稻种质，包括籼稻（XI, indica）和粳稻（GJ, japonica）。田间试验于2024年5月在中国四川省青莲进行，设置低氮（R1, 0 kg ha^?1）和中氮（R2, 60 kg ha^?1）两种处理。成熟期每份种质随机选取5株作为生物学重复，使用万深数字种子成像系统（SC-G）测量约50粒饱满均匀谷粒的粒长和粒宽。表型数据经SPSS（22.0）进行描述性统计和相关性分析。基因型数据来源于3K水稻基因组计划（3K-RGP），约480万个单核苷酸多态性（SNP）经PLINK（1.9）过滤（缺失率<10%，次要等位基因频率MAF >5%）。利用ADMIXTURE（1.3.0）分析群体结构，TASSEL（5.0）估算亲缘关系（K）并进行主成分分析（PCA）。GWAS采用混合线性模型（MLM），显著性阈值设为p < 1×10^?4，200 kb内连续显著SNP视为一个QTL。候选基因区间定义为峰值SNP侧翼100 kb，利用RAP-DB注释排除假基因和转座子，通过单倍型分析（Haploview 4.2）筛选与表型显著关联的基因。

结果

表型特征

两种氮环境下，粒长变幅为6.11–11.81 mm，粒宽为2.11–3.76 mm，变异系数相近，表明粒形性状受氮处理影响较小。偏度和峰度绝对值均低于1，符合正态分布，满足GWAS分析要求。相关性分析显示，粒长与粒宽在R1和R2环境下均呈极显著负相关。

群体结构分析

群体结构最优分组数K=6，ADMIXTURE聚类图显示231份种质可划分为6个祖先群体，遗传成分差异为后续关联分析提供基础。

粒形性状GWAS

共检测到4个稳定QTL，分布于3、4和11号染色体。其中，粒长相关QTL包括qGL3.1（Chr3: 16,729,325–16,832,795 bp，贡献率14.19%）、qGL3.2（Chr3: 16,878,251–16,893,486 bp，贡献率11.32%）和新位点qGL11（Chr11: 1,113,362–1,113,774 bp，贡献率15.07%）。粒宽相关QTL为qGW4（Chr4: 20,599,993–20,620,638 bp，贡献率13.03%）。qGL3.1和qGL3.2分别与已知基因GS3（Os03g0407400）和SMG3（Os03g0410000）共定位，qGW4与OsOFP14（Os04g0415100）共定位。

候选基因与单倍型分析

qGL11区间内筛选出两个候选基因：推测糖基转移酶家族成员LOC_Os11g03160和推测生长素外排载体蛋白LOC_Os11g02950。单倍型表型比较显示，LOC_Os11g02950的Hap3粒长显著高于其他单倍型，提示其可能为qGL11的候选基因。GS3基于功能SNP（C→A突变）分为Hap1和Hap2，Hap1（平均粒长8.62 mm）显著长于Hap2（8.33 mm），且Hap1主要分布于东亚籼稻。SMG3通过5′UTR和编码区变异区分出Hap1和Hap2，Hap1粒长（8.47 mm）显著大于Hap2。OsOFP14单倍型中，Hap3（TAA等位组合）对应较窄粒宽（约2.70 mm），多见于籼稻。

优异单倍型聚合分析

对GS3、SMG3和OsOFP14单倍型组合分析发现，Type 3（GS3^Hap1+ SMG3^Hap1+ OsOFP14^Hap3）种质粒长最大（平均8.98 mm），粒宽较小，呈细长粒表型，以北美籼稻（如Newbonnet、Skybonnet、Lemont）为主；Type 5（GS3^Hap2+ SMG3^Hap1+ OsOFP14^Hap2）粒长最短（7.19 mm）、粒宽最大（3.27 mm），呈短圆粒表型，多为亚洲粳稻（如Kosh、Kongyu 131）。共筛选出8份细长粒和22份短圆粒优异种质。

讨论

本研究鉴定的QTL中，qGL11为新颖粒长调控位点，qGL3.1（GS3）、qGL3.2（SMG3）和qGW4（OsOFP14）则进一步验证了已知基因在不同遗传背景下的功能。GS3通过竞争性结合G蛋白β亚基负调控粒长；SMG3编码泛素结合酶（E2），与DGS1互作参与油菜素内酯（BR）信号通路调控粒发育；OsOFP14与GS9互作调节颖壳细胞分裂方向。单倍型聚合分析表明，优异等位基因组合（如Type 3）可显著优化粒形，但其效应受遗传背景限制，且粒形最终受“库”（颖壳大小）和“源”（灌浆效率）共同制约。未来需结合基因编辑（如CRISPR/Cas9）和多环境表型鉴定，深化候选基因功能验证和等位基因互作机制研究。

结论

通过GWAS分析，本研究明确了4个水稻粒形关键QTL，发现新位点qGL11，并证实GS3、SMG3和OsOFP14单倍型聚合对粒形的协同调控。筛选的细长粒和短圆粒种质资源为水稻粒形精准改良提供了直接可利用的亲本材料和分子标记。