• <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)方能清洗校園“毒素”
    国产精品九九99| 成年人午夜在线观看视频| 欧美另类一区| 一级片'在线观看视频| 国产成人精品久久二区二区91| 老司机在亚洲福利影院| 国产黄频视频在线观看| 中文字幕人妻丝袜一区二区| 少妇粗大呻吟视频| 大片免费播放器 马上看| 别揉我奶头~嗯~啊~动态视频 | 亚洲精品久久成人aⅴ小说| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美国产精品va在线观看不卡| 国产成人精品久久二区二区免费| 国产成人啪精品午夜网站| 99久久国产精品久久久| 亚洲欧美日韩另类电影网站| 久久久久久免费高清国产稀缺| 免费黄频网站在线观看国产| 亚洲成av片中文字幕在线观看| 国产成人a∨麻豆精品| 大片电影免费在线观看免费| 精品一区二区三区四区五区乱码| 18在线观看网站| 国产精品自产拍在线观看55亚洲 | 亚洲,欧美精品.| 久久国产精品男人的天堂亚洲| 另类精品久久| 91九色精品人成在线观看| e午夜精品久久久久久久| 国产精品一区二区在线观看99| 91麻豆精品激情在线观看国产 | 亚洲伊人久久精品综合| 亚洲国产精品一区二区三区在线| 性高湖久久久久久久久免费观看| 99国产精品免费福利视频| 国产免费av片在线观看野外av| 久久精品亚洲熟妇少妇任你| 操出白浆在线播放| 久久久久精品人妻al黑| 天天躁狠狠躁夜夜躁狠狠躁| 午夜成年电影在线免费观看| 桃红色精品国产亚洲av| 久久久欧美国产精品| 国产精品久久久久久精品古装| 在线观看免费高清a一片| 成人国语在线视频| 丝袜美腿诱惑在线| av欧美777| 午夜两性在线视频| 人人妻人人添人人爽欧美一区卜| 夫妻午夜视频| 叶爱在线成人免费视频播放| 国产主播在线观看一区二区| 日韩 欧美 亚洲 中文字幕| 精品亚洲乱码少妇综合久久| 亚洲午夜精品一区,二区,三区| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美一级毛片孕妇| 两人在一起打扑克的视频| 国产国语露脸激情在线看| 最新的欧美精品一区二区| 久久天躁狠狠躁夜夜2o2o| 亚洲国产精品一区二区三区在线| 男女之事视频高清在线观看| 久久免费观看电影| 这个男人来自地球电影免费观看| av欧美777| 日韩 欧美 亚洲 中文字幕| 在线亚洲精品国产二区图片欧美| 精品亚洲乱码少妇综合久久| 美女午夜性视频免费| 国产亚洲av片在线观看秒播厂| 91九色精品人成在线观看| svipshipincom国产片| 丝袜美腿诱惑在线| 国产av一区二区精品久久| 在线观看www视频免费| 亚洲成人手机| 国产激情久久老熟女| 亚洲精品自拍成人| 久久国产精品大桥未久av| 亚洲精品久久成人aⅴ小说| 老熟妇仑乱视频hdxx| 日日夜夜操网爽| 91麻豆av在线| 中文字幕人妻熟女乱码| 国产成人精品无人区| 岛国毛片在线播放| 一个人免费看片子| 99久久人妻综合| 在线观看舔阴道视频| 欧美日韩国产mv在线观看视频| 亚洲色图综合在线观看| 免费在线观看日本一区| 亚洲欧洲日产国产| 成人国产av品久久久| 久久精品熟女亚洲av麻豆精品| 老司机靠b影院| 啦啦啦啦在线视频资源| 女人被躁到高潮嗷嗷叫费观| 天天躁夜夜躁狠狠躁躁| 国产深夜福利视频在线观看| 十八禁高潮呻吟视频| 国产在线一区二区三区精| 丁香六月天网| 精品亚洲成国产av| 欧美97在线视频| 午夜福利乱码中文字幕| 欧美午夜高清在线| 亚洲男人天堂网一区| 日本黄色日本黄色录像| 啦啦啦 在线观看视频| 久久久国产精品麻豆| 成在线人永久免费视频| 国产精品一区二区在线观看99| 日韩,欧美,国产一区二区三区| 国产欧美日韩精品亚洲av| 岛国毛片在线播放| 天天添夜夜摸| 男女之事视频高清在线观看| 女警被强在线播放| 超色免费av| av网站免费在线观看视频| 一级,二级,三级黄色视频| bbb黄色大片| 欧美另类亚洲清纯唯美| 亚洲精品国产色婷婷电影| 免费在线观看日本一区| 亚洲av片天天在线观看| 热99国产精品久久久久久7| 丰满少妇做爰视频| av电影中文网址| 黑人操中国人逼视频| 久久国产亚洲av麻豆专区| 人人妻人人添人人爽欧美一区卜| 亚洲av电影在线观看一区二区三区| 欧美少妇被猛烈插入视频| 美女扒开内裤让男人捅视频| 欧美日韩亚洲国产一区二区在线观看 | 成年美女黄网站色视频大全免费| 国产福利在线免费观看视频| 亚洲激情五月婷婷啪啪| 永久免费av网站大全| 久久久久久亚洲精品国产蜜桃av| 久久99热这里只频精品6学生| 不卡av一区二区三区| 国产在线一区二区三区精| 亚洲中文av在线| 亚洲熟女精品中文字幕| 久久国产精品人妻蜜桃| 麻豆乱淫一区二区| 亚洲国产精品成人久久小说| 一边摸一边做爽爽视频免费| 午夜影院在线不卡| www.av在线官网国产| 亚洲国产精品一区二区三区在线| 日韩人妻精品一区2区三区| 国产区一区二久久| 在线观看舔阴道视频| 国产在线一区二区三区精| 亚洲av日韩精品久久久久久密| 中文字幕人妻熟女乱码| 国产99久久九九免费精品| 在线 av 中文字幕| 久久狼人影院| 男女下面插进去视频免费观看| 国产成人精品久久二区二区91| 亚洲人成电影观看| 日韩中文字幕视频在线看片| 我要看黄色一级片免费的| 欧美 日韩 精品 国产| 日韩,欧美,国产一区二区三区| 制服诱惑二区| 久久精品人人爽人人爽视色| 亚洲天堂av无毛| 久久久久久久久免费视频了| 麻豆乱淫一区二区| 波多野结衣av一区二区av| 少妇裸体淫交视频免费看高清 | 日本欧美视频一区| 亚洲精品中文字幕一二三四区 | 777米奇影视久久| 午夜福利视频精品| 啦啦啦啦在线视频资源| 国产欧美日韩一区二区三 | 首页视频小说图片口味搜索| 啦啦啦中文免费视频观看日本| 国产精品国产三级国产专区5o| 999精品在线视频| 91精品三级在线观看| 18在线观看网站| 亚洲综合色网址| 啦啦啦中文免费视频观看日本| 麻豆乱淫一区二区| 久久天堂一区二区三区四区| 最近最新免费中文字幕在线| 在线av久久热| 丰满少妇做爰视频| 97在线人人人人妻| 人人妻人人澡人人爽人人夜夜| 久久这里只有精品19| 亚洲欧美日韩高清在线视频 | 性色av一级| 亚洲av美国av| 国产一区二区三区av在线| 欧美另类一区| 丰满饥渴人妻一区二区三| 人人妻人人爽人人添夜夜欢视频| 亚洲激情五月婷婷啪啪| 老司机福利观看| 亚洲精品一卡2卡三卡4卡5卡 | 51午夜福利影视在线观看| 日韩制服骚丝袜av| 国产一区二区 视频在线| 五月开心婷婷网| 国产伦人伦偷精品视频| 在线观看免费日韩欧美大片| 免费不卡黄色视频| 男女高潮啪啪啪动态图| 丝袜美足系列| 狂野欧美激情性bbbbbb| 欧美+亚洲+日韩+国产| 考比视频在线观看| 国产男女超爽视频在线观看| 久久免费观看电影| 免费不卡黄色视频| 在线观看免费高清a一片| 亚洲精品国产区一区二| 精品人妻熟女毛片av久久网站| 超色免费av| 国产精品一区二区在线不卡| 又黄又粗又硬又大视频| 91字幕亚洲| 人妻一区二区av| 永久免费av网站大全| 制服诱惑二区| 国产精品 欧美亚洲| 99九九在线精品视频| 激情视频va一区二区三区| 免费在线观看影片大全网站| 亚洲av电影在线进入| 午夜福利在线免费观看网站| 欧美日韩亚洲综合一区二区三区_| 亚洲av日韩在线播放| 日韩有码中文字幕| 国产成人精品久久二区二区免费| av天堂在线播放| 老司机影院毛片| 热99re8久久精品国产| 午夜日韩欧美国产| 精品亚洲乱码少妇综合久久| 大香蕉久久网| 香蕉国产在线看| 久久国产精品人妻蜜桃| 激情视频va一区二区三区| 国产精品二区激情视频| 精品熟女少妇八av免费久了| 国产成人a∨麻豆精品| 免费黄频网站在线观看国产| 亚洲一码二码三码区别大吗| 法律面前人人平等表现在哪些方面 | 又紧又爽又黄一区二区| 欧美激情 高清一区二区三区| 亚洲色图 男人天堂 中文字幕| 男男h啪啪无遮挡| 丝袜喷水一区| 亚洲精华国产精华精| 50天的宝宝边吃奶边哭怎么回事| 欧美 日韩 精品 国产| 男女免费视频国产| 91大片在线观看| 少妇粗大呻吟视频| 十八禁人妻一区二区| 日韩免费高清中文字幕av| 久久久久精品人妻al黑| 青青草视频在线视频观看| 一区福利在线观看| 桃花免费在线播放| 999久久久国产精品视频| 免费高清在线观看视频在线观看| 一级黄色大片毛片| 精品福利永久在线观看| 亚洲一区中文字幕在线| 国产免费福利视频在线观看| 免费在线观看视频国产中文字幕亚洲 | 老司机午夜福利在线观看视频 | 在线观看一区二区三区激情| 女人精品久久久久毛片| 18禁观看日本| 美女主播在线视频| 999精品在线视频| 国产精品 国内视频| 亚洲精品国产av成人精品| 美女高潮到喷水免费观看| 99久久精品国产亚洲精品| 欧美日韩精品网址| 免费女性裸体啪啪无遮挡网站| 美女扒开内裤让男人捅视频| 伊人久久大香线蕉亚洲五| 国产一区二区三区av在线| 欧美人与性动交α欧美软件| 久久久欧美国产精品| 亚洲精华国产精华精| 各种免费的搞黄视频| 97人妻天天添夜夜摸| 午夜福利一区二区在线看| 三级毛片av免费| 午夜激情av网站| 亚洲成av片中文字幕在线观看| 性少妇av在线| 动漫黄色视频在线观看| 99国产综合亚洲精品| 亚洲精品国产色婷婷电影| 日本五十路高清| 视频在线观看一区二区三区| 黄色a级毛片大全视频| 蜜桃国产av成人99| 90打野战视频偷拍视频| 日韩人妻精品一区2区三区| 欧美精品av麻豆av| 精品国产超薄肉色丝袜足j| 女性生殖器流出的白浆| 日韩中文字幕视频在线看片| 一级a爱视频在线免费观看| 久久亚洲国产成人精品v| 欧美人与性动交α欧美精品济南到| 日日夜夜操网爽| 男人添女人高潮全过程视频| 麻豆国产av国片精品| 一本久久精品| 日本精品一区二区三区蜜桃| 亚洲成国产人片在线观看| 国产亚洲一区二区精品| 国产av精品麻豆| 一级片'在线观看视频| 亚洲欧美激情在线| 女人久久www免费人成看片| 一区在线观看完整版| 亚洲欧美一区二区三区久久| 久久精品熟女亚洲av麻豆精品| 亚洲天堂av无毛| 亚洲欧美清纯卡通| 亚洲精品一二三| 亚洲一卡2卡3卡4卡5卡精品中文| 免费一级毛片在线播放高清视频 | 老司机亚洲免费影院| 日本精品一区二区三区蜜桃| 爱豆传媒免费全集在线观看| 蜜桃国产av成人99| 久久精品国产亚洲av高清一级| 国产在视频线精品| 黄频高清免费视频| a级毛片黄视频| 男女下面插进去视频免费观看| 热99国产精品久久久久久7| 黄频高清免费视频| 老熟妇乱子伦视频在线观看 | 777米奇影视久久| 飞空精品影院首页| 五月开心婷婷网| 高清av免费在线| 亚洲精品中文字幕在线视频| 日韩欧美免费精品| 黄色怎么调成土黄色| 日韩精品免费视频一区二区三区| bbb黄色大片| 少妇 在线观看| 国产亚洲欧美在线一区二区| 欧美中文综合在线视频| 侵犯人妻中文字幕一二三四区| 免费在线观看完整版高清| 在线观看免费午夜福利视频| 亚洲va日本ⅴa欧美va伊人久久 | 日本av免费视频播放| 午夜福利影视在线免费观看| 在线永久观看黄色视频| 五月开心婷婷网| 欧美另类亚洲清纯唯美| 免费在线观看影片大全网站| 亚洲七黄色美女视频| 精品少妇久久久久久888优播| 国产亚洲av片在线观看秒播厂| 黄片大片在线免费观看| 搡老熟女国产l中国老女人| av视频免费观看在线观看| 亚洲人成电影免费在线| 老司机影院成人| 最新在线观看一区二区三区| 动漫黄色视频在线观看| 狠狠狠狠99中文字幕| 精品欧美一区二区三区在线| 欧美 亚洲 国产 日韩一| 精品人妻一区二区三区麻豆| 性高湖久久久久久久久免费观看| 亚洲欧美成人综合另类久久久| 中文字幕高清在线视频| 亚洲av成人不卡在线观看播放网 | 法律面前人人平等表现在哪些方面 | 免费黄频网站在线观看国产| 91国产中文字幕| 亚洲国产精品成人久久小说| 精品国产一区二区三区久久久樱花| 中文字幕最新亚洲高清| 精品人妻1区二区| 国产成人影院久久av| 热99国产精品久久久久久7| 日韩欧美一区视频在线观看| 丝袜人妻中文字幕| 一本大道久久a久久精品| 亚洲,欧美精品.| 自拍欧美九色日韩亚洲蝌蚪91| 自线自在国产av| 香蕉国产在线看| 日韩 亚洲 欧美在线| 久久中文字幕一级| 免费在线观看日本一区| 久久天堂一区二区三区四区| 亚洲第一欧美日韩一区二区三区 | 亚洲午夜精品一区,二区,三区| 亚洲人成77777在线视频| 50天的宝宝边吃奶边哭怎么回事| 久热这里只有精品99| 爱豆传媒免费全集在线观看| 成年av动漫网址| 欧美日韩视频精品一区| 侵犯人妻中文字幕一二三四区| 男女之事视频高清在线观看| 99久久国产精品久久久| 99久久精品国产亚洲精品| 欧美日韩精品网址| 啦啦啦 在线观看视频| av又黄又爽大尺度在线免费看| 天堂俺去俺来也www色官网| 99精国产麻豆久久婷婷| 高清在线国产一区| av线在线观看网站| 1024视频免费在线观看| 国产av国产精品国产| 国产又爽黄色视频| 亚洲熟女精品中文字幕| 亚洲专区字幕在线| 精品人妻熟女毛片av久久网站| 可以免费在线观看a视频的电影网站| 亚洲欧美日韩另类电影网站| 午夜福利免费观看在线| 国产亚洲午夜精品一区二区久久| 国产一级毛片在线| 高清黄色对白视频在线免费看| 亚洲 欧美一区二区三区| 免费黄频网站在线观看国产| 中文欧美无线码| 亚洲va日本ⅴa欧美va伊人久久 | 色精品久久人妻99蜜桃| 人人妻人人添人人爽欧美一区卜| av线在线观看网站| 亚洲激情五月婷婷啪啪| 午夜91福利影院| 亚洲国产精品999| 欧美黄色淫秽网站| 动漫黄色视频在线观看| 国产av精品麻豆| 黑人欧美特级aaaaaa片| 久久久国产欧美日韩av| 午夜老司机福利片| 欧美日韩亚洲高清精品| 国产色视频综合| 一边摸一边做爽爽视频免费| 日韩大码丰满熟妇| 亚洲精品一卡2卡三卡4卡5卡 | 欧美黄色淫秽网站| 国产高清视频在线播放一区 | 亚洲欧美日韩另类电影网站| 香蕉丝袜av| 99九九在线精品视频| 精品人妻熟女毛片av久久网站| 两个人看的免费小视频| 成年动漫av网址| 一本—道久久a久久精品蜜桃钙片| 999久久久国产精品视频| 亚洲五月婷婷丁香| 黄网站色视频无遮挡免费观看| 国产成人精品无人区| 超碰成人久久| 免费观看av网站的网址| 天天躁日日躁夜夜躁夜夜| 免费高清在线观看日韩| 黑人操中国人逼视频| 国产av精品麻豆| 女人爽到高潮嗷嗷叫在线视频| 黑丝袜美女国产一区| av一本久久久久| 欧美乱码精品一区二区三区| 午夜免费鲁丝| 亚洲视频免费观看视频| 国产亚洲午夜精品一区二区久久| 国产91精品成人一区二区三区 | 国产成人免费观看mmmm| 波多野结衣av一区二区av| 欧美黄色淫秽网站| 国产成人系列免费观看| 午夜老司机福利片| 秋霞在线观看毛片| 国产精品自产拍在线观看55亚洲 | www.999成人在线观看| 日日摸夜夜添夜夜添小说| 亚洲专区国产一区二区| 欧美精品一区二区大全| 99热全是精品| 69精品国产乱码久久久| 三上悠亚av全集在线观看| 一区二区三区精品91| 国产亚洲av高清不卡| 最近最新免费中文字幕在线| 啦啦啦啦在线视频资源| 制服诱惑二区| 伊人久久大香线蕉亚洲五| 国产极品粉嫩免费观看在线| 999精品在线视频| 狠狠婷婷综合久久久久久88av| 国产成+人综合+亚洲专区| 99国产综合亚洲精品| 国产亚洲精品一区二区www | 亚洲精品第二区| 丝袜人妻中文字幕| 性色av一级| 亚洲天堂av无毛| 亚洲专区中文字幕在线| 热99re8久久精品国产| 乱人伦中国视频| 一本—道久久a久久精品蜜桃钙片| 国产一区有黄有色的免费视频| 首页视频小说图片口味搜索| 亚洲第一av免费看| 国产激情久久老熟女| 国产不卡av网站在线观看| 久久影院123| 高清视频免费观看一区二区| 国产99久久九九免费精品| 亚洲av日韩精品久久久久久密| netflix在线观看网站| 免费观看av网站的网址| 老司机深夜福利视频在线观看 | 午夜激情av网站| 电影成人av| 动漫黄色视频在线观看| 爱豆传媒免费全集在线观看| 国产欧美日韩精品亚洲av| 五月天丁香电影| 亚洲人成电影免费在线| 久久精品国产亚洲av香蕉五月 | 搡老乐熟女国产| 俄罗斯特黄特色一大片| 亚洲精品在线美女| 中文精品一卡2卡3卡4更新| 成年人免费黄色播放视频| 免费日韩欧美在线观看| 久久久国产成人免费| av欧美777| 丝袜美足系列| 一区二区三区精品91| 亚洲国产毛片av蜜桃av| 黄色a级毛片大全视频| 国产成人精品久久二区二区91| 18禁观看日本| 欧美精品av麻豆av| 国产亚洲精品久久久久5区| 美女视频免费永久观看网站| 免费久久久久久久精品成人欧美视频| 女性生殖器流出的白浆| 麻豆乱淫一区二区| 可以免费在线观看a视频的电影网站| 一二三四社区在线视频社区8| 美女国产高潮福利片在线看| 国产成人av教育| 热99久久久久精品小说推荐| 女人爽到高潮嗷嗷叫在线视频| 国产一区二区三区av在线| 青春草视频在线免费观看| 亚洲成人手机| 日韩大码丰满熟妇| 亚洲欧美色中文字幕在线| av在线播放精品| 亚洲精品自拍成人| av在线老鸭窝| 免费一级毛片在线播放高清视频 | 9色porny在线观看| 久久99一区二区三区| 自线自在国产av| 久久国产精品人妻蜜桃| 国产福利在线免费观看视频| 成年av动漫网址| 精品第一国产精品| 男女下面插进去视频免费观看| 日韩欧美国产一区二区入口| 国产精品一区二区免费欧美 | 999精品在线视频| √禁漫天堂资源中文www| 亚洲九九香蕉| 亚洲国产成人一精品久久久| 正在播放国产对白刺激| 亚洲精品中文字幕一二三四区 | 99国产综合亚洲精品| 国产一级毛片在线| 欧美少妇被猛烈插入视频| 国内毛片毛片毛片毛片毛片| 999精品在线视频| 丰满人妻熟妇乱又伦精品不卡| 国产亚洲精品久久久久5区| 亚洲专区中文字幕在线| 久久热在线av| 伊人亚洲综合成人网| 极品少妇高潮喷水抽搐| 亚洲 欧美一区二区三区|