關(guān)鍵詞：豐水期；長(zhǎng)蕩湖；水質(zhì)評(píng)價(jià)；微生物群落

中圖分類(lèi)號(hào)：X171.4 文獻(xiàn)標(biāo)志碼：B

前言

微生物是水生態(tài)系統(tǒng)的重要組成部分，在維持水生態(tài)系統(tǒng)的功能和健康中起著重要作用，同時(shí)微生物也是反映水環(huán)境健康的重要指標(biāo)。張建波等認(rèn)為后生浮游動(dòng)物可作為湖泊水環(huán)境變化的指示性生物。由于微生物能將對(duì)生態(tài)環(huán)境造成威脅的有機(jī)物和有毒物質(zhì)轉(zhuǎn)化為可以維護(hù)生態(tài)平衡的無(wú)機(jī)物。如枯草芽孢桿菌是一類(lèi)革蘭氏陽(yáng)性菌，其分泌的多種高活性胞外酶，能快速分解并吸收水環(huán)境中的有機(jī)物，以達(dá)到改善水質(zhì)的效果。因此全面調(diào)查長(zhǎng)蕩湖微生物群落組成是極有必要的。

長(zhǎng)蕩湖位于江蘇省常州市。長(zhǎng)蕩湖南北長(zhǎng)約15.5km，東西寬約9km，湖泊總面積85.31km²，平均水深1.1m。隨著長(zhǎng)蕩湖周邊農(nóng)業(yè)、漁業(yè)的快速發(fā)展以及過(guò)度城鎮(zhèn)化，導(dǎo)致湖泊水體嚴(yán)重富營(yíng)養(yǎng)化，水質(zhì)日漸惡化。

目前有關(guān)長(zhǎng)蕩湖的研究主要聚焦于藻類(lèi)、動(dòng)物的群落結(jié)構(gòu)特征以及沉積物污染風(fēng)險(xiǎn)等方面，但鮮有全面研究長(zhǎng)蕩湖微生物群落結(jié)構(gòu)組成。鑒于豐水期雨量充沛，水樣采集方便且豐水期長(zhǎng)蕩湖的水體污染更為嚴(yán)重。因此此研究基于2021年6月長(zhǎng)蕩湖采樣數(shù)據(jù)，采用2種水質(zhì)評(píng)價(jià)方法綜合評(píng)價(jià)豐水期長(zhǎng)蕩湖水質(zhì)現(xiàn)狀，并運(yùn)用高通量測(cè)序技術(shù)分析長(zhǎng)蕩湖微生物群落結(jié)構(gòu)特征，為控制長(zhǎng)蕩湖水體污染和生態(tài)修復(fù)提供理論依據(jù)。

1材料與方法

1.1采樣點(diǎn)設(shè)置與采樣時(shí)間

根據(jù)長(zhǎng)蕩湖地形輪廓，在長(zhǎng)蕩湖南部、中部、東部和北部共設(shè)置4個(gè)采樣點(diǎn)。4個(gè)采樣點(diǎn)編號(hào)依次為W156、W158、W160、W163，于2021年6月（豐水期）進(jìn)行水體樣品采集。采樣點(diǎn)的分布見(jiàn)圖1。

1.2水樣采集與測(cè)定方法

使用有機(jī)玻璃采水器采集長(zhǎng)蕩湖表層水，每個(gè)采樣點(diǎn)采集1L水樣，保存于已消毒滅菌并用樣本底水充分清洗的聚乙烯取樣瓶中，并在4℃條件下運(yùn)回實(shí)驗(yàn)室。其中500 mL水樣用于理化因子測(cè)定，500 mL在實(shí)驗(yàn)室無(wú)菌環(huán)境下使用0.45μm濾膜過(guò)濾后，將濾膜放入10mL凍存管并保存于-80℃冰箱，用于后續(xù)DNA測(cè)序。

此次測(cè)定理化指標(biāo)包括總有機(jī)碳（TOC）、總氮（TN）、總磷（TP）、水溫（WT）、pH值及溶解氧（DO）。TOC采用燃燒氧化-非分散紅外吸收法；TN采用堿性過(guò)硫酸鉀-紫外分光光度法；TP采用鉬酸銨分光光度法；采用YSI多參數(shù)水質(zhì)分析儀（YSI PR01020，USA）現(xiàn)場(chǎng)測(cè)定WT、pH和DO。

