CIP在环境中持续地迁移扩散, 当达到一定浓度后会对某些敏感微生物的生长产生抑制, 甚至直接杀死它们, 从而改变土著微生物区系和数量, 影响微生物群落结构分布.同时, 环境中残留的CIP等氟喹诺酮类抗生素能促进抗性基因的产生.目前, 在河流、湖泊、土壤、底泥等不同介质中就检出了由于抗生素作用而产生的抗性基因和耐药细菌.抗性基因的传播和扩散可能会加快抗药性菌群的大量繁殖, 且细菌抗药性一旦生成就很难在短时间内消失， 对微生物群落结构形成了潜在的威胁.抗性基因作为一种新型污染物, 其危害比抗生素本身要大, 这给人类的身体健康和生态环境安全带来了不可预测的风险.近年来, 围绕着CIP对水或土壤中微生物群落的影响开展了较多的研究, 但大多集中于CIP痕量水平, 关于高浓度CIP对活性污泥中微生物群落结构、功能微生物数量及抗性基因丰度的影响还少有报道.
Xie等针对某抗生素生产废水开展研究, 发现废水中CIP浓度很高, 达到5.06 mg ?L-1, 生化处理对CIP去除效果不佳, 处理出水中残留浓度仍高达3.23 mg ?L-1; 且该系统对有机物和氨氮去除效果较差, 推测是废水中的CIP对活性污泥中功能微生物造成了影响.此外, 本实验室前期研究也发现CIP对膜生物反应器运行效能会产生负面影响.
本研究旨在通过MBR处理CIP模拟废水实验, 考察不同CIP投加浓度下MBR中微生物群落结构和功能微生物丰度的变化; 探讨CIP-ARGs, 包括(bāo kuò)gyr
A、gyr
B、parC和parE在CIP选择压力下的变化情况, 以期为CIP生产废水的无害化处理提供数据支撑.
1 材料与方法1.1 实验用水
MBR进水中含有200 mg ?L-1葡萄糖、256.67 mg ?L-1无水乙酸钠(Sodium)、94.17 mg ?L-1硫酸铵和17.50 mg ?L-1磷酸二氢钾, 其CO
D、氨氮和磷(P)浓度分别约为400、20和4 mg ?L-1.整个实验分为4个工况, 共运行132 d, 运行工况1、2、3、4条件下分别向上述MBR进水中投加CIP浓度为0、5、10和15 mg ?L-1.所用CI
P、葡萄糖、无水乙酸钠、硫酸铵、磷酸二氢钾等药剂均为分析纯.
1.2 实验装置与运行条件
MBR主要由有机玻璃容器、膜组件和曝气装置组成.反应器中接种污泥取自某生活污水处理厂好氧膜池, 经稀释水洗和稀释后, 调整反应器中初始MLSS为6.2 g ?L-1.整个MBR装置和进水桶采取遮光处理, 避免CIP发生光解.反应器连续进水, 流速为25 mL ?min-1; 间歇抽吸出水, 每连续出水5 min后停止出水1 min, 以减缓膜污染, 出水流速为30 mL ?min-1.整个实验中, 反应器水力停留时间为8 h, 溶解氧为2~4 mg ?L-1, pH为7~8, 水温为25~30℃, COD容积负荷约1.2 kg ?-1, 除工况1视情况(Condition)周期性排泥外其余3个工况几乎没有排泥.
图 1
图 1 实验装置结构示意
每个工况下周期性取样分析常规水质指标和CIP的变化情况, 并取污泥样品密封避光保存于-70℃超低温冰箱内, 用于分析微生物群落及其抗性基因的变化.
1.3 分析项目与方法1.3.1 高通量测序
分别于第33
D、36 d和66
D、69 d和99
D、102 d和132 d采取污泥样品, 委托生工生物工程股份有限公司进行高通量测序.
泥水混合物3 000 r ?min-1离心去上清液, 沉淀（precipitation）固体部分使用E.Z.N.ATM Mag-Bind Soil DNA Kit试剂盒进行DNA提取, 之后用Qubit®2.0荧光计检测DNA浓度, 用琼脂糖凝胶检测DNA完整性, 用Qubit 2.0 DNA检测试剂盒对基因组DNA精确定量.PCR扩增采用融合了MiSeq测序平台的V3-V4通用引物341F CCTACGGGNGGCWGCAG)和805R.
1.3.2 抗性基因检测分析
采用PowerSoil® DNA Isolation Kit试剂盒并结合液氮冷冻研磨法对样品的基因组DNA进行提取.PCR的扩增引物如表 1所示.PCR反应体系:DNA模板1 μL, 2×Taq酶12.5 μL, 上下游引物各0.5 μL, 10.5 μL ddH2O, 反应体系总体积为25 μL.PCR反应程序:94℃预热5 min, 30个循环:94℃变性30 s, 退火温度30 s, 72℃延伸30 s, 最后72℃停留7 min.再采用Universal DNA Purification Kit试剂盒对CIP-ARGs的PCR产物进行纯化回收, 并委托中国科学院南京土壤研究所对纯化的PCR产物进行定量分析.
表 1 引物信息 Table 1 Primer information
1.3.3 其它指标（target aim)检测
CIP浓度采用高效液相色谱仪HPLC测定:色谱柱为ODS-2, 5 μm 4.6×250 mm; 检测器为紫外可见吸收检测器, 波长为277 nm; 流动相为色谱纯乙腈:水=20 :80; 流速为0.7 mL ?min-1; 温度为35℃; 进样体积为10 μL.CO
D、氨氮和MLSS分别采用重铬酸钾法、纳氏试剂分光光度法和重量法[30]测定; TOC采用TOC-VCSN总有机碳分析仪测定.
2 结果与讨论2.1 不同CIP投加浓度下MBR中微生物群落解析2.1.1 活性污泥的微生物群落结构
不同CIP投加浓度下MBR污泥中微生物群落在门、纲、目、科、属等分类水平上的平均种类分布如图 2所示.门、纲、目、科、属等5个分类水平微生物种类数分布大体上都随着CIP投加浓度和运行时间的增加呈现减少的趋势.
图 2
图 2 各分类水平微生物群落的数量分布
在工况1下MBR运行稳定时, 微生物群落分析中共检出了17个门、34个纲、60个目、104个科和232个属.
当运行至工况2, MBR进水中CIP浓度为5 mg ?L-1时, 微生物体系突然受到CIP的胁迫, 敏感微生物活性受到抑制而死亡, 使得微生物群落的门、纲和目种类减少; 继续运行至该阶段末, MBR运行逐渐趋于稳定, 此时科和属种类数相较于0 mg ?L-1时分别下降了21%和25%.
继续增加MBR进水中CIP浓度至10 mg ?L-1和15 mg ?L-1, 门和纲水平上微生物种类数基本趋于稳定; 但目、科、属水平上微生物种类数却有较大程度减少, 尤其是属水平, 工况4末的属水平种类数相较于工况1末下降了53%.结果显示, CIP的加入对微生物群落的门和纲数量影响不大, 但对目、科、属数量影响较大; 表明了微生物分类单位越小, CIP胁迫所引起的微生物种类变化越明显.
为了深入探究不同CIP浓度对微生物群落结构的影响, 接下来分别从门、科、属水平上进行群落结构分析.门水平上, MBR中微生物群落结构分布如图 3所示.
图 3
图 3 门水平上微生物群落结构分布
MBR在工况1中运行稳定时, 微生物群落共检出17个门, 优势菌群为变形菌门、拟杆菌门、浮霉菌门、酸杆菌门和硝化螺旋菌门, 其中变形菌门所占丰度比例最大, 为52%.工况2, MBR进水中开始投加5 mg ?L-1的CIP, 上述优势菌群受到CIP的突然冲击, 其相对丰度比例呈现出不同程度地降低; 厚壁菌门、放线菌门、Candidatus Saccharibacteria和绿弯菌门为该工况下新增的优势菌群.
随着运行时间的延长和CIP投加浓度的继续增加, 上述优势菌群的相对丰度比例呈现不同的变化趋势:变形菌门和拟杆菌门的相对丰度比例呈先升高后降低再升高再降低的趋势, 表明这两种菌群对CIP浓度的增加比较敏感; 浮霉菌门、酸杆菌门和硝化螺旋菌门的相对丰度比例均呈先增后减的趋势, 表明这3种菌群能适应一定浓度的CIP环境但在更高浓度的CIP环境中很难生存; 5 mg ?L-1时新检出的优势菌群均不能适应CIP的长期胁迫而持续减少, 甚至消失.
在高浓度CIP的长期选择性压力下变形菌门和拟杆菌门成为MBR污泥中的主要菌群, 工况4末二者相对丰度比例分别为57.5%和12.7%, 表明MBR中大部分抗CIP的微生物是属于这两种优势菌群, 这与前人的研究一致.
科水平上, MBR中微生物群落结构分布如图 4所示.
图 4
图 4 科水平上微生物群落结构分布
工况1条件下, MBR中检出的优势菌科有红环菌科、Kofleriacea
E、生丝微菌科、Chitinophagaceae和浮霉菌科, 所占比例分别为15.45%、5.80%、5.67%、16.10%和10.13%;前3者均属于变形菌门, 后两者分别隶属于拟杆菌门和浮霉菌门.工况2, 当MBR进水中开始投加CIP时, 这些优势菌科活性和生长受到CIP的抑制, 相对丰度比例均急剧降低; 而红杆菌科、消化链球菌科和叶杆菌科比例增大, 成为优势菌科.
当MBR在5~15 mg ?L-1的CIP环境中运行期间, 微生物群落结构产生不同的变化.其中红环菌科逐渐适应CIP环境, 相对丰度比例持续增加; Kofleriacea
E、Chitinophagaceae和浮霉菌科的相对丰度比例呈现出先升高后降低再升高再降低的过程, 表现出波动状态, 表明这3种菌科易受到CIP浓度增加的影响, 其中Chitinophagaceae的变化趋势与拟杆菌门完全一致; 生丝微菌科、红杆菌科、消化链球菌科和叶杆菌科相对丰度比例逐渐降低, 表明这4种菌科经不住CIP的长期胁迫.
面对高浓度CIP长时间的选择性压力, 在工况4末期选择出来的优势菌科为红环菌科、Chitinophagaceae和丛毛单胞菌科, 相对丰度比例分别为29.96%、5.44%和6.60%.
属水平上, MBR中微生物群落结构分布如图 5所示.MBR进水中未投加CIP时, 污泥中的优势菌（fungus)属为Azospir
A、Kofleri
A、Ferruginibacte
R、动胶菌属和生丝微菌属, 所占比例分别为13.07%、5.80%、10.41%、7.03%和4.48%.工况2, 当MBR进水中开始投加CIP时, 上述优势菌属相对丰度比例均急剧降低, 表明其活性和生长受到了CIP的抑制; 此时优势菌属为动胶菌属、Tissierell
A、Acetoanaerobium和Roseicitreum.
图 5
图 5 属水平上微生物群落结构分布
由图 5可知, 继续延长CIP作用时间并增加CIP的投加浓度, 微生物群落结构分布产生显著（striking）性差异.其中Azospir
A、Kofleri
A、Tissierell
A、生丝微菌属、Acetoanaerobium和Roseicitreum在CIP长期胁迫下很难生存而失去优势; 而Methyloversatili
S、Ferruginibacte
R、动胶菌属和丛毛单胞菌属在高浓度CIP的选择性压力下成为优势菌属, 工况4末各相对丰度比例分别为21.70%、7.56%、5.24%和4.15%.
门、科、属水平上的微生物群落结构分析(Analyse)表明, 在高浓度CIP的长期胁迫和选择性压力下, 敏感性微生物活性受到抑制而难以生存; 而变形菌门中红环菌科下的Methyloversatili
S、Ferruginibacter和动胶菌（fungus)属及丛毛单胞菌科下的丛毛单胞菌属最终成为了MBR中的优势菌属.
2.1.2 CIP投加浓度对MBR中微生物群落多样性的影响
考察CIP投加浓度对MBR中微生物群落丰富度与多样性的影响, 结果如表 2所示.Chao1和ACE指数用来估计微生物群落的丰富度, 而Simpson和Shannon指数代表群落多样性.由表 2可知, 随着MBR进水中CIP投加浓度的增加和反应时间的延长, Chao1指数呈下降趋势, 从0 mg ?L-1的1271.7降至15 mg ?L-1的469.7, 下降率为63%;ACE指数在整个运行中比较波动, 但整体上呈下降趋势; Simpson指数持续升高, 而Shannon指数持续降低.Chao1、AC
E、Shannon指数的降低和Simpson指数的升高表明了MBR污泥中微生物群落受到了CIP的抑制导致部分微生物难以生存而死亡, 使得污泥中微生物群落丰富度和多样性均呈下降趋势, 这与图 2显示的结果一致.
表 2 CIP投加浓度对微生物(Micro-Organism)群落Chao1、AC
E、Simpson和Shannon指数的影响
2.1.3 CIP投加浓度对部分功能微生物群落的影响
考察CIP投加浓度对TO
C、氨氮和CIP去除效果的影响, 浅析污染物去除率变化与微生物群落结构变化的关系.结合图 5和图 6可知, 随着CIP投加浓度从0 mg ?L-1增加到15 mg ?L-1, TOC平均去除率从97.56%下降至90.84%, 表明CIP对有机物去除有所影响但影响程度不大.众所周知, 菌胶团是活性污泥和生物膜的重要组成部分, 是有机物降解去除的主力军; 动胶菌属作为菌胶团的主要构成体, 在CIP长期胁迫下, 其相对丰度比例趋于稳定, 或许TOC去除效果较稳定就与之相关.
图 6
图 6 CIP投加浓度对污染物去除率的影响
图 5显示, 丛毛单胞菌属也逐渐成为MBR中优势菌属之一, 有研究报道该菌属能利用CIP作为碳源[32]; 但在整个运行期间CIP去除率较低, 除了5 mg ?L-1时主要靠污泥吸附去除59.03%的CIP外, 其余时段平均去除率仅为7.24%~24.90%, 表明大部分优势菌种虽可在CIP胁迫下生存和繁殖, 但对CIP利用率不高.
由图 6和图 7可知, 随着CIP投加浓度的增加和运行时间的延长, 亚硝化单胞菌属和硝化螺旋菌属相对丰度呈先降低后升高再降低的趋势; 产碱菌属和硝化杆菌属持续减少, 最后消失.这4种菌属的减少或消失对MBR中氨氮的筛除造成了负面影响, 平均去除率由0 mg ?L-1的96.67%降低至15 mg ?L-1的66.52%.虽然硝化微生物对低浓度CIP可能具有一定的抵抗性和适应性, 但在高浓度CIP环境中硝化微生物却很难成为优势菌群, 导致(cause)氨氮的去除效率降低.
图 7
图 7 CIP投加浓度对硝化菌属的影响
2.2 CIP对抗性基因的影响
由图 3~5可知, 在CIP投加浓度为5 mg ?L-1的阶段, MBR中微生物在门、科和属水平上的优势群落结构分布均存在显著（striking）性差异.因此选取该阶段为研究对象, 考察CIP-ARGs平均相对丰度随时间的变化情况, 判断ARGs水平是否也存在类似现象, 结果如表 3所示.本研究选取的抗性基因包括DNA促旋酶A亚单位与B亚单位和拓扑异构酶Ⅳ的C亚单位与E亚单位.为了最大化减少分析检测过程中产生的误差, 采用CIP-ARGs的相对丰度进行分析.表 3显示, CIP的加入对MBR微生物产生了较为显著的影响, 其中内参基因16S rRNA的绝对丰度随CIP暴露时间的增加总体上呈降低趋势, 由3.45×107 copies ?g-1降至33 d的1.68×107 copies ?g-1; 4种CIP-ARGs在运行3 d时, 相对丰度均呈现不同程度的降低, 以gyrA基因降低幅度最大, 原因可能是携带这4种抗性基因的部分菌种对高浓度CIP敏感而被杀死.之后随着运行时间的增加, 这4种CIP-ARGs的相对丰度呈波动状态, 但相较于CIP投加初期均有所增加, 原因可能在于细菌长期暴露在CIP环境中, 其DNA旋转酶和拓扑异构酶Ⅳ基因的抗性决定区发生基因突变导致细菌对CIP出现抗性作用, 抗药性细菌可能会逐渐被选择成为优势菌群, 这与微生物群落结构分布结果相一致.
表 3 CIP-ARGs在MBR运行中的波动变化
3 结论
通过(tōng guò)考察不同CIP投加浓度下MBR污泥中微生物(Micro-Organism)群落结构变化发现, 随着CIP投加浓度和运行时间的增加, 工况4末, 门水平上Proteobacteria和Bacteroidetes仍然为MBR中的优势（解释:能压倒对方的有利形势)菌群, 科水平上Rhodocyclacea
E、Chitinophagaceae和Comamonadaceae被选择成为优势菌科, 属水平上Methyloversatili
S、Ferruginibacte
R、Zoogloea和Comamonas被选择成为优势菌属.Chao1、AC
E、Shannon指数降低和Simpson指数升高表明MBR中微生物丰富度和多样性均降低.高浓度CIP对硝化微生物Nitrosomona
S、Nitrospir
A、Alcaligenes和Nitrobacter的活性具有抑制作用, 导致氨氮去除效果下降.
CIP-ARGs分析表明高浓度CIP的长期胁迫使得MBR中的gyr
A、gyrB和parC基因增加.
本研究仅为CIP对MBR中微生物群落结构特征和抗性基因的影响提供数据支撑,
