小麦白粉病菌全球种群结构与动态

《Molecular Ecology》：Global Population Structure and Dynamics of Wheat Powdery Mildew

　　

害虫与病原体的全球流动常可追溯到包括农作物迁徙在内的人类迁移活动。小麦白粉病的病原、子囊菌真菌Blumeria graminis f. sp. tritici（Bgt）呈国际分布，是适应最快的植病病原体之一。为研究种群遗传结构、多样性与迁移潜力，研究人员对来自16个国家、6大洲的653株Bgt分离株进行了简化代表性测序（reduced-representation sequencing）。研究人员开展了系统发育推断、聚类、种群分化分析，并基于隔离-迁移（isolation-with-migration）模型估算基因流。分析揭示全球Bgt种群在遗传上分化为两类：其一为源自欧亚大陆的较古老种群，内部遗传多样性高且彼此间存在混合；其二为遗传多样性较低的种群，定殖于前欧洲殖民地。值得注意的是，巴西菌株在系统发育上独立于具有共同英国殖民历史国家的菌株。低基因流表明在前殖民地中外源毒力分离株的自然迁移有限。相比之下，欧洲、中国与俄罗斯间较高的基因流估算值提示，迁移可能促进Bgt对欧亚大陆内新抗源源的适应。

研究背景与立项依据

小麦白粉病由专性活体营养病原真菌Blumeria graminis f. sp. tritici（Bgt，小麦白粉病菌）引起，是全球小麦减产七大病害之一。该病原菌凭借每年有性生殖周期与多循环流行中形成的较大有效种群大小，具备极高的重组与逆转录转座子驱动变异能力，能快速响应杀菌剂与宿主抗性选择。尽管已有研究指出Bgt随人类贸易与殖民活动全球扩散，但此前系统发育重建多聚焦更古老分化事件，且取样规模有限；Sotiropoulos等人2022年虽给出13国172株全基因组轮廓，仍不足以精细分辨近代奠基者瓶颈、迁移方向与有效种群大小变化。研究人员因而以653株跨16国6洲（1986—2019年收集）的分离株为材料，使用简化代表性测序获得的42020个单核苷酸多态性（single nucleotide polymorphism，SNP）标记，整合系统发育、群体遗传与贝叶斯人口学建模，回答“病原如何从起源中心扩散”及“当代结构是否映射人类小麦栽培史”两个核心问题。该文发表于《Molecular Ecology》。

主要技术方法概览

研究人员采用改良双酶切限制性位点关联DNA测序（ddRADseq）构建文库，Illumina NovaSeq 6000双端150 bp测序，比对至Bgt分离株96224更新参考基因组（GCA_900519115.1），FreeBayes与GATK单倍体调用共识SNP。群体框架上，RAxML最大似然系统树以两株1990年以色列分离株作外类群生根；智能PCA（smartpca）结合gap统计量定簇数k；STRUCTURE推断k=2与k=10 ancestries；SplitsTree4 Neighbor-Net呈网状；VariScan算π、Watterson θ_W、Tajima's D、Fu and Li's F*，vcftools算连锁不平衡（linkage disequilibrium，LD）。迁移与人口学方面，TreeMix推主要迁移边，IMa3（隔离-迁移模型）估分裂时间、有效种群大小与双向基因流，Pyabcranger做近似贝叶斯计算随机森林（ABC-RF）在M40/M41/M44三情景中选奠基顺序；样本队列含北美（美/加188）、欧洲（瑞/法/英/德71）、南半球（南非/巴西/澳90），近东与亚洲基干株不参与人口学模拟。

研究结果

3.1 系统发育推断

基于42020 SNP最大似然树（400次迭代）将全球Bgt裂为两支：一支为肥沃新月/亚洲组（土耳其、埃及、以色列、俄罗斯、尼泊尔、中国），各国亚群长枝高支持（100%）；阿达纳与梅内门土耳其局部株成无混合阶梯系。另一支含欧洲（英/法/瑞/德）及巴西、南非、澳大利亚、北美，后四者呈近期梳状拓扑；南非、澳大利亚、美/加自欧洲祖先后相继分出，巴西与该三洲间隔久远长内枝（100%）。美国蒙大拿株脱离其余美/加株独立成支。两株1990年以色列外类群未邻近日以列其余株而近土耳其株，暗示祖型或近今土耳其而非以色列。

3.2 PCA聚类

前三维解释35%、8%、6%变异；gap统计定k=11簇，ANOVA F与卡方极显著。聚类大体映地理：欧陆68%同簇，以色列75%，俄罗斯94%；巴西100%独簇，澳大利亚100%独簇但与蒙大拿株同簇（簇3）；南非98.3%独簇；土耳其核心簇占其75%以上并吸纳俄黑海株，阿达纳80%成私簇邻以色列，梅内门10株入“十字路口簇”与中国50%株、部分欧株混合，显示土/以/埃/欧间线性混合信号。

