IF50.5｜蝙蝠白鼻病病原实为两个宿主特化隐存种

蝙蝠白鼻病的元凶并非单一真菌，而是两个宿主特化的隐存种，且北美引入源可追溯至乌克兰。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

白鼻病（white-nose disease）是由真菌 Pseudogymnoascus destructans 引起的蝙蝠疾病，在北美已导致大量蝙蝠死亡，是哺乳动物中由疾病造成的最高死亡事件之一。此前研究普遍认为该病原为单一物种，且在北美为一次克隆引入后扩张。欧洲种群虽被报道具有遗传多样性，但受限于采样范围和生态数据，相关发现未得到重视。宿主因素和环境条件对疾病结局的影响已有较多探讨，但病原自身的遗传变异在疾病动态中的作用长期被忽视。

科学问题

卡在哪里： 此前研究基于 P. destructans 为单一物种的假设，未系统评估病原内部的遗传分化、宿主特异性及地理来源。有限的亚洲分离株虽显示遗传距离较远，但因缺乏生态背景和足够样本，其分类地位和致病性一直悬而未决。

本文要回答： 利用大规模洲际采样和基因组数据，厘清白鼻病病原的物种界限、宿主偏好、种群结构，并追溯北美引入的精确地理来源。

技术路线

作者主要用了：微卫星基因分型、长读长全基因组测序、Pool-seq、系统发育与种群遗传分析、SPASIBA 空间溯源。

作者首先对来自27个国家的5,479株真菌分离株进行18个微卫星位点基因分型，通过系统发育和主成分分析揭示出两个明显分化的进化支（Pd-1和Pd-2）。随后选取代表性菌株进行Nanopore长读长测序和Illumina短读长测序，组装高质量基因组，利用BUSCO基因和全基因组SNP分析量化两分支间的遗传分化。进一步通过Pool-seq扩大样本量，评估全基因组水平的FST和共享SNP比例。同时结合蝙蝠宿主种类和洞穴微环境数据，检验宿主特化。最后利用欧洲Pd-1的精细种群结构，通过DAPC和SPASIBA两种独立方法追溯北美引入源。

图 1　Fig. 1。图1展示了采样地点分布和基于微卫星的遗传结构。a图显示欧亚大陆255个采样点，b图基于1,866个独特多位点基因型构建的系统发育树清晰分为Pd-1和Pd-2两大支，北美分离株（星号标记）全部落在Pd-1支内。c图PCA进一步证实两分支完全分离。该图确立了隐存种的存在及北美来源属于Pd-1。

核心亮点

亮点 1｜两个隐存种同域共存且均致病

基于5,479株分离株的微卫星基因分型，系统发育树和PCA清晰分为两个进化支Pd-1和Pd-2，两者在欧亚大陆广泛同域分布，并在23个地点共同出现。北美分离株全部属于Pd-1。两分支均从具有白鼻病典型病变的蝙蝠上分离，证明二者都能致病。在实验室培养中，两分支的菌落颜色显著不同（t=-11.58, P<0.001），但生长速率无显著差异，且偏好相似的低温高湿环境。

图 2　Fig. 2。图2从基因组层面验证两分支分化。a图显示用于测序的18株分离株地理来源。b图基于664个BUSCO基因的系统发育树以100% bootstrap支持两分支单系，节点上的site concordance factors分别为92%和76%。c图箱线图显示分支间遗传距离（中位数0.14%）远大于分支内（0.022-0.026%）。d图全基因组FST分布中位数0.88，表明强烈分化。e图热图展示大量分支特异的syntenic clusters。

亮点 2｜宿主特化显著：Pd-2偏好Daubenton's bat

对欧洲16种蝙蝠的1,463份拭子分析显示，87%的Pd-1分离株来自Myotis myotis/blythii，而Pd-2主要来自M. daubentonii（39%）。贝叶斯层次模型估计M. daubentonii携带Pd-2的概率接近1，而携带Pd-1的概率极低，表明其几乎只被Pd-2感染。这种宿主偏好不能由蝙蝠分布解释，因为M. daubentonii在采样区广泛存在，且63%的Pd-1阳性洞穴也有该蝙蝠冬眠。

图 3　Fig. 3。图3揭示欧洲Pd-1的精细种群结构和遗传不连续性。a图EEMS分析显示三条主要遗传屏障：巴尔干与欧洲其他地区之间、波兰至斯洛文尼亚一线、法国/伊比利亚与欧洲其他地区之间。b图显示DAPC重分配距离分布，64%的分离株被正确分配到原采样点，平均距离仅52.6公里，而随机化后正确率仅0.55%，表明存在强烈的局部种群结构。

亮点 3｜基因组分化强烈且无基因流

对12株新测序和7株已发表基因组分析，664个BUSCO基因的系统发育树以100% bootstrap支持两分支单系。分支间FST中位数高达0.79（全基因组窗口），共享SNP比例不足2.7%。分支间序列分歧（0.14%）是分支内（0.022-0.026%）的5倍以上。近似贝叶斯计算支持严格隔离模型（无基因流），估计分歧时间约75万代。尽管两分支内均检测到重组和两种交配型，但同域共存下仍无遗传交换，表明存在生殖隔离。

图 4　Fig. 4。图4展示北美引入源的空间溯源结果。a图SPASIBA分析在超过500万平方公里的连续景观上计算33株北美分离株的分配对数似然，颜色从蓝（低）到红（高），最高似然区域集中在乌克兰西部Podillia地区（黄色点）。b图放大显示两个独立方法（DAPC和SPASIBA）均指向同一9公里半径区域，为该地区作为引入源提供了强有力证据。

亮点 4｜基因组结构差异与隐存种界定

基于微共线性网络分析，平均14.4%的syntenic clusters为分支特异，其中1,250个含单拷贝基因。由cluster存在/缺失矩阵重建的系统发育树以最大支持恢复两分支单系。结合互惠单系、分支内重组和分支间无基因流，符合真菌物种界定标准，支持将Pd-1和Pd-2视为两个隐存种。Pd-2目前仅见于欧亚大陆，未在北美发现，其引入可能对当地蝙蝠构成新威胁。

亮点 5｜北美引入源锁定乌克兰Podillia地区

利用欧洲Pd-1的2,261株分离株（225个地点）的微卫星数据，EEMS分析揭示三条主要遗传不连续带。DAPC和SPASIBA两种独立方法均将北美33株分离株的起源指向乌克兰西部Podillia地区，且定位在9公里半径内。该地区拥有世界最长洞穴之一，且自苏联解体后与美国纽约州洞穴探险界有密切交流，与北美首次发现地吻合。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： 微卫星基因分型、长读长全基因组测序、Pool-seq、系统发育与种群遗传分析、SPASIBA 空间溯源
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事真菌病原、种群基因组学、宿主-病原互作或生物入侵溯源研究的科研人员。

带走一句

病原体内部的隐存种和宿主特化可能被严重低估；大规模采样结合种群基因组学能揭示物种界限并精准溯源，为疾病监测和生物安全提供新范式。