《Plant, Cell & Environment》：Revealing the Key Genetic Basis of Late Leaf Spot Resistance in Peanut by Combining QTL Mapping, QTL-seq and RNA-seq
编辑推荐：
晚斑病（late leaf spot, LLS）是全球花生最具破坏性和分布最广的真菌病害，然而其遗传基础相对未知。本研究利用一张包含1327个重组bin（源自9533个SNP/InDel标记）的遗传图谱，结合多环境表型数据，在中花5号×ICGV 86699重组
晚斑病（late leaf spot, LLS）是全球花生最具破坏性和分布最广的真菌病害，然而其遗传基础相对未知。本研究利用一张包含1327个重组bin（源自9533个SNP/InDel标记）的遗传图谱，结合多环境表型数据，在中花5号×ICGV 86699重组自交系（recombinant inbred line, RIL）群体中重复检测到位于A02和B06染色体上的两个主效LLS抗性QTL。QTL-seq分析证实了qLLSA02，两种作图方法的整合将A02区间缩小至1.0–3.4?Mb区域。转录组分析鉴定出数千个差异表达基因。通过对该目标区域内变异、功能注释、表达模式、等位效应和比较基因组学的系统研究，研究人员筛选出10个LLS抗性候选基因，并进一步优先确定其中两个为最有希望的候选基因。基于这两个优先候选基因的SNP开发了两个竞争性等位基因特异性PCR（kompetitive allele-specific PCR, KASP）标记，显示基因型与LLS表型之间存在强相关性。通路分析表明，黄酮生物合成、植物激素信号转导和植物-病原互作途径参与了LLS感染的响应。本研究揭示了花生LLS抗性的遗传和分子机制，所开发的两个KASP标记可促进LLS抗性花生育种。
研究背景与科学问题
花生是全球重要的油料作物，但其生产受到多种病害的限制，其中晚斑病（late leaf spot, LLS）是最具破坏性和分布最广的真菌病害之一。LLS由病原真菌侵染叶片和茎秆引起，严重时导致落叶甚至植株死亡，通常造成约30%的产量损失，严重时超过50%。培育和种植抗病品种是防控LLS最经济、安全和环保的策略。然而，由于LLS抗性的遗传基础尚不清晰，抗性种质资源有限，目前生产中可利用的抗病品种很少。尽管已有大量研究致力于探索LLS抗性的遗传基础，但多数研究早于四倍体花生参考基因组的发布，其遗传图谱主要基于二倍体祖先种或无参考基因组的方法，无法提供精确的QTL信息。此外，大多数研究依赖单一方法，对抗性机制的研究仍显不足。因此，有必要采用合适的参考基因组并进行多方法联合分析，以深入解析LLS抗性的遗传基础。
研究内容与意义
本研究综合运用传统QTL定位、QTL-seq和RNA-seq技术，系统研究了花生LLS抗性的遗传与分子机制。通过QTL定位和QTL-seq分析，研究人员鉴定了LLS抗性的主效QTL；通过系统分析序列变异、功能注释、表达模式、等位效应和比较基因组学，发现了候选基因和优先候选基因；基于优先候选基因的变异开发了KASP标记；此外，转录组分析揭示了参与LLS感染响应的关键通路。这些发现加深了对LLS抗性的理解，并为花生抗性育种提供了资源和工具。
关键技术方法
- 1.
遗传连锁图谱构建与QTL定位：利用ddRAD测序数据，对中花5号×ICGV 86699 RIL群体（166个家系）进行基因分型，基于9533个高质量SNP/InDel标记构建包含1327个重组bin的高密度遗传图谱，采用复合区间作图法进行QTL定位。
- 2.
QTL-seq分析：从RIL群体中选取极端抗病和感病株系分别构建抗病池（resistant bulk, RB）和感病池（susceptible bulk, SB），连同双亲进行全基因组重测序，通过计算SNP-index和ΔSNP-index鉴定与LLS抗性相关的基因组区间。
- 3.
转录组测序与分析：对来自感病组三个家系在四个不同发病阶段（无、轻度、中度、重度）以及抗病组三个家系无症状样本共15个样本进行RNA-seq，鉴定差异表达基因（differentially expressed genes, DEGs）并进行KEGG富集分析。
- 4.
KASP标记开发：基于优先候选基因内的SNP位点设计竞争性等位基因特异性PCR（KASP）标记，用于区分抗感材料。
- 5.
共线性分析：利用minimap2和syri软件对栽培种花生Tifrunner与野生种A. cardenasii GKP10017的基因组序列进行比对和结构变异检测。
研究结果
3.1 遗传图谱构建与评估：利用938.52百万条测序读段，获得9533个高质量变异位点（8879个SNP和655个InDel），划分为1327个重组bin，构建了覆盖20个连锁群、总长1457.0 cM、平均bin间距1.11 cM的遗传图谱。各连锁群与其对应染色体的物理位置显示出高度共线性。
3.2 传统QTL定位鉴定LLS抗性QTL：共检测到13个分布于A01、A02、A04、A06、A08、B01、B03、B04、B06和B08染色体的LLS抗性QTL，LOD值范围3.40–11.98，表型变异解释率（phenotypic variance explained, PVE）为4.40%–16.90%。其中三个QTL（qLLSA02、qLLSA08-1和qLLSB06-1）在多个环境中被重复检测到，qLLSA02、qLLSA06、qLLSB01和qLLSB06-1为主效QTL。qLLSA02位于A02染色体2,415,578–3,434,671 bp区间，qLLSB06-1位于B06染色体6.78–18.35 Mb区间，抗性有利等位基因均来自ICGV 86699。
3.3 QTL-seq鉴定控制LLS抗性的主效QTL：对双亲和极端抗感池进行全基因组重测序后，通过计算ΔSNP-index，在A02染色体1.0–6.1 Mb区域鉴定出一个超过99%置信区间的显著基因组区间。在该候选区间内，抗病池和感病池的平均SNP-index分别为0.24和0.96，表明该区域存在一个强QTL信号。
3.4 转录组分析：构建并测序了15个RNA-seq文库，共组装出43,563个新转录本。在七个成对比较组中共鉴定出DEGs，数量从2,876到19,322不等。通过定义四个比较组（DEG1–DEG4），共获得10,107个非冗余DEGs。维恩图分析显示515个在所有四个比较组中共有的DEGs，代表了一组既组成型差异表达又响应感染的基因核心集。
3.5 综合候选基因分析与KASP标记开发：QTL定位与QTL-seq的重叠区间（1.0–3.4 Mb）被确定为靶标区域。该区域注释了222个基因和30,779个变异，包括39个R基因。通过表达谱分析，10个在感病材料中显著高表达的R基因被确定为候选基因。其中，arahy.FXRP5B和arahy.G3725Q表现最为突出：它们在抗病材料中表达极低，在感病材料中表达显著升高，且随病情发展呈递增趋势，符合隐性抗性特征，且在DEG1–DEG4所有比较组中均显著差异表达，被定为优先候选基因。两者均预测编码拟南芥中具有LRR和NB-ARC结构域的抗病蛋白同源物。基于这两个基因的SNP成功开发了两个KASP标记（A02-2511086和A02-2553959），能够有效区分抗感和感病RIL家系，基因型与LLS表型间存在强相关性。
3.6 KEGG分析：对10,107个DEGs进行KEGG富集分析，发现最显著富集的通路是“黄酮生物合成”（Q值：7.34E-19），涉及71个DEGs，其中关键酶基因（C4H、CHS、CHI、HCT）显著上调。“植物激素信号转导”通路富集基因数最多（136个DEGs），涉及水杨酸（salicylic acid, SA）、茉莉酸（jasmonic acid, JA）、乙烯（ethylene, ET）、生长素（auxin, IAA）和油菜素内酯（brassinosteroids, BR）介导的防御反应。SA通路组分（TGA、PR-1）上调，JA通路关键调控因子MYC2下调，ET受体ETR和关键信号因子EIN3上调，IAA响应基因（AUX/IAA、SAUR、GH3）普遍下调，BR信号中BAK1和BSK上调。“植物-病原互作”通路中，钙信号相关基因（CNGCs、CMLs、CDPKs、Rbohs）显著上调，暗示钙介导的CDPK-Rboh模块在促进活性氧（reactive oxygen species, ROS）爆发和抵御病原菌中发挥关键作用。
讨论与结论
讨论：本研究通过整合传统QTL定位、QTL-seq和转录组测序，成功鉴定了A02染色体上LLS抗性的基因组区间，筛选并优先确定了候选基因，证实了多方法整合策略在解析性状遗传基础方面的有效性。A02 1.0–3.4 Mb区间富含R基因（约16个/Mb，远高于全基因组平均~0.3个/Mb）。两个优先候选基因arahy.FXRP5B和arahy.G3725Q表现出隐性抗性特征，其表达水平在抗感材料间差异显著，且随病情发展递增，仅在它们身上同时满足组成型和动态感染响应特征。比较基因组分析表明，ICGV 86699的A02片段来源于野生种A. cardenasii，且两个优先候选基因恰好位于一个大的倒位结构变异区域内。KEGG分析揭示黄酮生物合成、SA和ET信号激活、钙信号介导的ROS爆发等通路协同参与抗病反应，而IAA和JA模块则受到抑制。
结论：通过整合QTL定位和QTL-seq，将A02染色体上的LLS抗性位点界定在1.0–3.4 Mb区域。对该目标区域内变异、功能注释、表达模式、等位效应和比较基因组学的系统研究揭示了十个候选基因，其中两个被优先确定。此外，还开发了两个用于LLS抗性的KASP标记。KEGG通路富集分析表明，黄酮生物合成、植物激素信号转导和植物-病原互作是响应LLS感染的关键通路。总之，本研究为理解花生LLS抗性的遗传基础和机制提供了见解。
