全球尺度寄主转换模式揭示立枯丝核菌（Rhizoctonia solani）内对比性特化水平

《Ecology and Evolution》：Global-Scale Patterns of Host Transitions Suggest Contrasting Levels of Specialization Within Rhizoctonia solani

　　

摘要

植物病原体感染多种寄主的能力对粮食安全和生态系统健康具有重大不利影响。立枯丝核菌（Rhizoctonia solani）是一种全球分布的土传真菌病原体，能够侵染多种具有重要经济价值的作物。其跨地理区域和寄主分类群的广泛分布伴随着高度的遗传多样性。历史上，根据菌丝吻合相容性，R. solani 种复合群被划分为14个吻合群（anastomosis groups，AGs）。尽管被广泛视为广谱病原体，R. solani 的寄主范围尚未在全面的跨AG框架下进行研究，以评估病原体多样性和寄主关联的全球模式。研究人员分析了过去70年间收集的4119条R. solani rDNA-ITS序列的全球数据集，包括来自内布拉斯加州的新采样分类群，并使用包含162条R. solani序列的较小的ITS和TEF-1α串联数据集来评估广泛的系统发育模式是否在不同标记间保持一致。系统发育重建表明，大多数AG形成了支持良好的分支，这些分支与已报道的AG命名大体一致。寄主科水平转换推断表明，广谱性AG可能经历更短暂的寄主关联，而专性谱系表现出更强的寄主相关遗传分化，这有助于观察到的遗传结构。祖先寄主状态重建分析进一步表明，AG-3谱系可能起源于一个具有更广泛寄主范围的祖先，随后主要与茄科（Solanaceae）关联，且在采样数据集中未推断出逆转。这些发现表明，寄主关联和寄主转换模式是广谱性-专性连续谱的重要特征，并有助于全球分布真菌病原体的种群结构形成。

论文解读：全球尺度下立枯丝核菌寄主转换模式及其特化水平差异研究

一、研究背景与科学问题

立枯丝核菌（Rhizoctonia solani）作为一类广泛分布的土传植物病原真菌种复合群，其寄主范围涵盖超过200种具有重要经济价值的植物种类。该病原体引起的立枯病、根腐病和纹枯病等病害，在全球范围内对多种农业和园艺作物造成严重影响。传统分类学依据菌丝融合（anastomosis）相容性，将R. solani划分为14个吻合群（anastomosis groups, AGs），包括AG-1至AG-13以及桥接群AG-BI。不同AG表现出差异化的寄主关联特征：例如，AG-4具有广泛的致病性，可侵染大豆、西兰花、番茄、甜瓜和甘蓝等多种作物；AG-1主要与禾本科作物的纹枯病相关；而AG-3则表现出较强的寄主限制性，主要与马铃薯、烟草等茄科（Solanaceae）植物关联。然而，以往研究多聚焦于单个AG与单一寄主间的互作，或少数AG在有限寄主范围内的致病性比较，缺乏在跨AG框架下对寄主范围进化模式与群体遗传结构的全球尺度综合分析。寄主转换（host switching）——即病原体持续感染并传播至新寄主种群的能力——在不同病原系统中的进化后果存在显著差异。因此，R. solani中寄主转换事件究竟是促进广谱性、专性还是两者兼有的混合模式，目前尚不清楚。该研究旨在探究R. solani各AG间寄主与生活方式转换是否随机发生，抑或呈现朝向更广泛广谱性或更强专性的定向模式。

二、主要技术方法

研究人员构建了包含4119条R. solani rDNA-ITS序列的全球数据集（涵盖NCBI GenBank公开数据及38条来自美国内布拉斯加州新采集的分离株序列），以及由162条ITS和TEF-1α串联序列构成的多位点验证数据集。采用最大似然法（Maximum Likelihood, ML）进行系统发育重建；使用BEAST v1.10.4开展马尔可夫跳跃（Markov-jump）分析推断寄主科间转换事件；运用判别主成分分析（Discriminant Analysis of Principal Components, DAPC）评估AG-3和AG-4的群体遗传结构；通过分子方差分析（AMOVA）和成对F_ST统计量化寄主相关的种群分化；并利用中位连接法构建单倍型网络。

三、研究结果

3.1 R. solani的系统发育重建

基于4119条ITS序列构建的全数据集最大似然系统发育树显示，分离株形成九个与报道AG相对应的独立分支。除AG-2和AG-4外，所有AG均形成支持率超过70%的单系群。AG-2被划分为两个支持度较高的独立分支，暗示该群内存在显著遗传分化。AG-4内部存在三个亚分支。值得注意的是，AG-BI聚在AG-2的某一分支内，该分支与AG-3和AG-10系统发育关系更近。此外，少数分离株在系统发育树上的聚类位置与原始报道AG不一致，涉及AG-2、AG-4、AG-5、AG-11和AG-BI。ITS单独与ITS+TEF-1α串联数据集的系统发育拓扑结构高度一致；不同比对修剪策略下主要AG水平拓扑关系稳健。

