Processing math: 100%

利用16S rRNA高通量测序技术考察温度对生物脱氮硝化过程中亚硝酸盐氧化菌代谢功能的影响

李奕燃, 董昆, 于雪, 李晓强, 王光杰, 孙洪伟. 利用16S rRNA高通量测序技术考察温度对生物脱氮硝化过程中亚硝酸盐氧化菌代谢功能的影响[J]. 环境工程学报, 2022, 16(3): 980-988. doi: 10.12030/j.cjee.202104086
引用本文: 李奕燃, 董昆, 于雪, 李晓强, 王光杰, 孙洪伟. 利用16S rRNA高通量测序技术考察温度对生物脱氮硝化过程中亚硝酸盐氧化菌代谢功能的影响[J]. 环境工程学报, 2022, 16(3): 980-988. doi: 10.12030/j.cjee.202104086
LI Yiran, DONG Kun, YU Xue, LI Xiaoqiang, WANG Guanjie, SUN Hongwei. Investigation on the effect of temperature on the metabolism and function of nitrite-oxidizing bacteria in the process of biological nitrification by 16S rRNA high-throughput sequencing technology[J]. Chinese Journal of Environmental Engineering, 2022, 16(3): 980-988. doi: 10.12030/j.cjee.202104086
Citation: LI Yiran, DONG Kun, YU Xue, LI Xiaoqiang, WANG Guanjie, SUN Hongwei. Investigation on the effect of temperature on the metabolism and function of nitrite-oxidizing bacteria in the process of biological nitrification by 16S rRNA high-throughput sequencing technology[J]. Chinese Journal of Environmental Engineering, 2022, 16(3): 980-988. doi: 10.12030/j.cjee.202104086

利用16S rRNA高通量测序技术考察温度对生物脱氮硝化过程中亚硝酸盐氧化菌代谢功能的影响

    作者简介: 李奕燃 (1995—),男,硕士研究生,1587381875@qq.com
    通讯作者: 孙洪伟(1978—),男,博士,教授,sunhw@ytu.edu.cn
  • 基金项目:
    省级大学生创新创业训练计划项目(2017107320315)
  • 中图分类号: X703.5

