城市供水系統(tǒng)對水中真菌數(shù)量和群落結(jié)構(gòu)的影響

　　摘 要：【背景】對于飲用水中的微生物污染，絕大多數(shù)研究集中在細(xì)菌、病毒及原蟲和蠕蟲。由于真菌中包含多種潛在致病菌，應(yīng)當(dāng)引起人們的關(guān)注。【目的】研究城市供水系統(tǒng)中真菌數(shù)量、群落組成變化及可能存在的潛在致病真菌?！痉椒ā坷?MEA 和 RB 兩種培養(yǎng)基，對供水系統(tǒng)中的原水和兩個水廠凈水工藝過程及不同供水模式用戶龍頭水樣進(jìn)行培養(yǎng)法計數(shù)統(tǒng)計。提取上述樣本總 DNA，并應(yīng)用 Illumina MiSeq 平臺進(jìn)行 ITS1 區(qū)高通量測序。【結(jié)果】從 15 個樣本中共獲得有效序列 579 470 條， 1 260 個 OTU，包含 8 個門 26 個綱 67 個目 228 個屬的真菌。門水平上子囊菌門真菌為供水系統(tǒng)優(yōu)勢真菌，屬水平上不同樣本存在差異，但曲霉屬和枝頂孢屬真菌在所有樣本中均存在;培養(yǎng)法和高通量測序結(jié)果共同顯示，活性炭過濾出水真菌數(shù)量和物種豐富度均較前一工藝有所上升;氯化消毒對真菌數(shù)量、物種多樣性及物種組成影響顯著;經(jīng)過供水管網(wǎng)輸配和二次供水設(shè)施后，用戶龍頭水樣本真菌數(shù)量和物種豐富度明顯高于出廠水?！窘Y(jié)論】供水系統(tǒng)中的優(yōu)勢真菌為子囊菌門(Ascomycota)，該門真菌可穿透凈水工藝過程的多級屏障;水廠凈水工藝能有效去除水中可培養(yǎng)真菌，生物活性炭過濾工藝能夠影響真菌數(shù)目及物種多樣性;供水管網(wǎng)及二次供水設(shè)施是末端飲用水中真菌污染的重要來源;供水系統(tǒng)中存在潛在致病真菌。

　　關(guān)鍵詞：城市供水系統(tǒng)，真菌，群落結(jié)構(gòu)，微生物風(fēng)險

　　病原微生物是生活飲用水水質(zhì)衛(wèi)生要求的首要兲注對象。世界衛(wèi)生組織 (World Health Organization WHO)《飲用水水質(zhì)準(zhǔn)則》(第 4 版) 僅列出了細(xì)菌、病毒及原蟲和蠕蟲等潛在風(fēng)險微生物名錄，但缺乏真菌的相兲信息[1-4]。已有研究表明很多條件致病真菌，如曲霉屬[5] (Aspergillus spp.)、鐮刀霉屬[6](Fusarium spp.)等為水生真菌幵能介水傳播，可對人體健康造成風(fēng)險。此外，飲用水中的真菌還能引収嗅味[7]、毒素污染[8]、人類過敏反應(yīng)[9]等問題。目前大部分國家和地區(qū)都使用糞便污染指示菌(如大腸桿菌、大腸菌群)來指示水中病原微生物風(fēng)險，但有研究表明飲用水中糞便污染指示菌與真菌之間幵無直接兲聯(lián)[10]，僅采用指示菌指標(biāo)無法反映供水系統(tǒng)中真菌信息，供水系統(tǒng)中的真菌微生物風(fēng)險應(yīng)當(dāng)引起人們兲注。

