黄河中游喹诺酮类抗生素（QNs）时空动态驱动细菌群落演替及洪水季节功能适应

《Environmental Technology & Innovation》：Spatiotemporal dynamics of quinolone antibiotics drive bacterial community succession in the Yellow River

　　

喹诺酮类抗生素（quinolones, QNs）的时空动态驱动黄河细菌群落演替。QNs降低物种均匀度（species evenness）和系统发育多样性（phylogenetic diversity）。细菌群落在较高QN胁迫下增强正负种间关联（interspecific association）以维持稳健性（robustness）。细菌群落通过春季洪水期间维持生存稳态（survival homeostasis）适应更重QN污染。夏季洪水期间群落功能向高效污染物降解转变。

研究背景：人类活动扩大使抗生素广泛存在，痕量水平（ng/L或μg/L）即可诱导非靶标毒性，并选择抗生素抗性细菌（antibiotic-resistant bacteria, ARB）和抗生素抗性基因（antibiotic resistance genes, ARGs）。喹诺酮类抗生素（quinolones, QNs）是全球消费前五的抗生素，约占市场17%，因广谱、药代动力学好、生物利用度高和成本低而广泛使用；其高水溶性促进环境迁移。QNs除诱导遗传突变外，还可触发非遗传代谢调节，如糖酵解、三羧酸循环（TCA cycle）或群体感应（quorum sensing, QS）系统，导致表型耐受。细菌群落可通过种间互作应对抗生素压力。河流是抗生素污染传播通道。黄河中游具有双洪水 regime：春汛（3–4月，冰坝驱动，低流量）和夏汛（季风驱动，高湍流，高泥沙，外源输入）。但洪水季节环境过滤与QN污染如何塑造细菌群落组装尚不清楚。研究人员提出假设：QNs动态季节不同，春汛更严重；重QN污染作为强确定性过滤改变生态位-中性平衡；群落功能可塑性从春季抗性稳态转向夏季降解代谢。

研究方法概括：研究人员在黄河中游河南段及支流设12个站点，包括污水处理厂（wastewater treatment plants, WWTPs）排放口LB和HY，于2023年3月和7月采集水样。样本队列来源为黄河中游河南段12站点。检测7种优先QNs：norfloxacin、lomefloxacin、ciprofloxacin、pefloxacin、fleroxacin、enrofloxacin和ofloxacin，采用固相萃取（solid-phase extraction, SPE）和超高效液相色谱-串联质谱（ultra-high-performance liquid chromatography-tandem mass spectrometry, UHPLC-MS/MS）。微生物用16S rRNA V3–V4扩增和Illumina MiSeq测序，共118样本。分析包括QIIME 2 alpha多样性、双因素方差分析（analysis of variance, ANOVA）、Bray-Curtis主坐标分析（principal coordinates analysis, PCoA）、PERMANOVA、LEfSe、中性群落模型（Neutral Community Model, NCM）、β最近分类群指数（β-nearest taxon index, βNTI）与Raup-Crick度量（RC_bray）、SparCC共现网络、NetCarto关键类群和PICRUSt2功能预测。

结果如下。

3.1 QN污染在春汛更严重。检测显示norfloxacin检出率最高（91.67%），平均17.29 ± 10.73 ng/L，SS春汛峰值42.68 ng/L；enrofloxacin检出率最低（29.17%），集中在WWTP排放口。lomefloxacin仅春汛检出，浓度1.98–91.81 ng/L，夏汛未检出。pefloxacin、ofloxacin和fleroxacin春汛检出率高于夏汛；ciprofloxacin检出率相似，但春汛浓度（12.19 ± 11.09 ng/L）高于夏汛（0.41 ± 0.63 ng/L）。结论为春汛QN污染更重，WWTP排放口和SS为热点。

3.2 细菌群落组成与多样性。春汛Proteobacteria（50.79 ± 10.63%）主导，Bacteroidota和Verrucomicrobiota为生物标志；夏汛Actinobacteriota（42.30 ± 10.32%）主导，Cyanobacteria、Bdellovibrionota和Gemmatimonadota为生物标志。科水平上，Sporichthyaceae和Ilumatobacteraceae春汛低于夏汛，Comamonadaceae、Rhodobacteraceae、Sphingomonadaceae和Flavobacteriaceae春汛较高。Observed species和Chao1两季无显著差异，Shannon、Pielou’s evenness和Faith’s系统发育多样性（phylogenetic diversity, PD）夏汛更高；HY、LB和SS春汛alpha多样性较低。

