城市污水處理系統真核微生物群落特性與地域性差異

　　城市污水是人類(lèi)活動(dòng)的必然產(chǎn)物, 其產(chǎn)生量隨著(zhù)城市與社會(huì )經(jīng)濟的發(fā)展不斷增加, 城市污水處理廠(chǎng)(wastewater treatment plant, WWTP)作為城市的“腸道”, 是控制水污染, 保證城市與人類(lèi)可持續發(fā)展的重要設施.我國城市污水處理廠(chǎng)普遍采用活性污泥法處理污水, 其中微生物是污染物去除的主體, 因此微生物群落是系統穩定運行的基礎.

　　在群落中, 除原核生物外, 真菌、原生動(dòng)物、藻類(lèi)等真核微生物作為不可或缺的一部分, 在保障系統穩定運行過(guò)程中發(fā)揮著(zhù)重要的作用.一方面, 原生動(dòng)物可以攝食處理水中游離細菌, 指示活性污泥性質(zhì); 通過(guò)優(yōu)勢種群變化、形態(tài)變化、繁殖方式變化來(lái)指示反應操作環(huán)境; 通過(guò)種群密度和種類(lèi)的變化預測出水水質(zhì); 通過(guò)對細菌、有機顆粒和藻類(lèi)等的強化和捕食作用強化凈化效果; 通過(guò)監測纖毛類(lèi)的多樣性可以預估污水處理廠(chǎng)的運行狀況.另一方面, 隱孢子蟲(chóng)和賈第鞭毛蟲(chóng)等病原性原生動(dòng)物能導致人和動(dòng)物腹瀉, 且具有個(gè)體小、潛伏期長(cháng)、致病性強且易暴發(fā)流行的特點(diǎn), 對下游生物造成極大的威脅.

　　活性污泥中真菌的出現一般也與水質(zhì)有關(guān).一方面, 絲狀真菌是形成活性污泥絮凝體的骨架, 使污泥具有良好的沉淀性能, 保持高的凈化效率、低的出水濃度和出水懸浮物濃度.另一方面, 絲狀真菌, 如Trichosporon屬, 若大量異常地增殖則會(huì )引發(fā)污泥膨脹現象, 常造成處理水質(zhì)下降和污泥的大量流失, 影響處理系統的穩定運行.因此, 深入理解污水處理系統中真核微生物的特性, 對提高出水水質(zhì)、保障污水處理系統的穩定運行具有重要作用.

　　傳統的微生物培養方法和形態(tài)學(xué)研究手段因自身存在的不足和問(wèn)題, 限制了人們對真核微生物的認識, Meng等利用密度梯度離心結合定量蛋白銀染色的方法在黃海48個(gè)站位的沉積物樣品中鑒定出198種底棲纖毛蟲(chóng); 代仁海利用同樣的方法在春季長(cháng)江口及東海的39個(gè)站位的沉積物樣品中檢獲278種底棲纖毛蟲(chóng), 但其中約37%的種類(lèi)由于數量過(guò)少等原因導致分類(lèi)信息不足, 難以準確鑒定.近年來(lái)高通量測序技術(shù)的研究手段極大地提高了微生物群落研究的深度與水平, 現已被廣泛用來(lái)解析農田、森林、土壤、海洋等環(huán)境中的微生物.在污水處理系統中, 相關(guān)研究大多集中于利用高通量測序技術(shù)結合其他生物信息學(xué)方法解析系統內細菌群落的結構和功能, 采用這一技術(shù)研究系統內真核微生物群落的報道還很少見(jiàn).此外, 我國不同區域的地理環(huán)境和飲食結構均存在差異, 且有研究顯示自然系統內真核微生物群落和污水處理系統內的細菌群落均具有地域性差異, 然而不同區域活性污泥系統內的真核微生物群落地域性差異研究尚未開(kāi)展.

　　因此, 本研究擬利用高通量測序技術(shù), 結合多種生態(tài)學(xué)及數理統計相關(guān)分析方法研究我國城市污水處理系統中真核微生物的群落特性與地域性差異, 深化對污水處理系統中真核微生物的認識, 以期為指導污水處理系統的優(yōu)化設計、提高出水水質(zhì)、保障系統穩定高效運行提供科學(xué)依據.

　　1 材料與方法1.1 樣品采集與DNA提取

