毒液相关而非化学感觉基因驱动寄生蜂 Asecodes parviclava 的宿主相关分化：全基因组重测序揭示生理屏障维持特化谱系

《Molecular Ecology》：Host-Associated Differentiation in a Parasitoid Wasp Connected to Venom-Related but Not Chemosensory Genes

寄生蜂（parasitoid wasps）巨大多样性（megadiversity）的驱动机制仍是进化生态学的核心问题。为探究寄生蜂宿主相关分化（host-associated differentiation）背后的基因组结构，研究人员构建了 de novo 基因组组装，并对寄生蜂 Asecodes parviclava（Eulophidae）的17个种群进行了全基因组重测序（whole-genome re-sequencing），这些种群跨越纬度梯度并涉及多个 Galerucella 宿主物种。尽管宿主种群存在 panmixia（随机交配/泛交），研究人员鉴定出三个遗传分化的寄生蜂谱系，其中最深刻的分化将 Galerucella tenella 相关寄生蜂与其他种群分开。位点特异性扫描（locus-specific scans）揭示选择作用高度不对称。研究人员检测到显著分化和拷贝数扩张（copy-number expansion）聚集于毒液相关基因家族（venom-related gene families），这与化学感觉基因（chemosensory genes）几乎完全缺乏分化形成鲜明对比。这种解耦进化（decoupled evolution）表明，宿主寻找可能依赖行为可塑性（behavioural plasticity），而生理兼容性（physiological compatibility）则通过多层免疫抑制（immune-suppression）机制被遗传硬连线（genetically hard-wired）。最后，北部种群中的不对称基因流（asymmetrical gene flow）表明，细尺度栖息地重叠（fine-scale habitat overlap）创造了生态连通性（ecological connectivity），但未侵蚀特化谱系，因为强烈的宿主介导选择清除了适应不良的杂种（maladapted hybrids）。最终，该研究表明，生理屏障而非行为识别可能维持特化寄生蜂谱系，为特化天敌中伴随基因流的物种形成（speciation with gene flow）的基因组基础提供了新见解。

寄生蜂（parasitoid wasp）是进化生物学中巨大多样性的重要类群。其宿主相关分化（host-associated differentiation, HAD）通常被认为由宿主多维生态位、宿主与寄生蜂之间的军备竞赛以及内寄生生活史共同驱动。内寄生蜂幼虫必须在宿主血腔中发育，而宿主主要防御是血细胞包囊（encapsulation）寄生蜂卵并将其窒息杀死，因此寄生蜂需进化出抑制宿主免疫的毒力机制。理论上，宿主相关分化应同时靶向化学感觉基因（chemosensory genes）以识别宿主，以及毒液相关基因（venom-related genes）以抑制宿主免疫；若识别与免疫抑制不匹配，可能在进化上致命。然而，种群水平上评估宿主相关谱系究竟由行为与生理耦合分化维持，还是二者可独立进化，仍十分缺乏。此前表型研究显示，Asecodes parviclava 的寄生率交互依赖 Galerucella 宿主种类和寄生蜂种群地理起源；南部瑞典寄生蜂毒力高于北部，且与 G. tenella 相关的寄生蜂具有不同于 G. calmariensis 或 G. pusilla 相关种群的毒力表型。Galerucella 宿主免疫能力存在差异，顺序为 G. pusilla > G. tenella > G. calmariensis；转录组分析还表明成功寄生涉及宿主免疫相关基因近乎完全沉默。因此，研究人员利用全基因组数据解析这些表型变化的基因组结构，重点检验毒液相关与化学感觉基因家族在宿主相关分化中的作用，并假设两类候选基因的共定位可能通过减少识别与毒力等位基因间重组来稳定宿主利用表型，同时假设分化最强的毒液位点可能功能上映射到宿主早期免疫反应级联的破坏和沉默。

研究人员开展了 A. parviclava 的 de novo 基因组组装和种群全基因组重测序，对来自17个种群、分为四组、覆盖南北纬度和三种 Galerucella 宿主的个体进行选择扫描，最终识别出三个遗传分化谱系。最显著分歧发生在从 G. tenella 羽化的寄生蜂（Gt-wasps）与从 G. calmariensis 或 G. pusilla 羽化的寄生蜂（Gc-wasps）之间，且该宿主相关分化在南部和北部地点一致。研究检测到毒液相关基因家族中显著分化和拷贝数扩张，而化学感觉基因几乎无分化。北部种群还存在从 Gc-wasps 到 Gt-wasps 的单向基因流，但未侵蚀特化谱系。该研究的意义在于提出生理屏障而非行为识别可能维持特化寄生蜂谱系，为伴随基因流的物种形成提供了基因组框架。

