基于Miseq測序技術(shù)分析黃顙魚不同養(yǎng)殖模式下池塘微生物群落結(jié)構(gòu)多樣性

王一亭 李 波 王厚紅 吳法龍 張 磊

(1.中國科學(xué)院水生生物研究所中國科學(xué)院水生生物多樣性與保護(hù)重點(diǎn)實(shí)驗(yàn)室，武漢 430072;2.中國科學(xué)院水生生物研究所淮安研究中心，淮安 223002;3.武漢市農(nóng)業(yè)科學(xué)院，武漢 430207;4.淮安市洪澤區(qū)三河鎮(zhèn)農(nóng)業(yè)技術(shù)服務(wù)站，淮安 223124)

微生物群落是養(yǎng)殖生態(tài)系統(tǒng)的重要組成部分[1],其在一定程度上決定養(yǎng)殖生態(tài)系統(tǒng)的穩(wěn)定性[2],并能夠較好地指示養(yǎng)殖環(huán)境健康狀況[3],因此,通過養(yǎng)殖環(huán)境微生物群落結(jié)構(gòu)的研究來比較不同養(yǎng)殖模式對養(yǎng)殖環(huán)境及魚類健康的影響受到越來越多學(xué)者的關(guān)注[4,5]。

黃顙魚Pelteobagrus fulvidraco(Richardson)是典型的底棲雜食性魚類,是我國重要的水產(chǎn)養(yǎng)殖品種。近年來,我國黃顙魚養(yǎng)殖產(chǎn)業(yè)發(fā)展迅速,據(jù)《2018年中國漁業(yè)統(tǒng)計(jì)年鑒》,2017年全國黃顙魚總產(chǎn)量超過48×107kg,產(chǎn)量位居淡水名優(yōu)魚類第三。目前,黃顙魚養(yǎng)殖模式以黃顙魚單養(yǎng)、黃顙魚混養(yǎng)草魚及黃顙魚混養(yǎng)草魚、鳙、鰱三種養(yǎng)殖模式最為常見。草魚是典型的草食性魚類,有助于清除池塘雜草;鰱、鳙是典型濾食性魚類,在控制水體生態(tài)環(huán)境及水質(zhì)健康中起重要作用[6,7]。通過比較黃顙魚不同養(yǎng)殖模式下池塘水體及底泥中微生物的群落結(jié)構(gòu)及差異,有助于掌握不同養(yǎng)殖模式環(huán)境微生物群落結(jié)構(gòu)的變化規(guī)律,從而有效調(diào)整或改進(jìn)現(xiàn)有養(yǎng)殖模式。高通量測序(HTS)技術(shù)由于其快速、準(zhǔn)確、數(shù)據(jù)量大的優(yōu)點(diǎn),被廣泛應(yīng)用到動(dòng)物個(gè)體微生物及環(huán)境微生物基因組學(xué)、功能基因組學(xué)的研究中[8,9]。對全面分析樣本的微生物群落多樣性等意義重大[10]。本研究基于構(gòu)建16S rDNA基因測序的方法,分析黃顙魚不同養(yǎng)殖模式下池塘養(yǎng)殖水體及底泥中微生物群落結(jié)構(gòu)及差異,旨在探討三種養(yǎng)殖模式下環(huán)境微生物組成及變化規(guī)律,為黃顙魚生態(tài)綠色健康養(yǎng)殖模式開發(fā)提供理論依據(jù)。

1 材料與方法

1.1 實(shí)驗(yàn)池塘及實(shí)驗(yàn)魚

實(shí)驗(yàn)池塘位于江蘇省淮安市淮陰區(qū)水產(chǎn)良種場,選擇9個(gè)標(biāo)準(zhǔn)化養(yǎng)殖池塘,每個(gè)池塘面積為6666.7 m2,水深2.2 m,進(jìn)排水方便。池塘水引自同一進(jìn)水渠(水源為運(yùn)河水),池底淤泥約30 cm,底質(zhì)一致性良好。在實(shí)驗(yàn)期間,水溫為自然水溫。實(shí)驗(yàn)所用魚苗由中國科學(xué)院水生生物研究所淮安研究中心自行繁殖的同批次魚苗,實(shí)驗(yàn)飼料由潛江加益加水產(chǎn)飼料有限公司生產(chǎn),飼料蛋白含量為40%。

