杜蕾, 李暢游, 李文寶, 史小紅, 楊旭, 劉晶晶

夏季達(dá)里諾爾湖浮游細(xì)菌群落表、底層結(jié)構(gòu)特征及其關(guān)鍵驅(qū)動因子

杜蕾, 李暢游*, 李文寶, 史小紅, 楊旭, 劉晶晶

內(nèi)蒙古農(nóng)業(yè)大學(xué)水資源保護(hù)與利用自治區(qū)重點實驗室, 呼和浩特 010018

浮游細(xì)菌是湖泊水體物質(zhì)循環(huán)和能量流動過程中的關(guān)鍵要素, 影響著整個生態(tài)系統(tǒng)的遷移轉(zhuǎn)化。為探索中國北方典型內(nèi)陸封閉型湖泊浮游細(xì)菌群落結(jié)構(gòu)對環(huán)境因子變動的響應(yīng), 以夏季內(nèi)蒙古達(dá)里諾爾湖(以下簡稱“達(dá)里湖”)為例, 采集表、底層水共32個樣品, 基于16S rRNA基因測序, 詳細(xì)分析了浮游細(xì)菌群落結(jié)構(gòu)的表、底層變化特征及其主要影響因素, 結(jié)果顯示: (1)夏季達(dá)里湖水體中浮游細(xì)菌群落多樣性和均勻度均呈現(xiàn)表層＞底層, 豐富度則呈底層＞表層。(2)浮游細(xì)菌群落以變形菌門(Proteobacteria)、放線菌門(Actinobacteria)、藍(lán)細(xì)菌門(Cyanobacteria)和擬桿菌門(Bacteroidetes)四大優(yōu)勢菌門為主要組成部分, 且變形菌門含量在表、底層水體中均占絕對優(yōu)勢: 在底層水中占比46.46 %, 在表層水中占比36.21 %。(3)基于不同環(huán)境因子和浮游細(xì)菌優(yōu)勢菌門間的Pearson相關(guān)性分析, 表層水體中總?cè)芙庑怨腆w(TDS)、溶解態(tài)無機磷(DIP)是多數(shù)優(yōu)勢菌門的主要驅(qū)動因子, 而底層水體中的總氮(TN)、葉綠素a(Chl)和水溫(WT)為主要驅(qū)動因子?？偟膩碚f, 達(dá)里諾爾湖表、底層浮游細(xì)菌群落結(jié)構(gòu)與不同環(huán)境因子間存在顯著相關(guān), 而水體環(huán)境因子的表、底層差異是影響多數(shù)優(yōu)勢菌門演化的重要因素。

浮游細(xì)菌; 表、底層變化; 夏季; 群落結(jié)構(gòu); 達(dá)里諾爾湖

0 前言

湖泊在區(qū)域生態(tài)系統(tǒng)中起著至關(guān)重要的作用, 而湖泊水體中的微生物則是這一生態(tài)系統(tǒng)中的關(guān)鍵組成部分[1-2]。浮游細(xì)菌作為主要分解者[3-4], 參與了湖泊水體中的各種物質(zhì)循環(huán)和能量流動過程[5-8]。此外, 水體中的浮游細(xì)菌群落結(jié)構(gòu)變化與其水質(zhì)也有著密切的聯(lián)系, 且對水體環(huán)境變化響應(yīng)快速: 不同形態(tài)營養(yǎng)元素的輸入、水體環(huán)境因子的變化都會導(dǎo)致細(xì)菌群落結(jié)構(gòu)的改變[9-10]。因此, 分析水體浮游細(xì)菌群落結(jié)構(gòu)及優(yōu)勢菌門含量變化可以很好地反映湖泊的水環(huán)境質(zhì)量狀況[11-14]。

目前, 關(guān)于湖泊水體生態(tài)系統(tǒng)中浮游細(xì)菌菌群結(jié)構(gòu)特征及其與環(huán)境因子之間關(guān)系的研究已經(jīng)廣泛展開。如Garcia, De Figueiredo等通過研究季節(jié)變化對湖泊水體環(huán)境的影響, 進(jìn)而討論了其中細(xì)菌群落結(jié)構(gòu)的變化[1,8]; 在國內(nèi), 姜學(xué)霞、呂明姬、沈烽、楊長明等分別對河、湖水體及沉積物中的浮游細(xì)菌群落結(jié)構(gòu)特征及影響因素進(jìn)行了分析[14-17]。結(jié)果顯示不同環(huán)境因子, 如水體溫度(WT)、溶解氧含量(DO)、溶解性總固體含量(TDS)、營養(yǎng)元素(如總磷-TP、總氮-TN等)等, 對浮游細(xì)菌的群落結(jié)構(gòu)影響明顯[17,18]。然而, 相對于這些研究區(qū)域, 中國內(nèi)蒙古高原內(nèi)陸封閉型湖泊不僅水文循環(huán)過程影響因素簡單[19], 而且氣候環(huán)境變化存在明顯的冬季冰封期(相對封閉環(huán)境)-夏季非冰封期(相對開放環(huán)境)交替的區(qū)域特點[20, 21]。因此, 在具有重要區(qū)域生態(tài)環(huán)境意義的內(nèi)陸湖泊開展菌群結(jié)構(gòu)變化研究, 既可以為深入探討區(qū)域生態(tài)環(huán)境演變提供理論基礎(chǔ), 又可以為不同氣候條件下的湖泊水體生態(tài)系統(tǒng)演變對比提供基礎(chǔ)材料。

