張弘杰,徐慧敏,過梓栩,何 斐,曾 巾,趙大勇

1 河海大學(xué)水文水資源與水利工程科學(xué)國家重點實驗室,南京 210098 2 生態(tài)環(huán)境部南京環(huán)境科學(xué)研究所,南京 210042 3 湖泊與環(huán)境國家重點實驗室,中國科學(xué)院南京地理與湖泊研究所,南京 210008

水產(chǎn)養(yǎng)殖業(yè)是農(nóng)業(yè)的重要組成部分,也是全球最大的食品產(chǎn)業(yè)之一[1]。近年來,克氏原螯蝦(Procambarusclarkii)、中華鱉(Pelodiscussinensis)等重要淡水物種的人工養(yǎng)殖發(fā)展迅速[2]?？耸显r隸屬于十足目(Decapoda)、螯蝦科(Astacidae)[3],是我國養(yǎng)殖最廣的淡水螯蝦品種[4]。中華鱉屬龜鱉目(Testudinata)、鱉科(Trionychidae)[5],在我國分布廣泛,藥用價值極高,是我國重要的淡水養(yǎng)殖物種之一[6]。

細菌是生態(tài)系統(tǒng)中的重要組成部分[7]。在水產(chǎn)養(yǎng)殖系統(tǒng)中,細菌既能作為分解者降解養(yǎng)殖動物無法利用的有機質(zhì),又可以作為養(yǎng)殖動物的直接或間接餌料[8]。此外,細菌群落對環(huán)境的變化十分敏感,其與周圍的環(huán)境因子相互作用、往復(fù)調(diào)控[9]。在水產(chǎn)養(yǎng)殖的水體中,細菌群落的結(jié)構(gòu)、多樣性等均隨水環(huán)境的變化而改變,水體細菌群落的變化又會對水環(huán)境造成影響[10],因此,對養(yǎng)殖環(huán)境中的水體細菌群落進行深入研究具有重要意義[11]。除養(yǎng)殖水體外,沉積物也是細菌附著的主要介質(zhì)[12]。大量的細菌富集于沉積物中,通過同化、異化等代謝過程來影響沉積物中營養(yǎng)鹽的分布、轉(zhuǎn)化和利用等[13]。因此,沉積物中的細菌群落可以作為水環(huán)境變化和演替的重要標志[14],探究水產(chǎn)養(yǎng)殖對沉積物中細菌群落的影響也具有重要意義。

近年來,快速發(fā)展的高通量測序技術(shù)由于其具有通量大、測試結(jié)果準確等優(yōu)點成為準確、便捷和全面地研究養(yǎng)殖系統(tǒng)中的細菌群落結(jié)構(gòu)和組成的研究手段之一[15]。胡常巨等人[16]利用高通量測序技術(shù)探究網(wǎng)箱養(yǎng)殖對水環(huán)境細菌群落的影響,表明養(yǎng)殖活動顯著改變了養(yǎng)殖水環(huán)境中浮游細菌群落結(jié)構(gòu)。陶玲[11]利用高通量測序技術(shù)研究發(fā)現(xiàn),稻田—池塘復(fù)合循環(huán)水養(yǎng)殖系統(tǒng)與常規(guī)池塘養(yǎng)殖系統(tǒng)中的細菌群落結(jié)構(gòu)存在差異。但目前的研究大多是對比水產(chǎn)養(yǎng)殖(或某種特殊的養(yǎng)殖系統(tǒng),如循環(huán)水養(yǎng)殖系統(tǒng))和天然水體細菌群落的差異[17—18],鮮有研究利用高通量測序技術(shù),對養(yǎng)殖不同物種的水環(huán)境細菌群落進行深入探究。而水產(chǎn)養(yǎng)殖不同物種對水環(huán)境中的細菌群落可能產(chǎn)生完全不同的影響[19]。不同物種的生活習(xí)性和管理方式等可能會對細菌群落產(chǎn)生不同的影響,而不同的細菌群落又可能會對養(yǎng)殖水環(huán)境或作物產(chǎn)生不同的影響[20]。因此,利用高通量測序技術(shù)對養(yǎng)殖不同物種的水體和沉積物中細菌群落進行深入研究,對探究水產(chǎn)養(yǎng)殖對水環(huán)境中細菌群落的影響有著重要意義?？耸显r和中華鱉是兩種在生理特性、生活習(xí)性及管理方式等方面具有明顯差異的淡水養(yǎng)殖物種?？耸显r是雜食性動物,生長發(fā)育較快,挖洞能力強,抗逆力很強,能生活在一些輕度污廢水中[21]。中華鱉以肉食為主,生長發(fā)育相對緩慢,繁殖力低,具有上岸曬背的習(xí)性[5,22]。因此,我們預(yù)期養(yǎng)殖這兩種作物會對水體和沉積物產(chǎn)生不同的影響。綜上,本研究以養(yǎng)殖不同物種(克氏原螯蝦和中華鱉)的水體和沉積物細菌群落為研究對象,以基于16S rRNA基因的高通量測序技術(shù)為研究方法,揭示了養(yǎng)殖不同物種的水體和沉積物間細菌群落α多樣性、群落組成和優(yōu)勢細菌的異同,以期為養(yǎng)殖環(huán)境的維護和改善提供參考。

