基于图像的水稻粒型数量性状位点（QTL）定位用于建立预测育种体系

《Plant Physiology and Biochemistry》：Image-based QTL mapping of grain size for establishment of predictive breeding in rice

　　

水稻是全球大部分人口的主食和热量主要来源。随着全球人口的持续增长，培育高产水稻品种对未来粮食安全至关重要。粒型是与水稻产量直接相关的关键性状。在这项研究中，研究人员利用游标卡尺收集表型数据和基于图像的表型分析进行了QTL定位。通过卡尺测量确定的所有QTL也均通过图像数据检测到，这使得能够以更高的LOD值（Logarithm of Odds，对数几率）实现更精确的定位。在染色体3、5、6和7上鉴定到了粒型相关QTL，包括GS3、qSW5和GW7等主要基因。研究人员在染色体6上的标记RM586和RM1163之间鉴定到一个新的QTL区域，该区域此前未见报道。该区域的渗入对粒长产生积极影响，且与qGL3结合时观察到加性效应。在RM586–RM1163区域内，注释了16个开放阅读框（ORF），基因本体（GO）分析表明它们在籽粒发育过程中调节细胞结构和细胞器功能。其中，OsGSq6从穗分化期到抽穗期的表达量显著增加。在该基因内鉴定到15个单核苷酸多态性（SNP），产生了11种不同的单倍型，其中几种主要在籼稻中发现。OsGSq6编码磷酸酪氨酸磷酸酶激活剂，表明其在籽粒发育中的作用。基于图像的表型分析还能检测品种混合物，有助于提高遗传纯度。该方法为提高水稻育种精度提供了一种有前景的策略。

该研究聚焦于水稻粒型这一决定产量的核心农艺性状，旨在通过结合高精度图像表型技术与传统的QTL定位方法，挖掘新的遗传位点并评估其在分子育种中的应用潜力。论文题为《基于图像的水稻粒型数量性状位点（QTL）定位用于建立预测育种体系》，发表于《Plant Physiology and Biochemistry》。研究背景源于全球人口增长对粮食安全的压力，以及传统育种中依赖人工测量的表型鉴定存在的效率低、误差大等问题。尽管GS3、GW2、qSW5等主要基因已被广泛利用，但仍有大量微效QTL及环境互作机制未被阐明。为此，研究人员构建了青城（Cheongcheong，长粒型 indica）与纳洞（Nagdong，短粒型 japonica）衍生的120个加倍单倍体（DH）群体（CNDH），综合运用游标卡尺与相机成像系统进行表型鉴定，并辅以分子标记连锁作图，以期精准定位调控粒长的遗传因子。

为实现上述目标，研究人员采用了几项关键技术手段。首先，利用NIKON D7200相机采集籽粒RGB图像，并通过ImageJ软件进行阈值分割、颗粒分析与椭圆拟合，提取粒长与粒宽数据。其次，构建了包含788个SSR标记的基因连锁图谱，利用IciMapping v4.2软件进行复合区间作图（CIM），以LOD > 3.0且表型变异率（PVE）> 10%作为QTL检测阈值。此外，还结合了RNA提取与实时荧光定量PCR（qRT-PCR）分析候选基因在穗分化期的表达模式，并利用RiceVarMap v2.0数据库进行了单倍型分析。

研究结果部分详细阐述了以下发现：

1. 1.
   基于游标卡尺与图像法的粒长、粒宽及长宽比调查
   通过对亲本及CNDH群体的表型分析发现，两种方法均显示青城的粒型显著大于纳洞，且CNDH群体表现出超亲分离现象。图像分析法相比传统卡尺测量，不仅数据一致性更高，且在两年间表现稳定。
2. 2.
   CNDH群体遗传图谱构建与QTL定位
   利用143个多态性SSR标记构建了覆盖12条染色体的遗传图谱。QTL定位结果显示，染色体3（RM15448–RM6266）和5（RM5311–RM440）上的QTL可通过两种方法共同检测到，而染色体6（RM586–RM1163）和7（RM8261–RM420）上的QTL仅通过图像表型数据检出。其中，染色体6上的qIGL6与qIRLW6为新增位点。
3. 3.
   QTL聚合对CNDH群体粒型的影响
   基因型组合分析表明，携带qIGL3^Cheongcheong+qIGW5^Cheongcheong+qIGL6^Cheongcheong+qIGW7^Cheongcheong组合的株系（如CNDH91）表现出最长的粒长、最大的千粒重及最高产量。相关性分析显示图像法与卡尺法测得的数据相关系数高达0.95以上。
4. 4.
   定位区域内候选基因的筛选
   在四个稳定的QTL区间内共筛选出52个候选基因。值得注意的是，染色体6的RM586–RM1163区间未发现已知主效基因，被确认为一个新的QTL区域，因此被选为进一步功能分析的重点。
5. 5.
   候选基因的遗传距离与关联分析
   系统发育树将候选基因分为四支，Venn图显示代谢与信号传导相关基因在各染色体区间普遍存在。GO富集分析表明，染色体6区间内的基因主要富集于细胞组分（Cellular Component）类别，如质体（Plastid）与细胞质（Cytoplasm）。
6. 6.
   候选基因表达的qRT-PCR分析
   对染色体6区间内16个候选基因的表达分析显示，Os06g0218500与OsGSq6在青城及长粒型DH系中表达量持续升高。特别是OsGSq6，其表达量从穗分化期至抽穗前5天逐渐增加，暗示其在籽粒发育中的关键作用。
7. 7.
   OsGSq6的序列变异与单倍型分析
   在OsGSq6基因内鉴定到15个SNP，将4649份水稻种质划分为11种单倍型。Hap2、Hap3、Hap5、Hap9和Hap11主要分布于籼稻，而Hap10为粳稻特有，显示了明显的亚种特异性。
8. 8.
   OsGSq6的序列同源性分析
   序列比对与蛋白互作网络分析表明，OsGSq6编码丝氨酸/苏氨酸蛋白磷酸酶2A激活剂，与小麦（Triticum aestivum）、大麦（Hordeum vulgare）等禾本科作物同源基因具有高度保守性，提示其可能通过调节磷酸酶信号通路参与籽粒发育调控。

在讨论与结论部分，研究人员指出，本研究证实了图像表型技术在QTL定位中具有比传统卡尺测量更高的分辨率，能够识别如染色体6上RM586–RM1163这类被常规方法遗漏的新位点。新鉴定的OsGSq6基因与已知的GS3基因存在加性效应，二者聚合可进一步增加粒长。OsGSq6编码的磷酸酪氨酸磷酸酶激活剂（PTPA）可能通过参与油菜素内酯（Brassinosteroid）信号途径中的去磷酸化过程，正向调控籽粒发育。该研究不仅为水稻粒型遗传改良提供了新的基因资源，也为建立基于图像的精准预测育种体系提供了理论依据和技术支撑。