1.2 實(shí)驗(yàn)設(shè)計(jì)

實(shí)驗(yàn)分為黃顙魚單養(yǎng)組A、黃顙魚草魚混養(yǎng)組B、黃顙魚草魚鰱鳙混養(yǎng)組C三種養(yǎng)殖模式,每種養(yǎng)殖模式分別由3個(gè)獨(dú)立的養(yǎng)殖池塘進(jìn)行實(shí)驗(yàn)(即每種養(yǎng)殖模式3個(gè)重復(fù))。每種魚類的放養(yǎng)平均質(zhì)量及放養(yǎng)密度為: 黃顙魚(8.20±0.10) g、15 ind./m2,草魚(10.03±0.04) g、0.05 ind./m2,鰱(9.85±0.11) g、0.04 ind./m2,鳙(11.56±0.32) g、0.03 ind./m2。實(shí)驗(yàn)時(shí)間為2017年8月至11月,共84d。實(shí)驗(yàn)期間每天定時(shí)投喂2次(9:00和16:00),每次投喂量約為黃顙魚魚體總重的4%,視黃顙魚吃食情況適量增減投喂量,并記錄。

1.3 樣品采集和處理

實(shí)驗(yàn)開始、結(jié)束時(shí)分別對池塘養(yǎng)殖水體及底泥取樣1次。水樣采用5 L采水器分別在池塘的中央和4個(gè)角水面下50 cm處采集,并將5個(gè)點(diǎn)的水樣混合后取2 L帶回實(shí)驗(yàn)室待測。在水樣沉淀后,取上清液通過0.22 μm微孔濾膜真空過濾,將過濾后的濾膜保存在5 mL離心管中,液氮速凍3—5min,然后轉(zhuǎn)移至-80℃冰箱中保存待測。泥樣采集用彼得森采泥器分別在池塘中央和4個(gè)角采集水底表層泥,混合后裝入無菌收集瓶。泥樣帶回實(shí)驗(yàn)室后分裝至 5 mL離心管中,液氮速凍3—4h,然后轉(zhuǎn)移至-80 ℃冰箱中保存待測。樣品命名見表1。

1.4 樣品DNA提取

水樣中微生物DNA提取,取出樣品濾膜充分剪碎(使樣品表面菌群更好的與裂解液接觸),置于1.5 mL滅菌離心管中,加入裂解液,使其釋放DNA,隨后使用E.Z.N.A.TMWater DNA Kit(Omega,美國)試劑盒提取;泥樣中微生物DNA提取使用E.Z.N.ATMMag-Bind Soil DNA Kit(Omega,美國)試劑盒,具體方法見試劑盒說明書。

1.5 PCR擴(kuò)增及測序分析

本研究選擇細(xì)菌16S rDNA V3-V4可變區(qū)作為目的擴(kuò)增片段,進(jìn)行后續(xù)高通量測序,采用16S rDNA V3-V4區(qū)通用引物341F/805R。引物序列為:341F: 5′-CCCTACACGACGCTCTTCCGATCTG-3′,805R: 5′-GACTGGAGTTCCTTGGCACCCGAGA ATTCCA-3′,同時(shí)在3′端添加Barcode區(qū)分序列樣品。

擴(kuò)增V3-V4可變區(qū)采用30 μL PCR體系。PCR體系為: 2×Taqmaster Mix 15 μL,primer F(10 μmol/L)1 μL,Primer R (10 μmol/L) 1 μL,DNA模板 10—20 ng,添加去離子水至30 μL。

PCR條件為: 95℃變性3min;94 ℃變性20s,55℃退火20s,72℃延伸30s,擴(kuò)增25個(gè)循環(huán);最后于72℃延伸5min。擴(kuò)增后PCR產(chǎn)物經(jīng)純化后,交生工生物工程有限公司(上海)進(jìn)行Miseq(Illumina)高通量測序。

表 1 水樣及泥樣具體命名表Tab.1 List of names of water samples and mud samples

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