主要关键技术方法（≤250字）：研究人员从瑞典采集寄生蜂 Asecodes parviclava，按纬度与宿主分为 North（3个 G. calmariensis、7个 G. tenella 种群）和 South（3个 G. calmariensis、1个 G. pusilla、3个 G. tenella 种群），最终74个个体来自17个种群、4组。采用单倍体雄蜂构建 de novo 基因组：10X linked-read 文库、Illumina HiSeq X 测序，Supernova 组装，purge_dups、ARCS、LINKS 改进，BUSCO 评估，MAKER 注释。种群重测序用 Illumina HiSeq X 2×151 bp；BWA-MEM 比对，GATK HaplotypeCaller 按倍性 call SNP，过滤得3,051,355个双等位SNP。用 PCA、ADMIXTURE、TreeMix 推断结构；10 kb窗口计算 F_ST、π 与 ΔAF，topGO 富集；毒液和化学感觉候选基因通过 eggNOG 注释与同源搜索鉴定。

研究结果按原文小标题概括如下。3.1 基因组评估和注释：研究人员组装出 A. parviclava 基因组，含11,196个 scaffold，N₅₀=302,549 bp，总大小384 Mbp；BUSCO 节肢动物数据库评估显示92.7%核心完整 BUSCO 组，其中91.9%单拷贝直系同源基因、0.8%重复、1.8%片段化、5.5%缺失。MAKER 注释预测22,290个蛋白编码基因位点，共33,884个 mRNA；18,747个基因获得功能注释，11,265个基因通过 BLAST 匹配已知蛋白。3.1.1 遗传多样性和种群结构：种群遗传多样性 π 为0.00152（Gc North）至0.00187（Gc South）。PCA 反映四个假设种群；ADMIXTURE 主要按宿主来源分离。PC1 分离 Gt 与 Gc/Gp，PC2 分离 Gc-wasps 南北，PC3 分离 Gt-wasps 南北。TreeMix 支持单次迁移事件（m=1），解释43%总遗传方差，优于无迁移模型（25%）；m=2 模型拟合下降（26%）。该模型识别出北部 Gc-wasps 到 Gt-wasps 的强单向迁移边，迁移权重约0.5；南部种群无种群间基因流证据，残差协方差可忽略。3.1.2 F_ST 离群 SNP 和候选选择区域鉴定：窗口 F_ST 呈双峰，采用0.25阈值。Gt-wasps 南 vs 北的 F_ST 较低（0.032），其余对比为 Gc-wasps vs Gt-wasps 南部0.059、Gc-wasps 南 vs 北0.064、Gt-wasps vs Gc-wasps 北部0.070；所有对比中82%–96%的 F_ST<0.1。Gt 南 vs 北离群 SNP 少，故从进一步分析中剔除。3.1.2.1 广泛 GO 比较：对每个种群对比 top 1%离群基因集进行 GO 富集，发现宿主利用对比（Gt vs Gc）的富集 GO 谱在南部和北部高度保守，表明宿主相关选择靶向独立于纬度的共享生物过程；Gc 种群内纬度对比（Gc-S vs Gc-N）产生独特功能特征。3.2 功能基因鉴定：毒液相关候选中，eggNOG 广泛搜索鉴定1,303个基因、34,956个 SNP；靶向同源搜索用107个外部参考序列得到190个独特毒液基因位点，全部多态，含20,417个 SNP。化学感觉候选方面，eggNOG 鉴定222个基因、5,723个 SNP；靶向搜索用91个参考序列得到104个独特化学感觉基因位点，全部多态，含6,428个 SNP。F_ST 离群窗口与功能基因交集显示，29个毒液基因分布在9个 scaffold 上呈高分化；仅2个化学感觉基因（gr28b_2 和 gr28b_3，均在 scaffold 5）在种群间分化。高分化毒液基因涉及蛋白酶抑制剂/pPO 调控、蛋白酶和细胞破坏、基因调控和表观遗传、信号和代谢等类别。仅 Gt-wasps 与 Gc/Gp-wasps 分化且不依赖纬度的基因中，除两个外共同位于 scaffold 4 的1 Mbp区段，包括 SRPN19、KLK13、mp1easpe、ING3、blm 和 PPP1CB；其他位于 scaffold 10和42。CLIPA3、ppaf2_2 和 KLK7 在南部 Gc-wasps 中唯一分化，bicra 呈类似趋势。五个基因在 scaffold 5和6上于北部 Gc-wasps 中唯一分化。scaffold 4 的离群基因簇被进一步定位。