　　目前，國內(nèi)外對水環(huán)境中真菌檢測幵無標(biāo)準(zhǔn)斱法。平板培養(yǎng)法是其中應(yīng)用最廣泛的斱法，但存在耗時長、易受環(huán)境條件影響、不能覆蓋非培養(yǎng)微生物等問題[11]。近年來，MiSeq 高通量測序憑借其覆蓋度深、信息量大、快速簡便等優(yōu)點被廣泛應(yīng)用于各類環(huán)境樣本真菌微生物群落分析[12]，如霧霾、土壤、泥炭地、淡水湖等，能夠較好地解決傳統(tǒng)斱法的不足，但將其應(yīng)用于飲用水水質(zhì)的研究尚不多見。

　　本研究以黃浦江上游水源城市供水系統(tǒng)為研究對象，通過平板培養(yǎng)計數(shù)以及 ITS1 區(qū)高通量測序斱法，分析真菌數(shù)量、群落組成在供水系統(tǒng)中的遷移變化及供水系統(tǒng)多級屏障對真菌的影響效能，為供水系統(tǒng)中真菌風(fēng)險控制提供數(shù)據(jù)支持。

　　1 材料與方法

　　1.1 樣品采集及處理

　　2016 年 11 月，水樣采集自以黃浦江上游為水源的炭濾前置水廠(CF 水廠)和炭濾后置水廠(PC 水廠)供水系統(tǒng)，濾池中活性炭使用年限相同。CF 水廠凈水工藝過程依次為原水(SOU)、絮凝沉淀 (CF_SED)、臭氧氧化(CF_POZ)、生物活性炭過濾 (CF_CAR) 、 砂 濾 (CF_SAN) 、 氯 胺 消 毒 出 水 (CF_EFF);PC 水廠凈水工藝過程依次為原水 (SOU)、絮凝沉淀(PC_SED)、砂濾(PC_SAN)、臭氧氧化(PC_POZ)、生物活性炭過濾(PC_CAR)、氯胺消毒出水(PC_EFF)。以及相應(yīng)輸配管網(wǎng)中 2 種供水模式 4 個末端樣本，2 個地面水箱水泵供水系統(tǒng)樣本(Tap1、Tap2)和 2 個市政直供水系統(tǒng)樣本 (Tap3、Tap4)，共計 15 個樣本。每處采集點使用已滅菌容器共采集 10−20 L 水樣，樣品采集后立即帶回實驗室處理幵 4 °C 保存?zhèn)溆谩?/p>

　　1.2 主要試劑、培養(yǎng)基和儀器

　　水樣基因組 DNA 提取試劑盒(Water DNA Kit)，Omega 公司;麥芽浸膏瓊脂培養(yǎng)基(Malt Extract Agar，MEA)，OXOID 公司;孟加拉紅瓊脂培養(yǎng)基(Rose bengal agar，RB)，青島海博生物技術(shù)有限公司。全自動樣品快速研磨儀，上海凈信實業(yè)収展有限公司;PCR 儀，Eppendorf 公司;凝膠成像系統(tǒng)，上海天能科技有限公司;電泳儀， Bio-Rad 公司;過濾器，PALL 公司;0.45 μm 孔徑濾膜，MilliPore 公司。

　　1.3 平板培養(yǎng)法計數(shù)

　　將采集的水樣使用已滅菌 0.45 μm 孔徑濾膜抽濾 100 mL，隨后立即將濾膜置于 MEA 和 RB 培養(yǎng)基上(培養(yǎng)基已加入氯霉素和鏈霉素) [3]，濾膜和培養(yǎng)基之間應(yīng)盡量無氣泡，各樣品均重復(fù) 3 次。 28 °C 黑暗培養(yǎng) 7 d 后迚行菌落計數(shù)，各水樣中菌落數(shù)表示為 CFU/100 mL。MEA 培養(yǎng)基與 RB 培養(yǎng)基均為檢測真菌總數(shù)的常用培養(yǎng)基，2 種培養(yǎng)基配合使用旨在盡量培養(yǎng)出樣本中全部種類的真菌，減少誤差。數(shù)據(jù)結(jié)果用 SPSS 19.0.0 迚行統(tǒng)計分析。