Investigation on the effect of temperature on the metabolism and function of nitrite-oxidizing bacteria in the process of biological nitrification by 16S rRNA high-throughput sequencing technology

    Corresponding author: SUN Hongwei, sunhw@ytu.edu.cn
  • 摘要: 为探究温度对生物脱氮硝化过程微生物代谢及功能基因的影响,采用16S rRNA高通量测序技术,以富集的亚硝酸盐氧化菌(NOB)为研究对象,考察了活性污泥系统硝化过程中菌属与比亚硝酸盐氧化速率(SNiOR)的相关性。通过LEfSe分析筛选出3种温度(25、35和40℃)条件下的生物标记物,并利用KEGG数据库解析了微生物的代谢功能、氮代谢相关酶类型及相对丰度。结果表明,在3种温度条件下,NOB富集系统内的Nitrospira是SNiOR的唯一影响因子,且两者呈现显著的正相关性,Spearman相关系数(ρ)为0.95。此外,在温度条件为25、35和40℃时,分别筛选出14种、6种和8种生物标记物,显著优势菌分别是Xanthomonadales目、Sphingomonadales目和Gammaproteobacteria纲。以上结果说明:温度对微生物的代谢功能具有一定程度的影响,低温有利于碳水化合物代谢、能量代谢;较高温度会抑制能量代谢,NOB丰度降低;随着温度的升高,亚硝酸盐氧化还原酶(nxr)类基因(nxrAnxrB)的数量基本保持不变。
  • 2021年,我国城市污水处理厂污水处理能力为2.1×109 m3·d−1,污水处理率达到97.89%。污水、污泥处理过程中产生并逸散的恶臭气体、挥发性有机物等气态污染物,影响环境空气质量以及污水处理厂内员工、周边居民的身体健康。2002年我国制定了《城镇污水处理厂污染物排放标准》(GB 18918-2002),严格限制硫化氢、氨、甲烷等气态污染物的排放。有研究[1-2]表明,城市污水处理厂的进水区(进水泵站、格栅间)和污泥处理区(污泥浓缩池、污泥脱水机房)是恶臭气体和挥发性有机物的主要逸散源。硫化氢和氨是常见的恶臭物质[3-4],。排放的VOCs包括烷烃、烯烃、芳香烃、卤代烃、含氧有机物、含氮有机物和含硫有机物等80余种[5-6]。以往的研究主要关注恶臭物质中的硫化氢和氨,较少涉及二硫化碳等有机硫化物以及磷化氢。随着我国污水的排放标准的逐步提高,以及人们对环境质量的要求也显著提高,使得城市污水处理厂排放的空气污染物受到越来越多的关注。A2O是一种典型的污水处理工艺,可常用于二级污水处理或三级污水处理以及中水回用,具有良好的脱氮除磷效果,在我国应用广泛。本研究在北方某座采用A2O工艺的城市污水处理厂设置采样点,监测主要处理单元空气中的硫化氢、氨、二硫化碳、磷化氢等恶臭物质,研究了恶臭随季节的变化特征,明确了其排放源及产生原因,评估了恶臭物质的嗅味影响和健康风险,以期为城市污水处理厂恶臭物质的有效削减和控制提供科学依据。

    本研究在京津冀地区某一污水处理厂定期采集样品,研究主要污水处理单元恶臭物质逸散的特征。该污水处理厂采用倒置A2O污水处理工艺,处理规模为1.0×105 m3·d−1,服务区域55.8 km2。污水经过粗格栅、细格栅后经水泵提升进入曝气沉砂池,随后进入生化池,处理后的水全部进入再生水厂深度处理。剩余污泥经离心式脱水机脱水后外运。

    图1所示,在污水处理厂主要处理单元共设置了8个臭味气体的采样点,包括进水(S1)、粗/中格栅(S2)、细格栅(S3)、曝气池进水(S4)、曝气池1(S5)、曝气池2(S6)、脱水间(S7)以及出水(S8)。所有采样点设置在距离地面1.5 m处。采样时间3—6月、11—12月,每月采集样品6~8次。

    图 1  采样点分布图
    Figure 1.  Distribution of sample points

    1)恶臭物质采集及分析。恶臭物质硫化氢、氨、二硫化碳和磷化氢的监测采用在线气体检测仪(PTM-600,中国)。监测时段为 8:30—10:30 am,每个采样点连续监测20 min,获得的恶臭物质的浓度为连续监测数据的平均值。同时采集对应时间和对应工艺段的污水、污泥样品带回实验室分析。

    2)水质及环境条件分析。污水的化学需氧量(COD)和无机离子分别采用多功能消解仪(ET3150B)/分光光度计(DR2800,HACH)和离子色谱仪(DIONEX ICS1000,美国)测定。环境温度和相对湿度采用手持式智能温湿度记录仪(179-TH,USA)监测,风速仪(DeltaOHM HD2303.0,Italy)记录风速。空气温度、相对湿度以及风速分别为6.7~29.1 ℃,23.8%~74.1% 以及0~1.7 m·s−1

    3)计算方法。典范对应分析(canonical correspondence analysis, CCA)方法解析逸散的恶臭气体与环境因素之间的相关性。通过计算气味活性值(odor activity value, OAV),即某种物质的化学浓度与气味检测阈值的比值,确定和分析检出的恶臭物质的嗅觉效应及其气味贡献百分比[7-8]。通过美国环保署公布的空气污染物风险评估方法计算吸入性慢性非致癌风险(hazard quotient,HQ),评估污水厂操作工人恶臭物质的暴露风险[9-10]

    图2所示,监测结果表明,恶臭物质在污水各处理单元的空气中均有检出,逸散水平有明显差异。进水阶段包括进水渠、中格栅、细格栅、曝气池进水等单元。在进水渠空气中检出的氨、硫化氢、二硫化碳和磷化氢的质量浓度分别为0~0.15、0.06~0.33、23.79~99.93和0.05~0.32 mg·m−3。中格栅间内,氨、硫化氢、二硫化碳和磷化氢的质量浓度平均值分别为0.56、1.62、20.88和1.65 mg·m−3。细格栅间空气中氨、硫化氢、二硫化碳和磷化氢的质量浓度分别为0~0.27、0~0.37、0.51~94.33和0~0.33 mg·m−3,平均值分别为0.11、0.16、38.82和0.13 mg·m−3。曝气池进水单元分别有0~0.12 mg·m−3氨、0.33~9.28 mg·m−3硫化氢、0~151.26 mg·m−3二硫化碳和0~1.29 mg·m−3的磷化氢被检出。污水生化处理阶段的采样点包括曝气池1和曝气池2。曝气池1空气中的氨、硫化氢和二硫化碳的质量浓度分别为0~0.11、0~0.14和0~7.14 mg·m−3,平均值分别为0.03、0.02和1.02 mg·m−3。曝气池2空气中只有少量的硫化氢检出,质量浓度为0~0.12 mg·m−3。污泥脱水间空气中氨、硫化氢和二硫化碳等气味物质均有检出,平均浓度分别为1.07、0.03和4.39 mg·m−3

    图 2  污水处理厂恶臭物质的逸散
    Figure 2.  The escape of odorous substances in wastewater treatment plant

    污水处理厂主要处理单元空气中的恶臭气体主要在进水区,平均质量浓度为0.19 mg·m−3(氨),1.07 mg·m−3(硫化氢),44.32 mg·m−3(二硫化碳)和0.58 mg·m−3(磷化氢),而在污水或污泥处理过程产生的较少。以往的恶臭排放特征的调查和研究的结果表明,城市污水处理厂臭味气体的产生源主要是进水池、格栅间、沉砂池、初沉池及污泥处理系统的储泥池、脱水机房等工艺段或构筑物[1, 11-12]

    污水处理厂预处理单元中的恶臭气体是污水在排水管网中长距离输送时产生的。污水在排水管网中长时间停留,溶解氧逐渐消耗,形成缺氧或厌氧环境。此时,厌氧微生物大量繁殖,污水中的含氮、含硫等物质会转化为氨、硫化氢等恶臭物质。有机物在污水管网中的生物转化研究[13-14]发现,从表面到沉积物中1~2 mm深度的沉积层中,硫化氢的质量浓度迅速增加,表明有硫化氢大量产生,同时,该区域的溶解氧值接近0,处于厌氧状态。管道生物膜中生长的硫酸盐还原菌,以乳酸或丙酮酸等有机物作为电子供体,在厌氧状态下,把硫酸盐、亚硫酸盐、硫代硫酸盐等还原为硫化氢。在无氧条件下,污水中的蛋白质,尿素、氨基酸等这些含氮有机物在微生物的作用下转化成氨。含硫、含磷有机物转化为硫醇、硫醚、二硫化碳(C2S)和磷化氢(PH3)等恶臭物质。当污水流经进水泵站和格栅时,排水管网中产生的恶臭气体大量逸出,形成该工艺段的恶臭污染。此外,污水处理厂格栅截留的漂浮物中含有较多的有机物,格栅渣的堆积也会导致有机物的厌氧发酵产生恶臭物质。另外,格栅间是相对密闭的空间,空气流动较少,恶臭物质积聚在格栅间内,导致格栅间空气中的恶臭浓度明显高于其他处理单元。剩余污泥在浓缩、堆置过程中,污泥内部形成厌氧环境,大量有机物进一步发酵、分解,产生NH3和挥发性有机硫等恶臭物质。污泥脱水过程中,恶臭物质从污泥中释放,使污泥处理单元也成为一个恶臭气体的排放源。污水生化处理阶段采用倒置A2O工艺,其中好氧段采用活性污泥法,微孔曝气方式提供生化反应所需的氧气,因此,曝气池内是好氧环境,水中的氨、硫化氢以及有机物等物质在该处理单元被氧化为硝酸盐、硫酸盐以及二氧化碳等物质。并且,生化池位于室外,空气流动加速气体扩散,使生化池的恶臭浓度值最低。

    在污水输送和处理过程中,恶臭气体的产生与排放与污水处理厂的进水水质以及环境条件(如光照强度、空气温度、相对湿度和风速)等因素密切相关[13]。本研究发现进水阶段的各处理单元是污水处理厂恶臭物质的主要逸散源,因此,进一步研究了进水水质(COD、氨氮、总氮、总磷、pH等),对各种检出的恶臭物质逸散水平的影响,并分析其相关性。结果表明,硫化氢的释放与进水中硫酸盐的含量呈显著正相关,这是由于输水管道内呈缺氧环境,污水中的硫酸盐在厌氧菌作用下被还原成硫化物,进一步产生硫化氢。以往的研究也报道了类似的结果[14-15]。此外,空气中氨的浓度与进水中的氨氮和总氮正相关(图3(a))。根据亨利定律,在一定温度下,当液面上的一种气体与溶液中所溶解的该气体达到平衡时,该气体在溶液中的浓度与其在液面上的平衡分压成正比[16]。污水从输送管道中流出进入进水渠、格栅等处理工艺段,气相的分压降低,气水平衡发生改变,水中的氨氮释放到周围空气中。因此,进水氨氮浓度高时,在进水区的空气中会有较多的氨检出。以往的研究表明,COD值高的污水,进入水厂后会释放较多的挥发性有机物。二硫化碳是挥发性有机硫化物,其逸散浓度与污水的COD呈现正相关(图3(b))。磷化氢微溶于水,其水溶液呈弱碱性。当pH增大,水溶液呈碱性时,有更多的磷化氢因溶解度降低而存在于气相中,因此呈现出其与pH呈正相关的现象。

    图 3  恶臭物质逸散的影响因素
    Figure 3.  Factors of malodorous substances emission

    恶臭物质的逸散浓度在各个季节也有明显差异。污水处理厂位于华北地区,每年的3~5月为春季、6~8月为夏季,9~11月为秋季,12~次年2月为冬季。恶臭物质的逸散受温度、相对湿度、风速等环境条件影响。相关性分析显示,氨、硫化氢、C2S 以及PH3的逸散浓度与风速均呈负相关关系,与环境温度、相对湿度呈正相关关系(图3(c))。在温度、相对湿度适宜且平均风速较小的夏季,污水处理各单元空气中恶臭物质的检出较多。

    恶臭物质影响分析包括恶臭物质对人体感官和健康两方面的影响分析。恶臭物质的嗅觉效应及贡献率通过计算其气味活性值(OAV)在总气味活性值中的占比进行评估,以识别主要的气味物质,计算结果如图4所示。主要的致臭物质为硫化氢,其在各工艺段的异味贡献率在78%~100%。在进水阶段的空气中,硫化氢和二硫化碳是引起异味的主要化合物,其异味贡献率分别为78.52%和21.46%。除硫化氢外,二硫化碳是污泥脱水间的主要异味物质,其异味贡献率为10.08%。

    图 4  恶臭物质的臭味贡献率
    Figure 4.  Odor contribution rate of malodorous substances

    氨、硫化氢、二硫化碳以及磷化氢均不是致癌物质,但是,长期接触会导致慢性致病。硫化氢是一种强烈的神经毒素,对粘膜有强烈刺激作用。低浓度的硫化氢对眼、呼吸系统及中枢神经都有影响。硫化氢是细胞色素氧化酶的强抑制剂,能与线粒体内膜呼吸链中的氧化型细胞色素氧化酶中的三价铁离子结合,而抑制电子传递和氧的利用,引起细胞内缺氧,造成细胞内窒息。通过吸入或皮肤接触进入人体后,氨在人体组织内遇水生成氨水,可以溶解组织蛋白质,与脂肪起皂化作用。氨水能破坏体内多种酶的活性,影响组织代谢。CS2是一种多系统亲和毒物,对神经、心血管、胃肠道等系统均有毒害作用。评价人体对这些恶臭物质的暴露风险,有助于相关部门明确污水厂气体污染物控制的优先次序、加强风险管理、保障现场操作人员和周边居民的身体健康。慢性致病风险值平均为0.33(氨),0.039(硫化氢),3.49(二硫化碳)和4.87(磷化氢)。其中,磷化氢的平均风险值最高。根据美国环保署的相关定义,对于慢性毒害,当HQ<1时,其风险可以忽略;当HQ>1时,长期暴露导致的危害应引起关注。磷化氢和二硫化碳是主要的具有潜在慢性致病风险的物质。

    污水处理厂的生化处理池和二沉池等污水处理池均为敞开式,未进行加盖密封处理。在这些工艺段逸散的恶臭气体,不仅在厂区弥漫,还会随风向周边扩散,影响污水厂厂外区域的空气质量。HYSPLIT模型是一种用于气溶胶和气体传输模拟的大气扩散模型,被广泛应用于空气质量监测、环境污染、气候变化等领域。利用HYSPLIT模型分别模拟了在曝气池进水(S4)、曝气池(S5和S6)以及出水(S8)等露天点位,硫化氢、二硫化碳和磷化氢的传输和扩散轨迹。图5反映了3种气体物质释放浓度最大时的模拟结果。硫化氢的扩散方向呈由南向北椭圆形分布(图5(a));10 min后,距离厂界为0.3 km时,其质量浓度降低至1.0×10−7 mg·m−3,影响范围为0.21 km2。二硫化碳和磷化氢的扩散趋势呈西北到东南椭圆形分布(图5(b)和图5(c)),10 min后,二硫化碳的质量浓度为1.1×10−6 mg·m−3 (距离厂界0.286 km),磷化氢的质量浓度为7.9×10−8 mg·m−3 (距离厂界0.25 km)。二硫化碳和磷化氢的影响范围分别为0.20 km2 和 0.14 km2。污水厂内产生的恶臭物质向周边随风扩散10 min之后,距离厂界300 m处空气中的硫化氢和二硫化碳的浓度均低于《恶臭污染物排放标准GB 14554-93》规定的周界恶臭污染物浓度限值。磷化氢慢性致病的风险值也降低至3.13×10−6(远低于1),属于可忽略的范围。

    图 5  气味化合物扩散模拟
    Figure 5.  Diffusion simulation of odor substances

    1)污水处理厂的进水区,尤其中格栅和曝气池进水单元是恶臭物质的主要逸散源,其逸散浓度与进水水质相关。周围空气中恶臭物质的浓度按污水处理单元的顺序依次降低。污水在该污水处理厂经过有效的处理,很少或几乎不释放恶臭物质。在温度、相对湿度适宜且平均风速较小的夏季,恶臭物质的检出相对较多。

    2)硫化氢和二硫化碳是臭味活性物质的主要贡献物质。二硫化碳和磷化氢是主要的具有潜在慢性致病风险的物质。由于空气稀释的作用,污水处理厂检出的恶臭物质对距离厂界300 m以外的区域,造成的影响极小。但对于长期在污水处理厂工作的职工,污染物的累积会增加暴露风险。因此,污水处理厂工作人员在上述污水处理工艺段操作时,需要做好相应的防护。

    3)未来,在污水处理过程中,建议采用适宜的方法如加盖密闭恶臭物质的主要逸散源、建立相应的处理设施等削减和控制污水处理产生的气体污染物。

  • 图 1  不同温度条件下微生物关联网络关系

    Figure 1.  Microbial association network at different temperatures

    图 2  基于分类等级树的组间差异分类单元展示

    Figure 2.  Inter-group difference classification units display based on classification level tree

    图 3  基于线性判别分析(LDA)的组间具有显著差异的分类单元

    Figure 3.  Linear discriminant analysis (LDA) based taxonomic units with significant differences between groups

    图 4  功能基因的KEGG统计

    Figure 4.  KEGG statistics of functional genes

    图 5  不同温度条件下典型代谢功能变化趋势

    Figure 5.  Trend of typical metabolic functions under different temperature conditions

    图 6  系统氮代谢特定物种代谢通路图

    Figure 6.  Map of species-specific metabolic pathways for systemic nitrogen metabolism

    图 7  不同氮转化过程对应的功能基因丰度

    Figure 7.  Functional genes abundance of different nitrogen conversion processes

  • [1] 马智明. 生物脱氮除磷理论与技术进展[J]. 化工管理, 2018(21): 179-180. doi: 10.3969/j.issn.1008-4800.2018.21.137
    [2] PENG Y Z, ZHU G B, GUO J H, et al. Biological nitrogen removal with nitrification and denitrification via nitrite pathway[J]. Applied Microbiology and Biotechnology, 2006, 73(1): 15-26. doi: 10.1007/s00253-006-0534-z
    [3] WANG T, ZHANG L G, LE B, et al. Simultaneous heterotrophic nitrification and aerobic denitrification at high concentrations of NaCl by Halomonas Bacteria[J]. IOP Conference Series:Earth and Environmental Science, 2019, 237: 052033. doi: 10.1088/1755-1315/237/5/052033
    [4] 史银银, 明红霞, 陈泉睿, 等. 亚硝酸盐氧化细菌的生态位以及对海洋环境变化的响应机制[J]. 微生物学通报. 2021, 48(8) : 2985-2901.
    [5] DAIMS H, LÜCKER S, WAGNER M, et al. A new perspective on microbes formerly known as nitrite-oxidizing bacteria[J]. Trends in Microbiology, 2016, 24(9): 699-712. doi: 10.1016/j.tim.2016.05.004
    [6] GOUDAR C T, GANJI S H, PUJAR B G. Substrate inhibition kinetics of phenol biodegradation[J]. Water Environment Research, 2000, 72(1): 50-55. doi: 10.2175/106143000X137103
    [7] SPIECK E, HARTWIG C, CORMACK I, et al. Selective enrichment and molecular characterization of a previously uncultured Nitrospira-like bacterium from activated sludge[J]. Environmental Microbiology, 2006, 8: 405-415. doi: 10.1111/j.1462-2920.2005.00905.x
    [8] 包鹏, 王淑莹, 马斌, 等. 不同溶解氧间歇曝气对亚硝酸盐氧化菌的影响[J]. 中国环境科学, 2016, 36(9): 2696-2702. doi: 10.3969/j.issn.1000-6923.2016.09.024
    [9] ZHU S M, CHEN S. The impact of temperature on nitrification rate in fixed film biofilters[J]. Aquacultural Engineering, 2002, 26(4): 221-237. doi: 10.1016/S0144-8609(02)00022-5
    [10] 陈翠忠, 高宇学, 王文迪, 等. 温度对SBR生物脱氮效能的长期影响[J]. 环境工程, 2018, 36(6): 68-72.
    [11] TAYLOR A E, MYROLD D D, BOTTOMLEY P J. Temperature affects the kinetics of nitrite oxidation and nitrification coupling in four agricultural soils[J]. Soil Biology and Biochemistry, 2019, 136: 107523-107523. doi: 10.1016/j.soilbio.2019.107523
    [12] 于雪, 孙洪伟, 李维维, 等. 温度对硝化杆菌(Nitrobacter)活性动力学影响[J]. 中国环境科学, 2019, 40(3): 1427-1430.
    [13] 吴鹏, 陆爽君, 徐乐中, 等. 温度对ABR-MBR复合工艺处理生活污水的影响及其微生物群落分析[J]. 环境科学, 2014, 35(9): 3466-3472.
    [14] WERTZ S, POLY F, LE ROUX X, et al. Development and application of a PCR-denaturing gradient gel electrophoresis tool to study the diversity of Nitrobacter-like nxrA sequences in soil[J]. FEMS Microbiology Ecology, 2008, 63(2): 261-271. doi: 10.1111/j.1574-6941.2007.00416.x
    [15] CHU Y J, COREY D R. RNA sequencing: Platform selection, experimental design, and data interpretation[J]. Nucleic Acid Therapeutics, 2012, 22(4): 271-274. doi: 10.1089/nat.2012.0367
    [16] 曾涛涛, 蒋小梅, 韩科昌, 等. 生活污水处理厂微生物群落结构解析[J]. 安全与环境学报, 2018, 18(2): 697-703.
    [17] FRIAS-LOPEZ J, SHI Y, TYSON G W, et al. Microbial community gene expression in ocean surface waters[J]. Proceedings of the National Academy of Sciences of the United States of America, 2008, 105(10): 3805-3810. doi: 10.1073/pnas.0708897105
    [18] 国家环境保护总局. 水和废水监测分析方法[M]. 4版. 北京: 中国环境科学出版社, 2002: 211-213
    [19] HEWSON I, PORETSKY R S, TRIPP H J, et al. Spatial patterns and light-driven variation of microbial population gene expression in surface waters of the oligotrophic open ocean[J]. Environmental Microbiology, 2010, 12(7): 1940-1956. doi: 10.1111/j.1462-2920.2010.02198.x
    [20] 张婷婷, 张建, 杨芳, 等. 温度对污水脱氮系统污染物去除效果及氧化亚氮释放的影响[J]. 环境科学, 2012, 33(4): 1283-1287.
    [21] LANGILLE M G I, ZANEVELD J, CAPORASO J G, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences[J]. Nature Biotechnology, 2013, 31(9): 814-821. doi: 10.1038/nbt.2676
    [22] WANG Q, YE L, JIANG G, et al. Side-stream sludge treatment using free nitrous acid selectively eliminates nitrite oxidizing bacteria and achieves the nitrite pathway[J]. Water Research, 2014, 55(10): 245-255.
    [23] MA B, YANG L, WANG Q, et al. Inactivation and adaptation of ammonia-oxidizing bacteria and nitrite-oxidizing bacteria when exposed to free nitrous acid[J]. Bioresource Technology, 2017, 245(8): 1266-1270.
    [24] RADAX R, RATTEI T, LANZEN A, et al. Metatranscriptomics of the marine sponge Geodia barretti: tackling phylogeny and function of its microbial community[J]. Environmental Microbiology, 2012, 14(5): 1308-1324. doi: 10.1111/j.1462-2920.2012.02714.x
    [25] LI Y, LI J, WANG Z, et al. Inhibition kinetics of nitritation and half-nitritation of old landfill leachate in a membrane bioreactor[J]. Journal of Bioscience and Bioengineering, 2016, 123(4): 482-488.
    [26] 杨浩, 张国珍, 杨晓妮, 等. 典型集雨人饮地区窖水微生物群落多样性及差异解析[J]. 环境科学, 2017, 38(4): 4733-4746.
  • 加载中
