《Molecular Ecology》：Host and Environmental Factors Outweigh Pathogen in Shaping Skin Microbiome of a Threatened Andean Amphibian Community
编辑推荐：
专业摘要翻译：真菌病原体蛙壶菌（Batrachochytrium dendrobatidis，Bd）已导致新热带地区两栖动物严重衰退，然而在流行后地方性病原体条件下，部分物种仍持续存在。研究人员调查了皮肤细菌群落与宿主黏蛋白组（mucosome）活性是否有助于
专业摘要翻译：真菌病原体蛙壶菌（Batrachochytrium dendrobatidis，Bd）已导致新热带地区两栖动物严重衰退，然而在流行后地方性病原体条件下，部分物种仍持续存在。研究人员调查了皮肤细菌群落与宿主黏蛋白组（mucosome）活性是否有助于一个曾为生物多样性极高的乔科-安第斯（Chocó-Andean）群落中七个两栖动物物种的持续存在，该群落曾经历与Bd相关的剧烈衰退。研究人员在旱季和雨季利用16S rRNA宏条形码技术（16S rRNA metabarcoding）表征皮肤细菌群落，并通过体外挑战实验（in vitro challenge assays）评估黏蛋白组介导的对当地Bd谱系的抑制作用。宿主物种身份是细菌群落组成的最强预测因子，而季节与繁殖模式（与栖息地利用密切相关）亦对微生物组结构有所贡献。宿主系统发育与大尺度细菌群落模式相关，而Bd感染状态与感染强度所解释的变异甚微。推测具有Bd抑制能力的细菌在所有宿主物种中均有检出，并在旱季及具有水生幼体的物种中相对丰度更高，但其丰富度与丰度均与Bd感染或载量无关。黏蛋白组实验显示多数物种对Bd生长几乎没有抑制作用；仅Pristimantis laticlavius（历史Bd流行后丰度有所增加的少数物种之一）表现出轻微抑制。总体而言，宿主身份、季节性与栖息地相关的繁殖模式，而非当前Bd感染状态，是两栖动物皮肤细菌群落组成的主要驱动因素。推测Bd抑制细菌的广泛存在及其在某些持续存在宿主中的富集，与潜在的微生物组介导反应一致，但有限的黏蛋白组活性表明，仅凭微生物组难以解释宿主持续存在。整合宿主免疫性状、病原体特征与环境条件，对于理解这一受威胁两栖动物群落中宿主-病原体的长期共存至关重要。
论文解读文章
一、研究背景与科学问题
蛙壶菌（Batrachochytrium dendrobatidis，Bd）是引发壶菌病（chytridiomycosis）的两栖动物皮肤致病真菌，已被视为全球两栖动物衰退的主要驱动因素之一。在新热带山地生态系统中，Bd曾导致大量两栖动物物种的种群崩溃甚至局部灭绝，即使在原始保护区内也不例外。然而，部分物种在经历初始流行（epizootic）后仍能持续存在，并在病原体转为地方性流行（enzootic）的条件下维持种群。这种差异性的疾病结局提示，宿主、病原体与环境之间的复杂互作可能决定物种的存续命运。
传统疾病三角模型（disease triangle）强调宿主、病原体与环境三者的相互作用，但宿主相关微生物组在疾病动态中的角色日益受到重视，由此催生了疾病金字塔（disease pyramid）框架，将宿主和环境相关微生物组作为第四个维度纳入疾病结局的调控网络。两栖动物皮肤表面覆盖着一层由宿主分泌物（如黏蛋白、抗菌肽）、常驻微生物及其代谢产物组成的黏蛋白组（mucosome），构成宿主先天免疫的重要界面。皮肤细菌群落中已鉴定出多种具有抗Bd活性的分类群，这些细菌可能通过产生抑真菌代谢物、竞争资源或调节宿主免疫等方式介导对病原体的抵抗。然而，在自然持续存在种群中，皮肤微生物组是否真正贡献于宿主防御，以及宿主身份、系统发育、环境因子与Bd感染状态如何共同塑造微生物组结构，仍是悬而未决的关键问题。
本研究以哥伦比亚拉普拉纳达自然保护区（Reserva Natural La Planada）为模型系统，该地区自20世纪80年代以来经历了由Bd驱动的剧烈两栖动物衰退，物种多样性下降约67%。研究人员选取七个在衰退后仍持续存在的蛙类物种，涵盖不同的繁殖模式、微栖息地利用和Bd感染动态，旨在检验以下预测：（1）低Bd感染流行率和载量的物种具有更高的细菌丰富度；（2）低感染物种富集更多Bd抑制细菌（符合适应性微生物组假说）；（3）微生物群落组成存在系统发育信号；（4）细菌丰富度存在季节性差异；（5）宿主黏蛋白组具有差异性的抗真菌活性。
二、研究设计与方法概述
研究人员于2022年旱季（8-9月，n=50）和2023年雨季（5-7月，n=62）在拉普拉纳达自然保护区（海拔1700-2100米，哥伦比亚纳里尼奥省）采集了七个物种共119只成年蛙的皮肤拭子样本，另于2025年7月采集了三个物种（Pristimantis hectus、P. laticlavius、Hyloscirtus alytolylax）的黏蛋白组样本。皮肤细菌群落通过16S rRNA V4区扩增子测序（Illumina MiSeq平台）进行表征，生成扩增子序列变异（ASV）并进行分类学注释。Bd感染状态与载量通过定量PCR（qPCR）检测。推测Bd抑制细菌通过将ASV序列与Woodhams等人（2015）的抗真菌分离株数据库及SILVA参考数据库联合比对进行鉴定。黏蛋白组对Bd生长的抑制作用通过体外挑战实验评估，使用2025年从拉普拉纳达分离的本地Bd谱系（属于高毒力全球泛流行谱系，BdGPL）。数据分析包括线性混合效应模型（LMM）、置换多元方差分析（PERMANOVA）、ANCOM-BC2与MaAsLin2差异丰度分析、Mantel检验及Blomberg's K系统发育信号检验等。
三、研究结果
3.1 Bd感染模式
qPCR检测显示，H. alytolylax（43%）、P. laticlavius（40%）和P. apiculatus（33%）具有最高的Bd流行率，而Centrolene peristicta（6019 zGE）和P. apiculatus（320 zGE）具有最高的平均感染载量。P. hectus和Sachatamia orejuela的流行率与载量均为最低（12%和10.3 zGE；14%和15 zGE）。这些结果与更大规模样本的研究一致，表明所选物种代表了不同的Bd易感性。
3.2 皮肤细菌微生物组的多样性与结构
共检测到14241个细菌ASV。物种间在纲、目、科、属水平均存在细菌组成差异，科和属水平的宿主物种特异性最为显著。S. orejuela和H. alytolylax拥有最高比例的特有ASV（分别为56.6%和29.5%），而多数ASV在3-5个物种间共享。旱季有23个特有ASV（41.1%），雨季有13个（23.2%），两季共享20个（35.7%）。具有水生幼体的物种有24个特有ASV（38.1%），直接发育物种有18个（28.6%），两繁殖模式共享21个（33.3%）。
3.3 皮肤细菌群落组成与结构的驱动因素
α多样性（ASV丰富度、Shannon指数、Faith系统发育多样性）随物种和季节变化，但不受繁殖模式、Bd感染状态或载量影响。S. orejuela的ASV丰富度和系统发育多样性显著较低，P. verecundus的Shannon多样性显著较低。仅考虑Bd阳性个体时，P. verecundus在所有三个指标上均显著降低，而H. alytolylax的ASV丰富度显著较高。雨季的Faith系统发育多样性显著高于旱季。β多样性分析显示，物种（未加权UniFrac R2=0.08，加权UniFrac R2=0.16）、季节（均为R2=0.02）和繁殖模式（未加权R2=0.02，加权R2=0.03）均显著影响群落组成，但Bd感染状态无显著效应。差异丰度分析鉴定出4个属在物种间差异显著（P. laticlavius富集Rhizobacter和Limnohabitans，后者为推测Bd抑制属；P. verecundus富集Roseateles），18个属在季节间差异显著（旱季以Flavobacterium、Acinetobacter、Bacteriovorax为主，前两者为推测Bd抑制菌），13个属在繁殖模式间差异显著（水生幼体物种富集Bordetella、Flavobacterium、Alkanindiges）。Mantel检验显示Bray-Curtis距离与宿主系统发育显著相关（r=0.46，p=0.03），Blomberg's K检验在PCoA轴1上检测到显著系统发育信号（K=1.05，p=0.04），但UniFrac距离未显示显著相关性。
3.4 推测Bd抑制细菌的存在
基于Woodhams等人（2015）数据库鉴定出73个推测Bd抑制ASV。其相对丰度在物种间差异显著（Kruskall-Wallis χ2=16.7，p=0.01），S. orejuela（7.7%）和C. peristicta（2.9%）最高，P. verecundus（0.5%）最低。旱季个体（3.2%）显著高于雨季（1%），具有水生幼体的物种（3.7%）显著高于直接发育物种（1.5%）。Stenotrophomonas是各分组中的优势推测抑制菌属。然而，推测Bd抑制细菌的丰富度和相对丰度均与Bd感染状态或载量无显著关联。
3.5 宿主皮肤黏蛋白组对Bd生长的抑制
黏蛋白组挑战实验显示，多数物种的Bd生长曲线与阳性对照高度相似，表明几乎没有抑制作用。仅P. laticlavius在实验后期表现出轻微抑制效应，其Bd承载能力（K）、生长速率（r）和曲线下面积（AUC）均较低。物种间Bd生长抑制率差异显著（χ2=61.7，p<0.001），H. alytolylax和P. hectus与P. laticlavius差异显著。值得注意的是，H. alytolylax的黏蛋白组显著增强了Bd生长（Mann-Whitney U，p=0.003）。实验结束时各处理组的Bd基因组当量（zGE）无显著差异。
四、讨论与结论
本研究表明，在拉普拉纳达这一经历Bd流行后仍存续的两栖动物群落中，宿主物种身份、系统发育、季节性和与栖息地密切关联的繁殖模式是皮肤细菌群落组成的主要驱动因素，而当前Bd感染状态和载量所解释的变异甚微。这一发现与部分研究一致，但也与另一些显示Bd驱动微生物组变化的报道形成对比。推测Bd抑制细菌在所有物种中均有检出但丰度较低（0.5%-7.7%），显著低于其他两栖动物群落的报道值（10%-20%）。尽管如此，这些细菌在未经历衰退的S. orejuela和C. peristicta中富集，而P. laticlavius——唯一丰度增加的物种——既富集了推测抑制菌属Limnohabitans，又表现出轻微的黏蛋白组抑制活性，提示特定有益细菌的存在和活性可能介导了微生物组相关的防御反应。H. alytolylax黏蛋白组促进Bd生长的趋势值得关注，因其在该地区疾病动态中扮演重要角色。
研究人员强调，当前数据无法确定Bd是否在历史上塑造了微生物组组装，但为未来评估微生物组的适应性变化提供了重要基线。研究结论指出，宿主免疫性状、病原体特征与环境条件之间的整合对于理解这一受威胁群落中宿主-病原体的长期共存至关重要。未来研究应整合黏蛋白组功能、宿主免疫基因活性与Bd谱系毒力相关基因表达，以阐明这些互作组分如何塑造长期共存格局。本研究为疾病金字塔框架在自然生态系统中的实证应用提供了范例，也为两栖动物保护策略的制定提供了微生物组层面的科学依据。
