市政污水氧化溝處理工藝對人類致病菌群落組成的影響以南京市某污水處理廠為例
——

吳璟涵，劉鵬，張徐祥，任洪強

(污染控制與資源化研究國家重點實驗室，南京大學環境學院，江蘇 南京 210023)

介水傳染病是由于飲用或接觸受病原體污染的水而引起的一類傳染病，其病原體主要包括細菌、病毒、原生動物等[1]。近年來，全球范圍內介水傳染病的暴發已經嚴重危害各國人民生命健康，據世界衛生組織報道，介水傳染病每年造成220多萬人死亡[1-2]。市政污水處理廠作為城市水循環的重要節點，既通過污水管網承接城市污水，又將處理后的污水排放進入自然水體，是人類致病菌(Human pathogenic bacteria，HPB)傳播的重要源和匯。

以活性污泥為代表的生化處理工藝是市政污水處理廠的核心，是污水中有機質和營養鹽去除的關鍵單元。隨著高通量測序技術被廣泛應用于污水和污泥細菌群落結構的解析[3]，生化處理工藝在塑造污水細菌群落結構尤其是HPB群落結構中的作用得到越來越多關注。污水中微生物種類繁多[4]，包括志賀氏菌[5]、銅綠假單胞菌[6]、沙門氏菌[7]、結核桿菌[8]、糞腸球菌[9]、鮑氏不動桿菌[10]等多種常見HPB，形成了復雜的HPB群落結構。Ye等[11]利用454焦磷酸測序技術在來自4個國家14座市政污水處理廠中污泥、進水、出水中檢測到34個不同屬的HPB。研究同時發現，部分HPB在生化處理前后的進水和出水樣品中豐度相差較大，如弓形菌屬在進水中相對豐度較高，而弧菌屬在出水中相對豐度較高。Li等[12]通過宏基因組學方法在香港某污水處理廠中檢測到113種HPB，其中79種在該污水處理廠首次檢測出，而活性污泥工藝對HPB的總體去除率達到98%，但不同種類的HPB去除效果存在較大差異。Huang等[13]利用454焦磷酸測序技術在某污水處理廠中檢測了243種潛在致病菌，活性污泥工藝顯著降低了污水HPB群落的多樣性。

氧化溝處理工藝是一種常見的活性污泥法變形工藝，具有工藝簡單，運行靈活，管理方便，成本較低的特點，被廣泛運用于中小型污水處理廠[14]，但氧化溝處理工藝對污水HPB群落組成的影響卻鮮有報道[15]?，F選擇南京市某采用氧化溝處理工藝的城市污水處理廠為研究對象，采集2016年1— 12月份污水處理廠的進水、污泥和生化處理出水樣品，基于16S rRNA基因的聚合酶鏈式反應(PCR)擴增和高通量測序技術，通過致病菌數據庫比對的方法解析了HPB群落組成特征，并研究了生化處理過程中水質理化指標與HPB群落組成的相關性，從而有助于加深對污水處理系統HPB群落的動態變化規律的理解，為氧化溝處理工藝去除HPB提供科學支撐。

1 材料與方法

1.1 樣品采集

選擇南京市某城市污水處理廠為研究對象，該污水處理廠除收集生活組團的污水外，同時處理部分工業企業污水，日處理水量約為4.5萬m3。該污水處理廠主體工藝采用改進型三槽式氧化溝為生化處理工藝，水力停留時間為18.3 h，污泥齡為18.1 d。于2016年1—12月，每月采集該污水處理廠進水、活性污泥與生化處理出水樣品。使用無菌聚乙烯桶采集進水樣品200 mL和生化處理出水樣品20 L，使用50 mL無菌離心管采集活性污泥樣品，并立即以1∶1(V/V)比例將污泥樣品與無水乙醇混合，以固定生物質。將上述樣品置于冰上并運送至實驗室保存于-20 ℃冰箱。

1.2 水質理化指標測定

1.3 樣品總DNA提取、PCR擴增與高通量測序