3.3 细菌群落的生态位-中性平衡。LEfSe识别春汛637个、夏汛884个生物标志。NCM显示生物标志中性分布比例春汛52.75%、夏汛53.96%，偏离中性，受环境过滤；非生物标志拟合较好，春汛97.55%、夏汛97.86%，背景群落可能随机。βNTI春汛显著高于夏汛。空模型显示春汛均质选择38.78%、均质扩散36.12%、生态漂变22.75%；夏汛均质选择65.20%，均质扩散24.07%、生态漂变10.45%。结论为春汛由均质选择和中性过程共同组装，夏汛由强均质选择驱动。

3.4 网络结构与关键类群。共现网络显示夏汛参与节点比例55.94%高于春汛34.17%；春汛网络连接度和聚类系数更高，模块性更低；正/负边比两季一致（1.21 ± 0.16）。8.08% taxa为关键类群，包括connectors、module hubs和network hubs。module hubs多由优势taxa组成，connectors多非优势。结论为春汛细菌关联更强，关键类群角色时空变化。

3.5 两个洪水季节潜在指示类群和功能通路。整合NCM、LEfSe、网络和功能预测，春汛指示类群与核心代谢（TCA循环、糖酵解、氨基酸/核苷酸/血红素生物合成）和应激响应（磷壁酸、decaprenyl phosphate、ppGpp）相关，可能关联外排、靶标保护或持留菌形成。夏汛指示类群与芳香化合物降解相关，如gallate、vanillin、toluene、catechol、salicylate、nitrobenzoate，Comamonadaceae、Moraxellaceae和Weeksellaceae正相关。结论为春汛偏被动抗性，夏汛偏合作降解。

3.6 细菌群落演替的潜在驱动因子。结构方程模型（structural equation modelling, SEM）显示春汛压力主要来自lomefloxacin，与Candidatus Planktophila强负相关（r = -0.64），可能降低稀有类群丰富度和Shannon多样性。夏汛由ofloxacin、pefloxacin、norfloxacin和fleroxacin共同驱动，增强Predibacter starrii和Aquaspirillum serpens竞争力，但稀有类群丰富度和Faith’s PD与ofloxacin、pefloxacin负相关。Mantel检验显示春汛核心代谢与总有机碳（total organic carbon, TOC）、溶解氧（dissolved oxygen, DO）、Cl^-等正相关；夏汛芳香化合物降解与总氮、温度、ofloxacin正相关。结论为春汛可能强化核心代谢或共代谢，夏汛高温促进代谢并驱动靶向降解。

讨论部分总结：4.1 春季lomefloxacin与夏季多QNs污染。春汛污染更重，因春季人兽用药多、温度约10℃、降解低；夏汛光解和微生物降解增强。主河道与支流浓度相当，除夏汛ofloxacin。4.2 生态位-中性动态改变多样性和种间关联。约90%群落成员符合NCM，但组装过程季节不同。春汛高QN下中性过程贡献超过50%，其他物理扰动和外源输入可能压倒环境过滤；夏汛低QN但其他理化因子驱动均质选择，促进远缘类群共存和PD。重lomefloxacin与种间关联增强相关，正负比不变，提示增强连接本身而非转向促进维持胁迫下稳定。4.3 功能从春季生存维持向夏季高效降解转变。春季功能模块化维持生存稳态，夏季转向芳香化合物等高效降解。4.4 候选指示类群潜在生态角色。春汛为Acinetobacter lwoffii和Ca. Planktophila；夏汛为Peredibacter starrii、Acinetobacter variabilis和Aquaspirillum serpens。各具网络角色和可能功能。这些通过生物信息学识别，不建立因果，需实验验证。

结论翻译：本研究调查QN污染对黄河浮游细菌群落的生态效应。关键结论为：春汛受更重QN污染，可能通过生态位-中性动态增强种间关联，同时降低物种均匀度和系统发育多样性。细菌群落可能通过春汛功能模块化维持生存稳态，并在夏汛转向高效污染物降解，从而季节适应QN胁迫。这些发现揭示抗生素驱动细菌群落季节演替机制，支持靶向风险管理。未来工作可关注候选指示类群间精确种间合作和拮抗，以及QS等机制在抗生素抗性发展中介导这些互作。实验室培养、宏基因组学和代谢组学可整合解析代谢通路并构建物种-功能网络。