环境科学  2023, Vol. 44 Issue (1): 158-168   PDF    
长江口近岸地区抗生素抗性基因与微生物群落分布特征
徐秋鸿1, 刘曙光1,2, 娄厦1,2, Larisa Dorzhievna Radnaeva3, Elena Nikitina3, Makhinov Aleksei Nikolavich4, Araruna José Tavares5, 金宇辰1, 李鑫1     
1. 同济大学水利工程系, 上海 200092;
2. 长江水环境教育部重点实验室, 上海 200092;
3. Laboratory of Chemistry of Natural Systems, Baikal Institute of Nature Management of Siberian Branch of the Russian Academy of Sciences, Republic of Buryatia, Republic of Buryatia 670000, Russian;
4. Institute of Water and Ecology Problems, Far East Branch, Russian Academy of Sciences, Khabarovsk 680000, Russia;
5. Pontifical Catholic University of Rio de Janeiro, Rio de Janeiro 999074, Brazil
摘要: 由于抗生素的大量使用, 环境中微生物对抗生素的抗性不断增加, 抗生素抗性基因(ARGs)问题越来越严重, 严重威胁生态安全和人类健康.为研究长江口近岸地区水体和底泥沉积物中的ARGs和微生物群落的分布特征, 通过野外调查采集了8个站点的水样和沉积物样本, 对2种磺胺类抗性基因(sul1sul2)、6种四环素类抗性基因(tetMtetCtetXtetAtetOtetQ)、1种整合子基因intI1、16S rRNA基因和微生物群落进行检测分析.结果表明, 长江口近岸地区10种抗性基因的检出率为100%.其中, 整合子基因intI1和水样中多种ARGs呈显著正相关关系.变形菌门(Proteobacteria)和拟杆菌门(Bacteroidota)是长江口近岸地区水环境中的优势菌门; Chloroplast为水体中的主要菌属, ChloroplastNitrospira为沉积物中的主要菌属.在水体中, 硝化螺旋菌门(Nitrospirota)是4种四环素类抗性基因(tetXtetAtetOtetQ)共同的潜在宿主; 在沉积物中, Sva0485是sul1intI1的共同潜在宿主.微生物群落的分布是长江口近岸地区ARGs迁移转化的重要影响因素.
关键词: 抗生素抗性基因(ARGs)      荧光定量PCR(qPCR)      16S扩增子测序      长江口近岸地区      微生物群落     
Distributions of Antibiotic Resistance Genes and Microbial Communities in the Nearshore Area of the Yangtze River Estuary
XU Qiu-hong1 , LIU Shu-guang1,2 , LOU Sha1,2 , Larisa Dorzhievna Radnaeva3 , Elena Nikitina3 , Makhinov Aleksei Nikolavich4 , Araruna José Tavares5 , JIN Yu-chen1 , LI Xin1     
1. Department of Hydraulic Engineering, Tongji University, Shanghai 200092, China;
2. Key Laboratory of Yangtze River Water Environment, Ministry of Education, Shanghai 200092, China;
3. Laboratory of Chemistry of Natural Systems, Baikal Institute of Nature Management of Siberian Branch of the Russian Academy of Sciences, Republic of Buryatia, Republic of Buryatia 670000, Russian;
4. Institute of Water and Ecology Problems, Far East Branch, Russian Academy of Sciences, Khabarovsk 680000, Russia;
5. Pontifical Catholic University of Rio de Janeiro, Rio de Janeiro 999074, Brazil
Abstract: Due to the extensive use of antibiotics, the resistance of microorganisms to antibiotics in the environment is increasing, and the problem of antibiotic resistance genes (ARGs) is becoming more and more severe, which seriously threatens ecological security and human health. In order to study the distribution characteristics of ARGs and the microbial community in different media in the coastal area of the Yangtze River Estuary, water and sediment samples from eight sites were collected through a field investigation. Two sulfonamide resistance genes (sul1, sul2) and six tetracycline resistance genes (tetM, tetC, tetX, tetA, tetO, and tetQ), one integrase gene intI 1, 16S rRNA gene, and the microbial community were detected and analyzed. The results showed that the detection rate of 10 resistance genes in the coastal area of the Yangtze River Estuary was 100%. intI1 was significantly positively correlated with various ARGs in the water samples. Proteobacteria and Bacteroidota were the dominant bacteria phyla in the water environment of the Yangtze River Estuary. Chloroplast was the main bacteria genus in water, and Chloroplast and Nitrospira were the main bacteria genera in sediment. In water, Nitrospirota was the common potential host of four tetracycline resistance genes (tetX, tetA, tetO, and tetQ). In sediments, Sva0485 was a potential host community shared by sul1 and intI1. The distribution of the microbial community was an important factor affecting the migration and transformation of ARGs in the nearshore area of the Yangtze River Estuary.
Key words: antibiotics resistance genes(ARGs)      quantitative-polymerase chain reaction(qPCR)      16S amplicon sequencing      nearshore area of the Yangtze River Estuary      microbial communities     

2006年, Pruden等[1]提出将抗生素抗性基因(antibiotic resistance genes, ARGs)作为一种新型的环境污染物, 世界卫生组织也于2015年明确提出将微生物抗性作为全球健康风险.全世界范围内越来越多的研究人员开始关注ARGs诱发的相关问题.近年来, 抗生素的滥用导致微生物(包括病原菌和非病原菌)对抗生素的抗性逐年增加, 环境中的细菌可以通过基因突变和代系遗传等方式获得抗性基因, 也可以通过转导、接合和转化等方式获得外源存在的抗性基因[2].也有研究表明, 环境中较低浓度的抗生素可以作为微生物种间或者种内的信号分子, 从而调控微生物群落对抗生素进行选择性适应[3].抗生素的存在会加速上述过程的发生, 增加环境中抗生素抗性出现的频率, 对环境安全和人类健康造成潜在威胁.近年来, 我国各大流域都已存在不同程度的ARGs污染[4~7].

