地埋式一体化污水处理设备工艺

抗生素的滥用诱导和加速了抗生素抗性基因(antibiotic resistance genes, ARGs)的产生、传播以及抗生素耐药菌(antibiotic resistant bacteria, ARB)形成的风险, 并且ARGs的存在是细菌表现耐药性的根本原因.ARGs作为一类新型的环境污染物, 可通过质粒、整合子以及转座子等可在细菌同种属菌株间和不同种属的菌株之间发生水平基因转移, 即使ARB死亡后, 携带ARGs的裸露DNA会在脱氧核苷酸酶的保护作用下长期存在, 故ARGs在环境中的持久性残留、传播和扩散比抗生素本身的危害还要大.目前, ARGs广泛存在于不同的生态环境中, 如水环境、土壤和大气, 其中作为水环境之一的污水厂被认为是ARB和ARGs巨大的储存地.以前针对ARGs的大多数研究主要集中在处理医院和生活污水的城市污水厂上.然而, 迄今为止, 很少有报道从医药化工废水中检测出ARB和ARGs.

 ARGs的检测方法有可培方法和非培方法这两大类.可培方法是通过培养基进行微生物分离, 而后根据所分离微生物的耐药性, 结合各种分子生物学方法或者对ARB进行全基因组.但是由于自然环境中**过99%的微生物不能用传统培养基方法进行分离培养, 故可培方法会遗漏大量ARB和ARGs.非培方法有PCR技术、qPCR技术、基因芯片技术及宏基因组学技术.无论是普通PCR技术还是qPCR技术只能局限于已知序列的单一耐药基因检测, 并且只能对少数的ARGs进行同时检测, 最后还需要测序分析来进一步确认.基因芯片技术虽然可以同时检测出大量的ARGs, 但仍然要依赖于根据已知的耐药基因序列设计出相应的特异性强的探针, 并且其制作复杂, 费用昂贵.而随着生物技术的不断发展, 宏基因组测序也越来越多地应用于医学研究及临床诊断领域, 如感染类型诊断、抗性基因的鉴定和传染病的防控等.宏基因组技术通过高通量测序, 打破了传统以培养为主的微生物研究方式, 直接提取遗传物质进行遗传操作, 不仅可以准确获得样品微生物物种组成及丰度等信息, 进入数据库预测微生物功能, 还能分析微生物代谢网络等多方面的信息, 并且可同时检测大量的ARGs, 为寻找新基因、研究微生物多样提供了便利的工具.

本研究借助宏基因组技术调查医药化工废水厂中微生物群落组成, 并分析ARGs的多样性及两者之间的关系, 对掌握医药化工废水中ARGs的污染现状及控制其传播, 保护微生态环境具有重要指导意义.

1 材料与方法

1.1 样品采集

本实验样品取自地埋式一体化污水处理设备工艺, 该废水厂主要功能是处理园区各医药化工企业排放的工业废水, 纳入废水厂管网的废水包括抗生素中间体废水及农药中间体废水等多种医药化工废水.该废水厂日处理量10万t, 收集的废水中含β-内酰胺类、氨基糖苷类、大环内酯类、喹诺酮类和四环素类等几大类抗生素生产过程中产生的高浓发酵废水、化学合成废水、其他制剂废水及混合废水(清洗水、冷却循环水和生活污水等).废水中含大量残留的抗生素原材料、副反应产物, 还含**溶剂二甲基甲酰胺、丙酮、乙腈、四氢呋喃和二等.废水厂主要采用厌氧-缺氧-好氧(A2O)工艺.本实验样品为活性污泥, 具体取自完全混合式好氧生物处理池.取样时先在池内前中后三段各取约500 mL泥水混合样品, 以保证样品具有代表性, 现场混合均匀后, 用无菌50 mL离心管分装后冰袋低温冷藏, 当日即进行基因组DNA提取.

2 结果与讨论

2.1 地埋式一体化污水处理设备工艺中微生物群落结构

在该好氧活性污泥样品中, 经过质量控制后, 共得到reads数目109 523 076条, 并由此组装得到55 255个scaftig, 并预测了105 892个unigene.利用DIAMOND软件将unigene在NCBI-NR数据库中进行比对, 后续用MEGAN软件进行注释, 统计样本在门水平及属水平上的前10个丰度较大的群落组成.