1.3微生物群落測(cè)序分析方法

先在無(wú)菌環(huán)境中將濾膜剪成碎片，然后使用PowerWater DNA（MOBIO，USA）試劑盒提取基因組DNA。通過(guò)1%瓊脂糖凝膠電泳檢測(cè)提取的基因組DNA。使用正向引物341F（5’-CCTAYGGGRBG-CASCAG -3’）和反向引物806R（5’-GGACTAC-NNGGGTATCTAAT-3’）對(duì)V3-V4可變區(qū)進(jìn)行PCR擴(kuò)增。使用2%瓊脂糖凝膠回收PCR產(chǎn)物，利用AxyPrep DNA Gel Extraction Kit（Axygen Biosci-ences，Union City，CA，USA）進(jìn)行純化，Tris-HCl洗脫，2%瓊脂糖電泳檢測(cè)。測(cè)序采用IlluminaMiSeq PE300高通量測(cè)序平臺(tái)進(jìn)行測(cè)序。PCR及測(cè)序均由上海凌恩生物科技公司完成。

1.4數(shù)據(jù)處理與分析

水質(zhì)評(píng)價(jià)方法分別采用水污染指標(biāo)法和內(nèi)梅羅指數(shù)法2種評(píng)價(jià)方法。

1.4.1水污染指數(shù)法

水污染指數(shù)法以單因子指數(shù)法為基礎(chǔ)，內(nèi)插計(jì)算出各污染指標(biāo)的WPI值，選取WPI值的最高值來(lái)判定水體污染等級(jí)。選取的環(huán)境參數(shù)為總氮、總磷、DO，每個(gè)參數(shù)的WPI值計(jì)算公式如式（1）：

基于Microsoft Excel處理的原始實(shí)驗(yàn)數(shù)據(jù)；微生物多樣性指數(shù)使用QIIME軟件（1.7.0版）計(jì)算，使用Origin 2018軟件繪制物種豐度柱狀圖。依托微信在線平臺(tái)繪制Circos圖分析各片區(qū)門(mén)分類(lèi)水平下微生物群落結(jié)構(gòu)特征。

2結(jié)果與討論

2.1長(zhǎng)蕩湖水質(zhì)現(xiàn)狀

長(zhǎng)蕩湖理化因子濃度與空間分布有一定關(guān)聯(lián)。如表1所示，長(zhǎng)蕩湖3個(gè)理化因子呈現(xiàn)從南到北濃度逐步降低的趨勢(shì)。其中TN濃度從1.14 mg/L降至0.68 mg/L；TP濃度從0.1 mg/L降至0.03 mg/L；TOC含量從5.4 mg/L降至3.8 mg/L。這是因?yàn)樨S水期降雨量大，且人湖河流在長(zhǎng)蕩湖北部的輸水量更多，因此造成了長(zhǎng)蕩湖南部污染物濃度高于北部。長(zhǎng)蕩湖各片區(qū)的WT、pH差距不大，長(zhǎng)蕩湖南部片區(qū)和北部片區(qū)DO濃度低于中部和東部片區(qū)。整體而言，除南部片區(qū)，長(zhǎng)蕩湖其他片區(qū)的理化指標(biāo)均達(dá)到地表Ⅲ類(lèi)水標(biāo)準(zhǔn)。

運(yùn)用水污染指數(shù)法、內(nèi)梅羅指數(shù)法分別對(duì)豐水期長(zhǎng)蕩湖進(jìn)行水質(zhì)評(píng)價(jià)。2種水質(zhì)評(píng)價(jià)方法結(jié)果顯示：長(zhǎng)蕩湖南部片區(qū)比其他3個(gè)片區(qū)水質(zhì)現(xiàn)狀嚴(yán)重，超標(biāo)因子為T(mén)N、TP。但采用水污染指數(shù)法評(píng)價(jià)的水質(zhì)現(xiàn)狀普遍低于內(nèi)梅羅指數(shù)法。因?yàn)樗廴局笖?shù)法是基于單因子指數(shù)法，即采用“最差規(guī)則”來(lái)確定水質(zhì)等級(jí)，雖然對(duì)單項(xiàng)指標(biāo)有較好的評(píng)價(jià)意義，但缺乏對(duì)水體環(huán)境綜合考慮?？傮w而言，長(zhǎng)蕩湖南部片區(qū)的部分水質(zhì)指標(biāo)（TN、TP）超標(biāo)，但長(zhǎng)蕩湖整體水質(zhì)達(dá)標(biāo)。

