-
食品加工、制药工业、皮革工业、垃圾渗滤液等废水不仅氨氮 (NH4+-N) 浓度高,还含有大量盐分[1-3]。生物法虽然作为一种经济有效的废水处理方法,但是高盐分可抑制微生物酶活,甚至引起细胞质壁分离导致微生物死亡[4-5],不利于微生物对NH4+-N的降解。因此,在对高盐废水的脱氮处理中,通常使用膜分离[6]、蒸发结晶[7]、芬顿[8]和电渗析技术[9]等物化法进行脱盐预处理,这不但使运行成本高、工艺流程长,还会造成二次污染。虽然经物化脱盐处理后的高盐含氮废水可进一步通过生物法脱氮,但去除效果并不稳定[10-11],即使在低盐度环境(<2%),传统硝化、反硝化菌和常规异养菌等常见脱氮微生物的生长代谢及活性也会受到明显抑制[12-13]。具备异养硝化-好氧反硝化(heterotrophic nitrification-aerobic denitrification, HN-AD)功能的好氧脱氮菌在好氧条件下因具备同时硝化反硝化、且无中间产物亚硝态氮 (NO2−-N) 和硝态氮 (NO3−-N) 积累的特质得到广泛关注[14]。不仅如此,部分HN-AD菌属还兼具耐盐性,已被广泛应用于高盐废水处理过程,如不动杆菌属 (Acinetobacter) 、盐单胞菌属 (Halomonas) 、副球菌属 (Paracoccus) [15]等。然而,在工程应用中,许多工业废水的盐度波动大,这使得包括HN-AD菌在内的众多微生物难以快速适应不同的盐度环境,菌群无法实现高效富集,进而影响废水处理效果[16-17]。为了增强微生物对高盐胁迫的抵抗能力,已有研究[18]通过增加盐浓度驯化功能菌属耐盐性能的方式以期获得高去污效率、高富集丰度的耐盐型脱氮功能微生物。然而,在实际应用中却发现,即使经过长期驯化,功能微生物对高盐胁迫的抵抗能力仍较弱[19-20],原因在于功能微生物易流失、难截留导致富集丰度较低,最终影响废水处理效果的稳定性[19, 21-22]。因此,为解决这一问题,亟待研发一种即可提高微生物对高盐胁迫的抵抗能力,又可提高功能菌属生长富集能力的新技术。
微生物在逆境中可通过积累或从环境中吸收相容性溶质 (compatible solutes, CS) 来抵抗外部条件带来的胁迫作用,在不需要细胞内酶的任何特殊适应下即可保护细胞免受不可逆转的损害[23-24]。所以CS在细胞中发挥了重要作用,它不仅不会阻碍正常的细胞代谢,还可以平衡细胞质内的渗透压以增强微生物的活性[23, 25]。在VYRIDES等[26]的研究中,盐度为4%时,假单胞菌株 (Pseudomonas sp strain ADP) 产生的相容性溶质主要为海藻糖;ZHANG等[27]在使用颗粒污泥处理含盐养殖废水的过程中,添加1 000 μmol·L−1的海藻糖,可使产甲烷活性和磷酸酶活性分别提高22%和27%。而海藻糖作为一种价格低、易采购的典型相容性溶质,在高盐废水处理中已被成功应用[28-29]。其作为渗透调节剂易被微生物吸收利用,通过调节细胞渗透压帮助微生物快速适应高盐环境,提升去污能力[26]。目前,已有海藻糖用于高盐废水的厌氧处理系统的相关报道[29-30]。BAI等[30]研究了在序批式反应器中添加海藻糖考察了厌氧氨氧化细菌处理富氮含盐废水的长期脱氮性能,在添加浓度为300 μmol·L−1时,就可使NH4+-N、NO2−-N的去除率分别提高为32.4%和42.2%;杨振琳等[28]通过添加海藻糖强化厌氧化耦合反硝化工艺,在盐度为3.2%的模拟废水中加入250 μmol·L−1的海藻糖,NH4+-N和NO2−-N的去除速率分别提高了81.25%和75%。可见海藻糖对提高长期脱氮效果有显著性作用,但其在高盐废水好氧生物处理中的研究却鲜有报道。盐度显著影响微生物的物理和生化特性,长期作用甚至能够改变整个环境的微生物群落结构[31]。
本研究利用膜曝气生物膜反应器 (membrane aerobic biofilm reactor, MABR) 富集好氧脱氮菌,通过调节外加海藻糖的浓度,采用高通量测序等分析方法,从微生物脱氮性能、群落结构及脱氮功能基因丰度等层面,研究海藻糖浓度梯度对高盐废水中微生物的物种丰度、菌群多样性及脱氮功能基因的影响,旨在为推动好氧脱氮菌在高盐废水好氧生物脱氮处理的实际应用提供理论基础及新技术思路。
-
在实验过程中,建立了5个MABR,工作体积为0.85 L,反应装置见图1。反应器的膜组件由128根中空纤维膜组成,用于空气运输和生物膜附着,中空纤维是用复合聚合物材料制造,只允许气体分子渗透。纤维的外径为300 μm,内径为200 μm,膜孔径为0.1 μm,有效比表面积为1.76 cm2·cm−3。膜组件一端通过流量计和压力表连接到空气泵,另一端连接到底端曝气管形成氧循环,曝气压力设定为10 kPa,进出水均采用蠕动泵完成,水力停留时间为48 h[32]。反应器顶部留有取样口,用于采样和溶解氧 (DO) 监测等操作。
-
合成的高盐高氮废水由无水CH3COONa (碳源) 、(NH4)2SO4 (氮源) 、NaCl和微量元素组成。进水COD为5 000 mg·L−1,NH4+-N为500 mg·L−1,盐度为2.5%,微量元素为50 mL·L−1。微量元素配方: MgSO4·7H2O 2.0 g·L−1,MnSO4·H2O 0.1 g·L−1,CaCl2 1.5 g·L−1,FeSO4·7H2O 0.1 g·L−1,K2HPO4 5.0 g·L−1。
该研究的接种菌液中主要为产碱杆菌(Alcaligenes)、不动杆菌(Acinetobacter sp.TAC-1)和假单胞菌 (Pseudomonas) 等具备HN-AD功能的好氧脱氮混合菌群。这些细菌主要分离于猪场沼液、垃圾渗滤液和化工废水,可耐受高浓度氨氮。为确保适当的菌液浓度和脱氮活性,首先进行菌液驯化,再分别接种于5个反应器中进行挂膜。挂膜完成后,分别加入不同浓度海藻糖 (0、40、120、360和1 080 μmol·L−1) 进行实验,对应反应器分别为C0、C40、C120、C360和C1 080。
-
反应器启动运行阶段的溶氧、盐度及pH值分别采用便携式溶氧仪、海水比重计及pH计进行直接测量。菌液浓度 (OD600 nm) 采用紫外分光光度法测定,进、出水NH4+-N、NO3−-N、NO2−-N、TN和COD的浓度均参照文献中的方法[33]进行检测。
利用扫描电镜 (scanning electron microscope, SEM) 技术表征反应器挂膜阶段的生物膜表观形态特征。将采集的生物膜样品经过固定、PBS清洗、乙醇梯度脱水 (30%、50%、70%、85%和95%质量分数的乙醇洗脱1次,100%乙醇洗脱2次(每次脱水时间为15~20 min)) 、冷冻干燥等步骤完成样品预处理,将制备完成的样品送至武汉铄思百检测技术有限公司并进行SEM观察。
-
在各反应器稳定运行阶段,采集微生物样品(每个反应器重复采集3次样品,每个样品约0.5 g),将样品短期保存于−80 ℃冰箱中。采用通用检测试剂盒 (MOBIO实验室,美国) 提取样品的总基因组DNA,采用引物338F和806R扩增V3~V4区域的l6S rRNA基因,测序引物序列为338F(5'-ACTCCTACGGGAGGCAGCA-3')和806R(5'-GGACTACHVGGGTWTCTAAT-3')。采用Trimmomatic和Flash软件对测序序列进行拼接,通过Usearch 软件对具有97%相似性的序列进行OTU聚类,经Usearch_global 获得OTU序列丰度统计表,进行后续微生物群落多样性分析;基于京都基因和基因组百科全书(http://www.genome.jp·KEGG)数据库及PICRUSt1软件,对样品中微生物脱氮功能基因进行预测分析。
-
SEM表征结果显示,反应器启动运行20 d,膜丝上明显附着了大量的以杆菌和球菌为主的微生物 (图2) ,说明反应器于20 d内成功完成挂膜启动。该研究反应器菌液挂膜时间较文献报道相关结果[32-34]略有增加,这可能是因为盐度抑制了微生物的生长繁殖,从而导致挂膜周期延长。反应器启动运行阶段的污染物去除情况如图3(a)~图3(e)中第1阶段所示。5个平行反应器的NH4+-N、TN及COD去除率分别为32%~45%、32%~40%和60%~71%,无明显有害中间产物NO2−-N和NO3−-N的积累。启动阶段反应器脱氮性能较低,究其原因主要在于高盐环境对反应器富集的脱氮微生物的生长及活性存在影响。高盐度会使微生物细胞处于一个高渗环境,影响细胞正常生长代谢,进而降低污染物的去除率[35]。
在4个平行反应器中添加不同浓度海藻糖 (40、120、360和1 080 μmol·L−1) ,设定为实验组 (C40、C120、C360和C1 080) ,对照组 (C0) 不添加海藻糖,连续运行80 d结束。添加海藻糖后的反应器运行性能如图3(a)~图3(e)中第2阶段所示。相较于对照组,4个实验组中NH4+-N去除率分别提高了10.70%、32.72%、27.36%、19.45%,TN去除率分别提高了8.32%、28.36%、22.53%、17.63%,COD去除率分别提高了12.09%、31.14%、25.27%、25.06%。这说明添加海藻糖对高盐给微生物带来的高渗胁迫具有缓解作用,提高了系统脱氮性能。C120组的NH4+-N、TN及COD去除率最高,C360组的TN去除率及C1 080组的NH4+-N、TN去除率相较C120组降低,3组中COD去除率无明显变化。可以看出,外源海藻糖浓度与污染物去除率无线性相关性,最适海藻糖浓度为120 μmol·L−1。这主要是因为微生物对海藻糖的最大渗透调节剂有一定的吸收范围[28],当体系内海藻糖浓度超过微生物作为渗透调节剂的吸收浓度范围时,过剩的海藻糖可充当常规异养反硝化菌 (heterotrophic denitrified bacteria, HDB) 的有机碳源[28]。因此,在本研究中,高浓度海藻糖 (360、1 080 μmol·L−1) 条件可能对体系内HDB生长更有利,进而影响好氧脱氮菌群的富集与活性,降低了反应器脱氮性能。
-
1) 微生物α多样性分析。通过微生物群落分析研究不同海藻糖浓度对反应器中物种丰富度及多样性的影响,结果如表1所示。通过MiSeq Illumina测序分析获得269~303个OTU,覆盖指数均大于99%,说明测序深度包含测序样品中的全部细菌数。Ace指数和Chao指数均为物种丰富度的指标,实验组较对照组有所增加,表明添加海藻糖提升了高盐环境下系统内的物种丰富度。Shannon指数反映微生物多样性特征,其数值同样在实验组中有所增长,说明添加海藻糖丰富了高盐环境下系统内的菌群多样性。C120组的Chao指数、Ace指数以及Shannon指数均最高,说明120 μmol·L−1的海藻糖浓度对增强高盐胁迫下微生物群落的抵抗力稳定性最有利。
2) 微生物群落组成分析。不同海藻糖浓度下微生物群落组成相似但相对丰度略有差异(图4)。如图4(a)所示,门水平相对丰度排名前3的优势菌门分别为变形菌(Proteobacteria)、拟杆菌(Bacteroidota)和厚壁菌(Firmicutes)。三者的总相对丰度在对照组C0中为82.60%,在各实验组 (C40、C120、C360和C1 080) 中分别提高了3.29%、5.64%、5.49%和3.84%,说明在高盐条件下添加海藻糖对系统中门水平的群落结构具有一定影响。有报道[36-37]指出,Proteobacteria和Bacteroidota在高盐、高氨氮废水处理中常占据主导地位,其包含着多样性最丰富的硝化和反硝化细菌 (如HN-AD菌) ;也有报道[38]指出,可在极端环境进行异养硝化的脱氮微生物 (如肠球菌属Enterococcus) 归属于Firmicutes。该3类菌门中微生物总相对丰度的增加说明添加海藻糖有利于高盐胁迫下脱氮菌群的生长富集。
由属水平相对丰度排名前50的微生物进行层级聚类Heatmap分析结果 (图4(c)) 可知,在C0组中,不动杆菌属 (Acinetobacter) 相对丰度为34.58%,是系统内的优势菌属,Acinetobacter在高氨氮浓度下是硝化和反硝化过程中生长的主要菌群,表现出高效的HN-AD能力[39],且具有耐受高盐、高氨氮和脱氮速率快的优点[40]。而添加海藻糖的4个实验组中,高丰度优势菌属仍为Acinetobacter且相对丰度保持稳定(分别为35.69%、34.95%、34.93%和31.04%)。这说明在高盐环境下,海藻糖对Acinetobacter生长富集的影响较弱,但Acinetobacter保证了高盐系统内的高脱氮效率。在C120组中,假黄褐藻属 (Pseudofulvimonas) 相对丰度明显高于对照组和其他实验组,并且其相丰度变化存在显著性差异 (P <0.05) (图4(b)) 。Pseudofulvimonas相对丰度在低浓度海藻糖时 (40、120 μmol·L−1),随海藻糖浓度的增加而增长,但在高浓度海藻糖时 (360、1 080 μmol·L−1) ,却随海藻糖浓度的增加而减少。Pseudofulvimonas为耐盐型HN-AD菌,是高盐废水处理过程中常见的脱氮功能微生物[41-42];同时,黄杆菌属 (Flavobacterium) [43]、红杆菌属 (Rhodobacteraceae) [44]、假单胞菌属 (Pseudomonas) [16]、副球菌属 (Paracoccus) [45]和产碱菌科 (Parapusillimonas) [46]等耐盐型HN-AD菌在海藻糖浓度作用下的相对丰度变化规律与Pseudofulvimonas一致并同样存在显著性差异 (P<0.05) (见图4(b)) 。Pseudofulvimonas通常与各种氮物质降解过程有关[47];Flavobacterium与生物絮体聚集和颗粒稳定性有关,其丰度增加有利于功能微生物的富集进而保证污染物稳定有效的去除[42, 48];Rhodobacteraceae是海水形成初始生物膜的主要定殖者[44],其具有较强的耐盐性能,能将NH4+-N转化为NO3−-N和NO2−-N,且可在曝气的条件下实现对NO3−-N的高效去除[49-50];Pseudomonas、Paracoccus具有亚硝酸盐还原和有机物降解等多种功能,在有机物去除方面起着关键作用[51-52]。当海藻糖浓度为120 μmol·L−1时,这些HN-AD菌属相对丰度均达到最高。说明低浓度海藻糖促进高盐废水中耐盐型HN-AD菌生长富集,使系统内主要发生异养硝化和好氧反硝化过程,可强化微生物对氮、碳污染物的去除性能。
在C360和C1 080组中,Enterococcus相对丰度由C0中的0.12%分别增至13.88%和21.22% (P <0.05) (图4(b)) 。Enterococcus是一种兼性脱氮菌[53],具有较强的脱氮耐盐性能及高效的产絮凝效果,在生长代谢过程中产生多糖、蛋白质、脂和核酸等成分,可促进水体中菌体细胞凝聚进而减缓功能微生物的流失[54-55]。Enterococcus相对丰度在C360和C1 080组得到明显提升,说明高浓度海藻糖的作用倾向于通过增强高盐条件下微生物凝聚能力,提高功能菌群富集。此外,硝化细菌属 (Nitrincola) 、特吕珀菌属 (Truepera) 和梭菌属 (Fusibacter) 3种HDB的相对丰度在高浓度海藻糖时 (360、1 080 μmol·L−1) 也呈现增长趋势。Nitrincola是一类耐盐碱的厌氧反硝化细菌[15],能够在极端环境中对地球化学元素循环起促进作用[56];Truepera[32]与厌氧反硝化作用密切相关[57],且能够在碱性、中等盐度和高温等多种极端条件下生长,有利于NO向N2O的转换[58-59];Fusibacter是一类耐盐、厌氧的棒杆细菌,可将硫酸盐和元素硫还原成硫化物,还能够广泛利用碳水化合物 (如乙酸、葡萄糖等) 生成CO2和H2等产物[60],且具有脱氮除磷的功能[61]。上述3种菌属均在高浓度氨氮环境的反硝化过程中发挥重要作用,其相对丰度在高海藻糖浓度作用下提高,也证实高浓度海藻糖促进高盐废水中HDB生长富集,使系统内主要发生常规异养反硝化过程。这主要是因为高浓度海藻糖促进Enterococcus相对丰度大幅增加,虽有利于微生物凝聚,但同时也产生大量菌胶团聚集堵塞膜孔,致使生物膜层出现缺氧环境,不利于好氧脱氮微生物的生长,继而引起HDB增殖。虽然HDB相对丰度在高海藻糖浓度下的增长促使C360和C1 080组的NH4+-N和TN去除率相比C0组有所提高,但其相比C40和C120组却呈下降趋势,推测耐盐型HN-AD菌的相对丰度变化对系统脱氮性能的影响相比HDB更为显著和重要,而HDB的相对丰度变化主要影响了系统的脱氮菌群结构。
3) 环境因子相关性分析。通过相关性Heatmap图进一步分析环境因子对系统内属水平群落组成的影响,结果如图5所示。与NH4+-N、COD及TN去除率和海藻糖浓度4个因子均显著关联的只有Pseudomonas的相对丰度,且呈负相关;NH4+-N、COD及TN去除率与Flavobacterium、Parapusillimonas 2种耐盐型HN-AD菌属相对丰度呈正相关,与浮游盐属 (Planktosalinus) 、新巨藻属(Neomegalonema)相对丰度呈负相关。Planktosalinus归属黄杆菌科 (Flavobacteriaceae) ,而Flavobacteriaceae是废水系统发生污泥膨胀的指示微生物[62-63];Neomegalonema同样属于污泥膨胀菌群,其丰度变化用于表征氮磷失衡下废水中的微生物群落健康状况[64-65]。海藻糖浓度与Enterococcus、球毛菌属 (Sphaerochaeta) 和脱硫微杆菌属 (Desulfomicrobium) 相对丰度呈正相关。可以看出,不同菌属相对丰度受到海藻糖浓度、污染物去除率影响的情况存在差异,可证明耐盐型HN-AD菌确为系统内主要脱氮贡献者,显著影响反应器的脱氮性能;以Enterococcus为代表的传统脱氮菌和其他功能微生物的相对丰度变化规律则主要反映了海藻糖浓度梯度对高盐废水中脱氮菌群组成的差异影响。
-
采用PICRUSt1对功能基因丰度进行分析,其相对丰度变化结果见图6(a)。该研究获得与硝化过程相关的羟胺氧化酶基因hao,是NH2OH氧化生成NO2−-N过程的关键酶基因[66],其在实验组 (C40、C120、C360和C1 080) 中相对丰度相较对照组 (C0) 分别提高了4.96倍、34.70倍、32.45倍和29.06倍,说明添加海藻糖明显加快系统内微生物的硝化进程。同时,还获得8个与反硝化过程相关的功能基因,包括硝酸盐还原酶(nasA、nasB、narG、narH、narI、napA、napB)、亚硝酸盐还原酶(nirK)、一氧化氮氧化酶(norB、norC)和一氧化二氮还原酶(nosZ)。4个实验组中反硝化功能基因总相对丰度较对照组分别提高了1.31倍、2.86倍、1.81倍和1.75倍,说明添加海藻糖同样促进了系统内微生物的反硝化进程。
硝酸还原酶催化细菌体内NO3−-N还原生成NO2−-N,是实现反硝化过程的关键酶之一。nasA是一种同化型硝酸盐还原酶基因,存在于细菌的细胞质中[67],也常见于HN-AD菌[68]。其在实验组中的相对丰度值相较对照组均有提高,低浓度海藻糖时 (40、120 μmol·L−1) ,nasA基因相对丰度随海藻糖浓度的增加而增长,于海藻糖浓度为120 μmol·L−1时达到最高值;但在高浓度海藻糖时 (360、1 080 μmol·L−1) ,其随海藻糖浓度增加出现降低。同样,异化型硝酸盐还原酶NAP的基因napA和napB相对丰度变化规律与nasA一致。NAP表达在好氧条件下占主导地位,常存于HN-AD菌中[69]。另一种异化型硝酸还原酶NAR的基因narG、narH和narI相对丰度却在海藻糖浓度为360 μmol·L−1和1 080 μmol·L−1时增加,而NAR表达在缺氧条件下占主导地位,常见于HDB[70]。以上结果说明,在高盐废水添加海藻糖可诱导增强以HN-AD菌和HDB为主的脱氮菌群的反硝化活性。低浓度海藻糖倾向于强化系统的异养硝化-好氧反硝化途径,而高浓度海藻糖则偏向于促进系统的异养反硝化途径 (图6(b)) 。
亚硝酸盐还原酶NIR是反硝化过程的另一类关键酶,驱动NO2−-N还原生成气体产物NO,存在于细胞膜外生物细胞质中,其基因有nirS和nirK 2种类型[71]。海藻糖的加入提高了nirK的相对丰度,且在海藻糖浓度为120 μmol·L−1时达到最高值,推进了系统的亚硝酸盐还原进程。NO还原酶基因 (norB和norC) 和nosZ分别是催化NO还原生成N2O、N2O进一步还原生成最终产物N2的关键酶基因。在污水处理中,反硝化细菌体内nosZ活性的降低或不具有nosZ是导致反硝化过程中温室气体N2O大量积累的主要因素[72]。该研究nosZ的相对丰度随着海藻糖浓度的增加而升高,说明添加海藻糖对控制好氧脱氮过程的N2O产生与积累具有一定作用。以上基因预测结果也进一步证明,海藻糖浓度为120 μmol·L−1时,反应器的NH4+-N、TN去除率最高 (图2) 。这主要是因为此海藻糖浓度条件激发HN-AD菌的硝化和反硝化进程达到最快,极大程度加速了NH4+-N氧化及NO3−-N和NO2−-N的还原,最终使反应器脱氮性能达到最优。
-
1) 添加海藻糖的实验组 (40、120、360和1 080 μmol·L−1) 中NH4+-N、TN和COD去除率相较对照组 (0 μmol·L−1) 显著提高,最适海藻糖浓度为120 μmol·L−1。外加海藻糖可优化高盐条件下的反应器脱氮性能。
2) 海藻糖浓度对脱氮菌群的相对丰度具有显著影响。低浓度海藻糖 (40、120 μmol·L−1) 有利于耐盐型HN-AD菌相对丰度的提高,而高浓度海藻糖 (360、1 080 μmol·L−1) 更倾向于促进HDB的富集;耐盐型HN-AD菌是系统的主要脱氮贡献者,其相对丰度变化与污染物去除率之间存在显著相关性,传统脱氮菌相对丰度变化则与高盐废水中脱氮菌群组成差异有关。耐盐型HN-AD菌对调控高盐条件下的反应器脱氮性能起着重要影响。
3) 海藻糖可有效提高以HN-AD菌和HDB为主的脱氮菌群的硝化和反硝化活性。在低浓度海藻糖 (40、120 μmol·L−1) 作用下,与HN-AD菌关联的硝化、反硝化基因(hao、nasA、napA和napB)丰度随海藻糖浓度的增加而增长,系统的异养硝化-好氧反硝化途径得到强化;在高浓度海藻糖 (360、1 080 μmol·L−1) 作用下,与HDB关联的反硝化基因 (narG、narH和narI) 相对丰度明显增加,系统的异养反硝化途径得到增强。