進水和生化處理出水樣品分別經0.45 μm微孔濾膜過濾，以富集生物質用于總DNA提取。采用FastDNA SPIN試劑盒[安倍醫療器械貿易(上海)有限公司]提取樣品總DNA，并采用NanoDrop?ND-2000超微量紫外分光光度計(美國基因有限公司)測定總DNA的濃度與純度。

總DNA樣品送至北京諾禾致源科技股份有限公司，擴增16S rRNA基因v4可變區，采用Hiseq 2500平臺進行高通量測序，測序策略為PE250。

1.4 生物信息學分析

生物信息學分析采用QIIME2平臺進行。首先，采用DADA2對原始測序序列降噪和去除冗余序列，得到擴增子序列變異(amplicon sequence variant, ASV)。其次，選擇每個ASV代表序列，采用Blast方法比對包括128種人類致病菌16S rRNA基因全長序列的數據庫，比對閾值E設置為10-3，將序列相似度超過97%的ASV鑒定為潛在HPB。最后，HPB的相對豐度和某一樣品中HPB總相對豐度按照公式(1)和(2)[16]計算：

(1)

(2)

式中：Ai——某一樣品中某一HPB的相對豐度；NHPB-like sequence——比對上的某HPB的序列條數；

N16S sequence——樣品中微生物總序列條數；AT——某一樣品中HPB的總相對豐度；n——某一樣品中對應檢出的HPB群落數量。

1.5 統計分析

采用R語言中vegan包進行HPB群落與水質指標的冗余分析(RDA)，并通過計算Spearman系數反映HPB群落豐度與水質指標的相關性，Spearman相關性系數絕對值|r|>0.80且P<0.05，即認為具有顯著相關性。

2 結果與討論

2.1 污水處理廠的水質理化指標

污水處理廠在2016年1—12月，進水水質呈現一定波動，但生化處理出水水質較為穩定，均符合《城鎮污水處理廠污染物排放標準》(GB 18918—2002)中一級A排放標準，水質理化指標情況見表1。

表1 污水處理廠進水、生化處理出水水質理化指標

2.2 進水、污泥和生化處理出水中HPB檢出種類數量與總相對豐度

利用基于16S rRNA基因的高通量測序技術，在所有進水、污泥和生化處理出水樣品中共檢測出28種HPB群落。3種樣品中HPB檢出種類數量與總相對豐度的年內變化見圖1(a)(b)。通過與數據庫中128種潛在HPB的信息對比發現，進水、生化處理出水中HPB的檢出種類數量較為平穩，基本穩定在14 ～ 20種。污泥中HPB的檢出種類數量變化較大，冬、春季節較少，其中1，3月最少，檢出12種；夏、秋季節較多，其中9月最多，檢出23種。

同時，進水和生化處理出水中HPB的總相對豐度變化不大，平均值為3%左右。而污泥中HPB的含量隨季節的變化差異明顯，2月污泥中HPB總相對豐度最低，為1.7%，5月最高，為15.1%。此外，污泥中HPB的總相對豐度普遍比進水、生化處理出水要高，顯示出活性污泥對HPB具有一定的吸附、蓄積潛力。

圖1 進水、污泥和生化處理出水中HPB檢出種類數量與總相對豐度的年內變化

2.3 進水、污泥和生化處理出水中HPB群落組成

本研究共檢出28種HPB，檢出類型與年內相對豐度變化見圖2，相對豐度值已通過對數標準化處理。其中，共有6種HPB在所有月份的進水、污泥和生化處理出水中均有檢出且具有較高的相對豐度，分別是鮑氏志賀菌(Shigellaboydii)(0.01%～0.28%)、腸道沙門氏菌(Salmonellaenterica)(0.05%～1.25%)、銅綠假單胞菌(Pseudomonasaeruginosa)(0.15%～2.46%)、結核分枝桿菌(Mycobacteriumtuberculosis)(0.01%～0.31%)、糞腸球菌(Enterococcusfaecalis)(0.07%～3.54%)、鮑氏不動桿菌(Acinetobacterbaumannii)(0.22%～9.80%)。這6種HPB在市政污水[17]與湖泊[18]、河流[19]等自然水體中均多有報道，是水環境中常見的HPB。

