城市污水处理厂微生物群落特征与群体感应调控

李松亚, 王林裴, 吴俊峰, 刘彪, 王乐, 赵琪彭, 李国琪, 黄善强, 费学宁. 城市污水处理厂微生物群落特征与群体感应调控[J]. 环境工程学报, 2022, 16(11): 3705-3715. doi: 10.12030/j.cjee.202207120
引用本文: 李松亚, 王林裴, 吴俊峰, 刘彪, 王乐, 赵琪彭, 李国琪, 黄善强, 费学宁. 城市污水处理厂微生物群落特征与群体感应调控[J]. 环境工程学报, 2022, 16(11): 3705-3715. doi: 10.12030/j.cjee.202207120
LI Songya, WANG Linpei, WU Junfeng, LIU Biao, WANG Le, ZHAO Qipeng, LI Guoqi, HUANG Shanqiang, FEI Xuening. Microbial community characteristics and their quorum sensing regulation in municipal wastewater treatment plants[J]. Chinese Journal of Environmental Engineering, 2022, 16(11): 3705-3715. doi: 10.12030/j.cjee.202207120
Citation: LI Songya, WANG Linpei, WU Junfeng, LIU Biao, WANG Le, ZHAO Qipeng, LI Guoqi, HUANG Shanqiang, FEI Xuening. Microbial community characteristics and their quorum sensing regulation in municipal wastewater treatment plants[J]. Chinese Journal of Environmental Engineering, 2022, 16(11): 3705-3715. doi: 10.12030/j.cjee.202207120

城市污水处理厂微生物群落特征与群体感应调控

    作者简介: 李松亚 (1990—) ,男,博士,讲师,20201013@hncj.edu.cn
    通讯作者: 李松亚(1990—),男,博士,讲师,20201013@hncj.edu.cn
  • 基金项目:
    河南省科技攻关项目(212102310277);国家自然科学基金资助项目(51778398);2021年河南城建学院大学生科技创新基金资助项目
  • 中图分类号: X703.1