讨论部分总结：研究人员认为，A. parviclava 的种群基因组分析揭示三个分化谱系，与先前表型数据一致，且该分化发生在宿主种群无相应结构的情况下。最显著基因组分歧位于 Gt-wasps 与 Gc-wasps 之间，并在南北一致。Gc-wasps 存在额外纬度分化，而 Gt-wasps 在所研究范围内基因组相似。差异可能与 Gc-wasps 宿主利用南北不同有关：南部为 Gc+Gp，北部仅 Gc。南部 Gc-wasp 种群克服宿主防御能力更强，可能源于 G. pusilla 免疫能力更强带来的选择压力；G. calmariensis 几乎不能包囊寄生蜂卵。南部 G. pusilla 的强烈选择压力似乎补偿了在较弱 G. calmariensis 上的较高存活率；来自两种宿主的南部寄生蜂均能攻击 G. pusilla，交叉表型显示二者功能毒力无差异，PCA 中 Gp-wasps 与 Gc-wasps 完全交叉，支持归为单一谱系。两个宿主种常同株共现，可能促进基因流并阻止进一步宿主相关分化。Gt-wasps 具有独特毒力表型。分析识别 scaffold 4–6 的进化热点，多个毒液相关离群基因成簇且多拷贝，在宿主利用对比和纬度对比中不同。结合先前基因表达研究，G. pusilla 攻击后39个经典免疫基因被诱导，而 G. calmariensis 仅1个被诱导；因此寄生蜂毒力可能沉默免疫级联早期基因，如宿主识别或信号基因。候选离群基因包括 ING3 和 atad2，可能通过重塑宿主染色质结构实现免疫相关转录的表观沉默；ythdf 作为 m⁶A 修饰主要阅读器，可能促进宿主免疫转录本选择性衰减；PPP1CB 样基因产物可能通过去磷酸化 MAPK 或 NF-κB 信号通路关键组分提供次级抑制；KLK7、CAPN7 和 chym2 可能直接干扰宿主血细胞。寄生蜂干扰宿主免疫可能复杂，涉及多基因产物，scaffold 4、5和6的复合体是未来功能研究关键靶点。其他基因可能因基因组碎片化被遗漏，基于同源和功能的基因列表需 RNAi 或基因编辑验证。北部种群间存在大量基因流，单向从 Gc-wasps 到 Gt-wasps，这与生态观察一致：G. calmariensis 存在地点 G. tenella 寄生增加；Filipendula ulmaria（Gt 宿主植物）比 Lythrum salicaria（Gc/Gp 宿主植物）更吸引寄生蜂，可能因其高芳香且产大量花蜜。该芳香差异导致寄生蜂切换宿主植物，主要在北部，因为两种植物常在海波罗的海岸边并排生长；南部生态分离，L. salicaria 占湿地，F. ulmaria 生于较干地点。虽然基因流可能干扰 Gt-wasp 毒力，研究人员认为杂种后代成功率高度可变，只有保留关键毒力基因的个体能存活宿主防御。化学感觉基因分化相对稀少，与毒液相关基因广泛分化形成鲜明对比；虽然两个气味受体显著分化，整体化学感觉工具包分歧最小，尽管已知宿主吸引行为差异。该差异可能表明该性状非遗传固定，而由关联学习或孵化时获得线索支配；若宿主寻找是学习而非硬连线，可解释毒液基因生理适应而无化学感觉系统遗传变化。

研究结论部分翻译：总之，该研究表明，Asecodes parviclava 的分化主要由宿主免疫景观与生态连通性之间复杂且不对称的互作驱动。寄生蜂种群间显著的基因组分化，尤其是 Gt-wasps 与 Gc-wasps 之间的深刻分歧，突出了宿主特异性选择在驱动寄生蜂进化中的强大作用。尽管宿主物种在广泛范围内保持 panmixia（随机交配/泛交），寄生蜂经历了强烈的局部选择并产生显著生理适应，主要映射到 scaffold 4、5和6上的基因组热点。这些高度分化区域包含多样且特化的毒液相关基因阵列。数据表明，这些基因可能通过多层机制促进宿主抑制，从宿主免疫位点的表观沉默和靶向 mRNA 衰减到直接信号干扰。此类复杂毒力策略似乎必要，以克服 Galerucella 物种间高度可变且强健的免疫能力。有趣的是，这种生理适应并未反映在寄生蜂的感觉装置中；化学感觉基因显著缺乏分化表明，宿主寻找可能更多依赖可塑性行为性状，如关联学习，而非固定遗传变化。这种“寻找”与“杀死”宿主机制之间的解耦进化表明，生理兼容性而非行为识别是该系统中宿主范围扩展的主要屏障。此外，北部种群中观察到的单向基因流强调，植物群落组成和细尺度栖息地重叠可影响特化寄生蜂的遗传完整性。研究结果明确表明，只要毒力选择足够强以清除适应不良杂种，即使面对显著迁移，维持不同宿主相关谱系仍可能。最终，这些结果为未来功能研究提供了具体基因组框架，以解析这些候选基因在寄生蜂-宿主军备竞赛中的具体作用。通过整合基因组、表型和生态数据，研究人员阐明了单一寄生蜂物种如何成功应对近缘宿主提供的不同防御环境。