传统多样性指标低估氟苯尼考暴露下尼罗罗非鱼（Oreochromis niloticus）肠道微生物组扰动：16S rRNA测序与共现网络拓扑分析

《Scientific Reports》：Conventional diversity metrics underestimate gut microbiome disruption under florfenicol exposure in Nile tilapia (Oreochromis niloticus)

肠道微生物组（gut microbiome）是鱼类健康和生长性能的重要组成部分，因此是水产养殖中的重要生产因素。抗生素（antibiotics）是主要应激源，不仅通过扰乱微生物多样性，还通过改变参与消化、定植抵抗（colonization resistance）和功能冗余（functional redundancy）的微生物相互作用。然而，当前风险评估主要依赖基于多样性的评价（diversity-based evaluations），很少考虑网络拓扑（network topology）的变化。在此，研究人员利用16S rRNA测序（16S rRNA sequencing），将传统多样性指标（conventional diversity metrics）与微生物共现网络建模（microbial co-occurrence network modelling）相结合，研究了氟苯尼考（florfenicol, FFC）对尼罗罗非鱼（Oreochromis niloticus）肠道微生物组的影响。尽管微生物群落在多样性水平上总体保持稳定，处理仅解释观测变异的一小部分，但网络分析（network analysis）揭示了显著碎片化（fragmentation）和较不稳健的网络结构。网络连通性（network connectivity）显著降低，氟苯尼考网络中中心枢纽类群（central hub taxa）几乎完全更替。研究人员的发现表明，仅基于多样性的评估可能大幅低估抗生素诱导的微生物组扰动（antibiotic-induced microbiome disturbance）程度，而网络分析能捕获群落内部对微生物组稳定性高度相关、若评价仅限于多样性指标否则仍无法检测到的变化。

商业水产养殖依赖集约化养殖条件下养殖鱼类的良好生长性能，以满足全球对水产品日益增长的需求。鱼类依赖与肠道微生物组（gut microbiome）的共生关系来实现最佳发育和应激条件下恢复力，这一点日益被认为是可持续养殖和兽医健康管理的关键。细菌通过参与营养消化，包括短链脂肪酸（short-chain fatty acid, SCFA）产生和必需维生素如钴胺素（维生素B₁₂）合成、免疫调节和针对病原体的定植抵抗（colonization resistance）来支持宿主生理。微生物相互作用包括合作性行为，如代谢交叉喂养或抗生素抗性传递，也包括竞争和拮抗相互作用。维持正负相互作用平衡对核心微生物组功能和群落适应扰动至关重要。氟苯尼考（florfenicol, FFC）是广泛使用的兽用抗生素，可治疗淡水养殖中鲴? 重要细菌病，如链球菌病和弧菌病。除抗生素抗性基因（antibiotic resistance genes, ARGs）传播风险外，研究还报道其对宿主和微生物组有组织病理学影响。抗生素诱导的微生物组扰动不能由单一通用失调标志概括。目前多数研究主要依赖α多样性、β多样性和差异丰度等标准多样性指标，难以捕获微生物-微生物相互作用层面的结构变化。由于丰度与网络影响并非严格耦合，稀有类群也可能作为关键类群对生态系统稳定性产生不相称影响。因此，仅关注组成分析可能低估抗生素处理后微生物组扰动程度。罗非鱼是全球养殖最广泛的鱼类之一，但在微生物组研究中相对关注不足。为此，研究人员将成年尼罗罗非鱼（Oreochromis niloticus）暴露于每日10 mg/kg bw的氟苯尼考处理10天，随后进行20天无干预恢复，并基于16S rRNA测序的V1-V2区，结合传统多样性指标和微生物共现网络分析，表征肠道微生物组动态。

关键技术方法简述：研究人员采用16S rRNA V1-V2区扩增子测序，使用DADA2和SILVA v138进行ASV（amplicon sequence variant，扩增序列变异）推断与分类注释；以Shannon指数评估α多样性，以Bray–Curtis距离和PERMANOVA评估β多样性，并用beta-binomial回归识别时间响应类群；在属水平采用SPIEC-EASI构建控制组和FFC组共现网络，计算度、介数、邻居中心性、模块性、枢纽（hubness）和IVI（Integrated Value of Influence，综合影响值），并通过随机和靶向节点去除模拟评估网络鲁棒性。样本队列来源：成年尼罗罗非鱼来自巴西Santa Bárbara, Zacarias, S?o Paulo商业鱼场André Luiz Scarano Camargo e Outro，未接受过抗菌治疗；实验在巴西圣保罗大学水产中心（Jaboticabal, SP, Brazil）开展。

