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

    海河天津段微囊藻及其毒素的空間分布及與水環(huán)境因子的關(guān)系

    2018-10-29 05:00:24王素炎陳裴裴李一鳴喬之怡
    中國(guó)環(huán)境科學(xué) 2018年10期
    關(guān)鍵詞:產(chǎn)毒海河微囊

    霍 達(dá),曹 琪,王素炎,陳裴裴,李一鳴,喬之怡*

    ?

    海河天津段微囊藻及其毒素的空間分布及與水環(huán)境因子的關(guān)系

    霍 達(dá)1,2,曹 琪1,王素炎1,陳裴裴1,李一鳴1,喬之怡1,2*

    (1.天津農(nóng)學(xué)院水產(chǎn)學(xué)院,天津 300384;2.天津市水產(chǎn)生態(tài)及養(yǎng)殖重點(diǎn)實(shí)驗(yàn)室,天津 300384)

    采取熒光定量PCR技術(shù)及ELISA酶聯(lián)免疫吸附法,以基因和-基因?yàn)榘谢驅(qū)ο募舅{(lán)藻暴發(fā)期間海河天津市區(qū)段的產(chǎn)毒微囊藻種群豐度和毒素含量進(jìn)行了研究.結(jié)果表明:夏季海河天津市區(qū)段微囊藻種群豐度具有明顯差異性:產(chǎn)毒微囊藻種群豐度為1.16×104~2.48×107copies/mL占總微囊藻種群的4.25%~28.59%.藻毒素含量最高點(diǎn)為8號(hào)采樣點(diǎn)天津站,每升水體中藻細(xì)胞共含毒素195.51μg,除去藻細(xì)胞水中毒素濃度為0.97μg/L.總的來(lái)說(shuō),海河天津市區(qū)段微囊藻總基因拷貝數(shù)較高,產(chǎn)毒微囊藻豐度在不同采樣點(diǎn)間差異較大,影響水體中微囊藻豐度的最主要環(huán)境因素是pH值.海河水體中毒素含量沒有超過(guò)安全閾值,但是單位體積水體中藻細(xì)胞內(nèi)毒素含量很高,具有較高的毒素釋放潛力,可能對(duì)下游河口生態(tài)系統(tǒng)造成潛在威脅.

    微囊藻;微囊藻毒素;熒光定量PCR技術(shù);ELISA酶聯(lián)免疫吸附法;海河

    微囊藻水華給淡水環(huán)境帶來(lái)了極大的影響.一些種類的微囊藻能產(chǎn)生次生代謝物微囊藻毒素(MC),該毒素是富營(yíng)養(yǎng)化水體中最主要的藍(lán)藻毒素類型,可引發(fā)水生生物中毒甚至死亡,直接威脅水中魚類及其它生物的生存.飲用水源中的微囊藻毒素會(huì)對(duì)人類的肝臟造成損傷,影響人類的健康和生存[1],同時(shí),微囊藻毒素是強(qiáng)致癌劑,與人群原發(fā)性肝癌的高發(fā)關(guān)系密切[2].為保護(hù)人類身體健康,世界衛(wèi)生組織規(guī)定的飲用水中微囊藻毒素安全濃度閾值<1.0μg/L[3].研究表明,微囊藻水華種群中同時(shí)含有有毒和無(wú)毒微囊藻基因型,兩者種群豐度及其比例的動(dòng)態(tài)變化直接影響了微囊藻水華毒素濃度,對(duì)微囊藻毒素濃度具有很好的指示作用[4-6].有毒和無(wú)毒微囊藻細(xì)胞形態(tài)一致,用傳統(tǒng)的顯微鏡形態(tài)觀察方法無(wú)法對(duì)兩者進(jìn)行區(qū)分和鑒定.有毒微囊藻細(xì)胞基因組中含有微囊藻毒素合成酶基因家族(mcy)[7-8],對(duì)毒素產(chǎn)生進(jìn)行調(diào)控,而無(wú)毒微囊藻株中則不含有mcy 基因家族.因此可以應(yīng)用PCR擴(kuò)增mcy基因單個(gè)或多個(gè)成員檢測(cè)和鑒別有毒和無(wú)毒藻株[9-11].

    近些年來(lái),隨著分子生物學(xué)技術(shù)的發(fā)展,越來(lái)越多的研究人員對(duì)于微囊藻的產(chǎn)毒相關(guān)基因進(jìn)行了定量的研究.Rinta-Kanto等[12]以基因和微囊藻16S rDNA基因?yàn)榘谢虿捎脽晒舛縋CR的方法分析了Western Lake Erie的產(chǎn)毒微囊藻構(gòu)成及其豐度,發(fā)現(xiàn)了該湖中微囊藻存在產(chǎn)毒和非產(chǎn)毒共存的現(xiàn)象,在國(guó)內(nèi)的研究中,學(xué)者們[13-15]對(duì)太湖及巢湖水華期間微囊藻種群豐度進(jìn)行了研究,發(fā)現(xiàn)了相同現(xiàn)象在中國(guó)南方水體中的廣泛存在.劉洋等[16]對(duì)南太湖入湖口的產(chǎn)毒微囊藻進(jìn)行了研究,證明了該方法可以快速準(zhǔn)確地檢測(cè)水體中微囊藻產(chǎn)毒藻數(shù)量.

    隨著全球氣溫變化的影響,藍(lán)藻暴發(fā)的范圍慢慢由南方水體延伸至北方.北方水體具有冬季結(jié)冰,四季氣溫變化大等特征,與南方水體差異較大.已有研究表明在溫度[17]、緯度[18]等因素影響下,微囊藻毒素的釋放具有不同特征.海河貫穿天津市區(qū),隨著近年來(lái)兩岸生活污水的排放,以及上游水質(zhì)的惡化,海河市區(qū)內(nèi)景觀河道出現(xiàn)了大面積的微囊藻水華,對(duì)生活用水造成了較大的威脅.本研究將天津海河作為模型,從毒素合成調(diào)控基因和毒素兩方面入手,研究北方緩速流動(dòng)型水體的毒素釋放特征,為水環(huán)境治理提供參考資料.

    1 材料與方法

    1.1 采樣位點(diǎn)及化學(xué)指標(biāo)測(cè)定方法

    表1 采樣點(diǎn)信息

    2015年7月,在海河天津市區(qū)段(坐標(biāo))設(shè)置8個(gè)采樣點(diǎn)(表1),分別為天津站、耳閘摩天輪、天津?yàn)彻珗@富民橋、獅子林橋.使用不銹鋼采水器采取1.5m水柱,混合均勻.同時(shí)利用YSI多功能水質(zhì)檢測(cè)儀測(cè)定采樣點(diǎn)水質(zhì)基本參數(shù)如溫度、pH值、氨氮、溶解氧等.TN、TP測(cè)定方法參照文獻(xiàn)[14].

    1.2 基因組DNA提取

    取200mL混勻后的水樣(根據(jù)微囊藻的量調(diào)整抽濾體積)使用Millipore玻璃纖維膜過(guò)濾.含藻濾膜置于-20℃保存.將膜剪碎后采用Rinta-Kanto等[12]的方法提取DNA.

    圖1 采樣點(diǎn)示意

    1.3 構(gòu)建熒光定量PCR 標(biāo)準(zhǔn)曲線

    熒光定量PCR所選用的引物如表2所示,其中,引物188F和254R用于特異性擴(kuò)增微囊藻藻藍(lán)蛋白轉(zhuǎn)錄間隔區(qū)PC-IGS[4],引物F2和R2用于特異性擴(kuò)增有毒微囊藻毒素合成基因[19].標(biāo)準(zhǔn)藻株為有毒微囊藻株FACHB 905,由中國(guó)淡水藻種庫(kù)(FACHB)提供,在實(shí)驗(yàn)室內(nèi)培養(yǎng),培養(yǎng)條件(光照強(qiáng)度為30mmol protons/(m2×s),光照周期12h/12h (L:D),藻細(xì)胞濃度培養(yǎng)到對(duì)數(shù)生長(zhǎng)期進(jìn)行采收).采用全細(xì)胞PCR的方法擴(kuò)增2個(gè)片段,將PCR產(chǎn)物進(jìn)行瓊脂糖電泳后回收.將純化后的產(chǎn)物與pEASY?克隆載體連接,并注入大腸桿菌感受態(tài)細(xì)胞Trans1-T1中,于LB液體培養(yǎng)基中振蕩培養(yǎng);離心收集菌液,混勻,取適量涂布在含有LB(Amp+)的固體培養(yǎng)皿中過(guò)夜培養(yǎng);挑取單個(gè)陽(yáng)性菌落按上述條件進(jìn)一步培養(yǎng),再對(duì)培養(yǎng)后菌落中的質(zhì)粒進(jìn)行測(cè)序.利用質(zhì)粒提取試劑盒提取質(zhì)粒,用超微量核酸蛋白測(cè)定儀(eppendorf)測(cè)定質(zhì)粒標(biāo)準(zhǔn)品的濃度.根據(jù)文獻(xiàn)[20]計(jì)算拷貝數(shù),將得到的質(zhì)粒標(biāo)準(zhǔn)品用1×TE梯度稀釋構(gòu)建標(biāo)準(zhǔn)曲線.

    表2 mcyD基因和PC-IGS基因qPCR反應(yīng)引物

    1.4 定量PCR反應(yīng)體系

    定量PCR反應(yīng)體系為20μL,其中含有10μL SYBR GREEN MIX,2pmoL引物,1μL模板,超純水補(bǔ)足20μL, PCR反應(yīng)在Bio-rad CFX Connet中進(jìn)行.-基因擴(kuò)增程序如下:95℃預(yù)變性2min, 95 ℃變性15s,58℃退火30s,72℃延伸30s,40個(gè)循環(huán).基因定量PCR反應(yīng)退火溫度為56℃,其它條件與PC-IGS相同.

    1.5 水中及細(xì)胞中微囊藻毒素檢測(cè)

    采用0.22μm一次性真空濾膜過(guò)濾5mL水樣,保留含藻濾膜進(jìn)行水中微囊藻毒素含量測(cè)定.細(xì)胞中的毒素使用GF/C膜過(guò)濾藻細(xì)胞并調(diào)節(jié)濃度后測(cè)定.所有毒素檢測(cè)樣本均使用中科院水生生物研究所研制的微囊藻毒素ELISA試劑盒進(jìn)行檢測(cè).

    1.6 數(shù)據(jù)統(tǒng)計(jì)

    數(shù)據(jù)統(tǒng)計(jì)在R語(yǔ)言環(huán)境下進(jìn)行,相關(guān)性計(jì)算采用Pearson相關(guān)性函數(shù)計(jì)算.數(shù)據(jù)展示采用平均值±標(biāo)準(zhǔn)差(mean±SD)形式.

    2 結(jié)果與分析

    2.1 海河水體理化參數(shù)

    測(cè)得采樣期間水體段理化指標(biāo)見表3,采樣期間海河水體的平均溫度為29.2℃,4號(hào)采樣點(diǎn)水溫最高為30.4℃.水體呈弱堿性,5號(hào)采樣點(diǎn)pH值最高.溶解氧變化范圍較廣,為4.77~17.15mg/L.總氮最高值為4號(hào)點(diǎn)3.55mg/L,最低值為1號(hào)點(diǎn)1.04mg/L.總磷最高點(diǎn)為2號(hào)點(diǎn)0.1mg/L,最低值為4號(hào)點(diǎn)0.05mg/L.

    表3 海河水體主要理化參數(shù)

    2.2 定量PCR標(biāo)準(zhǔn)曲線

    利用質(zhì)粒DNA為模板建立基因和基因定量PCR標(biāo)準(zhǔn)曲線.從圖2和圖3可以看出,循環(huán)閾值q與基因組拷貝數(shù)的對(duì)數(shù)呈顯著線性關(guān)系.定量PCR反應(yīng)中,序列定量PCR標(biāo)準(zhǔn)曲線方程為= -3.482+ 38.79 ,(2=0.993),擴(kuò)增效率為93.8%.-序列定量PCR標(biāo)準(zhǔn)曲線方程為= -3.41+ 39.97,(2=0.999),擴(kuò)增效率為96.5%.溶解曲線結(jié)果顯示和-的擴(kuò)增產(chǎn)物解鏈峰值溫度分別為82.49℃(變異系數(shù)CV=0.09%)和84.41℃(變異系數(shù)CV=0.26%)

    圖2 mcyD基因熒光定量PCR標(biāo)準(zhǔn)曲線

    圖3 PC-IGS基因熒光定量PCR標(biāo)準(zhǔn)曲線

    2.3 產(chǎn)毒微囊藻和總微囊藻種群豐度的分布

    定量PCR結(jié)果如圖4所示,海河微囊藻種群豐度具有顯著的空間分布差異.5號(hào)點(diǎn)天津之眼微囊藻種群豐度最高,每毫升水中微囊藻基因型拷貝數(shù)達(dá)到(1.53×108±1.13×107)copies/mL同時(shí)基因型拷貝數(shù)達(dá)到(2.48×107±7.12×106)copies/mL.3號(hào)點(diǎn)北洋園微囊藻種群豐度最低,微囊基因型拷貝數(shù)為(8.50×104±1.62×106) copies/mL,基因型拷貝數(shù)為(1.16×104±1.09×102)copies /mL.產(chǎn)毒微囊藻基因豐度與總微囊藻基因豐度比值為4.25%~28.59%平均比值為15.44%,各點(diǎn)之間差異較大.

    圖4 產(chǎn)毒微囊藻和總微囊藻種群豐度的空間分布

    2.4 微囊藻細(xì)胞和水體中微囊藻毒素含量

    ELISA酶聯(lián)免疫反應(yīng)檢測(cè)結(jié)果如圖5所示,藻細(xì)胞中毒素含量最高的點(diǎn)位為8號(hào)天津站為每升水中藻細(xì)胞含195.51μg毒素,最低3號(hào)為2.53μg.水中微囊藻毒素含量最高的也為8號(hào)天津站0.97μg/L.最低點(diǎn)3號(hào)0.18μg/L.4~8號(hào)點(diǎn)位顯著高于1~3號(hào)點(diǎn).通過(guò)線性回歸得到細(xì)胞產(chǎn)毒和水中微囊藻毒素相關(guān)性2=0.95.

    圖5 每升水體中細(xì)胞內(nèi)和細(xì)胞外的微囊藻細(xì)胞產(chǎn)毒素含量

    3 討論

    3.1 基于熒光定量PCR結(jié)果的微囊藻基因豐度差異

    近年來(lái),夏季中的海河備受藍(lán)藻困擾,藍(lán)藻暴發(fā)對(duì)居民生活和城市景觀帶來(lái)較大影響.當(dāng)前推測(cè)海河藍(lán)藻暴發(fā)的原因有:海河生態(tài)用水不足,水體流動(dòng)性差;海河水體中適合藻類生長(zhǎng)的營(yíng)養(yǎng)物質(zhì)含量較高;天津夏季氣候非常適應(yīng)藍(lán)藻生長(zhǎng)等[21],但均缺乏定量的分析.此外對(duì)水體中微囊藻是否產(chǎn)毒,產(chǎn)毒微囊藻的豐富度如何也不甚清楚,缺乏相應(yīng)的監(jiān)測(cè)和報(bào)道.

    近些年,熒光定量PCR在藍(lán)藻水華的研究中得到廣泛應(yīng)用.這項(xiàng)技術(shù)能夠?qū)Υ笈鷺悠吠瑫r(shí)進(jìn)行處理,對(duì)產(chǎn)毒微囊藻種群豐富度進(jìn)行定量化的研究,快速對(duì)藍(lán)藻水華是否產(chǎn)毒和毒性做出判斷,可以預(yù)測(cè)藍(lán)藻水華帶來(lái)的生態(tài)風(fēng)險(xiǎn)[22].基因qPCR法具有靈敏性高、成本較低、所用時(shí)間短、所需檢測(cè)樣本少、操作簡(jiǎn)單等特點(diǎn)[23].qPCR法對(duì)于微囊藻產(chǎn)毒潛能的預(yù)測(cè)與毒素的直接檢測(cè)結(jié)果具有高度一致性[24].研究中常選擇單個(gè)基因作為目的基因進(jìn)行qPCR擴(kuò)增片段[25-30],也可選取多個(gè)基因作為目標(biāo)基因進(jìn)行擴(kuò)增[31-33],本實(shí)驗(yàn)選取基因作為目的基因?qū)λw中產(chǎn)毒微囊藻進(jìn)行檢測(cè).-基因?yàn)樘禺愋詳U(kuò)增微囊藻的基因片段,可檢測(cè)出水體中總微囊藻基因豐富度.本實(shí)驗(yàn)通過(guò)建立基因和-基因的熒光定量PCR標(biāo)準(zhǔn)曲線,對(duì)海河水體中產(chǎn)毒微囊藻種群基因豐富度和總微囊藻種群基因豐富度進(jìn)行定量分析.通過(guò)熒光定量PCR的結(jié)果表明,海河干流不同水系的微囊藻種群豐富度具有較大差異,北洋園(3號(hào)位點(diǎn))產(chǎn)毒微囊藻豐富度與總微囊藻豐富度均最低;南運(yùn)河(1號(hào)位點(diǎn))和新開河(4號(hào)位點(diǎn))產(chǎn)毒微囊藻豐富度與總微囊藻豐富度均較高,因此,海河干流中的微囊藻可能主要來(lái)源于南運(yùn)河和新開河.

    3.2 水體理化指標(biāo)與微囊藻產(chǎn)毒基因及毒素的空間分布關(guān)系

    通過(guò)對(duì)8月份海河天津市區(qū)段水體理化指標(biāo)分析,可知此時(shí)的水體處于富營(yíng)養(yǎng)化,再根據(jù)qPCR結(jié)果印證了在富營(yíng)養(yǎng)化狀態(tài)下微囊藻的總豐度較高[34].在自然環(huán)境中,許多種環(huán)境因子都可能影響產(chǎn)毒微囊藻占總微囊藻種群的比例,在Yoshida等[35]研究Mikata湖泊中產(chǎn)毒微囊藻和非產(chǎn)毒微囊藻種群動(dòng)態(tài)隨時(shí)間變化規(guī)律時(shí),發(fā)現(xiàn)硝態(tài)氮濃度與產(chǎn)毒微囊藻占總微囊藻種群的比值呈顯著正相關(guān).這與本實(shí)驗(yàn)中南運(yùn)河(愛民橋 1號(hào)位點(diǎn))硝態(tài)氮濃度很高,產(chǎn)毒微囊藻占總微囊藻種群的比值也高的實(shí)驗(yàn)結(jié)論相符合.從研究結(jié)果可以看出,產(chǎn)毒微囊藻占總微囊藻的比值范圍4.25%~28.59%,平均值為15.50%.其中比值最高的點(diǎn)為獅子林橋(6號(hào)位點(diǎn)),比值最低的點(diǎn)為進(jìn)步橋(7號(hào)位點(diǎn)).相關(guān)報(bào)道表明不同湖泊產(chǎn)毒微囊藻占總微囊藻種群的比值的波動(dòng)范圍較大[36-37],獅子林橋與進(jìn)步橋的距離較為相近,但產(chǎn)毒藻株比例卻差異很大,導(dǎo)致這種結(jié)果的原因還需要進(jìn)一步研究.根據(jù)劉永定等關(guān)于中國(guó)藻類學(xué)的研究表明高pH (大于8)值有利于原核藻類生長(zhǎng)[38],海河水體8月pH值均在8以上,呈現(xiàn)弱堿性,非常有利于微囊藻生長(zhǎng),導(dǎo)致藍(lán)藻暴發(fā).6號(hào)位點(diǎn)pH值顯示為9.01,7號(hào)位點(diǎn)為8.01,推測(cè)藻類豐富度可能與pH值有關(guān).5號(hào)位點(diǎn)pH為9.3,是pH值最高的點(diǎn),5號(hào)位點(diǎn)也是基因豐富度和-基因豐富度最高的點(diǎn).產(chǎn)毒微囊藻基因豐富度與總微囊藻基因豐富度的比值最高的點(diǎn)為獅子林橋(6號(hào)位點(diǎn)),比值最低的點(diǎn)為進(jìn)步橋(7號(hào)位點(diǎn)),其中6號(hào)位點(diǎn)的溶解氧9.11mg/L大于7號(hào)位點(diǎn)的溶解氧4.77mg/L.新開河(4號(hào)位點(diǎn))溶解氧為最高17.5mg/L,產(chǎn)毒微囊藻基因豐富度也最高為2.08×106copies/mL,溶解氧最低的點(diǎn)為進(jìn)步橋(7號(hào)位點(diǎn)),產(chǎn)毒微囊藻基因豐富度的最低點(diǎn)為北洋園(3號(hào)位點(diǎn)),但溶解氧最低點(diǎn)不是產(chǎn)毒微囊藻基因豐富度最低點(diǎn).這與李大命等在巢湖夏季和冬季有毒微囊藻和無(wú)毒微囊藻種群豐度研究中的實(shí)驗(yàn)研究結(jié)果[39]出現(xiàn)了差異性,具體原因有待進(jìn)一步研究.

    3.3 水體理化指標(biāo)與微囊藻產(chǎn)毒基因及毒素的相關(guān)性分析

    表4 環(huán)境因子與微囊藻和有毒微囊藻豐度之間的相關(guān)系數(shù)(n=8)

    注: * :<0.05;**:<0.01.

    已有研究表明,水體理化因子的差異可以顯著的影響微囊藻毒素的釋放.其中觀點(diǎn)主要集中在溫度、總氮、總磷、磷酸鹽等[5,40-41]因素的影響,在本次研究中我們發(fā)現(xiàn)海河水體中微囊藻無(wú)論從總基因豐度到產(chǎn)毒基因豐度都與pH值顯著相關(guān)(<0.05),這與本文在海河的另一項(xiàng)基于高通量測(cè)序得到的結(jié)果是一致的[42].這說(shuō)明在高pH值情況下微囊藻無(wú)論在絕對(duì)豐度還是與其他原核生物共存的相對(duì)豐度水平上均能占據(jù)較大比例.在所有環(huán)境因子中,影響海河中微囊藻的豐度差異最重要的原因是pH值.然而目前的研究并不能給出pH值與微囊藻直接相關(guān)的明確原因,二者之間的因果關(guān)系尚不明晰.有一種推測(cè)是在微囊藻固碳過(guò)程中有一個(gè)中間形式為碳酸氫鹽[43].藍(lán)藻的集中暴發(fā)可能使碳酸氫根離子在水體中大量富集而后隨著藻細(xì)胞衰敗破裂而釋放到水體中.微囊藻產(chǎn)毒基因、產(chǎn)毒藻株比例與各點(diǎn)位實(shí)際觀測(cè)到的微囊藻毒素含量并沒有顯著相關(guān)性.這證明了雖然基因簇在微囊藻毒素合成中必不可少[44],但毒素的實(shí)際合成過(guò)程中受到轉(zhuǎn)錄因子調(diào)控或基因修飾影響顯著,造成產(chǎn)毒藻株的豐度并不與毒素含量有實(shí)際關(guān)聯(lián).微囊藻毒素合成基因與微囊藻總16S rRNA基因豐度極顯著相關(guān)(<0.01),這證明了在海河中微囊藻的產(chǎn)毒比例趨近一致,各采樣點(diǎn)位產(chǎn)毒株所占比例大致相同.單位水體內(nèi)胞外微囊藻毒素含量與胞內(nèi)微囊藻毒素含量極顯著相關(guān),這說(shuō)明在采樣期間微囊藻毒素釋放速率較為穩(wěn)定,并沒有發(fā)生大規(guī)模的胞內(nèi)毒素釋放事件.但值得注意的是,雖然藻細(xì)胞外水體中微囊藻毒素含量在各點(diǎn)位均未超過(guò)WHO規(guī)定的飲用水微囊藻毒素含量標(biāo)準(zhǔn)但是計(jì)算相同體積水體中藻細(xì)胞的毒素含量遠(yuǎn)超這一標(biāo)準(zhǔn).這說(shuō)明海河微囊藻具有較高的毒素釋放潛力,這意味著當(dāng)水體發(fā)生大規(guī)模藻細(xì)胞破裂事件時(shí),水環(huán)境中毒素含量會(huì)遠(yuǎn)超過(guò)現(xiàn)有水平.一旦被飲用會(huì)對(duì)人體健康造成巨大威脅.在海河及與其類似的通海河道水體中,滲透壓造成藻細(xì)胞破裂也可能會(huì)給河口生態(tài)系統(tǒng)造成巨大破壞.因此對(duì)于海河中微囊藻的控制是十分有必要的.

    4 結(jié)論

    4.1 夏季海河天津市區(qū)段微囊藻種群豐度具有明顯空間差異產(chǎn)毒微囊藻種群豐度為1.16′104~2.48′107copies/mL,占總微囊藻種群的4.25%~28.59%.

    4.2 藻毒素含量最高點(diǎn)為8號(hào)采樣點(diǎn)天津站,每升水體中藻細(xì)胞共含毒素195.51μg,除去藻細(xì)胞水中毒素濃度為0.97μg/L.水體中微囊藻毒素未達(dá)到WHO規(guī)定的閾值.

    4.3 產(chǎn)毒基因豐度及比例與環(huán)境中微囊藻毒素含量未呈現(xiàn)顯著相關(guān)關(guān)系.pH值的大小與微囊藻總豐度和產(chǎn)毒基因型豐度呈現(xiàn)顯著相關(guān)關(guān)系.

    4.4 夏季海河中微囊藻胞內(nèi)毒素含量遠(yuǎn)高于胞外含量,可能對(duì)下游的河口區(qū)域生態(tài)系統(tǒng)安全構(gòu)成潛在威脅.

    [1] Pouria S, Andrade A D, Barbosa J, et al. Fatal microcystin intoxication in haemodialysis unit in Caruaru, Brazil. [J]. Lancet, 1998,352(9121): 21-26.

    [2] Shun-Zhang Y U. Primary prevention of hepatocellular carcinoma [J]. Journal of Gastroenterology & Hepatology, 2013,95(4):427–434.

    [3] Organization W H. Cyanobacterial toxins : Microcystin-LR [J]. Guidelines for Drinking-Water Quality, 1998,2:83-127.

    [4] Kurmayer R, Kutzenberger T. Application of real-time PCR for quantification of microcystin genotypes in a population of the toxic cyanobacterium Microcystis sp [J]. Applied & Environmental Microbiology, 2003,69(11):6723-6730.

    [5] Davis T W, Berry D L, Boyer G L, et al. The effects of temperature and nutrients on the growth and dynamics of toxic and non-toxic strains of Microcystis during cyanobacteria blooms [J]. Harmful Algae, 2009,8(5):715-725.

    [6] Joung S H, Oh H M, Ko S R, et al. Correlations between environmental factors and toxic and non-toxic Microcystis dynamics during bloom in Daechung Reservoir, Korea [J]. Harmful Algae, 2011,10(2):188-193.

    [7] Dittmann E, Neilan B A, Erhard M, et al. Insertional mutagenesis of a peptide synthetase gene that is responsible for hepatotoxin production in the cyanobacterium Microcystis aeruginosa PCC 7806 [J]. Molecular Microbiology, 1997,26(4):779-787.

    [8] Tillett D, Dittmann E, Erhard M, et al. Structural organization of microcystin biosynthesis in Microcystis aeruginosa PCC7806: an integrated peptide-polyketide synthetase system. [J]. Chemistry & Biology, 2000,7(10):753-764.

    [9] Baker J A, Neilan B A, Entsch B, et al. Identification of cyanobacteria and their toxigenicity in environmental samples by rapid molecular analysis.[J]. Environmental Toxicology, 2001,16(6):472-482.

    [10] Hui P, Song L, Liu Y, et al. Detection of hepatotoxic Microcystis strains by PCR with intact cells from both culture and environmental samples. [J]. Archives of Microbiology, 2002,178(6):421-427.

    [11] Ouahid Y, Pérez-Silva G, Campo F F D. Identification of potentially toxic environmental Microcystis by individual and multiple PCR amplification of specific microcystin synthetase gene regions [J]. Environmental Toxicology, 2005,20(3):235-242.

    [12] Rinta-Kanto J, Ouellette A, Boyer G, et al. Quantification of toxic Microcystis spp during the 2003and 2004blooms in western lake Erie using Quantitative Real-Time PCR [J]. Environmental Science & Technology, 2005,39(11):4198-4205.

    [13] 李大命,葉琳琳,于 洋,等.太湖水華期間有毒和無(wú)毒微囊藻種群豐度的動(dòng)態(tài)變化[J]. 生態(tài)學(xué)報(bào), 2012,32(22):7109-7116.

    [14] 李大命,孔繁翔,張 民,等.太湖和巢湖夏季藍(lán)藻水華期間產(chǎn)毒微囊藻和非產(chǎn)毒微囊藻種群豐度的空間分布[J]. 應(yīng)用與環(huán)境生物學(xué)報(bào), 2011(4):480-485.

    [15] Li D, Kong F, Shi X, et al. Quantification of microcystin-producing and non-microcystin producing Microcystis populations during the 2009and 2010blooms in Lake Taihu using quantitative real-time PCR [J]. Journal of Environmental Sciences, 2012,24(2):284-90.

    [16] Liu Y, Peiru H U, Sisan M A, et al. Detection of microcystin- producing Microcystis cells at the entrance of rivers to southern Lake Taihu by real-time fluorescence quantitative PCR [J]. Journal of Lake Sciences, 2016,28(2):246-252.

    [17] Walls J T, Wyatt K H, Doll J C, et al. Hot and toxic: Temperature regulates microcystin release from cyanobacteria [J]. Science of the Total Environment, 2017,610-611:786.

    [18] Mantzouki E, Lürling M, Fastner J, et al. Temperature Effects Explain Continental Scale Distribution of Cyanobacterial Toxins [J]. Toxins, 2018,10(4):156.

    [19] Kaebernick M, Neilan B A, B?rner T, et al. Light and the Transcriptional Response of the Microcystin Biosynthesis Gene Cluster [J]. Applied & Environmental Microbiology, 2000,66(8):3387- 3392.

    [20] Neilan B A, Jacobs D, Del D T, et al. rRNA sequences and evolutionary relationships among toxic and nontoxic cyanobacteria of the genus Microcystis [J]. Int. J. Syst. Bacteriol., 1997,47(3):693-697.

    [21] 于洪利,杜慶有.天津海河2013年藍(lán)藻暴發(fā)過(guò)程及相關(guān)因素分析[J]. 海河水利, 2017,s1:30-33.

    [22] Baxa D V, Kurobe T, Ger K A, et al. Estimating the abundance of toxic Microcystis in the San Francisco Estuary using quantitative real-time PCR [J]. Harmful Algae, 2010,9(3):342-349.

    [23] 謝數(shù)濤,龍思思,韓博平,等.微囊藻毒素合成酶基因的PCR檢測(cè)方法 [J]. 熱帶亞熱帶植物學(xué)報(bào), 2005,13(4):315-318.

    [24] Fergusson K M, Saint C P. Molecular phylogeny of Anabaena circinalis and its identification in environmental samples by PCR. [J]. Applied & Environmental Microbiology, 2000,66(9):4145-4148.

    [25] Baker J A, Entsch BNeilan B A, Mckay D B. Monitoring changing toxigenicity of a cyanobacterial bloom by molecular methods [J]. Applied & Environmental Microbiology, 2002,68(12):6070-6076.

    [26] Tooming-Klunderud A, Mikalsen B, Kristensen T, et al. The mosaic structure of the mcyABC operon in Microcystis [J]. Microbiology, 2008,154(7):1886-1899.

    [27] Tillett D, Parker D L, Neilan B A. Detection of toxigenicity by a probe for the microcystin synthetase A gene(mcyA) of the cyanobacterial genus Microcystis, comparison of toxicitieswith 16S rRNA and phycocyanin operon (phycocyanin intergenic spacer) phylogenies [J]. Applied & Environmental Microbiology, 2001,67(6):2810–2818.

    [28] Bittencourt-Oliveira M D C. Detection of potential microcystin- producing in Brazilian reservoirs with a mcyB molecular marker [J]. Harmful Algae, 2003,2(1):51-60.

    [29] Kurmayer R, Dittmann E, Fastner J, et al. Diversity of microcystin genes within a population of the toxic cyanobacterium Microcystis spp. in Lake Wannsee (Berlin, Germany). [J]. Microbial Ecology, 2002, 43(1):107-118.

    [30] Jaana V, Anne R, Katrianna H, et al. Quantitative real-time PCR for determination of microcystin synthetase e copy numbers for microcystis and anabaena in lakes [J]. Applied & Environmental Microbiology, 2003,69(12):7289-97.

    [31] Youness O, Gonzalo P S, Campo F F D. Identification of potentially toxic environmental Microcystis by individual and multiple PCR amplification of specific microcystin synthetase gene regions [J]. Environmental Toxicology, 2005,20(3):235-242.

    [32] Oberholster P J, Myburgh J G, Govender D, et al. Identification of toxigenic Microcystis strains after incidents of wild animal mortalities in the Kruger National Park, South Africa [J]. Ecotoxicology & Environmental Safety, 2009,72(4):1177-1182.

    [33] Hotto A, Satchwell M, Boyer G. Molecular characterization of potential microcystin-producing cyanobacteria in lake Ontario embayments and nearshore waters [J]. Applied & Environmental Microbiology, 2007,73(14):4570-4578.

    [34] Ye C, Yang Y, Xu Q, et al. Volatile organic compound emissions from Microcystis aeruginosa under different phosphorus sources and concentrations [J]. Phycological Research, 2018,66:15-22.

    [35] Yoshida M, Yoshida T, Takashima Y, et al. Dynamics of microcystin- producing and non-microcystin-producing Microcystis populations is correlated with nitrate concentration in a Japanese lake [J]. Fems. Microbiology Letters, 2007,266(1):49-53.

    [36] Beversdorf L J, Chaston S D, Miller T R, et al. Microcystin mcyA and mcyE Gene Abundances Are Not Appropriate Indicators of Microcystin Concentrations in Lakes [J]. Plos. One, 2015,10(5): e0125353.

    [37] Hotto A M, Satchwell M F, Berry D L, et al. Spatial and temporal diversity of microcystins and microcystin-producing genotypes in Oneida Lake, NY [J]. Harmful Algae, 2008,7(5):671-681.

    [38] 劉永定,范 曉,胡征宇.中國(guó)藻類學(xué)研究 [M]. 武漢:武漢出版社, 2001:245.

    [39] 李大命,孔繁翔,陽(yáng) 振,等.巢湖夏季和冬季有毒微囊藻和無(wú)毒微囊藻種群豐度研究[J]. 環(huán)境科學(xué)學(xué)報(bào), 2011,12:2672-2680.

    [40] Vézie C, Rapala J, Vaitomaa J, et al. Effect of nitrogen and phosphorus on growth of toxic and nontoxic Microcystis strains and on intracellular microcystins concentrations [J]. Microbial Ecology, 2002, 43(4):443-454.

    [41] Rintakanto J M, Konopko E A, Debruyn J M, et al. Lake Erie Microcystis: relationship between microcystin production, dynamics of genotypes and environmental parameters in a large lake [J]. Harmful Algae, 2009,8(5):665-673.

    [42] Huo D, Chen Y X, Liu P, et al. Molecular detection of microbial communities associated with Microcystis vs Synechococcus dominated waters in Tianjin, China [J] Journal of oceanology and limnology, 2018,4(36):1145-1156.

    [43] Sandrini G, Ji X, Verspagen J M, et al. Rapid adaptation of harmful cyanobacteria to rising CO2[J]. Proc. Natl. Acad. Sci., 2016,113(33): 9315-9320.

    [44] Rantala A, Fewer D P, Hisbergues M, et al. Phylogenetic Evidence for the Early Evolution of Microcystin Synthesis [J]. Proceedings of the National Academy of Sciences of the United States of America, 2004,101(2):568-73.

    [45] 楊曉紅,蒲朝文,張仁平,等.水體微囊藻毒素污染對(duì)人群的非致癌健康風(fēng)險(xiǎn) [J]. 中國(guó)環(huán)境科學(xué), 2013,33(1):181-185.

    致謝:感謝天津農(nóng)學(xué)院基礎(chǔ)學(xué)院的劉萍老師在水化學(xué)分析中提供的幫助.

    The spatial distribution ofand microcystin and its relationship with environmental factors in Haihe Tianjin City.

    HUO Da1,2, CAO Qi1, WANG Su-yan1, CHEN Pei-pei1, LI Yi-ming1, QIAO Zhi-yi1,2*

    (1.Department of Fishery Science, Tianjing Agricultural College, Tianjing 300384, China;2.Tianjin Key Laboratory of Aqua-ecology and Aqua-culture, Tianjin 300384, China)., 2018,38(10):3897~3903

    The spatial distribution ofand microcystin in Haihe River during water bloom were investigated using quantitative real-time PCR (qPCR) and ELISA techniques. Microcystin synthetase gene () and-fragments were used as the target the total microcystin population and the potential microcystin-producing subpopulation. The results showed that microcystinexisted as a mixed population of potentially toxic and non-toxic genotypes in Haihe River. There was significantly spatial variation in the abundance of microcystin-producingThe abundance of microcystin-producingsubpopulation was ranged from 1.16×104to 2.48×107copies/mL and the ratio of its abundance to total microcystinvaried from 4.25% to 28.59%. The peakvalue of microcystinwas site point 8, which 1mL water contained 195.51μg microcystininalgae cells and 0.97μg/L in water. To sum up, total abundance ofremained a high level in Haihe Tianjin region. The toxin ofshowed various among different sampling sites. The most significant impact factors forabundance was pH value. However, the content of microcystin is not upon the value of WHO’s threshold. But the intercellular toxin in a given volume of water remains a high level, which may cause the potential risk to the ecosystem of down flow estuary.

    ;microcystin;quantitative real-time PCR;ELISA techniques;Haihe River

    X522

    A

    1000-6923(2018)10-3897-07

    霍 達(dá)(1993-),男,內(nèi)蒙古呼倫貝爾人,天津農(nóng)學(xué)院碩士研究生,主要從事分子生態(tài)學(xué)方向的研究.發(fā)表論文8篇.

    2018-03-23

    國(guó)家國(guó)際科技合作專項(xiàng)(2013DFA71340)

    * 責(zé)任作者, 副教授, qiaozhiyi7070@163.com

    猜你喜歡
    產(chǎn)毒海河微囊
    《海河水利》2021年總目次
    海河水利(2021年6期)2021-12-28 11:10:18
    一種臨近海河的超深地下室新型防水體系
    霉菌毒素什么條件下易產(chǎn)生
    孜然精油對(duì)產(chǎn)毒黃曲霉的抑制活性研究
    ε-聚賴氨酸鹽酸鹽對(duì)擴(kuò)展青霉生長(zhǎng)和產(chǎn)毒的影響
    破冰開河
    支點(diǎn)(2016年3期)2016-03-21 13:01:12
    微囊懸浮-懸浮劑和微囊懸浮劑不是同種劑型
    食品防腐劑對(duì)產(chǎn)毒黃曲霉孢子萌發(fā)及菌絲生長(zhǎng)的抑制作用研究
    微囊藻毒素-LR對(duì)秀麗線蟲精子形成的毒性作用
    馬錢子生物堿微囊的制備及評(píng)價(jià)
    а√天堂www在线а√下载| 在线a可以看的网站| 久9热在线精品视频| 欧美黑人欧美精品刺激| 91久久精品电影网| 一本综合久久免费| 尤物成人国产欧美一区二区三区| 国产成人aa在线观看| 国产黄色小视频在线观看| 亚洲精品色激情综合| 变态另类成人亚洲欧美熟女| 在线观看午夜福利视频| 亚洲精品日韩av片在线观看 | 国产综合懂色| 九色国产91popny在线| 欧美色视频一区免费| 亚洲无线在线观看| 97碰自拍视频| 在线观看舔阴道视频| 婷婷丁香在线五月| 久久九九热精品免费| 国内毛片毛片毛片毛片毛片| 国产一区二区在线观看日韩 | 亚洲中文日韩欧美视频| 在线观看午夜福利视频| 日韩欧美三级三区| 好男人电影高清在线观看| 特级一级黄色大片| 亚洲国产精品久久男人天堂| 午夜精品久久久久久毛片777| 亚洲,欧美精品.| 久久久久久九九精品二区国产| 一个人看视频在线观看www免费 | 国产高清视频在线播放一区| 日本撒尿小便嘘嘘汇集6| 成人无遮挡网站| 九九热线精品视视频播放| 88av欧美| 国内揄拍国产精品人妻在线| 日韩欧美 国产精品| 99久久99久久久精品蜜桃| 婷婷丁香在线五月| 久久久色成人| 一夜夜www| 久久久久久久久久黄片| 搡老岳熟女国产| 午夜福利高清视频| 亚洲专区国产一区二区| 无限看片的www在线观看| 精品久久久久久久人妻蜜臀av| 观看免费一级毛片| 亚洲久久久久久中文字幕| 夜夜夜夜夜久久久久| 精品一区二区三区视频在线 | 天堂网av新在线| 成年女人毛片免费观看观看9| 国产成人aa在线观看| 成年人黄色毛片网站| 精品久久久久久久末码| 男女之事视频高清在线观看| 亚洲精品粉嫩美女一区| 亚洲激情在线av| 91av网一区二区| 女人高潮潮喷娇喘18禁视频| 国产极品精品免费视频能看的| 久久精品国产亚洲av涩爱 | 国产成人系列免费观看| 国内少妇人妻偷人精品xxx网站| 国内精品一区二区在线观看| 91在线观看av| 亚洲av成人精品一区久久| 亚洲av成人精品一区久久| 日本一本二区三区精品| 日韩欧美在线二视频| 免费在线观看成人毛片| 1024手机看黄色片| 嫁个100分男人电影在线观看| 日韩欧美免费精品| 少妇的逼好多水| 深夜精品福利| 婷婷六月久久综合丁香| 午夜a级毛片| 欧美在线黄色| 亚洲精品成人久久久久久| 99国产极品粉嫩在线观看| 一进一出抽搐gif免费好疼| 亚洲中文字幕一区二区三区有码在线看| 成人特级黄色片久久久久久久| 亚洲精品一卡2卡三卡4卡5卡| 免费av观看视频| 国产av麻豆久久久久久久| 在线播放国产精品三级| 人妻久久中文字幕网| 精品无人区乱码1区二区| 日韩有码中文字幕| 久久九九热精品免费| 国产亚洲av嫩草精品影院| 亚洲av免费高清在线观看| 日本 欧美在线| 深爱激情五月婷婷| a级毛片a级免费在线| 18禁在线播放成人免费| 91在线精品国自产拍蜜月 | 婷婷丁香在线五月| 又爽又黄无遮挡网站| 国产中年淑女户外野战色| 免费看十八禁软件| 久久久久国内视频| 色精品久久人妻99蜜桃| 亚洲黑人精品在线| 亚洲精华国产精华精| 午夜福利在线观看吧| 成人无遮挡网站| 欧美成人一区二区免费高清观看| 岛国在线免费视频观看| 欧美性感艳星| 99国产综合亚洲精品| 噜噜噜噜噜久久久久久91| 久久久久久久午夜电影| 国产亚洲欧美在线一区二区| 床上黄色一级片| 天堂√8在线中文| 宅男免费午夜| 免费电影在线观看免费观看| 最新中文字幕久久久久| 国产午夜精品久久久久久一区二区三区 | www.色视频.com| 动漫黄色视频在线观看| 午夜两性在线视频| 欧美三级亚洲精品| 亚洲 欧美 日韩 在线 免费| 少妇熟女aⅴ在线视频| 亚洲成av人片在线播放无| 精品久久久久久久久久免费视频| 97超级碰碰碰精品色视频在线观看| aaaaa片日本免费| 亚洲一区高清亚洲精品| 欧美中文日本在线观看视频| 久久久国产成人精品二区| www国产在线视频色| 又爽又黄无遮挡网站| 精品午夜福利视频在线观看一区| 欧美另类亚洲清纯唯美| 人妻久久中文字幕网| 国产v大片淫在线免费观看| 亚洲欧美激情综合另类| 精品熟女少妇八av免费久了| 青草久久国产| 久久午夜亚洲精品久久| 午夜免费激情av| 亚洲五月天丁香| 噜噜噜噜噜久久久久久91| 嫩草影院精品99| 香蕉丝袜av| 久久久久久久精品吃奶| 午夜福利在线观看吧| 成人无遮挡网站| 美女cb高潮喷水在线观看| 真人一进一出gif抽搐免费| 精品熟女少妇八av免费久了| 国产麻豆成人av免费视频| 女警被强在线播放| 免费在线观看亚洲国产| 亚洲欧美精品综合久久99| 色av中文字幕| www国产在线视频色| 狂野欧美激情性xxxx| 熟女少妇亚洲综合色aaa.| 欧洲精品卡2卡3卡4卡5卡区| 亚洲人成网站高清观看| 成人欧美大片| 天堂av国产一区二区熟女人妻| 日本熟妇午夜| 国产精品久久久久久久电影 | 亚洲五月天丁香| 好男人在线观看高清免费视频| 国产日本99.免费观看| 色av中文字幕| 免费在线观看日本一区| 好看av亚洲va欧美ⅴa在| 少妇高潮的动态图| 亚洲色图av天堂| 全区人妻精品视频| 1000部很黄的大片| 最好的美女福利视频网| 国产精品久久久人人做人人爽| 黄片小视频在线播放| 90打野战视频偷拍视频| 波多野结衣高清作品| netflix在线观看网站| 好看av亚洲va欧美ⅴa在| 午夜两性在线视频| 欧美+亚洲+日韩+国产| 久久精品国产亚洲av涩爱 | 老汉色av国产亚洲站长工具| 欧美日韩中文字幕国产精品一区二区三区| 欧美一区二区国产精品久久精品| 欧美性感艳星| 午夜福利在线观看吧| 午夜精品久久久久久毛片777| 最新在线观看一区二区三区| 高清在线国产一区| 色综合婷婷激情| 日本熟妇午夜| 精品国产亚洲在线| 亚洲成人久久性| 精华霜和精华液先用哪个| 色吧在线观看| 舔av片在线| 欧美在线一区亚洲| 亚洲成av人片免费观看| e午夜精品久久久久久久| 一本精品99久久精品77| 丁香六月欧美| 免费无遮挡裸体视频| 一进一出抽搐gif免费好疼| 日本与韩国留学比较| 欧美av亚洲av综合av国产av| 亚洲av美国av| av在线天堂中文字幕| 亚洲精品久久国产高清桃花| 亚洲精品日韩av片在线观看 | 久久精品综合一区二区三区| 丝袜美腿在线中文| 国产精品嫩草影院av在线观看 | 真实男女啪啪啪动态图| 在线播放无遮挡| 国产乱人伦免费视频| 国产乱人视频| 国产伦精品一区二区三区视频9 | 青草久久国产| 亚洲成a人片在线一区二区| 精品久久久久久久久久久久久| 欧美高清成人免费视频www| 国产精品亚洲美女久久久| 少妇高潮的动态图| 三级国产精品欧美在线观看| 国产一区在线观看成人免费| 天堂动漫精品| 51午夜福利影视在线观看| 一本一本综合久久| 一进一出抽搐动态| 国产精品久久电影中文字幕| 久久久久久久亚洲中文字幕 | 一本精品99久久精品77| 9191精品国产免费久久| 在线a可以看的网站| 人人妻人人澡欧美一区二区| 搡老岳熟女国产| 99在线人妻在线中文字幕| 亚洲精品粉嫩美女一区| 一个人观看的视频www高清免费观看| 国产亚洲av嫩草精品影院| 亚洲精品456在线播放app | 国产免费av片在线观看野外av| 亚洲无线在线观看| 午夜免费激情av| 国产单亲对白刺激| 高潮久久久久久久久久久不卡| 日韩国内少妇激情av| 日韩免费av在线播放| 18禁国产床啪视频网站| 免费大片18禁| 国产美女午夜福利| 天堂网av新在线| 一个人看的www免费观看视频| 免费看a级黄色片| 久久婷婷人人爽人人干人人爱| 90打野战视频偷拍视频| 一区二区三区免费毛片| 欧美黄色淫秽网站| 脱女人内裤的视频| 中文亚洲av片在线观看爽| 一级a爱片免费观看的视频| a级毛片a级免费在线| 日韩亚洲欧美综合| 高清在线国产一区| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利高清视频| 日韩欧美三级三区| 老鸭窝网址在线观看| 女人被狂操c到高潮| 一区二区三区国产精品乱码| 日韩中文字幕欧美一区二区| 亚洲片人在线观看| 一个人免费在线观看电影| 狂野欧美白嫩少妇大欣赏| 国产精品一区二区免费欧美| 亚洲av二区三区四区| 亚洲国产精品999在线| 天堂√8在线中文| 午夜福利在线观看吧| 亚洲精品456在线播放app | 日本黄大片高清| 女人高潮潮喷娇喘18禁视频| 91麻豆精品激情在线观看国产| 舔av片在线| 好男人在线观看高清免费视频| 少妇的丰满在线观看| 日本熟妇午夜| 亚洲国产精品999在线| 免费高清视频大片| www.www免费av| 老司机午夜福利在线观看视频| 精品久久久久久久久久免费视频| 国产探花极品一区二区| 精品一区二区三区视频在线观看免费| 国产野战对白在线观看| 深夜精品福利| 亚洲精品久久国产高清桃花| 成年免费大片在线观看| 亚洲欧美日韩高清在线视频| 成人国产综合亚洲| 亚洲乱码一区二区免费版| 色吧在线观看| www.www免费av| 亚洲国产精品久久男人天堂| 特级一级黄色大片| 日韩人妻高清精品专区| 99在线人妻在线中文字幕| a在线观看视频网站| 亚洲精品乱码久久久v下载方式 | 亚洲,欧美精品.| 成人午夜高清在线视频| 久久久国产成人精品二区| 成人欧美大片| 悠悠久久av| 国产成年人精品一区二区| 午夜福利在线观看吧| 丰满人妻熟妇乱又伦精品不卡| 老司机在亚洲福利影院| 人人妻人人澡欧美一区二区| 久久香蕉国产精品| 男人和女人高潮做爰伦理| 亚洲在线观看片| 我的老师免费观看完整版| 97超视频在线观看视频| 久久久久九九精品影院| 黑人欧美特级aaaaaa片| 全区人妻精品视频| 99热只有精品国产| 亚洲最大成人手机在线| 久久国产精品人妻蜜桃| 国产一级毛片七仙女欲春2| 88av欧美| 国产亚洲精品久久久com| 亚洲av成人精品一区久久| 成年女人看的毛片在线观看| 美女黄网站色视频| 午夜老司机福利剧场| 一本精品99久久精品77| 小说图片视频综合网站| 久久99热这里只有精品18| 日韩欧美精品免费久久 | 成年免费大片在线观看| 成人永久免费在线观看视频| 无遮挡黄片免费观看| 免费人成视频x8x8入口观看| 校园春色视频在线观看| 成年版毛片免费区| 免费av毛片视频| 1024手机看黄色片| 国产欧美日韩精品一区二区| 日韩人妻高清精品专区| 亚洲国产欧洲综合997久久,| 日韩欧美精品免费久久 | 少妇人妻精品综合一区二区 | 国产精品乱码一区二三区的特点| 又粗又爽又猛毛片免费看| 久久香蕉国产精品| 亚洲av免费高清在线观看| 欧洲精品卡2卡3卡4卡5卡区| 在线十欧美十亚洲十日本专区| 国产乱人视频| 久久久国产精品麻豆| 国产视频内射| 亚洲国产欧美网| 久久国产乱子伦精品免费另类| 叶爱在线成人免费视频播放| 欧美中文日本在线观看视频| 又黄又粗又硬又大视频| 日韩 欧美 亚洲 中文字幕| 最近最新中文字幕大全免费视频| 法律面前人人平等表现在哪些方面| 国产乱人伦免费视频| 国产一区二区在线观看日韩 | 国产美女午夜福利| 脱女人内裤的视频| 国产高清视频在线播放一区| 变态另类丝袜制服| 在线观看日韩欧美| 国产激情欧美一区二区| 99久久久亚洲精品蜜臀av| 国产成人福利小说| 亚洲欧美激情综合另类| 少妇高潮的动态图| 深夜精品福利| 老司机午夜十八禁免费视频| 天天一区二区日本电影三级| 婷婷六月久久综合丁香| 99精品久久久久人妻精品| 中文字幕高清在线视频| 五月伊人婷婷丁香| 黄片大片在线免费观看| 国产激情欧美一区二区| av天堂在线播放| 欧美黑人巨大hd| 啦啦啦韩国在线观看视频| av中文乱码字幕在线| 岛国在线免费视频观看| 国产成人av激情在线播放| 国产三级在线视频| 最后的刺客免费高清国语| 天堂av国产一区二区熟女人妻| 日本与韩国留学比较| 日本 欧美在线| 脱女人内裤的视频| 国产成人系列免费观看| 桃红色精品国产亚洲av| 看黄色毛片网站| 久久这里只有精品中国| 成年人黄色毛片网站| 看黄色毛片网站| 搡女人真爽免费视频火全软件 | 午夜免费观看网址| 亚洲最大成人手机在线| 亚洲中文日韩欧美视频| 在线观看一区二区三区| 午夜精品在线福利| 国内精品一区二区在线观看| 亚洲国产欧洲综合997久久,| 桃色一区二区三区在线观看| 精品不卡国产一区二区三区| 哪里可以看免费的av片| 老司机福利观看| 亚洲精品在线美女| 欧美性感艳星| www.熟女人妻精品国产| 人人妻,人人澡人人爽秒播| 日本 欧美在线| 脱女人内裤的视频| 国产一区在线观看成人免费| 欧美色视频一区免费| 特级一级黄色大片| 国产欧美日韩一区二区三| 亚洲人成伊人成综合网2020| 网址你懂的国产日韩在线| 麻豆国产97在线/欧美| 亚洲avbb在线观看| 在线观看av片永久免费下载| 欧美黑人巨大hd| 亚洲国产欧美网| 欧美午夜高清在线| 免费在线观看亚洲国产| 99视频精品全部免费 在线| 欧美一区二区国产精品久久精品| 特大巨黑吊av在线直播| 老汉色∧v一级毛片| 欧美日韩乱码在线| 精品国内亚洲2022精品成人| 欧美色欧美亚洲另类二区| 国产激情欧美一区二区| 久久精品国产自在天天线| 操出白浆在线播放| 成人国产综合亚洲| 国产三级黄色录像| 色尼玛亚洲综合影院| 97碰自拍视频| 高清毛片免费观看视频网站| 国产午夜精品论理片| 国产三级黄色录像| 国产精品 欧美亚洲| 91字幕亚洲| 桃红色精品国产亚洲av| 搡老熟女国产l中国老女人| 国产成年人精品一区二区| 国产69精品久久久久777片| 欧美日韩瑟瑟在线播放| 高清日韩中文字幕在线| 毛片女人毛片| 久久国产精品人妻蜜桃| 久久久精品大字幕| 国产精品一区二区三区四区免费观看 | 老司机午夜十八禁免费视频| 国产探花极品一区二区| 国产欧美日韩精品亚洲av| 人妻久久中文字幕网| 国产爱豆传媒在线观看| 久久6这里有精品| 真人做人爱边吃奶动态| 男女做爰动态图高潮gif福利片| 亚洲人成电影免费在线| 在线观看免费视频日本深夜| 香蕉丝袜av| 亚洲在线自拍视频| 3wmmmm亚洲av在线观看| 一进一出抽搐动态| 成人特级黄色片久久久久久久| 超碰av人人做人人爽久久 | 亚洲精品一卡2卡三卡4卡5卡| 淫秽高清视频在线观看| 国产 一区 欧美 日韩| 97人妻精品一区二区三区麻豆| 午夜精品久久久久久毛片777| 一个人看的www免费观看视频| 国产精品自产拍在线观看55亚洲| 观看免费一级毛片| 一区福利在线观看| 中文字幕精品亚洲无线码一区| 国产精品亚洲一级av第二区| 国内毛片毛片毛片毛片毛片| 看片在线看免费视频| 久久婷婷人人爽人人干人人爱| 最新在线观看一区二区三区| 国内少妇人妻偷人精品xxx网站| 国模一区二区三区四区视频| 内射极品少妇av片p| 嫩草影院精品99| 夜夜躁狠狠躁天天躁| 亚洲片人在线观看| 久久精品国产99精品国产亚洲性色| 亚洲 国产 在线| 国产亚洲av嫩草精品影院| 一个人看的www免费观看视频| 亚洲av一区综合| 国产aⅴ精品一区二区三区波| 3wmmmm亚洲av在线观看| 美女高潮的动态| 一进一出抽搐gif免费好疼| 黄片大片在线免费观看| 国产精品国产高清国产av| 日本免费a在线| 欧美乱色亚洲激情| 亚洲性夜色夜夜综合| 国产欧美日韩精品一区二区| 99在线人妻在线中文字幕| 欧美国产日韩亚洲一区| 国产真人三级小视频在线观看| 两人在一起打扑克的视频| 国产三级中文精品| 欧美日韩一级在线毛片| 亚洲一区二区三区色噜噜| 国产亚洲av嫩草精品影院| 在线天堂最新版资源| 色精品久久人妻99蜜桃| 亚洲av五月六月丁香网| 狠狠狠狠99中文字幕| 内射极品少妇av片p| 极品教师在线免费播放| 免费av观看视频| 午夜福利免费观看在线| 国产欧美日韩精品一区二区| 少妇的逼好多水| 看免费av毛片| 亚洲av成人av| 日日夜夜操网爽| 国产三级黄色录像| 国产色婷婷99| 中文字幕人妻丝袜一区二区| 精品一区二区三区av网在线观看| 欧美午夜高清在线| 亚洲第一电影网av| www.色视频.com| 国产在视频线在精品| 亚洲美女视频黄频| 亚洲精品乱码久久久v下载方式 | 麻豆成人午夜福利视频| 国内精品美女久久久久久| 一级作爱视频免费观看| 波多野结衣高清无吗| 亚洲人与动物交配视频| 波多野结衣高清无吗| 一级黄片播放器| 一进一出抽搐gif免费好疼| 久久久国产成人免费| 老熟妇仑乱视频hdxx| 欧美成狂野欧美在线观看| 不卡一级毛片| 真人一进一出gif抽搐免费| 日本黄大片高清| 麻豆成人午夜福利视频| 久久久久久九九精品二区国产| 日韩有码中文字幕| 中文字幕av成人在线电影| 51国产日韩欧美| 精品人妻1区二区| 日本撒尿小便嘘嘘汇集6| 欧美又色又爽又黄视频| av女优亚洲男人天堂| 国产精品综合久久久久久久免费| 亚洲国产色片| 中文字幕av成人在线电影| 国产成人av教育| 国产av在哪里看| 性欧美人与动物交配| 欧美在线黄色| 亚洲不卡免费看| 一进一出好大好爽视频| 国产在线精品亚洲第一网站| 日韩欧美三级三区| 看免费av毛片| 亚洲精品国产精品久久久不卡| 一区二区三区国产精品乱码| 精品久久久久久久毛片微露脸| 亚洲国产欧美网| 少妇裸体淫交视频免费看高清| 国产精品一区二区三区四区免费观看 | 国产99白浆流出| 久久久国产成人免费| 91字幕亚洲| 97超级碰碰碰精品色视频在线观看| 午夜影院日韩av| 精品午夜福利视频在线观看一区| 男女做爰动态图高潮gif福利片| 国内久久婷婷六月综合欲色啪|