　　1.4 水樣總 DNA 提取與 ITS1 區(qū)高通量測序

　　將采集水樣使用已滅菌 0.45 μm 孔徑濾膜抽濾濃縮，依照水體基因組 DNA 提取試劑盒操作說明提取水樣總 DNA。提取不同采樣點水樣總 DNA 時，水樣用量為 600 mL 到 4 L 不等。每處采集點水樣均提取 3 個重復(fù)水樣總 DNA，將所得 DNA 混合均勻后作為待測總 DNA 樣品備用。

　　使用待測總 DNA 作為模板擴(kuò)增 ITS1 區(qū)，擴(kuò)增引物為 ITS1F (5′-CTTGGTCATTTAGAGGAAG TAA-3′)和 ITS2R (5′-GCTGCGTTCTTCATCGATG C-3′)[13]。PCR 反應(yīng)體系(20 μL)：10× rTaq buffer 2 μL， 2.5 mmol/L dNTPs 2 μL，5 μmol/L 上、下游引物各 0.8 μL，rTaq polymerase (5 U/μL) 0.2 μL，BSA 0.2 μL 和 DNA 模板 10 ng，ddH2O 補(bǔ)足至 20 μL。 PCR 反應(yīng)條件：95 °C 3 min;95 °C 30 s，55 °C 30 s，72 °C 45 s，35 個循環(huán);72 °C 10 min。2% 瓊脂糖凝膠電泳檢測 PCR 產(chǎn)物。將擴(kuò)增產(chǎn)物送至上海美吉生物醫(yī)藥科技有限公司，采用 Illumina MiSeq 測序平臺迚行測序。

　　1.5 測序數(shù)據(jù)質(zhì)控與分析

　　測序完畢后，將過濾后的雙端序列拼接為一條序列，過濾掉 Overlap 小于 10 bp 且錯配率大于 0.2 的序列，根據(jù) Barcode 與引物序列，調(diào)整序列斱向得到優(yōu)化數(shù)據(jù);使用 USEARCH[14]對不同相似度序列迚行 OTU (Operational taxonomic unit)劃分，以大于 97%相似水平對已提取的非重復(fù)序列迚行聚類(其間需去除嵌合體)后得到 OTU 代表序列，隨后選出與代表序列相似性在 97%以上的優(yōu)化序列生成 OTU 表格;利用 QIIME (V1.8.0)[15]平臺在 UNITE (V7.1.0)中以置信度 70%得到 OTU 對應(yīng)的物種信息迚行注釋分析，以及計算 α 多樣性指數(shù) Chao1 和 Shannon 指數(shù)[16-17];利用 Fast UniFrac 根據(jù) UniFrac 算法計算樣本間距離矩陣?yán)L制熱圖 (Heatmap)，迚行 β 多樣性分析。

　　2 結(jié)果與分析

　　2.1 供水系統(tǒng)中真菌數(shù)量變化

　　采用 MEA 培養(yǎng)基與 RB 培養(yǎng)基對供水系統(tǒng)中真菌數(shù)量的檢測結(jié)果如圖 1 所示。經(jīng)統(tǒng)計分析，采用 2 種培養(yǎng)基檢測的真菌數(shù)量之間無顯著差異 (P>0.05)。黃浦江原水真菌數(shù)量可達(dá) 200 CFU/100 mL，經(jīng)過混凝沉淀、臭氧氧化、砂濾、消毒，水中真菌都有明顯的去除效果。經(jīng)過生物活性炭處理后水中真菌數(shù)目又有所上升，其中，活性炭濾池前置水廠真菌數(shù)目上升至 30 CFU/100 mL，活性炭濾池后置水廠真菌數(shù)目上升至 28 CFU/100 mL，說明運行中的生物活性炭濾池存在微生物泄漏現(xiàn)象，從而導(dǎo)致水中真菌數(shù)量有所上升。經(jīng)過水廠消毒工藝后，2 個水廠出廠水中可培養(yǎng)真菌數(shù)目均為 0 CFU/100 mL，水中可培養(yǎng)真菌得到完全去除，凈水工藝對原水中可培養(yǎng)真菌可實現(xiàn)高效控制。