结果一：氟苯尼考对肠道微生物组多样性的影响有限。通过16S rRNA测序共检测10,028个ASV，代表31门854属；对照组和FFC组共享相似优势属，包括Cetobacterium、Romboutsia、Streptococcus、Diaminobutyricimonas和Bacteroides。Shannon多样性随时间显著变化，处理主效应不显著，时间×处理交互显著；β多样性中时间解释最多，处理主效应显著但很小，交互反映时间轨迹差异。差异丰度分析显示，854属中对照组31属、FFC组37属随时间变化，25属共享；43属在特定时间点组间显著但无贯穿实验的一致变化。最大响应数出现在处理第5天和处理后7天。FFC与Diaminobutyricimonas、Sarcina、Lactobacillus、Romboutsia和Pseudomonas等随时间更强耗竭相关；Aeromonas、Pantoea、Methylobacterium-Methylorubrum和Bacteroides初期下降后反弹；Bacillus、Streptococcus、Phreatobacter、Sphingomonas、Rhodobacter和Staphylococcus暂时升高。核心群落在两组均稳定，差异主要限于低丰度和低流行类群。

结果二：氟苯尼考暴露下网络稳定性和相互作用降低。研究人员将各时间点样本按处理组汇合，构建属水平共现网络，控制组n=30、FFC组n=33。两组节点数相似（261和268），正关联比例约74%，但FFC网络边数更少（409 vs. 516）。模块性差异较小（0.57 vs. 0.64），但度（degree）和邻居中心性显著降低，导致网络更稀疏、更碎片化；介数无显著差异，但FFC网络节点值分布更异质。FFC网络平均路径长度增加，模块内边比例升高。节点去除模拟中，FFC网络相对最大连通分量（largest connected component, LCC）下降更快、AUC（area under the curve，曲线下面积）更低，尤其靶向去除高连接节点时，说明结构鲁棒性降低、更易受节点丢失影响。阶段分层网络敏感性分析提示FFC网络在恢复期更强碎片化，但因样本量小且不等、边选择稳定性较低，未支持稳健的阶段比较。

结果三：氟苯尼考处理下中心枢纽节点转移。基于hubness评分并以IVI二次排序，研究人员在控制组识别23个中心枢纽，在FFC组识别22个。所有节点中相对丰度与hubness评分无相关。控制组枢纽主要属于Firmicutes、Proteobacteria和Bacteroidota；FFC处理后，除Clostridium sensu stricto 2和Proteiniclasticum外，所有控制组枢纽失去中心网络位置。Dechloromonas和Actinobacillus不再符合枢纽标准，但从外周位置转为连接不同子模块的桥节点。FFC组枢纽主要分布于Firmicutes、Actinobacteriota和Proteobacteria。控制组枢纽仅4个与核心属重叠，即Cetobacterium、Clostridium sensu stricto 1、Rhizobium和Lactobacillus；FFC组仅Streptococcus同时为枢纽和核心属。丰度变化未一致转化为网络重要性变化；Streptococcus丰度升高且网络重要性升高，而Lactobacillus和Hyphomicrobium丰度降低且网络影响剧降。

结果四：氟苯尼考诱导关联伙伴强更替。基于两组共享的209个节点，FFC处理下79个关联丢失、54个新关联形成。边丢失主要影响高连接枢纽如Aggregatibacter、Lactobacillus、Macellibacteroides和Hyphomicrobium，也影响非枢纽类群如Murdochiella和Bacteroides。关联增加最多者多为FFC网络中新枢纽，如Erysipelothrix、Erysipelatoclostridium、Frisingicoccus和Porphyromonas。保留中心位置的Clostridium sensu stricto 2和Proteiniclasticum，其直接关联伙伴完全替换。原有功能连贯簇被替代，例如Clostridium sensu stricto 2与Succiniclasticum、Alloprevotella的连接在FFC网络中替换为Proteus和Eubacterium nodatum group；Lactobacillus失去与Limosilactobacillus、Aurantimicrobium的关联，转而与Porphyromonas关联；Bacteroides与Fusobacterium、Aeromonas、Sutterella的中心连接缺失；Cetobacterium从与Macellibacteroides关联转为与Firmicutes Family XIII AD3011 group关联。还出现共生菌与潜在病原体属之间的连接变化，如Streptococcus与Macellibacteroides负关联。新增关联包括Erysipelothrix与Succiniclasticum、Pasteurella、Alloprevotella，以及Proteiniclasticum与Butyricicoccus、Prevotella 7。上述拓扑重排提示FFC网络中功能互补性可能降低。