Illumina Miseq測序得到的原始圖像數(shù)據(jù)文件經(jīng)堿基識別(Base Calling)分析轉(zhuǎn)化為原始測序序列(Raw Reads),根據(jù)barcode區(qū)分每個(gè)樣本,并對序列質(zhì)量進(jìn)行質(zhì)控和過濾;根據(jù)overlap關(guān)系進(jìn)行拼接,拼接后的序列再次進(jìn)行質(zhì)控和過濾,最后得到優(yōu)化序列,進(jìn)行下游分析。優(yōu)化序列進(jìn)行OTU 聚類分析和物種分類學(xué)分析。

OTU(操作分類單元,Operational Taxonomic Units)是在系統(tǒng)發(fā)生學(xué)或群體遺傳學(xué)研究中,為了便于進(jìn)行分析,人為給某一個(gè)分類單元(品系,屬,種、分組等)設(shè)置的同一標(biāo)志。根據(jù)不同的相似度水平,對所有序列進(jìn)行OTU劃分,本研究在97%的相似水平下的OTU進(jìn)行生物信息統(tǒng)計(jì)分析。

2 結(jié)果

2.1 序列豐度

對測序所得的序列進(jìn)行拼接、質(zhì)控、去除低質(zhì)量序列和嵌合體,12個(gè)樣本共獲得655100條合格16S rDNA序列,每個(gè)樣本的有效序列數(shù)目41256—77508條(表2)。細(xì)菌的有效序列長度主要分布在400—440 bp(圖1)。通過97%相似度的OTU進(jìn)行稀釋性曲線分析,結(jié)果發(fā)現(xiàn)各個(gè)樣本曲線均趨向平坦,說明測序數(shù)據(jù)量合理,更多的數(shù)據(jù)量只會產(chǎn)生少量新的OTU,因此測序所得序列能夠完全展示樣品微生物群落多樣性(圖2)。

表 2 樣品序列數(shù)統(tǒng)計(jì)及菌群多樣性分析Tab.2 Numbers of reads and cliversity indes of different samples

2.2 細(xì)菌群落的α-多樣性分析

利用測序劃分的OTU對不同養(yǎng)殖模式下水體和底泥的α-多樣性分析結(jié)果如表2所示: 所有樣品的覆蓋度均在97%以上,在三種養(yǎng)殖模式下底泥樣本在實(shí)驗(yàn)前后OTU數(shù)量變化不大,水體樣本OTU數(shù)量在養(yǎng)殖后期明顯增加。

Chao1和ACE指數(shù)反映群落分布豐度,指數(shù)越大,豐度越高,Shannon和Simpson指數(shù)指示群落分布多樣性,Shannon指數(shù)越大、Simpson指數(shù)越小,說明群落多樣性越高。在三種養(yǎng)殖模式下底泥樣本的Chao1、ACE和Shannon均高于水體樣本,說明底泥的群落豐度及多樣性明顯高于水體。養(yǎng)殖初期與末期樣本群落豐度比較,底泥樣本的Chao1和A C E 指數(shù)前后變化不明顯,而水體末期樣本(WAE、WBE和WCE)的Chao1和ACE指數(shù)明顯高于初期樣本,其中WBE最高,WCE最低。養(yǎng)殖初期與末期樣本群落多樣性比較,底泥樣本的Shannon和Simpson指數(shù)前后變化不明顯,而水體末期樣本(WAE、WBE和WCE)的Shannon指數(shù)明顯高于初期樣本,其中WCE最高,WAE最低;水體末期樣本的Simpson指數(shù)明顯低于初期樣本,其中WCE最低,WAE較高。因此,在三種養(yǎng)殖模式下水體中末期的群落豐度和多樣性均高于養(yǎng)殖初期,并且黃顙魚草魚混養(yǎng)組水體群落豐度最高,黃顙魚草魚鰱鳙混養(yǎng)組水體群落多樣性較高。

2.3 水樣與底泥樣品菌群之間的關(guān)系

基于OTU豐度的樣本聚類圖可知: 底泥和水體樣本分別聚為一簇,底泥樣本間的空間距離較短,相似性更高。水體和底泥初期樣本與末期樣本各自具有較高的相似性。三種養(yǎng)殖模式對比,A、B模式與C模式下的樣本表現(xiàn)出更大的空間距離(圖3)。

為進(jìn)一步探究每個(gè)樣本中細(xì)菌多樣性的相互關(guān)系,通過維恩分析觀察水體和底泥樣本OTU數(shù)目組成相似性及重疊情況。