Microbial community characteristics and their quorum sensing regulation in municipal wastewater treatment plants

    Corresponding author: LI Songya, 20201013@hncj.edu.cn
  • 摘要: 微生物群落在活性污泥中发挥有重要作用,为了解采用活性污泥法的城市污水处理厂微生物群落特征和群体感应调控作用,采集天津市某城市污水处理厂活性污泥样品,采用高通量测序对其微生物群落特征进行解析,并分析其群体感应调控作用。结果表明,该污水处理厂一年间活性污泥样品的微生物多样性和微生物群落结构存在一定差异性,主要是由于冬春季节发生了污泥膨胀。Ferruginibacterf_Saprospiraceae_unclassified、Terrimonasf_Blastocatellaceae_unclassified和metagenome等是该污水处理厂中主要优势菌属,冬春季节引发污泥膨胀的主要菌属是Candidatus Microthrix。群体感应调控作用分析表明C6-HSL能够抑制污泥膨胀,3OC12-HSL能够调控促进多糖的产生,f_Saprospiraceae_unclassified、f_Caldilineaceae_unclassified、SWB02、OLB12是主要的C6-HSL产生菌,f_Bacteroidetes_vadinHA17_unclassified是3OC12-HSL的群体淬灭菌。基于高通量测序的Tax4fun功能预测分析表明5月份污水处理厂的脱氮除磷性能均略逊于其他月份。
  • 加载中
  • 图 1  基于PCoA 的微生物群落差异分析

    Figure 1.  PCoA-based analysis of microbial community differences

    图 2  门水平微生物群落相对丰度

    Figure 2.  Relative abundance of microbial communities at phylum level

    图 3  主要微生物群落在属水平上的相对丰度

    Figure 3.  Relative abundances of major bacterial community at genus level

    图 4  信号分子浓度变化

    Figure 4.  Changes in signaling molecules concentration

    图 5  群体感应调控相关性分析

    Figure 5.  Correlation analysis for quorum sensing regulation

    图 6  AHLs与菌属之间相关性

    Figure 6.  Correlations of AHLs and genera

    图 7  代谢功能基因相对丰度的预测分析

    Figure 7.  Predictive analysis of relative abundance of metabolic function genes

    表 1  样品中微生物多样性指数

    Table 1.  Microbial diversity index in samples

    样品ACEChao1ShannonSimpson覆盖率/%
    S1613.16613.507.0150.98399.8
    S3588.12600.846.9200.98399.8
    S5576.22577.287.0090.97899.9
    S7646.34644.047.1670.98299.8
    S9643.03644.966.8260.97499.8
    S11647.57650.177.1180.98299.9
    样品ACEChao1ShannonSimpson覆盖率/%
    S1613.16613.507.0150.98399.8
    S3588.12600.846.9200.98399.8
    S5576.22577.287.0090.97899.9
    S7646.34644.047.1670.98299.8
    S9643.03644.966.8260.97499.8
    S11647.57650.177.1180.98299.9
    下载: 导出CSV

    表 2  功能菌属相对丰度

    Table 2.  Relative abundance of functional bacteria %

    功能菌分类功能菌属相对丰度参考文献
    S1S3S5S7S9S11
    GAOCandidatus Competibacter0.170.320.50.470.481.50[3]
    PAOTetrasphaera0.871.320.480.190.030.02[3, 23]
    Candidatus_Accumulibacter0.020.040.290.110.050.11[3]
    DPBf_Saprospiraceae_unclassified3.735.994.910.567.235.82[19]
    Paracoccus1.352.390.570.580.150.27[24]
    AOBNitrosomonas0.180.030.040.950.930.61[3]
    NOBNitrospira0.760.10.041.420.631.25[3]
    DNBFerruginibacter9.268.2213.369.9910.815.19[25]
    Terrimonas5.358.017.483.604.692.79[26]
    f_Blastocatellaceae_unclassified3.270.970.966.3910.478.74[18]
    Amphiplicatus1.363.135.461.610.710.48[27]
    Ottowia1.383.411.190.590.490.46[28]
    Defluviimonas1.402.220.871.070.600.91[29]
    Thauera0.480.370.391.841.952.03[3]
    Dokdonella0.240.791.20.761.030.33[11]
    Hyphomicrobium0.860.560.790.540.610.99[30]
    f_Rhodobacteraceae_unclassified0.872.070.680.330.140.19[30]
    Rhodobacter0.841.850.610.200.090.06[3]
    Novosphingobium0.360.70.510.220.160.18[30]
    Devosia0.620.720.280.170.10.12[30]
    Thiobacillus1.280.280.290.030.010.02[31]
    EPS分泌菌Terrimonas5.358.017.483.604.692.79[32]
    Thauera0.480.370.391.841.952.03[30]
    Devosia0.620.720.280.170.10.12[30]
    水解菌Ferruginibacter9.268.2213.369.9910.815.19[17]
    f_Saprospiraceae_unclassified3.735.994.910.567.235.82[33]
    Ottowia1.383.411.190.590.490.46[28]
    功能菌分类功能菌属相对丰度参考文献
    S1S3S5S7S9S11
    GAOCandidatus Competibacter0.170.320.50.470.481.50[3]
    PAOTetrasphaera0.871.320.480.190.030.02[3, 23]
    Candidatus_Accumulibacter0.020.040.290.110.050.11[3]
    DPBf_Saprospiraceae_unclassified3.735.994.910.567.235.82[19]
    Paracoccus1.352.390.570.580.150.27[24]
    AOBNitrosomonas0.180.030.040.950.930.61[3]
    NOBNitrospira0.760.10.041.420.631.25[3]
    DNBFerruginibacter9.268.2213.369.9910.815.19[25]
    Terrimonas5.358.017.483.604.692.79[26]
    f_Blastocatellaceae_unclassified3.270.970.966.3910.478.74[18]
    Amphiplicatus1.363.135.461.610.710.48[27]
    Ottowia1.383.411.190.590.490.46[28]
    Defluviimonas1.402.220.871.070.600.91[29]
    Thauera0.480.370.391.841.952.03[3]
    Dokdonella0.240.791.20.761.030.33[11]
    Hyphomicrobium0.860.560.790.540.610.99[30]
    f_Rhodobacteraceae_unclassified0.872.070.680.330.140.19[30]
    Rhodobacter0.841.850.610.200.090.06[3]
    Novosphingobium0.360.70.510.220.160.18[30]
    Devosia0.620.720.280.170.10.12[30]
    Thiobacillus1.280.280.290.030.010.02[31]
    EPS分泌菌Terrimonas5.358.017.483.604.692.79[32]
    Thauera0.480.370.391.841.952.03[30]
    Devosia0.620.720.280.170.10.12[30]
    水解菌Ferruginibacter9.268.2213.369.9910.815.19[17]
    f_Saprospiraceae_unclassified3.735.994.910.567.235.82[33]
    Ottowia1.383.411.190.590.490.46[28]
    下载: 导出CSV

    表 3  功能基因相对丰度

    Table 3.  Relative abundance of functional gene %

    代谢功能基因相对丰度
    S1S3S5S7S9S11
    氨基酸代谢13.9413.8812.8213.6312.9613.52
    碳水化合物代谢12.5812.6611.7712.3311.8712.26
    膜转运12.4212.0411.0712.0111.5811.79
    辅酶因子和维生素代谢7.377.357.197.337.277.34
    能量代谢7.357.307.327.377.437.45
    信号转导6.626.747.476.947.387.05
    核苷酸代谢5.445.455.645.545.625.61
    转译4.344.454.694.444.574.52
    外来生物的生物降解和代谢4.684.614.334.564.344.45
    复制和修复3.803.884.484.014.334.10
    脂质代谢3.693.693.723.673.633.61
    其他氨基酸代谢2.982.972.952.952.902.89
    萜类和多酮类化合物代谢2.572.572.402.502.382.44
    折叠、分类和降解2.242.262.342.282.332.31
    细胞运动2.032.002.642.212.582.29
    聚糖的生物合成和代谢2.022.112.552.152.412.23
    细胞生长和死亡1.661.621.921.741.931.79
    其他次生代谢物生物合成0.910.920.900.890.870.88
    运输和分解代谢0.350.360.360.350.340.34
    细胞群落—原核生物0.320.320.360.320.350.33
    转录0.190.190.260.210.240.21
    信号分子和相互作用0.0020.0010.020.0070.0160.009
    代谢功能基因相对丰度
    S1S3S5S7S9S11
    氨基酸代谢13.9413.8812.8213.6312.9613.52
    碳水化合物代谢12.5812.6611.7712.3311.8712.26
    膜转运12.4212.0411.0712.0111.5811.79
    辅酶因子和维生素代谢7.377.357.197.337.277.34
    能量代谢7.357.307.327.377.437.45
    信号转导6.626.747.476.947.387.05
    核苷酸代谢5.445.455.645.545.625.61
    转译4.344.454.694.444.574.52
    外来生物的生物降解和代谢4.684.614.334.564.344.45
    复制和修复3.803.884.484.014.334.10
    脂质代谢3.693.693.723.673.633.61
    其他氨基酸代谢2.982.972.952.952.902.89
    萜类和多酮类化合物代谢2.572.572.402.502.382.44
    折叠、分类和降解2.242.262.342.282.332.31
    细胞运动2.032.002.642.212.582.29
    聚糖的生物合成和代谢2.022.112.552.152.412.23
    细胞生长和死亡1.661.621.921.741.931.79
    其他次生代谢物生物合成0.910.920.900.890.870.88
    运输和分解代谢0.350.360.360.350.340.34
    细胞群落—原核生物0.320.320.360.320.350.33
    转录0.190.190.260.210.240.21
    信号分子和相互作用0.0020.0010.020.0070.0160.009
    下载: 导出CSV
  • [1] SHEIK A R, MULLER E E L, WILMES P. A hundred years of activated sludge: time for a rethink[J]. Frontiers in Microbiology, 2014, 5: 47.
    [2] WU L, NING D, ZHANG B, et al. Global diversity and biogeography of bacterial communities in wastewater treatment plants[J]. Nature Microbiology, 2019, 4(7): 1183-1195. doi: 10.1038/s41564-019-0426-5
    [3] 鞠峰, 张彤. 活性污泥微生物群落宏组学研究进展[J]. 微生物学通报, 2019, 46(8): 2038-2052.
    [4] 马切切, 袁林江, 牛泽栋, 等. 活性污泥微生物群落结构及与环境因素响应关系分析[J]. 环境科学, 2021, 42(8): 3886-3893.
    [5] 程珂珂, 曾艳华, 蔡中华, 等. 微生物的交流信号[J]. 生物化学与生物物理进展, 2022, 49(6): 960-974.
    [6] LI S, FEI X, CHI Y, et al. Impact of the acetate/oleic acid ratio on the performance, quorum sensing, and microbial community of sequencing batch reactor system[J]. Bioresource Technology, 2020, 296: 122279. doi: 10.1016/j.biortech.2019.122279
    [7] HU H, LIU Y, LUO F, et al. Stable and rapid partial nitrification achieved by boron stimulating autoinducer-2 mediated quorum sensing at room & low temperature[J]. Chemosphere, 2022, 304: 135327. doi: 10.1016/j.chemosphere.2022.135327
    [8] WANG N, GAO J, LIU Y, et al. Realizing the role of N-acyl-homoserine lactone-mediated quorum sensing in nitrification and denitrification: A review[J]. Chemosphere, 2021, 274: 129970. doi: 10.1016/j.chemosphere.2021.129970
    [9] LIU Q, FENG X, SHENG Z, et al. Enhanced wastewater treatment performance by understanding the interaction between algae and bacteria based on quorum sensing[J]. Bioresource Technology, 2022, 354: 127161. doi: 10.1016/j.biortech.2022.127161
    [10] MADDELA N R, SHENG B, YUAN S, et al. Roles of quorum sensing in biological wastewater treatment: A critical review[J]. Chemosphere, 2019, 221: 616-629. doi: 10.1016/j.chemosphere.2019.01.064
    [11] LI S, FEI X, CAO L, et al. Insights into the effects of carbon source on sequencing batch reactors: Performance, quorum sensing and microbial community[J]. Science of the Total Environment, 2019, 691: 799-809. doi: 10.1016/j.scitotenv.2019.07.191
    [12] AßHAUER K P, WEMHEUER B, DANIEL R, et al. Tax4Fun: predicting functional profiles from metagenomic 16S rRNA data[J]. Bioinformatics, 2015, 31(17): 2882-2884. doi: 10.1093/bioinformatics/btv287
    [13] CHU H, LIU X, MA J, et al. Two-stage anoxic-oxic (A/O) system for the treatment of coking wastewater: Full-scale performance and microbial community analysis[J]. Chemical Engineering Journal, 2021, 417: 129204. doi: 10.1016/j.cej.2021.129204
    [14] FENG Z, LI T, LIN Y, et al. Microbial communities and interactions in full-scale A2/O and MBR wastewater treatment plants[J]. Journal of Water Process Engineering, 2022, 46: 102660. doi: 10.1016/j.jwpe.2022.102660
    [15] WANG P, YU Z, QI R, et al. Detailed comparison of bacterial communities during seasonal sludge bulking in a municipal wastewater treatment plant[J]. Water Research, 2016, 105: 157-166. doi: 10.1016/j.watres.2016.08.050
    [16] 贺赟, 李魁晓, 王佳伟, 等. 不同季节城市污水处理厂微生物群落特性[J]. 环境科学, 2021, 42(3): 1488-1495. doi: 10.13227/j.hjkx.202007015
    [17] 黎镛, 袁辉洲, 柯水洲, 等. 微生物载体对MBBR工艺性能及微生物群落结构的影响[J]. 环境工程, 2021, 39(12): 100-106. doi: 10.13205/j.hjgc.202112015
    [18] ZHANG L, HAO S, WANG Y, et al. Rapid start-up strategy of partial denitrification and microbially driven mechanism of nitrite accumulation mediated by dissolved organic matter[J]. Bioresource Technology, 2021, 340: 125663. doi: 10.1016/j.biortech.2021.125663
    [19] 顾晓丹, 黄勇, 丁永伟, 等. 改良型UNITANK工艺冬季运行特性及微生物群落分析[J]. 环境工程学报, 2021, 15(7): 2480-2487. doi: 10.12030/j.cjee.202102134
    [20] 尚越飞, 王申, 宗倪, 等. 污水生物处理工艺低温下微生物种群结构[J]. 环境科学, 2020, 41(10): 4636-4643. doi: 10.13227/j.hjkx.202004116
    [21] FAN N, QI R, HUANG B, et al. Factors influencing Candidatus Microthrix parvicella growth and specific filamentous bulking control: A review[J]. Chemosphere, 2020, 244: 125371. doi: 10.1016/j.chemosphere.2019.125371
    [22] WANG J, QI R, LIU M, et al. The potential role of ‘Candidatus Microthrix parvicella’ in phosphorus removal during sludge bulking in two full-scale enhanced biological phosphorus removal plants[J]. Water Science and Technology, 2014, 70(2): 367-375. doi: 10.2166/wst.2014.216
    [23] LIU H, ZENG W, MENG Q, et al. Phosphorus removal performance, intracellular metabolites and clade-level community structure of Tetrasphaera-dominated polyphosphate accumulating organisms at different temperatures[J]. Science of the Total Environment, 2022, 842: 156913. doi: 10.1016/j.scitotenv.2022.156913
    [24] 叶丽红, 李冬, 张杰, 等. 亚硝化-反硝化除磷技术研究进展[J]. 北京工业大学学报, 2016, 42(4): 585-593. doi: 10.11936/bjutxb2015070009
    [25] 赖永恒, 刘敏. A/O-MBBR工艺处理生活污水性能及菌群结构[J]. 中国环境科学, 2022, 42(5): 2120-2128. doi: 10.3969/j.issn.1000-6923.2022.05.016
    [26] CHEN X, XU P, YANG C, et al. Study of enhanced nitrogen removal efficiency and microbial characteristics of an improved two-stage A/O process[J]. Environmental Technology, 2021, 42(27): 4306-4316. doi: 10.1080/09593330.2020.1754924
    [27] 张晶晶, 杨翠春, 丁鹏霖, 等. MBBR工艺用于唐山某污水厂提标改造效能分析[J]. 中国给水排水, 2020, 36(15): 78-85.
    [28] LI J, DU Q, PENG H, et al. Optimization of biochemical oxygen demand to total nitrogen ratio for treating landfill leachate in a single-stage partial nitrification-denitrification system[J]. Journal of Cleaner Production, 2020, 266: 121809. doi: 10.1016/j.jclepro.2020.121809
    [29] WANG L, ZHAN H, WU G, et al. Effect of operational strategies on the rapid start-up of nitrogen removal aerobic granular system with dewatered sludge as inoculant[J]. Bioresource Technology, 2020, 315: 123816. doi: 10.1016/j.biortech.2020.123816
    [30] HE Q, CHEN L, ZHANG S, et al. Hydrodynamic shear force shaped the microbial community and function in the aerobic granular sequencing batch reactors for low carbon to nitrogen (C/N) municipal wastewater treatment[J]. Bioresource Technology, 2019, 271: 48-58. doi: 10.1016/j.biortech.2018.09.102
    [31] ZHANG R, CHEN C, XU X, et al. The interaction between Pseudomonas C27 and Thiobacillus denitrificans in the integrated autotrophic and heterotrophic denitrification process[J]. Science of the Total Environment, 2022, 811: 152360. doi: 10.1016/j.scitotenv.2021.152360
    [32] ZHAO Y, JIANG B, TANG X, et al. Metagenomic insights into functional traits variation and coupling effects on the anammox community during reactor start-up[J]. Science of the Total Environment, 2019, 687: 50-60. doi: 10.1016/j.scitotenv.2019.05.491
    [33] 郑向阳, 罗晓, 袁立霞, 等. AO工艺处理淀粉污水效能及微生物群落解析[J]. 环境工程学报, 2018, 12(3): 804-814. doi: 10.12030/j.cjee.201708047
    [34] 王萍, 余志晟. 污水处理厂污泥膨胀和污泥发泡的比较分析[J]. 微生物学通报, 2019, 46(8): 1971-1981.
    [35] HUANG J, SHI Y, ZENG G, et al. Acyl-homoserine lactone-based quorum sensing and quorum quenching hold promise to determine the performance of biological wastewater treatments: An overview[J]. Chemosphere, 2016, 157: 137-151. doi: 10.1016/j.chemosphere.2016.05.032
    [36] 罗景阳, 张乐, 张腾, 等. 不同工艺对污水处理效能及细菌代谢功能的影响[J]. 中国给水排水, 2021, 37(23): 61-67. doi: 10.19853/j.zgjsps.1000-4602.2021.23.011
    [37] FANG D, ZHAO G, XU X, et al. Microbial community structures and functions of wastewater treatment systems in plateau and cold regions[J]. Bioresource Technology, 2018, 249: 684-693. doi: 10.1016/j.biortech.2017.10.063
  • 加载中