海岸带尤其是河口地区, 是连接陆地与海洋的重要过渡区域[8], 接纳了大量来自河流径流和污水排放的污染物, 包括多环芳烃[9]、重金属[10]、抗生素[11]和ARGs[12].含有ARGs的污水最终会通过河口排入海洋环境, 并在人类活动和水文气象等多种因素的耦合影响下[13, 14], 发生一系列复杂的物理-化学-生物作用过程.抗生素抗性基因导致细菌产生耐药性, 并可通过水平转移在细菌间传播扩散[15], 严重威胁生态系统安全和人类健康[16].因此, 厘清河口区域水环境中ARGs和微生物群落的分布特征至关重要.

本文选取长江口近岸地区作为研究区域, 采用荧光定量PCR(quantitative-polymerase chain reaction, qPCR)分析水体和沉积物中共10种ARGs的分布规律, 包括2种磺胺类抗性基因(sul1sul2)、6种四环素类抗性基因(tetMtetCtetXtetAtetOtetQ)以及1种整合子基因intI1和16S rRNA基因, 并采用16S扩增子测序技术分析水环境中微生物群落的分布特征, 揭示长江口近岸地区ARGs与微生物群落间的影响机制, 以期为控制长江口近岸地区水环境中ARGs污染提供数据支持和科学依据, 并为长江口近岸地区ARGs生态风险评价提供参考.

1 材料与方法 1.1 样品采集

本文采集了长江口近岸地区8个采样点的表层水和沉积物样品, 具体采样点分布如图 1所示.自北向南分别是浏河口(1)、石洞口(2)、吴淞口(3)、竹园排污口(4)、三甲港水闸(5)、朝阳农场(6)、大治河(7)和南汇东滩(8).其中浏河口(1)是浏河与长江的交汇地, 吴淞口(3)是黄浦江与长江的交汇地, 大治河(7)是大治河与长江交汇地, 石洞口(2)和竹园(4)位于城市污水处理厂的排污口, 朝阳农场(6)属于潮滩环境, 三甲港水闸(5)位于九段沙湿地国家自然保护区的下游断面, 南汇东滩(8)属于湿地保护区.

图 1 采样点分布示意 Fig. 1 Location of sampling sites

有研究表明, 温度对ARGs的传播扩散起着重要的作用[17].因此本次采样时间为2021年7月14日.每个采样点采集表层水(2 L)并装入超纯水洗过的聚乙烯瓶中, 同时采集100 g沉积物(0~5 cm), 并转移至灭菌的塑料袋中.每个采样点采集3个平行样品, 所有样品均用冰袋保存直至运送到实验室.

1.2 DNA提取

1 L水样经0.22 μm孔径的玻璃纤维滤膜过滤并收集滤膜, 沉积物称取5 mg湿样, 然后根据DNA提取试剂盒(Fast DNA Spin Kit for Soil, MP Biomedicals, 美国)说明书提取水样和沉积物样的总DNA.提取的DNA质量和完整性通过1.5%的琼脂糖凝胶电泳进行验证, 并通过超微量分光光度计(NanoDrop2000, 美国)进一步定量DNA的浓度.所有DNA样品存放于-20℃冰箱备用.

1.3 荧光定量PCR

本文通过荧光定量PCR技术(qPCR)定量分析了2类8种目标抗性基因和Ⅰ类整合子基因.本文使用的ARGs片段对应的引物序列、片段大小、退火温度和文献如表 1所示.荧光定量PCR反应体系(20 μL)为:2×Taq Pro Universal SYBR qPCR Master Mix 10 μL, 上下游引物各0.4 μL, DNA模板2 μL, ddH2 O 7.2 μL.基因扩增的热循环反应程序为:95℃预变性30 s; 40个循环(每个循环95℃ 10 s, 退火30 s, 退火温度见表 1); 溶解曲线温度设置为由60℃增加到95℃.

表 1 目的基因的qPCR引物信息和反应条件 Table 1 Primers, product lengths, and annealing temperatures of qPCR for target genes

为了生成用于测定提取DNA中每ng基因丰度的qPCR标准曲线, 将目标基因克隆到质粒(E. coli DH5a)中.使用NanoDrop2000测量质粒浓度, 并根据质粒长度和靶序列计算每毫升溶液中靶基因的丰度.通过对携带目的基因的质粒进行10倍的连续梯度稀释, 构建qPCR标准曲线.目的基因的标准曲线和扩增效率如表 2所示.

表 2 ARGs定量所需标准曲线和扩增效率1) Table 2 Standard curve and amplification efficiency of target genes

1.4 16S扩增子测序

为了研究长江口近岸地区水环境的细菌群落多样性, 使用上游引物341F(ACTCCTACGGGAGGC AGCA)和下游引物805R(GGACTACHVGGGTWTC TAAT)对16S rRNA基因V3/V4区域进行扩增.扩增产物送至上海伟寰生物科技有限公司, 采用Illumina NovaSeq PE250平台进行16S扩增子测序.

1.5 数据分析

