• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2016-10-27 06:16:07朱小奕陳瑾張建英
    生態(tài)毒理學(xué)報(bào) 2016年3期
    關(guān)鍵詞:類毒素球藻毒素

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

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

    2.浙江省有機(jī)污染過(guò)程與控制重點(diǎn)實(shí)驗(yàn)室,杭州310058

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

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

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

    2.浙江省有機(jī)污染過(guò)程與控制重點(diǎn)實(shí)驗(yàn)室,杭州310058

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

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

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

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

    目前,對(duì)微囊藻毒素的生態(tài)風(fēng)險(xiǎn)研究已從生物個(gè)體水平發(fā)展到到生態(tài)系統(tǒng)水平,但除此之外,水環(huán)境中節(jié)球藻毒素、柱孢藻毒素、類毒素和貝毒素等的研究?jī)H限于生物個(gè)體毒性,尚未報(bào)道其在系統(tǒng)尺度上的生態(tài)風(fēng)險(xiǎn)。本研究通過(guò)挖掘文獻(xiàn)中4種藻毒素的急性毒性試驗(yàn)數(shù)據(jù),構(gòu)建ICE-SSD復(fù)合模型,評(píng)價(jià)不同濃度水平藻毒素污染的水生態(tài)物種敏感性風(fēng)險(xiǎn),計(jì)算藻毒素對(duì)水生生物的風(fēng)險(xiǎn)閾值。擬提供在系統(tǒng)尺度上評(píng)價(jià)藻毒素污染風(fēng)險(xiǎn)的一個(gè)保證率為95%的可靠評(píng)價(jià)模型,從而為藻毒素污染水體的生態(tài)風(fēng)險(xiǎn)管理提供方法和依據(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曲線,用以評(píng)價(jià)其水生態(tài)風(fēng)險(xiǎn),數(shù)據(jù)來(lái)源于Web of Science(http:// www.webofknowledge.com/)收錄的相關(guān)文獻(xiàn)資料,其篩選原則參照歐盟生態(tài)風(fēng)險(xiǎn)評(píng)估標(biāo)準(zhǔn)“風(fēng)險(xiǎn)評(píng)估技術(shù)綱領(lǐng)”等[25]:(1)標(biāo)準(zhǔn)毒性試驗(yàn)方法,如經(jīng)濟(jì)合作與發(fā)展組織(Organization for Economic Cooperation and Development,OECD)或美國(guó)材料與試驗(yàn)協(xié)會(huì)(A-merican Society for Testing and Materials,ASMT)等發(fā)布方法;(2)實(shí)驗(yàn)條件的設(shè)定,如暴露周期、介質(zhì)等,具體篩選條件見(jiàn)表1[9],本研究中節(jié)球藻毒素、柱孢藻毒素和類毒素以淡水分布為主,選取其對(duì)淡水生物毒性試驗(yàn)結(jié)果,而貝毒素主要來(lái)源于海洋藻釋放,評(píng)價(jià)其對(duì)海洋生物的急性毒性;(3)急性毒性終點(diǎn)選取,以半數(shù)效應(yīng)濃度EC50表征,其效應(yīng)包括運(yùn)動(dòng)抑制、生長(zhǎng)抑制、胚胎發(fā)育抑制和致死率等中毒癥狀;(4)評(píng)價(jià)物種數(shù)不少于5個(gè),且盡量分布于3個(gè)以上營(yíng)養(yǎng)級(jí)。為彌補(bǔ)柱孢藻毒素、類毒性、貝毒素適當(dāng)毒性數(shù)據(jù)的缺失,本文進(jìn)一步采用EPA的ICE模型(http:// www.epa.gov/ceampubl/fchain/webice/),通過(guò)選取2個(gè)具有代表性的物種(大型蚤Daphnia magna,Daphnia pulex)為替代物種外推同類生物急性毒性數(shù)據(jù)。

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

    1.2 SSD擬合

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

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

    當(dāng)k趨于無(wú)窮大時(shí),BurrIII分布變化為Re-Weibull分布:

    當(dāng)c趨于無(wú)窮大時(shí),變化為RePareto分布:

    2 結(jié)果(Results)

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

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

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

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

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

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

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

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

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

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

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

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

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

    圖1 全部物種對(duì)5種藻毒素的物種敏感性分布(SSD)曲線圖

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

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

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

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

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

    3 討論(Discussion)

    圖2 脊柱動(dòng)物(I)、無(wú)脊椎動(dòng)物(II)對(duì)微囊藻毒素、柱孢藻毒素和類毒素的SSD曲線比較

    圖3 脊椎動(dòng)物、無(wú)脊椎動(dòng)物對(duì)柱孢藻毒素(I)、類毒素(II)的SSD曲線比較

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

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

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

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

    (References):

    [1] 謝平.微囊藻毒素對(duì)人類健康影響相關(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)用物種敏感性分布評(píng)估DDT和林丹對(duì)淡水生物的生態(tài)風(fēng)險(xiǎn)[J].環(huán)境科學(xué)學(xué)報(bào),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)用物種敏感性分布評(píng)估有機(jī)磷農(nóng)藥對(duì)淡水生物的急性生態(tài)風(fēng)險(xiǎn)[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)用物種敏感性分布評(píng)估重金屬對(duì)淡水生物的急性生態(tài)風(fēng)險(xiǎn)[J].中國(guó)環(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] 蔣丹烈,胡霞林,尹大強(qiáng).應(yīng)用物種敏感性分布法對(duì)太湖沉積物中多環(huán)芳烴的生態(tài)風(fēng)險(xiǎn)分析[J].生態(tài)毒理學(xué)報(bào),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)險(xiǎn)評(píng)估[J].應(yīng)用生態(tài)學(xué)報(bào), 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方法評(píng)價(jià)硝基芳烴的綜合毒性[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] 江天久,陳菊芳,鄒迎麟,等.中國(guó)東海和南海有害赤潮高發(fā)區(qū)麻痹性貝毒素研究[J].應(yīng)用生態(tài)學(xué)報(bào),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

    簡(jiǎn)介:張建英(1967-),女,環(huán)境工程博士,副教授,主要研究方向?yàn)樗廴旧鷳B(tài)化學(xué)、環(huán)境風(fēng)險(xiǎn)控制技術(shù),發(fā)表學(xué)術(shù)論文70余篇。

    國(guó)家自然科學(xué)基金項(xiàng)目(21177106,21477103)

    朱小奕(1990-),女,碩士研究生,研究方向?yàn)樗鷳B(tài)毒理學(xué),E-mail:21414020@zju.edu.cn

    *通訊作者(Corresponding author),E-mail:zjy@zju.edu.cn

    猜你喜歡
    類毒素球藻毒素
    真菌Simplicillium lanosoniveum DT06 對(duì)雨生紅球藻生長(zhǎng)與脂類合成的影響
    可再生能源(2022年8期)2022-08-17 06:37:52
    What Makes You Tired
    一種艱難梭菌類毒素A的制備方法
    QuEChERS-液相色譜-四極桿飛行時(shí)間質(zhì)譜法測(cè)定毒蘑菇中5種鵝膏肽類毒素
    一類具有毒素的非均勻chemostat模型正解的存在性和唯一性
    球藻沉浮的秘密植物
    大自然探索(2019年2期)2019-03-01 02:23:30
    毒蘑菇中毒素的研究進(jìn)展
    聚球藻硅質(zhì)化作用初探
    不同血液凈化模式對(duì)維持性血液透析患者蛋白結(jié)合類毒素的清除效果
    嚴(yán)苛標(biāo)準(zhǔn)方能清洗校園“毒素”
    美女大奶头黄色视频| cao死你这个sao货| 亚洲国产成人一精品久久久| 国产精品一区二区在线不卡| 国产精品偷伦视频观看了| 精品少妇黑人巨大在线播放| 日韩中文字幕欧美一区二区| 国产人伦9x9x在线观看| 久久久久国内视频| 精品一区二区三卡| 国产亚洲av高清不卡| 久久午夜综合久久蜜桃| 亚洲人成电影免费在线| 啦啦啦 在线观看视频| 午夜免费成人在线视频| 男男h啪啪无遮挡| 久久精品久久久久久噜噜老黄| 亚洲人成77777在线视频| 亚洲成人免费av在线播放| 老司机深夜福利视频在线观看 | 男女无遮挡免费网站观看| 黄色a级毛片大全视频| 黄片大片在线免费观看| 亚洲精品久久成人aⅴ小说| 欧美激情久久久久久爽电影 | 国产精品av久久久久免费| 人人妻,人人澡人人爽秒播| 国产片内射在线| 日韩大码丰满熟妇| 99re6热这里在线精品视频| 纯流量卡能插随身wifi吗| 少妇猛男粗大的猛烈进出视频| 精品熟女少妇八av免费久了| 捣出白浆h1v1| 天天影视国产精品| 亚洲专区中文字幕在线| 免费在线观看日本一区| 后天国语完整版免费观看| 欧美日本中文国产一区发布| 好男人电影高清在线观看| 在线观看免费高清a一片| 两人在一起打扑克的视频| 免费女性裸体啪啪无遮挡网站| 丁香六月天网| 国产精品久久久久成人av| 91麻豆av在线| 亚洲激情五月婷婷啪啪| 99re6热这里在线精品视频| 亚洲国产欧美日韩在线播放| 咕卡用的链子| 国产国语露脸激情在线看| 90打野战视频偷拍视频| 国产亚洲av片在线观看秒播厂| 曰老女人黄片| 国产99久久九九免费精品| 可以免费在线观看a视频的电影网站| 亚洲一区二区三区欧美精品| 午夜福利在线免费观看网站| 他把我摸到了高潮在线观看 | 超碰97精品在线观看| 欧美人与性动交α欧美软件| 高潮久久久久久久久久久不卡| 成人国语在线视频| 中国美女看黄片| 亚洲av日韩精品久久久久久密| av视频免费观看在线观看| 国产深夜福利视频在线观看| 国产精品免费视频内射| 大片免费播放器 马上看| 国产人伦9x9x在线观看| 亚洲免费av在线视频| 成人手机av| 精品国产乱码久久久久久男人| 午夜福利,免费看| 国产亚洲av片在线观看秒播厂| 热re99久久精品国产66热6| 欧美人与性动交α欧美软件| av国产精品久久久久影院| 成人三级做爰电影| 一区在线观看完整版| 永久免费av网站大全| 国产欧美亚洲国产| 欧美在线黄色| 国产一级毛片在线| av不卡在线播放| 黄色怎么调成土黄色| 国产男人的电影天堂91| 91国产中文字幕| 在线看a的网站| 亚洲精品国产精品久久久不卡| 日韩精品免费视频一区二区三区| 久久久久久人人人人人| 欧美精品亚洲一区二区| 咕卡用的链子| 黄色怎么调成土黄色| √禁漫天堂资源中文www| bbb黄色大片| 亚洲精品自拍成人| 亚洲精品中文字幕一二三四区 | 亚洲男人天堂网一区| av欧美777| 天天躁狠狠躁夜夜躁狠狠躁| 国产高清国产精品国产三级| 在线观看免费午夜福利视频| 大香蕉久久成人网| 亚洲精品日韩在线中文字幕| 性色av乱码一区二区三区2| 国产精品一区二区在线观看99| www.999成人在线观看| 亚洲成人国产一区在线观看| 久热这里只有精品99| 欧美精品高潮呻吟av久久| 亚洲精品一区蜜桃| 18禁观看日本| 欧美亚洲日本最大视频资源| 日韩一区二区三区影片| 午夜激情av网站| 亚洲午夜精品一区,二区,三区| 欧美人与性动交α欧美软件| 午夜91福利影院| 美女福利国产在线| 久久99热这里只频精品6学生| 亚洲 欧美一区二区三区| 淫妇啪啪啪对白视频 | 看免费av毛片| 国产成人一区二区三区免费视频网站| 欧美激情极品国产一区二区三区| 黑人巨大精品欧美一区二区mp4| 久久性视频一级片| 岛国毛片在线播放| 国产精品久久久av美女十八| 老司机在亚洲福利影院| 国产一区二区 视频在线| 一进一出抽搐动态| 天堂中文最新版在线下载| 考比视频在线观看| 18禁观看日本| 欧美性长视频在线观看| 国产精品一区二区在线不卡| 操出白浆在线播放| 欧美在线一区亚洲| 老司机午夜十八禁免费视频| 啪啪无遮挡十八禁网站| 在线天堂中文资源库| 国产日韩欧美亚洲二区| 一边摸一边抽搐一进一出视频| 免费av中文字幕在线| 国产欧美日韩综合在线一区二区| 欧美激情高清一区二区三区| 日韩精品免费视频一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 欧美性长视频在线观看| 少妇的丰满在线观看| 男女下面插进去视频免费观看| 日韩 欧美 亚洲 中文字幕| 无限看片的www在线观看| 多毛熟女@视频| 国产一区二区三区av在线| 美女高潮喷水抽搐中文字幕| av在线app专区| 亚洲精品一卡2卡三卡4卡5卡 | 久久午夜综合久久蜜桃| 亚洲av国产av综合av卡| av有码第一页| 一区福利在线观看| 国产一区有黄有色的免费视频| 天堂俺去俺来也www色官网| 真人做人爱边吃奶动态| 狂野欧美激情性xxxx| 又大又爽又粗| 女人久久www免费人成看片| 国产亚洲av片在线观看秒播厂| 天堂俺去俺来也www色官网| 热re99久久国产66热| 18禁裸乳无遮挡动漫免费视频| 久久久精品免费免费高清| 蜜桃国产av成人99| 欧美黑人精品巨大| 国产一级毛片在线| www日本在线高清视频| 中文字幕最新亚洲高清| 少妇裸体淫交视频免费看高清 | 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品国产av在线观看| 国产无遮挡羞羞视频在线观看| 9191精品国产免费久久| 99香蕉大伊视频| 美女高潮喷水抽搐中文字幕| www.自偷自拍.com| 人人妻,人人澡人人爽秒播| 性色av一级| 91精品国产国语对白视频| 色综合欧美亚洲国产小说| 久久九九热精品免费| 国产伦人伦偷精品视频| 热99久久久久精品小说推荐| 老汉色∧v一级毛片| 在线观看免费高清a一片| 成年人黄色毛片网站| 精品免费久久久久久久清纯 | 丝袜脚勾引网站| 国产1区2区3区精品| 国产精品 欧美亚洲| 日韩制服骚丝袜av| 美女大奶头黄色视频| 在线精品无人区一区二区三| 天天影视国产精品| 久久久久久人人人人人| 又紧又爽又黄一区二区| 黄片播放在线免费| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品亚洲av国产电影网| 永久免费av网站大全| 两个人免费观看高清视频| 国产视频一区二区在线看| 欧美性长视频在线观看| 亚洲精品久久午夜乱码| 精品久久蜜臀av无| 亚洲人成电影观看| 俄罗斯特黄特色一大片| 他把我摸到了高潮在线观看 | 亚洲中文av在线| 日韩欧美一区二区三区在线观看 | 免费看十八禁软件| 久久精品人人爽人人爽视色| 亚洲第一青青草原| 国产视频一区二区在线看| 丰满少妇做爰视频| 看免费av毛片| 女人精品久久久久毛片| 亚洲人成电影观看| 天堂中文最新版在线下载| 精品国产一区二区三区四区第35| 悠悠久久av| 欧美精品亚洲一区二区| 欧美乱码精品一区二区三区| 久久综合国产亚洲精品| 精品久久蜜臀av无| 久久久欧美国产精品| 黑人巨大精品欧美一区二区mp4| 18禁裸乳无遮挡动漫免费视频| 国产精品偷伦视频观看了| cao死你这个sao货| 久久精品国产亚洲av香蕉五月 | 男女午夜视频在线观看| 亚洲精品成人av观看孕妇| 女性被躁到高潮视频| 99久久综合免费| xxxhd国产人妻xxx| 欧美大码av| 一本色道久久久久久精品综合| 国产老妇伦熟女老妇高清| 免费av中文字幕在线| 欧美黄色淫秽网站| 日本撒尿小便嘘嘘汇集6| 欧美日韩成人在线一区二区| 欧美日韩亚洲综合一区二区三区_| 啦啦啦 在线观看视频| 亚洲男人天堂网一区| 天天躁夜夜躁狠狠躁躁| 久久久久久久国产电影| 国产无遮挡羞羞视频在线观看| netflix在线观看网站| 9色porny在线观看| 可以免费在线观看a视频的电影网站| 黄色片一级片一级黄色片| 久久狼人影院| 十八禁高潮呻吟视频| 少妇粗大呻吟视频| 国产成+人综合+亚洲专区| 国产成人av教育| 国产不卡av网站在线观看| 亚洲精品av麻豆狂野| 日本wwww免费看| 精品人妻熟女毛片av久久网站| 久久久久久免费高清国产稀缺| 久久影院123| 男女边摸边吃奶| 久久中文看片网| 妹子高潮喷水视频| 久久午夜综合久久蜜桃| 母亲3免费完整高清在线观看| 日本欧美视频一区| netflix在线观看网站| 婷婷丁香在线五月| av网站免费在线观看视频| 18禁国产床啪视频网站| 一区福利在线观看| 亚洲视频免费观看视频| h视频一区二区三区| 在线观看人妻少妇| 91成人精品电影| 亚洲全国av大片| 国产精品影院久久| 99久久人妻综合| 欧美精品啪啪一区二区三区 | 国产成+人综合+亚洲专区| 又黄又粗又硬又大视频| 热re99久久国产66热| 欧美日韩国产mv在线观看视频| 精品久久久精品久久久| 热re99久久国产66热| 国产1区2区3区精品| 色播在线永久视频| www.精华液| 男女高潮啪啪啪动态图| 黄片播放在线免费| 成人国语在线视频| 国产极品粉嫩免费观看在线| 国产有黄有色有爽视频| 国产福利在线免费观看视频| 高清黄色对白视频在线免费看| 欧美精品一区二区大全| 俄罗斯特黄特色一大片| 精品免费久久久久久久清纯 | 亚洲专区中文字幕在线| 日韩欧美国产一区二区入口| 一边摸一边做爽爽视频免费| 一二三四在线观看免费中文在| 午夜两性在线视频| 久久国产精品影院| 亚洲三区欧美一区| 水蜜桃什么品种好| 人人澡人人妻人| 电影成人av| 国产亚洲欧美在线一区二区| 午夜影院在线不卡| 丝袜喷水一区| 三上悠亚av全集在线观看| 777久久人妻少妇嫩草av网站| 一本大道久久a久久精品| 欧美少妇被猛烈插入视频| 精品国产一区二区三区久久久樱花| av在线老鸭窝| a 毛片基地| 啦啦啦视频在线资源免费观看| 国产麻豆69| 岛国毛片在线播放| av欧美777| 99国产极品粉嫩在线观看| 国产亚洲精品第一综合不卡| 黄色 视频免费看| 丁香六月天网| 欧美黑人精品巨大| 久久中文看片网| 另类亚洲欧美激情| 亚洲,欧美精品.| 黄片小视频在线播放| 亚洲中文av在线| 久久中文字幕一级| 欧美精品一区二区大全| 纵有疾风起免费观看全集完整版| 亚洲人成电影观看| 精品熟女少妇八av免费久了| 国产成人免费观看mmmm| 夜夜夜夜夜久久久久| 一本—道久久a久久精品蜜桃钙片| 777米奇影视久久| 18禁黄网站禁片午夜丰满| 日本一区二区免费在线视频| 99热网站在线观看| 亚洲第一欧美日韩一区二区三区 | 欧美黑人精品巨大| 国产亚洲欧美在线一区二区| 性色av乱码一区二区三区2| 91av网站免费观看| 亚洲av国产av综合av卡| 国产福利在线免费观看视频| 丝袜人妻中文字幕| 免费在线观看黄色视频的| 18在线观看网站| 欧美日韩福利视频一区二区| 免费在线观看黄色视频的| 老司机午夜福利在线观看视频 | 欧美人与性动交α欧美精品济南到| 色婷婷av一区二区三区视频| 夫妻午夜视频| 一级毛片女人18水好多| 韩国高清视频一区二区三区| 涩涩av久久男人的天堂| 80岁老熟妇乱子伦牲交| 久久久国产精品麻豆| 69av精品久久久久久 | 国产亚洲午夜精品一区二区久久| 午夜久久久在线观看| 精品第一国产精品| 国产日韩一区二区三区精品不卡| 色播在线永久视频| 久久久久国产精品人妻一区二区| 亚洲成人免费电影在线观看| 日本91视频免费播放| 久久中文看片网| 亚洲精品美女久久av网站| 人人妻人人爽人人添夜夜欢视频| 中文字幕高清在线视频| 亚洲成国产人片在线观看| 一进一出抽搐动态| 99香蕉大伊视频| 如日韩欧美国产精品一区二区三区| 高清黄色对白视频在线免费看| 久久人妻福利社区极品人妻图片| 91精品三级在线观看| 国产不卡av网站在线观看| 亚洲伊人色综图| www.自偷自拍.com| 91字幕亚洲| 国产高清视频在线播放一区 | 亚洲午夜精品一区,二区,三区| 在线观看免费高清a一片| 久久久久网色| 国产精品影院久久| 国产欧美日韩一区二区三 | 人人妻,人人澡人人爽秒播| 一个人免费看片子| 99精品欧美一区二区三区四区| 中文字幕另类日韩欧美亚洲嫩草| av福利片在线| 精品第一国产精品| a级毛片黄视频| 国产亚洲精品第一综合不卡| 免费在线观看完整版高清| 天天操日日干夜夜撸| 日韩制服骚丝袜av| 啦啦啦在线免费观看视频4| 韩国精品一区二区三区| 久久av网站| 狠狠狠狠99中文字幕| 99热网站在线观看| 免费一级毛片在线播放高清视频 | 日韩中文字幕欧美一区二区| 久久av网站| 一级,二级,三级黄色视频| 婷婷丁香在线五月| 欧美乱码精品一区二区三区| 黄色 视频免费看| 国产又爽黄色视频| 国产精品av久久久久免费| 国产一区二区在线观看av| 午夜免费鲁丝| 亚洲中文字幕日韩| 精品视频人人做人人爽| 黑丝袜美女国产一区| 久9热在线精品视频| 正在播放国产对白刺激| 欧美精品啪啪一区二区三区 | 9191精品国产免费久久| 美女视频免费永久观看网站| 亚洲精品久久成人aⅴ小说| 亚洲欧美精品自产自拍| 韩国高清视频一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 伊人久久大香线蕉亚洲五| 亚洲精品一二三| 欧美xxⅹ黑人| 午夜久久久在线观看| 黄色a级毛片大全视频| a 毛片基地| 亚洲国产看品久久| 国产在线一区二区三区精| 韩国高清视频一区二区三区| 国产在线视频一区二区| 国产伦人伦偷精品视频| 国产人伦9x9x在线观看| 国产一区二区 视频在线| 悠悠久久av| 黑丝袜美女国产一区| 亚洲一码二码三码区别大吗| 99热网站在线观看| 亚洲专区中文字幕在线| 少妇裸体淫交视频免费看高清 | 亚洲精品第二区| 精品国产一区二区三区久久久樱花| 成人影院久久| 日韩电影二区| 爱豆传媒免费全集在线观看| 精品亚洲乱码少妇综合久久| 国产一区二区激情短视频 | 日韩,欧美,国产一区二区三区| 窝窝影院91人妻| 十八禁网站网址无遮挡| 精品熟女少妇八av免费久了| 国产高清视频在线播放一区 | 老熟妇仑乱视频hdxx| 天天操日日干夜夜撸| 多毛熟女@视频| 久久国产精品男人的天堂亚洲| 男女之事视频高清在线观看| 成人三级做爰电影| 国产亚洲av高清不卡| 国产精品国产av在线观看| 国产伦人伦偷精品视频| 亚洲美女黄色视频免费看| 丰满少妇做爰视频| 另类亚洲欧美激情| 桃红色精品国产亚洲av| av天堂久久9| 天天躁夜夜躁狠狠躁躁| 五月开心婷婷网| 美女高潮到喷水免费观看| 精品久久蜜臀av无| 每晚都被弄得嗷嗷叫到高潮| 精品免费久久久久久久清纯 | 欧美少妇被猛烈插入视频| 考比视频在线观看| 精品一区在线观看国产| 久久天堂一区二区三区四区| 多毛熟女@视频| 国产免费现黄频在线看| 男人操女人黄网站| 国产不卡av网站在线观看| 国产精品99久久99久久久不卡| 纯流量卡能插随身wifi吗| 黄色a级毛片大全视频| 在线观看一区二区三区激情| 国产精品99久久99久久久不卡| 免费在线观看视频国产中文字幕亚洲 | 69av精品久久久久久 | 9191精品国产免费久久| 桃花免费在线播放| 国产熟女午夜一区二区三区| 日本av免费视频播放| 又大又爽又粗| 满18在线观看网站| 国产精品久久久久久人妻精品电影 | 欧美久久黑人一区二区| 两性夫妻黄色片| 交换朋友夫妻互换小说| 天天操日日干夜夜撸| 视频在线观看一区二区三区| 人成视频在线观看免费观看| 久久天躁狠狠躁夜夜2o2o| 国产伦人伦偷精品视频| 亚洲精品久久成人aⅴ小说| 亚洲欧美激情在线| 久久人人爽人人片av| 天堂俺去俺来也www色官网| 三上悠亚av全集在线观看| 80岁老熟妇乱子伦牲交| 一级黄色大片毛片| 在线观看www视频免费| 精品一区二区三卡| 天天操日日干夜夜撸| 亚洲精品国产一区二区精华液| 亚洲五月色婷婷综合| 91精品三级在线观看| 国产一区二区在线观看av| 亚洲av成人一区二区三| 亚洲 国产 在线| 色94色欧美一区二区| 久久ye,这里只有精品| 中文字幕制服av| 亚洲人成77777在线视频| 欧美另类亚洲清纯唯美| 国产精品1区2区在线观看. | 黄色视频在线播放观看不卡| 亚洲男人天堂网一区| 桃花免费在线播放| 久久人人爽人人片av| av有码第一页| 国产成人免费观看mmmm| 97精品久久久久久久久久精品| 他把我摸到了高潮在线观看 | 日韩精品免费视频一区二区三区| 久久天躁狠狠躁夜夜2o2o| 久久99热这里只频精品6学生| 天天躁夜夜躁狠狠躁躁| 亚洲av欧美aⅴ国产| 亚洲av片天天在线观看| 另类亚洲欧美激情| 一区二区av电影网| 精品福利观看| 午夜福利在线免费观看网站| 国产精品 欧美亚洲| 91字幕亚洲| 精品高清国产在线一区| 少妇精品久久久久久久| 欧美亚洲 丝袜 人妻 在线| 午夜91福利影院| h视频一区二区三区| 日本wwww免费看| 少妇 在线观看| 日韩视频一区二区在线观看| 国产精品一区二区在线不卡| 国产深夜福利视频在线观看| 色精品久久人妻99蜜桃| 午夜91福利影院| kizo精华| 亚洲伊人久久精品综合| 亚洲专区国产一区二区| 久久精品人人爽人人爽视色| 亚洲国产精品一区三区| 丝袜喷水一区| 欧美人与性动交α欧美精品济南到| a 毛片基地| 成年美女黄网站色视频大全免费| 欧美久久黑人一区二区| 搡老岳熟女国产| 国产精品秋霞免费鲁丝片| 国产老妇伦熟女老妇高清| 岛国毛片在线播放| 黑人猛操日本美女一级片| 十八禁高潮呻吟视频| a在线观看视频网站| 精品国产国语对白av| 欧美黑人精品巨大| 两人在一起打扑克的视频| 国产亚洲一区二区精品| 香蕉丝袜av| 国产精品熟女久久久久浪| 大陆偷拍与自拍| 在线观看舔阴道视频| 黄片小视频在线播放| 老司机午夜福利在线观看视频 | 香蕉丝袜av| 午夜福利视频精品|