图( 7)
计量
  • 文章访问数:  4713
  • HTML全文浏览数:  4713
  • PDF下载数:  93
  • 施引文献:  0
出版历程
  • 收稿日期:  2021-04-12
  • 录用日期:  2021-08-16
  • 刊出日期:  2022-03-10
李奕燃, 董昆, 于雪, 李晓强, 王光杰, 孙洪伟. 利用16S rRNA高通量测序技术考察温度对生物脱氮硝化过程中亚硝酸盐氧化菌代谢功能的影响[J]. 环境工程学报, 2022, 16(3): 980-988. doi: 10.12030/j.cjee.202104086
引用本文: 李奕燃, 董昆, 于雪, 李晓强, 王光杰, 孙洪伟. 利用16S rRNA高通量测序技术考察温度对生物脱氮硝化过程中亚硝酸盐氧化菌代谢功能的影响[J]. 环境工程学报, 2022, 16(3): 980-988. doi: 10.12030/j.cjee.202104086
LI Yiran, DONG Kun, YU Xue, LI Xiaoqiang, WANG Guanjie, SUN Hongwei. Investigation on the effect of temperature on the metabolism and function of nitrite-oxidizing bacteria in the process of biological nitrification by 16S rRNA high-throughput sequencing technology[J]. Chinese Journal of Environmental Engineering, 2022, 16(3): 980-988. doi: 10.12030/j.cjee.202104086
Citation: LI Yiran, DONG Kun, YU Xue, LI Xiaoqiang, WANG Guanjie, SUN Hongwei. Investigation on the effect of temperature on the metabolism and function of nitrite-oxidizing bacteria in the process of biological nitrification by 16S rRNA high-throughput sequencing technology[J]. Chinese Journal of Environmental Engineering, 2022, 16(3): 980-988. doi: 10.12030/j.cjee.202104086