基于此, 本文將以內(nèi)蒙古赤峰市克什克騰旗境內(nèi)達(dá)里諾爾湖(簡稱“達(dá)里湖”)為研究區(qū)域, 通過對2018年夏季湖水進(jìn)行垂向采樣, 分析不同空間部分環(huán)境因子參數(shù)(如WT、DO、TDS、水深、不同形態(tài)營養(yǎng)元素含量等)的變化特征; 同時, 利用16S rRNA基因測序?qū)Σ煌顦悠分械母∮渭?xì)菌群落結(jié)構(gòu)進(jìn)行測試, 分析浮游細(xì)菌群落結(jié)構(gòu)組成及垂向變化特征; 進(jìn)一步結(jié)合Pearson相關(guān)性分析方法研究浮游細(xì)菌優(yōu)勢菌門和不同環(huán)境因子間的關(guān)系, 討論影響浮游細(xì)菌群落組成及演變的關(guān)鍵因素, 為進(jìn)一步研究達(dá)里湖水環(huán)境變化、水生態(tài)演化過程提供理論基礎(chǔ), 以期為達(dá)里諾爾湖水體環(huán)境保護(hù)提供一定參考。

1 研究區(qū)概況

達(dá)里湖(43o12′—43o24′ N, 116o24′—116o56′ E)位于內(nèi)蒙古自治區(qū)赤峰市克什克騰旗境內(nèi), 為封閉式蘇達(dá)型半咸水湖, 屬高原內(nèi)陸湖, 湖水無外瀉, 補給源為大氣降水, 河水(亮子河、沙里河、浩來河與貢格爾河)及地下水 (圖1), 其中亮子河與貢格爾河為主要補給河流, 達(dá)里湖年入湖水量存在差異, 主要由于年均降水量不同[19]。

達(dá)里湖水面面積平均約為189 km2, 年均面積變化在1—2 km2內(nèi)[19]。湖泊的面積變化造成不同區(qū)域水深出現(xiàn)一定程度波動, 具體體現(xiàn)在從西南到東北逐漸變淺, 湖水最深處約13 m[21]。其中, 西南側(cè)水深變化明顯, 形成了顯著的水下陡坡。而其他區(qū)域湖岸均比較平緩, 湖水深度變化也較平緩[22, 23]。

注: (a) 取樣點分布示意圖。取樣點依具體位置由北到南依次編號。左上小圖代表內(nèi)蒙古自治區(qū)的年均降水量區(qū)域變化特征; (b)達(dá)里湖水深變化概況示意圖, 修改自肖舉樂, 等[22, 23]. 圖中白色數(shù)值代表當(dāng)時水深。由于近年達(dá)里湖水位變化, 具體湖水深度存在一定程度的波動。

Figure 1 Distribution of sampling points and water depth profile of Dali-nor Lake

2 材料與方法

2.1 樣品采集

根據(jù)達(dá)里湖的水面面積、水深以及入湖河口的位置, 本文采用GPS全球定位系統(tǒng), 在達(dá)里湖內(nèi)共布設(shè)17個取樣點(圖1), 于2018年7月進(jìn)行了湖水樣品采集工作。樣品采集主要劃分為表層、底層2層, 其中DL-9水深較淺, 僅采集表層水樣, 由于實驗操作失誤丟失DL-5底層水樣, 共采集水樣32個。采集完成的水樣裝入經(jīng)滅菌處理的1 L的聚乙烯取樣瓶中, 并立即密封避光保存在低溫保溫箱中盡快運回實驗室進(jìn)行過濾處理, 其中500 mL用于環(huán)境因子的測定, 500 mL用于收集細(xì)菌, 首先水樣先經(jīng)過5 μm孔徑濾膜過濾, 以去除顆粒雜質(zhì)(濾膜可用以測試附著細(xì)菌)[24, 25], 然后在經(jīng)過0.22 μm孔徑濾膜過濾[24, 25], 過濾完成后, 濾膜置于凍融管中保存在-20℃條件下, 用于后續(xù)實驗分析。

2.2 水體環(huán)境因子測定

利用YSI Pro多參數(shù)水質(zhì)分析儀現(xiàn)場測定DO、WT、TDS、EC, 塞氏圓盤測定SD、測深鉛錘測定水深。同時, TN、TP、溶解態(tài)無機磷(DIP)、溶解態(tài)總磷(DTP)、葉綠素(Chl)等環(huán)境因子的測定則是在內(nèi)蒙古農(nóng)業(yè)大學(xué)水資源利用與保護(hù)自治區(qū)重點實驗室完成, 標(biāo)準(zhǔn)參照《水和廢水監(jiān)測分析方法》(第四版)[26]。

2.3 DNA提取和高通量測序

