TH118FLM 来源黄瓜霜霉病抗性位点精细定位锁定乙烯相关气孔防御候选基因 CsACO、CsCA 和 CsSTPK

《Journal of the Science of Food and Agriculture》：Fine mapping of cucumber downy mildew resistance loci identifies candidate genes linked to ethylene-associated stomatal defence

　　

背景：由 Pseudoperonospora cubensis 引起的霜霉病（downy mildew, DM）是影响黄瓜的最具破坏性病害之一，仍难以通过依赖杀菌剂的防治策略实现可持续治理。尽管已报道若干抗性位点，但许多先前定位的区域仍过宽，难以用于高效分子标记辅助选择（marker-assisted selection, MAS）或可靠的候选基因优先排序。本研究旨在通过整合高代重组群体、QTL-seq（quantitative trait locus sequencing，数量性状位点测序）、亲本–混池一致性分析（parent–bulk concordance analysis）和候选基因评价，精细定位三个 TH118FLM 来源的霜霉病抗性位点，并鉴定与抗性相关的生物学上合理的候选基因。

结果：研究人员从 TH118FLM 与 WMEJ 杂交中构建了两个重组自交系（recombinant inbred line, RIL）群体 dm2.1-F₇ 和 dm2.2-F₇，以及一个近等基因/回交衍生群体 dm5.1-BC₄F₂。三个群体的表型分布均显示感病疾病指数（disease index, DI）类别频率较高，支持涉及多个抗性相关数量性状位点（quantitative trait loci, QTL），而非单个显性抗性因子。QTL-seq 将 dm2.1、dm2.2 和 dm5.1 区域分别缩小至 333、296 和 694 kb。随后对单核苷酸多态性（single nucleotide polymorphism, SNP）和插入/缺失（insertion/deletion, INDEL）进行亲本–混池一致性分析，进一步将这些位点缩小至 61.7、84.7 和 80.0 kb，分别包含 7、10 和 5 个注释基因。通过整合位置、序列、转录和预测结构证据，CsACO、CsCA 和 CsSTPK 分别被优先列为 dm2.1、dm2.2 和 dm5.1 的候选基因。表达和生理分析进一步显示，这些候选基因响应病原侵染和乙烯利（ethephon, ETH）处理，而抗性背景在乙烯利处理后表现出气孔开度（stomatal aperture）降低和气孔导度（stomatal conductance）改变，支持一个可检验的乙烯相关气孔防御假说。

结论：本研究将三个黄瓜霜霉病抗性位点精细定位为紧凑的、可用于标记开发的区间，并鉴定 CsACO、CsCA 和 CsSTPK 为生物学上可信的候选基因。研究结果为分子标记辅助黄瓜抗性育种提供了实用基因组资源，并提供了一个可检验的生物学框架，提示 TH118FLM 来源的抗性可能涉及乙烯相关气孔防御。

研究背景与问题：黄瓜（Cucumis sativus L.）是全球广泛种植的蔬菜作物，但由 Pseudoperonospora cubensis 引起的霜霉病（downy mildew, DM）是露地和保护地生产中最具破坏性的叶部病害之一，常导致黄化、坏死、落叶和严重产量损失。依赖杀菌剂防治面临成本、环境、残留以及病原群体对杀菌剂不敏感进化等问题，因此寄主抗性被认为是更可持续的防治策略。然而，黄瓜对 P. cubensis 的持久抗性难以实现，原因包括病原遗传多样性高、抗性常为数量性状而非单主效基因控制。已有研究报道了多个抗性位点及候选基因，如 CsSGR、CsLRK10L2 和 CsAAP2A，但许多数量性状位点（quantitative trait loci, QTL）区间仍过宽，难以进行标记辅助选择（marker-assisted selection, MAS）和因果基因鉴定。TH118FLM 是稳定抗霜霉病的韩国黄瓜自交系，前期研究定位了 dm2.1、dm2.2 和 dm5.1 三个抗性位点，但区间仍不适合精细候选基因分析。因此，研究人员开展本研究，旨在精细定位这三个位点，整合证据优先筛选候选基因，并检验是否支持乙烯相关气孔防御假说。