1 材料與方法

1.1 樣品采集

本次研究的水體和沉積物樣品均采集于江蘇省省級精品漁業(yè)園(119°01′48″E,32°13′51″N)的多年輪作稻田(Oryzasativa：亞洲栽培稻)濕地。在實驗開始前,利用不銹鋼板將該濕地劃分為兩個區(qū)域,分別投放克氏原螯蝦(Procambarusclarkii)和中華鱉(Pelodiscussinensis),并根據(jù)養(yǎng)殖物種分別命名為PC和PS區(qū)。該濕地首次嘗試稻-蝦/鱉共作的模式。PC和PS區(qū)均于2019年5月上旬進行水稻播種；6月底水稻插秧；7月上旬水稻進入返青分蘗期,此時投放幼蝦(20—30 kg/hm2)和幼鱉(1500只/hm2)；11月底收割水稻(共產(chǎn)出水稻360 kg)；12月底—2020年1月初打撈克氏原螯蝦和中華鱉。為保證PC和PS區(qū)中的本底微生物多樣性和結(jié)構(gòu)相似,兩塊稻田濕地的面積(5328 m2)、種植密度(46%)、管理方式和灌溉用水等均保持一致。為探究養(yǎng)殖不同作物對水體和沉積物細菌群落的影響,于養(yǎng)殖周期結(jié)束后(2020年1月)在上述2個區(qū)域中,各選擇6個采樣點,分別利用不銹鋼采水器(HAD-WB-SS)和抓斗式采泥器(CN-100)采集水體(水面以下30—50 cm)和表層沉積物(0—3 cm)樣品。為避免水稻根系微生物對沉積物樣品的影響,在采集沉積物樣品時,我們盡量避開水稻生長的地方,并確保所采集的樣品中無水稻根系。收集到樣品后,利用保溫箱(4 ℃)冷藏,快速運回實驗室。利用YSI 6600(YSI 6600,Yellow Springs,美國)對水體的溫度(T)、pH和溶解氧濃度(DO)分別進行原位測定。

1.2 樣品預(yù)處理

樣品運回實驗室后,利用真空冷凍干燥機(LABCONCO,美國)凍干沉積物樣品,并在樣品凍干后進行研磨和過篩(35目土壤篩)。取部分水體樣品,用0.22 μm的硅膠微孔濾膜進行過濾,分別冷凍保存濾膜和濾液于超低溫冰箱中(-70 ℃)。

1.3 樣品理化指標測定

沉積物的總氮(TN)和總磷(TP)濃度也采用堿性過硫酸鉀消解紫外分光光度法(HJ 636—2012)進行測定[23]。使用重鉻酸鉀氧化-分光光度法(HJ 615—2011)[26]測定沉積物的有機碳濃度(OC)。

各區(qū)域水體和沉積物的理化指標測定結(jié)果如表1所示。

表1 2個養(yǎng)殖區(qū)域水體和沉積物的環(huán)境變量Table 1 Environmental variables of water and sediment samples in the two aquacultural zones

1.4 DNA提取及高通量測序

水體濾膜和沉積物樣品分別采用E.Z.N.A.?水體DNA提取試劑盒(Omega Biotek,Doraville,GA,美國)和PowerSoil? DNA提取試劑盒(MoBio Laboratories,Inc.,Carlsbad,CA,美國)進行水體和沉積物總DNA的提取[27]。使用PowerCleanTMDNA純化試劑盒(MoBio Laboratories,Carlsbad,CA,美國)對DNA樣品進行純化[28]。采用1%的瓊脂糖凝膠電泳驗證所提取的DNA樣品。