按照FastDNA? Spin Kit for Soil試劑盒的具體實驗操作步驟方法, 提取湖水中微生物基因組DNA。隨后利用微量紫外分光光度計(NanoDrop-2000)測定提取出的DNA濃度和純度后, 選用通用引物515F(5`- GTGCCAGCMGCC-GCGGTAA-3`)和806R (5`-GGACTACVSGGG TATCTAAT-3`)[27]對V4-V5 可變區(qū)進(jìn)行PCR 擴(kuò)增。

擴(kuò)增程序[28]為95℃預(yù)變性3 min, 27 個循環(huán)(95 ℃變性30 s, 55 ℃退火30 s, 72 ℃延伸30 s), 最后72 ℃延伸10 min (PCR 儀: ABI GeneAmp? 9700 型)。擴(kuò)增體系為20 uL, 4 uL 5*FastPfu 緩沖液, 2 uL 2.5 mM dNTPs, 0.8 uL 引物(5 uM), 0.4 uL FastPfu 聚合酶; 10 ng DNA 模板。每個樣品3 個重復(fù), 將不同樣品的PCR 擴(kuò)增產(chǎn)物均一化至10 nmol·L-1后等體積混合, 利用上海美吉生物醫(yī)藥科技有限公司的MiSeq PE300測序儀完成序列測定。

2.4 數(shù)據(jù)分析

高通量測序所得數(shù)據(jù)保存在NCBI數(shù)據(jù)庫內(nèi), 其生物信息學(xué)分析采用 Qiime軟件進(jìn)行, 根據(jù)序列的相似度, 將序列歸為多個OTU, OTU 產(chǎn)出后, 統(tǒng)計各個樣品中含有OTU的情況及每個OTU中所含序列的數(shù)目, 得到每個OTU 的分類學(xué)信息[29]。選取相似度在97%條件下的OTU生成預(yù)期的稀釋曲線, 并利用Mothur軟件計算豐富度指數(shù)Chao和Ace, 文庫覆蓋度, 均勻度指數(shù)Heip以及多樣性指數(shù)Simpson 和Shannon指數(shù)進(jìn)行Alpha 多樣性分析[30]。利用非度量多維尺度分析(NMDS)方法, 對不同深度湖水中浮游細(xì)菌群落結(jié)構(gòu)的垂向差異進(jìn)行分析。利用Pearson相關(guān)分析進(jìn)行了環(huán)境因子與優(yōu)勢菌群間的相關(guān)性分析。

3 結(jié)果與分析

3.1 達(dá)里湖水體環(huán)境因子表、底層差異變化

為探究達(dá)里湖水體環(huán)境因子表、底層差異, 對各取樣點的理化指標(biāo)計算均值, 并計算各指標(biāo)表、底均值的變異系數(shù)。Chl的變化最大, 變異系數(shù)為10.889%; pH值的變化最小, 變異系數(shù)為0.007%。除Chl外, 各理化指標(biāo)的變異系數(shù)均較小, 可見達(dá)里湖表、底層水體交換過程相對緩慢, 水體理化指標(biāo)變化相對穩(wěn)定[31]。

表1 表、底層水體環(huán)境因子的平均值及變異系數(shù)*

注: “*”平均值為每層水體所有取樣點的平均值, 黑體字代表平均值最高。

在夏季開放條件下的達(dá)里湖, 表層水更易受太陽輻射、大氣降水及河流輸入等外界因素的直接影響, 因此導(dǎo)致達(dá)里湖表層水的溫度、TDS值較高。夏季適宜的水體溫度導(dǎo)致藻類大量繁殖, 造成表層水體的Chl值高于底層, 而底層水體中TN、DIP、COD環(huán)境因子參數(shù)的平均值高于表層。

3.2 達(dá)里湖浮游細(xì)菌多樣性分析

表、底層水體樣品測序結(jié)果的文庫覆蓋率均高于99.7%(表2), 測序條數(shù)達(dá)到10000條以上時稀釋性曲線趨向平穩(wěn)(圖2), 表明測序數(shù)據(jù)量合理, 能夠代表物種的豐富度。采用群落豐富度指數(shù)Chao和Ace、群落多樣性指數(shù)Shannon 和Simpson、均勻度指數(shù)Heip和文庫覆蓋率等對達(dá)里湖浮游細(xì)菌群落進(jìn)行評估, 結(jié)果表明樣品具有豐富的群落組成, 其群落多樣性均較高(表2)。多樣性指數(shù)Shannon、Simpson、均勻度指數(shù)Heip均體現(xiàn)表層水體浮游細(xì)菌群落多樣性最高, 豐富度指數(shù)Chao 和Ace為底層水體最高, 文庫覆蓋率為表層高于底層。從標(biāo)準(zhǔn)差對比來看, 達(dá)里湖表、底層水體的浮游細(xì)菌群落豐富度較多樣性差異明顯。

3.3 浮游細(xì)菌群落結(jié)構(gòu)組成的表、底層分布特征

為具體探究達(dá)里湖表、底層水體中浮游細(xì)菌群落結(jié)構(gòu)的差異性特征, 采用非度量多維尺度分析(NMDS)方法, 繪制出樣點的二維分布圖(圖3), 其中不同樣點在二維圖中的距離可以反映群落結(jié)構(gòu)的相似程度, 樣點之間的距離越近則表示群落組成結(jié)構(gòu)越相似。