對水體初期和末期各樣本中只有一條序列的OTU進(jìn)行統(tǒng)計(jì)分析,發(fā)現(xiàn)其數(shù)量已經(jīng)達(dá)到總有效OTU數(shù)量的60%以上,表明各樣本細(xì)菌種類多但大多數(shù)細(xì)菌豐度較小。養(yǎng)殖初期各水樣共有OTU數(shù)為151條(圖4a),占三者OTU總數(shù)的8.36%,養(yǎng)殖末期各水樣共有OTU數(shù)為476條(圖4b),占三者OTU總數(shù)的14.88%,表明隨著養(yǎng)殖活動(dòng)的進(jìn)行,水體樣本間細(xì)菌群落相似度有所增加。

從底泥樣本的菌群相似性來看,養(yǎng)殖初期各泥樣共有OTU數(shù)為1568條(圖4c),占三者OTU總數(shù)的31.17%,養(yǎng)殖末期各泥樣共有OTU數(shù)為1104條(圖4d),占三者OTU總數(shù)的21.64%,養(yǎng)殖前后底泥樣本之間共有OTU的比例均較高,說明在不同養(yǎng)殖時(shí)期、不同養(yǎng)殖模式下底泥群落相似性均比較高,但養(yǎng)殖后期三者共有OTU數(shù)占比下降,說明不同的養(yǎng)殖魚類使底泥樣本微生物群落結(jié)構(gòu)逐漸分化。

PcoA分析各樣品間的差異和距離,從圖5a可以看出,6個(gè)水體樣本可以分為2個(gè)相似群落,其中初期樣本和末期樣本分別聚集在一起,而6個(gè)泥樣樣本全部聚集在一起,這與聚類分析的結(jié)果一致。對水樣進(jìn)行PcoA分析發(fā)現(xiàn)(圖5b),沿PC1坐標(biāo)方向(貢獻(xiàn)率為37%),初期樣本與末期樣本分別趨于收斂,沿PC2坐標(biāo)方向(貢獻(xiàn)率為20%),初期樣本得以分開,說明黃顙魚的養(yǎng)殖使水體微生物群落構(gòu)成相似度逐漸升高。對泥樣進(jìn)行分析發(fā)現(xiàn)(圖5c),養(yǎng)殖初期各樣本在PC1(貢獻(xiàn)率為32%)及PC2(貢獻(xiàn)率為24%)方向均高度收斂,距離很近,說明初期樣本相似度較高,而末期各樣本在2個(gè)方向上都趨于分化,三種養(yǎng)殖模式下底泥的微生物群落構(gòu)成相似度降低,這與聚類分析及維恩分析的結(jié)果一致。

圖1 樣品中有效序列長度分布Fig.1 Distribution of bacterial sequence lengths in the samples

圖2 各樣品高通量測序結(jié)果的稀釋性曲線Fig.2 The rarefaction curve basing on high-through sequencing results

圖3 基于高通量測序結(jié)果的聚類分析結(jié)果Fig.3 Clustering analysis results basing on high-through sequencing results of different samples

2.4 各樣本微生物菌群構(gòu)成

高通量測序結(jié)果表明水樣和泥樣的核心微生物群明顯不同,養(yǎng)殖前后樣本之間也存在較大差異。所有樣本的序列分屬43個(gè)門,其中共有門有8個(gè),分別為變形菌門(Proteobacteria)、綠彎菌門(Chloroflexi)、厚壁菌門(Firmicutes)、酸桿菌門(Acidobacteria)、擬桿菌門(Bacteroidetes)、疣微桿菌門(Verrucomicrobia)、浮霉菌門(Planctomycetes)、放線菌門(Actinobacteria)。除此之外,還檢測到廣古菌門(Euryarchaeota)、梭桿菌門(Fusobacteria)、芽單胞菌門(Gemmatimonadetes)、硝化螺旋菌門(Nitrospirae)、衣原體門(Chlamydiae)、螺旋體門(Spirochaetes)、綠菌門(Chlorobi)、互養(yǎng)菌門(Synergistetes)等不常見的類群,在樣本中檢測量極低,這表明基于Illumina Miseq平臺的高通量測序可作為全面了解黃顙魚養(yǎng)殖環(huán)境微生物結(jié)構(gòu)多樣性的有效研究手段。