本文使用Origin软件对所有数据进行汇总分析, 绘制柱状图、相关性热图和聚类分析热图, 用Spearman法进行相关性分析(P < 0.05表示显著相关), 使用R语言包“Hmisc”计算ARGs与微生物群落的相关系数.ARGs丰度变化的主成分分析采用免费的在线数据分析平台Tutools进行分析(https://www.cloudtutu.com), ARGs与微生物群落网络分析图的绘制采用生科云在线数据分析平台(https://www.bioincloud.tech/).

2 长江口近岸地区ARGs的赋存 2.1 ARGs的分布特征

图 2为长江口近岸地区水环境中ARGs相对丰度分布柱状图, W1~W8表示长江口近岸地区1~8号采样点的水样, S1~S8表示1~8号采样点的沉积物样.本文利用16S rRNA基因的丰度对其余目的基因的丰度进行标准化处理, 用来表征单一细菌ARGs的相对丰度(基因拷贝数与16S rRNA拷贝数的比值)[18, 25], 以减少DNA提取效率以及水环境中微生物背景值干扰所造成的误差[26, 27].本文选择的10种ARGs在长江口近岸地区8个采样点的水样和沉积物中的检出率是100%.

图 2 水样和沉积物样中ARGs的相对丰度 Fig. 2 Relative abundance of ARGs in water and sediment samples

图 2所示, 在水样中, 9种目标抗性基因中tetC的总相对丰度最高, 其相对丰度范围在5.01×10-2~1.76之间, 为长江口近岸地区表层水中的优势基因.在沉积物中, tetMtetA的总相对丰度处于较高的水平, 相对丰度范围分别在1.51×10-3~1.96×10-1及6.19×10-3~2.19×10-1之间, 为长江口近岸地区沉积物中的优势基因.不同采样点中ARGs的总相对丰度大小分布规律为:W2>W1>W4>W3>W5>W8>W7>W6(水样中); S1>S2>S4>S3>S5>S6>S8>S7(沉积物中).在水样中, W2和W1的ARGs相对丰度显著高于其他采样点; 在沉积物中, S1、S2和S4的ARGs相对丰度显著高于其他采样点.石洞口和竹园采样点位于大型污水处理厂及其邻近水域, 而浏河口采样点位于城市内河与长江下游交汇处.本文分析结果表明ARGs在长江口近岸地区水环境中的富集与污水处理厂的污水排放以及受人类活动影响的城市内河径流有关.该结论与前人研究该区域的结论相似[28, 29].

本文长江口近岸地区水样和沉积物中ARGs的绝对丰度平均值分别为2.32×104~7.66×104 copies ·mL-1和1.41×107~6.59×107 copies ·g-1. Lu等[30]研究了ARGs在渤海和黄海的空间分布, 结果表明在渤海, 目的基因的总绝对丰度范围在2.05×102~7.25×103 copies ·mL-1之间, 而在黄海中, 目的基因的总绝对丰度范围在21.1~8×103 copies ·mL-1之间.在渤海和黄海的河口区域, 目的基因的总绝对丰度范围分别在1.23×104~3.94×105copies ·mL-1和8.91×104~3.67×105copies ·mL-1之间.Chen等[31]对东海海湾沉积物中的27个目的基因进行了研究, 在杭州湾和象山湾河口区域, 目的基因的总绝对丰度范围分别在8.39×106~2.33×107copies ·g-1和1.25×107~2.08×107 copies ·g-1之间.Lu等[25]对辽河口和大辽河口磺胺类抗性基因(sul1sul2sul3)进行了表征, 辽河和大辽河口水样ARGs的绝对丰度平均值 <1×103 copies ·mL-1. 杨颖[32]调查了北江流域四环素抗性基因的污染情况, 发现在所检出的四环素抗性基因中, 含量最高的是tetC, 与16S rRNA拷贝数比值范围在8.3×10-2~13.2之间.范长征等[5]分析了湘江长沙河段水体及底泥中4种常见四环素类抗生素的含量特征及其季节变化, 其优势基因tetCtetA的相对丰度范围在8.2×10-5~7.3×10-3之间.长江口近岸地区ARGs的污染情况在全国处于较高的水平.

为了更深入地研究长江口近岸地区水环境中ARGs的分布特征以及水体和沉积物ARGs分布的相关性, 本文采用主成分分析法(PCA)[33]对长江口近岸地区8个采样点不同介质中ARGs的变化特征进行聚类分析.如图 3所示, X轴表示第一主成分(PC1), Y轴表示第二主成分(PC2).点与点之间的距离表示差异程度.PC1和PC2总共解释了基因丰度变异的85.2%, 其中, PC1的值为54.6%.可以得出, 表层水样品在PC1上聚集, 沉积物样品在PC2上聚集.说明不同采样点不同介质中的ARGs分别聚类在一起, 而且水样与沉积物样之间没有关联性.只有个别采样点的ARGs分布超出了95%置信椭圆的范围, 分别是水样中的W2和沉积物样中的S1.水样中ARGs的分布变化可能是由于石洞口污水处理厂的污水排放提供了大量的营养物质, 促进了当地的细菌生长.沉积物样中ARGs的分布变化可能是由于浏河口接纳了大量来自城市内河径流和人类活动引起的陆源排放的污染物, 包括重金属、抗生素、抗性细菌等. 浏河口和石洞口样点中ARGs的相对丰度平均值分别为:2.62×10-1和4.50×10-1(水样中), 5.27×10-2和4.32×10-2(沉积物中).这两个样点的ARGs污染水平处于同一数量级, 其差异很小, 而且这两个样点以及吴淞口和竹园排污口的污染水平显著高于其他4个样点, 所以污水处理厂的污水排放和城市内河径流带来的陆源污染物是ARGs污染的主要来源.Proia等[34]通过研究比利时Zenne河中ARGs的分布规律, 发现污水处理厂处理水的排放会增加ARGs的丰度.Guo等[35]研究了ARGs在长江口的赋存特征, 结果表明, 污水影响区域和内河交汇区域会对ARGs的丰度造成较大的影响.

图 3 不同介质中ARGs丰度的聚类分析 Fig. 3 Cluster analysis of ARGs abundance in different media

2.2 intI1整合子基因与ARGs的相关性

水平转移是ARGs在环境介质以及微生物中扩散传播的主要途径, 而整合子(integrons)是促进水环境中ARGs水平转移的重要移动基因元件[36~38].

长江口近岸地区8个采样点的intI1的相对丰度从高到低为:W2>W3>W1>W5>W8>W4>W7>W6(水样中); S3>S6>S1>S4>S5>S2>S8>S7(沉积物中).采用Spearman等级相关性检验[39]分析了intI1与其他ARGs的相关关系(如表 3), 说明intI1与多种ARGs都存在显著的相关关系, 尤其是在水样中. sul1基因在水样和沉积物中都与intI1呈显著相关性(P < 0.05), 这可能是由于sul1的基因序列通常与intI1的基因序列重叠所导致, 前人的研究也证实了这一点[40~43].然而, 在沉积物中, 四环素类抗性基因和intI1无显著相关性(P>0.05), 与前人的研究结论不符[44]. Yang等[44]以长江中下游的15个淡水湖泊为研究区域, 对18种不同的ARGs进行检测分析.本研究结果证明四环素类抗性基因tetGintI1之间存在显著的相关关系(P < 0.05).这说明不同种抗性基因的抗性机制不同, 导致其传播扩散的影响因子不同.虽然intI1是促进水环境中ARGs水平转移的重要移动基因元件, 但自然界中还广泛存在其他可移动遗传元件, 包括整合子、转座子、噬菌体、质粒和基因岛等[45, 46].而且, 其他的环境因素包括抗生素、重金属和多环芳烃等也可能影响表层沉积物中ARGs的传播扩散[47, 48].

表 3 表层水和沉积物中intI1与其余ARGs的相关性分析1) Table 3 Correlation analysis between intI1 and other ARGs in water and sediment samples

