斑点钝口螈（Ambystoma maculatum）皮肤微生物组群落动态变化支持对蛙壶菌（Batrachochytrium dendrobatidis, Bd）地方流行性和外来流行性菌株的抵抗作用

《Molecular Ecology》：Changes to Skin Microbiome Community Dynamics Support Resistance to Enzootic and Epizootic Strains of Batrachochytrium dendrobatidis (Bd) in the Spotted Salamander (Ambystoma maculatum)

　　

摘要翻译：动物病原体的全球性传播已导致物种数量下降和灭绝。在特定疾病呈地方流行性的区域，病原体抑制可能通过宿主相关微生物组提供的保护作用而实现。两栖动物的皮肤微生物组能够抑制真菌病原体蛙壶菌（Batrachochytrium dendrobatidis, Bd）的生长，通过一系列微生物介导的抗真菌机制防止疾病发生，使宿主能够抵抗Bd感染。然而，皮肤微生物组在感染不同Bd菌株类型后，其群落组成或结构如何变化仍不清楚。研究人员评估了具有Bd抗性的两栖动物（斑点钝口螈，Ambystoma maculatum）在实验性暴露于Bd-GPL地方流行性和外来流行性菌株后的感染动态，并通过16S rRNA宏条形码技术追踪了从暴露到恢复期间的病原体载量和细菌皮肤微生物组群落响应。研究人员发现，微生物组群落在暴露后发生转变，已知具有Bd抑制功能的微生物的多样性、优势度、丰度和总比例均增加，表明感染期间存在微生物救援效应。研究人员还观察到恢复期间宿主内多样性变异降低，表明宿主种群存在共同的功能响应，并广泛指示微生物群落恢复力。暴露于地方流行性Bd的螈类随时间具有更高的病原体载量，并表现出更持久的群落变化和更多推定具有保护性的微生物组，而外来流行性Bd感染在抑制性微生物短暂增加后被更快清除。总体而言，这些结果表明皮肤微生物组可能为具有Bd抗性的两栖动物提供抵御真菌疾病的关键屏障，病原体暴露诱导的微生物群落结构变化有利于宿主，可能由感染动态的局部共进化变化所驱动。该研究阐明了复杂的宿主-病原体相互作用如何通过皮肤微生物组介导，其机制是微生物群落动态变化有利于抗病原体微生物。

论文解读文章

研究背景与目的：动物病原体的全球性传播对种群、群落和生态系统产生了深远影响，理解疾病从流行性暴发向地方流行性转变过程中的宿主抗性机制至关重要。宿主相关微生物组（包括皮肤和肠道微生物组）是病原体抗性的重要组成部分，在感染过程中，微生物间的相互作用可为宿主提供保护作用，即通过微生物救援效应（microbial rescue）驱动有益抗性功能并限制病原体对宿主的影响。两栖动物皮肤微生物组可通过产生抗真菌次级代谢产物、抗真菌酶以及竞争性排斥等方式抑制蛙壶菌（Batrachochytrium dendrobatidis, Bd）的生长，从而构成对壶菌病的抵抗机制。然而，皮肤微生物组在不同Bd菌株类型感染后其群落组成或结构如何变化，尤其是在具有Bd抗性的宿主中，仍不清楚。已有研究表明，在Bd地方流行性种群中，皮肤微生物组常含有更高比例的Bd抑制性微生物，且群落组成和结构变异较小，但这些动态是否由局部共进化驱动，以及宿主微生物组是否在功能上响应感染，尚缺乏实验证据。此外，对于抗性宿主在感染期间是否发生微生物失调（dysbiosis）也未有深入探究。因此，本研究旨在通过实验暴露于地方流行性和外来流行性Bd菌株，评估斑点钝口螈（Ambystoma maculatum）皮肤微生物组的群落动态变化，验证以下假设：感染后皮肤中具有Bd抑制功能的微生物丰度和优势度会增加，且这种响应具有菌株特异性，对本地菌株暴露产生更强的保护性响应，而外来菌株可能导致更严重的群落扰动和更高的病原体载量。

研究方法与主要技术：研究人员于2019年5月从加拿大安大略省巴克霍恩附近的一个池塘采集了10个斑点钝口螈卵块，在室外中宇宙（mesocosms）中饲养至变态，随后转入陆地环境。实验中使用两个Bd全球泛流行谱系（Bd-GPL）菌株：本地菌株JEL 261（分离自加拿大蒙特利尔）和外来菌株JEL 423（分离自巴拿马El Cope，已知高毒力）。将28只螈随机分配到本地Bd暴露组、外来Bd暴露组和假暴露对照组，每组15只用于微生物组分析（另有18只野生成年个体作为生态背景参考）。通过将1×10⁶个游动孢子直接滴加至螈腹部并浸泡24小时进行感染，在感染前（Day 0）、感染期（Day 5）和恢复期（Day 30）采集皮肤拭子样本。病原体载量通过Taq-man qPCR定量（ITS拷贝数），微生物群落通过16S rRNA基因V4高变区扩增子测序（Illumina MiSeq平台）获得。序列处理使用Qiime2（deblur流程）生成ASV（Amplicon Sequence Variants），并通过抗真菌分离物数据库（Antifungal Isolates Database version 2023.2）以100%序列相似性鉴定推定Bd抑制性ASV。群落分析包括β多样性（加权和未加权Unifrac距离）、群落离散度（BETADISP）、α多样性（Shannon指数和物种丰富度）、优势类群（>1%平均丰度）以及差异丰度分析（DESeq2和ANCOM-BC共识方法）。所有动物实验均获得Trent大学动物护理委员会批准（协议#25635），Bd培养和接种在生物安全许可（#29065）下进行。