讨论部分总结：FFC在常规多样性指标下对成年尼罗罗非鱼肠道微生物组影响中等，时间是最强影响因素。差异主要来自Diaminobutyricimonas、Lactobacillus、Romboutsia、Bacteroides和Pseudomonas等常见淡水鱼肠道成员短暂耗竭，这些类群可能参与复杂碳水化合物降解、脂质消化和SCFA产生。尽管这些共生类群在实验结束时恢复丰度，暂时失衡仍可能损害肠道屏障完整性和代谢，使宿主更易受后续应激影响。Han et al.曾报道锦鲤FFC处理后尽管微生物组多样性恢复，肠道损伤持续存在。机会病原富集是抗生素失调常见特征；本研究中仅Streptococcus sp.和Aeromonas sp.短暂增加，且不贯穿实验。二者包含高经济关注鱼类病原体，如S. iniae、S. agalactiae和A. hydrophila，其耐FFC可能与多重耐药决定因子和R策略有关。Citrobacter和Plesiomonas未检出。全长16S rRNA可改善种水平分辨率，但确认病原株需shotgun宏基因组或靶向培养。研究结果与尼罗罗非鱼、channel catfish、salmon和European seabass等研究一致；差异可能来自鱼龄、抗生素种类、温度依赖药代动力学、肠道区室和养殖系统。FFC及代谢物经粪便和未吃饲料进入水体，在封闭或循环系统中可能延长暴露。实验室条件下FFC是唯一主要应激源，商业系统存在多种压力；且因缺乏缸级独立重复，不能完全排除缸效应。丰度持续并不说明生态关系保留。尽管多样性变化小，FFC使网络更碎片化、更区室化，网络连通性降低；类似效应在鼠和人中有报道，但鱼中此前未如此详细。网络复杂性降低与稳态维持能力降低和功能效力下降有关。丢失的关联涉及营养网络互补底物利用，如Alloprevotella产琥珀酸、Succiniclasticum将琥珀酸转化为丙酸、Clostridium sensu stricto产丁酸；Lactobacillus与Limosilactobacillus的互养关系；Streptococcus与Macellibacteroides负连接消失可能提示病原抑制减弱。Streptococcus成为FFC网络中连接度最高的新枢纽，Erysipelothrix等机会病原体也类似。新形成的关联涉及Proteiniclasticum与Butyricicoccus、Prevotella，以及Clostridium sensu stricto 2与Proteus，多为低丰度类群，可能不能功能等价补偿丢失的代谢协同。共现关系不一定代表直接生态相互作用，需共培养验证。这些观察可能解释抗生素处理罗非鱼和海鲈中饲料利用、脂质消化和生长降低。枢纽类群大量更替；Cetobacterium、Clostridium和Lactobacillus等控制组枢纽有支持宿主生理的作用。Cetobacterium可合成钴胺素（维生素B₁₂），补充Cetobacterium可增强对A. hydrophila的定植抵抗和肠道屏障完整性。Aureispira可能通过ixotrophic predation抑制病原体，已在Vibrio spp.中证明，且在Japanese eel和swan mussels中检测到，低丰度高连接，是候选关键类群。Hyphomicrobium和Dechloromonas为环境相关短暂低丰度类群。枢纽丢失或替换可能致网络崩溃。节点丢失模拟显示FFC网络失去连接更快，靶向去除枢纽最明显。现实中重复抗生素、共同治疗、继发感染或环境应激可造成二次节点丢失。枢纽身份和网络位置与丰度不严格相关，因此多样性分析会漏检生态重要类群。Cetobacterium对抗生素响应不一致，本研究无变化；先前brown trout研究显示FFC后丰度恢复但中心性仍降低。当前肠道失调定义主要依赖多样性，很少纳入网络拓扑；网络构建无金标准、计算需求大且需较大样本量。每时间点重复不足，限制阶段或时间点分层网络推断。

研究结论翻译：当仅考虑多样性时，研究人员的结果提示，FFC处理未在成年尼罗罗非鱼中诱导明确的肠道微生物组失调（dysbiosis）模式或持续性微生物组转变。然而，对微生物共现结构的分析显示，尽管多样性相似，FFC使微生物组转向一种失稳且更脆弱的群落状态，其功能也可能不同于初始状态。从产业角度看，上述协同相互作用丧失导致的营养利用受损和恢复力降低，可能导致生长下降和死亡率增加，最终降低产量。多样性指标和网络拓扑捕获微生物组稳定性的不同维度；一个维度的恢复力不能直接外推至另一个。因此，仅评估多样性可能低估抗生素暴露后微生物组扰动程度。不应将组成稳定性和网络稳定性视为彼此分离的概念，而应将其视为互补层面，二者均为全面理解微生物组响应和可靠风险评估所必需。未来研究整合这两个层面，还可能有助于更稳健地定义健康鱼类微生物组基线，促进区分水产养殖中真生态（eubiotic）和失调（dysbiotic）状态的合适生物标志物识别，并支持开发新的益生菌候选物以替代反复化学干预。