3 长江口近岸地区微生物群落分布 3.1 微生物群落多样性分析

本文采用16S扩增子测序技术对长江口近岸地区8个采样点的微生物群落多样性进行分析.细菌群落丰富度指数(Chao1)以及多样性指数(Shannon和Simpson)[49]计算公式如下:

(1)

式中, Sobs为观测到的OTU(操作分类单元)数, n1为只有一条序列的OTU数, n2为只有两条序列的OTU数.

(2)
(3)

式中, Sobs为观测到的OTU数, ni为第i个OTU中所含的序列数, N为所有序列数.

表 4所示, 沉积物的Chao1、Shannon和Simpson指数都高于表层水, 这说明长江口近岸地区沉积物的细菌丰富度和生物多样性要远高于表层水.竹园采样点表层水的Shannon指数在8个采样点中最高(6.80), 沉积物中的Shannon(7.36)指数仅次于吴淞口采样点(7.38), 这可能是因为污水处理厂的污水排放促进了细菌的繁殖, 增加了竹园周围水环境的微生物群落多样性.吴淞口采样点表层水中Simpson指数远远低于其他采样点, 这说明吴淞口水样中的群落物种多样性和均匀度偏低, 其他采样点的Simpson指数都比较接近, 表明物种多样性程度较均匀.

表 4 细菌群落丰富度指数和多样性指数 Table 4 Bacterial community richness and diversity index

3.2 微生物群落结构分析

本文对长江口近岸地区8个采样点的微生物群落结构特征进行分析.总体而言, 长江口水环境中的微生物群落由81个门、231个纲、569个目、961个科、1 960个属和4 788个种组成.在门分类水平上, 不同介质中微生物群落的相对丰度如图 4所示.

图 4 长江口近岸地区水环境中微生物群落在门水平上的相对丰度 Fig. 4 Relative abundance of microbial communities in the nearshore area of the Yangtze River Estuary at the phylum level

在水样中, 变形菌门(Proteobacteria)、放线菌门(Actinobacteriota)、拟杆菌门(Bacteroidota)和蓝菌门(Cyanobacteria)为主要菌群, 其总相对丰度平均值为85.87%, 其相对丰度平均值占比分别为46%、13.03%、13.00%和13.84%.经过石洞口污水处理厂处理水排放后的下游吴淞口采样点中变形菌门(Proteobacteria)占比显著提高, 占总比例的83.27%, 为绝对优势菌门, 污水处理厂对变形菌门(Proteobacteria)有着显著的影响[50].在沉积物中, 变形菌门(Proteobacteria)、拟杆菌门(Bacteroidota)、疣微菌门(Verrucomicrobiota)和酸杆菌门(Acidobacteriota)为主要菌群, 其总相对丰度平均值为64.06%, 其相对丰度平均值占比分别为36.03%、10.72%、9.23%和8.08%.并且, 变形菌门(Proteobacteria)和拟杆菌门(Bacteroidota)作为长江口水环境中的优势菌群, 可能参与众多物质的生物化学循环过程[51].

本文在属分类水平上对样品及其菌属进行了聚类分析, 并根据各样品中微生物菌属的相对丰度绘制热图, 如图 5所示.在水体中, W1和W5采样点之间的菌属差异较小, 其中相对丰度较高的菌属是ChloroplastFlavobacteriumRheinheimera.而W2和W3以及W4和W5采样点之间的菌属差异较大.W2采样点相对丰度较高的是RheinheimeraFlavobacterium, 分别占总量的31.62%和19.87%; W3采样点的主导菌属为Hydrogenophaga, 其相对丰度占总量的87.17%, 且显著高于其他采样点.在W4采样点相对丰度较高的是ChloroplastNitrospira, 分别占总量的14.84%和12.83%; 而W5采样点主导菌属是Chloroplast, 其相对丰度占总量的45.50%. W4和W5采样点之间的微生物菌属存在显著差异, 说明污水处理厂的污水排放对水环境中微生物群落结构变化有较大的影响.有研究发现, 污水处理厂可以直接将微生物释放到环境中[52, 53].在沉积物中, S2和S4采样点微生物分布的相似性较高, 差异较小, 其主导菌属均是硝化螺菌属(Nitrospira), 相对丰度分别占总量的12.46%和11.75%, 这可能是由于污水处理厂的厌氧生物处理工艺增加了水环境中污泥富集硝化螺菌属的风险.硝化螺菌属(Nitrospira)通常是污水处理厂活性污泥体系中的优势菌属[54].