研究结果：
1. Bd感染动态：所有暴露的螈在接种后均检测到Bd阳性，但载量较低。排除Day 1后，暴露于本地Bd的螈在所有时间点的病原体载量均显著高于外来Bd（p<0.001）。本地菌株感染持续时间更长，Day 30仍有25%个体阳性；外来菌株感染迅速清除，Day 5时71.4%个体已清除，Day 30仅3.5%阳性。所有螈均无症状，生长速率与对照组无显著差异（本地菌株质量增长略低）。
2. 微生物组群落组成和结构：Bd感染显著改变微生物组β多样性，所有处理组和采样时间点间群落结构均有差异（PERMANOVA p<0.001）。群落离散度在感染期增加，恢复期降低。本地菌株处理中，感染期离散度显著增加，恢复期降至低于感染前水平，表明微生物群落恢复力；外来菌株处理离散度变化不显著，恢复期未显著低于感染前。对照组群落结构随时间变化但离散度无显著变化。
3. α多样性：全微生物组α多样性（Shannon指数）在感染过程中基本稳定，外来菌株在恢复期略有增加（p=0.051）。仅考虑Bd抑制性ASV时，外来菌株感染期和恢复期α多样性均显著增加，且恢复期高于本地菌株；本地菌株无显著变化。对照组全微生物组α多样性稳定，但抑制性ASV多样性随时间下降。
4. 优势类群：感染前两种处理均有单一Bd抑制性ASV（Caulobacter）占优势，感染和恢复期出现额外抑制性微生物优势类群。本地菌株在恢复期优势类群包括Acinetobacter和Comamonadaceae（属未知），而外来菌株在感染期出现Pseudomonas、Acinetobacter和Comamonadaceae的优势，恢复期仅Acinetobacter保持优势。对照组优势类群中Bd抑制性ASV随时间减少。
5. 差异丰度：从感染前到恢复期，本地菌株处理中4个抑制性ASV丰度增加，外来菌株处理中8个抑制性ASV增加，对照组中17个抑制性ASV丰度下降。所有增加的抑制性ASV均属于变形菌门（Proteobacteria）。Pseudomonas和Neorhizobium在两种处理中均增加，但Neorhizobium在对照组也增加。本地菌株特有的增加ASV包括Klebsiella和Comamonas，外来菌株特有的包括Bosea、Stenotrophomonas、Sphingobium和两种Brevundimonas。对照组中两种Pseudomonas和多种抑制性ASV丰度下降。
6. 与野生微生物组比较：野生螈皮肤微生物组中Bd抑制性ASV平均比例11.3%，实验期间检测到的抑制性ASV约54%也存在于野生个体。所有实验室优势Bd抑制性ASV均在野生微生物组中发现，且增加的ASV（如Neorhizobium、Klebsiella、Bosea等）在野生中亦存在。野生个体Bd阳性率56%，但载量极低（平均1.2拷贝）。

讨论与结论：研究表明，具有Bd抗性的斑点钝口螈在实验感染地方流行性和外来流行性Bd菌株后，皮肤微生物组发生了推定具有保护性的群落变化，符合微生物救援效应。暴露于本地菌株的螈产生了更持久和更强的保护性微生物组响应，包括抑制性微生物总比例和优势度增加，以及恢复期群落离散度降低，表明微生物群落恢复力。这可能反映了局部共进化适应，即地方流行性病原体毒力较低，宿主微生物组通过招募有益微生物来调节感染。相比之下，外来菌株虽毒力较高，但感染迅速被清除，可能依赖宿主组成性防御（如抗菌肽）与微生物组变化的协同作用，短期抑制性微生物增加后恢复至接近基线水平，未发生失调。研究还发现，对照组中抑制性微生物随时间下降，而感染组中无抑制性ASV减少，表明感染诱导了特异性响应。这些结果说明皮肤微生物组动态是宿主-病原体相互作用的重要媒介，通过改变微生物群落结构有利于抗病原体微生物，从而增强对真菌疾病的抵抗，且这种机制在地方流行性感染中尤为重要。研究强调，微生物组群落恢复力和抑制性功能增加共同构成了宿主对感染的有机体可塑性组成部分。