3.3 STRUCTURE聚类

LnP(D)选k=10，除Fertile Crescent/亚洲六亚群（埃及、中/尼、俄、土、阿达纳、梅内门）外，欧裔分四亚群：欧洲、巴西、南非、澳/美加。巴西与南非极少渗其他国背景；澳/蒙大拿株小比南非背景但主背景仍归美加；中国株裂两亚支——甘青川西支含33%–52%欧洲背景，陕渝黔豫冀东支52%–98%本土背景；以色列株背景排序为尼< />< />< />

3.4 SplitsTree网络

Neighbor-Net显至少四高度分化谱系（南非、巴西、美加澳、欧洲）的长内边；南非KwaZulu-Natal单株i04_bcG07置巴西与欧洲间，或为南半球×欧洲杂交；盒状不兼容分裂提示重组/杂交网态演化。

3.5 国家间遗传差

smartpca算F_ST：欧洲四国间≤0.1，美加0.037，土-以/埃/俄<0.15；欧-土<1.5（矩阵值<0.15）；南半球与前欧陆均>0.25（美-巴0.785，美-南非0.412，巴-南非0.752）。Mantel与MEM分析全量无显著地理-遗传距离相关，仅Fertile Crescent子集微调R²=0.06，说明跨洋殖民地隔离由奠基事件而非距离衰减主导。

3.6 种群多样与人口学

Tajima's D与Fu and Li's F全负；以/土/俄/瑞具强扩张信号（TajD<-1，π>0.0005，LD r²≤0.07）。澳（π=0.000214）、巴（0.000216）、加（0.000292）近中性、低π、高LD（≥0.08）符合新大陆强瓶颈；美（π=0.000295，F=-2.55，r²=0.02）与南非（π=0.000294）低π但深扩张信号，美LD极低示规律重组。MAT1-1-1与MAT1-2-1两交配型等位频率各国近1:1（美/土/俄0.47），二项式检验不偏离0.5，证有性重组频繁。

3.7 迁移与人口学史（TreeMix+IMa3）

OptM定6迁移边；TreeMix Residual正残示以→土/尼、入中国/巴西/南非潜在流。IMa3选PCA/STRUCTURE/Treemix混合对：多数分裂220–360年（俄-中~140年）；欧→美加、欧→南非显著（p<0.01）单向；巴→欧显著反向（p<0.01，与TreeMix入巴向异）；中俄间IMa3检得中→俄流。分裂时间90%HPD宽，但一致支持近数百年串行分裂。

3.8 ABC人口学建模

ABC-RF分类票：M44（欧先裂→北美→南半球）380/500（后验0.763），M41（北美先）93/500，M40（近同时）27/500。M44参数：鬼魂祖群→欧T₃=0.55千年前（KYA）（0.19–1.07），欧→北美T₁=0.31 KYA（0.09–1.19），北美→南半球T₂=0.04 KYA（0.001–0.20）；祖N_e≈6.53k，欧≈10.17k，北美≈10.51k，南半球≈17.07k；奠基严重度S₃=0.294、S₁=0.445、S₂=0.355（越小瓶颈越强）。迁移参数解释力低，证近期串行奠基为主。

讨论与结论翻译浓缩

研究人员在讨论中指出，Bgt与小麦伴生约600万年，但近500年欧殖民带作物与病原至南半球/美洲/澳，全球扩展呈阶跃式。ABC-RF与IMa3一致支持欧裔最早裂、次北美、最晚南半球，绝对时间落在近600年内，与Sotiropoulos等全球近期扩张论合但遗传-地理相关判断异（因取样/标记谱不同）。北欧单流行单元高相似与Jigisha等2025结论合；跨洋殖民地低π小N_e映人类介导奠基子集抽样。肥沃新月早农迁移带感染单粒小麦等Triticeae促Bgt随连续宿主跨欧亚扩散，故欧亚内混合高；跨洋株因宿主断带与遗传漂变孤立。阿达纳/梅内门私簇、中以西支欧背景、以株多源背景皆映区域间历史渗流。结论上，全球Bgt分“欧亚古老高多样混合系”与“前殖民地低多样奠基系”两类；巴西独立于英殖系；欧-中-俄基因流较高助新抗源适应，殖民地外源毒力自然迁入受限；小麦白粉病全球史是人类栽培史与真菌生物学约束（分生孢子脆、短距初侵染）的共同投影。