太湖流域地表水E. coli中高风险抗生素耐药基因（ARGs）的流行、空间热点与共现网络

《iScience》：High-risk ARGs in environmental E. coli from the Taihu Lake watershed

　　

全球抗菌药物耐药性（AMR，antimicrobial resistance）危机因环境细菌储库而加剧。本研究调查了太湖流域118份地表水样本中E. coli携带的抗生素耐药基因（ARGs，antibiotic resistance genes）。研究人员采用mTEC琼脂和uidA PCR分离E. coli，并针对覆盖15类耐药类别的92种ARGs对分离株进行筛查。来自湖州的样本ARG检出率高于其他地区样本，以氨基糖苷类、磺胺类和糖肽类耐药基因占优势。城市采样点构成ARG热点。共现分析揭示喹诺酮类与氨基糖苷类耐药基因之间存在关联，提示水平基因转移（HGT，horizontal gene transfer）。风险评估将25.3%的样本归类为高风险。这些观察结果说明当地水环境中存在人为驱动的ARG污染，并支持针对性监测和面向One Health（同一健康）的缓解策略。

研究背景与问题：全球AMR危机因环境储库中耐药细菌和耐药基因的传播而放大，地表水可接收污水处理厂出水、农业径流和城市径流，成为ARGs和抗生素耐药细菌（ARBs，antibiotic-resistant bacteria）的储库与传播通道。大肠杆菌（E. coli，Escherichia coli）既是粪便指示菌，也是ARGs的重要储库和监测哨兵。湖州位于长江三角洲、毗邻太湖，受城市化、工业、农业和污水收集处理不足等多重压力；但该地区地表水中E. coli相关ARGs的流行与分布仍缺乏系统本地数据，限制了对公共健康和生态风险的评估。因此，研究人员开展本研究，旨在分离鉴定湖州地表水E. coli，刻画其ARGs丰度与多样性、空间分布及与污染源邻近关系，分析ARGs与移动遗传元件（MGEs，mobile genetic elements）的关联，并用风险指数、关键ARG检出频率和聚类进行风险评估；研究结果可为针对性监测和One Health导向的缓解策略提供依据。

关键技术方法（不超过250字）：研究于2023年采集中国太湖流域湖州及非湖州地表水样本118份，来源于河流、湖泊、溪流、山泉和自来水等。研究人员用mTEC琼脂分离推测E. coli，以uidA PCR确认；每份水样所有确认菌落DNA合并，代表该样本进行靶向PCR筛查92种ARGs，并检测MGEs和质粒标记。数据分析采用描述统计、主成分分析（PCA，principal component analysis）、k-means聚类、风险指数（每样本检出ARG数）和检出频率分析。

