趙 瓊，張峰峰，周 可，孫海波，趙玉潔，謝鳳行

（天津市農(nóng)業(yè)生物技術(shù)研究中心，天津300384）

微生物具有豐富的物種和遺傳多樣性，在眾多的微生物多樣性的研究方法中，Biolog Eco法是以研究微生物群落水平生理特征（community level physiological profiling，CLPP）為主的一種方法，由美國(guó)BIOLOG公司于1989年建立[1].Biolog Eco方法的原理是微生物在利用碳源過程中產(chǎn)生的自由電子與四唑染料發(fā)生還原顯色反應(yīng)，顏色深淺可以反映微生物群落對(duì)碳源的利用程度.由于微生物對(duì)不同碳源的利用能力很大程度上取決于微生物的種類和固有性質(zhì)，因而Biolog Eco板培養(yǎng)后測(cè)定各孔的顏色變化即可表征樣品微生物群落的代謝功能多樣性與活性[2].Biolog Eco方法數(shù)據(jù)量十分豐富，每一孔的吸光度均可作為樣品的1個(gè)變量進(jìn)行分析，且具有方便、快捷等優(yōu)點(diǎn)，在表征細(xì)菌群落動(dòng)態(tài)變化的時(shí)空尺度上具有顯著優(yōu)勢(shì)[3-4]，因此被廣泛應(yīng)用于土壤、水體等環(huán)境領(lǐng)域[5-8].由于Biolog Eco方法能提供微生物群落碳代謝潛力的相關(guān)信息，因此還可以和高通量測(cè)序方法相結(jié)合，進(jìn)行微生物群落結(jié)構(gòu)和代謝功能的相關(guān)研究[9-11].

Biolog Eco方法最初是被應(yīng)用于土壤微生物群落的研究，近年來在海洋、湖泊、河口領(lǐng)域的應(yīng)用逐年增多[12].Biolog Eco方法的難點(diǎn)在于數(shù)據(jù)分析，目前進(jìn)行數(shù)據(jù)分析的方式主要有3種：①選用特定的培養(yǎng)時(shí)間點(diǎn)去分析Biolog Eco的數(shù)據(jù)[13]；②采用曲線擬合的方法，兼顧Biolog Eco反應(yīng)前后的變化[14]；③為規(guī)避接種密度對(duì)分析結(jié)果的影響采用固定顏色平均變化率（average well color development，AWCD）時(shí)的吸光度值進(jìn)行分析[15].由于選擇特定培養(yǎng)時(shí)間點(diǎn)的方法更簡(jiǎn)單易行，目前被國(guó)內(nèi)外學(xué)者普遍采用[13].根據(jù)不同環(huán)境下樣品微生物群落菌體密度和活性的不同，常用于數(shù)據(jù)分析的培養(yǎng)時(shí)間有72[16]、96[17]、120[18]、168[19]和240 h[20].賈夏等[19]研究土壤樣品Biolog Eco板培養(yǎng)時(shí)間對(duì)分析結(jié)果的影響，認(rèn)為應(yīng)盡量選擇240 h進(jìn)行連續(xù)讀數(shù)全記錄，再根據(jù)AWCD曲線“拐點(diǎn)”處或之后的時(shí)間進(jìn)行分析.

養(yǎng)殖水體的異養(yǎng)菌菌體濃度比土壤環(huán)境低很多，一般在103～105cfu/mL之間[21-22].Biolog Eco方法在養(yǎng)殖水體中的應(yīng)用研究也有少量報(bào)道，如Kurten等[23]認(rèn)為，利用處于指數(shù)增長(zhǎng)期的72 h或96 h進(jìn)行結(jié)果分析較為合適，因?yàn)殡S著時(shí)間的延長(zhǎng)群落底物利用特征（community substrate utilization profile，CSUP）的差異程度降低，AWCD曲線進(jìn)入穩(wěn)定期后提供的新信息有限.目前缺乏對(duì)養(yǎng)殖水體微生物群落代謝功能隨時(shí)間變化的相關(guān)研究.因此，本研究以鰱魚單養(yǎng)池塘和對(duì)蝦精養(yǎng)池塘水樣為研究對(duì)象，考察水樣微生物群落的AWCD值、6類碳源利用情況、多樣性指數(shù)和主成分分析達(dá)到?jīng)]有顯著差異的時(shí)間點(diǎn)，以及不同時(shí)間點(diǎn)對(duì)5種水樣的區(qū)分情況.本研究可為利用Biolog Eco板開展養(yǎng)殖水體微生物群落功能多樣性研究提供依據(jù)，也為Biolog Eco板在其他水體環(huán)境中的應(yīng)用提供參考.