通過表、底層各個取樣點中細(xì)菌群落結(jié)構(gòu)的分析, 可以發(fā)現(xiàn)表層水體各取樣點的NMDS1值主要分布在-0.4—0.2之間, NMDS2值主要分布在-0.2—0.45之間; 底層水各取樣點的NMDS1值主要分布在-0.4—0.3間, NMDS2值主要分布在-0.15—0.35間(圖3)。顯然, 夏季開放條件下達(dá)里湖水體浮游細(xì)菌菌群結(jié)構(gòu)差異性程度呈現(xiàn)表層＞底層的特征。

基于表、底層水體中浮游細(xì)菌群落結(jié)構(gòu)的差異, 對表層和底層水中的浮游細(xì)菌在門類水平上的變化進(jìn)行了統(tǒng)計, 具體分析了達(dá)里湖表、底層水體中浮游細(xì)菌群類結(jié)構(gòu)特征及差異(圖4, 表3), 變形菌門(Proteobacteria)、放線菌門(Actinobacteria)、藍(lán)細(xì)菌門(Cyanobacteria) 3種浮游細(xì)菌菌門豐度在表、底層水體中均占優(yōu)勢比例(表3), 其中, 變形菌門的標(biāo)準(zhǔn)差遠(yuǎn)高于其他菌門: 變形菌門(Proteobacteria)的豐度最高, 在底層水(47.65%)中高于表層水(36.21%)中; 而放線菌門的豐度在底層水(24.89%)高于表層水(24.66%); 藍(lán)細(xì)菌門在表層水(11.22%)中高于底層水(8.99%)中。相對地, 擬桿菌門(Bacteroidetes)、軟壁菌門(Tenericutes)、異常球菌-棲熱菌門(Deinococcus-Thermus)、疣微菌門(Verrucomicrobia)、硬壁菌門 (Firmicutes)、浮霉?fàn)罹T(Plancto-mycetes)、酸桿菌門(Acidobacteria)和綠菌門(Chlorobi)等8種門類浮游細(xì)菌同樣在表層、底層湖水中的比例存在差異(表3)。

表2 浮游細(xì)菌表、底層分布多樣性分析

注: 黑體字代表平均值最高。

圖2 浮游細(xì)菌稀釋性曲線

Figure 2 Rarefaction curve of bacterioplankton community

圖3 達(dá)里湖表、底層水體浮游細(xì)菌群落結(jié)構(gòu)的NMDS分析

Figure 3 NMDS analysis of phytoplankton community structure in surface and bottom waters of Dali-nor Lake

表3 表、底層湖水門類水平上浮游細(xì)菌群落結(jié)構(gòu)比例(%)

圖4 表、底層水體門類水平浮游細(xì)菌群落結(jié)構(gòu)組成

Figure 4 Composition of phytoplankton community in surface and bottom water bodies on phylum level

3.4 環(huán)境因子與浮游細(xì)菌優(yōu)勢菌門之間的關(guān)系

浮游細(xì)菌作為湖泊水體中主要的微生物種群, 其群落結(jié)構(gòu)變化會對水體環(huán)境因子的垂向變化出現(xiàn)一定程度的響應(yīng)。因此, 利用Pearson相關(guān)分析對達(dá)里湖表、底水體中門類水平上的浮游細(xì)菌物種豐富度和環(huán)境因子變化間的相互關(guān)系進(jìn)行了詳細(xì)分析, 結(jié)果顯示表、底層水體中優(yōu)勢細(xì)菌門類受不同環(huán)境因子的影響程度也不相同(圖5)。

注: 橫向坐標(biāo)為水體環(huán)境因子, 縱向為物種信息, 中間不同顏色代表Pearson 相關(guān)分析的相關(guān)系數(shù)r值, r值范圍為-1—1. 當(dāng)r值小于零, 則顯示兩者負(fù)相關(guān), 當(dāng)r值大于零, 則顯示兩者正相關(guān). "*"代表顯著性檢驗P值<0.05有顯著性差異, "**"代表有極顯著差異。

Figure 5 Pearson analyzing result of planktonic bacterial flora and environmental factors in Dali-nor Lake

表層水體中, TDS、DIP作為與浮游細(xì)菌優(yōu)勢菌門群落結(jié)構(gòu)相關(guān)性最顯著的環(huán)境因子, 在和不同優(yōu)勢門類浮游細(xì)菌間的相關(guān)性結(jié)果上也有體現(xiàn)。異常球菌-棲熱菌門(Deinococcus- Thermus)、硬壁菌門(Firmicutes)與TDS呈極顯著正相關(guān), 而與DIP呈顯著負(fù)相關(guān); 變形菌門(Proteobacteria)與這TN、SD、TDS 3個環(huán)境因子呈顯著負(fù)相關(guān); pH與擬桿菌門(Bacteroidetes)、軟壁菌門(Tenericutes)呈顯著負(fù)相關(guān); 綠菌門(Chlorobi)與DIP呈顯著負(fù)相關(guān) (圖5a)。