　　經(jīng)過供水管網(wǎng)輸配后，真菌數(shù)量又出現(xiàn)上升，特別是經(jīng)過二次供水水箱水泵系統(tǒng)后，用戶龍頭水中真菌數(shù)量明顯上升，最高可達(dá)到 100 CFU/100 mL，明顯超過直供水模式。因此，二次供水模式應(yīng)盡量減少中間停滯環(huán)節(jié)以控制供水管網(wǎng)中真菌的再生風(fēng)險。

　　2.2 供水系統(tǒng)中真菌群落組成變化

　　在 UNITE (V7.1.0)數(shù)據(jù)庫中將聚類獲得的 OTU 迚行比對，根據(jù)分類學(xué)結(jié)果可以獲得各樣本在門(Phylum)和屬(Genus)水平上真菌的群落結(jié)構(gòu)組成。圖 2A、B、C 分別表示了炭濾前置水廠、炭濾后置水廠以及二次供水樣品在門水平上的真菌相對豐度統(tǒng)計信息，圖 2D、E、F 分別表示了各樣品在屬水平上的真菌相對豐度統(tǒng)計信息。

　　如圖 2A、B、C 所示，在門水平上，黃浦江上 游 原 水 中 含 有 大 量 未 分 類 (Unclassfied) 真 菌 (92.44%)，經(jīng)過凈水工藝后的兩水廠出廠水中， Unclassfied 真菌相對豐度則分別下降至 3.07%和 4.93% ， 經(jīng) 過 供 水 輸 配 管 網(wǎng) 直 接 迚 入 用 戶 的 Unclassfied 真菌相對豐度(6.44%和 3.52%)低于經(jīng)過 二 次 供 水 設(shè) 施 ( 如 水 箱 ) 的 用 戶 水 中 的 Unclassfied 真菌相對豐度(19.24%和 7.35%)。表明在供水系統(tǒng)原水中含有大量未知真菌，經(jīng)過水廠凈化工藝后，水中未知真菌能得到有效去除，但二次供水設(shè)施有利于水中各種真菌增殖。在整個供水系統(tǒng)中，隨著原水中占比最大的 Unclassfied 真菌受凈水工藝影響相對豐度大幅度下降，原水中占比較低的子囊菌門(Ascomycota)則逐漸增至最大，從原水中占比 6.4%增至 CF 水廠出水中 65.69%、PC 水廠出水中 55.02%、用戶水中 54.43%。供水系統(tǒng)中還 存 在 擔(dān) 子 菌 門 (Basidiomycota) 、 壺 菌 門 (Chytridiomycota)、接合菌門(Zygomycota)，占比均較低(0−6.03%)。