1 材料與方法

1.1 養(yǎng)殖池塘的選擇及管理

本研究共采集了5個(gè)水樣.WS1和WS2水樣取自天津市西青區(qū)津靜公路北側(cè)魚塘，東經(jīng)117°05′，北緯39°10′，取樣時(shí)間為2018年5月31日.WS1池塘為鰱魚單養(yǎng)池塘，WS2池塘為南美白對(duì)蝦精養(yǎng)池塘，2個(gè)池塘面積均為4 hm2.2018年5月20日WS1池塘放養(yǎng)鰱魚4 500尾/hm2，WS2池塘放養(yǎng)蝦苗60萬尾/hm2.養(yǎng)殖期間WS1池塘不投喂飼料，WS2池塘投喂蝦飼料.WS3、WS4和WS5水樣均取自天津市西青區(qū)畜牧水產(chǎn)局養(yǎng)殖基地南美白對(duì)蝦精養(yǎng)池塘，東經(jīng)117°04′，北緯39°11′，取樣時(shí)間為2018年9月10日.3個(gè)池塘面積均為0.5 hm2，2018年6月20日每個(gè)池塘放養(yǎng)蝦苗30萬尾，養(yǎng)殖期間投喂蝦飼料.表1為5種水樣池塘水質(zhì)的基本情況和異養(yǎng)菌菌落總數(shù).

表1 池塘水體理化性質(zhì)和異養(yǎng)菌菌數(shù)Tab.1 Physical chemical characteristics and heterotrophic plate count of the water sample in ponds

1.2 樣品采集

用有機(jī)玻璃采水器在池塘四周和中心分別采集等量表層水樣（水面以下0.5 m），混和后裝入無菌取樣瓶，每個(gè)池塘取2份水樣，低溫保存，帶回實(shí)驗(yàn)室立即開展實(shí)驗(yàn).

1.3 測(cè)定方法

異養(yǎng)菌總數(shù)的測(cè)定：將所取池塘水樣混和，利用0.85%滅菌生理鹽水進(jìn)行梯度稀釋，取合適稀釋度的0.1 mL稀釋液涂布營(yíng)養(yǎng)瓊脂平板，28℃倒置培養(yǎng)3 d計(jì)算菌落總數(shù).

Biolog Eco板操作：用0.85%的生理鹽水對(duì)水樣進(jìn)行稀釋，用八道移液器將樣液加入預(yù)熱的Biolog Eco板中，每孔加150μL.將加好樣的Biolog Eco板加蓋放入保鮮盒中，25℃培養(yǎng)，每隔24 h利用酶標(biāo)儀（Elx808，美國(guó)Bio Tek公司）讀取各孔在750 nm和590 nm下的吸光度值[24].每個(gè)池塘的2份水樣分別加入到Biolog Eco板中，每個(gè)Biolog Eco板3個(gè)重復(fù)，共6個(gè)重復(fù).

Biolog Eco板培養(yǎng)時(shí)間的優(yōu)化：將所取水樣以原液加入Biolog Eco板中，記錄24～192 h的吸光度值，計(jì)算不同時(shí)間點(diǎn)微生物群落的AWCD值、6類碳源利用情況、多樣性指數(shù)并進(jìn)行主成分分析，確定各項(xiàng)指標(biāo)達(dá)到?jīng)]有顯著差異的時(shí)間點(diǎn)，同時(shí)研究不同時(shí)間點(diǎn)對(duì)5種水樣的區(qū)分.

1.4 計(jì)算方法

1.4.1 水樣微生物群落平均活性