　　污泥樣品采集于2014年7~8月, 選取了北京, 深圳和無(wú)錫的14個(gè)污水處理廠(chǎng), 為保證后續數據統計學(xué)分析結果的可靠性, 樣品采集遵循每個(gè)城市至少采集4個(gè)WWTPs, 每個(gè)WWTP至少采集3個(gè)曝氣池活性污泥樣品(分別在曝氣池前、中、后端)以及一個(gè)曝氣池上清液樣(曝氣池中端), 若同一WWTP具有不同的處理工藝, 則在不同處理工藝的曝氣池中均進(jìn)行了活性污泥樣品和曝氣池上清液樣品的采集.同時(shí), 原位測定曝氣池混合液的溫度、pH、溶解氧濃度和電導率等參數, 收集WWTPs的地理信息、城市氣候信息, WWTPs采用的處理工藝、設計參數、進(jìn)出水水質(zhì)情況和操作參數等基本信息.

　　采用統一的方法, 每個(gè)樣品采集30 mL, 并迅速置于冰上24 h內帶回實(shí)驗室, 充分搖勻后于15 000g離心10 min后棄去上清液, 取沉淀相保存于-80℃用于微生物DNA提取.采用PowerSoil®DNA提取試劑盒提取污泥樣品的DNA, 然后用ND-2000分光光度計和NanoDrop測定其濃度和純度, 保證每個(gè)樣品滿(mǎn)足測序的DNA質(zhì)量要求.

　　1.2 18S rDNA高通量測序