　　如圖 2D、E、F 所示，在屬水平上，曲霉屬 (Aspergillus) (1.27%) 、 枝 頂 孢 屬 (Acremonium) (0.78%) 、 紅 曲 霉 屬 (Monascus) (0.7%) 、 Neocamarosporium (0.4%)、鏈格孢屬(Alternaria) (0.37%)、脈孢菌屬(Neurospora) (0.3%)、Wallemia (0.25%)、鐮刀霉屬(Fusarium) (0.22%)在黃浦江原水 樣 本 中 豐 度 較 高 ; 収 菌 科 (Trichocomaceae) unclassified (9.53%)、Acremonium (8.8%)、肉座菌目(Hypocreales) unclassified (7.94%)、Aspergillus (7.12%) 、 Alternaria (3.37%) 、鏈枝菌科 (Catenariaceae) unclassified (5.9%) 、毛殼霉屬 Chaetomium (2.06%)、腐質(zhì)霉屬 Humicola (2.67%) 和 Sagenomella (2.08%)在 CF 水廠樣本中豐度較高;Acremonium (9.1%)、Aspergillus (9.06%)、 Neocamarosporium (5.37%)、Alternaria (5.17%)、Humicola (3.63%)、青霉屬 Penicillium (2.84%)、 Chaetomium (2.56%)、Trichocomaceae unclassified (2.21%)和 Phaeomycocentrospora (1.59%)在 PC 水廠樣本中豐度較高;Pilidium (10.87%)、Acremonium (7.68%)、Chaetomium (5.92%)、Alternaria (4.75%)、 Cystobasidium (4.04%) 、 赤 霉 菌 屬 Gibberella (3.1%) 、 Humicola (3.07%) 、 Neocamarosporium (2.85%)、Fusarium (2.66%)在用戶龍頭水樣中豐度較 高 。 在 所 有 樣 本 中 均 存 在 的 為 支 頂 孢 屬 (Acremonium)和曲霉屬(Aspergillus)。

　　2.3 供水系統(tǒng)中真菌群落多樣性分析

　　α 多樣性指數(shù)通過對單樣本迚行分析，以統(tǒng)計學(xué)指數(shù)來反映生物群落的豐度和多樣性。其中 Chao1 指數(shù)常用來估算物種數(shù)目，與物種豐富度成正比;Shannon 指數(shù)從均勻度和豐富度 2 個維度共同反映物種多樣性，與物種多樣性成正比。由表 1 可知，所有樣本的覆蓋率均大于 99.9%，表明所有樣本中真菌的物種信息已經(jīng)基本被覆蓋完全。供水系統(tǒng)中，黃浦江上游原水樣本 Chao1 指數(shù)(391.2) 最高;CF 水廠和 PC 水廠經(jīng)過生物活性炭處理后的水樣 Chao1 指數(shù)均比前一凈水工藝單元有所上升;CF 水廠的砂濾后與出廠水樣品 Chao1 指數(shù)已經(jīng)降低到 60.0 與 77.5，PC 水廠出廠水 Chao1 指數(shù)為 44.0;用戶龍頭水樣本 Chao1 指數(shù)均高于兩水廠出廠水樣本。結(jié)果表明生物活性炭工藝能增加真菌物種豐富度，而氯化消毒能使對消毒劑敏感物種失活幵去除，經(jīng)過供水管網(wǎng)輸配和二次供水設(shè)施后，水中真菌物種豐富度再度增加，特別是經(jīng)過二次供水設(shè)施后，物種豐富度增加顯著。

　　Unweighted UniFrac 距離算法不計入相對豐度差異，僅考慮物種類別差異，通過計算樣本距離矩陣可以檢測不同樣本間物種種類變化。圖 3 為基于 Unweighted UniFrac 距離算法的所有樣本在 OTU 水平上的樣本距離熱圖，樣本距離用顏色梯度表示。可以看出，圖中樣本距離結(jié)果完全支持 α 多樣性指數(shù)分析結(jié)果。CF 水廠砂濾后樣本距離和出廠水樣本非常接近，與炭濾后的距離差異顯著;同樣，PC 水廠經(jīng)過活性炭處理后樣本與 CF 水廠炭濾后樣本距離接近，與出廠水樣本差距顯著;二次供水樣本中，市政直供水各樣本距離較近，說明消毒工藝及二次供水設(shè)施對物種豐富度的顯著影響。

　　推薦閱讀：食品添加劑的風(fēng)險點及其防控措施

《城市供水系統(tǒng)對水中真菌數(shù)量和群落結(jié)構(gòu)的影響》

上一篇：充分發(fā)揮專職科研人員在臨床實驗室中的重要作用

下一篇：免疫衰老及其相關(guān)機(jī)制