朱小奕，陳瑾，張建英,*

1.浙江大學(xué)環(huán)境科學(xué)研究所，杭州310058

2.浙江省有機污染過程與控制重點實驗室，杭州310058

水環(huán)境中藻毒素生態(tài)風(fēng)險的物種敏感性分布評價

朱小奕1,2，陳瑾1,2，張建英1,2,*

1.浙江大學(xué)環(huán)境科學(xué)研究所，杭州310058

2.浙江省有機污染過程與控制重點實驗室，杭州310058

藻毒素對人體的健康風(fēng)險已受普遍關(guān)注，然而其對水生態(tài)物種敏感性分布的影響尚不明朗。本研究采集已有實驗數(shù)據(jù)，利用種間相關(guān)性分析(ICE)模型和物種敏感性分布評估(SSDs)方法，篩選了64個水生生物物種的71組急性毒性數(shù)據(jù)(EC50)，構(gòu)建水生生物對肝毒素(節(jié)球藻毒素、柱孢藻毒素)和神經(jīng)毒素(類毒素、貝毒素)的SSD方程。在此基礎(chǔ)上，計算不同暴露濃度下的潛在影響比例(PAF)以及保護95%物種基礎(chǔ)上藻毒素對水生生物的生態(tài)風(fēng)險閾值(HC5)，比較不同類別生物對藻毒素的敏感性以及藻毒素對水生生物的生態(tài)風(fēng)險。結(jié)果表明：(1)在95%物種保護保證率下，節(jié)球藻毒素、柱孢藻毒素、類毒素和貝毒素對全部物種的HC5值分別為74.96、205.39、194.39、0.3 μg·L-1，貝毒素水生態(tài)風(fēng)險最高，柱孢藻毒素最低。(2)柱孢藻毒素、類毒素對無脊椎動物的HC5值分別為122.93、95.19 μg·L-1，低于全部物種的HC5值，無脊椎動物受柱孢藻毒素、類毒素影響較其他物種大。(3)物種潛在影響比例可明確表征敏感性，柱孢藻毒素、類毒素在各濃度暴露情景下對無脊椎動物的PAF值均高于脊椎動物，顯示無脊椎動物的敏感性較高。

藻毒素；水生生物；物種敏感性分布；生態(tài)閾值

朱小奕,陳瑾,張建英.水環(huán)境中藻毒素生態(tài)風(fēng)險的物種敏感性分布評價[J].生態(tài)毒理學(xué)報，2016,11(3):131-139

Zhu X Y,Chen J,Zhang J Y.Assessing ecological risk of cyanotoxins based on interspecies correlation estimation and species sensitivity distributions[J]. Asian Journal of Ecotoxicology,2016,11(3):131-139(in Chinese)

湖泊藻華爆發(fā)遍布全球，其中25%～70%的藻華會產(chǎn)生有毒的二次代謝產(chǎn)物藻毒素。根據(jù)毒素對動物的致毒方式，可分為肝毒素和神經(jīng)毒素等，肝毒素包括微囊藻毒素(microcystins,MCs)、節(jié)球藻毒素(nodularins, NODs)和柱孢藻毒素(cylindrospermopsin,CYN)，神經(jīng)毒素包括類毒素(anatoxin-a,antx-a)、貝毒素(saxitoxins, STXs)等[1-2]。1990年至2014年間，中國已有60%的淡水湖泊出現(xiàn)不同程度的富營養(yǎng)化現(xiàn)象，引發(fā)嚴(yán)重藻華污染，釋放藻毒素以微囊藻毒素為主，濃度可達35.3 μg ·L-1[3]，近年柱孢藻毒素、類毒素和貝毒素等污染也日趨嚴(yán)重，如柱孢藻毒素濃度達8.25 μg·L-1[4]。藻毒素不但會引起肝、神經(jīng)毒性等健康風(fēng)險，還可導(dǎo)致魚類的死亡、浮游動物與水生植物生長和繁殖抑制等物種危害[5-6]，從而造成水生態(tài)物種多樣性破壞風(fēng)險[7-8]。然而，多種藻毒素在生態(tài)系統(tǒng)尺度上的風(fēng)險尚不清楚，其生態(tài)風(fēng)險閾值亟待確定。