图 5 长江口近岸地区属水平下的聚类分析热图 Fig. 5 Cluster analysis heat map of the nearshore area of the Yangtze River Estuary at the genus level

3.3 长江口近岸地区ARGs与微生物相关性分析

选取相对丰度前20的菌门与ARGs进行Speraman相关性分析, 相关性热图如图 6所示.20种菌门与10种ARGs之间有200种对应关系.在水样中, 菌门与部分ARGs间(29/200)呈显著相关(P < 0.05), 放线菌门(Actinobacteriota)、拟杆菌门(Bacteroidota)和蓝菌门(Cyanobacteria)与大部分ARGs呈正相关.硝化螺旋菌门(Nitrospirota)、髌骨门(Patescibacteria)、厚壁菌门(Firmicutes)、Latescibacterota、Bacteria、Desulfobacterota、MBNT15和NB1.j与ARGs呈显著负相关.在沉积物中, 菌门与部分ARGs间(25/200)呈显著相关(P < 0.05), 变形菌门(Proteobacteria)、Campilobacterota、Desulfobacterota、硝化螺旋菌门(Nitrospirota)、Spirochaetota、Sva0485和MBNT15与多种ARGs呈显著正相关关系.拟杆菌门(Bacteroidota)、浮霉菌门(Planctomycetota)、芽单胞菌门(Gemmatimonadota)和放线菌门(Actinobacteriota)与多种ARGs呈显著负相关关系.当微生物菌门与目的ARGs呈显著正相关关系时, 这些菌门可能是ARGs的潜在宿主.Wang等[55]研究了洪湖中微生物群落和6种目的ARGs之间的相关关系, 发现厚壁菌门(Firmicutes)、变形菌门(Proteobacteria)和芽单胞菌门(Gemmatimonadota)丰度与ARGs丰度呈显著正相关, 表明所属这些门的细菌可能是ARGs的宿主.微生物菌门与目的ARGs呈负相关关系可能意味着该菌门的丰度变化会影响ARGs的丰度和水平转移.Chen等[56]的研究表明, 硫杆菌属(Thiobacillus) 对土霉素具有一定的生物降解作用, 从而可能间接抑制ARGs的水平基因转移, 降低ARGs的相对丰度.

1. sul1 , 2. sul2, 3. tetM, 4. tetC, 5. tetX, 6. tetA, 7. tetO, 8. tetQ, 9. intI1, 10.16S rRNA;
*表示P < 0.05, **表示P < 0.01, ***表示P < 0.001
图 6 ARGs与微生物菌门相关性分析 Fig. 6 Correlation analysis between ARGs and microbial phylum

为了分析长江口近岸地区ARGs潜在的宿主信息, 本文采用网络共现图挖掘了ARGs与微生物群落的共现模式(图 7).基于Spearman相关性系数计算结果(r>0.7, P < 0.05), 选取ARGs与相对丰度排名前20的微生物菌门做网络共现分析.结果如图 7所示, 其中的连接线表示节点之间的相关系数大于0.7, 显著性小于0.05的相关关系, 红线表示正相关, 蓝线表示负相关.水样中共检测到5个潜在宿主, 沉积物中共检测到3个潜在宿主.在水样中, 磺胺类抗性基因sul1未检测到宿主信息, sul2只有1个宿主信息(放线菌门), 四环素类抗性基因(tetXtetAtetOtetQ)的潜在宿主较多.在水样中, 硝化螺旋菌门(Nitrospirota)是tetXtetAtetOtetQ共同的潜在宿主.微生物氮循环的频繁发生会产生硝化螺旋菌门(Nitrospirota)[8], 因此, 微生物氮的循环可能会间接影响水体中四环素类抗性基因的分布.在沉积物样品中, 只有sul1tetCintI1存在潜在宿主菌门.磺胺类抗性基因sul1的潜在宿主较多(2个), 说明在沉积物中sul1受微生物群落变化的影响较大. sul1的潜在宿主分别为Spirochaetota和Sva0485.硝化螺旋菌门(Nitrospirota)是tetC的潜在宿主.Sva0485是sul1intI1的共同潜在宿主.由于不同介质中微生物群落结构不同, ARGs的分布和诱导的潜在机制可能有所不同, 从而导致ARGs与微生物群落的相关性以及ARGs的潜在宿主菌存在差异.共同潜在宿主的存在可能会对ARGs的分布特征和水平转移产生影响, 会导致细菌获得多重耐药性, 在一定程度上加速细菌的产生, 生态安全和人类健康可能遭到威胁[57].长江口近岸地区多重耐药菌对ARGs的影响机制值得进一步研究和探索.

图 7 不同介质中ARGs与微生物群落网络分析 Fig. 7 Network analysis of the relationship among ARGs and microbial communities in different media

4 结论

(1) 2种磺胺类抗性基因(sul1sul2)、6种四环素类抗性基因(tetMtetCtetXtetAtetOtetQ)以及1种整合子基因intI1和16S rRNA基因在长江口近岸地区水环境中均有检出.城市内河径流和污水处理厂的污水排放对长江口近岸地区ARGs的分布有较大的影响.

(2) 在水体中, 整合子基因intI1与7种ARGs(sul1、sul2、tetMtetCtetXtetAtetO)存在显著正相关关系(P < 0.05); 在沉积物中整合子基因intI1只与一个磺胺类抗性基因sul1存在显著正相关关系.整合子基因intI1可能已经不是沉积物中调控ARGs迁移转化的重要影响因子.

(3) 微生物群落多样性分析上, 长江口近岸地区沉积物中细菌的丰富度和多样性要远高于表层水体.微生物群落结构上, 变形菌门(Proteobacteria)和拟杆菌门(Bacteroidota)是长江口近岸地区水环境中的优势菌门; Chloroplast为水体中的主要菌属, ChloroplastNitrospira为沉积物中的主要菌属.微生物群落多样性以及群落结构变化与污水处理厂的污水排放有关.

