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

    17β-雌二醇對(duì)巖溶稻田土壤微生物群落的影響

    2017-01-10 12:56:48黃雅丹靳振江李強(qiáng)劉暢

    黃雅丹,靳振江,李強(qiáng),劉暢

    1.桂林醫(yī)學(xué)院研究生學(xué)院,廣西 桂林 541004;2.中國地質(zhì)科學(xué)院巖溶地質(zhì)研究所,廣西 桂林 541004;

    3.聯(lián)合國教科文組織國際巖溶研究中心,廣西 桂林 541004;4.桂林理工大學(xué)環(huán)境科學(xué)與工程學(xué)院,廣西 桂林 541004

    17β-雌二醇對(duì)巖溶稻田土壤微生物群落的影響

    黃雅丹1,靳振江2,4,李強(qiáng)2,3*,劉暢2

    1.桂林醫(yī)學(xué)院研究生學(xué)院,廣西 桂林 541004;2.中國地質(zhì)科學(xué)院巖溶地質(zhì)研究所,廣西 桂林 541004;

    3.聯(lián)合國教科文組織國際巖溶研究中心,廣西 桂林 541004;4.桂林理工大學(xué)環(huán)境科學(xué)與工程學(xué)院,廣西 桂林 541004

    為明確17β-雌二醇對(duì)巖溶區(qū)稻田土壤生態(tài)系統(tǒng)微生物多樣性-穩(wěn)定性關(guān)系的影響,通過17β-雌二醇室內(nèi)添加模擬實(shí)驗(yàn),采用熒光定量聚合酶鏈反應(yīng)和聚合酶鏈反應(yīng)-變性梯度凝膠電泳法研究有氧和厭氧條件下細(xì)菌和真菌豐度及多樣性變化。結(jié)果表明,(1)0.05 mg?kg-1和0.20 mg?kg-117β-雌二醇對(duì)巖溶稻田土壤微生物豐度具有一定的表觀刺激效應(yīng),這種表觀刺激效應(yīng)在厭氧條件下尤為顯著。巖溶稻田土壤細(xì)菌和真菌總拷貝數(shù)及歸一化結(jié)果進(jìn)一步顯示,17β-雌二醇的添加量為0.05 mg?kg-1時(shí),巖溶稻田土壤微生物豐度最高,這可能是巖溶稻田土壤微生物對(duì)17β-雌二醇產(chǎn)生應(yīng)急反應(yīng)的一個(gè)重要應(yīng)答點(diǎn)。(2)在厭氧培養(yǎng)條件下,巖溶稻田土壤細(xì)菌和真菌均勻度指數(shù)與其豐度呈現(xiàn)相反的變化趨勢。盡管17β-雌二醇能夠促進(jìn)少數(shù)厭氧微生物種屬及其數(shù)量的增加,但整體上降低了巖溶稻田土壤微生物系統(tǒng)的多樣性和穩(wěn)定性。(3)培養(yǎng)條件與微生物群落的雙因素方差分析結(jié)果進(jìn)一步證實(shí),土壤通氣狀況(有氧和厭氧)是影響巖溶稻田土壤微生物群落穩(wěn)定性的極顯著因素。這說明,定期對(duì)巖溶稻田進(jìn)行翻耕、曬田或水-旱交替耕作能夠促進(jìn)巖溶稻田土壤中17β-雌二醇的降解,降低厭氧微生物活性,緩解17β-雌二醇對(duì)巖溶稻田土壤微生物帶來的生態(tài)風(fēng)險(xiǎn),從而為巖溶區(qū)稻作農(nóng)業(yè)生產(chǎn)管理提供技術(shù)支持。

    17β-雌二醇;土壤通氣狀況;微生物豐度;微生物多樣性

    具有雌激素效應(yīng)的17β-雌二醇能夠干擾生物體內(nèi)分泌物質(zhì)的合成、釋放、運(yùn)輸、結(jié)合、代謝等過程,以此激活或抑制生物體的內(nèi)分泌系統(tǒng)(Bechi et al.,2013)。加之17β-雌二醇在水中的溶解度很低,極性相對(duì)較低,極易被水體底泥和土壤吸附,且很難再移動(dòng),因此,即使17β-雌二醇質(zhì)量濃度低至ng?L-1級(jí)時(shí)仍對(duì)生態(tài)環(huán)境具有極大的風(fēng)險(xiǎn),為此《斯德哥爾摩公約》將其列為禁止使用和限制使用的持久性有機(jī)污染物(聯(lián)合國環(huán)境規(guī)劃署,2001)。

    17β-雌二醇廣泛存在于市政污水及家畜糞便中,直接或間接排入河流的17β-雌二醇將會(huì)對(duì)環(huán)境造成不同程度的污染。譬如,在湖北、黑龍江等畜禽養(yǎng)殖大省,部分河流中的17β-雌二醇質(zhì)量濃度高達(dá)29 ng?L-1(邵曉玲等,2008),相對(duì)而言,珠江水域沉積物中的17β-雌二醇質(zhì)量分?jǐn)?shù)變化范圍則在0.9~2.6 ng?g-1之間(Gong et al.,2011;龔劍等,2011)。在國外,17β-雌二醇對(duì)水環(huán)境造成的污染形勢也不可小覷,Kolpin et al.(2002)通過調(diào)查美國30個(gè)州139個(gè)地表水樣,發(fā)現(xiàn)將近40%的地表水受到17β-雌二醇不同程度的污染,而Kuch et al.(2001)在德國南部的一些飲用水中檢測到的17β-雌二醇質(zhì)量濃度為0.2~2.1 ng?L-1。大量毒理實(shí)驗(yàn)表明,環(huán)境中的17β-雌二醇能夠影響某些魚類的內(nèi)分泌系統(tǒng)及其生殖能力,甚至引起雄魚雌性化(陳家長等,2012),并且也能顯著抑制水體微生物的反硝化作用(吳偉恒等,2014)。此外,趙穎等(2014)12-14通過模擬實(shí)驗(yàn)證明當(dāng)17β-雌二醇在環(huán)境中的質(zhì)量濃度不超過0.5 ng?L-1時(shí),能提高水體產(chǎn)甲烷微生物的活性,導(dǎo)致水體甲烷釋放強(qiáng)度增加;而當(dāng)17β-雌二醇質(zhì)量濃度為1.0~10.0 ng?L-1時(shí)則對(duì)水體產(chǎn)甲烷微生物活性產(chǎn)生抑制效應(yīng)。

    稻作農(nóng)業(yè)是保證我國糧食安全的主要支撐,由于水稻整個(gè)生育期都需要維持一定的水位,因此在水資源分配不均勻的巖溶區(qū)不得不對(duì)稻田進(jìn)行污水灌溉(陳曉鋒等,2009),然而含17β-雌二醇污水對(duì)稻田土壤微生物生態(tài)系統(tǒng)的影響研究目前尚未涉及。為此,本實(shí)驗(yàn)在有氧及厭氧培養(yǎng)條件下,采用熒光定量聚合酶鏈反應(yīng)(Quantitative Real-time Polymerase Chain Reaction,qPCR)和聚合酶鏈反應(yīng)-變性梯度凝膠電泳(Polymerase Chain Reaction-Denaturing Gradient Gel Electrophoresis,PCR-DGGE)法研究不同添加量的17β-雌二醇對(duì)巖溶稻田土壤微生物群落的影響,以期為評(píng)價(jià)環(huán)境殘留的17β-雌二醇對(duì)巖溶稻田土壤微生物造成的生態(tài)風(fēng)險(xiǎn)提供基礎(chǔ)數(shù)據(jù)。

    1 材料與方法

    1.1 實(shí)驗(yàn)土壤

    培養(yǎng)實(shí)驗(yàn)所需稻田土壤采自桂林會(huì)仙巖溶濕地0~20 cm土層,利用總有機(jī)碳分析儀TOC-V(日本島津)測得土壤有機(jī)碳質(zhì)量分?jǐn)?shù)為22.33 mg?kg-1,利用inoLab? pH 7110實(shí)驗(yàn)室臺(tái)式pH測試儀(德國WTW)測得土壤pH為7.91。土樣采集期,該區(qū)晝夜氣溫21~30 ℃。采集的新鮮土樣過10 mm篩,用濕紗布蓋在土表于室溫下放置3 d,待微生物活化后,加入配好的17β-雌二醇溶液,添加到土壤中的17β-雌二醇分別為0(CK)、0.05、0.10、0.20、0.40 mg?kg-1,并用去離子水調(diào)節(jié)土壤濕度至40%,分別模擬稻田水、旱條件進(jìn)行有氧和厭氧(厭氧培養(yǎng)箱)處理,在25 ℃下恒溫培養(yǎng)30 d后進(jìn)行后續(xù)分析。同一處理?xiàng)l件設(shè)有3個(gè)平行樣。

    1.2 測定方法

    1.2.1 土壤微生物DNA的提取

    稱取0.25 g土壤樣品,采用美國MOBIO土壤強(qiáng)力DNA提取試劑盒提取土壤微生物總DNA,利用Quawell 5000微量紫外分光光度計(jì)測定其濃度和純度,保存于-80 ℃冰箱。

    1.2.2 熒光定量聚合酶鏈反應(yīng)

    為獲取16S rRNA基因拷貝數(shù),將提取的土壤微生物DNA作為qPCR模板,采用對(duì)大多數(shù)細(xì)菌具有特異性的引物對(duì)338F/518R(5′-CCT ACG GGA GGC AGC AG-3′,5′-ATT ACC GCG GCT GCT GG-3′)進(jìn)行擴(kuò)增(Muyzer et al.,1993)696-699。25 μL反應(yīng)體系為:5 ng?μL-1土壤微生物DNA模板1 μL,Green-2-Go qPCR Mastermix(上海生工)12.5 μL,去離子無菌水9.5 μL,10 μmol?μL-1的338F/518R引物各1 μL。qPCR反應(yīng)程序設(shè)定如下:95 ℃預(yù)變性3 min,95 ℃變性1 min,56 ℃退火1 min,72 ℃延伸1 min,39個(gè)循環(huán),最后72 ℃延伸5 min,擴(kuò)增效率90%~105%,R2>0.99。

    為獲取18S rRNA基因拷貝數(shù),將提取的土壤微生物DNA作為qPCR模板,采用對(duì)真菌具有特異性的引物對(duì)Fungi/NSI(5′-ATT CCC CGT TAC CCG TTG-3′,5′-GTA GTC ATA TGC TTG TCT C-3′)進(jìn)行擴(kuò)增(May et al.,2003)。25 μL反應(yīng)體系為:5 ng?μL-1土壤微生物DNA模板1 μL,Green-2-Go qPCR Mastermix(上海生工)12.5 μL,去離子無菌水9.5 μL,10 μmol?μL-1的Fungi/NSI引物各1 μL。qPCR反應(yīng)程序設(shè)定如下:95 ℃預(yù)變性15 min,95 ℃變性1 min,56 ℃退火1 min,72 ℃延伸1 min,39個(gè)循環(huán),最后72 ℃延伸5 min,擴(kuò)增效率90%~105%,R2>0.99。

    每克干土中的16S rRNA和18S rRNA基因拷貝數(shù)即為細(xì)菌和真菌豐度。

    1.2.3 土壤微生物群落分析

    利用細(xì)菌F968/R1401引物對(duì)(5′-CGC CCG GGG CGC GCC CCG GGC GGG GCG GGG GCA CGG GGG GAA CGC GAA GAA CCT TAC-3′,5′-CGG TGT GTA CAA GAC CC-3′)對(duì)土壤微生物DNA進(jìn)行PCR擴(kuò)增(靳振江等,2014b)4286。25 μL反應(yīng)體系為:5 ng?μL-1土壤微生物DNA模板1 μL,Taq PCR Master Mix(上海生工)12.5 μL,去離子無菌水9.5 μL,10 μmol?μL-1的F968/R1401引物各1 μL;擴(kuò)增條件為:94 ℃預(yù)變性3 min,94 ℃變性30 s,56 ℃退火30 s,72 ℃延伸1 min,30個(gè)循環(huán),最后72 ℃延伸10 min。

    在真菌正向引物Fungi 5′端添加GC夾子(5′-CCG CCG CGC GGC GGG CGG GGC GGG GGC ACG GGG-3′)后,將NSI作為反向引物對(duì)土壤DNA進(jìn)行PCR擴(kuò)增(靳振江等,2014b)4286。25 μL反應(yīng)體系為:5 ng?μL-1土壤微生物DNA模板1 μL,Taq PCR Master Mix(上海生工)12.5 μL,去離子無菌水9.5 μL,10 μmol?μL-1的正向引物和反向引物各1 μL;擴(kuò)增條件為:95 ℃預(yù)變性15 min,91 ℃變性1 min,57 ℃退火1 min,72 ℃延伸2 min,共35個(gè)循環(huán),然后在68 ℃完全延伸10 min。

    以上PCR擴(kuò)增產(chǎn)物用1.0%瓊脂糖凝膠進(jìn)行電泳檢測,檢測條帶清晰無拖尾表示擴(kuò)增結(jié)果良好,將PCR擴(kuò)增產(chǎn)物進(jìn)行變性梯度凝膠電泳分析。細(xì)菌和真菌PCR擴(kuò)增產(chǎn)物均采用8%的聚丙烯酰胺凝膠,變性劑濃度自下至上分別為55%~70%和15%~35%,并分別在60 ℃、100 V和60 ℃、90 V恒定電壓下電泳9 h,經(jīng)銀染后利用Gel-Doc-2000凝膠影像分析系統(tǒng)(美國Bio-Rad)拍照。

    1.3 數(shù)據(jù)處理與統(tǒng)計(jì)

    細(xì)菌和真菌豐度歸一化參照文獻(xiàn)中的方法(李強(qiáng)等,2014),將不同培養(yǎng)條件下的16S rRNA和18S rRNA基因拷貝數(shù)加和后分別進(jìn)行歸一化處理,計(jì)算公式為:

    式中,xi為細(xì)菌和真菌在17β-雌二醇不同添加量下的基因拷貝數(shù)。

    16S rRNA、18S rRNA的變性梯度凝膠拍照后,采用美國Bio-Rad公司的Imange lab分析軟件捕獲圖像進(jìn)行多樣性分析。多樣性指數(shù)的計(jì)算方法如下(靳振江等,2014b)4286:

    香農(nóng)-維納多樣性指數(shù):

    Margalef豐富度指數(shù):

    均勻度指數(shù):

    式中,S為OTUs總數(shù),即條帶數(shù);N為泳道中OTU的總亮度峰值;ni為第i個(gè)OTU在所有克隆中的比值。

    圖1 不同培養(yǎng)條件下16S rRNA、18S rRNA基因豐度Fig.1 Abundances of paddy soil bacteria and fungus under different treatment

    此外,利用SPSS軟件對(duì)巖溶稻田土壤微生物數(shù)據(jù)進(jìn)行雙因素方差分析。

    2 結(jié)果與討論

    2.1 不同培養(yǎng)條件對(duì)巖溶稻田土壤細(xì)菌和真菌豐度的影響

    在有氧培養(yǎng)條件下,巖溶稻田土壤細(xì)菌豐度為每克土壤1.180~2.240×1011copies,真菌豐度為每克土壤2.672~6.313×107copies;而在厭氧培養(yǎng)條件下,巖溶稻田土壤細(xì)菌豐度為每克土壤4.092~10.396×1011copies,真菌豐度為每克土壤0.204~0.395×107copies(圖1)。具有較高有機(jī)碳質(zhì)量分?jǐn)?shù)、pH值的巖溶水稻土對(duì)17β-雌二醇具有良好的吸附作用(Lai et al.,2000;靳振江等,2014a;李沛辰等,2015),為此吸附在土壤顆粒表面的17β-雌二醇在有氧條件下的降解速率高于厭氧條件(20 ℃時(shí),有氧條件下為0.120 g?h-1?L-1,厭氧條件下為0.057 g?h-1?L-1)(何芳等,2012),降解速率亦與微生物數(shù)量密切相關(guān)。此外,鑒于一定量的17β-雌二醇對(duì)厭氧微生物的生長具有促進(jìn)作用(趙穎等,2014)12-14,因此,受17β-雌二醇在不同供氧條件下降解速率的影響,在同一質(zhì)量分?jǐn)?shù)的17β-雌二醇暴露條件下,巖溶稻田土壤細(xì)菌豐度在有氧培養(yǎng)條件下普遍低于厭氧培養(yǎng)條件(圖1a)。而巖溶稻田土壤真菌豐度在有氧培養(yǎng)條件下普遍高于厭氧培養(yǎng)條件(圖1b),加之絕大多數(shù)真菌都是專性好氧菌,這一研究結(jié)果與熊鴻焰等(2009)認(rèn)為稻田水作造成真菌質(zhì)量分?jǐn)?shù)顯著降低的結(jié)果相一致。

    此外,通過圖1(a)還可看出,巖溶稻田土壤細(xì)菌豐度在17β-雌二醇的添加量為0.05 mg?kg-1和0.10 mg?kg-1時(shí)均高于對(duì)照組,并在厭氧培養(yǎng)條件下具有顯著性差異;而當(dāng)添加量為0.40(有氧)、0.20(厭氧)、0.40(厭氧)mg?kg-1時(shí)其細(xì)菌豐度則低于對(duì)照組。這是因?yàn)?,厭氧條件下,低量的17β-雌二醇能夠提高厭氧微生物產(chǎn)甲烷活性,而高量的17β-雌二醇則對(duì)厭氧微生物產(chǎn)甲烷活性產(chǎn)生抑制效應(yīng)(趙穎等,2014)12-14。加之17β-雌二醇作為雌激素對(duì)真菌菌絲體定向生長具有一定的促進(jìn)作用(郭堅(jiān)華等,1995),因此巖溶稻田土壤在面臨不同質(zhì)量分?jǐn)?shù)的17β-雌二醇暴露風(fēng)險(xiǎn)時(shí)其真菌豐度均高于對(duì)照組(圖1b)。由此說明,改善巖溶稻田土壤通氣狀況能夠降低17β-雌二醇對(duì)厭氧微生物的影響,并能緩解甲烷釋放速率。

    通常土壤細(xì)菌占土壤微生物總數(shù)的70%~90%(陳文新,1996),并對(duì)土壤微生物數(shù)量具有絕對(duì)控制作用。有氧和厭氧培養(yǎng)條件下巖溶稻田土壤細(xì)菌和真菌總拷貝數(shù)及歸一化結(jié)果顯示:17β-雌二醇的添加量為0.05 mg?kg-1時(shí),巖溶稻田土壤微生物豐度最高(表1),這可能是巖溶稻田土壤微生物對(duì)17β-雌二醇產(chǎn)生應(yīng)急反應(yīng)的一個(gè)重要應(yīng)答點(diǎn)。

    2.2 不同培養(yǎng)條件下巖溶稻田土壤細(xì)菌和真菌群落多樣性特征

    DGGE技術(shù)能夠直觀反映自然界微生物群落遺傳多樣性,對(duì)評(píng)價(jià)土壤環(huán)境變化具有重要的指導(dǎo)作用(Muyzer et al.,1993)696-699。根據(jù)電泳圖譜中每條條帶的信息,通過綜合分析各實(shí)驗(yàn)組中細(xì)菌和真菌的香農(nóng)-維納指數(shù)、豐富度指數(shù)和均勻度指數(shù),得知:厭氧培養(yǎng)條件下,細(xì)菌的香農(nóng)-維納指數(shù)和豐富度指數(shù)均高于對(duì)照組,均勻度指數(shù)則低于對(duì)照組,而真菌的香農(nóng)-維納指數(shù)、豐富度指數(shù)和均勻度指數(shù)均低于對(duì)照組;有氧培養(yǎng)條件下,細(xì)菌和真菌的香農(nóng)-維納指數(shù)除在17β-雌二醇的添加量為0.40 mg?kg-1時(shí)高于對(duì)照組外,其它條件下均低于對(duì)照組,并且其豐富度指數(shù)和均勻度指數(shù)均沒有明顯的變化規(guī)律(表2)。由此說明通氣狀況(有氧或厭氧)以及17β-雌二醇的添加量能影響巖溶稻田土壤生態(tài)系統(tǒng)的微生物多樣性。

    在厭氧培養(yǎng)條件下,針對(duì)不同量的17β-雌二醇添加實(shí)驗(yàn),巖溶稻田土壤細(xì)菌豐度與其香農(nóng)-維納指數(shù)和豐富度指數(shù)總體上呈現(xiàn)一致的變化趨勢,但與其均勻度指數(shù)呈現(xiàn)相反的變化趨勢;巖溶稻田土壤真菌豐度與其3種多樣性指數(shù)總體上呈現(xiàn)相反的變化趨勢(圖1,表3)。由于均勻度與物種數(shù)目無關(guān),只與個(gè)體數(shù)目在各個(gè)物種中分布的均勻程度有關(guān),故群落中的常見種與稀少種的差距逐漸變小,整個(gè)群落向著物種均勻化方向發(fā)展(賀紀(jì)正等,2013)414-416。盡管17β-雌二醇能夠促進(jìn)少數(shù)微生物(細(xì)菌,特別是產(chǎn)甲烷菌)種屬及其數(shù)量的增加,但厭氧培養(yǎng)條件下巖溶稻田土壤細(xì)菌和真菌均勻度的降低則表明17β-雌二醇對(duì)巖溶稻田土壤微生物產(chǎn)生了抑制作用,這種抑制作用不但降低了巖溶稻田土壤微生物生態(tài)系統(tǒng)的穩(wěn)定性(賀紀(jì)正等,2013)415,而且將加劇巖溶稻田在淹水期間的甲烷排放量(秦曉波等,2012),進(jìn)而影響全球變化。

    表1 不同培養(yǎng)條件下巖溶稻田土壤微生物豐度Table 1 Abundances of karst paddy soil microbes under different treatment

    表2 不同培養(yǎng)條件下基于巖溶稻田土壤16S rRNA、18S rRNA DGGE圖譜的微生物多樣性指數(shù)Table 2 The karst paddy soil microbial diversity based on the DGGE profile of 16S rRNA and 18S rRNA under different treatment

    表3 土壤通氣狀況、17β-雌二醇添加量及其交互作用與16S rRNA、18S rRNA基因豐度雙因素方差分析Table 3 Two-way analysis of variance relating to experimental condition and the abundances of kart paddy soil bacteria and fungus

    在有氧條件下,17β-雌二醇通常被微生物脫氫氧化為雌激素酮(何芳等,2008)1962-1964,而雌激素酮具有較長的半衰期并更容易在微生物體內(nèi)積累(何芳等,20081963-1965;王代懿等,2014)。因此,在有氧條件下,17β-雌二醇添加量為0.05、0.10、0.20 mg?kg-1時(shí)的巖溶稻田土壤細(xì)菌和真菌具有較低的香農(nóng)-維納指數(shù),而細(xì)菌和真菌豐度與其香農(nóng)-維納指數(shù)呈現(xiàn)相反的變化趨勢。由于土壤中高添加量(0.40 mg?kg-1)的17β-雌二醇在有氧條件下降解不完全,因此巖溶稻田土壤細(xì)菌和真菌在17β-雌二醇和雌激素酮的交互作用下具有較高的香農(nóng)-維納指數(shù)。

    土壤中絕大多數(shù)真菌都是專性好氧菌(靳振江等,2014b)4289,受氧氣體積分?jǐn)?shù)的限制以及17β-雌二醇的抑制作用,巖溶稻田土壤真菌的3種多樣性指標(biāo)在有氧培養(yǎng)條件下均高于厭氧培養(yǎng)條件(表2)。此外,由于目前已知能夠降解17β-雌二醇的微生物大部分為細(xì)菌(Li et al.,2012;方會(huì)等,2012),加之17β-雌二醇在厭氧條件下能夠提高厭氧微生物的活性,因此,除香農(nóng)-維納指數(shù)外,巖溶稻田土壤細(xì)菌的其它兩種多樣性指標(biāo)總體上表現(xiàn)出有氧條件>厭氧條件(表2)。

    綜上所述,巖溶稻田土壤微生物的均勻度在維持巖溶稻田土壤生態(tài)系統(tǒng)微生物多樣性-穩(wěn)定性關(guān)系上具有重要的指示作用。因此,本研究結(jié)果表明,在17β-雌二醇暴露環(huán)境中,均勻度指數(shù)能夠作為水、旱交替條件下巖溶稻田土壤微生物群落特性變化的評(píng)價(jià)指標(biāo)。

    2.3 不同培養(yǎng)條件與巖溶稻田土壤微生物群落的雙因素方差分析

    采用F檢驗(yàn)對(duì)本實(shí)驗(yàn)所涉及的土壤通氣狀況(有氧、厭氧)、17β-雌二醇添加量、細(xì)菌和真菌豐度進(jìn)行雙因素方差分析,結(jié)果如表3所示。巖溶稻田土壤通氣狀況(有氧/厭氧)和17β-雌二醇添加量對(duì)巖溶稻田土壤細(xì)菌豐度有極顯著影響,而其交互作用對(duì)巖溶稻田土壤細(xì)菌豐度有顯著影響;土壤通氣狀況(有氧/厭氧)對(duì)巖溶稻田土壤真菌豐度也產(chǎn)生極顯著影響,但17β-雌二醇添加量及其與土壤通氣狀況的交互作用對(duì)巖溶稻田土壤真菌豐度影響不顯著。

    土壤通氣狀況(有氧、厭氧)能夠?qū)r溶稻田土壤細(xì)菌和真菌豐度產(chǎn)生極顯著影響,可能與17β-雌二醇的生物降解特性有關(guān)。在有氧條件下,17β-雌二醇的降解速率高于厭氧條件(何芳等,2008)1963-1965并符合一級(jí)反應(yīng)動(dòng)力學(xué)(反應(yīng)速率常數(shù)0.23~4.79)(Li et al.,2005)。此外,17β-雌二醇在有氧條件下能夠?yàn)橥寥乐械腗n(Ⅳ)或Fe(Ⅲ)提供電子,自身被氧化成雌激素酮(Oscarson et al., 1981),最終被微生物降解。而Lee et al.(2002)在有氧和厭氧條件下利用活性污泥處理17β-雌二醇,發(fā)現(xiàn)水體中17β-雌二醇質(zhì)量濃度為200 μg?L-1時(shí)幾乎能完全被降解掉。因此,本實(shí)驗(yàn)結(jié)果說明改善巖溶區(qū)稻田土壤通氣狀況,定期對(duì)稻田進(jìn)行翻耕、曬田或水-旱輪作能夠降低17β-雌二醇對(duì)土壤微生物及生態(tài)環(huán)境帶來的危害。

    3 結(jié)論

    (1)17β-雌二醇對(duì)巖溶稻田土壤微生物的生長具有一定的表觀刺激效應(yīng),這種表觀刺激效應(yīng)在厭氧條件下極顯著。

    (2)17β-雌二醇在促進(jìn)巖溶稻田厭氧微生物活性的同時(shí),降低了巖溶稻田土壤微生物系統(tǒng)的抵抗力和恢復(fù)力,進(jìn)而降低了巖溶稻田土壤生態(tài)系統(tǒng)微生物的多樣性及其穩(wěn)定性。

    (3)定期對(duì)巖溶區(qū)的稻田進(jìn)行翻耕、曬田或水-旱輪作能夠促進(jìn)17β-雌二醇的降解性能,進(jìn)而緩解17β-雌二醇對(duì)巖溶稻田土壤微生物造成的生態(tài)風(fēng)險(xiǎn)。

    BECHI N, SORDA G, SPAGNOLETTI A, et al.2013 Toxicity assessment on trophoblast cells for some environment polluting chemicals and 17β-estradiol [J].Toxicology in Vitro, 27(3): 995-1000.

    GONG J, RAN Y, CHEN D Y, et al.2011.Occurrence of endocrine-disrupting chemicals in riverine sediments from the Pearl River Delta, China [J].Marine Pollution Bulletin, 63(5): 556-563.

    KOLPIN D W, FURLONG E T, MEYER M T, et al.2002.Pharmaceuticals, hormones, and other organic wastewater contaminants in US streams, 1999—2000: a national reconnaissance [J].Environmental Science & Technology, 36(6): 1202-1211.

    K?RNER W, HANF V, SCHULLER W, et al.1999.Development of a sensitive E-screen assay for quantitative analysis of estrogenic activity in municipal sewage plant effluents [J].Science of the Total Environment, 225(1): 33-48.

    KUCH H M, BALLSCHMITER K.2001.Determination of endocrinedisrupting phenolic compounds and estrogens in surface and drinking water by HRGC-(NCI)-MS in the picogram per liter range [J].Environmental Science & Technology, 35(15): 3201-3206.

    LAI K M, JOHNSON K L, SCRIMSHAW M D, et al.2000.Binding of waterborne steroid estrogens to solid phases in river and estuarine systems [J].Environmental Science & Technology, 34(18): 3890-3894.

    LEE H B, LIU D.2002.Degradation of 17 β-estradiol and its metabolites by sewage bacteria [J].Water, Air, and Soil Pollution, 134(1-4): 351-366.

    LI F, YUASA A, OBARA A, et al.2005.Aerobic batch degradation of 17-β estradiol (E2) by activated sludge: Effects of spiking E2 concentrations, MLVSS and temperatures [J].Water Research, 39(10): 2065-2075.

    LI Z, NANDAKUMAR R, MADAYIPUTHIYA N, et al.2012.Proteomic analysis of 17β-estradiol degradation by Stenotrophomonas maltophilia [J].Environmental Science & Technology, 46(11): 5947-5955.

    MAY L A, SMILEY B, SCHMIDT M G.2003.Comparative denaturing gradient gel electrophoresis analysis of gungal communities associated with whole plant com silage [J].Canadian Journal of Microbiology Ecology, 43(1): 1-11.

    MUYZER G, DE WAAL E C, UITTERLINDEN A G.1993.Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA [J].Applied and Environmental Microbiology, 59(3): 695-700.

    OSCARSON D W, HUANG P M, DEFOSSE C, et al.1981.Oxidative power of Mn (IV) and Fe (III) oxides with respect to As (III) in terrestrial and aquatic environments [J].Nature, 291(5810): 50-51.

    陳家長, 王澤镕, 瞿建宏, 等.2012.17β-雌二醇與1-萘酚對(duì)雄性羅非魚(GIFT Oreochromis niloticus) 雌激素效應(yīng)的比較[J].生態(tài)環(huán)境學(xué)報(bào), 21(4): 754-759.

    陳文新.1996.土壤和環(huán)境微生物學(xué)[M].北京: 北京農(nóng)業(yè)大學(xué)出版社.

    陳曉鋒, 常志州, 黃紅英, 等.2009.養(yǎng)殖污水回用對(duì)水稻產(chǎn)量及土壤肥力的影響研究[J].中國土壤與肥料, (2): 39-42, 70.

    方會(huì), 陳梅雪, 秦德韜, 等.2012.集約化養(yǎng)豬廢水SBR中17β-雌二醇高效降解菌的分離鑒定及其降解特性[J].環(huán)境工程學(xué)報(bào), 6(10): 3361-3366.

    龔劍, 冉勇, 陳迪云, 等.2011.珠江三角洲兩條主要河流沉積物中的典型內(nèi)分泌干擾物污染狀況[J].生態(tài)環(huán)境學(xué)報(bào), 21(4): 754-759.

    郭堅(jiān)華, 曹以勤.1995.綿霉屬 (Achlya) 真菌性激素的研究進(jìn)展[J].生物學(xué)雜志, (3): 11-13.

    何芳, 李富生, AKIRA Y, 等.2012.湖泊沉積物對(duì)17 β-雌二醇的降解效能[J].土木建筑與環(huán)境工程, 34(4): 125-130.

    何芳, 李富生, AKIRA Y.2008.湖泊底泥中17 β-雌二醇的生物代謝機(jī)制研究[J].環(huán)境科學(xué), 29(7): 1961-1966.

    賀紀(jì)正, 李晶, 鄭袁明.2013.土壤生態(tài)系統(tǒng)微生物多樣性-穩(wěn)定性關(guān)系的思考[J].生物多樣性, 21(4): 412-421.

    靳振江, 程亞平, 李強(qiáng), 等.2014a.會(huì)仙喀斯特溶洞濕地、稻田和旱田土壤有機(jī)碳含量及其與養(yǎng)分的關(guān)系[J].濕地科學(xué), 12(4): 485-490.

    靳振江, 湯華峰, 李敏, 等.2014b.典型巖溶土壤微生物豐度與多樣性及其對(duì)碳循環(huán)的指示意義[J].環(huán)境科學(xué), 35(11): 4284-4290.

    李沛辰, 毋偉, 張豐松, 等.2015.秸稈生物碳的結(jié)構(gòu)特征及其對(duì)17β-雌二醇的吸附性能[J].環(huán)境科學(xué)研究, 28(8): 1260-1266.

    李強(qiáng), 胡清菁, 張超蘭, 等.2014.基于土壤酶總體活性評(píng)價(jià)鉛鋅尾礦砂坍塌區(qū)土壤重金屬污染[J].生態(tài)環(huán)境學(xué)報(bào), 23(11): 1839-1844.

    聯(lián)合國環(huán)境規(guī)劃署.2001.關(guān)于持久性有機(jī)污染物的斯德哥爾摩公約[Z].斯德哥爾摩: 聯(lián)合國環(huán)境規(guī)劃署.

    秦曉波, 李玉娥, 石生偉, 等.2012.稻田溫室氣體排放與土壤微生物菌群的多元回歸分析[J].生態(tài)學(xué)報(bào), 32(6): 1811-1819.

    邵曉玲, 馬軍.2008.松花江水中13種內(nèi)分泌干擾物的初步調(diào)查[J].環(huán)境科學(xué)學(xué)報(bào), 28(9): 1910-1915.

    王代懿, 梁漢東, 張豐松, 等.2014.水分對(duì)典型天然性激素在我國東北黑土中降解的影響[J].安全與環(huán)境學(xué)報(bào), 14(4): 232-237.

    吳偉恒, 阮愛東, 戴韻秋, 等.2014.雌激素污染對(duì)水體微生物反硝化作用的影響[J].環(huán)境科學(xué)與技術(shù), 37(12): 1-5, 15.

    熊鴻焰, 李廷軒, 余海英, 等.2009.水旱輪作條件下免耕土壤微生物特性研究[J].植物營養(yǎng)與肥料學(xué)報(bào), 15(1): 145-150.

    趙穎, 阮愛東, 劉忱瀟, 等.2014.雌激素污染對(duì)水體微生物產(chǎn)甲烷功能的影響[J].環(huán)境科學(xué)與技術(shù), 37(12): 11-15.

    Responses of Karst Paddy Soil Microbial Community Structure to 17β-estradiol in the Simulated Experiments

    HUANG Yadan1, JIN Zhenjiang2,4, LI Qiang2,3, LIU Chang2
    1.Graduate School of Guilin Medical University, Guilin 541004, China; 2.Institute of Karst Geology, Chinese Academy of Geological Sciences, Guilin 541004, China; 3.International Research Center on Karst under the Auspices of UNESCO, Guilin 541004, China; 4.Environmental Science and Engineering College, Guilin University of Technology, Guilin 541004, China

    To understand the risk of 17β-estradiol on the microbial diversity-stability relationship in the karst paddy soil, the exogenous 17β-estradiol with the content of 0, 0.05, 0.10, 0.20, 0.40 mg?kg-1was added to the karst paddy soil.Moreover, the soil was cultured in the aerobic or anaerobic condition.At the end of the culture, Quantitative Real-time Polymerase Chain Reaction and Polymerase Chain Reaction-Denaturing Gradient Gel Electrophoresis methods were applied to analyze the karst paddy soil microbial abundance and diversity.17β-estradiol (0.05 mg?kg-1and 0.10 mg?kg-1) has the apparent stimulatory effects on the abundances of karst paddy soil bacteria and fungus under different treatment, which is quite significant in the anaerobic condition.Moreover, according to the normalized karst paddy soil bacteria and fungus abundances, it was inferred that one of the critical and emergency response point from 17β-estradiol to karst paddy soil microbes is 0.05 mg?kg-1.In the anaerobic condition, the species evenness of karst paddy soil bacteria and fungus is lower than the control group and also has the opposite change with their abundance.In the aerobic condition, 17β-estradiol can be degraded.The karst paddy soil microbial diversity has no clear variation in the aerobic condition.Though some anaerobic microorganisms were increased, the karst paddy soil microbial abundance and diversity were depressed.Moreover, the two-way analysis of variance relating to experimental condition and the abundances of karst paddy soil bacteria proved that the soil aeration condition is the controlling factor for karst paddy soil microbial community.From the results, it can be inferred that tillage, sunning the fields or waterlog-drought rotation can improve the degradation rate of 17 β-estradiol in the karst paddy soil, which decreased the anaerobic microbe activity and their methanogenic function.Therefore, the eco-environmental risk for karst paddy soil microbe from 17β-estradiol will be reduced.

    17β-estradiol; soil aeration condition; microbial abundance; microbial diversity

    10.16258/j.cnki.1674-5906.2016.10.018

    X172; P593

    A

    1674-5906(2016)10-1721-06

    黃雅丹, 靳振江, 李強(qiáng), 劉暢.2016.17β-雌二醇對(duì)巖溶稻田土壤微生物群落的影響[J].生態(tài)環(huán)境學(xué)報(bào), 25(10): 1721-1726.

    HUANG Yadan, JIN Zhenjiang, LI Qiang, LIU Chang.2016.Responses of Karst paddy soil microbial community structure to 17β-estradiol in the simulated experiments [J].Ecology and Environmental Sciences, 25(10): 1721-1726.

    國家重點(diǎn)研發(fā)計(jì)劃課題(2016YFC0502501);國家自然科學(xué)基金項(xiàng)目(41641026;41361054);廣西自然科學(xué)基金項(xiàng)目(2015GXNSFGA139010;2014GXNSFCA118012);廣西科學(xué)研究與技術(shù)開發(fā)計(jì)劃項(xiàng)目(桂科合14123001-13);中國地質(zhì)科學(xué)院項(xiàng)目(YYWF201505);桂林市科學(xué)研究與技術(shù)開發(fā)計(jì)劃項(xiàng)目(20140122-1)

    黃雅丹(1980年生),女,助理研究員,碩士,主要從事環(huán)境生物學(xué)研究。*通信作者。李強(qiáng),E-mail: glqiangli@163.com

    2016-08-12

    女性被躁到高潮视频| 欧美日韩乱码在线| 91九色精品人成在线观看| 国产亚洲欧美98| 91大片在线观看| 国产黄a三级三级三级人| 欧美黑人精品巨大| 日韩精品青青久久久久久| 夜夜看夜夜爽夜夜摸 | 老司机亚洲免费影院| 在线国产一区二区在线| 看片在线看免费视频| 丁香六月欧美| 久久国产亚洲av麻豆专区| 怎么达到女性高潮| 国产高清国产精品国产三级| 欧美丝袜亚洲另类 | 亚洲精品国产一区二区精华液| 黄色片一级片一级黄色片| 免费人成视频x8x8入口观看| 露出奶头的视频| av有码第一页| 岛国视频午夜一区免费看| 男女下面进入的视频免费午夜 | 高清黄色对白视频在线免费看| 亚洲av电影在线进入| 久久久水蜜桃国产精品网| 免费在线观看完整版高清| 精品一品国产午夜福利视频| 天天影视国产精品| 免费看a级黄色片| 韩国精品一区二区三区| 一边摸一边做爽爽视频免费| 成年版毛片免费区| 中文字幕精品免费在线观看视频| 国产精品久久久久成人av| 窝窝影院91人妻| 在线观看免费高清a一片| 搡老乐熟女国产| 国产成人一区二区三区免费视频网站| 首页视频小说图片口味搜索| 桃色一区二区三区在线观看| 成人国语在线视频| www国产在线视频色| 99热国产这里只有精品6| 欧美成狂野欧美在线观看| 亚洲五月婷婷丁香| 午夜免费激情av| 久久99一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 久久精品国产亚洲av香蕉五月| 男人舔女人下体高潮全视频| 国产一区在线观看成人免费| 精品国产一区二区三区四区第35| 免费少妇av软件| 国产xxxxx性猛交| 成人三级黄色视频| 成人三级黄色视频| 久久久国产欧美日韩av| 欧美日本亚洲视频在线播放| 亚洲美女黄片视频| 精品国产一区二区三区四区第35| 岛国视频午夜一区免费看| www.999成人在线观看| 亚洲第一青青草原| 亚洲av成人av| 国产av精品麻豆| 18禁美女被吸乳视频| 一级片免费观看大全| 中文字幕人妻丝袜一区二区| 身体一侧抽搐| 悠悠久久av| 欧美一区二区精品小视频在线| 悠悠久久av| 亚洲美女黄片视频| 99re在线观看精品视频| 国产高清激情床上av| 在线观看一区二区三区激情| 999久久久国产精品视频| 波多野结衣高清无吗| 亚洲成人精品中文字幕电影 | 成年人免费黄色播放视频| 亚洲熟妇熟女久久| 国产亚洲精品久久久久5区| 99精国产麻豆久久婷婷| 国产成人精品无人区| 亚洲自拍偷在线| 麻豆久久精品国产亚洲av | 亚洲九九香蕉| av有码第一页| 免费女性裸体啪啪无遮挡网站| 亚洲黑人精品在线| 麻豆国产av国片精品| 亚洲成a人片在线一区二区| av网站在线播放免费| 国产单亲对白刺激| 久久人妻av系列| 日韩大尺度精品在线看网址 | 精品国产国语对白av| 亚洲国产欧美网| 777久久人妻少妇嫩草av网站| 老司机亚洲免费影院| 国产亚洲av高清不卡| 99久久久亚洲精品蜜臀av| 新久久久久国产一级毛片| 国产高清激情床上av| 欧美黄色片欧美黄色片| 成人免费观看视频高清| 日本免费a在线| 国产精品 国内视频| 一级作爱视频免费观看| 亚洲欧美精品综合久久99| www.熟女人妻精品国产| 精品第一国产精品| xxxhd国产人妻xxx| 国产亚洲精品综合一区在线观看 | 欧美国产精品va在线观看不卡| 亚洲午夜精品一区,二区,三区| 美女国产高潮福利片在线看| 国产乱人伦免费视频| 妹子高潮喷水视频| 韩国精品一区二区三区| 丁香欧美五月| 日韩欧美一区二区三区在线观看| 亚洲精品美女久久av网站| 五月开心婷婷网| 黄网站色视频无遮挡免费观看| 变态另类成人亚洲欧美熟女 | 两个人看的免费小视频| 亚洲情色 制服丝袜| 美国免费a级毛片| 国产精品1区2区在线观看.| 久久精品91蜜桃| 午夜免费成人在线视频| 两个人免费观看高清视频| 亚洲av美国av| 啦啦啦 在线观看视频| 午夜亚洲福利在线播放| 午夜精品在线福利| 亚洲成a人片在线一区二区| 免费av中文字幕在线| 91精品国产国语对白视频| 久久精品亚洲精品国产色婷小说| 欧美日韩av久久| 亚洲 欧美 日韩 在线 免费| 欧美日本亚洲视频在线播放| 亚洲视频免费观看视频| 国产91精品成人一区二区三区| 天堂√8在线中文| 一二三四社区在线视频社区8| 久久久精品欧美日韩精品| 欧美激情高清一区二区三区| 69av精品久久久久久| 欧美色视频一区免费| 国产av一区二区精品久久| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品久久成人aⅴ小说| 88av欧美| 久久久久久久午夜电影 | 精品第一国产精品| 看片在线看免费视频| 国产精品国产高清国产av| 久久久久久免费高清国产稀缺| 亚洲全国av大片| 国产蜜桃级精品一区二区三区| 日韩三级视频一区二区三区| 久久久久久久精品吃奶| 国产欧美日韩一区二区精品| 国产成人欧美在线观看| 法律面前人人平等表现在哪些方面| 国产在线精品亚洲第一网站| netflix在线观看网站| 亚洲精品成人av观看孕妇| 乱人伦中国视频| 两个人免费观看高清视频| 日韩三级视频一区二区三区| а√天堂www在线а√下载| 中文字幕人妻熟女乱码| 亚洲男人的天堂狠狠| 国产成人精品久久二区二区免费| 国产精品98久久久久久宅男小说| 国产欧美日韩一区二区三区在线| 欧美乱码精品一区二区三区| 真人做人爱边吃奶动态| 午夜激情av网站| 丝袜在线中文字幕| 波多野结衣av一区二区av| 成年人免费黄色播放视频| а√天堂www在线а√下载| 国产欧美日韩一区二区精品| 国产成人欧美在线观看| 丁香六月欧美| 国产精品亚洲一级av第二区| 亚洲专区中文字幕在线| 久久精品国产亚洲av香蕉五月| 女生性感内裤真人,穿戴方法视频| 亚洲国产精品一区二区三区在线| 免费人成视频x8x8入口观看| 视频区欧美日本亚洲| www日本在线高清视频| 热re99久久国产66热| 国产激情久久老熟女| 国产精品秋霞免费鲁丝片| 成人18禁在线播放| 精品久久久精品久久久| 又黄又爽又免费观看的视频| 日本免费一区二区三区高清不卡 | 色在线成人网| 亚洲成人国产一区在线观看| 99国产综合亚洲精品| 麻豆国产av国片精品| 法律面前人人平等表现在哪些方面| 9191精品国产免费久久| 热re99久久精品国产66热6| 国产精品香港三级国产av潘金莲| 国产免费现黄频在线看| av超薄肉色丝袜交足视频| 欧美色视频一区免费| 亚洲avbb在线观看| 热re99久久国产66热| 成人特级黄色片久久久久久久| 亚洲国产毛片av蜜桃av| 美女高潮喷水抽搐中文字幕| 亚洲一码二码三码区别大吗| 久9热在线精品视频| 欧美久久黑人一区二区| 国产免费男女视频| 国产精品98久久久久久宅男小说| 国产高清视频在线播放一区| 亚洲第一青青草原| 亚洲黑人精品在线| av片东京热男人的天堂| 最新美女视频免费是黄的| 亚洲 欧美一区二区三区| 激情在线观看视频在线高清| 一进一出抽搐gif免费好疼 | 国产精品美女特级片免费视频播放器 | 操出白浆在线播放| 亚洲国产中文字幕在线视频| 亚洲自偷自拍图片 自拍| 免费在线观看影片大全网站| 久热这里只有精品99| 久久久久久大精品| 18禁裸乳无遮挡免费网站照片 | 亚洲国产欧美网| 99国产精品一区二区三区| 黄色 视频免费看| 国产免费av片在线观看野外av| 女生性感内裤真人,穿戴方法视频| 欧美在线一区亚洲| 精品国产一区二区久久| 国产99久久九九免费精品| 亚洲一区二区三区色噜噜 | 国产精品 国内视频| 久久草成人影院| 丝袜美足系列| 天堂俺去俺来也www色官网| 国产高清激情床上av| 亚洲欧美一区二区三区黑人| 日本vs欧美在线观看视频| 国产激情欧美一区二区| 黄网站色视频无遮挡免费观看| 两个人免费观看高清视频| 亚洲黑人精品在线| 大码成人一级视频| 亚洲av片天天在线观看| 国产成人欧美| 1024香蕉在线观看| 丝袜人妻中文字幕| 亚洲精品国产一区二区精华液| 色婷婷久久久亚洲欧美| 日本欧美视频一区| 欧美日韩亚洲综合一区二区三区_| 欧美日韩一级在线毛片| 啦啦啦在线免费观看视频4| 麻豆久久精品国产亚洲av | 两个人看的免费小视频| 激情在线观看视频在线高清| 国产又色又爽无遮挡免费看| 国产成人精品久久二区二区91| 两性午夜刺激爽爽歪歪视频在线观看 | 男女之事视频高清在线观看| 欧美成狂野欧美在线观看| 在线播放国产精品三级| 国产视频一区二区在线看| 国产不卡一卡二| 国产亚洲精品久久久久5区| 亚洲免费av在线视频| 精品卡一卡二卡四卡免费| 国产深夜福利视频在线观看| 亚洲熟妇中文字幕五十中出 | 亚洲欧美日韩另类电影网站| 久久人人精品亚洲av| 大陆偷拍与自拍| 国产亚洲欧美98| 老司机深夜福利视频在线观看| 久久婷婷成人综合色麻豆| 成年人免费黄色播放视频| 日本一区二区免费在线视频| 高潮久久久久久久久久久不卡| www日本在线高清视频| 亚洲国产精品sss在线观看 | 一区二区三区精品91| 欧美日韩亚洲国产一区二区在线观看| 国产伦一二天堂av在线观看| 91字幕亚洲| 国产精品九九99| 青草久久国产| 国产成人啪精品午夜网站| 一区二区三区国产精品乱码| 99久久久亚洲精品蜜臀av| 免费在线观看完整版高清| 国产av又大| 成人18禁高潮啪啪吃奶动态图| 9热在线视频观看99| 精品人妻在线不人妻| 国产成人精品无人区| 久久久国产一区二区| 手机成人av网站| 亚洲国产精品sss在线观看 | 日本wwww免费看| 91九色精品人成在线观看| 最好的美女福利视频网| 欧美激情高清一区二区三区| 欧美中文日本在线观看视频| 757午夜福利合集在线观看| 日韩 欧美 亚洲 中文字幕| 亚洲aⅴ乱码一区二区在线播放 | 亚洲人成电影免费在线| 日日摸夜夜添夜夜添小说| 国产精品影院久久| 亚洲黑人精品在线| xxx96com| 亚洲一区中文字幕在线| 涩涩av久久男人的天堂| 亚洲九九香蕉| 精品久久蜜臀av无| 色老头精品视频在线观看| 日韩欧美三级三区| 后天国语完整版免费观看| 50天的宝宝边吃奶边哭怎么回事| 午夜精品国产一区二区电影| 亚洲精品美女久久av网站| 国产熟女午夜一区二区三区| 亚洲中文日韩欧美视频| 免费看a级黄色片| 国产成年人精品一区二区 | 亚洲情色 制服丝袜| 每晚都被弄得嗷嗷叫到高潮| a级片在线免费高清观看视频| 中亚洲国语对白在线视频| 搡老乐熟女国产| 纯流量卡能插随身wifi吗| 极品教师在线免费播放| 中出人妻视频一区二区| 97人妻天天添夜夜摸| 久久精品国产亚洲av香蕉五月| 成年人黄色毛片网站| 精品卡一卡二卡四卡免费| 国产精品久久久av美女十八| 中国美女看黄片| 日本三级黄在线观看| 精品国产超薄肉色丝袜足j| 久久狼人影院| 高清黄色对白视频在线免费看| 一级片'在线观看视频| 欧美不卡视频在线免费观看 | 动漫黄色视频在线观看| 黄片小视频在线播放| 人人妻人人澡人人看| av电影中文网址| 日韩 欧美 亚洲 中文字幕| 日韩人妻精品一区2区三区| xxxhd国产人妻xxx| 亚洲 国产 在线| 亚洲aⅴ乱码一区二区在线播放 | 一级片'在线观看视频| 亚洲欧美激情在线| 国产三级黄色录像| x7x7x7水蜜桃| 在线十欧美十亚洲十日本专区| av在线播放免费不卡| 欧美乱妇无乱码| 女人爽到高潮嗷嗷叫在线视频| 精品一品国产午夜福利视频| tocl精华| 欧美乱色亚洲激情| 在线观看免费视频网站a站| 久久久国产成人免费| 老司机深夜福利视频在线观看| 国产免费av片在线观看野外av| 别揉我奶头~嗯~啊~动态视频| 女人爽到高潮嗷嗷叫在线视频| av欧美777| 国产熟女午夜一区二区三区| 不卡av一区二区三区| 欧美黑人欧美精品刺激| 国产欧美日韩一区二区精品| 热re99久久国产66热| 天天躁狠狠躁夜夜躁狠狠躁| 男人的好看免费观看在线视频 | 男人的好看免费观看在线视频 | 日韩大尺度精品在线看网址 | www.自偷自拍.com| 日本免费a在线| 国产激情久久老熟女| 亚洲精品中文字幕在线视频| 国产精品免费一区二区三区在线| 18禁黄网站禁片午夜丰满| 欧美激情久久久久久爽电影 | 老汉色av国产亚洲站长工具| 欧美成狂野欧美在线观看| 国产亚洲精品综合一区在线观看 | 免费高清视频大片| 久99久视频精品免费| 老熟妇乱子伦视频在线观看| 亚洲久久久国产精品| 国产区一区二久久| 国产精品一区二区精品视频观看| 国产精品一区二区三区四区久久 | 我的亚洲天堂| 在线观看www视频免费| 久久久久久久午夜电影 | 老汉色∧v一级毛片| 亚洲av电影在线进入| 色综合欧美亚洲国产小说| 一区二区三区激情视频| 国产精品自产拍在线观看55亚洲| 国产av一区在线观看免费| 不卡一级毛片| 一个人观看的视频www高清免费观看 | 伦理电影免费视频| 免费在线观看影片大全网站| 俄罗斯特黄特色一大片| 9色porny在线观看| 99久久精品国产亚洲精品| 久久久久国产精品人妻aⅴ院| 男女床上黄色一级片免费看| 女人被狂操c到高潮| 成人特级黄色片久久久久久久| 中文字幕另类日韩欧美亚洲嫩草| 9色porny在线观看| 亚洲国产精品999在线| 91国产中文字幕| netflix在线观看网站| 亚洲视频免费观看视频| 首页视频小说图片口味搜索| 亚洲成人免费电影在线观看| 大陆偷拍与自拍| 精品一区二区三区av网在线观看| 国产真人三级小视频在线观看| 别揉我奶头~嗯~啊~动态视频| 黑人巨大精品欧美一区二区mp4| 久久人妻福利社区极品人妻图片| av天堂久久9| 国产xxxxx性猛交| 久99久视频精品免费| 国产精品亚洲一级av第二区| 99国产综合亚洲精品| 制服人妻中文乱码| 亚洲国产毛片av蜜桃av| 精品熟女少妇八av免费久了| 不卡av一区二区三区| 国产精品爽爽va在线观看网站 | 999精品在线视频| 亚洲国产欧美一区二区综合| 欧美日韩国产mv在线观看视频| 在线天堂中文资源库| 人妻丰满熟妇av一区二区三区| 国产精品国产av在线观看| 首页视频小说图片口味搜索| 天天躁夜夜躁狠狠躁躁| 久久性视频一级片| 欧美日韩瑟瑟在线播放| 欧美激情极品国产一区二区三区| 夫妻午夜视频| 日本免费一区二区三区高清不卡 | 免费人成视频x8x8入口观看| 国产成人精品久久二区二区免费| 国产一区在线观看成人免费| 黄色片一级片一级黄色片| 精品国产超薄肉色丝袜足j| 国产精品久久久久久人妻精品电影| 欧美乱码精品一区二区三区| 18禁国产床啪视频网站| 婷婷丁香在线五月| av天堂在线播放| 一级作爱视频免费观看| 又紧又爽又黄一区二区| 久久精品国产亚洲av高清一级| 看黄色毛片网站| 日韩三级视频一区二区三区| 丰满迷人的少妇在线观看| 啦啦啦免费观看视频1| 日韩一卡2卡3卡4卡2021年| 99久久久亚洲精品蜜臀av| 桃红色精品国产亚洲av| 手机成人av网站| 男女之事视频高清在线观看| 亚洲av片天天在线观看| 91精品三级在线观看| www.熟女人妻精品国产| 国产精品综合久久久久久久免费 | 午夜福利一区二区在线看| 欧美一级毛片孕妇| 国产亚洲精品第一综合不卡| 后天国语完整版免费观看| 一级黄色大片毛片| 欧美成人性av电影在线观看| 亚洲 欧美一区二区三区| 宅男免费午夜| av网站在线播放免费| 久久久久久人人人人人| 久久天躁狠狠躁夜夜2o2o| 激情在线观看视频在线高清| 国产伦人伦偷精品视频| 91精品国产国语对白视频| 动漫黄色视频在线观看| a在线观看视频网站| 亚洲欧美激情在线| 又大又爽又粗| 国产成+人综合+亚洲专区| 黑丝袜美女国产一区| 日韩欧美在线二视频| 国产在线观看jvid| 久久久久久亚洲精品国产蜜桃av| 欧美+亚洲+日韩+国产| 亚洲一卡2卡3卡4卡5卡精品中文| 国产一区二区三区视频了| 99riav亚洲国产免费| 国产精品亚洲一级av第二区| 亚洲精华国产精华精| 国产精品1区2区在线观看.| 少妇裸体淫交视频免费看高清 | 国产97色在线日韩免费| 99国产极品粉嫩在线观看| 操出白浆在线播放| 国产亚洲av高清不卡| 一边摸一边抽搐一进一小说| 欧美日韩精品网址| 亚洲aⅴ乱码一区二区在线播放 | 999精品在线视频| 国产色视频综合| 欧美亚洲日本最大视频资源| 国产精品成人在线| 亚洲一区二区三区色噜噜 | 亚洲欧美日韩高清在线视频| 88av欧美| 国产99久久九九免费精品| 日本黄色日本黄色录像| 国产人伦9x9x在线观看| 久久中文看片网| 精品久久久久久久久久免费视频 | 中国美女看黄片| 麻豆成人av在线观看| 免费搜索国产男女视频| 水蜜桃什么品种好| 亚洲国产欧美一区二区综合| 97人妻天天添夜夜摸| 久久亚洲真实| 亚洲人成网站在线播放欧美日韩| 18禁裸乳无遮挡免费网站照片 | 大陆偷拍与自拍| 一本综合久久免费| 久久人妻福利社区极品人妻图片| 香蕉久久夜色| 午夜亚洲福利在线播放| www日本在线高清视频| 狠狠狠狠99中文字幕| 亚洲色图av天堂| 久久精品91无色码中文字幕| 在线观看免费视频日本深夜| 99riav亚洲国产免费| 久久久久久免费高清国产稀缺| 欧美精品一区二区免费开放| 妹子高潮喷水视频| 国产极品粉嫩免费观看在线| 嫩草影视91久久| 很黄的视频免费| 亚洲国产毛片av蜜桃av| www国产在线视频色| 精品一品国产午夜福利视频| 国产av一区在线观看免费| 欧美黄色片欧美黄色片| 日韩大尺度精品在线看网址 | 久久人妻福利社区极品人妻图片| 亚洲va日本ⅴa欧美va伊人久久| 99国产精品一区二区蜜桃av| 亚洲欧美日韩另类电影网站| 新久久久久国产一级毛片| 最新在线观看一区二区三区| 久久国产乱子伦精品免费另类| 一级片'在线观看视频| 久久亚洲精品不卡| 在线播放国产精品三级| 亚洲视频免费观看视频| 人人妻人人添人人爽欧美一区卜| 在线永久观看黄色视频| 两个人看的免费小视频| 亚洲精品美女久久久久99蜜臀| 国产在线精品亚洲第一网站| 色综合站精品国产| 免费在线观看视频国产中文字幕亚洲| 丝袜美腿诱惑在线| 欧美乱色亚洲激情| 最近最新中文字幕大全免费视频| 黄色成人免费大全| 真人一进一出gif抽搐免费| 三上悠亚av全集在线观看| 免费不卡黄色视频| 日韩 欧美 亚洲 中文字幕| 正在播放国产对白刺激| 国产成年人精品一区二区 | 亚洲自拍偷在线| 精品久久久久久成人av|