物種敏感性分布評估方法(species sensitivity distributions,SSDs)是基于不同物種對某一環(huán)境脅迫的敏感度服從一定累積概率分布假設(shè)，以統(tǒng)計分布模型來描述不同物種樣本對脅迫因素的敏感性差異，實現(xiàn)將單一物種的測試結(jié)果外推至生態(tài)系統(tǒng)的風(fēng)險評估方法[9]。該方法被廣泛應(yīng)用于生態(tài)風(fēng)險評價[10-11]和水質(zhì)基準(zhǔn)推導(dǎo)[12-13]，目前國內(nèi)學(xué)者也嘗試運用SSD方法對農(nóng)藥[14-15]、重金屬[16]、多環(huán)芳烴[17]、微囊藻毒素[18]等進行生態(tài)風(fēng)險評估。SSD擬合物種數(shù)要求不少于5個[19-20]，為彌補不同物種，特別是稀缺物種毒性數(shù)據(jù)缺失對SSD方法結(jié)果不確定性的影響，美國EPA提出種間相關(guān)性分析(interspecies correlation estimations, ICE)模型。在ICE模型中，假定替代物種和預(yù)測物種的毒性值符合最小二乘法線性回歸，因此可通過替代物種毒性數(shù)據(jù)，外推預(yù)測物種的毒性效應(yīng)[21]。ICE模型彌補了SSD評估中毒性數(shù)據(jù)的斷層，已被用于水生無脊椎動物和魚類的毒性外推研究，如ICE得出硝基芳烴對不同生物毒性數(shù)據(jù)與其實驗數(shù)據(jù)擬合結(jié)果在1%的置信水平上呈顯著的線性正相關(guān)[22-23]；對鋅水質(zhì)基準(zhǔn)研究顯示基于ICE和實驗數(shù)據(jù)得出的SSD估算HC5值無顯著差異[24]，表明應(yīng)用ICE方法預(yù)測缺失毒性數(shù)據(jù)用于水生態(tài)風(fēng)險評估是可行的。

目前，對微囊藻毒素的生態(tài)風(fēng)險研究已從生物個體水平發(fā)展到到生態(tài)系統(tǒng)水平，但除此之外，水環(huán)境中節(jié)球藻毒素、柱孢藻毒素、類毒素和貝毒素等的研究僅限于生物個體毒性，尚未報道其在系統(tǒng)尺度上的生態(tài)風(fēng)險。本研究通過挖掘文獻中4種藻毒素的急性毒性試驗數(shù)據(jù)，構(gòu)建ICE-SSD復(fù)合模型，評價不同濃度水平藻毒素污染的水生態(tài)物種敏感性風(fēng)險，計算藻毒素對水生生物的風(fēng)險閾值。擬提供在系統(tǒng)尺度上評價藻毒素污染風(fēng)險的一個保證率為95%的可靠評價模型，從而為藻毒素污染水體的生態(tài)風(fēng)險管理提供方法和依據(jù)。

1 材料與方法（Materials and methods）

1.1 毒性數(shù)據(jù)獲取與物種分組

表1 毒性數(shù)據(jù)篩選條件Table 1 Criteria for selecting toxicity data