Biolog Eco板具有96個(gè)孔31種碳源，每種碳源設(shè)3個(gè)重復(fù).用（OD590-OD750）的值表示微生物代謝活性，OD590和OD750的值為加樣孔吸光度值減去相應(yīng)對(duì)照孔的吸光度值，當(dāng)上述差值為負(fù)或小于0.06時(shí)修正為0[24].微生物的平均活性采用Biolog Eco孔的AWCD來描述，計(jì)算公式如下：

式（1）中：Ci為每個(gè)碳源孔測(cè)定的兩波段吸光度值的差值；R為對(duì)照孔的吸光度值；31為碳源種類數(shù).

1.4.2 多樣性指數(shù)分析

采用Shannon-Wiener指數(shù)（H′）、Pielou均勻度指數(shù)（J）、Simpson指數(shù)（D）、McIntosh指數(shù)（Dmc）和McIntosh均勻度（E）來描述微生物群落的功能多樣性，計(jì)算公式如下[25-26]：

式（2）～式（6）中：ni為每個(gè)孔的相對(duì)吸光度值即（Ci-R）；N為Biolog Eco板上每個(gè)重復(fù)區(qū)域（32個(gè)孔）相對(duì)吸光度值的總和即∑（Ci-R）；pi為第i孔相對(duì)吸光度值與其相應(yīng)重復(fù)的相對(duì)吸光度值總和的比值，即pi=為顏色變化的孔的數(shù)目，即Ci-R＞0的孔的數(shù)目；在計(jì)算Simpson指數(shù)時(shí)，數(shù)據(jù)擴(kuò)大了1 000倍以防止出現(xiàn)負(fù)數(shù)[25].

1.4.3 碳源利用情況分析

Biolog Eco板含有的31種碳源可以分為6大類，其中氨基酸（amino acids）6種、碳水化合物類（carbohydrates）12種（D-葡萄糖胺酸和D-半乳糖醛酸列入碳水化合物類）、羧酸（carboxyl acids）5種、聚合物（polymers）4種、胺類（amines）2種和酚類2種[5]，計(jì)算每類碳源對(duì)應(yīng)孔（OD590-OD750）的平均值，以此表征微生物群落對(duì)此類碳源的利用情況.

1.4.4 主成分分析

對(duì)5個(gè)池塘水樣72～192 h微生物群落的底物利用特征進(jìn)行主成分分析（principal components analysis，PCA），同時(shí)利用主成分分析比較72～192 h時(shí)間點(diǎn)對(duì)5種水樣的區(qū)分.采用協(xié)方差矩陣進(jìn)行主成分分析[27-28].

1.4.5 統(tǒng)計(jì)分析

用Excel 2016和SPSS19.0進(jìn)行數(shù)據(jù)統(tǒng)計(jì)，每種水樣的6組吸光度值數(shù)據(jù)互為平行[24].AWCD值、6類碳源的利用和多樣性指數(shù)均采用單因子方差分析進(jìn)行差異統(tǒng)計(jì)學(xué)意義檢驗(yàn)，P＜0.05為差異具有統(tǒng)計(jì)學(xué)意義.ANOSIM方法（analysis of similarities）用于主成分分析后微生物群落底物利用特征的差異統(tǒng)計(jì)學(xué)意義檢驗(yàn)，SIMPER分析（similarity percentages）用于不同時(shí)間群落底物利用特征的差異度分析，2種方法通過PAST Version 3軟件完成[29]，應(yīng)用Bray-Curtis距離指數(shù).

2 結(jié)果與分析

2.1 不同培養(yǎng)時(shí)間微生物群落的AWCD值

5種水樣微生物群落AWCD值隨培養(yǎng)時(shí)間的變化如圖1所示.

圖1 微生物群落AWCD值的變化曲線Fig.1 Change curve of AWCD value of microbial community

AWCD曲線在經(jīng)過一定時(shí)間的延遲后迅速升高，進(jìn)入指數(shù)增長(zhǎng)期，隨著培養(yǎng)時(shí)間的延長(zhǎng)趨于穩(wěn)定.其中，WS1水樣最早進(jìn)入穩(wěn)定期，96 h后各時(shí)間點(diǎn)AWCD值沒有顯著差異.WS2、WS3和WS5水樣進(jìn)入指數(shù)增長(zhǎng)的時(shí)間比WS1水樣延遲，AWCD值趨于穩(wěn)定的時(shí)間分別為168、168和120 h.WS4水樣AWCD值顯著低于其他水樣，在144 h后各時(shí)間點(diǎn)的AWCD值沒有顯著差異.