2.2 地埋式一体化污水处理设备工艺中微生物群落功能注释

用DIAMOND软件将unigene在KEGG数据库进行比对, 共注释了54564个unigene, 结果从中可知, 该微生物菌群的功能通路相关的基因数量依次为：代谢>遗传信息处理>环境信息处理>人类疾病>细胞过程>生物系统, 其中较主要的功能通路是代谢相关, 其次是遗传信息处理, *三是环境信息处理.属于代谢类别的功能分类中丰度依次为：碳代谢>氨基酸代谢>能量代谢>核苷酸代谢>辅因子和维生素的代谢>异生素生物降解与新陈代谢>脂质代谢>萜类和聚酮化合物代谢>其他氨基酸代谢>其他次生代谢物生物合成>多糖生物合成与代谢, 其中碳代谢、氨基酸代谢和能量代谢的基因数量尤为丰富, 分别为7 788、6 155和5 067个.属于遗传信息处理类别的功能分类丰度依次为：复制与修复>折叠、分类与降解>翻译>转录, 其中复制与修复是较主要部分, 基因数量为2 535个.该微生物菌群功能通路中还有一部分的环境信息处理类别的功能, 其分类类别丰度依次为：膜转运>信号转导>信号分子与相互作用, 其中较主要的是膜转运, 基因数据为2 678个, 其他两个类别的基因合计100个, 数目甚少.属于人类疾病类别丰度依次为：传染性疾病>神经退行性疾病>内分泌和代谢疾病>癌症>物质依赖>*疾病>心血管疾病, 其中细菌传染性疾病相关的基因数目尤其**, 为517个, 是其他人类疾病基因总数的约1.6倍.而在传染病疾病中, 细菌感染性疾病相关的基因占据了大部分, 约80%.

在该废水厂中细菌传染疾病涉及到9类, 其丰度依次为军团杆菌病>结核病>沙门氏菌感染>幽门螺杆菌感染的上皮细胞信号传导>金黄色葡萄球菌感染>上皮细胞的细菌入侵>百日咳>致病性大肠杆菌感染>霍乱弧菌感染, 其中较主要的传染病为军团杆菌病, 其基因数目为136个;其次是结核病, 基因数目为108个;沙门氏菌感染和幽门螺杆菌感染的上皮细胞信号传导, 基因数目均在60个左右;其他疾病的基因数目在15个以下.

根据功能注释结果, 该废水厂菌群中较主要的功能通路是代谢相关, 并且还存在许多与人类疾病相关的基因, 其中主要是细菌感染性疾病.到目前为止, **新出现的传染病事件的发生率有所增加, 废水厂等特殊环境可以作为许多病原体的**宿主对人类健康存在潜在威胁.

2.3 地埋式一体化污水处理设备工艺中ARGs及耐药机制

用BLAST软件将unigene在抗性基因数据库CARD中进行比对, 结果.经过氨基酸大于50%标准控制后, 共注释到74类抗生素功能基因, 涉及包括β-内酰胺类、头孢菌素类、氨基糖苷类、大环内酯类、多肽类、喹诺酮类及磺胺类等多种抗生素抗性基因.

在这74类ARGs中sav1866、dfrE和mfd是较主要ARGs.sav1866基因属于多重耐药基因, 其编码的蛋白是ABC膜转运蛋白, 其耐药机制属于ABC外排泵.在细菌中, 许多ABC家族转运蛋白都与抗生素耐受性有关.sav1866被检测为丰度较大的ARGs, 推测这与该废水厂处理多种抗生素中间体废水及多种农药中间体废水密切相关.dfrE基因**被发现在粪肠球菌中, 其编码的蛋白是一种耐甲氧苄啶二氢叶酸还原酶(dihydrofolate reductases, DHFRs), 该种DHFRs具有高的甲氧苄啶抗性, 其耐药机制是抗生素靶位点改变.mfd基因的产物Mfd是转录-修复偶联因子蛋白, 属于DNA移位酶, 也是DNA修复的必须因子, 其耐药机制是抗生素靶位点保护.Mfd主要功能是促进DNA修复, 但在某些情况下, 如在复制-转录冲突区域和在稳定阶段的突变中, Mfd反而增加了突变, 如Mfd促使空肠弯曲杆菌产生了对喹诺酮药物的抗性.在新的研究结果中, Mfd被认为是一种促进突变的进化因子, 加速了抗生素耐药性的发展, 是不同类群的细菌快速获得各种抗生素耐药性所必须的[35].由此推测, mfd基因的大量存在及积累可能会促进废水厂中微生物产生更多抗生素抗性.