底層水體中, 浮游細(xì)菌優(yōu)勢菌門群落結(jié)構(gòu)受TN、Chl和WT因子影響最顯著。TN與擬桿菌門(Bacteroidetes)、綠菌門(Chlorobi)、軟壁菌門(Tenericutes)呈顯著負(fù)相關(guān); DIP與軟壁菌門(Tenericutes)呈顯著正相關(guān); Chl與放線菌門(Actinobacteria)、擬桿菌門(Bacteroidetes)、綠菌門(Chlorobi)呈顯著正相關(guān); WT與酸桿菌門(Acidobacteria)、異常球菌-棲熱菌門(Deinococcus- Thermus)和硬壁菌門(Firmicutes)呈顯著正相關(guān); 水深與酸桿菌門(Acidobacteria)、異常球菌-棲熱菌門(Deinococcus-Thermus)呈顯著負(fù)相關(guān)(圖5b)。

4 討論

本文利用16S rRNA基因測序技術(shù)闡明了達(dá)里湖浮游細(xì)菌群落的空間變化, 多樣性指數(shù)表現(xiàn)為表層高, 底層低, 而豐富度指數(shù)為表層低, 底層高, 每個樣品在菌群組成上表、底層差異顯著, 這與其中變形菌門(Proteobacteria)、放線菌門(Actinobacteria)、藍(lán)細(xì)菌門(Cyanobacteria)是達(dá)里湖的主要菌門有關(guān), 這一分布特征也和大部分的淡水湖泊浮游細(xì)菌的組成相似[32,33]。

在達(dá)里湖夏季水體中, Proteobacteria作為優(yōu)勢菌中數(shù)量最高的菌門, 在達(dá)里湖表、底層水體中均占據(jù)明顯優(yōu)勢, 其物種多樣性豐富, 生存方式多樣[34], 這也導(dǎo)致在表、底層湖水中Proteobacteria 豐度變化與湖水理化指標(biāo)參數(shù)的相互關(guān)系復(fù)雜多變; 而相對豐度也相對較高的 Actinobacteria, 其豐度變化往往會受到水體富營養(yǎng)化等因素影響[35], 隨著水深的增加, Actinobacteria與營養(yǎng)元素的相關(guān)性產(chǎn)生差異; 其次是Cyanobacteria, 其生長方式為在表層湖水進(jìn)行光合作用, 消耗表層水中的營養(yǎng)物質(zhì), 當(dāng)營養(yǎng)物質(zhì)被消耗盡時, 他們便會隨細(xì)胞體下沉到底層水體或沉積物中繼續(xù)吸收營養(yǎng)物質(zhì), 進(jìn)而造成達(dá)里湖表層水體到底層水體與氮、磷營養(yǎng)鹽的相關(guān)關(guān)系發(fā)生改變。

達(dá)里湖屬典型內(nèi)陸封閉型湖泊, 常年湖水無外泄, 在夏季湖泊水體處于開放條件下, 表層水體更易受太陽輻射、大氣降水及河流輸入等外界因素的直接影響, 因此導(dǎo)致達(dá)里湖表層水的溫度、TDS值較高, 并且由于夏季適宜的水體溫度導(dǎo)致藻類大量繁殖, 造成表層水體的Chl值也較底層偏高, 進(jìn)而顯示出表層水體的理化性質(zhì)主要受外界因素輸入的影響; 而底層水體中TN、DIP、COD等環(huán)境因子參數(shù)高于表層, 這是由于在湖泊底部, 底層水體與底泥接觸, 而底泥中的氮、磷等污染物釋放造成污染物濃度偏高。隨著環(huán)境因子在表、底層水體的變化, 其浮游細(xì)菌群落結(jié)構(gòu)也相應(yīng)發(fā)生改變。總的來說, 湖泊水體中微生物群落的結(jié)構(gòu)特征與其環(huán)境理化性質(zhì)緊密相關(guān)并相互限制影響[36]。本研究發(fā)現(xiàn), 與達(dá)里湖浮游細(xì)菌呈顯著相關(guān)的環(huán)境因子分別為TN、DIP、Chl、PH、WT、SD、TDS、水深。TN、DIP及水深因子在表、底層均為主要環(huán)境因子。氮、磷元素循環(huán)是各種元素循環(huán)的中心, 微生物參與下的氮、磷循環(huán)是生態(tài)系統(tǒng)中的主要初級生產(chǎn)力, 在水生態(tài)環(huán)境中, 過量的氮、磷會導(dǎo)致水體發(fā)生富營養(yǎng)化。Haukka等[37]在對韋西湖微生物群落結(jié)構(gòu)的分析中發(fā)現(xiàn), 氮濃度與浮游細(xì)菌的組成直接相關(guān), 其中Verrucomicrobia能夠?qū)Φ谞I養(yǎng)鹽進(jìn)行有效的利用, 而Actinobacteria則適宜低濃度氮磷營養(yǎng)環(huán)境。在本研究中達(dá)里湖浮游細(xì)菌對氮磷營養(yǎng)鹽的響應(yīng)并不明顯, 氮濃度與多數(shù)浮游細(xì)菌的關(guān)系為不顯著相關(guān), 且DIP作為對達(dá)里湖部分浮游細(xì)菌有顯著影響的磷營養(yǎng)鹽, 對Chlorobi與Tenericutes有顯著促進(jìn)作用, 對Deinococcus-Thermus、Firmicutes呈抑制作用; 李玉華等[38]在對松花湖水體的研究中得出水深為主要環(huán)境因子, 這與本研究中水深變化影響水體環(huán)境因子改變進(jìn)而影響菌群結(jié)構(gòu)產(chǎn)生差異的結(jié)果相同, 顯示了松花湖與達(dá)里湖浮游細(xì)菌菌群結(jié)構(gòu)均受到水深變化的影響。夏季開放條件下, 達(dá)里湖表層懸浮大量藻類, 使得Chl含量偏高且與底層存在顯著差異, 從Pearson相關(guān)分析中Chl對多數(shù)優(yōu)勢菌門的生長呈顯著促進(jìn)作用, 除此之外, TDS也對表層的優(yōu)勢菌門具有顯著影響, 這與王歡等對碧流河水庫的研究存在一定差異[39]。顯然, 湖泊水體微生物結(jié)構(gòu)特征除受水體環(huán)境因子的調(diào)控外, 水域所處區(qū)域氣候條件及自身的區(qū)域特征也會對水體浮游細(xì)菌的菌群結(jié)構(gòu)及優(yōu)勢菌門的分布存在影響。

