《Crop Science》：Unraveling genetic diversity and loci conferring leaf blast resistance in East and Southern African rice germplasm
编辑推荐：
为解决非洲粮食安全受水稻叶瘟病严重威胁的问题，研究人员对来自东非和南部非洲(ESA)的261份水稻种质进行了遗传多样性与全基因组关联分析(GWAS)。研究揭示了ESA水稻独特的遗传结构，并鉴定出多个与叶瘟病抗性显著相关的标记-性状关联(MTAs)，特别是染色体4上贡献20.6%表型变异的主要数量性状位点qBL4及其候选基因LOC_Os04g51030。该研究为培育持久、广谱抗稻瘟病的水稻品种提供了宝贵的基因资源和分子育种靶点。
水稻(Oryza sativa L.)是全球超过一半人口的主要粮食作物。然而，在非洲，水稻生产面临着严峻的挑战，特别是由稻瘟病菌(Pyricularia oryzae)引起的稻瘟病，它能在许多非洲国家造成高达48%的产量损失，严重威胁着地区的粮食安全和农民生计。长久以来，对抗这种病害主要依赖于化学防治和利用抗性基因，但病原菌的快速进化常常导致抗性“失效”，抗性品种很快失去保护作用。更为关键的是，目前水稻育种中使用的大多数抗性基因资源来源于非非洲地区的种质，它们可能不完全适应当地的环境条件，存在所谓的“适应性差距”。为了填补这一空白，挖掘和利用本地适应种质中蕴藏的宝贵抗性基因，就成为培育持久、广谱抗病品种的关键。那么，在东非和南部非洲地区广泛栽培的水稻地方品种和改良品种中，蕴藏着怎样的遗传多样性？它们与全球性的水稻遗传资源相比有何特点？更重要的是，能否从中找到控制稻瘟病抗性的关键遗传位点，为育种家提供“本土化”的武器？发表在《Crop Science》上的这项研究，正是为了回答这些问题。
为了揭示上述科学问题，研究人员设计并实施了一项系统性的研究。他们首先构建了一个由261份来自东非和南部非洲(ESA)不同国家的水稻种质组成的核心群体(ESA panel)，并将其与包含2987份种质的全球水稻3K基因组计划(3K panel)数据集进行整合比较，以评估ESA种质的遗传多样性和在全球水稻基因组中的位置。在布隆迪、肯尼亚、乌干达和科特迪瓦的非洲水稻中心(AfricaRice)Mbe研究站等多个病害高发点，对261份种质进行了叶瘟病(BL)抗性的自然发病评价，采用国际水稻研究所(IRRI)的0-9级病情严重度评分标准，并利用线性混合模型计算最佳线性无偏估计值(BLUEs)用于后续分析。在基因型分析方面，研究采用了基因分型测序技术(Genotyping-by-sequencing)对ESA种质进行基因分型，获得了26,383个高质量的单核苷酸多态性(SNP)标记。通过与3K panel共享的479个SNP，研究人员将ESA种质投射到3K panel定义的主成分分析(PCA)空间中，以量化其遗传关系。最后，利用全基因组关联分析(GWAS)方法，结合单座位模型SUPER和多座位模型BLINK，在控制了群体结构(使用前5个主成分PCs作为协变量)后，系统地扫描了与叶瘟病抗性显著相关的标记-性状关联(MTAs)，并对显著区域进行了连锁不平衡(LD)分析和候选基因(CG)挖掘。
3.1 Variation of BL severity among ESA accessions
对叶瘟病病情严重度的表型分析表明，261份种质的平均病情严重度(BLUE)为4.94，表现出中度的感病水平。三次重复评分间的相关性显著且强烈，表明评价结果稳定可靠。更重要的是，该性状的广义遗传力(H2)估计值高达0.84，说明叶瘟病抗性主要受遗传因素控制，适合进行遗传解析。
3.2 Clustering patterns of ESA accessions in relation to the broader set of 3K varieties
通过系统发育树和主成分分析，研究发现ESA种质在遗传上形成了一个独特的簇，与全球3K panel中的大多数种质(主要来自亚洲)有明显分化。然而，它们与3K panel中的aus(奥萨)亚群以及部分indica(籼稻)亚群(特别是来自非洲的品种)关系最为密切。遗传距离计算进一步确认，ESA种质与aus亚群的遗传距离最近，而与温带粳稻(temp)、热带粳稻(trop)等亚群距离最远。这表明ESA水稻具有独特的遗传背景，其遗传多样性可能未被全球性的资源库充分代表。
3.3 Population structure and genetic diversity within ESA panel
对ESA群体内部的遗传结构分析显示，261份种质可划分为5个遗传簇(Cluster I-V)。这些簇的划分与地理来源和品种类型有一定关联，例如Cluster I主要由来自坦桑尼亚大陆和桑给巴尔的“Supa”系列品种组成，而Cluster V则包含了大量来自非洲水稻中心(AfricaRice)的品种，包括一些非洲栽培稻(O. glaberrima)及其衍生材料。不同簇间的平均叶瘟病严重度存在显著差异。特别值得注意的是，核苷酸多样性(π)分析表明，Cluster V(以AfricaRice材料为主)在整个基因组上表现出最高的遗传多样性(π ≥ 0.3)，而Cluster I的多样性最低。这提示AfricaRice的育种材料为ESA水稻基因库注入了丰富的遗传变异。
3.4 Significant MTAs for blast resistance within ESA panel
利用SUPER和BLINK两种GWAS模型，研究在染色体1、4、6、7、11和12上鉴定出多个与叶瘟病抗性显著相关的MTAs。其中，由两个模型共同检测到、位于染色体4上的标记rs50247和rs16398最为显著，它们分别解释了20.62%的表型变异(pG)，其有利等位基因能使病情严重度降低约1.51个单位(评分0-9分)。
3.5 Linkage block and major QTL on chromosome 4
对染色体4上显著MTAs的连锁不平衡分析揭示了一个紧密的连锁区块(r2≥ 0.8)，该区块定义了一个主要的数量性状位点(QTL)，命名为qBL4。qBL4位于染色体4上30.06 Mb至30.23 Mb的基因组区间内。在该区间内，共挖掘到29个候选基因(CG)，其中包含LOC_Os04g51030，该基因编码一个细胞壁关联激酶(Wall-associated kinase)。根据与LOC_Os04g51030连锁的MTAs的等位基因组合，种质被划分为三个单倍型组(A、B、C)。携带有利等位基因纯合子的A组(63份种质)平均病情严重度显著低于携带不利等位基因纯合子的C组(187份种质)，证明了该位点等位基因的抗病效应。
3.6 Transcription analysis of LOC_Os04g51030, a major CG
为探究候选基因LOC_Os04g51030的功能，研究人员分析了已发表的转录组数据。数据显示，在接种稻瘟病菌后，抗病品种Hui1586中LOC_Os04g51030的表达量在12小时显著高于感病品种日本晴(Nipponbare)，上调了约2.3倍。这进一步支持了LOC_Os04g51030在调控水稻对稻瘟病抗性中可能发挥重要作用。
4 DISCUSSION
研究的结论与讨论部分深刻阐述了本项工作的意义。首先，本研究首次将东非和南部非洲(ESA)的代表性水稻种质系统地整合到全球水稻3K基因组框架中进行评估，揭示了ESA种质独特的遗传结构，它们与全球资源库(尤其以亚洲种质为主的3K panel)存在显著遗传分化，但与aus亚群关系密切。这种分化可能源于历史上的引入路径、长期的本土化适应，以及非洲水稻中心育种项目中非洲栽培稻(O. glaberrima)基因的渗入。这凸显了当前全球种质库中ESA遗传资源代表性不足的问题，强调了加强该地区种质收集与评价的重要性。其次，研究在ESA种质中鉴定出了丰富的遗传多样性，特别是来源于AfricaRice育种项目的种质(Cluster V)展现了最高的核苷酸多样性，是未来育种可资利用的宝贵遗传库。最为关键的是，通过GWAS分析，研究成功定位到一个主效抗叶瘟病QTL——qBL4，并精确定位到候选基因LOC_Os04g51030。单倍型分析和已发表的表达数据共同支持了该基因在抗病中的潜在功能。这一发现为分子标记辅助选择(MAS)提供了直接可用的分子靶标。携带qBL4位点有利单倍型的本地适应种质，可以作为培育广谱、持久抗稻瘟病水稻品种的理想供体亲本，从而有可能克服当前依赖外源抗性基因所带来的适应性差和抗性易丧失的困境。总之，这项研究不仅增进了我们对非洲水稻遗传资源的科学认识，更重要的是为面向非洲水稻生产的、可持续的稻瘟病抗性育种提供了坚实的遗传学基础和具体的育种策略，对保障非洲粮食安全具有重要的实践意义。