图( 7) 表( 3)
计量
  • 文章访问数:  7711
  • HTML全文浏览数:  7711
  • PDF下载数:  189
  • 施引文献:  0
出版历程
  • 收稿日期:  2022-07-24
  • 录用日期:  2022-10-25
  • 刊出日期:  2022-11-30
李松亚, 王林裴, 吴俊峰, 刘彪, 王乐, 赵琪彭, 李国琪, 黄善强, 费学宁. 城市污水处理厂微生物群落特征与群体感应调控[J]. 环境工程学报, 2022, 16(11): 3705-3715. doi: 10.12030/j.cjee.202207120
引用本文: 李松亚, 王林裴, 吴俊峰, 刘彪, 王乐, 赵琪彭, 李国琪, 黄善强, 费学宁. 城市污水处理厂微生物群落特征与群体感应调控[J]. 环境工程学报, 2022, 16(11): 3705-3715. doi: 10.12030/j.cjee.202207120
LI Songya, WANG Linpei, WU Junfeng, LIU Biao, WANG Le, ZHAO Qipeng, LI Guoqi, HUANG Shanqiang, FEI Xuening. Microbial community characteristics and their quorum sensing regulation in municipal wastewater treatment plants[J]. Chinese Journal of Environmental Engineering, 2022, 16(11): 3705-3715. doi: 10.12030/j.cjee.202207120
Citation: LI Songya, WANG Linpei, WU Junfeng, LIU Biao, WANG Le, ZHAO Qipeng, LI Guoqi, HUANG Shanqiang, FEI Xuening. Microbial community characteristics and their quorum sensing regulation in municipal wastewater treatment plants[J]. Chinese Journal of Environmental Engineering, 2022, 16(11): 3705-3715. doi: 10.12030/j.cjee.202207120