利用16S rRNA高通量测序技术考察温度对生物脱氮硝化过程中亚硝酸盐氧化菌代谢功能的影响

    通讯作者: 孙洪伟(1978—),男,博士,教授,sunhw@ytu.edu.cn
    作者简介: 李奕燃 (1995—),男,硕士研究生,1587381875@qq.com
  • 1. 兰州交通大学环境与市政工程学院,兰州 730070
  • 2. 烟台大学环境与材料工程学院,烟台 264000
  • 3. 山东同济测试科技股份有限公司,烟台 264000
基金项目:
省级大学生创新创业训练计划项目(2017107320315)

摘要: 为探究温度对生物脱氮硝化过程微生物代谢及功能基因的影响,采用16S rRNA高通量测序技术,以富集的亚硝酸盐氧化菌(NOB)为研究对象,考察了活性污泥系统硝化过程中菌属与比亚硝酸盐氧化速率(SNiOR)的相关性。通过LEfSe分析筛选出3种温度(25、35和40℃)条件下的生物标记物,并利用KEGG数据库解析了微生物的代谢功能、氮代谢相关酶类型及相对丰度。结果表明,在3种温度条件下,NOB富集系统内的Nitrospira是SNiOR的唯一影响因子,且两者呈现显著的正相关性,Spearman相关系数(ρ)为0.95。此外,在温度条件为25、35和40℃时,分别筛选出14种、6种和8种生物标记物,显著优势菌分别是Xanthomonadales目、Sphingomonadales目和Gammaproteobacteria纲。以上结果说明:温度对微生物的代谢功能具有一定程度的影响,低温有利于碳水化合物代谢、能量代谢;较高温度会抑制能量代谢,NOB丰度降低;随着温度的升高,亚硝酸盐氧化还原酶(nxr)类基因(nxrAnxrB)的数量基本保持不变。