總之, 在開放條件下的達(dá)里湖, 表層和底層水體受外界環(huán)境變化影響的程度存在差異, 這進(jìn)一步導(dǎo)致了水體中環(huán)境因子的含量出現(xiàn)波動, 最終引起浮游細(xì)菌群落結(jié)構(gòu)及優(yōu)勢菌門的豐度出現(xiàn)變化。

5 結(jié)論

1) 夏季開放條件下, 湖水中變形菌門(Proteo-bacteria)、放線菌門(Actinobacteria)、藍(lán)細(xì)菌門(Cyanobacteria)和擬桿菌門(Bacteroidetes)為優(yōu)勢菌門, 且Proteobacteria含量在表、底層水體中中均占絕對優(yōu)勢: 底層水占比46.46 %, 表層水中占比36.21 %。表、底層水體中優(yōu)勢菌的占比存在一定差異, 這與水體環(huán)境差異有關(guān)。

2)表層水體中浮游細(xì)菌的群落多樣性、分布均勻程度最高。不過, 底層水體的整體細(xì)菌豐富度要高于表層水。

3) 水深變化引起的不同水體中環(huán)境因子含量的波動成為夏季達(dá)里湖浮游細(xì)菌菌群結(jié)構(gòu)和優(yōu)勢菌門豐度變化的關(guān)鍵因子, 即不同水深處浮游細(xì)菌群落結(jié)構(gòu)及優(yōu)勢菌門與不同水體環(huán)境因子間的相關(guān)性存在差異, 不同環(huán)境因子促進(jìn)或抑制了相應(yīng)優(yōu)勢菌門的發(fā)育。其中, 表層水體中TDS、DIP對浮游細(xì)菌優(yōu)勢菌門影響最顯著, 底層水體中TN、Chl和WT為主要影響因子。整體上, 水體環(huán)境因子的表、底層差異成為影響浮游細(xì)菌的優(yōu)勢菌門群落結(jié)構(gòu)的重要因素。

[1] Garcia S L, Salka I, Grossart H P, et al. Depth-discrete profiles of bacterial communities reveal pronounced spatio-temporal dynamics related to lake stratification[J]. Environmental Microbiology Reports, 2013, 5(4): 549–555.

[2] Li M, Wei G, Shi W,et al. Distinct distribution patterns of ammonia-oxidizing archaea and bacteria in sediment and water column of the Yellow River estuary[J]. Scientific Reports, 2018, 8(1): 1584.

[3] 陳禎, 何聃, 任麗娟. 溫度和營養(yǎng)鹽水平對淡水浮游細(xì)菌群落結(jié)構(gòu)的潛在影響[J].南京大學(xué)學(xué)報(自然科學(xué)), 2016, 52(4): 590–600.

[4] Newton R J, Jones S E, Eiler A, et alA guide to the natural history of freshwater lake bacteria[J]. Microbiology and Molecular Biology Reviews, 2011, 75(1): 14–49.

[5] Hanson C A, Fuhrman J A, Horner-Devine M C, et al. Beyond biogeographic patterns: processes shaping the microbial landscape[J]. Nature Reviews Microbiology, 2012, 10(7): 497–506.

[6] 閆冰,李俊生,陸晴,等. 北京城市公園綠地土壤微生物群落碳源代謝活性特征[J].環(huán)境科學(xué)研究, 2019, 32(9): 1567–1574.

[7] Yuan G, He J, Min C, et al. Factors dominating bacterioplankton abundance and production in the Nordic seas and the Chukchi Sea in summer 2012[J]. Acta Oceanologica Sinica, 2017, 36: 153–162.

[8] De Figueiredo D R, Pereira M J, Correia A. Seasonal modulation of bacterioplankton community at a temperate eutrophic shallow lake[J]. World Journal of Microbiology Biotechnology, 2010, 26(6): 1067–1077.

[9] Chen K, Allen J, Lu J R. Community Structures of Phytoplankton with Emphasison Toxic Cyanobacteria in an Ohio Inland Lake during Bloom Season[J]. Journal of Water Resource and Protection, 2017, 9(11): 1–29.