對水體樣本在門水平的物種相對豐度(圖6a)分析發(fā)現(xiàn),各樣本中豐度占比最大的是變形菌門,但其在養(yǎng)殖末期占比明顯下降,3種養(yǎng)殖模式的水體樣本分別從初期的75%、91.39%和80.17%,下降到養(yǎng)殖末期的39.5%、42.64%和30.14%。同樣明顯下降的還有厚壁菌門。豐度上升的有放線菌門、浮霉菌門、疣微桿菌門、綠彎菌門、酸桿菌門。其中,A模式水樣擬桿菌門豐度明顯上升(1.95%到10.98%),而B模式變化不大(2.11%到2.39%),C模式則略有下降(1.63%到0.76%)。從屬水平(圖6b)分析水樣的組成及豐度發(fā)現(xiàn),養(yǎng)殖末期水體較養(yǎng)殖前期水體的細(xì)菌屬更豐富。部分水體樣本中包含了一些致病菌或條件致病菌,如弧菌屬(Vibrio)、假單胞菌屬(Pseudomonas)、氣單胞菌屬(Aeromonas)等。其中,A模式下水體樣本中黃顙魚條件致病菌氣單胞菌屬(Aeromonas)在養(yǎng)殖末期豐度占優(yōu)勢(5.81%),而其B、C模式下檢出量極少甚至未檢出。

對底泥樣本在門水平的物種組成豐度進(jìn)行分析(圖6a)發(fā)現(xiàn),底泥樣本相似度較高,豐度占比最大的門同樣為變形菌門,但其在養(yǎng)殖后期豐度均略有上升。其他豐度較高的還有厚壁菌門、綠彎菌門、浮霉菌門等。其中A、B模式酸桿菌門豐度下降,而C模式則有所上升。從屬水平分析(圖6b),三種養(yǎng)殖模式泥樣的假單胞桿菌屬在養(yǎng)殖末期豐度顯著增加。A、C模式不動(dòng)桿菌屬(Acinetobacter)豐度下降,而B模式豐度升高。A、B模式的Litorilinea屬豐度有所下降,而C模式下豐度升高。

圖4 維恩圖Fig.4 Bacterial Venn diagram at distance 0.03

3 討論

本研究提取了黃顙魚不同養(yǎng)殖模式下養(yǎng)殖水體和底泥的微生物宏基因組DNA,采用Illumina Miseq高通量測序平臺分析了微生物多樣性。得到有效序列41256—77508條,可歸納為608—3686個(gè)分類操作單元?；跍y序數(shù)據(jù),評估了黃顙魚不同養(yǎng)殖模式下養(yǎng)殖水體和底泥的微生物多樣性及群落結(jié)構(gòu)差異。

3.1 養(yǎng)殖水體及底泥微生物群落結(jié)構(gòu)及多樣性分析

維恩圖分析及主坐標(biāo)分析表明,隨著養(yǎng)殖活動(dòng)的進(jìn)行,不同養(yǎng)殖模式下的水體樣本共有OTU比例上升,水體樣本在主坐標(biāo)分析圖中趨于收斂,表明水體樣本間相似性逐漸增加。底泥樣本雖然初期樣本相似度較高,但后期樣本共有OTU比例下降,樣本在主坐標(biāo)分析圖中距離增大,表明樣本細(xì)菌群落分布趨于分化,這可能與黃顙魚的底棲特性相關(guān)。由于池塘底泥與水接觸層中微生物含量豐富[11],不同養(yǎng)殖模式造成的微生物群落分布差異性可能通過黃顙魚對底層水體及底泥的攪動(dòng)而在底泥樣本中被放大,顯示出分化現(xiàn)象[12]。

樣本聚類分析結(jié)果表明,黃顙魚單養(yǎng)及黃顙魚、草魚混養(yǎng)模式下的水體和底泥的菌群組成及群落結(jié)構(gòu)相似性較高,黃顙魚、草魚、鰱、鳙混養(yǎng)模式下樣本單獨(dú)形成一類,即樣本聚類按照有無鰱鳙聚在一起,表明鰱鳙對混養(yǎng)水體菌群組成和群落結(jié)構(gòu)有更明顯的影響。這與唐永濤等[6]關(guān)于混養(yǎng)鰱鳙對池塘養(yǎng)殖水體微生物群落結(jié)構(gòu)影響的結(jié)果一致。