English Abstract

  • 在短程硝化反硝化、同步硝化反硝化、厌氧氨氧化等[1]非传统生物脱氮工艺的生物脱氮系统中,氨化菌会在好氧条件下将有机氮转化为氨氮(NH+4-N),然后氨氧化菌(ammonia-oxidizing bacteria,AOB)和亚硝酸盐氧化菌(nitrite-oxidizing bacteria,NOB)会协同作用将氨氮(NH+4-N)转化为硝态氮(NO3-N)。该协同作用的过程即AOB先将NH+4-N氧化为亚硝态氮(NO2-N),NOB再将NO2-N进一步氧化为NO3-N[2]。其中,NOB为革兰氏阴性菌,是自养微生物,其适宜的生长温度为25~30 ℃[3]。在好氧条件下,NOB以CO2为碳源参与硝化反应,以氧气和硝酸盐为电子受体进行亚硝酸盐的氧化。此外,NOB还可利用其他能源参与微生物的代谢活动,如NOB作为专性氧化菌以甲酸、硫化物或其他有机化合物为底物发生代谢反应[4]。根据微生物的形态特征和细胞内膜的排列方式,NOB可划分为7个属:硝化杆菌属(Nitrobacter)、硝化球菌属(Nitrococcus)、硝化刺菌属(Nitrospina)、硝化螺菌属(Nitrospira)、Nitrotoga属、Candidatus Nitromaritima属和Nitrolancea[5]。其中,NitrobacterNitrospira是污水生物处理系统中的主导NOB[6]。在高溶解氧、高底物浓度条件下,Nitrobacter含量远高于Nitrospira,而在低溶解氧和低基质底物浓度条件下,Nitrospira为优势菌属[7-8]

    温度是硝化反应过程重要的影响因素之一。以溶解氧(dissolved oxygen,DO)和总氨氮(total ammonia nitrogen,TAN)为限制条件,当温度每升高1℃(在20 ℃基础上)时,硝化速率可分别提高1.108%和4.275%[9]。陈翠忠等[10]采用缺氧-好氧SBR处理合成废水时发现,当温度从15 ℃升至25 ℃时,氨氮去除率随之升高;当温度从25 ℃升至35 ℃时,氨氮去除率却逐渐降低。TAYLOR等[11]在测定NO2氧化动力学参数时发现,反应温度为30 ℃时的最大反应速率(Vmax)大于17 ℃时的Vmax。于雪等[12-13]发现:当温度较低时,微生物细胞膜呈凝胶状,水分和营养物质跨膜运输受阻,此时的细胞因缺乏营养而停止生命活动,反应速率随之降低;随着温度的升高,生化反应速率逐渐恢复并加快,且当温度超过一定范围后,细菌体内的蛋白质和核酸等物质对温度敏感的细胞组分发生变性,并使得其生长停止、最终死亡。

    近年来,16S rRNA基因Illumina MiSeq高通量测序技术被广泛应用于海洋、湖泊、水库、污水处理系统等水体微生物多样性、群落结构及功能研究领域[14-15]。众多学者基于该技术研究了温度对硝化菌群组成的影响,发现温度会影响NOB菌群的多样性和结构[16-19],但鲜有研究涉及NOB菌群的微生物代谢功能及相关功能基因。

    本研究基于16S rRNA高通量测序技术,探究在3种温度(25 ℃、35 ℃和40 ℃)条件下,人工模拟废水中富集的NOB菌群间的关联网络关系、生物特征标记物及代谢功能,并利用KEGG数据库分析硝化过程微生物代谢途径,以期进一步揭示温度对硝化过程中的微生物的影响,为生物脱氮工艺的优化调控提供参考。

    • 在实验采用的人工模拟废水中,NO2-N的质量浓度通过投加NaNO2来调节(5.0~140 mg·L−1);初始pH利用浓度为1 mol·L−1的盐酸溶液和浓度为1 mol·L−1的氢氧化钠溶液来调节。为满足微生物自身增殖的需要,向每升模拟废水中添加的微量元素为:0.4 g NaHCO3;1 g KH2PO4;1.31 g K2HPO4;1.25 g EDTA;0.55 g ZnSO4·7H2O;0.40 g CoCl2·6H2O;1.275 g MnCl2·4H2O;0.40 g CuSO4·5H2O;0.05 g Na2MoO4·2H2O;1.375 g CaCl2·2H2O;1.25 g FeCl3·6H2O;44.4 g MgSO4·7HO。

      实验接种污泥取自兰州市西固区兰炼污水处理厂。该污泥具有良好的生物脱氮性能,其混合液挥发性悬浮固体浓度(MLVSS)为3 300 mg·L−1

      实验装置为SBR(有效容积为3 L,排水比为1∶3)。反应器的进出水NO2-N质量浓度及MLVSS采用标准方法测定[20-21]。采用WTW Multi-3420多参数测定仪对温度、pH及DO进行实时监测。

    • 在实验开始前对接种污泥进行了为期30 d的培养驯化,以获得稳定的硝化性能。在SBR内,不同阶段泥水混合液的温度分别控制为25 ℃(阶段Ⅰ)、35 ℃(阶段Ⅱ)和40 ℃(阶段Ⅲ)。每个阶段SBR分别运行100个周期,以获得相应温度条件下成熟的活性污泥。SBR的各运行周期由持续进水(1 min)、硝化(通过实时监测pH和DO控制硝化终点,硝化的同时曝气)、沉淀(30 min)、排水(5 min)和闲置5个阶段组成。此外,为避免pH和DO成为亚硝酸盐氧化的限制因素,pH维持在(7.5±0.1),DO维持在4.0~4.5 mg·L−1。利用比亚硝态氮氧化速率(SNiOR)来反映NOB的活性,其计算公式为式(1)。

      式中:SNiOR为比亚硝态氮氧化速率,g·(g·d)−1(以每日每克可挥发性悬浮物氧化的亚硝态氮的质量计);c0(NO2-N)为曝气开始时NO2-N的质量浓度,mg·L−1ct(NO2-N)为曝气结束时NO2-N的质量浓度,mg·L−1;MLVSS为污泥的质量浓度,mg·L−1tN为曝气时间,min[22]

    • 在25 ℃(第99周期)、35 ℃(第199周期)和40 ℃(第299周期)时,取SBR中的活性污泥样品(每个温度包括3个平行样品)并提取其DNA。将样品中提取到的DNA,对16S rRNA V3-V4区进行扩增,以活性污泥样本所提取到的DNA原液为PCR模板,通用引物序列为:F:520F(5-ACTCCTACGGGCAGCA-3)R:802R:(5-GGACTACHVGGGTWTCTAAT-3)。PCR反应条件为:98 ℃时预变性2 min;98 ℃时变性15 s;55 ℃时退火30 s;72 ℃时延伸30 s后再保温5 min,提取后于4 ℃下保存。委托上海派森诺生物科技有限公司进行Illumina MiSeq高通量测序,通过16S rRNA高通量进行测序分析并获得微生物的结构组成。在此基础上,利用KEGG数据库,借助PICRUSt工具预测分析微生物菌群的代谢功能,并对氮转化的相关功能基因进行比对。其中,相关性分析基于SPSS统计分析软件完成,LEfSe分析借助派森诺云分析平台完成。

    • 基于微生物的关联网络图可分析系统微生物群落之间、微生物与环境因子之间的相互关系。通过关联网络分析的方法,可探究系统内与微生物种群相互“协作”或“竞争”的模式。图1为不同温度条件下各分类水平(门phylum、纲class、目order、科family和属genus)的微生物种群与SNiOR的关联网络特征。图中相同颜色的节点为同一门,且节点越大,该物种在微生物类群中所占的相对丰度越大。而与节点相连的线越多,则表明该物种与其他微生物联系越密切。两节点间的连线表示微生物两两之间存在相关性(Spearman相关系数)。其中,红线表示正相关,绿线表示负相关。在各分类水平上,各微生物间具有错综复杂的交互关系,且呈正相关的微生物多于呈负相关的微生物。此外,在亚硝酸盐氧化过程中,基质底物为亚硝酸盐,故SNiOR是NOB活性的直接体现。关联网络图表明,SNiOR与硝化螺菌门Nitrospirae硝化螺菌纲Nitrospirae鞘脂杆菌目sphingobacterials、鞘脂杆菌科sphingobacterials及硝化螺菌属Nitrospira均表现出正相关性。Nitrospira是污水生化处理系统中NOB菌群的代表菌属[23-25]。在本研究中,SNiOR与Nitrospira呈现正相关关系(Spearman相关系数ρ为0.95)。此外,对于生化系统内的Nitrospira,与其呈正相关的分别为黑罗河杆菌属(Rhodanobacter)(ρ=0.88)和独岛菌属(Dokdonella)(ρ=0.88),而未发现与Nitrospira呈负相关关系的微生物菌属。在各个分类水平上,SNiOR仅与TM6门呈负相关关系(ρ=−0.87)。这表明,在生化系统中亚硝酸盐的氧化是基于多种微生物协同作用来完成的,且较多的微生物种群间以相互“协作”作用共存,仅较少的微生物种群间存在相互“竞争”的生化反应。

    • LefSe是一种基于线性判别分析(linear discriminant analysis,LDA)效应量的分析方法,其本质是将线性判别分析与非参数的Kruskal-Wallis及Wilcoxon秩和检验相结合,从而筛选关键生物标记物(即关键群落成员)。在3种温度条件下,采用LefSe筛选出不同温度条件下污泥样品的关键生物,结果如图2所示。活性污泥样品的分类等级树展示了微生物群落中从门到属(从内圈到外圈依次排列)所有微生物的等级关系。其中,节点大小对应于该物种的平均相对丰度,黄色圆点代表未体现出显著组间差异的微生物物种,而其他颜色(如绿色和红色)则表明这些物种呈现显著组间差异,且在该色所代表分组样品中丰度较高。

      此外,通过线性判别分析(LDA值)方法进行了生物标记分析,以获得不同温度样品间具有显著差异的菌群。图3为在3种温度条件下污泥样品组间LDA值大于3.0的具有显著差异的微生物种群,柱状图的长度代表差异物种的显著性(即为LDA值)。在25℃、35℃和40℃温度条件下,共筛选出28种具有显著差异的微生物标记物。每种温度条件下的生物标记物分别有14种、6种和8种。以35℃活性污泥样品为例,SBR生化系统中存在Sphingomonadales(4.1%)、Rhodanobacter(4.3%)、Sphingomonadaceae(3.5%)、Novosphingobium(2.7%)、Zoogloea(1.8%)、Sulfuricella(0.02%)等6种生物标记物。在该温度下,这些微生物在系统内具有相对较高的丰度,是微生物种群的优势菌属。

    • 将上述研究得到的基因组数据与KEGG数据库进行了比对,以推测系统中微生物代谢通路并分析其代谢功能。在分类水平上,KEGG数据库包括4个等级:第一等级为生物代谢通路;第二等级为代谢通路的子功能;第三等级为代谢通路图;第四等级为KO 编号(对应代谢通路上各个KO的具体注释信息)[26]图4为基于KEGG数据库的3种温度条件下活性污泥微生物代谢功能分析结果。通过数据比对获得了7类生物代谢通路功能,即细胞进程、环境信息处理、遗传信息处理、人类疾病、新陈代谢、生物体系统和未分类这7个功能模块。在3种温度下,KEGG一级代谢通路中7个功能模块丰度变化趋势一致。即代谢功能模块在KEGG一级代谢通路中基因丰度最高,环境信息处理、遗传信息处理以及未分类的功能模块相应的基因丰度次之,生物体系统和人类疾病模块相关基因丰度较低。在代谢功能模块中,氨基酸代谢、碳水化合物代谢和能量代谢是最大的功能类群,这三者分别占总生物代谢通路的10.2%、9.9%和6.3%(在3种温度下取平均值,以下同理);其次是辅助因子和维生素的代谢、脂质代谢和核苷酸代谢等;环境信息处理模块中膜运输占总生物代谢通路的10.5%,信号转导占2.4%;遗传信息处理模块中,对遗传信息的复制和修复、翻译占据主要优势,其相对丰度分别为7.1%和4.6%;未分类功能模块中,细胞加工和信号传递、遗传信息处理以及代谢功能分别占总生物代谢通路的4.0%、2.5%和2.5%。

      以上结果表明,低温有利于微生物对碳水化合物的代谢、能量代谢;较高温度则会抑制能量代谢,即影响微生物的新陈代谢。因此,应避免在高温下运行以硝化反应占主导的脱氮工艺系统,以避免降低其中微生物的活性。

      为进一步研究温度对微生物在脱氮工艺中代谢功能的影响,本研究建立了KEGG二级代谢模型(图5),即在3种温度下6种代谢功能的变化规律。随着温度的升高(25 ℃→35 ℃→40 ℃),系统中微生物氨基酸代谢、辅助因子和维生素的代谢、脂质代谢和核苷酸代谢的相对丰度变化较小,而对碳水化合物及能量的代谢在35 ℃升至40 ℃的过程中呈下降趋势。此外,NOB氧化亚硝酸盐的过程属于能量代谢,且随着温度升高,能量代谢的相对丰度表现为先升高、后降低的趋势。

      在3种温度条件下,一级代谢通路中7个功能模块丰度变化趋势一致,代谢功能模块基因丰度最高。但在不同温度条件下,KEGG二级代谢相对丰度发生变化。NOB 降解NO2-N的代谢功能随着温度的升高,呈现先升高、后降低的趋势。KO一级水平层次上的7类代谢功能和KO二级水平层次上的 41 类二级代谢功能均具有显著的丰度差异。这种丰度差异在编码新陈代谢的基因序列模块中体现尤为突出,表明在活性污泥样品中存在大量具有特定生态功能的微生物菌群,其新陈代谢活动较为活跃。

    • 基于16S rRNA基因测序结果和PICRUSt工具可获得KEGG直系同源基因簇(KO)的具体注释信息。将注释得到的所有KO值映射到KEGG代谢通路数据库的“氮代谢通路图”中,可比对出与氮代谢相关的功能基因及其相对丰度。温度批次实验的所有样品均存在5种氮代谢途径(硝化、异化硝酸盐还原、同化硝酸盐还原、反硝化和固氮)和24个功能基因(图6)。硝化过程相关的功能酶为氨单加氧酶(amoABC)和羟胺氧化酶(hao)/还原酶(hcp)及亚硝酸盐氧化还原酶(nxrAnxrB)。在3种温度下,亚硝酸盐氧化还原酶(nxrAnxrB)丰度(9.9%、10.5%、10.3%)均远远高于氨单加氧酶(AmoCAB)和羟胺脱氢酶(hao)丰度总和(0.3%、0.3%、0.1%)。这表明亚硝酸盐氧化作用在系统硝化过程中占据主导地位。

      对比在3种温度条件下氮转化功能基因的丰度可发现,相关酶存在于活性污泥样品的氮转化过程中,但其功能基因的相对丰度存在一定差异。5种氮转化过程对应的功能基因如图7所示。在25 ℃、35 ℃和40 ℃条件下,亚硝酸盐氧化还原酶(nxr)类基因(nxrAnxrB)相对丰度分别为9.9%、10.5%和10.3%,随着温度的升高,基因数量基本保持不变。此外,在硝化过程中存在催化氨为羟胺的氨单加氧酶类基因(AmoABC)、将羟胺转化为亚硝酸盐的羟胺脱氢酶(hao),以及将亚硝酸盐氧化为硝酸盐的亚硝酸盐氧化还原酶(nxrAnxrB),其中nxrAnxrB分别为NOB菌属中硝化杆菌属(Nitrobactor)和硝化螺菌属(Nitrospira)的功能基因。这2种菌属为硝化过程中的关键菌属,可将NO2-N氧化为NO3-N。硝化过程参与的酶相对丰度大小为:nxrB>nxrA>hao>AmoCAB。这表明系统中存在大量将NO2-N氧化为NO3-N的基因,亦说明整个系统内的亚硝酸盐氧化效果显著。

    • 1)在3种温度条件下,多种生物的协同作用是实现亚硝酸盐氧化的生物反应机制。SNiOR仅与Nitrospira呈现正相关性(Spearman系数ρ=0.95),而Nitrospira相对丰度是SNiOR的唯一影响因子。通过LefSe分析表明,微生物种群的群落结构也存在显著差异,分别获得了14种、6种和8种具有显著差异(P<0.05)的微生物标记物。

      2)温度对微生物的代谢功能有一定程度的影响。碳水化合物代谢、能量代谢、氨基酸代谢、核苷酸代谢及酶类代谢共存于活性污泥系统中,而能量代谢、核苷酸代谢和酶类代谢是微生物新陈代谢的主导功能模块,且低温有利于碳水化合物代谢、能量代谢。较高温度会抑制能量代谢,NOB丰度降低。

      3)温度对氮代谢有关基因丰度的影响并不显著。随着温度的升高,亚硝酸盐氧化还原酶(nxr)类基因(nxrAnxrB)的数量基本保持不变。此外,对硝化过程的基因统计结果表明,具有亚硝酸盐氧化功能的基因相对丰度远大于具有氨氧化作用的基因相对丰度,系统的亚硝酸盐氧化效果显著。

    参考文献 (26)

返回顶部

目录

/

返回文章
返回