使用引物對515F(5′-GTGCCAGCMGCCGCGG-3′)和806R(5′-GGACTACHVGGGTWTCTAAT-3′)對16S rRNA基因的高可變區(qū)V4區(qū)進行擴增[29]。PCR擴增體系為20 μL(0.8 μL上下游引物、10 ng DNA模板、2 μL dNTP、0.4 μL FastPfu酶,及無菌蒸餾水),PCR擴增程序為95 °C預(yù)變性2 min,25次的擴增循環(huán)(95 °C變性30 s,55 °C退火30 s,72 °C延伸30 s),并最終在72 °C下延伸5 min[30]。PCR產(chǎn)物合并后通過2%的瓊脂糖凝膠電泳檢驗。PCR產(chǎn)物經(jīng)純化和定量后,由諾禾致源生物科技有限公司(北京)利用Illumina Miseq測序平臺進行測定。本研究所獲得的所有細菌基因序列均已上傳到National Center for Biotechnology Information(NCBI)數(shù)據(jù)庫中,序列登記號為PRJNA698288。

測序結(jié)果使用QIIME (v 1.9.1)進行處理[31],使用“trim”功能篩選出截斷讀數(shù)超過10 bp且配對成功的序列[32]。然后使用FLASH (v 1.0.0)拼接序列[33],再使用QIIME (v 1.9.1)中的“vsearch”功能查找并去除嵌合體[34],使用“pick_de_novo_otus.py”命令,根據(jù)3%的差異原則劃分操作分類單元(Operational Taxonomic Units,OTU)[35]。對照 SILVA 16S rRNA 基因數(shù)據(jù)庫(SILVA123_QIIME_release,2016)[36],匹配分類學(xué)信息。使用“filter_taxa_from_otu_table.py”命令去除非細菌序列和序列數(shù)不高于總序列數(shù)0.0005%的稀有OTU[37],使用“parallel_align_seqs_pynast.py”命令計算進化樹,最后使用“multiple_rarefactions_even_depth.py”命令對最終完善的文件進行重采樣,并導(dǎo)出最終的OTU表。

基于OTU表,在QIIME (v 1.9.1)中使用“alpha_diversity.py”計算Observed OTUs和Shannon指數(shù)來反映細菌群落的α多樣性。再利用QIIME (v 1.9.1)中的“beta_diversity.py”計算Weighted Unifrac距離矩陣來反應(yīng)細菌群落的β多樣性。

1.5 多元統(tǒng)計分析

本研究利用軟件GraphPad Prism (v 8)繪制箱型圖來反映細菌群落α多樣性的高低?；贠TU表,利用R(v 3.5.3)中的程序包“Vennerable”進行Venn圖的繪制。根據(jù)β多樣性矩陣,使用R(v 3.5.3)中的“vegan”程序包進行非度量多維尺度分析(Nonmetric multidimensional scaling,NMDS),來可視化不同樣品間群落結(jié)構(gòu)的差異和相似性。為探究細菌群落與環(huán)境因子的相關(guān)性關(guān)系,利用R(v 3.5.3)中的“vegan”程序包進行冗余分析(redundancy analysis,RDA)。將平均相對豐度大于0.5%的細菌門作為優(yōu)勢菌門,并利用軟件GraphPad Prism (v 8)繪制堆疊圖。篩選平均相對豐度大于1%的OTUs作為優(yōu)勢OTUs,使用軟件IBM SPSS Statistics (v 22.0)分析優(yōu)勢OTUs與環(huán)境因子的相關(guān)關(guān)系(Pearson相關(guān)),并利用R(v 3.5.3)中的“Heatmap”程序包分別繪制熱圖。采用獨立樣本T檢驗對水體樣品和沉積物樣品指標進行差異顯著性檢驗。

2 結(jié)果與分析2.1 2個養(yǎng)殖區(qū)域的水體和沉積物環(huán)境變量

2.2 水產(chǎn)養(yǎng)殖不同物種的水體和沉積物中細菌群落的α多樣性特征

