2021年8月10日，中国科学院分子植物科学卓越创新中心韩斌研究组在Molecular Plant在线发表了题为“Cytoplasmic and nuclear genome variations of rice hybrids and their parents inform trajectory and strategy of hybrid rice breeding”的研究论文。该论文基于群体基因组学分析，系统描绘了杂交稻/不育系细胞器基因组图谱，并整合细胞核与细胞质变异信息，对水稻不育系进行分类，精细描绘了不育系的多态性特征并讨论了不同类型不育系的育种轨迹和育种策略。该研究为水稻杂种优势的形成和利用提供了新见解，为未来的杂交稻设计育种，即为高效创制优良的杂交组合提供了重要的遗传资源。

　　自上世纪70年代以来，随着袁隆平先生和石明松先生相继发现天然败育水稻材料并成功创制三系和两系杂交水稻制种体系，杂种优势被成功用于水稻育种，大幅提高了水稻的产量。三系和两系杂交制种体系分别基于细胞质雄性败育（cytoplasmic male sterility, CMS）系和环境敏感型雄性败育（environment sensitive genic male sterility, EGMS）系创制而成。目前，引起这两种败育类型的分子机制研究得较为清晰(Ding et al., 2012; Fan et al., 2016; Luo et al., 2013; Wang et al., 2006; Yi et al., 2002; Zhou et al., 2012; Zhou et al., 2014)。且细胞质雄性不育系和环境敏感型不育系的核基因组特征以及对杂交育种的贡献也有报道(Huang et al., 2016; Lin et al., 2020; Lv et al., 2020)。但是，细胞器基因组作为细胞质雄性不育的重要参与者，却常在杂种优势研究中被忽视，基于群体水平的杂交稻细胞器基因组多态性分析鲜有报道，且质核互作限制三系不育系核基因组多态性拓展的程度也未见评估。

　　细胞器基因组由于其母系遗传特征，是系统发育和溯源分析的重要工具，可被用于材料精细分类和血统溯源。本研究以细胞器基因组为切入点，首先分析了1495份杂交稻的细胞器基因组，发现杂交稻/不育系拥有五种主要的细胞器类型，且细胞器基因组与不育系核基因组存在关联，分别属于五种生产上广泛使用的不育系：含败育基因WA352c的三系不育系（以野败型不育系为代表），含orf79的三系不育系（以包台型不育系为代表）、含orfH79的三系不育系（以红莲型不育系为代表），含粳型细胞器的两系不育系（如农垦58S）和含籼型细胞器的两系不育系（如株1S）。通过细胞器溯源和基于核基因组的Kinship分析，研究还发现不育系细胞器携同核基因组分化的主要原因是：在不育系育种过程中，不同类型的核心不育系作为败育资源供体和底盘材料被广泛使用，改良产生一系列亲缘关系较近的衍生材料。基于此，可以将水稻不育系分为五种主要的类型（图1）。此外，该研究还利用二三代混合组装的方法，对5个代表性细胞器基因组（包括叶绿体和线粒体基因组）进行组装，鉴定并评估了类群特异性序列变异。

　　后续分析基于不同的不育系类群探讨不育系基因组特征和改良育种的分子印迹。该研究以D2基因为例，揭示了不育系对育性相关基因的固定可带来遗传累赘效应，从而限制育种家对育种偏好基因的使用（图2）。XP-CLR分析和FST鉴定了不同类型不育系中趋同和分化选择的位点，并基于这些位点分析了不同类型不育系与对应恢复系的配组模式。在多数位点上，杂交稻双亲能够实现遗传互补，但仍有少量位点，如与稻瘟病抗性相关的Pi21基因，与抽穗期和产量相关的Hd3a基因以及与分蘖角度和株型相关的TAC1基因，在含粳型细胞器的两系不育系和对应恢复系中都为育种劣势基因型（图3）。此外，还有少量落于选择信号内的基因，其劣势基因型在不育系群体中受到选择，如与结实率相关的PTB1基因，与柱头外露相关的qSTL3基因。该研究推测可能的原因是：参与不育系育种的遗传资源较为狭窄，少量核心种质资源被作为育种底盘材料而广泛使用，其携带的育种劣势基因型便在不育系群体中广泛传播。因此，发展籼粳杂交稻、挖掘野生稻资源，或可成为打破杂交稻育种瓶颈的重要手段。

　　中国科学院分子植物科学卓越创新中心韩斌院士为该论文的通讯作者，分子植物卓越中心顾周琳博士为论文第一作者，上海师范大学黄学辉教授对本研究提供了建设性意见。本项目得到了国家自然科学基金和中国科学院战略性先导科技专项的资助。

　　论文链接：https://www.cell.com/molecular-plant/fulltext/S1674-2052(21)00320-8#relatedArticles

图1. 杂交稻5种类型细胞器以及其与不育系核基因组的关联

图2. 株型相关基因D2与育性相关基因Rf3紧密连锁

图3. 三种类型不育系的XP-CLR分析以及在相关位点上与对应恢复系的配组模式