分別利用節(jié)球藻毒素、柱孢藻毒素、類毒素與貝毒素的EC50急性毒性數(shù)據(jù)構(gòu)建SSD曲線，用以評價其水生態(tài)風(fēng)險，數(shù)據(jù)來源于Web of Science(http:// www.webofknowledge.com/)收錄的相關(guān)文獻資料，其篩選原則參照歐盟生態(tài)風(fēng)險評估標(biāo)準(zhǔn)“風(fēng)險評估技術(shù)綱領(lǐng)”等[25]：(1)標(biāo)準(zhǔn)毒性試驗方法，如經(jīng)濟合作與發(fā)展組織(Organization for Economic Cooperation and Development,OECD)或美國材料與試驗協(xié)會(A-merican Society for Testing and Materials,ASMT)等發(fā)布方法；(2)實驗條件的設(shè)定，如暴露周期、介質(zhì)等，具體篩選條件見表1[9]，本研究中節(jié)球藻毒素、柱孢藻毒素和類毒素以淡水分布為主，選取其對淡水生物毒性試驗結(jié)果，而貝毒素主要來源于海洋藻釋放，評價其對海洋生物的急性毒性；(3)急性毒性終點選取，以半數(shù)效應(yīng)濃度EC50表征，其效應(yīng)包括運動抑制、生長抑制、胚胎發(fā)育抑制和致死率等中毒癥狀；(4)評價物種數(shù)不少于5個，且盡量分布于3個以上營養(yǎng)級。為彌補柱孢藻毒素、類毒性、貝毒素適當(dāng)毒性數(shù)據(jù)的缺失，本文進一步采用EPA的ICE模型(http:// www.epa.gov/ceampubl/fchain/webice/)，通過選取2個具有代表性的物種(大型蚤Daphnia magna,Daphnia pulex)為替代物種外推同類生物急性毒性數(shù)據(jù)。

根據(jù)文獻和ICE模型預(yù)測，共獲得64個水生生物物種的71組急性毒性數(shù)據(jù)，包含36個甲殼類、17個魚類、6個昆蟲類、2個軟體類、2個棘皮類和1個原生動物。節(jié)球藻毒素、柱孢藻毒素、類毒素與貝毒素的毒性數(shù)據(jù)分別含7、24、24和9個物種。同一物種擁有多個毒性數(shù)據(jù)時，以幾何平均值作為該物種的數(shù)據(jù)點[26]。根據(jù)獲得物種毒性數(shù)據(jù)分布，柱孢藻毒素和類毒素在對全部物種進行分析的基礎(chǔ)上，進一步細分為脊椎動物和無脊椎動物以比較不同類別生物受到污染物危害的風(fēng)險大小。全部物種包括脊椎動物和無脊椎動物；脊椎動物包括魚類和兩棲類；無脊椎動物包括甲殼類、昆蟲類、軟體動物、蠕蟲類及其他無脊椎動物。

1.2 SSD擬合

將生物的毒性數(shù)據(jù)濃度值(或其對數(shù)值)對以濃度排列的分位數(shù)作圖，并對其進行參數(shù)擬合得到SSD物種敏感性分布曲線。對比不同參數(shù)擬合方法，BurrIII模型分布靈活，可根據(jù)參數(shù)值的大小轉(zhuǎn)換成ReWeibull和RePareto兩種模式，且數(shù)據(jù)量要求較少(最小數(shù)據(jù)量為7個)，對本研究篩選獲得的數(shù)據(jù)擬合度較好，除貝毒素對全部物種的SSD曲線擬合R2為0.89，其余曲線擬合R2均大于0.93。因此文中采用澳大利亞聯(lián)邦科學(xué)和工業(yè)研究組織CSIRO提供的BurrlizO(1.0.14)軟件計算BurrIII型分布模式[18]，其方程為：

式中：x為污染物濃度值(μg·L-1)；b、c、k為函數(shù)的3個參數(shù)，其中b為尺度參數(shù)，c為形狀參數(shù)，曲線位置則由c、k共同決定。

當(dāng)k趨于無窮大時，BurrIII分布變化為Re-Weibull分布：

當(dāng)c趨于無窮大時，變化為RePareto分布：

2 結(jié)果（Results）

2.1 水生生物對藻毒素的敏感響應(yīng)

不同類別生物對節(jié)球藻毒素、柱孢藻毒素、類毒素、貝毒素的急性毒性響應(yīng)值見表2和表3。水生生物對貝毒素、節(jié)球藻毒素、柱孢藻毒素、類毒素的對數(shù)EC50值分別為1.81、3.11、3.48、3.59 μg·L-1，敏感響應(yīng)濃度依次升高，顯示貝毒素的毒性最大。柱孢藻毒素、類毒素對無脊椎動物的對數(shù)EC50值均低于脊椎動物，表明這2種藻毒素對無脊椎動物的毒性更大。

表2 全部物種的毒理數(shù)據(jù)統(tǒng)計特征（對數(shù)變換后）Table 2 Log-transformed toxicity data statistics of all species

表3 脊椎動物和無脊椎動物的毒理數(shù)據(jù)統(tǒng)計特征（對數(shù)變換后）Table 3 Log-transformed toxicity data statistics of vertebrates and invertebrates

表4 利用BurrlizO計算SSD參數(shù)結(jié)果Table 4 SSD parameters calculated by BurrlizO

表5 不同暴露濃度下藻毒素的PAF預(yù)測值（單位：%）Table 5 Predicted PAF values under different concentrations of cyanotoxins(unit:%)

2.2 藻毒素水生態(tài)風(fēng)險的物種敏感性分布規(guī)律

表4是使用BurrlizO軟件計算得出的柱孢藻毒素、類毒素與節(jié)球藻毒素、貝毒素SSD方程參數(shù)結(jié)果。其中柱孢藻毒素、類毒素對全部物種和無脊椎物種的敏感性分布曲線符合Burr III分布，對脊椎動物符合ReWeibull分布；節(jié)球藻毒素對全部物種適用Burr III分布擬合；貝毒素對全部物種適用Re-Pareto分布擬合。

2.3 不同暴露濃度下藻毒素的PAF值

不同濃度暴露情景下，柱孢藻毒素、類毒素與節(jié)球藻毒素、貝毒素對水生生物的PAF預(yù)測值見表5。表5顯示了由不同濃度值得到的PAF值的大小，反映其對不同物種的損害程度。在1 000 μg·L-1下，柱孢藻毒素和類毒素對全部物種和無脊椎動物的損害比例為20.3%～31.6%，超過了5%的生態(tài)風(fēng)險閾值，對水生態(tài)系統(tǒng)造成一定危害；在各濃度暴露情景下，柱孢藻毒素、類毒素對無脊椎動物的PAF值均高于脊椎動物，顯示無脊椎動物的敏感性較高。節(jié)球藻毒素濃度達到1 000 μg·L-1時，對全部物種的PAF值為38.9%，危害明顯。貝毒素對生態(tài)系統(tǒng)的影響最為顯著，當(dāng)貝毒素濃度為100 μg·L-1時，對全部物種的PAF值即高達43.2%。

2.4 藻毒素對水生生物的HC5值及生態(tài)風(fēng)險比較

柱孢藻毒素、類毒素與節(jié)球藻毒素、貝毒素對不同類別生物的HC5值見表6，與文獻[18]中微囊藻毒素急性毒性數(shù)據(jù)HC5值(19.22 μg·L-1)相比較可得，5種藻毒素對全部物種的HC5值從小到大依次為：貝毒素<微囊藻毒素<節(jié)球藻毒素<類毒素<柱孢藻毒素，表現(xiàn)出依次降低的生態(tài)風(fēng)險，貝毒素水生態(tài)風(fēng)險最高，柱孢藻毒素最低。柱孢藻毒素、類毒素對無脊椎動物的HC5值(122.93、95.19 μg·L-1)均低于脊椎動物(1 134.77、1 662.73 μg·L-1)，表明無脊椎動物受到損害的風(fēng)險更高。

HC5值反映低濃度下不同類別生物對藻毒素的耐受性閾值，而通過SSD曲線可以直觀比較不同濃度范圍內(nèi)物種對污染物的敏感性差異。對4種藻毒素的毒性數(shù)據(jù)構(gòu)建SSD曲線，并與文獻[18]中對微囊藻毒素構(gòu)建的SSD曲線相比較，全部物種對5種藻毒素的SSD曲線如圖1所示。貝毒素曲線較陡峭，在低濃度時即對水生生物產(chǎn)生較大的生態(tài)風(fēng)險，并且隨著濃度增加，物種敏感性增大的趨勢迅速，在濃度高于932.60 μg·L-1時，累計概率即達到100%，生態(tài)風(fēng)險極大。節(jié)球藻毒素、柱孢藻毒素和類毒素的SSD曲線斜率差距較小，從曲線整體位置來看，低濃度時從左到右依次為：節(jié)球藻毒素、類毒素、柱孢藻毒素，生態(tài)風(fēng)險依次降低，與比較HC5值得到的結(jié)果一致；而隨著濃度的升高，柱孢藻毒素對水生生物的生態(tài)風(fēng)險逐漸大于類毒素。微囊藻毒素的SSD曲線較緩，在低濃度時，生物敏感性高于節(jié)球藻毒素、柱孢藻毒素和類毒素，但隨著濃度增加，敏感性增大的趨勢較緩慢，在高濃度時，微囊藻毒素對水生生物的生態(tài)風(fēng)險逐漸低于節(jié)球藻毒素、柱孢藻毒素和類毒素。

圖1 全部物種對5種藻毒素的物種敏感性分布（SSD）曲線圖

不同類別生物(I脊柱動物，II無脊椎動物)對微囊藻毒素、柱孢藻毒素和類毒素的SSD曲線如圖2，其中，微囊藻毒素的SSD曲線借鑒于文獻[18]。圖I中，脊椎動物對柱孢藻毒素和類毒素的敏感性差異不顯著，柱孢藻毒素對脊椎動物的生態(tài)風(fēng)險略大于類毒素；脊椎動物對微囊藻毒素的SSD曲線位于最左側(cè)，表明微囊藻毒素對脊椎動物的生態(tài)風(fēng)險大于柱孢藻毒素和類毒素。從圖II中可見，無脊椎動物對柱孢藻毒素和類毒素的敏感性差異同樣不顯著，依然是柱孢藻毒素的生態(tài)風(fēng)險略大；無脊椎動物對微囊藻毒素的SSD曲線變化較緩，在低濃度時，無脊椎動物對微囊藻毒素的敏感性高于柱孢藻毒素和類毒素，但隨著濃度的增加，敏感性逐漸低于柱孢藻毒素和類毒素。

對比脊椎動物、無脊椎動物對不同藻毒素(柱孢藻毒素，類毒素)的SSD曲線(圖3)可看出，柱孢藻毒素與類毒素的情況相似，在低濃度時，無脊椎動物與脊椎動物的敏感性差異較顯著，無脊椎動物的敏感性高于脊椎動物；隨著濃度增大，脊椎動物的敏感性增大趨勢更迅速，因此在到達高濃度時，無脊椎動物與脊椎動物的敏感性差異較小。

表6 藻毒素對水生生物的HC5值（單位：μg·L-1）Table 6 HC5values of cyanotoxins to aquatic organisms(unit:μg·L-1)

2.5 自然水體中藻毒素生態(tài)風(fēng)險評估

運用本文擬合的SSD曲線對我國典型水體藻毒素污染水平進行生態(tài)風(fēng)險評估[4,27]，其物種損害比例見表7。在廣東省東莞市水濂山水庫和河北省秦皇島市洋河水庫，分別檢測到柱孢藻毒素濃度和類毒素濃度為8.25 μg·L-1、0.29 μg·L-1，均低于本研究HC5值，對全部物種的影響比例為0.1%和0.01%，顯示水生態(tài)風(fēng)險不顯著。節(jié)球藻毒素和貝毒素在水體中的濃度水平較低，但在我國東海和南海海域的貝類中毒素的檢出率分別為14%和30%[28]，表明水生動物以藻類為食亦可造成藻毒素體內(nèi)富集和產(chǎn)生毒害效應(yīng)，進而威脅水生態(tài)安全。

3 討論（Discussion）

圖2 脊柱動物（I）、無脊椎動物（II）對微囊藻毒素、柱孢藻毒素和類毒素的SSD曲線比較

圖3 脊椎動物、無脊椎動物對柱孢藻毒素（I）、類毒素（II）的SSD曲線比較

表7 國內(nèi)水體柱孢藻毒素和類毒素對淡水生物的PAF預(yù)測值（單位：%）Table 7 Predicted PAF values of cylindrospermopsin and anatoxin-a to freshwater species in studied waters(unit:%)

通過對4種藻毒素水生態(tài)風(fēng)險的評價，可見不同藻毒素的水生態(tài)風(fēng)險差異顯著，同為神經(jīng)毒素的類毒素與貝毒素對全部物種的HC5值相差近650倍，可歸因于其致毒分子機制的差異。研究表明，類毒素主要與乙酰膽堿存在煙堿受體的競爭性結(jié)合作用，在對虹鱒(Oncorhynchus mykiss)魚苗的暴露實驗中，魚體肌肉組織的乙酰膽堿酯酶(AChE)活性隨著類毒素劑量的增加而顯著增強[29]；而麻痹性貝類毒素(paralytic shellfish toxins,PSTs)之一的貝毒素則通過阻止神經(jīng)細胞鈉離子通道抑制神經(jīng)沖動的正常動作電位傳輸，從而引起神經(jīng)肌肉的麻痹，被認(rèn)為是最有效的和快速的神經(jīng)毒素之一[30]。此外，同種藻毒素對不同種類生物毒性影響亦差異顯著，無脊椎動物對柱孢藻毒素與類毒素暴露的敏感性均明顯高于脊椎動物?？赡芘c其對藻毒素降解機制的差異有關(guān)，谷胱甘肽轉(zhuǎn)移酶(glutathione s-transferase,GST)通過催化內(nèi)源性谷胱甘肽(glutathione,GSH)與毒素結(jié)合而起到降毒作用，是生物體毒素降解的重要指標(biāo)之一。10 μg·L-1柱孢藻毒素對羅非魚(Oreochromis niloticus)暴露7 d，魚體內(nèi)GST酶活性增強，誘導(dǎo)率約為35%[31]；而648.6 ng·L-1柱孢藻毒素對大型蚤(Daphnia magna)暴露24 h，蚤體內(nèi)GST酶活性即被誘導(dǎo)，誘導(dǎo)率達25.14%[32]。顯示無脊椎動物大型蚤(Daphnia magna)在低濃度藻毒素暴露下短時間內(nèi)即產(chǎn)生應(yīng)急毒性響應(yīng)，敏感性更高。

物種敏感性分布是一種基于單物種測試的外推方法，它把可能發(fā)生的風(fēng)險依靠統(tǒng)計模型以概率的方式表達出來。影響SSD評估不確定性的因素包括毒性數(shù)據(jù)的選擇和統(tǒng)計模型的選擇[33]。目前的研究結(jié)果表明數(shù)據(jù)的選擇比統(tǒng)計模型的選擇對HC5值影響更大[9]。SSD方法對毒理數(shù)據(jù)的選擇暫沒有統(tǒng)一規(guī)定，經(jīng)濟合作與發(fā)展組織(OECD)推薦的最小數(shù)據(jù)量為5個，美國環(huán)保署(US EPA)建議至少涉及3個營養(yǎng)級8個分類群。本研究中柱孢藻毒素和類毒素的毒理數(shù)據(jù)量足夠覆蓋3個營養(yǎng)級，而節(jié)球藻毒素和貝毒素的毒理數(shù)據(jù)量不足3個營養(yǎng)級，但超過OECD推薦的5個數(shù)據(jù)點限值。此外，物種敏感性分布假定有限的生物物種是從整個生態(tài)系統(tǒng)中隨機取樣的，故認(rèn)為評估有限物種的損害概率是適合整個生態(tài)系統(tǒng)的。但在實際應(yīng)用中，物種樣本的數(shù)據(jù)是根據(jù)數(shù)據(jù)的有效性來采用的，而不是隨機樣本，因此在很大程度上增加了不確定性。因此針對毒理數(shù)據(jù)缺乏帶來的風(fēng)險評估不確定性，需進一步開展水生受試生物的毒理學(xué)實驗，豐富毒理學(xué)研究數(shù)據(jù)，提高采用SSD方法評估的準(zhǔn)確性與可靠性。

綜上所述，藻毒素因致毒機制等差異可對水生生物造成不同程度毒性影響，有必要對多種藻毒素進行水生態(tài)風(fēng)險評估，以獲得其生態(tài)安全濃度閾值。本研究基于不同營養(yǎng)級水生生物對肝毒性和神經(jīng)毒性藻毒素的急性毒性數(shù)據(jù)運用SSD概率統(tǒng)計分布方法評估4種藻毒素風(fēng)險，指出其生態(tài)風(fēng)險閾值HC5為0.3～205.39 μg·L-1不等，尤以神經(jīng)毒素貝毒素的水生態(tài)風(fēng)險最高，為藻毒素污染的水生態(tài)風(fēng)險管理提供依據(jù)。但在SSD模型評估中，毒性數(shù)據(jù)缺失及不同地區(qū)敏感物種組成的差異性，導(dǎo)致風(fēng)險評估不確定性增加，有必要開展進一步的毒理學(xué)實驗或野外研究，獲得多營養(yǎng)級更多分類群的毒理數(shù)據(jù)，進一步修正本文的相關(guān)研究結(jié)果。

（References）：

[1] 謝平.微囊藻毒素對人類健康影響相關(guān)研究的回顧[J].湖泊科學(xué),2009,21(5):603-613

Xie P.A review on the studies related to the effects of microcystins on human health[J].Journal of Lake Sciences, 2009,21(5):603-613(in Chinese)

[2] Merel S,Walker D,Chicana R,et al.State of knowledge and concerns on cyanobacterial blooms and cyanotoxins [J].Environment International,2013,59:303-327

[3] Otten T G,Xu H,Qin B,et al.Spatiotemporal patterns and ecophysiology of toxigenic microcystis blooms in Lake Taihu,China:Implications for water quality management[J].Environmental Science&Technology,2012, 46(6):3480-3488

[4] Lei L,Peng L,Huang X H,et al.Occurrence and dominance ofCylindrospermopsis raciborskiiand dissolved cylindrospermopsin in urban reservoirs used for drinking water supply,South China[J].Environmental Monitoring and Assessment,2014,186(5):3079-3090

[5] Ferrao-Filho A S,Kozlowsky-Suzuki B.Cyanotoxins: Bioaccumulation and effects on aquatic animals[J].Marine Drugs,2011,9(12):2729-2772

[6] Saqrane S,Ghazali I E,Ouahid Y,et al.Phytotoxic effects of cyanobacteria extract on the aquatic plantLemna gibba: Microcystin accumulation,detoxicationandoxidative stress induction[J].Aquatic Toxicology,2007,83(4):284-294

[7] Bittencourt-Oliveira M C,Hereman T C,Cordeiro-Araujo M K,et al.Phytotoxicity associated to microcystins:A review[J].Brazilian Journal of Biology,2014,74(4):753-760

[8] Landsberg J H.The effects of harmful algal blooms on aquatic organisms[J].Reviews in Fisheries Science,2002, 10(2):113-390

[9] Wheeled J R,Grist E P M,Leung K M Y,et al.Speciessensitivity distributions:Data and model choice[J].Marine Pollution Bulletin,2002,45:192-202

[10] Raimondo S,Vivian D N,Delos C,et al.Protectiveness of species sensitivity distribution hazard concentration for acute toxicity used in endangered species risk assessment [J].Environmental Toxicology and Chemistry,2008,27 (12):2599-2607

[11] Garner K L,Suh S,Lenihan H S,et al.Species sensitivity distributions for engineered nanomaterials[J].Environmental Science&Technology,2015,49(9):5753-5759

[12] Australian and New Zealand Environment and Conservation Council and Agriculture and Resource Management Council of Australia and New Zealand.Australian and New Zealand Guidelines for Fresh and Marine Water Quality[R].Canberra:ANZECC and ARMCANZ,2000

[13] National Institute for Public Health and the Environment. Guidance for the Derivation of Environmental Risk Limits within the Framework of‘International and National Environmental Quality Standards for Substances in the Netherlands’(INS)[R].BA Bilthoven,Netherlands:RIVM, 2007

[14] 王印,王軍軍,秦寧,等.應(yīng)用物種敏感性分布評估DDT和林丹對淡水生物的生態(tài)風(fēng)險[J].環(huán)境科學(xué)學(xué)報,2009,29(11):2407-2414

Wang Y,Wang J J,Qin N,et al.Assessing ecological risks of DDT and lindane to fresh water organisms by species sensitivity distributions[J].Acta Scientiae Circumstantiae,2009,29(11):2407-2414(in Chinese)

[15] 徐瑞祥,陳亞華.應(yīng)用物種敏感性分布評估有機磷農(nóng)藥對淡水生物的急性生態(tài)風(fēng)險[J].湖泊科學(xué),2012,24 (6):811-821 Xu R X,Chen Y H.Assessing acute ecological risks of organophosphorus pesticides to freshwater organisms by species sensitivity distributions[J].Journal of Lake Sciences,2012,24(6):811-821(in Chinese)

[16] 孔祥臻,何偉,秦寧,等.應(yīng)用物種敏感性分布評估重金屬對淡水生物的急性生態(tài)風(fēng)險[J].中國環(huán)境科學(xué), 2011,31(9):1555-1562

Kong X Z,He W,Qin N,et al.Assessing acute ecological risks of heavy metals to freshwater organisms by species sensitivity distributions[J].China Environmental Science, 2011,31(9):1555-1562(in Chinese)

[17] 蔣丹烈,胡霞林,尹大強.應(yīng)用物種敏感性分布法對太湖沉積物中多環(huán)芳烴的生態(tài)風(fēng)險分析[J].生態(tài)毒理學(xué)報,2011,6(1):60-66

Jiang D L,Hu X L,Yin D Q.Ecological risk assessment on polycyclic aromatic hydrocarbons of sediment in Taihu Lake using species sensitivity distributions[J].Asian Journal of Ecotoxicology,2011,6(1):60-66(in Chinese)

[18] 陳瑾,劉奕梅,張建英.基于物種敏感性分布的微囊藻毒素與氮污染水體生態(tài)風(fēng)險評估[J].應(yīng)用生態(tài)學(xué)報, 2014,25(4):1171-1180

Chen J,Liu Y M,Zhang J Y.Aquatic ecological risk assessment of microcystins and nitrogen pollution based on species sensitivity distribution[J].Chinese Journal of Applied Ecology,2014,25(4):1171-1180(in Chinese)

[19] Vanstraalen N M,Denneman C A J.Ecotoxicological evaluation of soil quality criteria[J].Ecotoxicology and Environmental Safety,1989,18(3):241-251

[20] Australian and New Zealand Environment and Conservation Council and Agriculture and Resource Management Council of Australia and New Zealand.Australian and New Zealand Guidelines for Fresh and Marine Water Quality[R].Canberra:Environment and Conservation Council,Agriculture and Resource Management Council, 2000

[21] Raimondo S,Jackson C R,Barron M G.Web-Based Interspecies Correlation Estimation(Web-ICE)for Acute Toxicity:User Manual Version 3.2[DB].Washington DC: US Environmental Protection Agency,Office of Research and Development,National Health and Environmental Effects Research Laboratory,2013

[22] Raimondo S,Jackson C R,Barron M G.Influence of taxonomic relatedness and chemical mode of action in acute interspecies estimation models for aquatic species[J].Environmental Science&Technology,2010,44(19):7711-7716

[23] 王斌,余剛,黃俊,等.應(yīng)用ICE和PCA方法評價硝基芳烴的綜合毒性[J].環(huán)境科學(xué),2007,28(8):1774-1778 Wang B,Yu G,Huang J,et al.Evaluation of integrated toxicity of nitroaromatic compounds by the combination of ICE and PCA[J].Environmental Science,2007,28(8): 1774-1778(in Chinese)

[24] Feng C L,Wu F C,Dyer S D,et al.Derivation of freshwater quality criteria for zinc using interspecies correlation estimation models to protect aquatic life in China[J]. Chemosphere,2013,90(3):1177-1183

[25] European Chemicals Bureau.Technical Guidance for Deriving Environmental Quality Standards.Common Implementation Strategy for the Water Framework Directive (2000/60/EC)[R].Luxembourg:European Chemicals Bureau,European Commission Joint Research Center,European Communities,2011

[26] Hose G C,Van den Brink P J.Confirming the speciessensitivity distribution concept for endosulfan using laboratory,mesocosm,and field data[J].Archives of Environ-mental Contamination and Toxicology,2004,47(4):511-520

[27] Li Z L,Yu J W,Yang M,et al.Cyanobacterial population and harmful metabolites dynamics during a bloom in-Yanghe Reservoir,North China[J].Harmful Algae,2010, 9(5):481-488

[28] 江天久,陳菊芳,鄒迎麟,等.中國東海和南海有害赤潮高發(fā)區(qū)麻痹性貝毒素研究[J].應(yīng)用生態(tài)學(xué)報,2003, 14(7):1156-1160

Jiang T J,Chen J F,Zou Y L,et al.Paralytic shellfish toxins in shellfish from the coast of high frequent harmful algae blooms occurrence areas in East China Sea and South China Sea[J].Chinese Journal of Applied Ecology, 2003,14(7):1156-1160(in Chinese)

[29] Osswald J,Carvalho A P,Guimaraes L,et al.Toxic effects of pure anatoxin-a on biomarkers of rainbow trout, Oncorhynchus mykiss[J].Toxicon,2013,70:162-169

[30] Zagatto P A,Buratini S V,Aragao M A,et al.Neurotoxicity of twoCylindrospermopsis raciborskii(cyanobacteria) strains to mice,daphnia,and fish[J].Environmental Toxicology and Chemistry,2012,31(4):857-862

[31] Guzman-Guillen R,Prieto A I,Vasconcelos V M,et al. Cyanobacterium producing cylindrospermopsin cause oxidative stress at environmentally relevant concentrations in sub-chronically exposed tilapia(Oreochromis niloticus) [J].Chemosphere,2013,90(3):1184-1194

[32] Nogueira I C G,Saker M L,Pflugmacher S,et al.Toxicity of the cyanobacteriumCylindrospermopsis radborskii toDaphnia magna[J].Environmental Toxicology,2004, 19(5):453-459

[33] Aldenberg T,Jaworska J S.Uncertainty of the hazardous concentration and fraction affected for normal species sensitivity distributions[J].Ecotoxicology and Environmental Safety,2000,46(1):1-18◆

Assessing Ecological Risk of Cyanotoxins Based on Interspecies Correlation Estimation and Species Sensitivity Distributions

Zhu Xiaoyi1,2,Chen Jin1,2,Zhang Jianying1,2,*

1.Environmental Science Institute,Zhejiang University,Hangzhou 310058,China
2.Zhejiang Provincial Key Laboratory of Organic Pollution Process and Control,Hangzhou 310058,China

15 July 2015 accepted 17 November 2015

The risk of cyanotoxins to human health has been widely concerned.However,their effect on aquatic species sensitivity distribution was not well understood.This study aimed to assess the ecological risks of two typical cyanobacterial hepatotoxins(nodularins,cylindrospermopsin)and two typical cyanobacterial neurotoxins(anatoxin-a,saxitoxins)by employing the interspecies correlation estimation(ICE)and species sensitivity distributions (SSDs)models.SSDs were constructed based on 71 available sets of acute toxicity data covering 64 aquatic organisms exposure to cyanotoxins.The potentially affected fractions(PAF)under various concentrations and the hazardous concentrations for 5%of the species(HC5)were calculated and compared among different aquatic species.Re-sults show that the estimated HC5values for nodularins,cylindrospermopsin,anatoxin-a and saxitoxins were 74.96 μg·L-1,205.39 μg·L-1,194.39 μg·L-1and 0.3 μg·L-1,respectively,suggesting that the hazardous risk of species is assessed at higher and lower level for saxitoxins and cylindrospermopsin.The invertebrates exposure to cylindrospermopsin and anatoxin-a presented lower HC5than that of the all aquatic organisms with 122.93 μg·L-1and 95.19 μg ·L-1,respectively,which indicate that cylindrospermopsin and anatoxin-a cause more damage to invertebrates than other species.The PAF of 10-1000 μg·L-1cylindrospermopsin and anatoxin-a to invertebrates were higher than that of vertebrates,suggesting that invertebrates are more sensitive to cylindrospermopsin and anatoxin-a than vertebrates.

cyanotoxins;aquatic organisms;species sensitivity distributions;ecological threshold

2015-07-15 錄用日期：2015-11-17

1673-5897（2016）3-131-09

X171.5

A

10.7524/AJE.1673-5897.20150715005

簡介：張建英(1967-)，女，環(huán)境工程博士，副教授，主要研究方向為水污染生態(tài)化學(xué)、環(huán)境風(fēng)險控制技術(shù)，發(fā)表學(xué)術(shù)論文70余篇。

國家自然科學(xué)基金項目(21177106，21477103)

朱小奕(1990-)，女，碩士研究生，研究方向為水生態(tài)毒理學(xué)，E-mail:21414020@zju.edu.cn

*通訊作者（Corresponding author），E-mail:zjy@zju.edu.cn