2.2 不同培養(yǎng)時(shí)間6類碳源利用情況

5種水樣72～192 h的碳源利用情況如圖2所示.由圖2可以看出，WS1水樣微生物群落對(duì)氨基酸類碳源的利用在96 h后沒有顯著差異，而其他各類碳源的利用在72 h后沒有顯著差異，即6類碳源在96 h后沒有顯著差異；WS2水樣微生物群落對(duì)碳水化合物類、氨基酸類、羧酸類、多聚物類、酚酸類和胺類的利用沒有顯著差異的時(shí)間點(diǎn)分別為96、120、120、192、144、120 h，即6類碳源的利用在192 h后沒有顯著差異；WS3水樣微生物群落對(duì)6類碳源的利用沒有顯著差異的時(shí)間點(diǎn)分別為144、168、144、144、120、120 h，即6類碳源的利用在168 h后沒有顯著差異；WS4水樣微生物群落對(duì)6類碳源的利用沒有顯著差異的時(shí)間點(diǎn)分別為120、144、120、144、144、168 h，即6類碳源的利用在168 h后沒有顯著差異；WS5水樣微生物群落對(duì)6類碳源的利用沒有顯著差異的時(shí)間點(diǎn)分別為96、72、72、120、72、72 h，即6類碳源的利用在120 h后沒有顯著差異.

對(duì)比不同時(shí)間點(diǎn)水樣碳源利用的大小順序，發(fā)現(xiàn)144 h后WS1水樣碳源利用順序穩(wěn)定為胺類＞羧酸類＞多聚物類＞碳水化合物類＞氨基酸類＞酚酸類；192 h后WS2水樣穩(wěn)定為多聚物類＞氨基酸類＞酚酸類＞碳水化合物類＞羧酸類＞胺類；WS3、WS4和WS5水樣6類碳源的利用順序相同，分別在144 h、168 h和120 h穩(wěn)定為多聚物類＞碳水化合物類＞羧酸類＞氨基酸類＞酚酸類＞羥酸類胺類.即除WS1水樣外其他4種水樣6類碳源的利用達(dá)到?jīng)]有顯著性差異之后，6類碳源的利用順序都趨于穩(wěn)定.

圖2 不同培養(yǎng)時(shí)間下微生物群落的6類碳源利用情況Fig.2 Utilization of six grouped carbon sources by microbial communities at different incubation time

2.3 不同培養(yǎng)時(shí)間多樣性指數(shù)的變化

5種水樣72～192 h微生物群落的多樣性指數(shù)如表2所示.由表2可以看出，WS1水樣微生物群落的Shannon-Wiener指數(shù)、Pielou均勻度指數(shù)、Simpson指數(shù)、McIntosh指數(shù)和McIntosh均勻度分別在144、72、120、72、72 h后沒有顯著差異，即5個(gè)多樣性指數(shù)在144 h后沒有顯著差異；WS2水樣微生物群落的各項(xiàng)多樣性指數(shù)分別在144、96、120、96、96 h后沒有顯著差異，即5個(gè)多樣性指數(shù)在144 h后沒有顯著差異；WS3水樣微生物群落5個(gè)多樣性指數(shù)均在120 h后沒有顯著差異；WS4水樣微生物群落的5個(gè)多樣性指數(shù)分別在120、72、120、96、72 h后沒有顯著差異，即5個(gè)多樣性指數(shù)在120 h后沒有顯著差異；WS5水樣微生物群落除McIntosh指數(shù)在120 h后沒有顯著差異外，其他各項(xiàng)指數(shù)均在72 h后沒有顯著差異.

表2 不同培養(yǎng)時(shí)間下微生物群落的多樣性指數(shù)Tab.2 Diversity indices for microbial communities at different incubation time

2.4 不同培養(yǎng)時(shí)間水樣主成分的變化

