下载: 
点击查看大图全文HTML
--> --> --> 科技的高速发展在提高人类生活水平的同时对环境也会造成极大影响，近年来，国内外水体污染问题日益严峻[1-3]。2015年，中国废水排放总量为7.353×1010 t。其中，城镇生活污水排放量达到了5.352×1010 t[4]。城镇生活污水的处理问题备受关注，因此，AO、A2O、OAO、氧化沟等优良工艺应运而生[5-6]。缺氧好氧工艺即AO工艺是我国目前广泛应用的一种生物脱氮工艺，由于其脱氮效果较好、原理简单被国内外****和工程师持续研究和改进[7-9]。OAO工艺即是AO工艺的一种改良工艺，其在传统AO工艺的缺氧池前端增设一级曝气池，具有良好的废水处理效果[10-11]。
微生物是污水生物处理的核心部分，利用分子生物学技术对污水处理过程中活性污泥的微生物菌群进行动态跟踪和功能群种鉴定，有助于筛选并培养出利于污水处理的优势菌种，构建更加合理的活性污泥生态系统，提高污水处理效率[12]。MATSUDA等[13]对日本大阪8个城市污水厂的活性污泥微生物组成进行了分析，发现不同工艺在不同时期微生物群落差异很大。DING等[14]和YE等[15]分别应用AO工艺和A2O工艺处理混合废水和酱油废水，发现Proteobacteria是污水处理中的第一优势菌门。HAN等[16]、曾涛涛等[17]研究氧化沟工艺处理生活污水时发现，大部分细菌属于Proteobacteria和Bacteroidetes这2大门类，其优势菌属有Pseudomonas、Phycisphaera、Methylocystis等。但是，关于OAO工艺处理城镇生活污水过程中活性污泥的微生物特性，目前还少有研究。本研究拟采用16S rDNA高通量测序技术对OAO和AO工艺处理城镇生活污水过程中各反应池的活性污泥中的微生物进行测序分析，为OAO工艺处理城镇生活污水中优势菌种的筛选及培养提供参考。
1.1. 实验装置
实验装置如图1所示，AO和OAO工艺所用一体化装置均为有机玻璃制作。AO工艺装置长、宽、高分别为0.78、0.4、0.6 m，其中，缺氧池有效容积为30 L，好氧池有效容积为90 L。OAO工艺则在AO工艺基础上，前置1个容积为46 L的曝气池及1个一级沉淀池，装置尺寸为1.16 m×0.4 m×0.6 m。
1.2. 实验用水
本实验采用实验室配置的模拟城镇生活污水，水质指标如下：COD为340~360 mg·L?1、NH3-N为50.2~72.3 mg·L?1、TN为54.2~75.5 mg·L?1、TP为3.3~8.0 mg·L?1、pH为7.3~8.0。1.3. 污泥的驯化与反应器的运行
本研究从郑州市某生活污水处理厂的污泥脱水机房获取污泥，其含水率约为80%，将反应器平均污泥浓度设为3 500 mg·L?1，以此计算出OAO工艺需投加的污泥量为2 905 g，AO工艺为2 100 g。同时，从该污水厂硝化液回流管道中取部分混合液，将混合液与污泥一并加入反应器内。而后，开始进水并搅拌混合液，使泥水混合均匀。同时，开启风机进行闷曝，闷曝2 d后，污泥基本恢复生物活性。闷曝工作完成后，开始连续不间断地向反应器内进水以驯化污泥，AO工艺进水流量约为12 L·h?1，OAO工艺约为16.6 L·h?1。AO工艺污泥回流比为100%，混合液回流比为200%，OAO工艺2段污泥和混合液回流比均为100%。AO工艺为单点进水，新鲜污水均进入缺氧池中。考虑到缺氧池碳源问题，OAO工艺采用两点进水，1/3的污水进入缺氧池，2/3进入曝气池。驯化期间，每天对进出水水质进行检测，COD采用重铬酸钾法(HJ 828-2017)、NH3-N采用纳氏试剂分光光度法(HJ 535-2009)、TN采用碱性过硫酸钾消解紫外分光光度法(HJ 636-2012)、TP采用钼酸铵分光光度法(GB 11893-1989)、溶解氧采用电化学探头法(HJ 506-2009)、pH采用玻璃电极法(GB 6920-1986)、温度使用水温计法(GB 13195-1991)进行测量[18-24]。进水检测大约持续30 d，AO工艺在第9天达到稳定状态，OAO工艺在第12天时达到稳定，稳定后，2种工艺对COD、NH3-N、TN、TP的平均去除率见图2。2种工艺稳定后的出水COD、NH3-N指标基本符合《城镇污水处理厂污染物排放标准》中的一级A标准[25]。
1.4. 样品的采集与DNA获取
在连续进水第30天时进行取样，此时2种工艺各反应池均已稳定运行。污泥样品取自AO工艺的好氧池、AO工艺的缺氧池、OAO工艺的好氧池、OAO工艺的缺氧池、OAO工艺的曝气池，分别记为AO.H、AO.Q、OAO.H、OAO.Q、OAO.B。每个反应池取100 mL泥水混合液，静置30 min，除去上清液，所得沉淀混合均匀后，分装于瓶中，保存于?30 ℃环境中，以便后续DNA提取。污泥样品中的DNA采用PowerSoil?DNA提取试剂盒，按照试剂盒使用方法提取。1.5. 16S rDNA高通量测序
采用NanoDrop 1000分光光度法测定样品中提取的DNA的浓度和纯度[26]，将符合要求的DNA样品使用无菌水稀释至1 ng·μL?1后，用作模板DNA。以上述基因组DNA为模板，在16S rDNA基因的保守区V4区域，通过融合引物515F(5′-GTGCCAGCMGCCGCGG-3′)和806R(5′-GGACTACHVGGGTWTCTAAT-3′)进行PCR扩增，扩增反应体系体积为30 μL，扩增流程为：95 ℃预变性3 min；95 ℃变性30 s，55 ℃退温30 s，72 ℃延伸40 s，进行30个循环；72 ℃延伸10 min；最后将样品保持在10 ℃恒温环境中。扩增产物采用1%琼脂糖凝胶电泳检测，选择明亮主条带的样品用于进一步的实验。使用AxyPrep DNA凝胶抽提试剂盒(AXYGEN)对PCR产物切胶纯化。之后，通过Illumina建库试剂盒和相应的质检方案完成文库的构建和质检。将合格的文库送往Illumina平台(Miseq)进行测序，用下机得到的数据进行相应的生物信息分析。2.1. 微生物群落多样性分析
为了反映样品中微生物菌群的丰富度及多样性，对每个样本分别进行微生物菌群的Alpha多样性分析[27]，结果见表1。覆盖率均在99%以上，证明本次测序数据可靠度极高。在97%相似度下，将样本中的优质Tags进行聚类，共得到1 237个OTUs，其中OAO工艺的OTUs为1 135个，AO工艺为1 063个。OAO工艺中的3个反应器都比AO工艺中的OTU数目高，反映了本实验中OAO工艺的细菌数目要大于AO工艺。从Chao指数和Ace指数来看，2种工艺各反应器的物种丰富度为OAO曝气池<AO好氧池<AO缺氧池<OAO缺氧池<OAO好氧池。OAO工艺独有的曝气池丰富度较低，初步分析其可能的原因是，该曝气池的相对进水流量为AO工艺和OAO工艺后续反应池的2/3，污水水量较少，营养物质浓度较低，且没有混合液回流，水源单一，进水水质的差异造成了其丰富度的差异[28]。同时，曝气池后的初沉池污泥直接内回流于曝气池，处于一种相对封闭的状态，故OAO工艺的一级曝气池丰度要略低于其他反应器。从Shannon和Simpson指数来看，OAO工艺中的微生物群落多样性显著高于AO工艺。2种工艺好氧池物种多样性均大于缺氧池，其原因可能是，在缺氧和好氧2种不同运行条件下，微生物的生态系统体系有所差异[29]。总体来看，OAO工艺的微生物物种丰度和多样性均优于AO工艺。2.2. 门、纲水平的微生物物种分布及分析
本次测序共检测出已知门38个，结果见图3。OAO工艺中的细菌门类较AO多。Proteobacteria和Bacteroidetes为2种工艺的第1和第2优势菌门，这与HAN等[16]、曾涛涛等[17]、郑向阳等[30]在氧化沟和AO工艺研究中的研究结果一致。Proteobacteria和Bacteroidetes是自然界细菌中2个较大的门类，一些反硝化细菌和大多数固氮细菌也属于这2种门类[31]。康晓荣[32]研究发现总氮和总磷的去除率与Proteobacteria和Bacteroidetes的丰度成正相关。本实验中，AO工艺中的Proteobacteria丰度要大于OAO工艺，但OAO工艺中的Bacteroidetes丰度要远大于AO工艺。
在纲水平下，对Proteobacteria门和Bacteroidetes门的细菌在各反应器中的分布进行统计分析，结果如表2所示。γ-Proteobacteria除在OAO的曝气池为第2优势菌纲外，在其他各反应器均为第1优势菌纲。在Proteobacteria门里，除γ-Proteobacteria外，β-Proteobacteria、α-Proteobacteria、δ-Proteobacteria也为2种工艺的优势菌纲。β-Proteobacteria和γ-Proteobacteria的富集可以大大提高污水中有机物和氮、磷去除率[33]；α-Proteobacteria作为一种成员变异性极大、共通点极少的菌纲[34]，在OAO的各反应器内的占比均明显高于AO工艺，其成员有根瘤菌、立克次氏体、光合细菌等，在某些光生物反应器实验中表现为优势菌种[33, 35]。在Bacteroidetes门中，Saprospirae、Sphingobacteria、Flavobacteriia、Cytophagia在2种工艺里占比较高。Cytophagia在AO工艺中的占比几乎是OAO工艺的2倍，Flavobacteriia在OAO缺氧池占比较小，其他该门菌纲在OAO工艺中的各反应器占比明显高于AO工艺。OGONOWSKI等[36]研究发现，Cytophagia在天然纤维素基材上富集较多，该菌纲在OAO工艺较少且曝气池>缺氧池>好氧池，这可能是因为部分纤维素作为部分微生物的碳源在前置曝气池中被转化。Flavobacteriia在许多海洋环境中都很丰富，可以在极端环境中降解水中不稳定的有机碳[37]。Saprospirae、Sphingobacteria在一些其他AO改良工艺和厌氧、好氧工艺中也表现为优势菌纲[38-40]。
除上述2种菌门外，在2种工艺中均存在且丰度较高的菌门还有Nitrospirae、Planctomycetes、Firmicutes、Verrucomicrobia、OD1等，Nitrospirae、Planctomycetes、Firmicutes中的某些细菌具有较好的碳氮去除能力，短程硝化反硝化中的厌氧氨氧化菌就来自Planctomycetes。TM7、Actinobacteria、Chlamydiae、Acidobacteria等菌门在OAO工艺中的丰度较高，但在AO工艺中则相反。Actinobacteria是OAO好氧池中的第3优势菌门，具有多样化的次级代谢，在生态环境中有许多潜在用途[41]；Acidobacteria是OAO缺氧池中的第4优势菌门，在土壤生态中占细菌总量的30%~50%，具有纤维素降解的功能，并且参与土壤中的铁循环和光合作用[42]。除Verrucomicrobia外，AO工艺中相对丰度大于0.5%的菌门在OAO中的相对丰度也均在0.5%以上。Verrucomicrobia作为一种在氮、磷、钾含量较低的情况下依然能利用碳源生长的菌门[43-44]，在AO工艺中的平均占比几乎是OAO工艺中的2倍，这与OAO工艺前置曝气池消耗掉部分有机物不无关系。
2.3. 属水平的微生物物种分布及分析
在本次测序的5个样本中，共检出菌属347个，丰度大于0.5%的已知菌属为23个，各池微生物群落相对丰度见图4。未知菌属(Unclassified)在各池的占比为35.07%~65.13%，在OAO工艺内的相对丰度要远大于AO工艺。已知属中来自Proteobacteria门的Thiothrix、Thermomonas、Dechloromonas相对丰度较高，三者总丰度占比为24%。与何嘉鹏[45]的实验结果不同，作为硝化细菌的重要菌属Nitrosomonas、Nitrospira虽均在2种工艺中检出，但含量皆小于0.5%，在本实验中，被归类为其他菌属(others)。Thiothrix、Pseudomonas、Hydrogenophaga是AO工艺好氧池中的优势菌属；Thiothrix、Dechloromonas、Leadbetterella是AO工艺缺氧池中的优势均属；Thermomonas、Flavobacterium、Dechloromonas是OAO工艺曝气池中的优势菌属；Thermomonas、Dechloromonas是OAO工艺好氧池中的优势菌属；Thermomonas、Dechloromonas、Rhodobacter是OAO工艺缺氧池中的优势菌属。
WANG等[46]研究发现，Thermomonas与硝酸盐还原酶基因呈正相关，对同时厌氧氨氧化-反硝化系统中的硝酸盐还原起重要作用；WANG等[47]研究发现，Dechloromonas在反硝化过程时可以进行磷的去除。在本实验中，上述2种菌属在OAO工艺各部分的相对丰度均大于AO工艺，这可能是OAO工艺对TN和TP去除效果显著优于AO工艺的原因之一。OAO工艺里其他优势菌属中，Flavobacterium为好氧反硝化菌[48]，Rhodobacter具有固氮的功能，也属于好氧反硝化菌[49]。AO工艺里，Pseudomonas为第1优势菌属，其主要功能为异养硝化-好氧反硝化菌[50]和反硝化聚磷菌[48]，其在AO工艺和OAO工艺好氧池中的大量存在可能是2种工艺高效脱氮的重要原因之一。Thiothrix作为一种能以硫化物为唯一能源生长的菌属，被认为是造成活性污泥膨胀的主要原因之一[51]。本次实验虽未见明显的污泥膨胀现象，但2种工艺在HRT、温度、水质等条件相同的情况下，Thiothrix在AO工艺中的相对丰度要远大于OAO工艺，这使得AO工艺发生污泥膨胀的可能性要远大于OAO工艺。Nitrosomonas在2种工艺的好氧池中的相对丰度分别为AO.H 0.05%、OAO.H 0.47%；Nitrospira则为AO.H 0.02%、OAO.H 0.45%。MA等[52]研究发现，污水厂活性污泥中的AOB和NOB相对丰度在0.1%~1%时，脱氮效率依旧很高，这与本次实验水质检测结果基本相符。2种硝化菌属在2种工艺中好氧池的丰度均为OAO好氧池>AO好氧池，这可能是因为OAO工艺中的第1段曝气池的作用使得工艺整体的污泥BOD负荷较AO工艺小，从而减轻硝化细菌与其他异养型细菌的竞争压力，为硝化细菌提供了更佳的生长环境。
2.4. 2种工艺微生物群落的相似性与差异性
图5为2种工艺各反应池中微生物在属水平下的物种丰度热图，该图将微生物属水平的聚类结果按照其在该样品中所占相对丰度数值用不同梯度的颜色表示，用颜色矩阵来直观地表示不同样品中物种的丰度。微生物群落组成在不同反应池间的相似性和差异性由纵向聚类表示，距离越近，枝长越短，说明样品的物种组成及丰度越相似。由于本次测序结果显示不同物种间相对丰度差异较大，颜色梯度难以满足要求，因此，将原始相对丰度以10为底进行对数转化，若某物种相对丰度为0，则取样品中物种相对丰度最小(大于0)的值除以2然后进行对数转化。
由图5可知，OAO和AO工艺间微生物群落具有较大的差异性，OAO工艺内部的缺氧池、曝气池和好氧池3个反应池间的差异性要略大于AO工艺内部的2个反应池间的差异性。这进一步说明OAO工艺中活性污泥微生物生态体系更为复杂，多样性更高。Planctomyces、Dokdonella、Mycoplana等在2种工艺中的相对丰度相近，初步推测其受2种工艺间差异的影响较小，采用OAO工艺并不能使该类菌属在原AO工艺中发挥的作用发生改变。Thiothrix、Pseudomonas、Tolumonas、Zoogloea、Leadbetterella、Hydrogenophaga、Janthinobacterium在AO工艺的丰度明显大于OAO工艺，同时，Thermomonas、Runella、Flavobacterium、Rhodobacter在OAO工艺中的丰度又显著大于AO工艺。采用OAO工艺会增加或减少上述菌属在原AO工艺反应器内的丰度，某些具有脱氮除磷功能的菌属的增减使污水处理效果发生变化。
5个反应池内的优势菌属中部分功能菌属的相对丰度统计分析结果见表3。Hydrogenophaga、Pseudomonas作为异养硝化-好氧反硝化菌[53]，在BOD负荷较高有机碳源充足的AO工艺好氧池中大量存在，其可在同一个空间内进行硝化和反硝化，在本次AO工艺生物脱氮过程中起到了重要作用。Dechloromonas、Flavobacterium、Rhodobacter、Thermomonas均为具有反硝化功能的菌属，在OAO工艺中的相对丰度大于AO工艺，尤其是Thermomonas作为一种自养反硝化菌[54]，在BOD负荷较小的OAO工艺中的相对丰度远高于AO工艺。这解释了OAO工艺在对TN的去除效果上优于AO工艺的原因。
2)微生物群落的Alpha多样性分析表明，AO工艺和OAO工艺二者间的丰度和多样性有显著差异，OAO工艺的丰度和多样性皆优于AO工艺。
3)微生物物种结构分析结果表明，2种工艺的优势菌门均为Proteobacteria和Bacteroidetes，且Proteobacteria门中的β-Proteobacteria和γ-Proteobacteria均为2种工艺的优势菌纲。2种工艺在优势菌属上差异较大，OAO工艺的优势菌属为Thermomonas、Dechloromonas、Rhodobacter，AO工艺的优势菌属为Pseudomonas、Thiothrix、Dechloromonas。
4)属水平下物种丰度热图表明，OAO和AO工艺间微生物群落具有较大的差异性，Thiothrix、Pseudomonas、Thermomonas等受2工艺间差异的影响，相对丰度差异性较大。
5)由于OAO工艺前置曝气池的设置，OAO工艺污泥BOD负荷较AO工艺小，导致Verrucomicrobia门、Cytophagia纲等对碳源需求较高的微生物在OAO中丰度较小，Thermomonas属、Dechloromonas属、Nitrospira属、Nitrosomonas属等具有脱氮除磷功能的细菌有了合适的生长环境，其在OAO的丰度明显高于AO工艺。据此推测，若将OAO工艺用于处理中高浓度的污水，其脱氮效率要远大于AO工艺。
参考文献