(4) 在水体中, 硝化螺旋菌门(Nitrospirota)是4种四环素类抗性基因(tetXtetAtetOtetQ)共同的潜在宿主; 在沉积物中, Sva0485是sul1intI1的共同潜在宿主.微生物群落的丰度变化对长江口近岸地区大部分ARGs的迁移转化有较大影响.

参考文献
[1] Pruden A, Pei R T, Storteboom H, et al. Antibiotic resistance genes as emerging contaminants: studies in northern Colorado[J]. Environmental Science & Technology, 2006, 40(23): 7445-7450.
[2] Alekshun M N, Levy S B. Molecular mechanisms of antibacterial multidrug resistance[J]. Cell, 2007, 128(6): 1037-1050. DOI:10.1016/j.cell.2007.03.004
[3] Aminov R I. The role of antibiotics and antibiotic resistance in nature[J]. Environmental Microbiology, 2009, 11(12): 2970-2988. DOI:10.1111/j.1462-2920.2009.01972.x
[4] 杨继平, 邱志刚, 袁兆康, 等. 天津海河流域抗生素抗性基因分布特征及与指示微生物的关系[J]. 环境与健康杂志, 2017, 34(4): 313-316.
Yang J P, Qiu Z G, Yuan Z K, et al. Distribution of antibiotic resistance genes and their relationship with indicator microorganisms in Haihe River basin, Tianjin[J]. Journal of Environment and Health, 2017, 34(4): 313-316.
[5] 范长征, 王聪, 鲁伦慧, 等. 湘江四环素及抗性基因含量特征及其季节变化[J]. 湖南大学学报(自然科学版), 2015, 42(6): 107-112.
Fan C Z, Wang C, Lu L H, et al. Occurrence and seasonal changes of tetracycline antibiotics and antibiotic resistance genes in the Xiang River[J]. Journal of Hunan University (Natural Sciences), 2015, 42(6): 107-112. DOI:10.3969/j.issn.1674-2974.2015.06.018
[6] 房平, 代鹤峰, 庄僖, 等. 东江下游典型饮用水源地抗生素抗性基因分布研究[J]. 生态环境学报, 2019, 28(3): 548-554.
Fang P, Dai H F, Zhuang X, et al. Distribution of antibiotic resistance genes in typical drinking water sources in the lower Dongjiang River[J]. Ecology and Environmental Sciences, 2019, 28(3): 548-554.
[7] 沈群辉, 冀秀玲, 傅淑珺, 等. 黄浦江水域抗生素及抗性基因污染初步研究[J]. 生态环境学报, 2012, 21(10): 1717-1723.
Shen Q H, Ji X L, Fu S J, et al. Preliminary studies on the pollution levels of antibiotic and antibiotic resistance genes in Huangpu River, China[J]. Ecology and Environmental Sciences, 2012, 21(10): 1717-1723. DOI:10.3969/j.issn.1674-5906.2012.10.015
[8] 郭行磐. 长江口滨岸水环境中抗生素抗性基因的赋存特征[D]. 上海: 华东师范大学, 2019.
Guo X P. Occurrence and abundance of antibiotic resistance genes in aquatic environments along the Yangtze estuary[D]. Shanghai: East China Normal University, 2019.
[9] Liu M, Hou L J, Yang Y, et al. Distribution and sources of polycyclic aromatic hydrocarbons in intertidal flat surface sediments from the Yangtze estuary, China[J]. Environmental Geology, 2001, 41(1-2): 90-95. DOI:10.1007/s002540100347
[10] Lin S, Hsieh I J, Huang K M, et al. Influence of the Yangtze River and grain size on the spatial variations of heavy metals and organic carbon in the East China Sea continental shelf sediments[J]. Chemical Geology, 2002, 182(2-4): 377-394. DOI:10.1016/S0009-2541(01)00331-X
[11] Hao S, Yi Y, Min L, et al. Occurrence and distribution of antibiotics in the surface sediments of the Yangtze Estuary and nearby coastal areas[J]. Marine Pollution Bulletin, 2014, 83(1): 317-323. DOI:10.1016/j.marpolbul.2014.04.034
[12] Xu T, Zhao W T, Guo X P, et al. Characteristics of antibiotics and antibiotic resistance genes in Qingcaosha Reservoir in Yangtze River Delta, China[J]. Environmental Sciences Europe, 2020. DOI:10.21203/rs.2.19989/v2
[13] Zhang J, Zhang Z F, Liu S M, et al. Human impacts on the large world rivers: would the Changjiang (Yangtze River) be an illustration?[J]. Global Biogeochemical Cycles, 1999, 13(4): 1099-1105. DOI:10.1029/1999GB900044
[14] Jiao N Z, Zhang Y, Zeng Y H, et al. Ecological anomalies in the East China Sea: impacts of the three gorges dam?[J]. Water Research, 2007, 41(6): 1287-1293. DOI:10.1016/j.watres.2006.11.053
[15] Jutkina J, Marathe N P, Flach C F, et al. Antibiotics and common antibacterial biocides stimulate horizontal transfer of resistance at low concentrations[J]. Science of the Total Environment, 2018, 616-617: 172-178. DOI:10.1016/j.scitotenv.2017.10.312
[16] Liu X H, Lu S Y, Guo W, et al. Antibiotics in the aquatic environments: a review of lakes, China[J]. Science of the Total Environment, 2018, 627: 1195-1208. DOI:10.1016/j.scitotenv.2018.01.271
[17] Pu C J, Yu Y, Diao J X, et al. Exploring the persistence and spreading of antibiotic resistance from manure to biocompost, soils and vegetables[J]. Science of the Total Environment, 2019, 688: 262-269. DOI:10.1016/j.scitotenv.2019.06.081
[18] Luo Y, Mao D Q, Rysz M, et al. Trends in antibiotic resistance genes occurrence in the Haihe River, China[J]. Environmental Science & Technology, 2010, 44(19): 7220-7225.
[19] Pei R T, Kim S C, Carlson K H, et al. Effect of River Landscape on the sediment concentrations of antibiotics and corresponding antibiotic resistance genes (ARG)[J]. Water Research, 2006, 40(12): 2427-2435. DOI:10.1016/j.watres.2006.04.017
[20] Ng L K, Martin I, Alfa M, et al. Multiplex PCR for the detection of tetracycline resistant genes[J]. Molecular and Cellular Probes, 2001, 15(4): 209-215. DOI:10.1006/mcpr.2001.0363
[21] Huang M H, Zhang W, Liu C, et al. Fate of trace tetracycline with resistant bacteria and resistance genes in an improved AAO wastewater treatment plant[J]. Process Safety and Environmental Protection, 2015, 93: 68-74. DOI:10.1016/j.psep.2014.04.004
[22] Schmitt H, Stoob K, Hamscher G, et al. Tetracyclines and tetracycline resistance in agricultural soils: microcosm and field studies[J]. Microbial Ecology, 2006, 51(3): 267-276. DOI:10.1007/s00248-006-9035-y
[23] Aminov R I, Chee-Sanford J C, Garrigues N, et al. Development, validation, and application of PCR primers for detection of tetracycline efflux genes of gram-negative bacteria[J]. Applied and Environmental Microbiology, 2002, 68(4): 1786-1793. DOI:10.1128/AEM.68.4.1786-1793.2002
[24] Koike S, Krapac I G, Oliver H D, et al. Monitoring and source tracking of tetracycline resistance genes in lagoons and groundwater adjacent to swine production facilities over a 3-year period[J]. Applied and Environmental Microbiology, 2007, 73(15): 4813-4823. DOI:10.1128/AEM.00665-07
[25] Lu Z H, Na G S, Hui G, et al. Fate of sulfonamide resistance genes in estuary environment and effect of anthropogenic activities[J]. Science of the Total Environment, 2015, 527-528: 429-438. DOI:10.1016/j.scitotenv.2015.04.101
[26] Guo X P, Li J, Yang F, et al. Prevalence of sulfonamide and tetracycline resistance genes in drinking water treatment plants in the Yangtze River Delta, China[J]. Science of the Total Environment, 2014, 493: 626-631. DOI:10.1016/j.scitotenv.2014.06.035
[27] 赵婉婷, 黄智峰, 郭雪萍, 等. 太湖周边饮用水处理厂中抗生素抗性基因污染分布特征[J]. 环境化学, 2020, 39(12): 3271-3278.
Zhao W T, Huang Z F, Guo X P, et al. Pollution and distribution characteristics of antibiotic resistance genes in drinking water treatment plants around Taihu Lake[J]. Environmental Chemistry, 2020, 39(12): 3271-3278.
[28] Guo X P, Liu X R, Niu Z S, et al. Seasonal and spatial distribution of antibiotic resistance genes in the sediments along the Yangtze Estuary, China[J]. Environmental Pollution, 2018, 242: 576-584. DOI:10.1016/j.envpol.2018.06.099
[29] Lin L, Yuan K, Liang X M, et al. Occurrences and distribution of sulfonamide and tetracycline resistance genes in the Yangtze River Estuary and nearby coastal area[J]. Marine Pollution Bulletin, 2015, 100(1): 304-310. DOI:10.1016/j.marpolbul.2015.08.036
[30] Lu J, Zhang Y X, Wu J, et al. Occurrence and spatial distribution of antibiotic resistance genes in the Bohai Sea and Yellow Sea areas, China[J]. Environmental Pollution, 2019, 252: 450-460. DOI:10.1016/j.envpol.2019.05.143
[31] Chen J Y, Su Z G, Dai T J, et al. Occurrence and distribution of antibiotic resistance genes in the sediments of the East China Sea bays[J]. Journal of Environmental Sciences, 2019, 81: 156-167. DOI:10.1016/j.jes.2019.01.016
[32] 杨颖. 北江水环境中抗生素抗性基因污染分析[D]. 广州: 中山大学, 2010.
Yang Y. Characterizing the pollution by the representative antibiotic resistance genes (ARGs) in the Beijiang River, South China[D]. Guangzhou: Sun Yat-Sen University, 2010.
[33] Guo X P, Zhao S, Chen Y R, et al. Antibiotic resistance genes in sediments of the Yangtze Estuary: from 2007 to 2019[J]. Science of the Total Environment, 2020, 744. DOI:10.1016/j.scitotenv.2020.140713
[34] Proia L, Anzil A, Subirats J, et al. Antibiotic resistance along an urban river impacted by treated wastewaters[J]. Science of the Total Environment, 2018, 628-629: 453-466. DOI:10.1016/j.scitotenv.2018.02.083
[35] Guo X P, Yang Y, Lu D P, et al. Biofilms as a sink for antibiotic resistance genes (ARGs) in the Yangtze Estuary[J]. Water Research, 2018, 129: 277-286. DOI:10.1016/j.watres.2017.11.029
[36] Stalder T, Barraud O, Casellas M, et al. Integron involvement in environmental spread of antibiotic resistance[J]. Frontiers in Microbiology, 2012, 3. DOI:10.3389/fmicb.2012.00119
[37] Stedtfeld R D, Stedtfeld T M, Waseem H, et al. Isothermal assay targeting class 1 integrase gene for environmental surveillance of antibiotic resistance markers[J]. Journal of Environmental Management, 2017, 198: 213-220. DOI:10.1016/j.jenvman.2017.04.079
[38] Stokes H W, Hall R M. A novel family of potentially mobile DNA elements encoding site-specific gene-integration functions: integrons[J]. Molecular Microbiology, 1989, 3(12): 1669-1683. DOI:10.1111/j.1365-2958.1989.tb00153.x
[39] 程森, 路平, 冯启言. 渔业复垦塌陷地抗生素抗性基因与微生物群落[J]. 环境科学, 2021, 42(5): 2541-2549.
Cheng S, Lu P, Feng Q Y. Distribution of antibiotic resistance genes and microbial communities in a fishery reclamation mining subsidence area[J]. Environmental Science, 2021, 42(5): 2541-2549.
[40] Bennett P M. Plasmid encoded antibiotic resistance: acquisition and transfer of antibiotic resistance genes in bacteria[J]. British Journal of Pharmacology, 2008, 153(S1): S347-S357. DOI:10.1038/sj.bjp.0707607
[41] Chee-Sanford J C, Mackie R I, Koike S, et al. Fate and transport of antibiotic residues and antibiotic resistance genes following land application of manure waste[J]. Journal of Environmental Quality, 2009, 38(3): 1086-1108. DOI:10.2134/jeq2008.0128
[42] Nigro S J, Farrugia D N, Paulsen I T, et al. A novel family of genomic resistance islands, AbGRI2, contributing to aminoglycoside resistance in Acinetobacter baumannii isolates belonging to global clone 2[J]. Journal of Antimicrobial Chemotherapy, 2013, 68(3): 554-557. DOI:10.1093/jac/dks459
[43] 江月, 陈奕涵, 何义亮. 长江下游某水源型水库抗生素抗性基因污染研究[J]. 环境科学学报, 2018, 38(3): 884-891.
Jiang Y, Chen Y H, He Y L. Study on antibiotic resistance genes contamination in a reservoir located at Yangtze River downstream[J]. Acta Scientiae Circumstantiae, 2018, 38(3): 884-891.
[44] Yang Y Y, Liu W Z, Xu C, et al. Antibiotic resistance genes in lakes from middle and lower reaches of the Yangtze River, China: effect of land use and sediment characteristics[J]. Chemosphere, 2017, 178: 19-25. DOI:10.1016/j.chemosphere.2017.03.041
[45] van Elsas J D, Bailey M J. The ecology of transfer of mobile genetic elements[J]. FEMS Microbiology Ecology, 2002, 42(2): 187-197. DOI:10.1111/j.1574-6941.2002.tb01008.x
[46] Zhang M, He L Y, Liu Y S, et al. Variation of antibiotic resistome during commercial livestock manure composting[J]. Environment International, 2020, 136. DOI:10.1016/j.envint.2020.105458
[47] Nguyen C C, Hugie C N, Kile M L, et al. Association between heavy metals and antibiotic-resistant human pathogens in environmental reservoirs: a review[J]. Frontiers of Environmental Science & Engineering, 2019, 13(3). DOI:10.1007/s11783-019-1129-0
[48] Das N, Kotoky R, Maurya A P, et al. Paradigm shift in antibiotic-resistome of petroleum hydrocarbon contaminated soil[J]. Science of the Total Environment, 2021, 757. DOI:10.1016/j.scitotenv.2020.143777
[49] Finotello F, Mastrorilli E, Di Camillo B. Measuring the diversity of the human microbiota with targeted next-generation sequencing[J]. Briefings in Bioinformatics, 2018, 19(4): 679-692.
[50] 罗晓, 张文丽, 袁立霞, 等. 纳污河流抗性基因和微生物群落相关性[J]. 中国环境科学, 2019, 39(6): 2606-2613.
Luo X, Zhang W L, Yuan L X, et al. Correlation between resistance genes and microbial community in polluted rivers[J]. China Environmental Science, 2019, 39(6): 2606-2613. DOI:10.3969/j.issn.1000-6923.2019.06.044
[51] Guo X P, Lu D P, Niu Z S, et al. Bacterial community structure in response to environmental impacts in the intertidal sediments along the Yangtze Estuary, China[J]. Marine Pollution Bulletin, 2018, 126: 141-149.
[52] Kim S, Aga D S. Potential ecological and human health impacts of antibiotics and antibiotic-resistant bacteria from wastewater treatment plants[J]. Journal of Toxicology and Environmental Health, Part B, 2007, 10(8): 559-573.
[53] Yu M D, Liu S J, Li G W, et al. Municipal wastewater effluent influences dissolved organic matter quality and microbial community composition in an urbanized stream[J]. Science of the Total Environment, 2020, 705. DOI:10.1016/j.scitotenv.2019.135952
[54] Yao Q, Peng D C. Nitrite oxidizing bacteria(NOB) dominating in nitrifying community in full-scale biological nutrient removal wastewater treatment plants[J]. AMB Express, 2017, 7(1). DOI:10.1186/s13568-017-0328-y
[55] Wang Z, Han M Z, Li E H, et al. Distribution of antibiotic resistance genes in an agriculturally disturbed lake in China: their links with microbial communities, antibiotics, and water quality[J]. Journal of Hazardous Materials, 2020, 393. DOI:10.1016/j.jhazmat.2020.122426
[56] Chen J F, Yang Y W, Liu Y Y, et al. Bacterial community shift and antibiotics resistant genes analysis in response to biodegradation of oxytetracycline in dual graphene modified bioelectrode microbial fuel cell[J]. Bioresource Technology, 2019, 276: 236-243.
[57] 王龙飞, 程逸群, 胡晓东, 等. 江苏省代表性水源地抗生素及抗性基因赋存现状[J]. 环境科学, 2021, 42(2): 749-760.
Wang L F, Cheng Y Q, Hu X D, et al. Occurrence of antibiotics and antibiotic resistance genes in representative drinking water resources in Jiangsu province[J]. Environmental Science, 2021, 42(2): 749-760.