圖2 進水、污泥和生化處理出水中HPB組成的年內變化

在上述6種HPB中，鮑氏志賀菌、糞腸球菌和鮑氏不動桿菌均呈現污泥中相對豐度高于進水和生化處理出水的分布特征，表明這3種HPB環境存活能力較強，具有在活性污泥中富集的潛力。腸道沙門氏菌、銅綠假單胞菌和結核分枝桿菌則在進水中相對豐度較高，經氧化溝工藝去除后，活性污泥和生化處理出水中相對豐度則均出現下降。

由圖1可見，5和9月的進水與生化處理出水樣品檢出的HPB種類個數相同，但對比發現，2類樣品中檢出的HPB種類并不完全相同。首先，在5月的樣品中，巴頓氏菌(Bartonellahenselae)、布魯氏菌(Brucellasuis)、白喉桿菌(Corynebacteriumdiphtheriae)在進水中檢出，經過氧化溝生化處理后，未在出水中檢出。而產氣莢膜梭菌(Clostridiumperfringens)、杰氏棒桿菌(Corynebacteriumjeikeium)、肺炎支原體(Mycoplasmapneumoniae)只在生化處理出水樣品中檢出。其次，在9月的樣品中，流感嗜血桿菌(Haemophilusinfluenzae)、嗜肺軍團菌(Legionellapneumophila)、肺炎鏈球菌(Streptococcuspneumoniae)在進水中檢出，經過氧化溝生化處理后，未在生化處理出水中檢出。而蠟樣芽孢桿菌(Bacilluscereus)、多形擬桿菌(Bacteroidesthetaiotaomicron)、無乳鏈球菌(Streptococcusagalactiae)只在生化處理出水中檢出。2個月份中不同的HPB種類并無重合，說明不同季節進水的多變性和復雜性往往直接影響生化處理的效果以及出水中的HPB群落變化。

同時，在生化處理出水中，HPB的相對豐度還表現出一定的月份差異。糞腸球菌在4—7月，其相對豐度達到1.5% ～ 3.5%，遠超出年度平均值1.0%。而某些HPB則在秋、冬季節的生化處理出水中相對豐度較高。鮑氏志賀菌在9—11月，相對豐度分別為0.08% ～ 0.13%，也高于年內其余月份的相對豐度。結核分枝桿菌則在8—10月相對豐度較高。此外，鮑氏不動桿菌在4—5月，以及8—11月出現了2個高相對豐度的時間段。

在生化處理出水中，蠟樣芽孢桿菌(Bacilluscereus)、鮑氏志賀菌(Shigellaboydii)、脆弱類桿菌(Bacteroidesfragilis)、副百日咳博德特菌(Bordetellaparapertussis)等10種HPB在超過6個月的樣品中檢出。此外，還有5種HPB僅在少數月份的生化處理出水中檢出，主要包括多形擬桿菌(Bacteroidesthetaiotaomicron)、杰氏棒桿菌(Corynebacteriumjeikeium)等。

2.4 生化處理出水理化指標與HPB群落組成的相關性

圖3 生化處理出水理化指標與HPB群落組成的冗余分析

表2 出水HPB豐度與理化指標的Spearman相關性系數

3 結論

(1)進水、生化處理出水中HPB的檢出種類數量較為平穩，總相對豐度變化不大。污泥中HPB的檢出種類數量變化較大，夏、秋季多于冬、春季；HPB總相對豐度隨季節變化有明顯差異，2月總相對豐度最低，5月最高。污泥中HPB的總相對豐度普遍高于進水、生化處理出水。

(2)進水、污泥、生化處理出水3種樣品中檢出的HPB具有相似性。其中，優勢HPB包括鮑氏志賀菌、腸道沙門氏菌、銅綠假單胞菌、結核桿菌、糞腸球菌和鮑氏不動桿菌。