3.2 AG间的寄主关联

寄主关联分析显示，多数AG寄主范围广泛，符合R. solani的广谱性特征。然而，AG-3表现出对茄科的强烈关联，GenBank来源的ITS数据集中约97%的AG-3序列来自茄科寄主。AG-2中与苋科（Amaranthaceae）关联的亚分支并未与十字花科（Brassicaceae）关联的亚分支紧密聚类，提示AG-2亚群内可能存在寄主相关的特化。

3.3 马尔可夫跳跃分析揭示寄主转换的模式

马尔可夫跳跃推断显示，苋科是寄主转换事件的重要受体。禾本科（Poaceae）是寄主转换的主要供体，向所有其他寄主科均有显著数量的跳跃事件。十字花科和茄科的寄主转换次数相对较低。所有统计显著的寄主转换事件均具有负z分数，表明观察到的寄主科间转换频率低于随机预期，即寄主转换受到显著约束。祖先寄主状态重建显示，AG-3与AG-2的共同祖先（most recent common ancestor, MRCA）具有多样化寄主关联（以茄科为主），随后AG-3谱系向茄科的关联概率显著增加，并在较近谱系中趋于专一。

3.4 AG-3与AG-4的群体结构

DAPC分析表明，AG-3的遗传结构与寄主关联高度对应，K-means聚类提示约三个遗传簇：簇1主要为马铃薯相关基因型，簇2主要为烟草相关基因型，簇3主要为非茄科寄主基因型。大陆尺度上未见明显的地理聚类。AG-4的BIC曲线缺乏明显拐点，DAPC散点图中不同寄主科个体显著重叠，表明其寄主相关遗传结构较弱。

3.5 AG-3与AG-4的群体遗传学参数

AG-3总体单倍型多样性较低（Hd_Total=0.269±0.019），核苷酸多样性中等（π_Total=0.005）。马铃薯相关群体具有74个多态位点，单倍型多样性极低（Hd=0.125±0.016）；烟草相关群体具有60个多态位点，单倍型多样性较高（Hd=0.989±0.006）。马铃薯群体Tajima's D为-2.716（p<0.001），呈现显著负值。AMOVA显示，寄主分组解释了AG-3中91.8%的ITS变异（Φ_ST=0.918, p<0.001），而仅解释AG-4中6.4%的变异（Φ_ST=0.064, p<0.001）。AG-3中马铃薯与烟草相关群体间F_ST=0.702（p=0.001），烟草与非茄科群体间F_ST=0.570（p=0.001）。AG-4的成对F_ST值范围为0.003至0.378，分化较弱且异质性较高。

3.6 AG-3与AG-4的单倍型网络

AG-3的单倍型网络显示两个主要单倍型中心：Hap_2构成“马铃薯簇”的核心，Hap_4构成“烟草簇”的核心，两簇之间无直接连接或共享单倍型。95.59%的单倍型为各寄主群体所独有。AG-4的单倍型网络呈高度互联结构，五个主要单倍型与多个寄主科相关联，缺乏清晰的寄主科特异性聚类。

四、讨论与结论

该研究提供了令人信服的证据，表明ITS系统发育信号与传统AG分类大致一致。AG-3与AG-4的比较分析揭示了广谱性与专性种群间寄主相关遗传结构的对比性模式：AG-3表现出明显的寄主相关遗传分化，马铃薯与烟草相关群体间存在强分化；而AG-4则保留更广泛的寄主关联和更弱的分化。AG-3的祖先寄主状态重建表明其从多寄主祖先分化而来，逐渐转向茄科优势关联，呈现出寄主相关特化的进化轨迹。寄主转换推断表明，R. solani的寄主转换并非随机，而是受到显著约束，禾本科作为主要供体寄主科，苋科作为主要受体寄主科。这些发现表明，寄主关联与寄主转换模式是广谱性-专性连续谱的重要特征，寄主相关选择可能塑造了部分R. solani种群的遗传结构。马铃薯相关AG-3中观察到的低遗传多样性可能与集约化单一种植导致的瓶颈效应有关，但需注意核糖体DNA重复序列可能受致同进化（concerted evolution）影响。该项研究代表了全球尺度上首个跨AG的系统发育与群体结构分析，为理解土传病原真菌的寄主范围进化提供了重要资源，并为制定针对特定R. solani谱系的病害管理策略提供了理论依据。该论文发表在《Ecology and Evolution》。