[10] Hahn M W. The microbial diversity of inland waters[J]. Current Opinion in Biotechnology, 2006, 17(3): 256-261.

[11] Sims A, Zhang Y, Gajaraj S, et al. Toward the development of microbial indicators for wetland assessment[J]. Water Research, 2013, 47(5): 1711–1725.

[12] Jiang J G, Shen Y F. Development of the microbial communities in Lake Donghu in relation to water quality[J]. Environmental Monitoring and Assessment, 2007, 127(1-3): 227–236.

[13] 程豹, 望雪, 徐雅倩, 等. 瀾滄江流域浮游細(xì)菌群落結(jié)構(gòu)特征及驅(qū)動因子分析[J]. 環(huán)境科學(xué), 2018, 39(8): 3649–3659.

[14] 姜學(xué)霞, 朱木得, 陳雪曼, 等. 肇慶波海湖細(xì)菌群落多樣性研究[J]. 環(huán)境科學(xué)導(dǎo)刊, 2019, 38(5): 20–23.

[15] 呂明姬, 汪杰, 范錚, 黃藝. 滇池浮游細(xì)菌群落組成的空間分布特征及其與環(huán)境因子的關(guān)系[J]. 環(huán)境科學(xué)學(xué)報, 2011, 31(2): 299–306.

[16] 沈烽, 趙大勇, 黃睿, 等. 南京玄武湖浮游細(xì)菌群落結(jié)構(gòu)的季節(jié)變化及其與環(huán)境因子的關(guān)系[J]. 湖泊科學(xué), 2017, 29(3): 662–669.

[17] 楊長明, 吳亞瓊, 王育來, 等. 南淝河表層沉積物細(xì)菌群落結(jié)構(gòu)特征及驅(qū)動因素[J]. 中國環(huán)境科學(xué), 2018, 38(9): 3552–3561.

[18] 李文寶, 劉志嬌, 楊旭, 等. 內(nèi)蒙古高原達(dá)里諾爾湖夏季水體穩(wěn)定同位素變化特征[J]. 湖泊科學(xué), 2019, 31(2): 539–550.

[19] 劉曉旭, 李暢游, 李文寶, 等. 冰封期達(dá)里諾爾湖同位素與營養(yǎng)鹽分布特征及關(guān)系的定量分析[J]. 湖泊科學(xué), 2015, 27(6): 1159–1167.

[20] 李衛(wèi)平. 高原典型湖泊營養(yǎng)元素地球化學(xué)循環(huán)與重金屬污染研究[D]. 呼和浩特: 內(nèi)蒙古農(nóng)業(yè)大學(xué), 2012.

[21] 王旭陽. 基于3S技術(shù)的達(dá)里諾爾湖水深反演研究[D]. 呼和浩特: 內(nèi)蒙古農(nóng)業(yè)大學(xué), 2017.

[22] 李文寶, 李暢游, 劉曉旭, 等. 達(dá)里諾爾湖水體穩(wěn)定氫、氧同位素組成變化對結(jié)冰過程的響應(yīng)[J]. 地球科學(xué)(中國地質(zhì)大學(xué)學(xué)報), 2015, 40(12): 2081–2090.

[23] 范佳偉, 肖舉樂, 溫銳林, 等. 內(nèi)蒙古達(dá)里湖全新世有機碳氮同位素記錄與環(huán)境演變[J]. 第四紀(jì)研究, 2015, 35(4): 856–870.

[24] Yang C Y, W Q, Simon P N, et al. Distinct Network Interactions in Particle-Associated and Free-Living Bacterial Communities during aBloom in a Plateau Lake[J]. Frontiers in Microbiology, 2017, 8: 1202–1227.

[25] Angelika R, Herlemann D P R. Particle-Associated Differ from Free-Living Bacteria in Surface Waters of the Baltic Sea[J]. Frontiers in Microbiology, 2015, 6: 1297–1310.

[26] 國家環(huán)境保護(hù)總局. 水和廢水監(jiān)測分析方法(第四版)[M]. 北京: 中國環(huán)境科學(xué)出版社, 2016, 35(1): 139–147.

[27] Bates S T, Berg-Lyons D, Caporaso J G, et al. Examining the global distribution of dominant archaeal populations in soil[J]. The International Society for Microbial Ecology Journal, 2011, 5(5): 908–917.

[28] Marti E, Balcázar, José Luis. Use of pyrosequencing to explore the benthic bacterial community structure in a river impacted by wastewater treatment plant discharges[J]. Research in Microbiology, 2014, 165(6): 468–471.

[29] Caporaso J G, Kuczynski J, Stombaugh J, et al. QIIME allows analysis of high-throughput community sequencing data [J]. Nature Methods, 2010, 7(5): 335–336.

[30] Schloss P D, Westcott S L, Ryabin T, et al. Introducing mothur: open-source, platform-independent, community- supported software for describing and comparing microbial communities[J]. Applied and Environmental Microbiology, 2009, 75(23): 7537–7541.

[31] 劉志嬌. 達(dá)里諾爾湖水動力條件及氫氧穩(wěn)定同位素試驗研究[D]. 呼和浩特: 內(nèi)蒙古農(nóng)業(yè)大學(xué), 2015.