研究人员开展的研究、结论与意义：研究人员以 TH118FLM 为抗性亲本、WMEJ 为感病亲本，构建 dm2.1-F₇、dm2.2-F₇ 两个重组自交系（recombinant inbred line, RIL）群体和一个 dm5.1-BC₄F₂ 近等基因/回交衍生群体，结合重复表型评价、QTL-seq、亲本–混池一致性过滤、候选基因表达、乙烯利（ethephon, ETH）响应、气孔行为和结构预测，将 dm2.1、dm2.2 和 dm5.1 分别缩小至 61.7、84.7 和 80.0 kb，并优先选择 CsACO、CsCA 和 CsSTPK 作为候选基因。研究结论认为，TH118FLM 来源抗性可能涉及乙烯相关气孔防御，但该模型仍需功能验证。其意义在于提供了可用于标记辅助育种的紧凑基因组区间，并提出了可检验的生物学框架。

主要关键技术方法：以 TH118FLM×WMEJ 为亲本，发展高代重组自交系（RIL）dm2.1-F₇、dm2.2-F₇ 和近等基因/回交衍生群体 dm5.1-BC₄F₂；这些群体于 2013–2019 年在韩国 Asia Seed Company（Icheon, Republic of Korea）温室通过标记辅助选择和重复表型筛选发展，最终精细定位阶段评估 100、131 和 146 个高代重组个体。采用 1–8 疾病指数（disease index, DI）在两季自然发病条件下评价抗性，从表型极端个体构建抗性和感病混池，进行全基因组重测序、QTL-seq 和 Δ(SNP-index) 分析；随后实施亲本–混池 SNP/INDEL 一致性过滤。用 RT-qPCR、乙烯利处理、气孔孔径与气孔导度测量及 SWISS-MODEL 结构预测进行候选基因评价。

研究结果：
1. 霜霉病抗性的表型变异与年际一致性（Phenotypic variation and interannual consistency of downy mildew resistance）：通过两个季节的自然发病表型评价，三个群体 DI 值范围为 1–8，感病 DI 类别频率较高，抗性个体较少；抗性混池与感病混池表型明显分离。2019 与 2020 年季节平均 DI 的 Spearman 秩相关分别为 ρ=0.944、0.924 和 0.971，均 P<0.05，说明基因型相对排名高度一致。结论是 TH118FLM 来源抗性为数量性状，可能由多个抗性相关位点控制，而非单个显性因子。

2. 抗性与感病混池的全基因组重测序（Whole-genome resequencing of resistant and susceptible bulks）：通过对抗性混池和感病混池进行全基因组重测序，获得约 90 百万条过滤 reads，比对到黄瓜参考基因组的覆盖率为 99.22%–99.93%。结论是测序数据质量足以支持 QTL-seq 和下游变异分析；相关数据已存入 NCBI Sequence Read Archive。

3. QTL-seq 显著缩小原始抗性区间（QTL-seq substantially reduced the original resistance intervals）：通过基于混池分离的 QTL-seq 和 Δ(SNP-index) 分析，将原始 dm2.1、dm2.2 和 dm5.1 区间分别缩小至 333、296 和 694 kb；基于 Chinese Long 9930 v3.0 注释，这些区间分别含 31、35 和 55 个注释基因。结论是 QTL-seq 大幅降低候选空间，但区间仍不足以直接确定候选基因。

4. 亲本–混池一致性挖掘将每个位点精炼至近基因分辨率（Parent–bulk concordance mining refined each locus to near-gene resolution）：通过保留在双亲间固定且对比、并在抗性混池和感病混池中相互固定的 SNP 与 INDEL，dm2.1、dm2.2 和 dm5.1 进一步缩小至 61.7、84.7 和 80.0 kb，分别含 7、10 和 5 个推定基因。结论是该方法将宽 QTL 转化为紧凑、可用于标记开发的区间。

5. 整合证据优先选择三个精细位点内的候选基因（Integrated evidence prioritised candidate genes within the three refined loci）：通过整合位置、序列、转录、ETH 响应和结构证据，dm2.1 中 CsaV3_2G012920（CsACO）因 I91V 变异、接种后抗性背景表达较高而被优先选择；dm2.2 中 CsaV3_2G030320（CsCA）因 4 个替换、接种后抗性背景表达较高而被优先选择；dm5.1 中 CsaV3_5G015170（CsSTPK）因 C653R 变异、ETH 响应和最强结构差异而被优先选择。结论是 CsACO、CsCA 和 CsSTPK 分别是三个位点的领先候选基因。