2個區(qū)域細菌群落的α多樣性如圖1所示,對于水體和沉積物樣品,PS區(qū)的Observed OTUs和Shannon指數(shù)均顯著(P<0.05)高于PC區(qū)。即PS區(qū)中水體和沉積物細菌群落的α多樣性均顯著高于PC區(qū)。

圖1 2個養(yǎng)殖區(qū)域中水體和沉積物的α多樣性Fig.1 α-diversity of bacterial communities within water and sediment samples between the two zonesObserved OTUs:OTU數(shù),Shannon:Shannon指數(shù);PC:克氏原螯蝦,PS:中華鱉;圖中*表示顯著性差異(獨立樣本T檢驗),*:P<0.05,***:P<0.001

對于養(yǎng)殖水體,2個區(qū)域共有2639個相同的OTUs,PS區(qū)的特有OTU數(shù)為2588遠大于PC區(qū)的195(圖2)。對于沉積物,PC和PS區(qū)的OTU數(shù)基本一致,絕大多數(shù)OTUs為兩區(qū)域共有,但PS區(qū)的特有OTU數(shù)(408)仍明顯高于PC區(qū)(143)。由此可見,沉積物中的共有OTU數(shù)占總OTU數(shù)的比重大于水體。并且,PS區(qū)中水體和沉積物的特有OTU數(shù)均高于PC區(qū)。

對于水體,PC區(qū)的特有OTUs共分類為15個細菌門,PS區(qū)的特有OTUs共分類為16個細菌門。兩個區(qū)域的特有OTUs均是屬于擬桿菌門(Bacteroidetes)的最多(PC：33個,PS：280個),并且SR1僅在PS區(qū)中觀測到(圖2)。對于沉積物,PC區(qū)的特有OTUs共分類為13個細菌門,其中屬于擬桿菌門(Bacteroidetes)和β-變形菌綱(Betaproteobacteria)的特有OTUs均為19個明顯高于其他門(圖2)。PS區(qū)的特有OTUs共分類為14個細菌門,其中屬于擬桿菌門(Bacteroidetes)的特有OTUs最多(70個)。

圖2 2個養(yǎng)殖區(qū)域中水體和沉積物細菌群落共有和特有OTU數(shù)的Venn圖及2個養(yǎng)殖區(qū)域中水體和沉積物細菌群落特有OTUs分類學(xué)信息(門/亞門水平)Fig.2 Venn plots showing of the number of OTUs of bacterial communities shared and specialized between the two zones within water and sediment samples.Taxonomy information of OTUs (phyla/subphyla)of bacterial communities specialized between the two zones within water and sediment samplesPC:克氏原螯蝦,PS:中華鱉

2.3 水產(chǎn)養(yǎng)殖不同物種的水體和沉積物中細菌群落結(jié)構(gòu)

圖3 基于Weighted Unifrac距離矩陣的非度量尺度分析(NMDS)及細菌群落結(jié)構(gòu)與環(huán)境因子間的冗余分析(RDA)Fig.3 Nonmetric multidimensional scaling analysis (NMDS)based on Weighted Unifrac dissimilarity and redundancy analysis (RDA)between bacterial community structure and environmental variablesPC:克氏原螯蝦,PS:中華鱉;圖中*表示顯著相關(guān),*:P<0.05,**:P<0.01,***:P<0.001

2.4 水產(chǎn)養(yǎng)殖不同物種的水體和沉積物中細菌群落組成

本研究共獲得1487129條序列、7247個OTUs,它們屬于細菌的34門、114綱、258目、504科和955屬。2個區(qū)域的優(yōu)勢(相對豐度>0.5%)細菌門/亞門如圖4所示。水體樣品中共有11個優(yōu)勢菌門,這11個門的相對豐度和占到所有細菌的91.46%—98.54%。沉積物樣品中共有13個優(yōu)勢菌門,其相對豐度和占到總細菌的77.52%—84.98%。

圖4 2個養(yǎng)殖區(qū)域中水體和沉積物細菌群落的優(yōu)勢菌門/亞門(相對豐度>0.5%)相對豐度Fig.4 Relative abundance of dominant phyla/subphyla (relative abundance>0.5%)of bacterial communities within water and sediment samples between the two zonesPC:克氏原螯蝦,PS:中華鱉