[32] 錢瑋, 范竟成, 王金虎, 等. 東太湖地區(qū)濕地水質(zhì)與浮游細(xì)菌群落結(jié)構(gòu)關(guān)聯(lián)性研究[J]. 生物技術(shù), 2017, 1: 90–96.

[33] 劉晨宇, 祝亞楠, 王繼華, 等. 不同季節(jié)尾水排放對河道細(xì)菌群落結(jié)構(gòu)的影響[J]. 環(huán)境科學(xué)研究, 2018, 31(5): 893–900.

[34] 丁軼睿, 李定龍, 張毅敏, 等. 滆湖底泥細(xì)菌群落結(jié)構(gòu)及多樣性[J]. 環(huán)境科學(xué)學(xué)報, 2017, 37(5): 1649–1656.

[35] Allgaier M, Grossart H P. Diversity and Seasonal Dynamics of Actinobacteria Populations in Four Lakes in Northeastern Germany[J]. Applied and Environmental Microbiology, 2006, 72(5): 3489–3497

[36] 萬甜, 何夢夏, 任杰輝, 等. 渭河流域水體細(xì)菌群落的環(huán)境響應(yīng)及生態(tài)功能預(yù)測[J]. 環(huán)境科學(xué), 2019, 40(08): 3588–3595.

[37]Haukka K, Kolmonen E, Hyder R, et al. Effect of Nutrient Loading on Bacterioplankton Community Composition in Lake Mesocosms[J]. Microbial Ecology, 2006, 51(2): 137– 146.

[38] 李玉華, 許其功, 趙越, 等. 松花湖水體中不同空間分布的細(xì)菌群落結(jié)構(gòu)分析[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2013, 32(4): 764–770.

[39] 王歡, 趙文, 謝在剛, 等. 碧流河水庫細(xì)菌群落結(jié)構(gòu)特征及其關(guān)鍵驅(qū)動因子[J]. 環(huán)境科學(xué), 2018, 39(8): 3660– 3669.

Surface and bottom characteristics of bacterioplankton community in summer Dali-nor Lake and its key driving factors

Du Lei, Li Changyou*, Li Wenbao, Shi Xiaohong, Yang Xu, Liu Jingjing

IMAR Key Laboratory of Water Resources Protecting and Utilization, Inner Mongolia Agricultural University, Hohhot 010018, China

Bacterioplankton is the key element in the material cycle and energyflow in lake, affecting the migration and transformation of the ecosystem. In this paper, we collect 32 water samples from the surface and bottom water in summer Dali-nor Lake, which locates in Inner Mongolia Plateau and is an inland closed lake. Based on the T-RFLP technology and 16S rRNA sequencing,the community structure of bacterioplankton and its vertical variation characteristics were analyzed in detail. In summer, the diversity and evenness of the bacterioplankton in surface water are higher than bottom water. Meanwhile, bottom water has the highest species richness.On the whole, the Proteobacteria, Actinobacteria, Cyanobacteria, and Bacteroidetes are the dominant phyla in summer Dali-nor Lake. The Proteobacteria is the most dominant phylum in all water layers, with the content of about 46.46% in bottom water, and about 36.21% in surface water, respectively. Furthermore, TDS, DIP in surface water are correlated with most dominant bacteria, TN, Chland WT in the bottom water are the main environmental factors influencing community distribution. In summary, there is a significant correlation between the surface and bottom planktonic bacterial community structure of Dali-norl Lake and different environmental factors, while the water environmental factors in surface and bottom waters have become important factors affecting the evolution of most dominant bacteria.

bacterioplankton; change of surface and bottom water layers; summer; community structure characteristics; Dali-nor Lake

10.14108/j.cnki.1008-8873.2021.06.002

X172

A

1008-8873(2021)06-013-08

杜蕾, 李暢游, 李文寶, 等. 夏季達(dá)里諾爾湖浮游細(xì)菌群落表、底層結(jié)構(gòu)特征及其關(guān)鍵驅(qū)動因子[J]. 生態(tài)科學(xué), 2021, 40(6): 13–20.

Du Lei, Li Changyou, Li Wenbao, et al. Surface and bottom characteristics of bacterioplankton community in summer Dali-nor Lake and its key driving factors[J]. Ecological Science, 2021, 40(6): 13–20.

2020-05-11;

2020-06-20

國家自然科學(xué)基金地區(qū)項目(52160021, 51869020); 國家重點研發(fā)計劃—國際合作專項(2017YFE0114800); 內(nèi)蒙古自治區(qū)自然科學(xué)基金項目(2019MS05032); 內(nèi)蒙古自治區(qū)科技攻關(guān)項目(2020GG0009); 內(nèi)蒙古自治區(qū)高等學(xué)?！扒嗄昕萍加⒉胖С钟媱潯表椖?NJYT-20-A14)

杜蕾(1996—), 女, 內(nèi)蒙古赤峰市人, 碩士研究生, 主要研究方向為湖泊水環(huán)境演變及修復(fù), E-mail: 18847163393@163.com

通信作者:李暢游, E-mail: nndlichangyou@163.com