近年来，分子生物学作为有效手段用来研究污水处理过程中微生物群落结构特性，如变形梯度凝胶电泳、克隆文库和高通量测序等. MiSeq高通量测序以Illumina的测序技术为基础，通过可逆终止试剂方法对数百万个基因片段同时进行大规模平行测序，具有分析结果精确、高速等特点，被广泛(extensive)应用于污水处理过程中微生物结构和多样性研究.
本文以颗粒硫磺和黄铁矿的硫自养反硝化反应器为研究对象，探究反应器在经过30 d饥饿忍耐后的恢复情况，并利用MiSeq高通量测序对饥饿忍耐前后反应器中细菌群落的变化情况进行分析，以期为污水处理厂脱氮保障工艺的季节性、间歇性运行提供技术参考.
1 材料与方法1.1 试验装置
硫自养反硝化工艺采用降流式生物滤池，试验装置示意图如图 1所示.反应器为密封的有机玻璃材质柱，内径140 mm，高度1 170 mm，填充高度500 mm，有效容积5.94 L.进水口距柱底600 mm，沿反应器不同高度处分别设置取填料口，取填料口距填料顶部的距离分别是100 mm和300 mm.
图 1 试验装置示意
试验共设2个反应器，1号反应器加入硫磺和白云石的混合物，2号反应器加入黄铁矿和白云石的混合物. 2个反应器中的填料均按
1:1的体积比混合装填，硫磺、黄铁矿和白云石的粒径均为5~10 mm.其中，硫磺/白云石的填充区域孔隙率为45.9%，黄铁矿/白云石为44.1%.
1.2 进水水质和接种污泥
本试验用水为人工模拟污水厂尾水，水质指标为: NO3--N浓度为30 mg?L-1，NH4+-N浓度为2 mg?L-1，TP浓度为1 mg?L-1，COD浓度为18~23 mg?L-1，pH为6.8~7.2.接种污泥取自高碑店污水处理厂二沉池回流污泥.
1.3 反应器运行方式
反应器在低温下运行85 d ，主要分为稳定（解释:稳固安定；没有变动)期、饥饿期和恢复期:稳定期，正常进水，反应器稳定运行，平均温度为12℃;饥饿期，停止进水，反应器处于饥饿状态，平均温度12.8℃;恢复期，恢复进水，反应器进入恢复期，平均温度为14℃.饥饿期结束后及恢复期采集反应器中的污泥样品 进行MiSeq高通量测序分析，分析反应器饥饿期及稳定期细菌群落的变化情况.
1.4 测试指标和方法1.4.1 水质指标的测定
反应器运行期间，定期采样进行水质指标的测定，检测项目包括氨氮水(Nitric acid)盐氮、亚硝氮、总氮、总磷(P) 等.其中NO3--N采用紫外分光光度法测定;NO2--N采用N--乙二胺分光光度法测定;TN采用碱性过硫酸钾消解紫外分光光度法测定;TP采用钼酸盐分光光度测定.试验所用分光光度计为HACH DR5000紫外可见分光光度计.生物量采用脂磷法测定.
1.4.2 微生物多样性的测定
DNA的提取与PCR的扩增:为了探究系统的特殊结构：荚膜、鞭毛、菌毛群落，在反应的饥饿期和恢复期进行取样.采用E. Z. N. A. Soil DNA试剂盒，按照试剂盒说明书提供的操作控制步骤提取.提取的DNA用1%琼脂凝胶电泳进行检测. PCR的扩增区域为16S rRNA的V3-V4区，细菌16S rRNA扩增引物采用通用引物.引物名称和引物序列分别是338F 和806R . PCR反应体系: Forward primer ，2 μL;Reverse primer ，2 μL;dNTPs ，4 μL;10×PCR buffer，5 μL;Pyrobest DNA Polymerase ，0.3 μL;补充ddH2O至50 μL.反应程序:先95℃预热5 min;然后进行25个循环，最后72℃延伸10 min.用AxyPrepDNA凝胶回收试剂盒 回收PCR产物，经Tris-HCl洗脱后，用2%琼脂糖凝胶电泳检测.根据电泳结果.将PCR产物用QuantiFluorTM -ST 进行定量，按照每个样品的测序量要求，进行相应比例的混合.
MiSeq文库构建与测序:用高通量测序平台Illumina MiSeq对扩增产物进行MiSeq高通量测序.测序得到的PE reads根据overlap关系进行拼接，同时过滤掉低质量的reads，区别样品后进行OUT聚类分析和物种分类学分析.
生物多样性和分类学分析(Analyse):首先将序列按照彼此的相似性归为操作分类单元.按照97%相似性进行OUT聚类，采用RDP classifier对97%相似水平的OUT代表序列进行分类学分析.利用MOTHUR软件对样品覆盖率，ACE，Chao丰富指数以及Shannon多样性指数进行计算.同时为了保证分析的高可信度，根据SILVA106库中的参考序列对OUT进行种属鉴定，种属比对的可信度阈值设定为80%.
2 结果与讨论2.1 反应器运行情况
稳定期、饥饿期和恢复期的硫自养1号和2号反应器进出水水质的变化情况如图 2所示.进水NO3--
N、NO2--
N、TN和TP平均浓度分别为30.00、0.01、35.88和1.00 mg?L-1.稳定期，反应器在低温下 运行30 d后，1号反应器出水NO3--
N、NO2--
N、TN和TP浓度分别为1.78、0.03、2.87和0.91 mg?L-1，TN去除率为91.9%，1号反应器具有较高的脱氮效果. 2号反应器出水NO3--
N、NO2--
N、TN和TP浓度分别为11.32、0.11、13.10和0.14 mg?L-1，TN和TP去除率分别为63.49%和86.41%. 2号反应器具有同步脱氮除磷(P)的效果.
图 2 饥饿前后反应器NO3--
N、NO2--
N、T
N、TP的变化
反应器稳定运行后，停止进水，进入饥饿期.经过30 d的饥饿忍耐后，反应器在低温下 恢复进水，进入恢复期.由图 2可以看出，饥饿忍耐后，1号反应器出水NO3--N增加到27.87 mg?L-1，2号反应器出水NO3--N增加到26.56 mg?L-1.但随着反应器的恢复，出水NO3--N浓度逐渐降低，1号反应器经5 d的恢复，出水NO3--N和TN浓度分别为0.38 mg?L-1和2.37 mg?L-1，去除率分别为98.8%和93.6%;2号反应器经11 d的恢复，出水NO3--N和TN浓度分别为10.51 mg?L-1和12.91 mg?L-1，去除率为66.40%和65.57%.其次，在恢复期2个反应器的出水NO2--N浓度均呈现先增加后降低的趋势，恢复初期均出现NO2--N的积累.由图 2可知，1号反应器恢复期第3 d 出水NO2--N浓度为2.17 mg?L-1;2号反应器恢复期第5 d ，出水NO2--N浓度为0.37 mg?L-1.随着反应器稳定运行，最终1号反应器NO2--N浓度下降(descend)到0.1 mg?L-1以下，2号反应器有0.1 mg?L-1 NO2--N的积累量.此外，饥饿忍耐未对2号反应器TP的去除效果产生明显的影响.原因如下，2号反应器对磷(P)的去除主要是通过黄铁矿中的铁与磷结合形成铁磷沉淀物及铁的氧化（oxidation)物和氢氧化物对磷有吸附作用，饥饿条件下并未明显影响这一物理化学过程.结果表明，与1号反应器相比，2号反应器的恢复时间略长，但2个反应器都能在低温、短期内恢复到正常处理(chǔ lǐ)效果.
饥饿期 及恢复期 反应器不同高度处微生物量的变化如表 1所示.从中可知，1号和2号反应器饥饿忍耐后不同高度处微生物量均低于恢复期.分析原因，主要是饥饿期间反应器中部分微生物对饥饿环境的忍耐能力差，裂解死亡;部分细菌因营养物质缺乏而进入休眠状态，不能进行生长繁殖.反应器恢复进水后，生物膜的活性和生长可得到不同程度的恢复，使得微生物量有所增加，同时反硝化性能得以恢复.
表 1 饥饿忍耐前后反应器不同高度处生物量的变化/nmol?g
2.2 微生物多样性的分析
利用MiSeq平台对A1、A2、B1、B2这4个污泥样品进行高通量测序，分别获得39 989、39 302、45 227、61 572条优化序列.将优化序列截齐后与SILVA106库比对后进行聚类，在97%的相似性下分别获得3 188、3 412、1 734、1 976个OUTs.并且4个污泥样品的覆盖率 均在95%以上，表明本研究中构建的序列库可以覆盖细菌群落的多样性.
表 2 样品的物种丰富度和多样性分析
Chao指数和Shannon指数表示微生物群落结构的变化.群落丰富度用Chao指数描述，其值越高表明群落物种的丰富度越高;Shannon指数可以反映样品的多样性程度，其值越高表明群落物种的多样性越高. 表 2显示，4个污泥样品的Chao和Shannon指数均为A2>A1>B2>B1.结果表明，1号反应器中的物种丰度和多样性均高于2号反应器.分析原因，主要是黄铁矿/白云石反应器 中所用硫源为黄铁矿，其硫含量低于硫磺.这使得硫磺/白云石反应器 中的硫自养反硝化菌群因得到充足基质而更好的生长.即由于基质含量的不同，黄铁矿/白云石反应器中的自养反硝化菌种的生长受到一定的限制.另一方面，经过30 d的饥饿忍耐后，硫磺/白云石和黄铁矿白云石2个反应器中微生物丰富度和多样性明显低于恢复期.分析原因，主要是饥饿期部分微生物对饥饿环境的忍耐能力差，裂解死亡，部分细菌因营养物质缺乏而进入休眠状态;恢复期生物膜的活性和生长可得到不同程度的恢复.这与上述反应器中微生物量饥饿期低于恢复期的变化是一致的.
2.3 群落结构分析
基于SILVA数据库的分类信息，对饥饿期及恢复期污泥样品的高通量测序数据进行了门、纲、属水平上的分类分析.门分类水平上，4个污泥样品共检测出39个类群，门水平上的大量类群 如图 3所示.从中可知，A1和A2样品中Proteobacteria为优势（解释:能压倒对方的有利形势)菌群，丰度分别为89.70%和89.77%，其次为Bacteroidetes 和Chlamydiae ，其他菌类比例相对较低. B1和B2样品中Proteobacteria为优势菌群，比例分别为51.82%和55.42%.此外，Bacteroidete
S、Chlamydia
E、Chloroflex
I、Actinobacteri
A、Acidobacteri
A、Chlorob
I、Planctomycete
S、Verrucomicrobi
A、Gemmatimonadetes也是B1和B2样品中主要的门类.结果表明，1号和2号反应器在生物群落结构和丰度上具有比较明显的差异.两组反应器最主要的优势菌门虽均为Proteobacteria，但其相对丰度存在明显的差异，且2号反应器在门水平上的主要类群较为多样.此外，每组反应器饥饿期和稳定期丰度变化不大.
图 3 门水平上群落结构
纲分类水平上，1号和2号反应器中共检测出68个类群，其中7个类群为主要的纲，如图 4所示. A1和A2样品中主要的纲类为β-GOOGLE PRoteobacteria，丰度为73.42%和69.15%，且β-Proteobacteria中的一些细菌是与污泥反硝化相关的[20~22];其次为α-Proteobacteri
A、γ-Proteobacteri
A、Sphinggobacteri
A、Chlamydiae . B1和B2样品中主要的纲类为α-Proteobacteri
A、β-Proteobacteri
A、γ-Proteobacteria.此外B1和B2污泥样品中还含有的δ-Proteobacteri
A、Ignavibacteri
A、Chlorobi
A、Actinobacteri
A、Anaerolinea
E、Acidobacteri
A、Gemmatimonadete
S、Phycisphaerae. 图 4表明，不同样品中类群结构及丰度差异性较大，1号反应器中β-Proteobacteria的相对丰度高于2号反应器，2号反应器中α-Proteobacteria和γ-Proteobacteria的相对丰度则高于1号反应器.
图 4 纲水平上群落结构
属分类水平上，1号和2号反应器大量类群如图 5所示.在属水平上，2个反应器中共检测出296个细菌（fungus)类群.其中A1和A2样品中主要属类为Sulfobacillu
S、Variovorax以及Thiobacillus，其丰度分别35.67%、8.82%、19.54% 和31.98%、13.25%、11.86%.其中，Sulfobacillus是噬酸且耐热的硫化物矿石氧化菌，也是亚铁离子氧化菌. Thiobacillus为革兰氏阴性细菌，是目前被报道最多的用于还原NO3--N的硫氧化细菌，可用于硫自养反硝化处理(chǔ lǐ)市政污水和地下水中的NO3--N.此外，A1和A2样品中的主要类群 还包括Thermomona
S、Ferruginibacte
R、Denitratisom
A、Sulfurimona
S、Rhizobium.
图 5 属水平上群落结构
与A1和A2样品相比，B1和B2样品中的主要类群包括Sulfobacillu
S、Variovora
X、Thiobacillu
S、Denitratisom
A、Ferruginibacte
R、Neochlamydi
A、Woodshole
A、Blastocatell
A、Bradyrhizobiu
M、Sulfuritalea等菌属.并且不同时期反应器中某些菌属的相对丰度有一定的差别，如B1和B2样品中Sulfobacillus和Thiobacillus在饥饿期和恢复期的相对丰度分别为0.06%、4.88% 和3.89%、0.70%。中空纤维膜纺丝机通过膜技术进行水处理，应用于制药、酿造、餐饮、化工、市政污水回佣、医院、小区污水会用、造纸等生产生活污水处理。膜分离技术是一种广泛应用于溶液或气体物质分离、浓缩和提纯的分离技术。膜壁微孔密布，原液在一定压力下通过膜的一侧，溶剂及小分子溶质透过膜壁为滤出液，而大分子溶质被膜截留，达到物质分离及浓缩的目的。膜分离过程为动态过滤过程，大分子溶质被膜壁阻隔，随浓缩液流出，膜不易被堵塞，可连续长期使用。膜生物反应器在污水处理，水资源再利用领域，MBR又称膜生物反应器（Membrane Bio-Reactor ）,是一种由膜分离单元与生物处理单元相结合的新型水处理技术。中空纤维膜纺丝机外形像纤维状，具有自支撑作用的膜。它是非对称膜的一种，其致密层可位于纤维的外表面／如反渗透膜，也可位于纤维的内表面(如微滤膜和超滤膜)。对气体分离膜来说，致密层位于内表面或外表面均可。具体参见污水宝商城资料或
3 结论
颗粒硫磺/白云石和黄铁矿/白云石反应器经过30 d的饥饿期后，分别需要5 d和11 d就可使NO3--N去除率恢复饥饿前的水平，且恢复期初期均出现NO2--N的积累现象.两个反应器在短时间内均能恢复稳定运行，在低温条件仍能保持良好的脱氮效果.
高通量测序结果表明，黄铁矿/白云石反应器的细菌群落的丰度和多样性低于硫磺/白云石反应器.经过30 d的饥饿忍耐后，两个反应器的细菌群落的丰度和多样性均略低于恢复期.
微生物群落结构的分析表明，颗粒硫磺/白云石和黄铁矿/白云石硫自养反应器的优势菌门均为GOOGLE PRoteobacteria，主要的纲类均为β-Proteobacteria.颗粒硫磺/白云石反应器中存在起主要反硝化作用的Thiobacillus.