將相對豐度>1%作為篩選優(yōu)勢OTUs的閾值,水體細菌群落中共篩選出15個優(yōu)勢OTUs,而沉積物細菌群落中僅篩選出7個優(yōu)勢OTUs(圖5)。水體細菌群落優(yōu)勢OTUs中,6個屬于β-變形菌綱(Betaproteobacteria)、5個屬于擬桿菌門(Bacteroidetes)、3個屬于放線菌門(Actinobacteria)和1個屬于α-變形菌綱(Alphaproteobacteria)。PS區(qū)中OTU3394(叢毛單胞菌科,Comamonadaceae)、OTU12347(腐螺旋菌科,Saprospiraceae)、OTU12385(叢毛單胞菌科,Comamonadaceae)、OTU15668(嗜甲基菌屬,Methylophilus)、OTU16336(叢毛單胞菌科,Comamonadaceae)、OTU20544(叢毛單胞菌科,Comamonadaceae)、OTU38268(黃桿菌屬,Flavobacterium)、OTU41298(放線菌目,Actinomycetales)、OTU42336(沉積桿菌屬,Sediminibacterium)、OTU52126(多核桿菌屬,Polynucleobacter)和OTU61118(黃桿菌屬,Flavobacterium)的相對豐度顯著低于PC區(qū)(P<0.01),OTU94966(腐螺旋菌科,Saprospiraceae)的相對豐度顯著高于PC區(qū)(P<0.001)。

圖5 2個養(yǎng)殖區(qū)域中水體和沉積物的優(yōu)勢OTUs(相對豐度>1%)相對豐度熱圖Fig.5 Heatmap showing the relative abundance of dominant OTUs (relative abundance>1%)of bacterial communities within water and sediment samples between the two zonesPC:克氏原螯蝦,PS:中華鱉;圖中*表示組間存在顯著性差異(獨立樣本T檢驗),**:P<0.01,***:P<0.001;最低分類學(xué)等級中,不同小寫字母表示不同分類水平,p:門,c:綱,o:目,f:科,g:屬

沉積物優(yōu)勢OTUs中,5個屬于擬桿菌門(Bacteroidetes)、1個屬于β-變形菌綱(Betaproteobacteria)、1個屬于綠彎菌門(Chloroflexi)。這7個優(yōu)勢OTUs的相對豐度在兩個區(qū)域之間均存在顯著差異(P<0.001)。PS區(qū)中OTU11213(Dechloromonas屬)和OTU59770(Prolixibacteraceae科)的相對豐度均顯著高于PC區(qū)(P<0.001)。而OTU1482(擬桿菌門,Bacteroidetes)、OTU89828(擬桿菌門,Bacteroidetes)、OTU97626(擬桿菌門,Bacteroidetes)、OTU103576(擬桿菌門,Bacteroidetes)和OTU103578(厭氧繩菌科,Anaerolineaceae)的相對豐度則是PC區(qū)顯著更高(P<0.001)。

圖6 2個養(yǎng)殖區(qū)域中水體和沉積物的優(yōu)勢OTU(相對豐度>1%)與環(huán)境因子的相關(guān)性熱圖Fig.6 Heatmap showing correlation between the relative abundance of the dominant OTUs (relative abundance>1%)and environmental factors in water and sediment samples between the two zones*表示顯著相關(guān),*:P<0.05,**:P<0.01

對于沉積物,TN、TP和OC濃度均與優(yōu)勢OTUs有較強的相關(guān)性。整體上,TP與優(yōu)勢OTUs的相關(guān)性較TN和OC弱。OTU11213和OTU59770與TN和OC均顯著負相關(guān)(P<0.05),OTU1482、OTU89828、OTU97626、OTU103576和OTU103578與TN和OC均顯著正相關(guān)(P<0.01)。

3 討論

3.1 水產(chǎn)養(yǎng)殖不同物種對水體和沉積物中細菌群落多樣性及結(jié)構(gòu)的影響

