澳大利亚小麦叶枯病菌（Parastagonospora nodorum）种群对寄主改变的适应

《Physiological and Molecular Plant Pathology》：Adaptation of the Australian Parastagonospora nodorum population to alterations in the host

　　

小麦叶枯病（SNB）是由Parastagonospora nodorum引起的一种经济上重要的小麦真菌病害。该病主要通过培育抗性小麦品种进行防控，但过去50年的经验表明，新的病原种群会迅速进化，对当前流行品种表现出更强毒力。研究人员组装了360株P. nodorum分离菌面板，该集合可划分为八个亚群，包括一个核心种群和七个瞬时种群，它们在空间与时间分布、交配型、效应子单倍型以及名为Molly的Tc-1 mariner转座子的完整与降解拷贝模式上具有对比特征。Molly可在真菌基因组中增殖并随机插入，其在有性种群中的扩增可能触发重复诱导点突变（RIP），这在一定程度上解释了广泛的遗传多样性。Molly活动与有性繁殖及不同选择压力共同作为进化驱动力，塑造了澳大利亚P. nodorum种群结构，并阐释了该重要小麦病原菌新谱系形成的模型。在小麦上测试时，近期涌现的群体在现代精英品种上表现出更强致病性，与此前观察到的低幅涨落周期（LABAB）一致。该研究识别并表征了可能导致病原更好适应与寄主专化的潜在机制，并提出了提高SNB抗性育种效率与持久性的实用措施。

研究背景与问题提出

小麦叶枯病（Septoria nodorum blotch, SNB）由Parastagonospora nodorum引起，是澳大利亚及全球部分麦区的重要叶部与颖壳病害，造成产量与品质损失。生产上依赖抗性品种与农药，但病原反复进化突破抗性，呈现不同于锈病与白粉病经典“boom-and-bust”的“低幅涨落周期（low-amplitude boom-and-bust, LABAB）”。前期基于SSR标记发现西澳种群存在核心与瞬时组，但全基因组层面转座子（transposable element, TE）尤其Molly（Tc-1 mariner类）如何与有性繁殖、RIP（repeat-induced point mutation，重复诱导点突变）共同塑造种群并驱动寄主适应仍不清楚。研究人员因此整合50年跨度360株澳大利亚分离菌全基因组数据，重审种群结构、Molly动态与致病性演化。

主要技术方法概括

研究人员构建三组澳大利亚P. nodorum分离菌队列（历史151株、西澳DPIRD 97株、2019—2021年Northam苗圃112株，合计360株，采样1968—2021年、覆盖西澳38个地点）。采用Illumina短读长测序全体360株并组装，另选3株做Oxford Nanopore长读长端粒到端粒组装；以SN15为参考行GATK流程SNP calling，筛选18,208个SNP用STRUCTURE与DAPC定群；用MetaEuk判交配型（MAT1-1/MAT1-2）与效应子（SnToxA等）单倍型；以Mauve比对核心与瞬时组基因组、用RIPper与OcculterCut分析LRAR（Large RIP Affected Region）与GC分区；以read-coverage与BLAST定Molly拷贝数及完整/降解分类；用poppr/PopGenome算rbarD、重组率、F_ST、π、Tajima’s D；以幼苗接种 assay 测Era I–III小麦品种（Halberd、Eradu、Calingiri、Wyalkatchem、Calibre、Mace、Scepter）响应。

研究结果

种群结构

基于18,208 SNP的PCA、STRUCTURE（K=8）与DAPC一致将360株分为1个核心组（Group 1，n=248，1968—2021，33地点）和7个瞬时组（Group 2—8，n=11—24）。核心组MAT1-1:MAT1-2≈1:1（χ² P=0.703），私有等位基因最高、rbarD最低（0.0012）、重组率高；瞬时组F_ST对核心组0.24—0.57，组内效应子单倍型单调，证实瞬时组源自核心组并受限时空分布。

群体特征

核心组为有性主导、高多样性的储库；Group 3、8每株重组率最高。Tajima’s D滑动窗口揭示Group 3受选区域方向异于其他，SnToxA所在Chromosome_04、SnTox267与交配型区邻近选泽热点。92.8%分离菌携5个效应子（SnToxA/1/3/5/267）活性等位，瞬时组多为单一单倍型。

基因组含量差异

Mauve比对边合成块内散布序列与独有序列，4524个差异位点中约94%在P. nodorum参考基因池中存在，Molly、Alternaria alternata来源及Pyrenophora teres f. teres（Ptt）水平转移片段占少数，证明瞬时组分化伴随局部基因内容获得/丢失而非全基因组重构。

Molly与RIP

3个长读长基因组（15FG38、16FG168、PN315）示Molly全基因组随机插入，71%落入LRAR；Group 1 Molly拷贝0—45（均数完整3.65/降解3.08），Group 5、7拷贝激增（均数总41—45、34—42）且完整拷贝多，Group 6无Molly。Molly序列中91—100% SNP具RIP特征（C:G?T:A），伴随AT富集区，证实Molly爆发→RIP泄漏至非重复区→贡献多样性的链条。

致病性差异

幼苗接种显示：Era I老品种对最老瞬时Group 6最感；Era II（Wyalkatchem）与Era III（Calibre、Mace、Scepter）上，新近Group 7（2019单点）与Group 8（2019—2021）毒力显著更高，符合LABAB——核心组维持变异，新品种大面积推广后短时涌现更适配瞬时谱系。

讨论与结论翻译

讨论指出，SNB在澳大利亚年损估约1.12亿澳元。研究以360株全基因组确认八组结构（一核心七瞬时），与前人SSR/pangenome研究衔接；瞬时组源自核心有性库，偶发涌现系强寄主专化或短暂环境选择所致。Molly（Tc-1 mariner）可随机插入基因/非基因区，其增殖在有性世代触发RIP，RIP从转座子富集区泄漏至编码区附近（48% Molly插入距基因<1 kb，个别<10 bp，16FG168中邻近SnToxA），与负选择平衡非同义突变，共同生成可适应变异。近期Group 7/8在现代品种上更强致病性验证LABAB模型。结论：核心组为有性变异储库，Molly活动+RIP+有性重组+选育压塑造澳大利亚P. nodorum结构；新谱系通过转座子扰动与效应子单倍型固定实现寄主适应；育种应跨Era混合选取核心组与代表性新瞬时组隔离菌评价抗性，以延缓LABAB式毒力回升。