城市污水处理厂微生物群落特征与群体感应调控

    通讯作者: 李松亚(1990—),男,博士,讲师,20201013@hncj.edu.cn
    作者简介: 李松亚 (1990—) ,男,博士,讲师,20201013@hncj.edu.cn
  • 1. 河南城建学院市政与环境工程学院,平顶山 467036
  • 2. 河南省水体污染防治与修复重点实验室,平顶山 467036
  • 3. 天津城建大学理学院,天津 300384
基金项目:
河南省科技攻关项目(212102310277);国家自然科学基金资助项目(51778398);2021年河南城建学院大学生科技创新基金资助项目

摘要: 微生物群落在活性污泥中发挥有重要作用,为了解采用活性污泥法的城市污水处理厂微生物群落特征和群体感应调控作用,采集天津市某城市污水处理厂活性污泥样品,采用高通量测序对其微生物群落特征进行解析,并分析其群体感应调控作用。结果表明,该污水处理厂一年间活性污泥样品的微生物多样性和微生物群落结构存在一定差异性,主要是由于冬春季节发生了污泥膨胀。Ferruginibacterf_Saprospiraceae_unclassified、Terrimonasf_Blastocatellaceae_unclassified和metagenome等是该污水处理厂中主要优势菌属,冬春季节引发污泥膨胀的主要菌属是Candidatus Microthrix。群体感应调控作用分析表明C6-HSL能够抑制污泥膨胀,3OC12-HSL能够调控促进多糖的产生,f_Saprospiraceae_unclassified、f_Caldilineaceae_unclassified、SWB02、OLB12是主要的C6-HSL产生菌,f_Bacteroidetes_vadinHA17_unclassified是3OC12-HSL的群体淬灭菌。基于高通量测序的Tax4fun功能预测分析表明5月份污水处理厂的脱氮除磷性能均略逊于其他月份。