本研究中,PS區(qū)的水體和沉積物細菌群落α多樣性均顯著高于PC區(qū)(P<0.05,圖1)。細菌群落多樣性受宿主影響較大[38],如張瓊瓊等人研究發(fā)現(xiàn),微生物多樣性與宿主呈現(xiàn)專一性,即宿主會在周圍環(huán)境中形成獨特的微環(huán)境并對細菌具有趨化和富集的作用,進而使得不同宿主對應(yīng)了不同大小的細菌群落多樣性[39]。因此,養(yǎng)殖克氏原螯蝦(PC)和中華鱉(PS)的水體和沉積物中極有可能分別形成了獨特的微環(huán)境,如2種不同物種的排泄物、皮膚分泌物等對細菌分別具有獨特的趨化和富集作用,導(dǎo)致PC和PS區(qū)的水體和沉積物間細菌群落種類和數(shù)量發(fā)生變化,進而使細菌群落多樣性與宿主呈現(xiàn)專一性,最終導(dǎo)致了PC和PS區(qū)細菌群落的α多樣性呈現(xiàn)顯著差異。Venn圖的結(jié)果顯示,細菌群落α多樣性更高的PS區(qū)中擁有更多的特有OTUs(圖2),這些OTUs可能就是在養(yǎng)殖中華鱉形成的特定微環(huán)境下定植并富集的。中華鱉還具有上岸曬背的習(xí)性,這一習(xí)性也可能將陸地中的細菌帶入到養(yǎng)殖水環(huán)境中,使特有OTUs數(shù)量增加,進而導(dǎo)致細菌群落多樣性升高[5,22]。本研究統(tǒng)計了PC和PS區(qū)中特有OTUs的分類學(xué)信息(圖2),但由于目前認知所限,這些特有OTUs大多只能獲得較為粗糙的分類學(xué)信息,因此這些OTUs的特性大部分還不清楚,需要在后續(xù)的研究中進一步探究,以探明養(yǎng)殖作物的哪些生理特征或生活習(xí)性會導(dǎo)致哪些特有OTUs的定殖和富集。此外,環(huán)境中的營養(yǎng)濃度對細菌群落多樣性也有著較大影響[40]。研究表明,在營養(yǎng)濃度較高的養(yǎng)殖環(huán)境中,更低的營養(yǎng)負荷反而可能會造成更高的多樣性[41]。而本研究中細菌群落α多樣性更高的PS區(qū)中,水體TN、TP和COD均略低于PC區(qū),沉積物TN、TP和OC均顯著低于PC區(qū)(P<0.05,表1)。這可能是克氏原螯蝦和中華鱉的生理特征、生活習(xí)性及管理方式的不同造成的。由于克氏原螯蝦的雜食性,在養(yǎng)殖過程中存在餌料的過量投放和產(chǎn)生大量排泄物等問題[42]。另一方面,中華鱉生長發(fā)育緩慢,餌料投放較少,并且其代謝緩慢[43]?？耸显r挖洞刨沙的習(xí)性會導(dǎo)致沉積物中的營養(yǎng)物質(zhì)更容易釋放到水體中[21]。因此,養(yǎng)殖中華鱉的水環(huán)境更容易保持較低的營養(yǎng)負荷,使細菌群落α多樣性更高。綜上,由于PC和PS區(qū)宿主不同分別形成了獨特的微環(huán)境(排泄物、皮膚分泌物等對細菌獨特的趨化和富集作用),及生理特征、生活習(xí)性和管理方式不同導(dǎo)致的水環(huán)境營養(yǎng)鹽負荷的差異,二者共同作用導(dǎo)致了PS區(qū)的細菌群落α多樣性顯著高于PC區(qū)。

3.2 養(yǎng)殖不同物種的水體和沉積物中的優(yōu)勢OTUs

本研究篩選出的7個沉積物優(yōu)勢OTUs中,OTU11213、OTU59770和103578均與去除沉積物營養(yǎng)濃度有關(guān)。在本研究中,OTU11213是沉積物中平均相對豐度最高的(圖5),其屬于Dechloromonas屬,Dechloromonas廣泛存在于水產(chǎn)養(yǎng)殖系統(tǒng)中并參與氮循環(huán)[56]。OTU59770屬于Prolixibacteraceae科,Prolixibacteraceae主要參與氮循環(huán),能在缺氧的條件下,將硝酸鹽還原為亞硝酸鹽[57]。OTU103578屬于厭氧繩菌科(Anaerolineaceae),是沉積物和濕地土壤中常見的優(yōu)勢類群,主要參與沉積物中氮素的去除[58]。

4 結(jié)論

(1)由于宿主的不同和營養(yǎng)濃度的差異,使得PC和PS區(qū)的水體和沉積物中細菌群落的α多樣性呈現(xiàn)顯著(P<0.05)差異。并且在水體和沉積物中,細菌群落α多樣性均為PS區(qū)顯著高于的PC區(qū)。