统计这些ARGs的耐药机制, 以调查该废水厂的耐药机制情况, 在该医药化工废水好氧活性污泥处理系统中, 微生物存在耐药的机制有5种, 分别是抗生素失活(8%)、抗生素靶位点改变(4%)、抗生素靶位点保护(7%)、抗生素靶位点替换(11%)和抗生素特异的外排泵(70%), 其中抗生素特异的外排泵是主要的耐药机制, 占比50%以上.在所研究的菌群中的外排泵组成， 主要由ABC外排泵(18%)、MFS外排泵(10%)、RND外排泵(34%)和SMR外排泵(3%)组成, 其中RND外排泵是该微生物菌群抗生素外排泵耐药机制的较主要类型, 其次是ABC外排泵.据报道, 细菌细胞膜上的主动外排泵可将抗菌药物排至细菌外, 从而阻碍抗菌药物与细菌内的靶点结合, 引起耐药.并且, 外排泵在细菌中广泛存在, 在介导细菌多重耐药中起着重要作用.本研究表明, 氯霉素、四环素、磺胺类、氟喹诺酮类和大环内酯类抗生素均能通过细菌的主动外排机制而被排出细胞外.而该研究的医药化工污水系统所具有的外排泵中, RND外排泵是主要机制.RND外排泵的外排底物比较广泛, 包括多种抗菌药物、消毒剂和染料等, 由此可推测该菌群对多种抗生素具有抗性. 属于RND外排泵的ARGs有smeD、mdtB、mexJ、smeR、cmeB和mexQ等, 涉及氨基抗生素、四环素抗生素、大环内酯类、氨基糖苷抗生素、三氯沙抗生素和氟喹诺酮类抗生素、头孢菌素等ARGs;属于ABC外排泵的ARGs有lmrD、otrC、novA、macB和sav1866, 涉及林可酰胺类抗生素、四环素类抗生素和大环内酯类抗生素的耐药.

2.4 废水厂中ARGs与微生物类群的关联

对废水厂中ARGs与微生物类群的共存模式进行网络分析, 结果显示有44类ARGs与污水厂中属级的微生物分类群高度相关, 并且主要类群的属都是多重耐药菌株, 同时携带多种ARGs, 涉及到4大类ARGs耐药机制.Hyphomicrobium属和Thauera属均可能携带有近20个ARGs, Paracoccus属和Methyloversatilis属可能携带有近10个ARGs.如Hyphomicrobium属中可能携带的抗性基因有dfrA26、mdsB、qepA、tetT、mexQ、mexF、mexL、rosB、sul2、dfrE、mfd和bacA等;Thauera属则可能携带有mdtD、bacA、mexJ、tet36、oprM、acrB和novA等.

(1) 该地埋式一体化污水处理设备工艺的菌群组成结构中, 细菌域是主要微生物类类群, 并且较优势的菌门是Proteobacteria, 其次是Chloroflexi.较优势属是Hyphomicrobium, 较优势种为Hyphomicrobium zavarzinii.

(2) 该医药化工废水厂中微生物群落较主要的功能通路是代谢相关, 并还存在许多与人类疾病相关的基因, 其中主要是细菌感染性疾病.

(3) 该医药化工废水厂中ARGs共检测到74类ARGs, 涉及包括β-内酰胺类、头孢菌素类、氨基糖苷类、大环内酯类、多肽类、喹诺酮类、磺胺类等多种抗生素抗性基因, 其中sav1866、dfrE和mfd是较主要ARGs, 而较主要的Hyphomicrobium属也同时携带这3种主要的ARGs.该废水厂菌群的耐药机制主要是抗生素特异的外排泵, 其中RND外排泵是主要机制.

(4) 该废水厂活性污泥中微生物属级分类与ARGs高度相关, 并且主要类群的属都是多重耐药菌株.