2.2長(zhǎng)蕩湖微生物群落分布特征

2.2.1微生物群落Alpha多樣性指數(shù)分析

長(zhǎng)蕩湖微生物群落Alpha多樣性指數(shù)分析結(jié)果見(jiàn)表2。

長(zhǎng)蕩湖4個(gè)采樣點(diǎn)的Coverage指數(shù)均在97%以上，說(shuō)明此次測(cè)序結(jié)果可靠。Ace、Chaol指數(shù)代表微生物豐富度。北部片區(qū)Ace指數(shù)最高，達(dá)5 462；中部片區(qū)最低，為2 590。相同地，北部片區(qū)Chaol指數(shù)最高（5248），而中部片區(qū)最低（2528）。Shannon、Simpson指數(shù)表示微生物群落的多樣性。Shannon指數(shù)最高的是北部片區(qū)，其次為東部片區(qū)、南部片區(qū)、中部片區(qū)。Simpson指數(shù)與Shannon指數(shù)反映結(jié)論一致?？傮w而言，長(zhǎng)蕩湖微生物群落的豐富度與多樣性呈現(xiàn)出北部片區(qū)gt;東部片區(qū)gt;南部片區(qū)gt;中部片區(qū)的規(guī)律。

2.2.2門(mén)分類(lèi)水平微生物群落結(jié)構(gòu)

長(zhǎng)蕩湖微生物群落在門(mén)分類(lèi)水平上具有較高的多樣性。4個(gè)采樣點(diǎn)中共檢測(cè)出32種已知微生物菌門(mén)，長(zhǎng)蕩湖微生物在門(mén)分類(lèi)水平上組成如圖2（a）所示，其中相對(duì)豐度排在10名之后的菌門(mén)和其他未知物種歸為其它（ others）。長(zhǎng)蕩湖相對(duì)豐度大于10_10的主要菌門(mén)有放線菌門(mén)（Actinobacteria）、變形菌門(mén)（Proteobacteria）、厚壁菌門(mén)（Firmicutes）、擬桿菌門(mén)（Bacteroidetes）、疣微菌門(mén)（Verrucomicrobia）、藍(lán)細(xì)菌門(mén)（Cyanobacteria）、綠彎菌門(mén)（Chloroflexi）、酸桿菌門(mén)（Acidobacteriota）。