6. 候选基因响应乙烯利处理并与气孔行为改变相关（Candidate genes responded to ethephon treatment and were associated with altered stomatal behaviour）：通过 ETH 处理和表达分析，三个候选基因在抗性背景中于处理后 6–24 h 表达升高；ETH 处理使抗性基因型 TH118FLM 气孔孔径显著小于模拟处理，并使气孔导度发生改变，包括处理后 21 h 导度显著降低。结论是 ETH 响应信号可能与抗性背景气孔行为改变相关，但未建立乙烯信号、气孔调控与抗病性之间的直接因果关系。

7. 计算机模拟结构比较突出 CsSTPK（In silico structural comparison highlighted CsSTPK）：通过 SWISS-MODEL 预测蛋白结构，CsACO 的 I91V 未造成整体结构大变化；CsCA 的 4 个替换后结构整体相似；CsSTPK 的 C653R 位于预测激酶区域，抗性等位基因模型预测额外 β-折叠及 C 端构象改变。结论是 CsSTPK 显示最显著的预测结构差异。

8. 由整合证据推导的假设工作模型（Hypothetical working model derived from integrated evidence）：基于位置、转录、结构和生理证据，作者提出工作模型，认为 TH118FLM 来源的 dm2.1、dm2.2 和 dm5.1 等位基因可能贡献乙烯相关气孔防御；抗性基因型中乙烯积累和下游防御激活可能促进气孔关闭并限制病原侵染，而感病基因型响应较弱、气孔开放增加。该模型被定位为未来功能验证框架，而非因果通路证明。

讨论部分总结：讨论认为，利用高代重组材料显著提高了定位分辨率，将三个位点缩小至小于 100 kb 的紧凑区间；但这种群体主要经过目标位点标记辅助选择和重复表型筛选，并非全基因组背景选择，因此不能正式检验位点间加性或聚合效应，残余连锁背景也不能完全排除，标记相关效应大小仍未确定。亲本–混池一致性挖掘被总结为 QTL-seq 后的有效精炼策略，可将仍含数十个基因的区间压缩到近基因分辨率，适用于重组率中等、参考资源不完整或表型成本高的作物。关于候选基因，CsACO 编码 1-氨基环丙烷-1-羧酸氧化酶（1-aminocyclopropane-1-carboxylate oxidase, ACO），催化乙烯生物合成最终限速步骤；CsCA 编码推定碳酸酐酶（carbonic anhydrase, CA），可能参与 CO₂/HCO₃^? 转换、保卫细胞信号、氧化还原平衡和免疫；CsSTPK 编码预测丝氨酸/苏氨酸蛋白激酶（serine/threonine protein kinase, STPK），属于免疫信号网络核心组分，其抗性背景中表达较低可能反映等位特异催化活性或信号阈值，而非简单转录上调。讨论还指出，ETH 处理诱导三个候选基因并伴随气孔行为改变，气孔作为动态免疫界面，可能影响 P. cubensis 这种专性活体营养卵菌的侵染效率、产孢环境或病斑扩展；但这些观察不能证明内源乙烯直接介导抗性，也不能证明气孔关闭是抗性因果基础。育种意义上，持久抗性可能需要聚合机制部分独立的位点；本研究精细区间可用于设计育种友好标记并降低连锁累赘，但候选基因仍需互补、靶向突变、启动子比较、等位特异表达、瞬时功能分析以及乙烯生物合成/信号抑制剂和基因编辑验证。

研究结论翻译：本研究将三个 TH118FLM 来源的黄瓜霜霉病抗性位点 dm2.1、dm2.2 和 dm5.1 精细定位到分别 61.7、84.7 和 80.0 kb 的紧凑区间。通过整合位置、序列、转录、ETH 响应、生理和计算机模拟结构证据，CsACO、CsCA 和 CsSTPK 被优先作为领先候选基因。结果支持一个可检验的工作假说，即 TH118FLM 来源的抗性可能涉及乙烯相关气孔防御，尽管需要功能验证来确定因果基因–机制关系。总体而言，这些精细区间和候选基因为分子标记辅助黄瓜抗性育种提供了实用基因组资源，并为未来功能研究提供了聚焦框架。