　　利用引物TAReuk454FWD1(5′-CCAGCASC YGCGGTAATTCC-3′)和TAReukREV3(5′-ACTTTCG TTCTTGATYRA-3′)通過(guò)Illumina MiSeq測序平臺完成活性污泥系統內真核微生物18S rDNA測序, 然后在Galaxy平臺(http://zhoulab5.rccc.ou.edu:8080/root)上完成原始數據的整合, 物種信息分類(lèi), 生成操作分類(lèi)單元(operational taxonomic unit, OTU)表格, .

　　1.3 數據分析

　　本研究使用了非度量多維尺度法(nonmetric multidimensional scaling, NMDS)和相異性分析非參數檢驗方法ADONIS(permutational multivariate analysis of variance)、ANOSIM(analysis of similarities)和MRPP(multi-response permutation procedure)檢驗北京、深圳和無(wú)錫這3個(gè)城市微生物群落結構的差異, 采用多元回歸矩陣分析(multiple regression on distance matrix, MRM)和partial mantel test分析確定影響群落結構的重要因素.以上分析通過(guò)R語(yǔ)言vegan、ape、ecodist程序包和Excel完成.

　　2 結果與討論2.1 OTU及其相對豐度的分布特征

　　為了保證樣品間數據的均等性, 本研究按照每個(gè)樣品8 190條序列的標準對測序結果進(jìn)行隨機抽平處理, 61個(gè)樣品共得到499 590條序列.然后按照97%相似度的標準進(jìn)行序列的OTU劃分, 共獲得3 794個(gè)OTUs.

　　數據分析結果表明, OTU的相對豐度分布極其不均勻, 只有少數的OTU相對豐度較高, 其中, 相對豐度大于0.10%的OTU僅有120個(gè)[圖 1(a)], 這也是常被觀(guān)測到的一種生態(tài)學(xué)現象.此外, 前84個(gè)OTUs(1.07%)的累積相對豐度就已達到總豐度的80%[圖 1(b)].同樣地, OTU在樣品中的分布也十分不均勻[圖 1(c)], 只有4個(gè)OTUs被所有樣品共有, 存在于50或更多個(gè)樣品中的OTU只有37個(gè), 85.40%的OTU僅存在于不超過(guò)10個(gè)樣品中.另外, 從圖 1(c)還可以看出：活性污泥系統真核微生物群落中絕大多數OTU的分布頻率與相對豐度都很低, 但分布頻率較高的OTU傾向于具有較高的相對豐度, 即存在少量分布廣泛且相對豐度較高的OTU, 這與污水處理系統中細菌群落的分布規律相似.

　　圖 1

(a)OTU相對豐度排名-相對豐度; (b)OTU相對豐度排名-累積相對豐度; (c)OTU分布頻率-相對豐度圖 1 OTU豐度-頻率分布示意

　　2.2 活性污泥系統內真核微生物群落的α多樣性

　　α多樣性能夠反映樣品內的微生物群落的豐富度和均勻度.本研究中, 活性污泥系統內每個(gè)樣品中真核微生物群落的OTU數目為212~689, 估算的OTU豐富度(Chao1)為370~1192, 香農-威納指數為1.58~4.80, 辛普森相關(guān)指數為1.8~38.65, 均勻度指數為0.29~0.73.城市水平來(lái)看, 3個(gè)城市活性污泥系統內真核微生物群落多樣性不同(表 1), 其中無(wú)錫群落的物種豐富度(實(shí)際值和估測值)最高, 深圳次之, 北京群落的物種豐富度則明顯低于無(wú)錫和深圳, 香農-威納指數在3個(gè)城市間也呈現相同的規律, 各城市群落物種均勻度差異不大.

表 1 城市樣品水平活性污泥系統內真核微生物群落α多樣性指數

 　　與傳統的研究方法相比, 本研究利用高通量測序技術(shù)檢獲的真核微生物多樣性更高. Matsunaga等通過(guò)構建18S rRNA基因克隆文庫的方法研究了活性污泥中的真核微生物多樣性, 在獲取的834個(gè)克隆中, 60%的克隆很難鑒定, 最終僅獲取了80個(gè)OTUs; Matsubayashi等利用同樣的方法研究了厭氧污泥中的真核微生物多樣性, 發(fā)現系統內85%的真核微生物都是未知的, 樣品中估算的OTU數量(Chao1)僅為34~77個(gè).盡管不同的數據處理軟件、參數的選擇等均可導致OTU數量的改變, 進(jìn)而影響分子多樣性的評價(jià), 且前期PCR擴增和測序等過(guò)程均可能導致微生物群落物種多樣性的高估, 但目前對其高估程度還難以準確定量.

　　與活性污泥系統內的細菌群落比較發(fā)現, 真核微生物多樣性較低. Zhang等使用454焦磷酸測序方法獲得的活性污泥系統內細菌群落OTU數目為1 183~3 567個(gè), Xia等采用Illumina測序方法分別研究了MBR和OD系統內12周的活性污泥樣品內的細菌群落多樣性, 結果表明MBR系統和OD系統樣品的OTU數目分別為2 297~3 330和1 802~3 517個(gè), 均高于真核微生物多樣性.具體聯(lián)系污水寶或參見(jiàn)http://www.sharpedgetext.com更多相關(guān)技術(shù)文檔。

　　2.3 不同空間尺度下的活性污泥系統內真核微生物群落物種組成

　　物種組成是微生物群落研究的一項重要內容, 是探究群落優(yōu)勢菌種或核心物種的基礎, 但不同空間尺度下活性污泥系統內的微生物群落組成可能有所差異.因此, 本研究主要從區域、城市和污水處理廠(chǎng)這3個(gè)空間尺度, Supergroup、Division兩個(gè)水平對OTU所屬的物種信息進(jìn)行分析, 從而全面地揭示活性污泥系統內真核微生物群落的組成特征.

　　在區域尺度上, 在Supergroup水平下, 3 794個(gè)OTUs分屬于10個(gè)Supergroup[圖 2(a)], 其中, Opisthokonta的相對豐度最高, 為54.06%, Alveolata次之, 為35.90%.在Division水平下, 3 794個(gè)OTUs分屬于26個(gè)Division, 其中, 相對豐度大于1.0%的僅有10個(gè)[圖 2(b)], 但它們的累積相對豐度已大于98%.具體地, 87.08%的OTUs屬于Ciliophora、Fungi和Metazoa, 其中, 相對豐度最高的是Ciliophora, 為34.35%, 其次是Fungi和Metazoa, 相對豐度分別為33.02%和19.71%.

　　圖 2

 

(a)整體尺度, Supergroup水平; (b)整體尺度, Division水平; (c)城市尺度, Devision水平; (d)污水處理廠(chǎng)尺度, Devision水平(相對豐度小于1.0%的統一劃歸為others)圖 2 不同尺度/水平下活性污泥系統內真核微生物群落組成 

　　在城市尺度、Division水平下[圖 2(c)], 3個(gè)城市的活性污泥真核微生物群落均以Ciliophora, Fungi和Metazoa為主要組成成分, 總相對豐度為86.22%~89.40%, 但是北京和無(wú)錫群落均以Ciliophora為最主要的組分, 深圳群落則以Fungi為最主要的組分.

　　在污水處理廠(chǎng)尺度、Division水平下[圖 2(d)], 除BJ6外的13個(gè)污水處理廠(chǎng)內的Ciliophora、Fungi和Metazoa的總相對豐度為83.21%~97.30%, BJ6較特殊, 三者總相對豐度僅為74.31%, 其他組分相對豐度較高.分析發(fā)現Perkinsea在各污水處理廠(chǎng)群落內分布極其不均勻, 在BJ6群落內的相對豐度高達19.99%, 但它們在BJ3和SZ4兩個(gè)群落內不存在, 在其他廠(chǎng)內的相對豐度也均很低(0.02%~2.53%). Perkinsea主要包括一些雙殼類(lèi)軟體動(dòng)物的寄生蟲(chóng), 常存活于淡水環(huán)境中.

　　2.4 活性污泥系統內真核微生物群落地域分布差異

　　采用基于Bray-Curtis距離的NMDS方法對來(lái)自于不同污水處理廠(chǎng)的樣品進(jìn)行了群落結構分析, 結果如圖 3所示：整體上, 來(lái)自北京、深圳和無(wú)錫這3個(gè)城市的活性污泥系統內真核微生物群落各自聚為一類(lèi), 表明3個(gè)城市的活性污泥系統內真核微生物群落結構存在差異.相異性分析非參數檢驗方法ADONIS、ANOSIM和MRPP的結果也為3個(gè)城市真核微生物群落結構存在顯著(zhù)差異提供了佐證(P < 0.01).具體看, 來(lái)自某一個(gè)城市不同污水處理廠(chǎng)的活性污泥真核微生物群落傾向于聚類(lèi)在一起, 即來(lái)自同一城市的活性污泥系統內真核微生物群落具有較強的相似性, 不同污水處理廠(chǎng)的活性污泥系統內真核微生物群落結構存在差異, 表明活性污泥系統內真核微生物存在地域性差異.這與自然系統內真核微生物群落和污水處理系統內的細菌群落均存在地域性差異的結論一致.

　　圖 3

 

圖 3 基于NMDS的活性污泥系統內真核微生物群落結構分析

　　2.5 活性污泥系統內真核微生物群落結構的影響因素

　　群落多樣性和結構主要受確定性因素和隨機因素的影響.確定性因素通�；�于生態(tài)位理論描述, 隨機因素由中性理論闡述, 現有研究主要關(guān)注確定性因素中的環(huán)境條件對活性污泥系統內微生物群落結構的影響.基于高通量方法的研究表明WWTPs的水質(zhì)參數、操作條件、地理位置和處理工藝等因素均會(huì )對污水處理系統微生物群落結構產(chǎn)生影響, 研究對象不同, 影響作用較顯著(zhù)的因素種類(lèi)不盡相同.

　　為了探明城市污水處理廠(chǎng)活性污泥系統內真核微生物群落結構的影響因子, 本研究利用partial mantel test分析了地理位置、進(jìn)出水水質(zhì)參數、操作條件等參數與群落結構的相關(guān)性, 結果表明, 在參與分析的20個(gè)因子中, 10個(gè)因子與群落結構顯著(zhù)相關(guān)(P<0.05), 其中相關(guān)系數最大、相關(guān)性最強的因子是出水總氮濃度(r=0.61), 其次是水力停留時(shí)間(r=0.50), 污泥齡(r=0.43)和出水總磷濃度(r=0.38).與出水污染物參數相比, 進(jìn)水污染物參數與群落相關(guān)性較弱[進(jìn)水BOD濃度(r=0.33), 進(jìn)水總磷濃度(r=0.26)], 此外, 與群落結構具有顯著(zhù)相關(guān)性的因素還有地理緯度(r=0.30), 曝氣池溶解氧濃度(r=0.22), 年平均氣溫(r=0.22)和曝氣池混合液濃度(r=0.22).

　　由于partial mantel test只能逐一分析各環(huán)境因子與群落結構的相關(guān)性或將多個(gè)環(huán)境因子作為一個(gè)整體考察其與群落結構的相關(guān)性, 為了探究不同環(huán)境因子對群落結構影響的相對重要性, 本研究使用相同的影響因子組合構建了MRM模型.結果表明, 模型顯著(zhù)且擬合度較高(R2=0.92, P<0.01).對于單個(gè)因子而言(表 2), 曝氣池混合液溫度是造成群落存在差異最重要的因素.環(huán)境溫度影響生物處理構筑物中生物化學(xué)過(guò)程的速度和趨向, 影響微生物的新陳代謝、進(jìn)化速率等生理活動(dòng), 適宜的溫度能夠促進(jìn)、強化微生物的生理活動(dòng), 反之, 能夠減弱甚至破壞微生物的生理活動(dòng), 導致微生物形態(tài)和生理特性的改變, 甚至可能導致微生物死亡, 進(jìn)而影響活性污泥的微生物種群的形成和發(fā)展.

　　表 2 MRM分析各環(huán)境因子對活性污泥系統內真核微生物群落結構影響的權重1)

　