對(duì)5種水樣不同培養(yǎng)時(shí)間的數(shù)據(jù)進(jìn)行主成分分析，分析不同培養(yǎng)時(shí)間微生物群落底物利用特征的差異，結(jié)果如圖3所示.圖3中展示了每一水樣6個(gè)平行樣從72～192 h各時(shí)間點(diǎn)的排序情況，從圖3中可以看出每個(gè)樣品72 h、96 h的各點(diǎn)分布與120、144、168和192 h的各點(diǎn)分布距離較遠(yuǎn)，反映出樣品微生物群落培養(yǎng)前期CSUP與培養(yǎng)后期差別較大.通過ANOSIM方法分析5種水樣微生物群落不同時(shí)間點(diǎn)CSUP的差異顯著性，結(jié)果發(fā)現(xiàn)5種水樣微生物群落的CSUP分別在96、120、144、120和96 h后沒有顯著差異，各水樣ANOSIM分析的R值和P值如表3所示.

圖3 水樣不同培養(yǎng)時(shí)間下微生物群落底物利用特征的主成分分析Fig.3 Community substrate utilization profile at different incubationtimebyprincipalcomponentanalysis（PCA）

表3 不同培養(yǎng)時(shí)間下微生物群落底物利用特征的ANOSIM分析Tab.3 ANOSIM analysis between CSUPs of microbial communities at different incubation time

2.5 不同時(shí)間點(diǎn)對(duì)5種水樣的區(qū)分

為了解不同時(shí)間點(diǎn)對(duì)5種水樣的區(qū)分情況，利用5種水樣在72～192 h的吸光度值數(shù)據(jù)進(jìn)行主成分分析，結(jié)果如圖4所示.由圖4可以看出，72、96和120 h時(shí)間點(diǎn)5種水樣在第一和第二主成分上的距離和分布可見明顯變化，144、168和192 h時(shí)間點(diǎn)5種水樣在第一和第二主成分上的分布基本一致，說明144 h后微生物群落的底物利用特征趨于穩(wěn)定.ANOSIM分析的R值和P值如表4所示.通過ANOSIM方法分析發(fā)現(xiàn)各時(shí)間點(diǎn)5種水樣微生物群的CSUP均存在顯著差異（P＜0.05），即各時(shí)間點(diǎn)數(shù)據(jù)均能區(qū)分5種水樣.對(duì)不同時(shí)間點(diǎn)5種水樣微生物群落的CSUP進(jìn)行SIMPER相似度分析，結(jié)果如表5所示.由表5可以看出，在72～192 h內(nèi)隨著培養(yǎng)時(shí)間的延長(zhǎng)兩兩水樣之間的差異百分比逐步降低.

圖4 不同培養(yǎng)時(shí)間對(duì)5種水樣的區(qū)分Fig.4 Differentiation of the five water samples by different incubation time

表4 不同培養(yǎng)時(shí)間5種水樣兩兩比較的ANOSIM分析Tab.4 ANOSIM analysis of paired comparison between the five water samples at different incubation time

表5 不同培養(yǎng)時(shí)間5種水樣兩兩比較的SIMPER分析Tab.5 SIMPER analysis of paired comparison between the five water samples at different incubation time%

3 討論

在研究微生物多樣性時(shí)，由于環(huán)境樣品的復(fù)雜性，在利用Biolog Eco板分析時(shí)要選擇眾多參數(shù)包括AWCD值、6類碳源利用情況、多樣性指數(shù)及主成分分析等進(jìn)行描述和比較不同樣品間的差異[30-31].本研究發(fā)現(xiàn)養(yǎng)殖水體微生物群落的上述指標(biāo)達(dá)到?jīng)]有顯著差異的時(shí)間不同，其中WS1水樣AWCD值、6類碳源利用情況、多樣性指數(shù)及主成分分析達(dá)到?jīng)]有顯著差異的時(shí)間點(diǎn)分別為96、96、144和96 h；WS2水樣為168、192、144和120 h；WS3水樣為168、168、120和144 h；WS4水樣為144、168、120和120 h；WS5水樣為120、120、120和96 h.這是因?yàn)橥粯悠犯髦笜?biāo)達(dá)到穩(wěn)定的時(shí)間不同與各指標(biāo)對(duì)時(shí)間的敏感度不同相關(guān).AWCD作為顏色平均變化率反映的是Biolog Eco板的平均活性，AWCD值穩(wěn)定后，各孔的顏色還可能不斷變化，因此各類碳源利用情況和多樣性指數(shù)還會(huì)發(fā)生變化.