研究结果：
1. 采样与E. coli分离概况：研究人员在2023年开展跨区域地表水调查，湖州核心区覆盖28个点位并含四季采样，另纳入非湖州样本，共形成118份样本；通过mTEC培养和uidA PCR确认E. coli，并将每样本确认菌落DNA合并用于ARG筛查。
2. E. coli丰度与时空格局：所有采样点均检出E. coli，总体回收成功率70%，浓度范围从<1至>4×10⁴ CFU/100 mL；城市影响点位最高（中位数4.0×10⁴ CFU/100 mL），郊区点位较低（中位数348.5 CFU/100 mL），夏季最高（中位数4.5×10⁴ CFU/100 mL），冬季最低（中位数265 CFU/100 mL），表明其受人为活动和季节气候驱动。
3. 湖州ARGs流行率更高且聚类区别于非湖州地区：PCR筛查显示湖州总体ARG检出率均值28.86%±9.91%，显著高于非湖州17.50%±10.85%；湖州优势类为庆大霉素（58.24%±28.68%，如ant(2″)-I）、磺胺甲噁唑（41.0%±28.46%，如sul1/2）和糖肽类（39.46%±27.92%，vanA达71.28%），多黏菌素类低（4.4%±9.31%）；非湖州高检出主要由庆大霉素、新霉素、链霉素和磺胺甲噁唑相关基因驱动。湖州ARGs层次聚类形成5簇，簇1（13基因，平均66.2%，如blaCTX-M-1）和簇2（16基因，52.9%，如blaFOX）为主，簇3（45基因，7.9%）为低丰度潜在新兴威胁；非湖州形成4簇，簇2（9基因，58.4%）为高风险组。比较散点图显示多数基因在湖州检出率更高，FDR校正后显著差异组差异更大。
4. 人为热点与基因特征的空间分布：按7类区域分析92种ARGs，大学平均检出率最高（33.48%±33.56%），城市公园次之（29.08%±23.83%），郊区公园最低（19.56%±26.57%）；大学与城市公园ARG流行率呈正相关但离散度大。热图显示社区富集gltA和ant(2″)-I，大学富集E. coli gyrA和parC，城市河流富集blaCMY/blaFOX，城市公园富集gyrA/blaCTX-M-1；sul2和emrE在部分地区低，gyrA、parC和blaCTX-M-9广泛高。雷达图显示gyrA、vanA和aac(3)-IV在大学及郊区公园高，ant(2″)-I、blaCTX-M-9和parC在社区和城市河流高，郊区湖多数基因低，表明空间异质性和位点特异性选择压力。
5. 共现网络与多重耐药（MDR，multidrug resistance）基因协同：湖州87份样本中多数携带多个ARGs，单样本检出ARG数从4.35%（4/92）至61.96%（57/92）；92种ARGs中89种至少检出一次，平均每种在25.96份样本中检出，ant(2″)-I最流行（68份）。共现网络前5强关系为E. coli gyrA与E. coli parC共现61次、ant(2″)-I与gyrA共现55次、ant(2″)-I与parC共现55次、aac(3)-IV与gyrA共现52次、aph(3′)-III与gyrA共现52次；gyrA与parC接近完全共现，提示共享HGT机制或相似选择压力。
6. 高风险样本、关键基因与聚类风险评估：风险指数（每样本检出ARG数）范围4—57，均值26.6，中位数26，第75百分位32.5；22份样本（25.3%）风险指数≥32.5，为高风险。检出频率最高为ant(2″)-I（78.2%），高频基因（>70%）包括E. coli gyrA和parC（均74.7%）及vanA（71.3%）；中等频率（50%—70%）包括aph(3′)-III、aac(3)-IV、blaVIM等。PCA前5个主成分累计解释29.2%方差，提示样本差异复杂；k-means聚类（k=3，轮廓系数0.0623）得到簇0（16份，平均风险指数32.98，高风险）、簇1（23份，18.05，低风险）和簇2（48份，32.68，高风险），簇间ARG谱不同。

讨论总结：作者指出，湖州ARG平均检出率高于非湖州，提示局部人为活动（污水、农业等）可能促进ARG积累；氨基糖苷类、磺胺类和糖肽类高流行反映人医和农业使用形成的持续选择压力，多黏菌素类低提示区域使用有限；ant(2″)-I和vanA等高检出说明环境储库作用，但其临床相关性仍需验证。聚类显示高丰度簇富集临床重要β-内酰胺酶和氟喹诺酮标记。空间梯度与土地利用强度、污水输入和人群密度相关，城市河流富集β-内酰胺酶基因，大学位点富集磺胺类基因和质粒标记，小尺度亦存在异质性，提示需本地化监测。共现关系提示HGT及质粒（如IncF、IncI）和1类整合子等MGEs驱动ARG共定位与传播；城市热点中最后防线抗生素耐药基因（如blaVIM、blaNDM、vanA/vanB、mcr变异体）的存在构成暴露风险。风险框架识别25.3%样本为高风险，聚类为管理优先级提供依据。作者同时指出局限：靶向PCR不量化丰度、表达或功能耐药，存在交叉反应可能；未进行基因组测序、抗生素残留定量和接合实验，ARG遗传背景与可移动性仍待研究。

研究结论部分翻译：总之，本研究利用靶向PCR对湖州地表水E. coli中ARGs的流行、空间分布和共现进行了初步表征。研究识别出显著的ARGs遗传负荷，城市地区流行率更高，并提示潜在的人为驱动因素和HGT热点。这些发现强调需要针对性监测和未来研究（如宏基因组学和测序）来验证风险并指导One Health策略。