　　與進(jìn)水污染物濃度(進(jìn)水BOD濃度)相比, 出水污染物濃度(出水總氮濃度、出水總磷濃度)與群落結構具有較強的相關(guān)性, 其中出水總氮濃度與群落結構的相關(guān)性強于出水總磷濃度(表 2), Wei等基于ITS數據對污水處理廠(chǎng)中真菌群落結構影響因素研究表明, 出水總氮濃度與群落結構有很強的相關(guān)性.氮是微生物生長(cháng)、繁殖、代謝不可缺少的重要元素, 另外, 出水總氮濃度也是出水水質(zhì)監測的一個(gè)重要指標, 因此這也側面反映了真核微生物作為出水水質(zhì)的指示生物, 其數量和密度與出水水質(zhì)有重要關(guān)聯(lián).齊嶸等采用活體鏡檢計數的方法對無(wú)機廢水硝化過(guò)程中原生動(dòng)物群落結構進(jìn)行了研究, 結果發(fā)現原生動(dòng)物優(yōu)勢種類(lèi)的數量變化與出水NO3--N和NO2--N的濃度存在一定聯(lián)系, 即隨著(zhù)NO2--N的大量累積, 表殼蟲(chóng)生長(cháng)繁殖旺盛, 密度最高可達575 ind.·mL-1, 成為系統中絕對優(yōu)勢的原生動(dòng)物; 但當NO2--N濃度小于10 mg·L-1時(shí), 幾乎未見(jiàn)表殼蟲(chóng)的存在, 而系統中的另一種典型的原生動(dòng)物累枝蟲(chóng)的數量則隨著(zhù)NO3--N濃度的持續增加而逐漸增加.