本實(shí)驗(yàn)中養(yǎng)殖水體微生物群落各功能指標(biāo)達(dá)到穩(wěn)定的時(shí)間不同，但是各時(shí)間點(diǎn)都在AWCD曲線的達(dá)到?jīng)]有顯著差異的時(shí)間點(diǎn)（即拐點(diǎn)）附近.其中除WS1水樣多樣性指數(shù)在拐點(diǎn)之后48 h，WS2和WS4水樣6類碳源的利用在拐點(diǎn)之后的24 h外，5種水樣其他各項(xiàng)指標(biāo)達(dá)到穩(wěn)定的時(shí)間都在拐點(diǎn)之前.說明當(dāng)水樣AWCD曲線達(dá)到拐點(diǎn)時(shí)，微生物群落的各項(xiàng)代謝功能指標(biāo)都趨于穩(wěn)定.這與賈夏等[19]研究土壤樣品Biolog Eco板培養(yǎng)時(shí)間對(duì)分析結(jié)果的影響時(shí)的結(jié)論相一致.

本研究同時(shí)討論了不同時(shí)間點(diǎn)對(duì)5種水樣的區(qū)分，ANOSIM分析發(fā)現(xiàn)72～192 h各個(gè)時(shí)間點(diǎn)都可以很好地區(qū)分5種水樣，但是SIMPER分析表明水樣之間的差異程度隨培養(yǎng)時(shí)間的延長(zhǎng)而降低.可能因?yàn)殡S著培養(yǎng)時(shí)間的延長(zhǎng)，Biolog Eco板上很多孔的顏色處于飽和，只能通過顯色較小或者沒有顏色變化的孔進(jìn)行樣品的區(qū)分，所以導(dǎo)致水樣之間差異程度減小.如果要充分研究水樣微生物群落CLPP各相關(guān)指標(biāo)，可以利用AWCD值達(dá)到?jīng)]有顯著差異的時(shí)間點(diǎn)進(jìn)行分析，這時(shí)可以兼顧生長(zhǎng)較慢的微生物，能夠全面地反映環(huán)境中的微生物群落[19].但建議不過多延長(zhǎng)培養(yǎng)時(shí)間，因?yàn)殡S著培養(yǎng)時(shí)間的延長(zhǎng)存在平行樣品標(biāo)準(zhǔn)差增大和光孔反應(yīng)飽和的問題[13].

如果為最大程度地區(qū)分樣品微生物群落CLPP各項(xiàng)指標(biāo)的差異，可以選擇AWCD值指數(shù)增長(zhǎng)期的時(shí)間點(diǎn)進(jìn)行分析，此時(shí)可以較好地反映出環(huán)境樣品中生長(zhǎng)較快的細(xì)菌種類.但是利用指數(shù)增長(zhǎng)期數(shù)據(jù)進(jìn)行分析，一般要求大致相同的接種密度[1].本實(shí)驗(yàn)中WS1水樣在接入Biolog Eco板后快速進(jìn)入指數(shù)增長(zhǎng)期，而WS2～WS5水樣的指數(shù)增長(zhǎng)較緩慢.目前養(yǎng)殖水體富營(yíng)養(yǎng)化問題突出，池塘水樣成分復(fù)雜，養(yǎng)殖水體理化指標(biāo)是否會(huì)對(duì)Biolog Eco板的顯色產(chǎn)生影響還有待于進(jìn)一步研究.

綜上所述，養(yǎng)殖池塘水樣接入Biolog Eco板的培養(yǎng)時(shí)間對(duì)分析微生物群落的功能多樣性具有顯著影響，微生物群落的各項(xiàng)功能指標(biāo)均在AWCD值達(dá)到?jīng)]有顯著差異的時(shí)間點(diǎn)趨于穩(wěn)定.為了更好地區(qū)分樣品，可以選擇指數(shù)增長(zhǎng)期的時(shí)間點(diǎn)，此時(shí)水樣之間的差異度高于穩(wěn)定期.