不同片區(qū)優(yōu)勢(shì)菌門(mén)的種類(lèi)相似。如放線菌門(mén)（Actinobacteria）為長(zhǎng)蕩湖各片區(qū)占比最高的優(yōu)勢(shì)菌門(mén)，據(jù)報(bào)道放線菌門(mén)（ Actinobacteria）是濕地、湖泊及海洋等各種水環(huán)境細(xì)菌群落中的優(yōu)勢(shì)門(mén)類(lèi)。相比其他細(xì)菌，放線菌能更好地躲避原生動(dòng)物的捕食，死亡率相對(duì)較低；在氮、磷等營(yíng)養(yǎng)物質(zhì)元素充足的湖泊中，碳源常會(huì)限制細(xì)菌生長(zhǎng)繁殖，而放線菌在利用碳源方面具有更強(qiáng)的競(jìng)爭(zhēng)力，所以水中放線菌門(mén)（Actinobacteria）的占比相對(duì)較高。但不同片區(qū)優(yōu)勢(shì)菌門(mén)的相對(duì)豐度卻有所差別。藍(lán)細(xì)菌門(mén)（Cyanobac-teria）在長(zhǎng)蕩湖北部片區(qū)和東部片區(qū)的相對(duì)豐度分別為13.27%和5.86%，高于長(zhǎng)蕩湖中部片區(qū)（2.93%）和南部片區(qū)（0.81%）。研究表明藍(lán)細(xì)菌門(mén)（Cyanobacteria）繁殖過(guò)多則會(huì)引發(fā)水華、使水體富營(yíng)養(yǎng)化且藍(lán)藻生長(zhǎng)潛力由營(yíng)養(yǎng)物質(zhì)濃度決定。但此研究中北部、東部片區(qū)的氮磷濃度低于南部和中部片區(qū)，猜測(cè)是人湖河流含有大量的藍(lán)細(xì)菌，進(jìn)而導(dǎo)致北部、東部片區(qū)藍(lán)細(xì)菌門(mén)（Cyanobacteria）的豐度急劇上升。長(zhǎng)蕩湖南部片區(qū)和中部片區(qū)厚壁菌門(mén)（Firmicutes）的相對(duì)豐度為16.86%、20.92%，遠(yuǎn)高于東部（0.25%）和北部片區(qū)（0.35%）。這與李建松等研究結(jié)論一致。研究發(fā)現(xiàn)厚壁菌門(mén)（Firmi-cutes）是藍(lán)藻水華降解末期的最優(yōu)勢(shì)菌門(mén)。原因可能是藍(lán)藻降解末期水體中溶解氧幾乎耗盡，引起了嗜陽(yáng)菌屬（Heliophilum）的大量增殖，進(jìn)而導(dǎo)致厚壁菌門(mén)（Firmicutes）在水體中的比例上升。

Circos圖可以直觀展現(xiàn)長(zhǎng)蕩湖不同片區(qū)的門(mén)分類(lèi)水平微生物群落組成情況。如圖2（b）所示，放線菌門(mén)（Actinobacteria）和變形菌門(mén)（Proteobacteria）在不同片區(qū)生物量分布均勻。如放線菌門(mén)（Actinobac-teria）是長(zhǎng)蕩湖各片區(qū)占比最高的菌門(mén)，其在南部片區(qū)、中部片區(qū)、東部片區(qū)及北部片區(qū)的相對(duì)豐度依次為36.6%、43.01%、39.56%、33.85%。疣微菌門(mén)（Verrucomicrobia）、藍(lán)細(xì)菌門(mén)（Cyanobacteria）相對(duì)豐度呈現(xiàn)由南到北逐漸遞增的趨勢(shì)。而厚壁菌門(mén)（Firmicutes）相對(duì)豐度呈現(xiàn)由南到北逐漸遞減的趨勢(shì)。特別是擬桿菌門(mén)（Bacteroidetes）在北部片區(qū)（13.72%）、南部片區(qū)（11.52%）、東部片區(qū)（11.85%）的相對(duì)豐度相近，但高于中部片區(qū)（3.94%）。由于人湖河流攜帶的污染物對(duì)長(zhǎng)蕩湖水質(zhì)有一定影響，猜測(cè)長(zhǎng)蕩湖北部、南部和東部片區(qū)的擬桿菌門(mén)相對(duì)豐度的提升來(lái)源于人湖河流的輸入?？梢?jiàn)長(zhǎng)蕩湖微生物群落結(jié)構(gòu)與空間位置、外源輸入有關(guān)。

2.2.3科分類(lèi)水平微生物群落結(jié)構(gòu)

豐水期長(zhǎng)蕩湖共檢測(cè)出378種微生物菌科，如圖3所示，相對(duì)豐度前5的菌科依次為：螢光桿菌科（Ilumatobacteraceae）、魚(yú)孢菌科（Sporichthyaceae）、動(dòng)球菌科（Planococcaceae）、叢毛單胞菌科（Co-mamonadaceae）。莫拉菌科（Moraxellaceae）。