對樣本的α-多樣性分析結(jié)果表明,不同養(yǎng)殖模式下的黃顙魚養(yǎng)殖均可顯著促進(jìn)養(yǎng)殖水體微生物群落豐度和多樣性增加,但水體樣本相似度也有所增加。不同養(yǎng)殖模式末期樣本之間相比較,黃顙魚草魚混養(yǎng)模式下水體的群落豐度最高,引入鰱、鳙混養(yǎng)后水體微生物群落豐度明顯下降,這與Xia等[13]研究發(fā)現(xiàn)草魚混養(yǎng)鰱魚可以有效減輕草魚養(yǎng)殖對環(huán)境的影響結(jié)論一致。草魚是一種食草性魚類,主要食用自然環(huán)境中的水生植物,但有研究表明草魚更喜歡養(yǎng)殖池塘中的顆粒飼料[14]。草魚與黃顙魚混養(yǎng),不僅會與黃顙魚爭奪飼料,其排泄的飼料殘?jiān)?、魚糞等廢物經(jīng)過發(fā)酵分解,為細(xì)菌提供適宜環(huán)境,使水體微生物大量生長,進(jìn)而培育出浮游生物,敗壞水質(zhì)[13]。而適量的濾食性魚類鰱鳙引入養(yǎng)殖模式后,鰱鳙可以將水體中富余的微生物和浮游生物作為餌料,有效減輕草魚養(yǎng)殖對環(huán)境的影響[15]。同時(shí),黃顙魚、草魚、鰱、鳙混養(yǎng)模式下的水體微生物群落多樣性最高,表明鰱鳙的引入還可以增加水體微生物群落多樣性。這與田相利等[16]發(fā)現(xiàn)的草魚、鰱、鯉(Cyprinus carpio,Linnaeus)混養(yǎng)使養(yǎng)殖環(huán)境中微生物的結(jié)構(gòu)組成與代謝功能更加多樣化的結(jié)論一致。

圖5 PcoA分析Fig.5 PcoA analysis

3.2 養(yǎng)殖水體及底泥微生物門水平和屬水平差異

在本研究中,門類水平上的測序結(jié)果表明,在三種養(yǎng)殖模式下水體及底泥樣本中,變形菌門的物種豐度明顯高于其他微生物豐度,尤其是在養(yǎng)殖初期占絕對地位,這與Hugenholtz等[17]、陳香碧等[18]認(rèn)為的變形菌門是自然環(huán)境中最常見菌群結(jié)論一致。同時(shí)變形菌門在很多養(yǎng)殖環(huán)境中也占據(jù)優(yōu)勢地位[19],趙曉偉等[20]對紅鰭東方鲀(Takifugu rubripes,Temminck &Schlegel)養(yǎng)殖水體和底泥的微生物群落分析結(jié)果表明,變形菌門為優(yōu)勢菌。但隨著養(yǎng)殖時(shí)間的延長,各種細(xì)菌的分布逐漸均勻,變形菌門細(xì)菌優(yōu)勢逐漸減弱。

放線菌門在物質(zhì)代謝上具有獨(dú)特的優(yōu)點(diǎn),在微生物修復(fù)水體環(huán)境中也有應(yīng)用[21],尤其以溶藻作用在防治赤潮方面具有重要作用[22]。在本研究中,放線菌門在養(yǎng)殖初期豐度不高,但到養(yǎng)殖末期,卻成為水體樣本中占比第二的優(yōu)勢門,尤其在黃顙魚草魚混養(yǎng)模式下,放線菌豐度占比超過三分之一。這可能是由于魚類養(yǎng)殖過程中,養(yǎng)殖水體中的殘餌、魚糞等逐漸富集,造成浮游植物生物量變大,為放線菌的生長繁殖提供營養(yǎng)。

擬桿菌門是溶解性有機(jī)物的主要消費(fèi)者。有研究表明,養(yǎng)殖環(huán)境中殘留餌料及養(yǎng)殖動(dòng)物糞便在水體中未被及時(shí)分解并沉積在水體及底泥中時(shí),會造成養(yǎng)殖池塘底部擬桿菌門細(xì)菌含量增加[23]。本研究黃顙魚單養(yǎng)及混養(yǎng)草魚時(shí),養(yǎng)殖末期水體及底泥擬桿菌門細(xì)菌含量均明顯增加,這可能是因?yàn)榛祓B(yǎng)草魚后餌料殘?jiān)汪~糞沉積在水體尤其是底泥中,導(dǎo)致擬桿菌門細(xì)菌大量繁殖。而黃顙魚、草魚、鰱、鳙混養(yǎng)模式下的環(huán)境樣本中擬桿菌門細(xì)菌含量大幅下降,指示濾食性魚類鰱鳙的引入有效消耗水體中的草魚糞便及殘余餌料[13],抑制擬桿菌門細(xì)菌生長。

