《Theoretical and Applied Genetics》：Dissection of local haplotype diversity at soybean rust loci reveals resistance-associated and context-dependent variation patterns in diverse germplasm
编辑推荐：
为解决大豆锈病(SBR)抗性遗传基础复杂、单基因抗性易被病原克服的问题，研究人员对2815份大豆种质进行了关联与单倍型分析。他们通过连锁不平衡(LD)引导的局部单倍型分析，揭示了染色体13和18上重要SBR位点的等位基因多样性，鉴定出与降低病害严重度相关的单倍型，并发现了染色体13上与大斑花叶病毒(SMV)抗性相关的特异性单倍型。研究结果展示了利用新型表型数据进行局部单倍型分型在识别连锁因果等位基因和育种应用中的潜力，为基于单倍型的抗性育种提供了见解。
大豆是全球重要的油料和蛋白作物，但其生产常常受到病虫害的严重威胁。其中，由专性寄生真菌豆薯层锈菌(Phakopsora pachyrhizi)引起的大豆锈病(Soybean Rust, SBR)是一种毁灭性的叶部病害，在适宜的环境条件下可导致高达50-90%的产量损失。面对气候变化可能加剧病原体毒力和传播的挑战，寻找和利用更持久的遗传抗性资源变得尤为紧迫。传统的抗性通常表现为不产生可见病斑的免疫反应或产生红褐色(Reddish-brown, RB)病斑，而感病反应则表现为黄褐色(Tan, TAN)病斑。然而，病斑颜色本身可能并非抗性的客观衡量标准，因为孢子的数量和产孢水平与病斑颜色不一定相关，且在RB和TAN类型内部都存在抗性程度的连续变异。此外，病原菌高度适应，其种群的选择性清除和遗传漂变可能导致能够克服已部署的单基因Rpp抗性的新菌株出现。因此，需要采用能够捕捉完整表型梯度而非将其简化为二元分类的基因组学方法，来识别与抗性相关的遗传变异模式。
过往研究已通过连锁作图和全基因组关联研究(Genome-wide association studies, GWAS)鉴定了八个主要的Rpp基因和众多与SBR相关的数量性状位点(Quantitative Trait Loci, QTL)。然而，常规的单点GWAS模型通常独立测试每个单核苷酸多态性(Single-Nucleotide Polymorphism, SNP)与表型的关联。对于大豆这种自花授粉物种，由于自交下的重组减少，以及驯化瓶颈和针对优良性状的强烈选择，存在广泛的连锁不平衡(Linkage Disequilibrium, LD)。高LD可能导致非因果SNP与性状出现关联，并且连锁变异的联合表型效应可能被低估或无法解析。局部单倍型分型作为一种GWAS的精细定位方法，可以分析紧密连锁、共遗传的SNP群组，以识别性状相关位点周围有生物学意义的功能变异模式。在SBR相关基因组区域的表征方面已有多次细化，例如Rpp1基因最初被鉴定为单个显性基因，后来被定位到染色体18上。染色体13上一个97.97-111.88 cM的区域(qSBR_Gm13)也被报道与SBR相关，该区域富含核苷酸结合位点-亮氨酸重复序列(Nucleotide-Binding-Site Leucine-Rich Repeat, NBS-LRR)家族的抗性基因，并且与包括大豆花叶病毒(Soybean Mosaic Virus, SMV)抗性位点Rsv1在内的多种病原抗性相关。
为了更精细地解析SBR抗性的遗传结构，一项发表于《Theoretical and Applied Genetics》的研究采用了基于连锁不平衡(LD)的局部单倍型分析方法，对大规模、表型多样的大豆种质资源进行了深入分析。
本研究主要应用了以下几项关键技术方法：首先，研究利用了来自美国农业部(USDA)种质资源库的2815份表型多样的大豆种质，其SBR病害严重度（病斑密度，1-5级）和病斑颜色（RB或TAN）数据来源于先前已发表的温室接种评估项目。其次，基因型数据来源于SoySNP50K芯片，并利用399份重测序的多样大豆种质作为参考面板，使用Beagle 5.4软件进行了基因型填充，获得了高密度的SNP数据集用于后续分析。再者，使用混合线性模型(Mixed Linear Model, MLM)和固定与随机模型循环概率统一模型(Fixed and random model Circulating Probability Unification, FarmCPU)进行了全基因组关联分析(GWAS)，以鉴定与SBR抗性相关的显著SNP位点。最后，也是核心方法，即利用R软件包crosshap对染色体13和18上GWAS鉴定的显著SNP所在基因组区域进行了局部单倍型分析，该方法通过将强连锁不平衡（r2）的SNP聚类成标记组(Marker Groups, MGs)，并根据这些标记组的等位基因组合将个体分配到不同的单倍型组中，从而解析局部的连锁结构和单倍型多样性。
结果
全基因组关联分析(GWAS)鉴定大豆锈病抗性相关位点
利用SoySNP50K数据和FarmCPU模型进行的GWAS在5条染色体上鉴定出8个显著关联位点。最显著的SNP位于染色体6上，位于之前报道的抗锈病QTL内。另一个显著SNP位于染色体18的56,100,116 bp位置(Gm18_56100116)，该区域与已知的Rpp1基因位点重叠。利用填充数据集进行的关联分析揭示了12个位点，包括染色体18上位于Rpp1位点内的最显著SNP (Gm18_56472200)。在染色体13上29,973,410 bp位置也发现了一个显著标记(Gm13_29973410)，该区域位于之前报道的与多种大豆病害相关的QTL富集区域。
连锁不平衡分析与单倍型窗口选择
围绕显著SNP的LD分析用于界定局部单倍型分析的基因组区域。在染色体18上，围绕Gm18_56100116的一个199.53 kb区域被界定用于分析。在染色体13上，围绕Gm13_29973410的一个359.85 kb区域被选定，该区域包含多个单倍型结构和注释的抗病基因。
染色体18上大豆锈病相关Rpp1位点的局部单倍型与标记组组成
在包含所有RB和TAN个体的全群体分析中，局部单倍型分析识别出由7个标记组(MG1-MG7)构成的12种单倍型组合(A-L)。其中，Gm18_MG5与降低的病害严重度相关（表型差异为-0.29）。携带已知Rpp1/Rpp1b等位基因的材料（如PI 587905）被归类到具有Gm18_MG1和Gm18_MG5交替等位基因的单倍型J中。单倍型G由与降低病害严重度相关的连锁标记组(MG1, MG2, MG3, MG5)构成。值得注意的是，Gm18_MG6和Gm18_MG7显示出正表型差异，表明其交替等位基因与增加的病害严重度相关，且主要出现在TAN个体中。专门针对RB群体的分析发现了6个紧密连锁的标记组(MG1-MG6)，其交替等位基因与这些个体中抗性的降低相关。这些标记组中的变异包括Glyma.18G280400（富含亮氨酸重复序列蛋白）的错义、同义和下游基因变异等。而另一个与改善抗性相关的标记组MG7中的变异则涉及Glyma.18G282000（一个ADP-核糖基化因子GTP酶激活蛋白）。
染色体13上大豆锈病相关位点的局部单倍型组成
在染色体13目标区域的全群体分析中，识别出由5个标记组构成的10种单倍型组合。其中MG3与降低的病害严重度相关。专门针对RB群体的分析发现了几个具有负表型差异的标记组，例如MG5与显著降低的病害严重度相关。MG5中的变异被注释为影响多个基因的内含子、上游或3'非翻译区变异，包括Glyma.13G187300（编码一个未表征的保守蛋白家族UPF0565）的错义变异。
从大豆锈病抗性(RB)材料到感病(TAN)材料的标记组转移
为了评估在RB群体中鉴定出的抗性相关标记组在TAN群体中是否具有相似效应，研究将RB单倍型分析得到的所有显著GWAS位点的标记组应用到TAN群体中相应的基因组区域。结果显示，六个位点的标记组在转移到TAN背景时效应减弱或发生逆转。值得注意的是，位于Gm07_43600726（具有自噬基因Glyma.07G261000的变异）的标记组在两种背景下都保持了负效应方向，尽管在TAN中幅度减小。而在Gm13_29973410位点，MG5在RB个体中显示负表型差异，但在TAN个体中则表现出效应逆转（正表型差异）。
染色体13上大豆花叶病毒(SMV-G1)抗性的局部单倍型分析
在相同的染色体13区域内，针对SMV株系G1的抗性进行了局部单倍型分析，识别出5个标记组。其中MG1、MG2和MG3显示出负表型差异，与抗性相关。由25个抗病材料和1个感病材料组成的单倍型D，具有MG1、MG2和MG3的交替等位基因，表现为一个独特的SMV抗性单倍型组合。携带单显性Rsv1抗SMV基因的材料PI 96983也被归入单倍型D。单倍型D中的变异包括多个基因，例如候选Rsv1基因Glyma.13G184900的内含子变异。
研究结论与讨论
本研究通过基于连锁不平衡的局部单倍型分析，在包含2815份表型多样的大豆种质中，精细解析了染色体13和18上重要大豆锈病抗性位点周围的等位基因多样性和单倍型复杂性。与之前主要关注已知抗源的研究不同，本研究在未经过预先选择的大规模群体中揭示了抗性位点周围存在的连续表型变异和复杂的单倍型结构。例如，在Rpp1位点，鉴定出了既包含与降低病害严重度相关（如Gm18_MG5）也包含与增加病害严重度相关（如Gm18_MG6/MG7）的标记组，并且某些单倍型在RB和TAN个体中均有出现，这凸显了表型变异的光谱性质。
跨群体标记组转移分析表明，在RB群体中鉴定出的抗性相关单倍型效应在转移到TAN群体背景时常常是背景依赖的，其效应可能减弱、消失甚至逆转。这强调了遗传背景、上位性互作或其他位点变异在决定最终表型中的重要性。然而，也发现了在两种背景下均保持稳定抗性效应的位点，如染色体7上的位点（涉及自噬基因Glyma.07G261000），这些位点可能蕴含更持久或背景独立性更强的抗性等位基因，对育种具有重要价值。
此外，在染色体13上富含NBS-LRR基因的区域，研究不仅定位了与SBR相关的单倍型，还成功鉴定出了与SMV抗性相关的特异性单倍型（如单倍型D）。这展示了局部单倍型分型在解析共享基因组区域内病原特异性抗性变异方面的能力。尽管该区域与多种病害抗性相关，但本研究结果表明，针对SBR和SMV的抗性由不同的单倍型所贡献，而非完全共享的遗传变异。
总之，这项研究证明了利用高密度填充SNP数据和完整表型梯度进行局部单倍型分析的强大能力。该方法能够超越单个SNP的关联，捕捉紧密连锁变异的累积效应，精细解析抗性位点的遗传结构，识别出有育种价值的候选单倍型以及携带优良等位基因的个体。研究结果为理解大豆病害抗性的遗传架构提供了新见解，并为实施基于单倍型的、更精准高效的抗性育种策略奠定了坚实基础。未来需要对鉴定出的候选因果等位基因和基因进行功能验证，以最终确认其生物学功能并应用于分子设计育种。