English Abstract

  • 近年来,我国污水排放总量不断增加,大量的污染物排放到水环境中将给水环境造成严重影响,城市污水处理是关系国计民生的重要问题。城市污水处理厂作为城市重要的基础设施之一,在水环境污染防治和生态文明建设方面发挥重要作用。活性污泥法因其高效、经济、环保等优点,成为我国乃至世界上大部分城市污水处理厂普遍采用的工艺[1]。活性污泥作为活性污泥工艺的主体,含有高度复杂的微生物群落,构成一种独特的微生物生态系统,并依靠这些微生物降解有机物、营养物质等,从而实现活性污泥系统稳定运行和良好处理性能[2-3]。马切切等[4]对活性污泥微生物群落结构及与环境因素响应关系研究发现,温度与NitrospiraAquihabitansTerrimonasDechloromona等多种关键的脱氮微生物呈负相关,当温度低于15 ℃时,污水处理厂对TN和NH4+-N的去除率仅有49.67%和63.19%;ZoogloeaArcobacterAcidovoraxAcinetobacter等优势功能性菌属相对丰度越高,污染物的去除率也越高。因此,污水处理厂的稳定运行在很大程度上依赖于系统微生物群落的特征,特别是一些功能菌的多样性、丰度。深入分析、研究污水处理厂系统中微生物群落组成、功能菌的生理特性,从而选择性调控微生物群落,对于提高污水处理厂活性污泥处理能力具有至关重要的意义。

    群体感应(quorum sensing,QS)是微生物间通过自发产生、释放并接收特定化学信号而建立起来的一种通讯机制,所依赖的化学信号分子就是他们的“通讯语言”[5]。近年来,群体感应在有机物降解、生物脱氮等废水生物处理过程中的调控作用受到越来越多的关注[6-8]。LI等[6]在不同乙酸/油酸比条件下对序批式反应器(sequencing batch reactor,SBR)进行了研究,发现N酰化高丝氨酸内脂(N-acyl homoserine lactones,AHLs)介导的群体感应对化学需氧量(COD)和总氮(TN)去除、污泥沉降性等均具有显著的调控作用。WANG等[8]总结了AHLs介导的群体感应在废水生物反应器氮代谢效率、颗粒聚集和生物膜的形成等方面的调控作用。大量研究[9-10]表明,群体感应在废水生物处理中发挥重要作用。因此,本研究对天津市某城市污水处理厂进行了为期1 a的监测研究,重点分析了其微生物群落结构特征和群体感应信号分子变化以探索污水处理厂活性污泥系统的群体感应调控作用,以期为提升污水处理厂处理能力及稳定性提供参考。

    • 对天津市某城市污水处理厂的曝气池活性污泥进行样本采集,该城市污水处理厂主体采用A/O处理工艺,并辅以化学法去除尾水中余磷。采样时间为2019年1月、3月、5月、7月、9月和11月中旬,分别命名为S1、S3、S5、S7、S9和S11。样品采集后于−80 ℃保存,待所有样品收集完毕后进行测序。每月的5、15、25日取活性污泥混合液用于群体感应信号分子的提取和检测。

    • AHLs类信号分子采用Bond Elut ENV固相萃取柱(500 mg,6 mL,Agilent)进行固相萃取,具体操作步骤如下:用0.45 μm滤膜过滤水样,分别用5 mL甲醇(Fisher Scientific,LC-MS级)和5 mL纯水活化Bond Elut ENV固相萃取柱,然后以5~10 mL·min−1的速度将过滤后的水样载入固相萃取柱,之后用5mL 5%甲醇水溶液淋洗,待固相萃取柱抽干后用10 mL甲醇(Fisher Scientific,LC-MS级)洗脱萃取的AHLs信号分子,并收集洗脱液。将洗脱液用旋转蒸发仪脱溶后,使用初始流动相定容至2 mL,采用高效液相色谱(Agilent 1290 Infinity II)-质谱(Agilent 6470)联用仪参考LI等[11]采用的方法对AHLs进行定量分析。

    • 待所有活性污泥样本收集完毕后送往苏州金唯智生物科技有限公司进行DNA提取。使用Magen Hipure Soil DNA Kit从污泥样品中提取DNA,并使用Qubit® dsDNA HS Assay Kit 检测DNA样品浓度。采用金唯智设计的包含“CCTACGGRRBGCASCAGKVRVGAAT”序列的上游引物和包含“GGACTACNVGGGTWTCTAATCC”序列的下游引物对V3和V4这2个高度可变区进行扩增。对扩增纯化后的DNA通过Illumina MiSeq测序平台进行2×250bp双端测序。

    • 对序列进行质量过滤、去除嵌合体后采用VSEARCH(1.9.6)进行OTU聚类(97%相似度),利用RDP classifier贝叶斯算法对OTU进行分类学注释,参考数据库为Silva 132。根据OTU结果进行α多样性指数分析,通过主坐标分析(PCoA)比较各样本间微生物群落差异,利用Tax4Fun进行功能预测分析[12]

    • 本研究中城市污水处理厂不同活性污泥样本微生物α多样性指数如表1所示。由表1可知,所有活性污泥样本的覆盖率均不低于99.8%。这表明样本中的绝大部分微生物均被检测到,测序结果可以反映样本的真实情况。S5样本的ACE指数和Chao指数是所有样本中最低的,表明5月份微生物群落的丰富度较其他时期较低。这是由于此时处于污泥膨胀恢复期,样本内微生物的种类较少。Shannon指数和Simpson指数均反映样本的多样性,Shannon指数可以反映物种均匀度,其值越大表示物种均匀性更高,Simpson指数越大,其多样性越大。而9月份污泥样本的Shannon指数和Simpson指数比其他时期内污泥样品低,表明9月份样品中微生物群落的均匀度和多样性最低。这可能是由于9月份系统运行温度及运行条件适宜某些功能微生物生长,其相对丰度较高,从而使整个系统微生物群落的均匀度降低(Shannon指数低),同时部分微生物无法适应该运行条件而被演替使整个系统微生物多样性降低(Simpson指数降低)。

      为分析不同样本间群落组成的差异性,基于Bray-Curtis距离算法,采用PCoA分析衡量样本间物种多样性差异,结果如图1所示。由图1可知,PCO1与PCO2分别解释方差为61.45%和17.66%,累计解释能力达79.11%,PCO1将样本S1、S3、S5的OTU与S7、S9和S11明显区分开,样本S1、S3、S5的污泥处于污泥膨胀状态和污泥膨胀恢复状态,样本相似性高;而S7、S9和S11均处于污泥正常状态,故具有高度相似性。同时,注意到S1、S3和S5之间明显分离,而S7、S9和S11之间虽然存在分离,但距离较小,表明在1—5月份污水处理厂由于发生污泥膨胀,微生物群落发生较大的变化,而7—11月份污水处理厂微生物群落差异变化不大,系统比较稳定。

    • 1)门水平上微生物群落结构分析。在门水平上,所有样本均检测到18~26个菌门,其中S3样本仅18个菌门,而S9和S11样本检测到26个菌门,相对丰度在前10的菌门见图2。由图2可见,S1、S3、S5样本中相对丰度较高的菌门为BacteroidetesProteobacteriaPatescibacteria;而S7、S9和S11中主要优势菌门为BacteroidetesProteobacteriaAcidobacteria。在活性污泥系统中,ProteobacteriaBacteroidetes作为优势菌门已被多次报道[13-14]。在废水生物处理系统中,绝大部分具有去除有机物和脱氮除磷等功能的微生物属于ProteobacteriaBacteroidetes菌门[14],其相对较高的丰度可保证该污水处理厂污染物的稳定去除。

      同时值得注意的是,不同样本的菌门相对丰度存在较大的差异,如ActinobacteriaChloroflexi,这可能主要是由于污水处理厂进水水质差异以及不同时间段环境温度差异所致。有研究表明,ActinobacteriaChloroflexi是活性污泥中的主要丝状菌门,其中属于Actinobacteria的丝状菌如Candidatus Microthrix常常引发污泥膨胀,而隶属Chloroflexi的丝状菌通常充当菌胶团骨架[15]。因此,由图2可以看到,在污水处理厂发生污泥膨胀的S1和S3样本中Actinobacteria相对丰度高于其他样本,而在其他样本中Chloroflexi相对丰度比在S1和S3样本中的相对丰度高。

      2)属水平微生物群落结构。属水平上微生物群落结构分析更有助于理解系统微生物群落的组成和演替,6个污水处理厂活性污泥样本中相对丰度前30的细菌属见图3。由图3可以看出,该污水处理厂活性污泥中主要菌属有Ferruginibacterf_Saprospiraceae_unclassified、Terrimonasf_Blastocatellaceae_unclassified和metagenome,相对丰度分别为5.19%~13.36%、3.73%~10.56%、2.79%~8.01%、0.96%~10.47%、1.06%~10.22%。贺赟等[16]对不同季节城市污水处理厂微生物群落研究也发现FerruginibacterBlastocatellaceae是污水处理厂的主要优势菌属。Ferruginibacter具有将复杂有机物分解为小分子有机物的能力[17]Blastocatellaceae能够降解复杂有机物同时可以完全还原NO- 3-N[18]Saprospiraceae是活性污泥系统中的核心菌属,其不仅能够降解有机物和蛋白质,还具有反硝化除磷能力[19]Terrimonas是一种典型的反硝化菌[20]。综上所述,这些优势菌属均为能够去除污染物的功能菌,其较高的相对丰度保证了污水处理厂的稳定运行。

      同时可以注意到,经常报道能够引发活性污泥膨胀的一种丝状菌属Candidatus Microthrix [21]有被检测到。由图3可以看到,污水处理厂S1和S3样本中Candidatus Microthrix相对丰度较高,分别为3.03%和4.24%,在S5和S7样本中其相对丰度降低至1.49%和0.28%,并在S9和S11中相对丰度进一步降低。在1—3月份,该污水处理厂发生污泥膨胀,污泥SVI均在150 mL·g−1以上,自3月下旬开始其逐渐降低,污泥沉降性能逐渐恢复。由此可以推断,该污水处理厂污泥沉降性能的恶化主要是由于Candidatus Microthrix过度生长所导致,这与春季低温条件利于Candidatus Microthrix的生长有关[21]。故Candidatus Microthrix是1—3月份引发污泥膨胀的主要菌属。这与WANG等[22]对北方城市污水处理厂研究发现是冬春季节引发污泥膨胀的主要菌属是Candidatus Microthrix的研究结果是一致的。

      3)功能菌群分析。功能菌群对污水处理厂的稳定运行具有重要作用,活性污泥系统中目前已报道的主要功能菌群有聚糖菌(GAO)、聚磷菌(PAO)、氨氧化菌(AOB)、亚硝酸盐氧化菌(NOB)、反硝化菌(DNB)、反硝化聚磷菌(DPB)、胞外聚合物(EPS)分泌菌等。根据对所研究污水处理厂样本在属水平上检测到的功能菌群进行分类,结果如表2所示。

      典型聚糖菌Candidatus Competibacter在该污水处理厂样本内检测到,其相对丰度(0.17%~1.5%)高于典型聚磷菌Candidatus Accumulibacter的相对丰度(0.02%~0.29%),较低丰度的Candidatus Accumulibacter则会影响系统的除磷效果。Candidatus Accumulibacter与Candidatus Competibacter间的竞争关系已被多次报道,系统中Candidatus Competibacter相对丰度较高,则会抑制Candidatus Accumulibacter,使系统除磷效果较差。此外,另一种具有除磷功能的聚磷菌Tetrasphaera[23]在该污水处理厂样品中检出,且其相对丰度(0.02%~1.32%)也显著高于Candidatus Accumulibacter。Tetrasphaera也是一种引起污泥膨胀的丝状菌属,该属的丝状菌能够利用疏水性的长链脂肪酸等,使其在低温条件下占据优势生态位[34]。可以注意到,在污水处理厂发生污泥膨胀期间(1月和3月)Tetrasphaera的相对丰度较高,而非膨胀期间(9月和11月)其相对丰度较低,故该污水处理厂中检测到的Tetrasphaera可能主要是引发污泥膨胀的丝状菌。

      此外,2种具有反硝化除磷功能的菌属f_Saprospiraceae_unclassified和Paracoccus也在该污水处理厂样品中检测到。对全球污水处理厂的调查研究表明,Saprospiraceae是活性污泥中的核心菌属[2],其具有代谢葡萄糖、半乳糖、特定蛋白质等能力[33]。同时,Saprospiraceae能够将亚硝酸盐氮作为电子受体,是反硝化除磷系统中的优势菌[19]Paracoccus在缺氧条件下能够以亚硝酸盐氮为电子受体氧化厌氧条件下合成的聚β-羟基丁酸盐(PHB)以产生能量用于吸磷[24]。2种反硝化除磷功能菌的存在表明系统中存在反硝化除磷现象,且f_Saprospiraceae_unclassified和Paracoccus的相对丰度(3.73%~10.56%和0.15%~2.39%)明显高于聚磷菌Candidatus Accumulibacter的相对丰度,故该污水处理厂可能以反硝化除磷为主。

      NitrosomonasNitrospira分别是活性污泥中经常出现的AOB和NOB,在该污水处理厂中其相对丰度分别为0.03%~0.95%和0.04%~1.42%。在3月份和5月份的污泥样品中,NitrosomonasNitrospira的相对丰度均为最低,表明NitrosomonasNitrospira受污泥膨胀影响较大,污泥膨胀发生期间,丝状菌的过度生长抑制了NitrosomonasNitrospira的生长。

      此外,该污水处理厂污泥样品中还检测到了多种DNB菌属,丰度相对较高的是Ferruginibacter (5.19%~13.36%)、Terrimonas(2.79%~8.01%)、f_Blastocatellaceae_unclassified (0.96%~10.47%)、Amphiplicatus (0.48%~5.46%)、Ottowia (0.46%~3.41%)等。这些相对丰度较高的反硝化菌是系统实现脱氮的保证。Ferruginibacter菌属能够水解尿素、明胶等有机物,具备将大分子有机物分解为小分子的能力[17],同时也具有反硝化能力[25]。1—9月的样品中Ferruginibacter均为优势反硝化菌,11月样品中优势反硝化菌却为f_Blastocatellaceae_unclassified,故Ferruginibacter对该污水处理厂有机物的降解和脱氮具有重要作用。f_Blastocatellaceae_unclassified能够代谢复杂有机聚合物,同时它也可以利用聚羟基脂肪酸酯(PHA)完全还原硝酸盐氮[18]。故1—9月Ferruginibacter与11月f_Blastocatellaceae_unclassified优势反硝化菌转变的主要原因可能是进水水质差异(如有机物种类、含量等)所致。可以注意到,在整个研究周期内活性污泥样品中f_Blastocatellaceae_unclassified相对丰度变化比较大,污泥膨胀发生期间(3月和5月)f_Blastocatellaceae_unclassified的相对丰度较其他月份显著降低,表明其受污泥膨胀影响较大。据报道,Terrimonas是一种好氧反硝化菌属,能够还原硝酸盐氮,是改良两级A/O和传统两级A/O工艺中的优势菌属[26]Amphiplicatus是严格好氧型细菌,具备反硝化脱氮性能,却不能水解纤维素、果胶、酪蛋白等基质[27]Ottowia不仅报道与反硝化作用有关,还是一种重要的具有水解功能细菌,所以能够同时发挥脱氮和有机物降解的作用[28]。此外,经常报道的优势反硝化菌如Thauera等,在该污水厂中相对丰度并不高,且可以注意到在1—5月其相对丰度较低,表明其受低温影响较大。这些反硝化菌的存在保证该污水处理厂脱氮性能的稳定。

      能够分泌EPS的絮凝菌也是活性污泥微生物的重要组成部分。有报道称,Terrimonas通过分泌EPS在微生物聚集体中起促进絮凝和聚集的重要作用[32]。HE等[30]总结了不同操作条件下能够产生EPS的关键功能菌,其中ThaueraDevosia能够分泌EPS并进行反硝化。这些EPS产生菌有助于维持活性污泥的聚集和稳定以及对部分营养物质的去除。

    • 1)群体感应信号分子。在该污水处理厂样品中检测到N-己酰高丝氨酸环内酯(C6-HSL)、3-羰基辛酰基高丝氨酸内酯(3OC8-HSL)和3-羰基十二酰基高丝氨酸内酯(3OC12-HSL)共3种AHLs,其每月平均浓度变化如图4所示。由图4可知,3OC8-HSL和3OC12-HSL的质量浓度相对比较稳定,分别维持在2.7 ng·L−1和0.6 ng·L−1(略有浮动)。而C6-HSL变化较大,1—7月其质量浓度逐渐由0.09 ng·L−1升高至7.8 ng·L−1,此后又降低至3.7 ng·L−1,而在11月份样品中未检出。可以注意到,较低气温下C6-HSL浓度也较低,而较高温度下其浓度也比较高,表明C6-HSL的产生可能与温度呈现一定的相关性。另外,污泥膨胀期间,C6-HSL浓度较低,而非膨胀期间,其浓度相对比较高,故C6-HSL可能与污泥膨胀或丝状菌生长存在调控关系。

      2)群体感应与污泥沉降性。为研究群体感应调控作用,进行Person相关性分析结果如图5所示。可见污泥容积指数(SVI)与C6-HSL呈显著负相关(P<0.05),与3OC12-HSL呈显著正相关(P<0.05),则C6-HSL、3OC12-HSL均与污泥沉降性有关。其中,C6-HSL能够调控抑制系统污泥膨胀的发生,使系统污泥沉降性能良好,而3OC12-HSL能够调控引发系统污泥膨胀的发生,使系统污泥沉降性能恶化。C6-HSL与污泥膨胀之间的调控作用与之前LI等[6]的研究结果是一致的,故鉴于C6-HSL对污泥膨胀的调控作用可以向发生污泥膨胀的系统中投加C6-HSL或者促进C6-HSL信号分子的产生以控制污泥膨胀。此外,3OC12-HSL与多糖(PS)含量呈显著正相关(P<0.05),表明3OC12-HSL主要调控多糖的产生。群体感应对EPS产生的调控作用在不同反应器和城市污水处理厂中均有报道,但在本研究中并没有观察到,这可能主要是工艺、运行条件等的差异[35]

      3)群体感应与微生物群落调控作用。为分析群体感应与微生物群落间的调控作用,将相对丰度前50的菌属与AHLs进行Person相关分析,与AHLs信号分子具有显著相关的细菌属只有f_Saprospiraceae_unclassified、f_Bacteroidetes_vadinHA17_unclassified、SWB02、OLB12、f_Caldilineaceae_unclassified(图6)。f_Saprospiraceae_unclassified、f_Caldilineaceae_unclassified、SWB02与C6-HSL呈显著正相关,其相关系数分别为0.913、0.864、0.904,OLB12与C6-HSL呈极显著正相关,其相关系数为0.982,表明f_Saprospiraceae_unclassified、f_Caldilineaceae_unclassified、SWB02、OLB12是主要的C6-HSL产生菌。而f_Bacteroidetes_vadinHA17_unclassified与3OC12-HSL呈显著负相关,相关系数为−0.819,表明f_Bacteroidetes_vadinHA17_unclassified是3OC12-HSL的群体淬灭菌。

    • 为进一步分析该污水处理厂微生物群落的关键功能基因表达情况,采用Tax4Fun基于该污水处理厂污泥样品高通量测序进行KEGG功能基因预测,分析结果见表3。氨基酸代谢(amino acid metabolism)、碳水化合物代谢(carbohydrate metabolism)、膜转运(membrane transport)、能量代谢(energy metabolism)和辅酶因子和维生素代谢(metabolism of cofactors and vitamins)为该污水处理厂中相对丰度较高的代谢功能。氨基酸代谢(amino acid metabolism)、碳水化合物代谢(carbohydrate metabolism)、能量代谢(energy metabolism)和辅酶因子和维生素代谢(metabolism of cofactors and vitamins)均属于新陈代谢(metabolism)功能模块,而膜转运(membrane transport)属于环境信息处理(environmental information processing)功能模块。碳水化合物代谢(carbohydrate metabolism)途径可以将大分子碳水化合物分解、稳定为小分子物质,氨基酸代谢(amino acid metabolism)过程可为反硝化脱氮提供电子,能量代谢(energy metabolism)途径在ATP代谢、电子传输等能量运输中起着关键作用[36],膜转运(membrane transport)在物质运输、群体感应调控方面发挥重要作用[30],这些代谢功能保证了该污水处理厂的正常运行。可以发现一些代谢功能的相对丰度在不同样品中差异较大,可能是由进水水质、运行条件等差异所导致,从而影响污水处理厂有机物的降解、脱氮除磷等性能。

      此外,通过深入分析脱氮除磷过程中相关基因的相对丰度表达情况,发现4种参与反硝化过程的nor(一氧化氮还原酶)、nos(一氧化二氮还原酶)、nir(亚硝酸盐还原酶)和nar(硝酸盐还原酶)[37],其中在硝酸盐氮还原为亚硝酸盐氮过程中起重要作用的nar(硝酸盐还原酶)相对丰度(0.104%~0.135%)明显高于其他还原酶((图7(a))。硝酸盐氮的还原是反硝化过程的第一步,是影响脱氮性能的关键过程,因此nar(硝酸盐还原酶)是影响脱氮性能的重要基因,对脱氮性能具有重要决定作用。S5样品的nar(硝酸盐还原酶)相对丰度是所有样品中最低的,同时该样品中nos(一氧化二氮还原酶)相对丰度也是最低的,表明5月份污水处理厂的脱氮性能较其他月份较差。

      发现4种参与除磷相关的功能基因,其中ppk(聚磷酸盐激酶)的相对丰度(0.125%~0.133%)明显高于phdp(聚-β-羟基丁酸解聚酶,相对丰度0.041%~0.051%)、hdh(羟丁酸脱氢酶,相对丰度0.036%~0.047%)和hh(羟基丁酸二聚体水解酶,相对丰度0.000 221%~0.000 297%)。ppk(聚磷酸盐激酶)与Poly-P形成相关,而其他3种均与聚-β-羟基丁酸降解有关[37],虽然具有较高的Poly-P形成能力,但受限于较低的聚-β-羟基丁酸降解能力,导致系统的除磷能力较低,故仍需借助化学除磷实现磷的去除。S5样品中与聚-β-羟基丁酸降解有关的3种功能基因相对丰度在所有样品中最低,表明5月份污水处理厂的磷去除能力略逊于其他月份。

    • 1)在1—3月污水处理厂发生污泥膨胀为微生物多样性和微生物群落结构存在差异的主要原因,Candidatus Microthrix过量生长是引发污泥膨胀的主要原因。

      2) BacteroidetesProteobacteriaAcidobacteria为该污水处理厂的主要优势菌门,Ferruginibacterf_Saprospiraceae_unclassified、Terrimonasf_Blastocatellaceae_unclassified和metagenome为该污水处理厂中主要优势菌属。Candidatus Accumulibacter、NitrosomonasNitrospiraFerruginibacterTerrimonas等相关功能菌的存在使系统维持相对稳定运行。

      3)污水处理厂中检测到C6-HSL、3OC8-HSL和3OC12-HSL3种AHLs,C6-HSL主要调控抑制系统污泥膨胀的发生,3OC12-HSL主要调控多糖的产生。群体感应与微生物群落间的调控作用分析表明f_Saprospiraceae_unclassified、f_Caldilineaceae_unclassified、SWB02、OLB12是主要的C6-HSL产生菌,f_Bacteroidetes_vadinHA17_unclassified是3OC12-HSL的群体淬灭菌。

      4)功能预测分析结果表明,5月份样品中nar(硝酸盐还原酶)、nos(一氧化二氮还原酶)、phdp(聚-β-羟基丁酸解聚酶)、hdh(羟丁酸脱氢酶)和hh(羟基丁酸二聚体水解酶)相对丰度均较其他样品低,故其脱氮除磷性能均略逊于其他月份。

    参考文献 (37)

返回顶部

目录

/

返回文章
返回