酸桿菌門是土壤中重要的一類微生物,是苯酚的主要降解菌[24]。有學(xué)者推測酸桿菌門的部分菌種對纖維素有一定的分解能力,而纖維素是植物細(xì)胞壁的主要組成成分[22]。在本研究中,底泥樣本中酸桿菌門占比較高,這可能是由于底泥中沉積的餌料殘?jiān)绕涫遣蒴~糞便中含有大量纖維素。

在本研究中,厚壁菌門也是樣本中豐度較高的優(yōu)勢門,但近年來針對該門的研究主要集中在其作為主要腸道菌群與肥胖的關(guān)系上[25]。對于本研究中養(yǎng)殖水體厚壁菌門豐度下降的原因還有待進(jìn)一步研究。

水產(chǎn)養(yǎng)殖動(dòng)物細(xì)菌性疾病的爆發(fā)往往與腸道細(xì)菌生態(tài)平衡異常相關(guān)[26]。養(yǎng)殖魚類腸道菌群組成與養(yǎng)殖池塘底泥及水體微生物群落結(jié)構(gòu)組成密切相關(guān)[27],尤其是淡水魚腸道內(nèi)優(yōu)勢菌與養(yǎng)殖水體優(yōu)勢菌組成一般具有相似性[28]。有研究表明,魚類腸炎的病發(fā)與腸道內(nèi)氣單胞菌大量繁殖相關(guān)[29]。周金敏等[26]對患細(xì)菌性出血病的黃顙魚與健康黃顙魚腸道中的氣單胞菌進(jìn)行比較,發(fā)現(xiàn)患病黃顙魚腸道內(nèi)的氣單胞菌數(shù)量大幅增加。在本研究中黃顙魚單養(yǎng)模式下水體樣本中氣單胞菌在養(yǎng)殖末期豐度占優(yōu)勢,極有可能是單養(yǎng)模式下的黃顙魚腸道內(nèi)氣單胞菌數(shù)量上升所致,該模式下黃顙魚患出血性腸炎的概率大大增加。因此,有針對性的增加一些套養(yǎng)魚類可以降低黃顙魚單養(yǎng)潛在的養(yǎng)殖風(fēng)險(xiǎn)。

圖6 樣品在門和屬水平物種相對豐度Fig.6 Relative abundance of species in phylum and genus level

4 結(jié)論

本研究結(jié)果顯示,不同養(yǎng)殖模式下的黃顙魚養(yǎng)殖環(huán)境微生物均體現(xiàn)了高多樣性和高豐度特征,這可能對池塘微生態(tài)系統(tǒng)穩(wěn)定和致病菌的抑制起到積極調(diào)控作用。然而,在黃顙魚單養(yǎng)模式下水體條件致病菌數(shù)量有所上升,存在潛在致病風(fēng)險(xiǎn);在黃顙魚草魚混養(yǎng)模式下,草魚排泄物造成水體中具有分解有機(jī)質(zhì)功能的部分微生物豐度顯著增大,同時(shí)整體群落豐度也較高,促進(jìn)浮游生物的生長(相關(guān)數(shù)據(jù)未發(fā)表),可能導(dǎo)致水體富營養(yǎng)化升高。只有在適量的濾食性魚類鰱、鳙引入養(yǎng)殖模式后,通過鰱、鳙的混養(yǎng)有效減輕黃顙魚養(yǎng)殖對水體及底泥菌群組成和群落結(jié)構(gòu)的負(fù)面影響,同時(shí)增加水體微生物多樣性,幫助恢復(fù)水質(zhì)健康。因此,從微生物群落結(jié)構(gòu)和多樣性的角度考慮,黃顙魚、草魚、鰱、鳙混養(yǎng)是一種值得推薦的養(yǎng)殖模式。