不同片區(qū)的科分類(lèi)水平下微生物相對(duì)豐度有一定差別。如動(dòng)球菌科（Planococcaceae）在長(zhǎng)蕩湖南部片區(qū)和中部片區(qū)的相對(duì)豐度遠(yuǎn)高于東部和北部片區(qū)，據(jù)報(bào)道動(dòng)球菌科（Planococcaceae）的生長(zhǎng)受多種環(huán)境因素共同影響，如pH、鹽度等。推測(cè)可能是人湖河流的輸入導(dǎo)致東部和北部片區(qū)水體理化性質(zhì)發(fā)生變化，從而造成動(dòng)球菌科（Planococcaceae）含量降低。長(zhǎng)蕩湖南部片區(qū)和中部片區(qū)莫拉菌科（Moraxellaceae）的含量同樣高于東部和北部片區(qū)。說(shuō)明科分類(lèi)水平下微生物與外源輸入有關(guān)。

2.2.4屬分類(lèi)水平微生物群落結(jié)構(gòu)

長(zhǎng)蕩湖4個(gè)采樣點(diǎn)共檢測(cè)出467種微生物菌屬；但所有采樣點(diǎn)均有未能確定分類(lèi)學(xué)地位的菌屬。如圖4所示，相對(duì)豐度排在10名之后的菌屬和其他未知物種歸為其它（others）。長(zhǎng)蕩湖相對(duì)豐度前5的菌屬依次為：CL500-29 marine group、hgcl clade、類(lèi)芽孢八疊球菌屬（Paenisporosarcina）、不動(dòng)桿菌屬（Acinetobacter）、Clade Ill_norank。

長(zhǎng)蕩湖不同片區(qū)優(yōu)勢(shì)菌屬種類(lèi)相似，但占比不同。CL500-29 marine group是長(zhǎng)蕩湖占比最高的菌屬，其在中部片區(qū)的占比高于其他片區(qū)。hgclclade分布情況與CL500-29 marine group一致。類(lèi)芽孢八疊球菌屬（Paenisporosarcina）、不動(dòng)桿菌屬（Acinetobacter）的相對(duì)豐度呈現(xiàn)出南部片區(qū)、中部片區(qū)高于東部片區(qū)和北部片區(qū)。研究表明不動(dòng)桿菌屬（ Acinetobacter）屬于變形菌門(mén)，是一類(lèi)重要的反硝化聚磷菌，其生長(zhǎng)繁殖受營(yíng)養(yǎng)物質(zhì)影響；由于長(zhǎng)蕩湖氮磷濃度自南到北逐步降低，因此不動(dòng)桿菌屬（Acinetobacter）的相對(duì)豐度也是自南到北逐漸下降。相反地，Clade Ill_norank在東部片區(qū)、北部片區(qū)的相對(duì)豐度略高于南部片區(qū)與中部片區(qū)?？赡苁侨牒恿鬏斎朐斐傻?。

3結(jié)論

兩種水質(zhì)評(píng)價(jià)方法結(jié)果顯示，豐水期長(zhǎng)蕩湖南部片區(qū)污染最為嚴(yán)重，其中超標(biāo)因子為T(mén)N、TP，但整體上長(zhǎng)蕩湖水質(zhì)達(dá)標(biāo)。高通量測(cè)序結(jié)果表明，長(zhǎng)蕩湖4個(gè)采樣點(diǎn)共檢測(cè)出微生物群落32門(mén)378科467屬。長(zhǎng)蕩湖各片區(qū)的微生物群落中優(yōu)勢(shì)菌門(mén)、菌科、菌屬種類(lèi)相似，但相對(duì)豐度卻有所差別。相對(duì)豐度前3的優(yōu)勢(shì)菌門(mén)、菌科分別為放線菌門(mén)（Actinobacte-ria）、變形菌門(mén)（Proteobacteria）、厚壁菌門(mén)（Firmi-cutes）與螢光桿菌科（Ilumatobacteraceae）、魚(yú)孢菌科（Sporichthyaceae）、動(dòng)球菌科（Planococcaceae）；占比前3的優(yōu)勢(shì)菌屬包含CL500-29 marine group、hgclclade、類(lèi)芽孢八疊球菌屬（Paenisporosarcina）。微生物群落Alpha多樣性分析結(jié)果顯示，長(zhǎng)蕩湖微生物群落多樣性和豐富度呈現(xiàn)北部片區(qū)gt;東部片區(qū)gt;南部片區(qū)gt;中部片區(qū)的規(guī)律。長(zhǎng)蕩湖微生物群落結(jié)構(gòu)與空間位置、外源輸入有關(guān)。