　　此外, 地理距離也是影響群落結構的重要因子(表 2), 這與上文3個(gè)城市群落結構存在顯著(zhù)差異的結果相一致, 說(shuō)明微生物種群的確存在地域性差異, 事實(shí)上, 地理距離的差異導致了環(huán)境的異質(zhì)性, 如氣溫、降水量的異質(zhì)性, 進(jìn)而導致了群落地域性差異. Bahram等報道了影響表層土壤真菌分布的最主要因素是降水量. Telford等的研究發(fā)現自然系統內的真核微生物(硅藻)會(huì )因為擴散限制而呈有限的地理分布, 其他學(xué)者的研究結果也同樣證明了微生物地域性差異的存在, 且地理距離造成的群落間的差異可能因空間限制而被長(cháng)期保留.

　　此外, 常規污水處理廠(chǎng)操作運行參數曝氣池pH、污泥齡和曝氣池混合液濃度對群落結構也存在重要影響, 它們也是活性污泥系統內細菌群落的重要影響因子.

　　3 結論

　　(1) 不同空間尺度下, 活性污泥系統內真核微生物群落優(yōu)勢種群組成相似.在Devision水平上, 3個(gè)城市真核微生物群落優(yōu)勢種群主要由真菌、原生動(dòng)物(纖毛亞門(mén))和后生動(dòng)物組成.

　　(2) 不同城市污水處理廠(chǎng)活性污泥系統內真核微生物多樣性存在差異, 無(wú)錫群落的多樣性(豐富度和香農-威納多樣性指數)最高, 北京的最低, 且3個(gè)城市的活性污泥系統內真核微生物群落結構存在顯著(zhù)的地域性差異.

　　(3) 利用Partial mantel test和MRM探討了影響活性污泥系統真核微生物群落結構的因素, 結果表明, 曝氣池混合液溫度和出水總氮濃度與活性污泥系統內真核微生物群落結構的相關(guān)性較強.(來(lái)源：環(huán)境科學(xué) 作者：秦文韜)