• <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)方能清洗校園“毒素”
    欧美zozozo另类| 久久精品久久久久久久性| av天堂中文字幕网| 纯流量卡能插随身wifi吗| 在线观看国产h片| 久久99热这里只有精品18| 国产一区二区三区av在线| 中文在线观看免费www的网站| 一本色道久久久久久精品综合| 天堂中文最新版在线下载| 欧美一级a爱片免费观看看| 亚洲av免费高清在线观看| 一级片'在线观看视频| 午夜日本视频在线| 久久国产精品男人的天堂亚洲 | 欧美xxⅹ黑人| 亚洲av在线观看美女高潮| 国产综合精华液| 成人毛片a级毛片在线播放| 80岁老熟妇乱子伦牲交| 欧美xxxx黑人xx丫x性爽| 中文欧美无线码| 国产午夜精品一二区理论片| 国产亚洲精品久久久com| 少妇人妻久久综合中文| 久久久久久久久久成人| 少妇熟女欧美另类| 一二三四中文在线观看免费高清| 国产精品av视频在线免费观看| 国产精品麻豆人妻色哟哟久久| 久久 成人 亚洲| 黄色日韩在线| 久热久热在线精品观看| 中国三级夫妇交换| 熟女电影av网| 中文字幕免费在线视频6| 国语对白做爰xxxⅹ性视频网站| 国产精品人妻久久久影院| 女人久久www免费人成看片| 国产久久久一区二区三区| 大片免费播放器 马上看| 欧美精品一区二区免费开放| 小蜜桃在线观看免费完整版高清| 成年免费大片在线观看| 国产综合精华液| 亚洲av成人精品一二三区| 精品酒店卫生间| 一个人免费看片子| 能在线免费看毛片的网站| 亚洲欧洲国产日韩| 伊人久久精品亚洲午夜| 亚洲精品一区蜜桃| 男人爽女人下面视频在线观看| 国产精品三级大全| 久久热精品热| 男女无遮挡免费网站观看| 欧美日韩国产mv在线观看视频 | 亚洲国产最新在线播放| 亚洲三级黄色毛片| av线在线观看网站| 国产亚洲一区二区精品| 日本欧美国产在线视频| 久久韩国三级中文字幕| 啦啦啦在线观看免费高清www| 国产精品一及| 搡女人真爽免费视频火全软件| 一区在线观看完整版| av免费观看日本| 日产精品乱码卡一卡2卡三| 久久精品人妻少妇| 亚洲av欧美aⅴ国产| 中文字幕av成人在线电影| 男人舔奶头视频| 3wmmmm亚洲av在线观看| 午夜福利视频精品| 欧美高清性xxxxhd video| 成人综合一区亚洲| 成人毛片60女人毛片免费| 五月伊人婷婷丁香| 校园人妻丝袜中文字幕| 美女高潮的动态| 久久久精品免费免费高清| 亚洲电影在线观看av| 久久精品久久久久久噜噜老黄| 成人18禁高潮啪啪吃奶动态图 | 在线观看三级黄色| 国产片特级美女逼逼视频| 乱码一卡2卡4卡精品| 国产精品.久久久| 免费观看a级毛片全部| 有码 亚洲区| 久久午夜福利片| 免费在线观看成人毛片| 日日摸夜夜添夜夜添av毛片| 午夜福利网站1000一区二区三区| 日韩av免费高清视频| 一区二区三区免费毛片| 亚洲精品一区蜜桃| 久久人妻熟女aⅴ| 亚洲欧美精品专区久久| 亚洲精品国产成人久久av| 最近手机中文字幕大全| 最近的中文字幕免费完整| 精品久久久久久久末码| 国产淫片久久久久久久久| 韩国av在线不卡| 国产 一区 欧美 日韩| 在线观看一区二区三区激情| 国产在线一区二区三区精| 直男gayav资源| 国产成人精品一,二区| 久久亚洲国产成人精品v| 大香蕉久久网| 国产 一区精品| 一区二区av电影网| 多毛熟女@视频| 大香蕉97超碰在线| 一级av片app| 久久精品国产a三级三级三级| 国产深夜福利视频在线观看| 插逼视频在线观看| 日韩伦理黄色片| 国产成人91sexporn| 日本黄大片高清| 日韩制服骚丝袜av| 黑人猛操日本美女一级片| 免费观看性生交大片5| 少妇裸体淫交视频免费看高清| av国产久精品久网站免费入址| 夜夜骑夜夜射夜夜干| 建设人人有责人人尽责人人享有的 | av国产久精品久网站免费入址| 水蜜桃什么品种好| 久久精品国产a三级三级三级| 成人一区二区视频在线观看| 三级经典国产精品| 欧美bdsm另类| 精品午夜福利在线看| 黄色怎么调成土黄色| 国产精品熟女久久久久浪| 亚洲国产毛片av蜜桃av| 99精国产麻豆久久婷婷| 亚洲欧美精品专区久久| 噜噜噜噜噜久久久久久91| 女性被躁到高潮视频| 欧美3d第一页| 日韩欧美一区视频在线观看 | 少妇高潮的动态图| 午夜福利网站1000一区二区三区| 亚洲美女视频黄频| 在线观看免费高清a一片| 国产精品熟女久久久久浪| 国产精品久久久久久精品古装| 亚洲av在线观看美女高潮| 三级经典国产精品| 日本午夜av视频| 日韩大片免费观看网站| 人妻少妇偷人精品九色| 乱码一卡2卡4卡精品| 97超碰精品成人国产| 人妻一区二区av| 黄色日韩在线| 日本av免费视频播放| 99久久精品热视频| 舔av片在线| 精品人妻熟女av久视频| 国产成人精品一,二区| 日韩国内少妇激情av| 久久精品久久久久久久性| 欧美人与善性xxx| 三级国产精品欧美在线观看| 国产精品福利在线免费观看| 日日啪夜夜爽| 直男gayav资源| 黄色视频在线播放观看不卡| 天天躁日日操中文字幕| 九色成人免费人妻av| 日日撸夜夜添| 亚洲欧美日韩无卡精品| 日本黄大片高清| 久久婷婷青草| 精品一区二区三卡| 精品亚洲成国产av| 国产无遮挡羞羞视频在线观看| 久久久久国产精品人妻一区二区| 日本av手机在线免费观看| 黑人猛操日本美女一级片| 美女xxoo啪啪120秒动态图| 18禁在线播放成人免费| 男女国产视频网站| 夜夜骑夜夜射夜夜干| 男人狂女人下面高潮的视频| 欧美一级a爱片免费观看看| 成人毛片a级毛片在线播放| 少妇被粗大猛烈的视频| 久久久精品94久久精品| 天天躁日日操中文字幕| 午夜免费观看性视频| 成人亚洲欧美一区二区av| 精品熟女少妇av免费看| 亚洲av成人精品一二三区| 国产综合精华液| 亚洲国产色片| 日韩欧美一区视频在线观看 | 国产成人午夜福利电影在线观看| 夜夜看夜夜爽夜夜摸| 中文欧美无线码| 美女内射精品一级片tv| 高清在线视频一区二区三区| 亚洲国产日韩一区二区| 一区二区av电影网| 美女xxoo啪啪120秒动态图| 国产精品久久久久成人av| 久久这里有精品视频免费| 日本-黄色视频高清免费观看| 黄片wwwwww| 国模一区二区三区四区视频| 久久精品国产鲁丝片午夜精品| 内射极品少妇av片p| 日本wwww免费看| 成人18禁高潮啪啪吃奶动态图 | 亚洲国产成人一精品久久久| 久久国产精品男人的天堂亚洲 | 国产一区二区三区av在线| 一本色道久久久久久精品综合| 日韩亚洲欧美综合| 国产精品国产av在线观看| 十八禁网站网址无遮挡 | 少妇人妻精品综合一区二区| 亚洲欧美清纯卡通| av视频免费观看在线观看| 九草在线视频观看| 你懂的网址亚洲精品在线观看| 我的老师免费观看完整版| 精品久久久久久电影网| 国产精品秋霞免费鲁丝片| 婷婷色综合www| 超碰97精品在线观看| 在线观看免费视频网站a站| 亚洲成色77777| 又粗又硬又长又爽又黄的视频| 国产亚洲精品久久久com| 国产亚洲av片在线观看秒播厂| av网站免费在线观看视频| 亚州av有码| 春色校园在线视频观看| 九九爱精品视频在线观看| 久久人妻熟女aⅴ| 免费不卡的大黄色大毛片视频在线观看| 国产综合精华液| 国产人妻一区二区三区在| av线在线观看网站| 看免费成人av毛片| 久久久久久九九精品二区国产| 哪个播放器可以免费观看大片| 亚洲精品国产色婷婷电影| 免费观看无遮挡的男女| 99热这里只有精品一区| 国产色婷婷99| av网站免费在线观看视频| 一本—道久久a久久精品蜜桃钙片| 久久精品久久精品一区二区三区| 免费黄网站久久成人精品| 亚洲久久久国产精品| 国产爱豆传媒在线观看| 精品国产三级普通话版| 嫩草影院入口| 精品亚洲成国产av| 日本免费在线观看一区| av免费观看日本| 亚洲成人手机| 亚洲图色成人| 精品国产一区二区三区久久久樱花 | 亚洲欧美日韩东京热| 免费播放大片免费观看视频在线观看| 老熟女久久久| 男人和女人高潮做爰伦理| 久久婷婷青草| 久久精品国产亚洲av涩爱| 99视频精品全部免费 在线| 日韩欧美 国产精品| 中国三级夫妇交换| 女人十人毛片免费观看3o分钟| 99热这里只有是精品在线观看| 一级毛片我不卡| 国产亚洲一区二区精品| 小蜜桃在线观看免费完整版高清| 我的老师免费观看完整版| 又粗又硬又长又爽又黄的视频| 免费高清在线观看视频在线观看| 精品久久久久久久久av| 国产精品一区www在线观看| 久久久久久久久久久丰满| 一个人看视频在线观看www免费| 国内揄拍国产精品人妻在线| 欧美一区二区亚洲| 国产一级毛片在线| 18禁动态无遮挡网站| 夫妻性生交免费视频一级片| 观看av在线不卡| 日韩av免费高清视频| 久久久色成人| 国产一区亚洲一区在线观看| 香蕉精品网在线| freevideosex欧美| 99久久精品一区二区三区| 久久99热这里只频精品6学生| 青青草视频在线视频观看| 久久久久网色| 嘟嘟电影网在线观看| av专区在线播放| 女性生殖器流出的白浆| 亚洲欧美精品专区久久| 高清欧美精品videossex| 国产视频首页在线观看| 少妇人妻久久综合中文| 亚洲美女视频黄频| 爱豆传媒免费全集在线观看| 2018国产大陆天天弄谢| 国产精品一区二区在线不卡| 天堂8中文在线网| 久久 成人 亚洲| 日韩欧美精品免费久久| 日日摸夜夜添夜夜爱| 久久 成人 亚洲| 久久影院123| 成人影院久久| 国内揄拍国产精品人妻在线| 午夜日本视频在线| 高清视频免费观看一区二区| 国产大屁股一区二区在线视频| 国产成人精品福利久久| 欧美精品一区二区免费开放| 一二三四中文在线观看免费高清| 直男gayav资源| 亚洲电影在线观看av| 一本色道久久久久久精品综合| 最近中文字幕高清免费大全6| 美女高潮的动态| 99国产精品免费福利视频| 亚洲天堂av无毛| 亚洲国产色片| 日韩视频在线欧美| 欧美精品人与动牲交sv欧美| 亚洲欧洲日产国产| 久久久久久人妻| 欧美区成人在线视频| 亚洲一级一片aⅴ在线观看| 最近2019中文字幕mv第一页| 91精品一卡2卡3卡4卡| 国产 一区精品| 色视频在线一区二区三区| 在线观看一区二区三区激情| 亚洲激情五月婷婷啪啪| 日本午夜av视频| 婷婷色麻豆天堂久久| 久久久久久九九精品二区国产| 自拍偷自拍亚洲精品老妇| 国产亚洲5aaaaa淫片| 亚洲激情五月婷婷啪啪| 自拍欧美九色日韩亚洲蝌蚪91 | 国产v大片淫在线免费观看| 97在线人人人人妻| 色吧在线观看| 人妻制服诱惑在线中文字幕| 我的老师免费观看完整版| 国产极品天堂在线| 亚洲精品色激情综合| 中文精品一卡2卡3卡4更新| 亚州av有码| 青春草亚洲视频在线观看| 最近手机中文字幕大全| 免费黄网站久久成人精品| 蜜桃在线观看..| 99视频精品全部免费 在线| av.在线天堂| 亚洲美女黄色视频免费看| 精品午夜福利在线看| 在线观看美女被高潮喷水网站| 久久午夜福利片| 国产91av在线免费观看| 国产日韩欧美在线精品| 欧美激情极品国产一区二区三区 | 国产成人免费无遮挡视频| 偷拍熟女少妇极品色| 日韩av免费高清视频| 又大又黄又爽视频免费| 国产真实伦视频高清在线观看| 亚洲美女视频黄频| 久久久国产一区二区| 99九九线精品视频在线观看视频| 精品99又大又爽又粗少妇毛片| 欧美日韩一区二区视频在线观看视频在线| 精品久久久久久电影网| 三级国产精品片| 性高湖久久久久久久久免费观看| 夫妻午夜视频| av在线老鸭窝| 亚洲精品自拍成人| 欧美zozozo另类| 国产久久久一区二区三区| 婷婷色综合大香蕉| 在线 av 中文字幕| 人妻制服诱惑在线中文字幕| 青青草视频在线视频观看| 高清在线视频一区二区三区| 婷婷色麻豆天堂久久| 美女内射精品一级片tv| 国产伦精品一区二区三区四那| 亚洲天堂av无毛| 最后的刺客免费高清国语| 国产色婷婷99| 蜜臀久久99精品久久宅男| 国产av国产精品国产| 国产日韩欧美在线精品| 国产黄频视频在线观看| 亚洲婷婷狠狠爱综合网| 久久毛片免费看一区二区三区| 视频区图区小说| 老司机影院成人| 尾随美女入室| 日韩不卡一区二区三区视频在线| 国内精品宾馆在线| 亚洲国产精品999| 成人国产麻豆网| 久久久久性生活片| 欧美区成人在线视频| 舔av片在线| 精品人妻一区二区三区麻豆| 男女边摸边吃奶| 一区二区三区免费毛片| 中文资源天堂在线| 久久久成人免费电影| 99热这里只有是精品在线观看| 亚洲成色77777| 国产精品免费大片| 国产爽快片一区二区三区| 国产免费福利视频在线观看| 啦啦啦中文免费视频观看日本| 国产大屁股一区二区在线视频| av播播在线观看一区| 夜夜骑夜夜射夜夜干| 亚洲av中文字字幕乱码综合| 日本色播在线视频| av不卡在线播放| 91精品一卡2卡3卡4卡| 亚洲av男天堂| 国产视频内射| 身体一侧抽搐| 久久ye,这里只有精品| 国产大屁股一区二区在线视频| 亚洲精品日韩在线中文字幕| 18禁在线无遮挡免费观看视频| 国产午夜精品一二区理论片| 国产精品国产三级国产av玫瑰| a级毛色黄片| 99视频精品全部免费 在线| 观看美女的网站| 国产日韩欧美在线精品| 成人国产麻豆网| 久久人人爽人人片av| 国产av精品麻豆| av一本久久久久| 新久久久久国产一级毛片| 中文欧美无线码| 日韩免费高清中文字幕av| 91精品一卡2卡3卡4卡| 国产精品人妻久久久影院| av黄色大香蕉| 熟妇人妻不卡中文字幕| 国产有黄有色有爽视频| 欧美日韩视频高清一区二区三区二| 免费少妇av软件| 久久精品国产自在天天线| 精品人妻熟女av久视频| 精品酒店卫生间| 激情五月婷婷亚洲| 亚洲精华国产精华液的使用体验| 久久国产精品大桥未久av | 欧美3d第一页| 日本欧美视频一区| 一级毛片黄色毛片免费观看视频| 国产成人免费无遮挡视频| 在线播放无遮挡| 亚洲内射少妇av| 又爽又黄a免费视频| 黄片无遮挡物在线观看| 亚洲综合精品二区| av在线蜜桃| 热re99久久精品国产66热6| 国产乱人偷精品视频| 国产中年淑女户外野战色| 精品一区二区三卡| 在线天堂最新版资源| 成人国产av品久久久| 亚洲精品久久久久久婷婷小说| 国产精品精品国产色婷婷| 青春草国产在线视频| 久久久久久久久久成人| 日韩伦理黄色片| 毛片一级片免费看久久久久| 国产又色又爽无遮挡免| 午夜福利高清视频| 一边亲一边摸免费视频| 色婷婷久久久亚洲欧美| 欧美一级a爱片免费观看看| av网站免费在线观看视频| 男女边吃奶边做爰视频| 乱系列少妇在线播放| 男人添女人高潮全过程视频| 22中文网久久字幕| 国产精品人妻久久久影院| 中文字幕精品免费在线观看视频 | 大香蕉97超碰在线| 国产一区二区三区综合在线观看 | 欧美 日韩 精品 国产| 午夜福利高清视频| av国产久精品久网站免费入址| 插逼视频在线观看| 国产高清有码在线观看视频| 国产av一区二区精品久久 | 青春草亚洲视频在线观看| 高清在线视频一区二区三区| 日韩三级伦理在线观看| 纯流量卡能插随身wifi吗| 国产精品一区二区性色av| 亚洲欧美日韩卡通动漫| 99精国产麻豆久久婷婷| 嫩草影院入口| 日日摸夜夜添夜夜添av毛片| 一边亲一边摸免费视频| 亚洲人成网站在线观看播放| 日韩av免费高清视频| 国产色爽女视频免费观看| 22中文网久久字幕| 建设人人有责人人尽责人人享有的 | 天天躁夜夜躁狠狠久久av| 少妇猛男粗大的猛烈进出视频| 蜜桃在线观看..| 爱豆传媒免费全集在线观看| 国产日韩欧美亚洲二区| 亚洲综合色惰| 亚洲国产成人一精品久久久| 国产一区二区三区av在线| a 毛片基地| 欧美精品一区二区免费开放| 在线免费观看不下载黄p国产| 精品国产一区二区三区久久久樱花 | av在线蜜桃| 久久久久性生活片| 国产黄色免费在线视频| 丝袜脚勾引网站| 女人久久www免费人成看片| 国产精品久久久久成人av| 国产一区二区三区综合在线观看 | 少妇 在线观看| 在线精品无人区一区二区三 | 老女人水多毛片| 欧美日韩视频精品一区| 直男gayav资源| 这个男人来自地球电影免费观看 | 观看美女的网站| 国产又色又爽无遮挡免| 成人特级av手机在线观看| 看免费成人av毛片| 啦啦啦中文免费视频观看日本| 小蜜桃在线观看免费完整版高清| 欧美zozozo另类| 不卡视频在线观看欧美| 天美传媒精品一区二区| 国产精品国产av在线观看| 一级a做视频免费观看| 日本欧美视频一区| 国产一区亚洲一区在线观看| 国产成人一区二区在线| 青春草亚洲视频在线观看| 精品视频人人做人人爽| 黑人猛操日本美女一级片| 在线观看免费视频网站a站| 丰满乱子伦码专区| 少妇的逼水好多| 五月玫瑰六月丁香| 在线观看免费视频网站a站| 一边亲一边摸免费视频| 看免费成人av毛片| 久久 成人 亚洲| 伦精品一区二区三区| 高清不卡的av网站| 九草在线视频观看| 国产高清有码在线观看视频| 国产熟女欧美一区二区| 久久久久人妻精品一区果冻| 日韩一区二区视频免费看| 国国产精品蜜臀av免费| 免费人成在线观看视频色| 国产精品国产三级国产av玫瑰| 秋霞伦理黄片| 夫妻午夜视频| 国产一区有黄有色的免费视频| 国产男女内射视频| 国产成人a区在线观看| 中文精品一卡2卡3卡4更新| 久久亚洲国产成人精品v| 青春草亚洲视频在线观看| 高清视频免费观看一区二区| 性色avwww在线观看| 肉色欧美久久久久久久蜜桃| 18禁裸乳无遮挡免费网站照片| 汤姆久久久久久久影院中文字幕| 国产精品一二三区在线看| 日产精品乱码卡一卡2卡三| 精品亚洲成a人片在线观看 | 成人影院久久| 免费大片黄手机在线观看| 亚洲国产精品成人久久小说|