《Ecosphere》：Differential skin microbiomes and disease outcomes in sympatric hosts despite similar infection burdens
编辑推荐：
摘要：两栖动物正面临全球性衰退，其主要原因之一是壶菌病（chytridiomycosis），该疾病由真菌病原体蛙壶菌（Batrachochytrium dendrobatidis，Bd）感染两栖动物皮肤所致。对Bd感染的抵抗力和耐受性可表现出种内和种间变异，这
摘要：两栖动物正面临全球性衰退，其主要原因之一是壶菌病（chytridiomycosis），该疾病由真菌病原体蛙壶菌（Batrachochytrium dendrobatidis，Bd）感染两栖动物皮肤所致。对Bd感染的抵抗力和耐受性可表现出种内和种间变异，这种变异可能部分取决于宿主皮肤微生物组。在加利福尼亚州内华达山脉，濒危的山黄腿蛙（MYL frogs，Rana sierrae 和 Rana muscosa）对壶菌病高度易感，而同域分布的太平洋雨蛙（Pseudacris regilla）则相对不易感，并可能充当Bd的储存宿主。为探究这两种同域宿主间不同的Bd感染动态是否与其皮肤微生物群落差异相关，研究人员（1）比较了23个内华达山脉样点中MYL蛙与P. regilla的Bd感染情况，并（2）表征了宿主皮肤微生物组（或细菌组，bacteriome）的细菌组分。两种蛙类的Bd感染流行率均约为80%，且P. regilla与MYL蛙的平均感染强度无显著差异。相比之下，细菌组组成存在显著差异。P. regilla的皮肤细菌组多样性高于MYL蛙，且仅MYL蛙的细菌组多样性与Bd感染强度相关。此外，P. regilla细菌组中推定抗真菌细菌的丰富度高于MYL蛙。这些同域宿主间细菌组的差异可能源于其生理或生活史的差异。未来研究对于确定P. regilla微生物组更高的多样性和表观稳定性是其固有Bd耐受性的原因还是结果至关重要。
论文解读：同域两栖动物皮肤微生物组差异与壶菌病感染结局的关联机制研究
一、研究背景与科学问题
两栖动物正面临全球范围内前所未有的种群衰退，壶菌病（chytridiomycosis）作为由真菌病原体蛙壶菌（Batrachochytrium dendrobatidis，Bd）引发的皮肤病，已被确认为导致数百种两栖动物种群崩溃和物种灭绝的关键驱动因素。Bd通过水生游动孢子感染两栖动物角化皮肤，尽管其毒力和传播效率极高，但不同宿主物种对Bd感染的响应存在显著差异：部分物种高度易感并发生大规模死亡，而另一些物种则能够抵抗或耐受感染而不表现明显病症。
在加利福尼亚州内华达山脉，濒危的山黄腿蛙（MYL蛙，Rana sierrae 和 Rana muscosa）对壶菌病高度易感，已从其历史分布区的约93%范围内消失，而与其同域分布的太平洋雨蛙（Pseudacris regilla）则种群稳定，能够在携带远超MYL蛙致死阈值的高Bd感染负荷下存活，因此被推测可能充当Bd的储存宿主（reservoir host）。然而，此前关于P. regilla皮肤细菌组及其与Bd感染关系的知识极为有限，同域宿主间Bd感染动态差异是否与皮肤微生物群落差异相关尚不清楚。本研究旨在通过景观尺度的自然变异，系统比较两种同域两栖动物的Bd感染状况和皮肤细菌组特征，探究皮肤微生物组在宿主疾病耐受性中的潜在作用。
二、研究设计与主要技术方法
研究人员于2023年7月至9月在内华达山脉中部和南部的23个样点（包括约塞米蒂国家公园、因约国家森林、塞拉国家森林和红杉-国王峡谷国家公园）采集了P. regilla（n=113）和MYL蛙（n=142）的皮肤拭子样本，样点涵盖单一物种占据和两物种同域共存的多种类型。Bd感染强度通过定量聚合酶链反应（qPCR）检测ITS1基因拷贝数进行量化；皮肤细菌组通过16S rRNA基因V4区域扩增子测序（Illumina MiSeq平台）进行表征，并使用QIIME2流程进行生物信息学分析，包括扩增子序列变异（ASV）鉴定、α多样性（香农多样性、Faith系统发育多样性、均匀度和丰富度）和β多样性（加权和非加权UniFrac距离）计算。推定抗真菌细菌通过将ASV序列与Amphibac数据库以99%序列相似性进行比对鉴定。统计分析采用Wilcoxon秩和检验、Spearman秩相关、线性混合效应模型、置换多元方差分析（PERMANOVA）和ANCOM-BC差异丰度分析等方法。
三、主要研究结果
Bd感染流行率和感染强度在P. regilla和R. sierrae中均较高**
在全部23个样点中，Bd在19个样点被检出。P. regilla的Bd感染流行率为72.6%，MYL蛙为57%，前者显著高于后者（p=0.0098）。然而，在排除四个预先确定的Bd阴性样点后，两种蛙类的Bd感染流行率几乎相同（P. regilla为80.4%，R. sierrae为80.2%）。感染个体的Bd负荷范围为每拭子115至2,021,118个ITS1基因拷贝，两种蛙类间的平均感染强度无显著差异（Welch's t检验，p=0.124）。在同域共存样点，感染P. regilla的平均Bd负荷略高于R. sierrae（p=0.0185）。
P. regilla皮肤细菌组比R. sierrae更多样化，且与Bd负荷无相关性**
P. regilla在所有α多样性指标（丰富度、Faith系统发育多样性和均匀度）上均显著高于MYL蛙（Wilcoxon秩和检验，p<0.001）。P. regilla的皮肤细菌多样性与Bd负荷无任何相关性（Spearman秩相关，所有p>0.5），也不随感染状态变化。相比之下，MYL蛙的细菌多样性与Bd负荷呈显著负相关（所有指标p<0.001），未感染个体的丰富度和均匀度高于感染个体。线性混合效应模型进一步证实，宿主物种与Bd负荷对香农多样性存在交互效应：MYL蛙的香农多样性与Bd负荷负相关，而P. regilla则无此关系。
皮肤细菌组组成在物种和样点间存在差异
加权和非加权UniFrac距离的PERMANOVA分析均显示，P. regilla与MYL蛙的皮肤细菌组组成显著不同（加权UniFrac伪F=75.91，p=0.001；非加权UniFrac伪F=12.45，p=0.001），且细菌组组成也因样点而异。P. regilla的细菌组组成不随Bd感染状态变化，而MYL蛙感染与未感染个体的细菌组组成显著不同。adonis2模型表明，样点和蛙类物种解释了细菌组变异的大部分，Bd负荷的效应较小但存在显著的物种×Bd负荷交互作用。ANCOM-BC分析鉴定出13个属在两种蛙类间差异丰度显著，其中Flavobacterium和未鉴定的Comamonadaceae在MYL蛙中富集，Salinivibrio和Enterococcus等在P. regilla中富集。
P. regilla相比MYL蛙具有更高的推定抗真菌类群丰富度
通过Amphibac数据库比对，在252个两栖动物皮肤拭子样本中共鉴定出214个推定抗真菌ASV。P. regilla皮肤细菌组中推定抗真菌细菌的丰富度显著高于MYL蛙（Wilcoxon检验，p<0.001）。MYL蛙的推定抗真菌细菌丰富度与Bd负荷呈显著但微弱的负相关（p<0.001），而P. regilla中该丰富度在不同感染强度间保持稳定。差异丰度分析显示，P. regilla中Enterobacteriaceae（两个ASV）、Microbacteriaceae、Yersiniaceae以及Aeromonas、Bacillus、Phyllobacterium、Polaromonas和Pseudomonas等属的推定抗真菌ASV丰度更高，而MYL蛙皮肤细菌组中两个未分类至属水平的Comamonadaceae ASV丰度更高。
四、讨论与结论
本研究发现，尽管P. regilla和MYL蛙在Bd感染流行率和感染强度上相似，但两者的疾病结局截然不同：MYL蛙因壶菌病持续衰退，而P. regilla种群保持稳定，这支持了P. regilla作为Bd储存宿主的假说。两种同域蛙类皮肤细菌组的显著差异可能源于其生活史差异——P. regilla变态后转向陆地栖息，而MYL蛙终身高度水生——这影响了它们从不同环境来源获取微生物的机会。MYL蛙细菌组多样性与Bd负荷的负相关表明Bd感染会扰动其皮肤细菌群落，而P. regilla细菌组对Bd感染的稳定性可能源于其固有的宿主耐受性，即Bd感染引起的皮肤病理较轻，从而对细菌组造成较少干扰。P. regilla更高的推定抗真菌细菌丰富度可能为其提供更强的抗Bd能力，但研究人员强调，基于16S序列相似性的推定功能分析仅是假设生成性的，需通过细菌培养和Bd挑战实验验证。研究结论指出，同域物种在相似病原体负荷下表现出不同的疾病结局和细菌组特征，凸显了未来通过实验手段验证这些关系因果机制的必要性。
