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

    冶煉企業(yè)周邊農(nóng)田土壤的多環(huán)芳烴污染及其細(xì)菌群落效應(yīng)

    2016-12-12 01:41:57吳宇澄林先貴朱清禾曾軍丁慶旻
    生態(tài)毒理學(xué)報(bào) 2016年2期
    關(guān)鍵詞:芳烴農(nóng)田群落

    吳宇澄,林先貴,,朱清禾,曾軍,丁慶旻

    1. 中國科學(xué)院南京土壤研究所 土壤與農(nóng)業(yè)可持續(xù)發(fā)展國家重點(diǎn)實(shí)驗(yàn)室,南京 210008 2. 中國科學(xué)院南京土壤研究所-香港浸會(huì)大學(xué)土壤與環(huán)境聯(lián)合開放實(shí)驗(yàn)室,南京 210008

    ?

    冶煉企業(yè)周邊農(nóng)田土壤的多環(huán)芳烴污染及其細(xì)菌群落效應(yīng)

    吳宇澄1,2,林先貴1,2,,朱清禾1,2,曾軍1,2,丁慶旻1,2

    1. 中國科學(xué)院南京土壤研究所 土壤與農(nóng)業(yè)可持續(xù)發(fā)展國家重點(diǎn)實(shí)驗(yàn)室,南京 210008 2. 中國科學(xué)院南京土壤研究所-香港浸會(huì)大學(xué)土壤與環(huán)境聯(lián)合開放實(shí)驗(yàn)室,南京 210008

    多環(huán)芳烴是一類持久性有機(jī)污染物,進(jìn)入土壤后可能產(chǎn)生多方面生態(tài)效應(yīng)。為研究多環(huán)芳烴對(duì)土壤微生物的影響,選取南京某冶煉企業(yè)周邊農(nóng)田樣品,在分析污染物含量基礎(chǔ)上,采用高通量測(cè)序、定量PCR等方法綜合評(píng)價(jià)了土壤細(xì)菌多樣性和組成以及多環(huán)芳烴降解細(xì)菌豐度等特征。17個(gè)土壤樣品中,多環(huán)芳烴總量為0.25 ~ 31.08 mg·kg-1,并具有隨污染源距離增加而降低的空間分布特征。與土壤理化性質(zhì)如pH相比較,多環(huán)芳烴污染對(duì)土壤細(xì)菌的總體多樣性和群落組成影響不顯著。進(jìn)一步分析發(fā)現(xiàn)多環(huán)芳烴與潛在降解微生物的相對(duì)豐度和降解功能基因(芳香環(huán)羥基化雙加氧酶,PAH-RHDα)拷貝數(shù)顯著正相關(guān)。污染較重樣品的克隆、測(cè)序分析表明,土壤中PAH-RHDα基因主要屬于革蘭氏陽性細(xì)菌nidA3/fadA1類群,且與分支桿菌相關(guān)序列較為接近。這些結(jié)果綜合評(píng)價(jià)了冶煉企業(yè)周邊農(nóng)田土壤多環(huán)芳烴污染對(duì)微生物群落的影響,提示土壤污染在多環(huán)芳烴潛在降解細(xì)菌中的富集作用,將為后續(xù)污染土壤生物修復(fù)提供重要科學(xué)依據(jù)。

    多環(huán)芳烴;土壤細(xì)菌;芳香環(huán)羥基化雙加氧酶;分枝桿菌;群落效應(yīng)

    Received 26 October 2015 accepted 16 November 2015

    多環(huán)芳烴(polycyclic aromatic hydrocarbons,PAHs)是一類具有稠合苯環(huán)結(jié)構(gòu)的有機(jī)污染物,主要來自化石燃料,特別是煤的燃燒、車輛尾氣以及一些生物質(zhì)燃燒過程[1]。工業(yè)排放是環(huán)境中多環(huán)芳烴的重要來源,可導(dǎo)致附近農(nóng)田土壤的污染[2-3],并進(jìn)而通過食物鏈傳遞導(dǎo)致人體健康風(fēng)險(xiǎn)。

    多環(huán)芳烴對(duì)土壤生態(tài)系統(tǒng)具有廣泛的影響,如引起土壤動(dòng)物死亡[4],抑制植物生長[5],損害土壤功能如硝化等[6]。微生物作為土壤生態(tài)系統(tǒng)中廣泛存在、最為活躍的組成部分,和PAHs污染物之間存在復(fù)雜的相互作用。一方面,PAHs污染具有微生物群落效應(yīng),體現(xiàn)為對(duì)土壤微生物多樣性、組成及生理功能的潛在影響[7-8],反映污染物對(duì)微生物的潛在毒害效應(yīng);另一方面,部分微生物可通過多種酶學(xué)機(jī)制如芳香環(huán)羥基化雙加氧酶(PAH ring-hydroxylating dioxygenase, PAH-RHD)等轉(zhuǎn)化PAHs[9],是控制土壤有機(jī)污染物的關(guān)鍵因素。但是,目前對(duì)PAHs土壤微生物群落效應(yīng)的認(rèn)識(shí),多基于極端污染情況[10],未考慮污染物的老化過程[7],不能反映一般污染條件下的微生物響應(yīng)規(guī)律。

    冶煉企業(yè)是重要的多環(huán)芳烴工業(yè)排放源,其周邊土壤往往存在PAHs污染[11]。本研究系統(tǒng)采集南京市郊某冶煉企業(yè)周邊農(nóng)田土壤,測(cè)定污染物含量,采用高通量測(cè)序和定量PCR方法深入解析土壤細(xì)菌群落多樣性和組成,并分析潛在PAHs降解細(xì)菌的響應(yīng)情況。研究的主要目的在于分析工業(yè)排放導(dǎo)致農(nóng)田土壤PAHs污染的微生物群落效應(yīng),其結(jié)果將為評(píng)價(jià)PAHs的微生物生態(tài)毒性以及污染土壤生物修復(fù)提供科學(xué)依據(jù)。

    1 材料與方法(Materials and methods)

    1.1 樣品采集

    南京市郊某大型冶煉企業(yè),其煉焦、熱電等生產(chǎn)部門產(chǎn)生的廢氣排入大氣,是潛在的多環(huán)芳烴排放源。該企業(yè)東、南方主要為農(nóng)田,受到大氣沉降的影響,存在一定污染風(fēng)險(xiǎn)。以廠區(qū)為起點(diǎn),在該區(qū)域盛行下風(fēng)向沿偏東、偏南兩個(gè)方向開展采樣,共設(shè)置17個(gè)樣點(diǎn)(MS-1~MS-17,圖1)。其中,MS-1、MS-2最接近潛在污染源(熱電廠煙囪,約500 m);偏南方向包括MS-3至MS-10,最遠(yuǎn)處距污染源約5 km;偏東方向包括MS-11至MS-17,最遠(yuǎn)處距污染源約4.5 km。各采樣點(diǎn)均為蔬菜地。采樣于2014年12月進(jìn)行,采用手動(dòng)土壤采樣器按5點(diǎn)采樣法收集5 m×5 m范圍內(nèi)表層農(nóng)田土壤(0~10 cm),充分混勻成為一個(gè)混合土樣。土樣一部分經(jīng)風(fēng)干、磨細(xì)后用于PAHs及理化性質(zhì)分析,另一部分置于-20 ℃保存,用于微生物生態(tài)研究。

    圖1 冶煉企業(yè)周邊采樣點(diǎn)分布Fig. 1 Soil sampling sites around the steel plant

    1.2 土壤理化性質(zhì)測(cè)定

    采用酸度計(jì)測(cè)定土壤的pH值(水土比5:1);標(biāo)準(zhǔn)方法分析土壤總碳、總氮、全磷、全鉀[12];采用2 mol·L-1KCl溶液提取土壤無機(jī)氮(銨態(tài)氮、硝態(tài)氮),流動(dòng)分析儀(SKALAR)測(cè)定。

    1.3 土壤多環(huán)芳烴的提取與分析

    測(cè)定土壤中15種USEPA優(yōu)先控制多環(huán)芳烴的含量,即萘、苊、芴、菲、蒽、熒蒽、芘、苯并[a]蒽、屈、苯并[b]熒蒽、苯并[k]熒蒽、苯并[a]芘、茚并[1,2,3-cd]芘、二苯并蒽和苯并[g, h, i]苝。具體方法為:稱取10 g干土,用60 mL二氯甲烷在索氏裝置上提取24 h,30 ℃旋轉(zhuǎn)蒸干。加入2 mL環(huán)己烷溶解,取0.5 mL溶解液過硅膠柱,加入正己烷/二氯甲烷混合液(1:1)洗脫。洗脫液經(jīng)氮?dú)獯蹈?,用乙腈溶解并定容? mL。

    采用超快速液相色譜(島津UFLC-20系統(tǒng))檢測(cè)土壤多環(huán)芳烴的含量。流動(dòng)相為乙腈/水混合液,流速為0.8 mL·min-1,柱溫40 ℃。分離柱為反相Waters C18-PAH專用柱(4.6 mm×250 mm,粒徑5 μm),熒光檢測(cè)器檢測(cè)。

    1.4 土壤DNA提取

    采用FastDNA SPIN Kit for Soil試劑盒(MP Biomedicals)提取土壤DNA,微量分光光度計(jì)(NanoDrop 2000)和電泳法檢測(cè)DNA的質(zhì)量。為避免共提取的土壤腐殖質(zhì)等雜質(zhì)干擾PCR反應(yīng),DNA提取液經(jīng)10倍稀釋后用于下游分析。

    1.5 細(xì)菌16S rRNA基因高通量測(cè)序及數(shù)據(jù)分析

    采用通用引物519F和907R擴(kuò)增細(xì)菌16S rRNA基因片段。正向引物序列中包含5 bp的條形碼(barcode)序列。PCR反應(yīng)體系為50 μL,包括25 μL Taq DNA聚合酶預(yù)混液(Takara), 1 μL模板(約50 ng基因組DNA),1 μL正向及反向引物,23 μL ddH2O。PCR擴(kuò)增條件為95 ℃ 3 min;95 ℃ 45 s, 56 ℃ 45 s, 72 ℃ 60 s,35個(gè)循環(huán);72 ℃ 7 min。PCR產(chǎn)物經(jīng)純化后,構(gòu)建測(cè)序文庫,采用Illumina MiSeq系統(tǒng)進(jìn)行雙向高通量測(cè)序。

    基于QIIME分析平臺(tái)進(jìn)行高通量數(shù)據(jù)分析。序列經(jīng)拼接、比對(duì)后在97%相似性水平劃分操作分類單元(operational taxonomic unit,OTU),在此基礎(chǔ)上計(jì)算多樣性指數(shù),包括observed species,Chao1,Phylogenetic Diversity (PD),Shannon和Simpson多樣性指數(shù)等;通過與Greengene database數(shù)據(jù)庫比較確定序列的系統(tǒng)學(xué)分類。

    1.6 PAH-RHDα基因定量PCR分析

    參照文獻(xiàn)報(bào)道方法[13],采用定量PCR方法測(cè)定各樣品革蘭氏陽性和陰性細(xì)菌PAH-RHDα基因的拷貝數(shù)。定量PCR標(biāo)準(zhǔn)品的制備方法為:克隆目標(biāo)基因(具體方法見1.7),提取質(zhì)粒后測(cè)定濃度并計(jì)算拷貝數(shù),梯度稀釋制備標(biāo)準(zhǔn)曲線(拷貝數(shù)范圍107~101μL-1)。定量PCR標(biāo)準(zhǔn)曲線R2>0. 99,擴(kuò)增效率>80%。

    定量PCR采用Sybr Green方法,總反應(yīng)體系為25 μL,包含12.5 μL TransStart Green qPCR SuperMix(全式金),2 μL土壤DNA和10 pmol引物。反應(yīng)程序?yàn)槿椒ǎ?5 ℃預(yù)變性3 min;94 ℃ 15 s,55℃退火 45 s,95℃延伸30 s并讀取熒光信號(hào),共進(jìn)行40個(gè)循環(huán);隨后進(jìn)行熔解曲線分析和電泳評(píng)價(jià)擴(kuò)增單一性。

    1.7 革蘭氏陽性細(xì)菌PAH-RHDα基因克隆、測(cè)序與序列分析

    PCR擴(kuò)增革蘭氏陽性細(xì)菌的PAH-RHDα基因,經(jīng)割膠純化后克隆至pEASY-T1載體,并轉(zhuǎn)化感受態(tài)大腸桿菌。根據(jù)菌落顏色隨機(jī)挑選陽性克隆并測(cè)序,所獲序列經(jīng)去除載體和引物部分后,用Mothur軟件進(jìn)行比對(duì)和劃分操作分類單元。采用MEGA6軟件,以Maximum likelihood (ML)方法重建PAH-RHDα基因的系統(tǒng)發(fā)育樹。

    1.8 數(shù)據(jù)分析

    用SPSS13.0進(jìn)行相關(guān)性分析,雙側(cè)檢驗(yàn)判斷顯著性;CANOCO 4.5軟件進(jìn)行典范對(duì)應(yīng)分析(canonical correspondence analysis, CCA)。

    2 結(jié)果(Results)

    2.1 土壤理化性質(zhì)

    土壤理化性質(zhì)分析結(jié)果如表1所示。由各因子變異系數(shù)可見,17個(gè)樣品的全磷、銨態(tài)氮、硝態(tài)氮和PAHs總量性質(zhì)存在較大差異,而總碳、總氮、碳氮比和全鉀總體變化不大。pH值的變化幅度為4.74~7.66,表現(xiàn)出較高的土壤酸堿性差異。

    2.2 土壤PAHs含量及組成特征

    所有土壤樣品中均檢出多環(huán)芳烴,15種PAHs總量(Σ15PAHs)為0.25~31.08 mg·kg-1(表1)。最靠近排放源的MS-1和MS-2中,Σ15PAHs接近7 mg·kg-1。隨著距離增加,土壤PAHs含量迅速降低。在偏東方向,Σ15PAHs由MS-11的1.72 mg·kg-1逐漸降至MS-15至MS-17的0.30 mg·kg-1左右;在偏南方向,除MS-9外,PAHs總量處于0.51~0.96 mg·kg-1之間。所有樣品中,MS-9的PAHs總量最高,達(dá)31.08 mg·kg-1。

    圖2 熒蒽/(熒蒽+芘)和茚并[1,2,3-cd]芘/(茚并[1,2,3-cd]芘+苯并[g,h,i]苝)的十字交叉圖 注:多環(huán)芳烴比值范圍參照文獻(xiàn)[14-15]。Fig. 2 Cross plot of Fla/(Fla+Pyr) and IcdP/(IcdP+BghiP) ratio Note: The diagnostic ratios are referred to [14-15].

    表1 各樣點(diǎn)土壤基本性質(zhì)及多環(huán)芳烴含量

    按環(huán)數(shù)統(tǒng)計(jì)PAHs組成,接近污染源樣品(MS-1和MS-2)和偏南方向樣點(diǎn)(MS-3至MS-10)中均以4環(huán)及以上PAHs為主,偏東部分樣點(diǎn)中三環(huán)及以下的低環(huán)組分比例較高(表1)。計(jì)算熒蒽/(熒蒽+芘)(Fla/(Fla+Pyr))和茚并[1,2,3-cd]芘/(茚并[1,2,3-cd]芘+苯并[g,h,i]苝)(IcdP/(IcdP+BghiP))比值,并作十字交叉圖(圖2)。圖中大部分樣點(diǎn)聚為一群,體現(xiàn)出煤炭及化石燃料燃燒的特征;偏南方向5個(gè)樣點(diǎn)(MS-5,MS-7至MS-10)明顯偏離,其中MS-9部分呈現(xiàn)石油源特征。

    2.3 土壤細(xì)菌群落多樣性與組成

    17個(gè)土壤樣品共獲得22.5萬余條16S rRNA基因序列,各樣品序列數(shù)為8 786~20 433條。基于OTU計(jì)算各樣品多樣性指數(shù)并與環(huán)境因子作相關(guān)分析,結(jié)果顯示PAHs與細(xì)菌多樣性指數(shù)間無顯著相關(guān)性,在土壤因子中僅有pH與各多樣性指數(shù)間存在顯著正相關(guān)(表2)。

    各土壤在門水平的細(xì)菌群落組成如圖3所示。該區(qū)域農(nóng)田土壤的優(yōu)勢(shì)細(xì)菌有變形菌門(Proteobacteria)、擬桿菌門(Bacteroidetes)、放線菌門(Actinobacteria)、酸桿菌門(Acidobacteria)等。CCA排序圖中,第一和第二排序軸共解釋32.3%的群落變異,pH和總磷兩個(gè)因子對(duì)細(xì)菌組成影響顯著(P<0.05)。MS-9與MS-5的細(xì)菌群落組成與其他樣點(diǎn)分異較大,pH中性的土壤MS-1、MS-2、MS-3、MS-4、MS-10聚在一起,位于CCA圖的右側(cè),而其他更偏酸性的土壤則聚集在圖中的左側(cè)。

    在屬的水平,與Σ15PAHs顯著正相關(guān)(P<0.05),且至少一個(gè)樣品中相對(duì)豐度在0.5%以上的細(xì)菌有12個(gè),其中已有文獻(xiàn)報(bào)道具有PAHs降解能力的是分枝桿菌(Mycobacterium)[16]、紅球菌(Rhodococcus)[17]、Actinomadura[18]、Nocardioides[19]等4個(gè)。如果排除PAHs總量異常偏高的MS-9,則僅有分枝桿菌的相對(duì)豐度與Σ15PAHs顯著正相關(guān)(圖5a)。

    表2 多樣性指數(shù)與土壤性質(zhì)及多環(huán)芳烴總量間的相關(guān)性

    注:*和**分別表示顯著(P<0.05)和極顯著相關(guān)(P<0.01)。

    Note: * and ** indicate significance at the level of P<0.05 and P<0.01, respectively.

    圖4 細(xì)菌和環(huán)境因子的CCA分析Fig. 4 Canonical correspondence analysis of bacterial community in PAHs-contaminated soils

    2.4 PAH-RHDα基因豐度和組成

    定量PCR檢測(cè)中,各樣品中革蘭氏陰性菌均低于檢測(cè)限。革蘭氏陽性菌PAH-RHDα基因擴(kuò)增產(chǎn)物經(jīng)熔解曲線和電泳分析為單一峰(條帶),表明實(shí)驗(yàn)結(jié)果不受引物二聚體或非特異擴(kuò)增的影響。17個(gè)土壤樣品中,革蘭氏陽性菌PAH-RHDα基因拷貝數(shù)范圍為每克干土0.39×104~1.6×105拷貝,并與土壤中的多環(huán)芳烴總量顯著正相關(guān)(圖5b)。

    選擇PAH-RHDα基因豐度較高的MS-1、MS-2 和MS-9樣品進(jìn)行基因克隆和測(cè)序分析,共獲得55個(gè)序列,經(jīng)序列比對(duì)后在95%相似性劃分為4個(gè)OTU,其中OTU-1和OTU-2占所獲序列的91%。在系統(tǒng)發(fā)育樹上,這些序列與分支桿菌M. vanbaalenii和M. rhodesiae的nidA3B3基因最為接近,同屬于雙加氧酶的nidA3/fadA1類群(圖6)。

    圖5 PAHs含量與(a)分支桿菌相對(duì)豐度和(b)革蘭氏陽性細(xì)菌PAH-RHDα基因豐度的關(guān)系Fig. 5 Correlations between PAHs concentration and (a) relative abundance of Mycobacterium and (b) copies of GP PAH-RHDα gene

    圖6 三個(gè)土壤中革蘭氏陽性細(xì)菌PAH-RHDα基因的系統(tǒng)發(fā)育樹 注:本研究所獲序列以RHDα OTU開頭,括號(hào)中數(shù)字為OTU包含序列數(shù)。Fig. 6 Maximum likelihood tree of GP bacterial PAH-RHDα gene recovered from MS-1, MS-2 and MS-9 Note: The sequences recovered in this study begin with RHDα OTU. The numbers in the parentheses represent the number of sequences comprising the OTU.

    3 討論(Discussion)

    3.1 冶煉企業(yè)周邊農(nóng)田土壤PAHs污染特征

    本研究涉及17個(gè)農(nóng)田樣點(diǎn),距離污染源(熱電廠煙囪)距離從數(shù)百米至5 km不等。所有受試土壤中均檢出PAHs,含量變異較大。根據(jù)加拿大、荷蘭等國家現(xiàn)行土壤環(huán)境質(zhì)量標(biāo)準(zhǔn),接近排放源樣品中有4個(gè)(MS-1,MS-2,MS-3和MS-11)苯并[a]芘含量超出限定值(>0.1 mg·kg-1);以總量計(jì)算,接近排放源的樣品(MS-1至MS4,MS-11和MS12)接近或超過1 mg·kg-1(表1),污染程度相對(duì)較高。隨著距離的增加,PAHs總量迅速降低至1 mg·kg-1以內(nèi),呈現(xiàn)由近及遠(yuǎn)污染程度下降的趨勢(shì)。高污染土壤的Fla/(Fla+Pyr)比值接近燃煤特征值[15](圖2)。由于研究區(qū)域內(nèi)沒有其他高能耗企業(yè)或大量民用煤炭消耗,因此,該冶煉企業(yè)的排放可能是鄰近農(nóng)田土壤PAHs污染的主要原因。根據(jù)IcdP/(IcdP+BghiP)比值,偏東及偏南方向部分樣品PAHs具有化石燃料燃燒特征[15],這可能與這部分樣點(diǎn)臨近公路,受到車輛排放影響有關(guān)。MS-9樣品除PAHs總量外,苯并[a]芘含量高達(dá)4.0 mg·kg-1,是調(diào)查中所發(fā)現(xiàn)污染最嚴(yán)重的樣點(diǎn),由于其遠(yuǎn)離冶煉企業(yè),具體污染源尚待進(jìn)一步調(diào)查。

    3.2 多環(huán)芳烴對(duì)土壤細(xì)菌群落的影響

    高濃度多環(huán)芳烴污染土壤可以顯著改變微生物的多樣性[20]、群落組成[7]和生理活性[8]。但在農(nóng)田生態(tài)系統(tǒng)中,PAHs污染程度一般遠(yuǎn)低于前述研究[21-22];PAHs與土壤的相互作用影響其賦存形態(tài),將導(dǎo)致生物有效性的變化[23],這些因素可能導(dǎo)致不同的微生物響應(yīng)特征。本研究對(duì)某冶煉廠周邊農(nóng)田土壤中細(xì)菌群落的分析表明,相對(duì)于pH等土壤因子,PAHs污染未能明顯影響細(xì)菌多樣性和群落組成(表2,圖4),這可能與研究區(qū)域土壤污染水平相對(duì)不高和PAHs老化后生物有效性低有關(guān)。

    由于部分細(xì)菌具有PAHs降解功能,本研究重點(diǎn)分析了土壤中潛在PAHs降解菌的變化情況。從高通量測(cè)序數(shù)據(jù)中可以發(fā)現(xiàn)分支桿菌、紅球菌、Actinomadura、Nocardioides等革蘭氏陽性PAHs降解細(xì)菌的相對(duì)豐度與Σ15PAHs顯著正相關(guān)。另一方面,PAH-RHD是細(xì)菌降解PAHs的關(guān)鍵酶,其基因(PAH-RHDα)拷貝數(shù)通常反映了PAHs降解細(xì)菌的豐度[13],而高濃度多環(huán)芳烴往往導(dǎo)致土壤中PAH-RHDα基因的富集[24]。本研究未能檢出革蘭氏陰性菌PAH-RHDα基因,提示土壤中該類降解菌的數(shù)量較少。對(duì)PAHs含量較高土壤中革蘭氏陽性菌PAH-RHDα的克隆測(cè)序結(jié)果表明其組成單一,主要屬于nidA3/fadA1類群[25],與分支桿菌等細(xì)菌所含序列較為接近。因此,以分支桿菌為代表的革蘭氏陽性菌可能是該區(qū)域農(nóng)田土壤的主要PAHs降解微生物。

    微生物降解是土壤PAHs消減的主要機(jī)制,以分支桿菌為代表的PAHs降解細(xì)菌[26-27]已被用于污染土壤的生物強(qiáng)化修復(fù)。但是,微生物的降解能力受到包括土壤特性在內(nèi)的多種因素制約[28],不適宜的環(huán)境條件可能影響微生物強(qiáng)化等修復(fù)策略的效果。本研究的意義在于考察PAHs污染對(duì)農(nóng)田細(xì)菌群落的影響,主要揭示了自然污染狀態(tài)下PAHs降解細(xì)菌的富集,這些微生物可能較為適應(yīng)該冶煉企業(yè)周邊污染農(nóng)田的條件,在污染物降解中具有潛在作用。這些結(jié)果將為實(shí)施污染土壤生物修復(fù)提供重要的科學(xué)依據(jù)。

    [1] 傅家謨, 盛國英. 環(huán)境有機(jī)地球化學(xué)初探[J]. 地學(xué)前緣, 1996, 3(1-2): 127-132

    Fu J M, Sheng G Y. Preliminary study on environmental organic geochemistry [J]. Earth Science Frontiers, 1996, 3(1-2): 127-132 (in Chinese)

    [2] 田靖, 朱媛媛, 楊洪彪, 等. 大型鋼鐵廠及其周邊土壤多環(huán)芳烴污染現(xiàn)狀調(diào)查、評(píng)價(jià)與源解析[J]. 環(huán)境化學(xué), 2013, 32(6): 1002-1008

    Tian J, Zhu Y Y, Yang H B, et al. Investigation, assessment and source analysis of polycyclic aromatic hydrocarbons (PAHs) pollution in soil from a large iron and steel plant and its surrounding areas in China [J]. Environmental Chemistry, 2013, 32(6): 1002-1008 (in Chinese)

    [3] 葛成軍, 安瓊, 董元華, 等. 南京某地農(nóng)業(yè)土壤中有機(jī)污染分布狀況研究[J]. 長江流域資源與環(huán)境, 2006, 15(3): 361-365

    Ge C J, An Q, Dong Y H, et al. Distribution of organic pollutants in agricultural soil in Nanjing City [J]. Resources and Environment in the Yangtze Basin, 2006, 15(3): 361-365 (in Chinese)

    [4] Sverdrup L E, Nielsen T, Krogh P H. Soil ecotoxicity of polycyclic aromatic hydrocarbons in relation to soil sorption, lipophilicity, and water solubility [J]. Environmental Science & Technology, 2002, 36(11): 2429-2435

    [5] Maliszewska-Kordybach B, Smreczak B. Ecotoxicological activity of soils polluted with polycyclic aromatic hydrocarbons (PAHs) - effect on plants [J]. Environmental Technology, 2000, 21(10): 1099-1110

    [6] Sverdrup L E, Ekelund F, Kroghet P H, et al. Soil microbial toxicity of eight polycyclic aromatic compounds: Effects on nitrification, the genetic diversity of bacteria, and the total number of protozoans [J]. Environmental Toxicology and Chemistry, 2002, 21(8): 1644-1650

    [7] Sawulski P, Clipson N, Doyle E. Effects of polycyclic aromatic hydrocarbons on microbial community structure and PAH ring hydroxylating dioxygenase gene abundance in soil [J]. Biodegradation, 2014, 25(6): 835-847

    [8] de Menezes A, Clipson N, Doyle E. Comparative metatranscriptomics reveals widespread community responses during phenanthrene degradation in soil [J]. Environmental Microbiology, 2012, 14(9): 2577-2588

    [9] Fuchs G, Boll M, Heider J. Microbial degradation of aromatic compounds - from one strategy to four [J]. Nature Reviews Microbiology, 2011, 9(11): 803-816

    [10] Mukherjee S, Juottonen H, Siivonen P, et al. Spatial patterns of microbial diversity and activity in an aged creosote-contaminated site [J]. The ISME Journal, 2014, 8(10): 2131-2142

    [11] 彭曉春, 吳彥瑜, 謝莉. 廣東省醫(yī)療廢物焚燒廠周圍土壤多環(huán)芳烴特性[J]. 中國環(huán)境科學(xué), 2013, 33(S1): 108-122

    Peng X C, Wu Y Y, Xie L. Distribution characteristics of polycyclic aromatic hydrocarbons (PAHs) around medical waste incineration plant in Guangdong Province [J]. China Environmental Science, 2013, 33(S1): 108-122 (in Chinese)

    [12] 魯如坤. 土壤農(nóng)業(yè)化學(xué)分析方法[M]. 北京: 中國農(nóng)業(yè)科技出版社, 1999: 107-190

    [13] Cébron A, Norini M-P, Beguiristain T, et al. Real-time PCR quantification of PAH-ring hydroxylating dioxygenase (PAH-RHDα) genes from gram positive and gram negative bacteria in soil and sediment samples [J]. Journal of Microbiological Methods, 2008, 73(2): 148-159

    [15] Yunker M B, Macdonald R W, Vingarzan R, et al. PAHs in the Fraser River basin: A critical appraisal of PAH ratios as indicators of PAH source and composition [J]. Organic Geochemistry, 2002, 33(4): 489-515

    [16] Zeng J, Lin X G, Zhang J, et al. Isolation of polycyclic aromatic hydrocarbons (PAHs)-degrading Mycobacterium spp. and the degradation in soil [J]. Journal of Hazardous Materials, 2010, 183(1-3): 718-723

    [17] Sun G D, Xu Y, Liu Y, et al. Microbial community dynamics of soil mesocosms using Orychophragmus violaceus combined with Rhodococcus ruber Em1 for bioremediation of highly PAH-contaminated soil [J]. Applied Microbiology and Biotechnology, 2014, 98(24): 10243-10253

    [18] Subramanian A, Menon S. Novel polyaromatic hydrocarbon (PAH) degraders from oil contaminated soil samples [J]. International Journal of Advanced Research, 2015, 3(7): 999-1006

    [19] Saito A, Iwabuchi T, Harayama S. A novel phenanthrene dioxygenase from Nocardioides sp. strain KP7: Expression in Escherichia coli [J]. Journal of Bacteriology, 2000, 182(8): 2134-2141

    [20] Bengtsson G, T?rneman N, De Lipthay J R, et al. Microbial diversity and PAH catabolic genes tracking spatial heterogeneity of PAH concentrations [J]. Microbial Ecology, 2013, 65(1): 91-100

    [21] 劉增俊, 滕應(yīng), 黃標(biāo), 等. 長江三角洲典型地區(qū)農(nóng)田土壤多環(huán)芳烴分布特征與源解析[J]. 土壤學(xué)報(bào), 2010, 47(6): 1110-1117

    Liu Z J, Teng Y, Huang B, et al. Distribution and sources analysis of PAHs in farmland soils in areas typical of the Yangtze River Delta, China [J]. Acta Pedologica Sinica, 2010, 47(6): 1110-1117 (in Chinese)

    [22] Ping L F, Luo Y M, Zhang H B, et al. Distribution of polycyclic aromatic hydrocarbons in thirty typical soil profiles in the Yangtze River Delta Region, East China [J]. Environmental Pollution, 2007, 147(2): 358-365

    [23] Alexander M. Aging, bioavailability, and overestimation of risk from environmental pollutants [J]. Environmental Science & Technology, 2000, 34(20): 4259-4265

    [24] Chen S C, Peng J J, Duan G L. Enrichment of functional microbes and genes during pyrene degradation in two different soils [J]. Journal of Soils and Sediments, 2015, 16(2): 417-426

    [25] Kweon O, Kim S J, Baek S, et al. A new classification system for bacterial Rieske non-heme iron aromatic ring-hydroxylating oxygenases [J]. BMC Biochemistry, 2008, 9(1): 11

    [26] Li X, Hou L, Liu M, et al. Abundance and diversity of polycyclic aromatic hydrocarbon degradation bacteria in urban roadside soils in Shanghai [J]. Applied Microbiology and Biotechnology, 2015, 99(8): 3639-3649

    [27] Leys N M, Ryngaert A, Bastiaens L, et al. Occurrence and community composition of fast-growing Mycobacterium in soils contaminated with polycyclic aromatic hydrocarbons [J]. FEMS Microbiology Ecology, 2005, 51(3): 375-388

    [28] Bogan B W, Sullivan W R. Physicochemical soil parameters affecting sequestration and mycobacterial biodegradation of polycyclic aromatic hydrocarbons in soil [J]. Chemosphere, 2003, 52(10): 1717-1726

    Polycyclic Aromatic Hydrocarbons (PAHs) Pollution and Their Effects on Bacterial Community in Agricultural Soils Near a Smelting Plant

    Wu Yucheng1,2, Lin Xiangui1,2,*, Zhu Qinghe1,2, Zeng Jun1,2, Ding Qingmin1,2

    1. State Key Laboratory of Soil and Sustainable Agriculture, Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210008, China 2. Joint Open Laboratory of Soil and the Environment, Institute of Soil Science, Chinese Academy of Sciences and Hong Kong Baptist University, Nanjing 210008, China

    Polycyclic aromatic hydrocarbons (PAHs) are a group of persistent organic pollutants that may pose significant risks on soil biota. In this study, 17 soil samples were collected from agricultural sites potentially polluted by a smelting plant in the outskirt of Nanjing, Jiangsu Province. The distances of the sampling sites to the potential PAHs source ranged from 500 m to 5 km. The total amounts of 15 PAHs (Σ15PAHs) in these samples ranged from 0.25~31.08 mg·kg-1, with most contaminated samples (> 1 mg·kg-1) found in closely adjacent sites to the plant. High molecular weight PAHs (≥ 4 rings) were dominant in most samples. PAH diagnostic ratios suggested these PAHs were largely related to coal and fossil fuel combustion. To reveal the bacterial community composition, 16S rRNA gene was amplified and analyzed using next generation sequencing. The Illumina’s MiSeq sequencing produced more than 225 000 reads, with averagely 13 257 reads obtained for each sample. The most abundant phyla across all samples were Proteobacteria, Bacteroidetes, Actinobacteria and Acidobacteria. Bacterial alpha diversity, as measured by observed species, Chao1, Phylogenetic Diversity (PD), Shannon and Simpson indices, were directly correlated to pH rather than PAHs. An ordination analysis indicated that the bacterial community composition was significantly influenced by pH and total phosphorus, while the contribution of PAHs was minimal. However, the PAHs levels were positively correlated to the relative abundance of a few potential PAHs degraders such as Mycobacterium, Rhodococcus, Actinomadura and Nocardioides. This trend was further confirmed by the quantitative PCR (qPCR) quantification of bacterial PAH ring-hydroxylating dioxygenase genes (PAH-RHDα). Although the gram negative (GN) bacterial PAH-RHDα gene abundance was below the detectable level, gram positive (GP) PAH-RHDα was recovered from all samples and its abundance was positively correlated with the PAHs pollution. The GP PAH-RHDα gene in three selected samples affiliated to the nidA3/fadA1 group, and their closest matches in Genbank were largely derived from Mycobacterium. Overall, these findings indicate the influence of industrial PAHs emission on the adjacent agricultural soils. The PAHs pollution may cause the enrichment of specific PAHs degraders such as Mycobacterium, although soil pH could be more significant in shaping total community of soil bacteria.

    polycyclic aromatic hydrocarbons; soil bacteria; PAH-ring hydroxylating dioxygenase; Mycobacterium; community effect

    10.7524/AJE.1673-5897.20151026001

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)(2014CB441106);國家自然科學(xué)基金(41371310,41201301);江蘇省自然科學(xué)基金(BK20131462);土壤與農(nóng)業(yè)可持續(xù)發(fā)展國家重點(diǎn)實(shí)驗(yàn)室優(yōu)秀青年人才項(xiàng)目(Y212000014)

    吳宇澄(1977-),男,副研究員,研究方向?yàn)槲廴就寥郎镄迯?fù)及污染生態(tài)學(xué),E-mail: ycwu@issas.ac.cn

    *通訊作者(Corresponding author), E-mail: xglin@issas.ac.cn

    2015-10-26 錄用日期:2015-11-16

    1673-5897(2016)2-484-08

    Q142;X53

    A

    簡介:林先貴(1955-),男,研究員,博士生導(dǎo)師。主要研究方向土壤微生物多樣性及其生態(tài)功能、環(huán)境微生物及其應(yīng)用,發(fā)表學(xué)術(shù)論文300余篇。

    吳宇澄, 林先貴, 朱清禾, 等. 冶煉企業(yè)周邊農(nóng)田土壤的多環(huán)芳烴污染及其細(xì)菌群落效應(yīng)[J]. 生態(tài)毒理學(xué)報(bào),2016, 11(2): 484-491

    Wu Y C, Lin X G, Zhu Q H, et al. Polycyclic aromatic hydrocarbons (PAHs) pollution and their effects on bacterial community in agricultural soils near a smelting plant [J]. Asian Journal of Ecotoxicology, 2016, 11(2): 484-491 (in Chinese)

    猜你喜歡
    芳烴農(nóng)田群落
    大學(xué)生牙齦炎齦上菌斑的微生物群落
    達(dá)爾頓老伯的農(nóng)田
    合成微生物群落在發(fā)酵食品中的應(yīng)用研究
    關(guān)于重芳烴輕質(zhì)化與分離的若干思考
    農(nóng)田創(chuàng)意秀
    農(nóng)田搞養(yǎng)殖需辦哪些證
    農(nóng)田制作所
    輪胎填充油中8種多環(huán)芳烴的檢測(cè)
    高芳烴環(huán)保橡膠填充油量產(chǎn)
    環(huán)保型橡膠填充油中芳烴及多環(huán)芳烴組成研究
    a级毛片a级免费在线| 成人无遮挡网站| 亚洲欧美激情综合另类| 看片在线看免费视频| 精品不卡国产一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 成在线人永久免费视频| 日本撒尿小便嘘嘘汇集6| 精品午夜福利视频在线观看一区| 欧美日韩福利视频一区二区| av国产免费在线观看| 91麻豆精品激情在线观看国产| 欧美日本亚洲视频在线播放| 免费av不卡在线播放| 在线观看舔阴道视频| 两个人看的免费小视频| 校园春色视频在线观看| www国产在线视频色| 久久久国产成人精品二区| 欧美中文日本在线观看视频| e午夜精品久久久久久久| bbb黄色大片| 美女 人体艺术 gogo| 在线观看一区二区三区| 身体一侧抽搐| 麻豆久久精品国产亚洲av| 亚洲专区中文字幕在线| 久久人人精品亚洲av| 亚洲乱码一区二区免费版| 在线观看一区二区三区| 免费看十八禁软件| 久久人妻av系列| 国语自产精品视频在线第100页| 国产真人三级小视频在线观看| 国产精品av久久久久免费| aaaaa片日本免费| 成人欧美大片| 大型黄色视频在线免费观看| 中文字幕熟女人妻在线| 国产日本99.免费观看| 麻豆一二三区av精品| 免费在线观看影片大全网站| 国产精品日韩av在线免费观看| 最近在线观看免费完整版| 国产在线精品亚洲第一网站| 两个人的视频大全免费| 国产亚洲欧美在线一区二区| xxxwww97欧美| 久久久久免费精品人妻一区二区| 亚洲欧美日韩无卡精品| 久久久久久久精品吃奶| 欧美国产日韩亚洲一区| 日韩中文字幕欧美一区二区| 一二三四社区在线视频社区8| 日本一二三区视频观看| 久久久久久久久免费视频了| 一区二区三区高清视频在线| 99在线人妻在线中文字幕| 听说在线观看完整版免费高清| 好男人在线观看高清免费视频| 在线观看免费午夜福利视频| 淫妇啪啪啪对白视频| 亚洲成人免费电影在线观看| 天堂√8在线中文| 亚洲成av人片免费观看| 天堂av国产一区二区熟女人妻| 12—13女人毛片做爰片一| 亚洲精品美女久久久久99蜜臀| 日韩高清综合在线| 中文字幕久久专区| 嫩草影视91久久| 欧美日韩中文字幕国产精品一区二区三区| 亚洲成a人片在线一区二区| 亚洲在线自拍视频| 一本精品99久久精品77| 天堂动漫精品| 久久久国产精品麻豆| 国产精品久久久久久人妻精品电影| 青草久久国产| 国产伦精品一区二区三区四那| 99久久久亚洲精品蜜臀av| 哪里可以看免费的av片| 99久久精品国产亚洲精品| 国产伦一二天堂av在线观看| 一区福利在线观看| 久久性视频一级片| 亚洲中文日韩欧美视频| 日本黄色视频三级网站网址| 少妇熟女aⅴ在线视频| 久久久久国内视频| 成人午夜高清在线视频| 亚洲欧美日韩高清专用| 精品午夜福利视频在线观看一区| 国产成人av激情在线播放| 精品久久久久久久久久免费视频| 老司机午夜福利在线观看视频| 最新在线观看一区二区三区| 老司机午夜十八禁免费视频| 色视频www国产| 美女扒开内裤让男人捅视频| 亚洲av成人精品一区久久| 久久久久久久久中文| h日本视频在线播放| 日本黄大片高清| 亚洲精品一区av在线观看| 热99在线观看视频| 变态另类丝袜制服| 午夜亚洲福利在线播放| 一本一本综合久久| 欧美中文综合在线视频| 超碰成人久久| 欧美黑人欧美精品刺激| xxxwww97欧美| 国产黄片美女视频| 天堂√8在线中文| 亚洲欧美激情综合另类| 无限看片的www在线观看| 啦啦啦免费观看视频1| 国产高清有码在线观看视频| 黑人欧美特级aaaaaa片| 午夜福利在线在线| 午夜福利在线观看吧| 国产 一区 欧美 日韩| 亚洲熟女毛片儿| www.999成人在线观看| 亚洲人与动物交配视频| 中文资源天堂在线| 少妇丰满av| 美女cb高潮喷水在线观看 | 亚洲 国产 在线| 嫩草影院精品99| 久久精品国产99精品国产亚洲性色| 亚洲av五月六月丁香网| 亚洲熟妇中文字幕五十中出| 桃色一区二区三区在线观看| 日本成人三级电影网站| 欧美性猛交黑人性爽| 老司机午夜福利在线观看视频| 国产精品一区二区精品视频观看| 日韩欧美三级三区| 国产单亲对白刺激| 九九在线视频观看精品| 免费大片18禁| 高潮久久久久久久久久久不卡| 日韩人妻高清精品专区| 一夜夜www| 国产av一区在线观看免费| 日本黄大片高清| 天天添夜夜摸| 国产亚洲精品综合一区在线观看| 久久久久久人人人人人| 性色av乱码一区二区三区2| 成人鲁丝片一二三区免费| 免费在线观看视频国产中文字幕亚洲| 日韩欧美精品v在线| 少妇裸体淫交视频免费看高清| 99久久综合精品五月天人人| 成年女人看的毛片在线观看| 视频区欧美日本亚洲| 夜夜夜夜夜久久久久| 国产精品乱码一区二三区的特点| 午夜福利欧美成人| 18禁裸乳无遮挡免费网站照片| 夜夜躁狠狠躁天天躁| 午夜视频精品福利| 少妇裸体淫交视频免费看高清| 色综合欧美亚洲国产小说| 婷婷精品国产亚洲av| 一进一出好大好爽视频| 久久这里只有精品19| 搞女人的毛片| 久久这里只有精品中国| 久久久色成人| 欧美高清成人免费视频www| 在线观看免费视频日本深夜| 无遮挡黄片免费观看| 亚洲性夜色夜夜综合| 欧美xxxx黑人xx丫x性爽| 一本精品99久久精品77| 成人av一区二区三区在线看| 91av网站免费观看| 国产精品一区二区免费欧美| 精品久久久久久,| 一夜夜www| 国内久久婷婷六月综合欲色啪| 午夜两性在线视频| 国产三级黄色录像| 成人三级做爰电影| 国产高清视频在线播放一区| 欧美在线一区亚洲| 国产精品野战在线观看| 午夜福利欧美成人| 成人av在线播放网站| 精品久久久久久久人妻蜜臀av| 人人妻,人人澡人人爽秒播| 精品免费久久久久久久清纯| 久久热在线av| 日本与韩国留学比较| 性色av乱码一区二区三区2| 日本在线视频免费播放| 不卡一级毛片| 在线观看美女被高潮喷水网站 | 日韩欧美精品v在线| 夜夜爽天天搞| 午夜日韩欧美国产| 久久精品国产亚洲av香蕉五月| 欧美成狂野欧美在线观看| 国产精品电影一区二区三区| av国产免费在线观看| 男插女下体视频免费在线播放| 国产1区2区3区精品| 男人舔女人下体高潮全视频| 亚洲成a人片在线一区二区| 亚洲熟女毛片儿| 国内毛片毛片毛片毛片毛片| 午夜精品一区二区三区免费看| 色在线成人网| 国产成人av教育| 免费无遮挡裸体视频| 亚洲熟妇熟女久久| 国产 一区 欧美 日韩| 亚洲美女黄片视频| 成人av在线播放网站| 曰老女人黄片| 女生性感内裤真人,穿戴方法视频| 婷婷丁香在线五月| av天堂在线播放| 男人的好看免费观看在线视频| 91在线观看av| 免费高清视频大片| 国产又黄又爽又无遮挡在线| 人人妻人人澡欧美一区二区| 99国产精品一区二区蜜桃av| 久久精品国产99精品国产亚洲性色| 看免费av毛片| 免费av不卡在线播放| 亚洲精品美女久久久久99蜜臀| 在线免费观看不下载黄p国产 | 黑人操中国人逼视频| 99视频精品全部免费 在线 | 97超级碰碰碰精品色视频在线观看| 国产精品久久电影中文字幕| 国产精品99久久99久久久不卡| 亚洲国产精品成人综合色| 国产精品美女特级片免费视频播放器 | 岛国在线免费视频观看| 在线免费观看不下载黄p国产 | 黄色 视频免费看| 啦啦啦免费观看视频1| 制服丝袜大香蕉在线| 夜夜夜夜夜久久久久| 免费看光身美女| 精品国产三级普通话版| 国产97色在线日韩免费| 国产爱豆传媒在线观看| 国产精华一区二区三区| 麻豆成人av在线观看| 观看美女的网站| 亚洲欧洲精品一区二区精品久久久| 久久人人精品亚洲av| 免费观看的影片在线观看| 757午夜福利合集在线观看| 波多野结衣高清作品| 亚洲成人久久爱视频| 亚洲成av人片免费观看| 91av网站免费观看| 在线观看日韩欧美| 老司机深夜福利视频在线观看| 国产欧美日韩一区二区三| 国产成人av教育| 国产精品1区2区在线观看.| 久久精品综合一区二区三区| 国产97色在线日韩免费| 亚洲欧美精品综合一区二区三区| www.999成人在线观看| 一进一出抽搐动态| 高清在线国产一区| 午夜福利高清视频| 中文字幕熟女人妻在线| 熟女少妇亚洲综合色aaa.| 国产精品一区二区精品视频观看| 亚洲av美国av| 国产视频内射| 亚洲美女视频黄频| or卡值多少钱| 又爽又黄无遮挡网站| 欧美不卡视频在线免费观看| 他把我摸到了高潮在线观看| 一进一出抽搐动态| 婷婷亚洲欧美| 老汉色av国产亚洲站长工具| e午夜精品久久久久久久| 中文字幕高清在线视频| 精品电影一区二区在线| 免费看美女性在线毛片视频| www国产在线视频色| 全区人妻精品视频| 久久久久久久午夜电影| 久久久久国产一级毛片高清牌| 国产一级毛片七仙女欲春2| 亚洲精品美女久久久久99蜜臀| 丝袜人妻中文字幕| 国产精品一区二区免费欧美| www.自偷自拍.com| 88av欧美| 久久午夜亚洲精品久久| 国产伦在线观看视频一区| 中文字幕人成人乱码亚洲影| 男女视频在线观看网站免费| 国产精品久久电影中文字幕| 此物有八面人人有两片| 偷拍熟女少妇极品色| 中出人妻视频一区二区| 757午夜福利合集在线观看| 国产91精品成人一区二区三区| 亚洲成av人片免费观看| 丁香六月欧美| 嫩草影视91久久| 免费看美女性在线毛片视频| 午夜成年电影在线免费观看| 亚洲中文字幕一区二区三区有码在线看 | 国产精品自产拍在线观看55亚洲| 午夜激情欧美在线| 欧美乱色亚洲激情| 热99在线观看视频| 久久午夜亚洲精品久久| 国产成人精品久久二区二区免费| 欧美又色又爽又黄视频| 男插女下体视频免费在线播放| 人人妻人人看人人澡| 国产高清激情床上av| 欧美黑人巨大hd| 天堂av国产一区二区熟女人妻| 最好的美女福利视频网| 丝袜人妻中文字幕| 一级毛片女人18水好多| 精品人妻1区二区| 老鸭窝网址在线观看| 99精品在免费线老司机午夜| 神马国产精品三级电影在线观看| 色精品久久人妻99蜜桃| 亚洲av片天天在线观看| 亚洲美女视频黄频| 一进一出抽搐动态| 黄色日韩在线| 久久国产乱子伦精品免费另类| 日韩高清综合在线| 欧美激情久久久久久爽电影| 天堂√8在线中文| 亚洲精品乱码久久久v下载方式 | 每晚都被弄得嗷嗷叫到高潮| 成人特级av手机在线观看| 特大巨黑吊av在线直播| 麻豆成人av在线观看| 五月玫瑰六月丁香| 99热6这里只有精品| 国产亚洲精品久久久久久毛片| 久久久水蜜桃国产精品网| 成人鲁丝片一二三区免费| 国产精品久久久久久久电影 | 国产精品一及| 可以在线观看的亚洲视频| 超碰成人久久| 热99在线观看视频| 亚洲人成电影免费在线| 免费在线观看视频国产中文字幕亚洲| 精品国产美女av久久久久小说| 久久久色成人| 成人高潮视频无遮挡免费网站| 床上黄色一级片| 成人18禁在线播放| 亚洲美女黄片视频| 国产亚洲av嫩草精品影院| 成人国产综合亚洲| 嫩草影院精品99| 天堂网av新在线| 国产主播在线观看一区二区| 日本撒尿小便嘘嘘汇集6| 身体一侧抽搐| 天天躁日日操中文字幕| 中文字幕高清在线视频| 在线免费观看的www视频| 无人区码免费观看不卡| 白带黄色成豆腐渣| 最近视频中文字幕2019在线8| 日韩欧美国产在线观看| 成人三级黄色视频| 狂野欧美白嫩少妇大欣赏| 麻豆成人av在线观看| 精品免费久久久久久久清纯| 一级毛片高清免费大全| 91麻豆精品激情在线观看国产| 亚洲第一欧美日韩一区二区三区| 亚洲一区二区三区色噜噜| 久久久久久大精品| 日本三级黄在线观看| 午夜福利高清视频| 听说在线观看完整版免费高清| 精品无人区乱码1区二区| 国语自产精品视频在线第100页| 他把我摸到了高潮在线观看| 久久久国产成人免费| 日本熟妇午夜| 老汉色av国产亚洲站长工具| 黄色丝袜av网址大全| 午夜福利18| 亚洲色图 男人天堂 中文字幕| 天天一区二区日本电影三级| 亚洲精品粉嫩美女一区| 午夜久久久久精精品| 国产一区二区在线观看日韩 | 国产成人精品无人区| 黑人操中国人逼视频| 老汉色av国产亚洲站长工具| www.熟女人妻精品国产| 91av网一区二区| 久久这里只有精品19| 色哟哟哟哟哟哟| 国模一区二区三区四区视频 | 两个人视频免费观看高清| 国产人伦9x9x在线观看| 国产成人系列免费观看| 精品无人区乱码1区二区| 久久精品91蜜桃| 国模一区二区三区四区视频 | 母亲3免费完整高清在线观看| 欧美日本视频| 欧美黄色片欧美黄色片| 老司机福利观看| 99热这里只有是精品50| 亚洲乱码一区二区免费版| 在线观看日韩欧美| 午夜福利在线观看免费完整高清在 | 国产日本99.免费观看| 成年女人看的毛片在线观看| 香蕉丝袜av| 国产综合懂色| 最好的美女福利视频网| 午夜福利在线观看吧| 亚洲国产精品合色在线| 国产亚洲欧美在线一区二区| 中出人妻视频一区二区| 黄色 视频免费看| 日日干狠狠操夜夜爽| 夜夜躁狠狠躁天天躁| 美女午夜性视频免费| 精品熟女少妇八av免费久了| 免费av不卡在线播放| av在线天堂中文字幕| 久久久精品大字幕| 啦啦啦韩国在线观看视频| 日韩精品青青久久久久久| 中亚洲国语对白在线视频| 国产精品一区二区三区四区免费观看 | 欧美极品一区二区三区四区| 日韩欧美精品v在线| 国产麻豆成人av免费视频| 久久精品91蜜桃| 午夜福利在线观看吧| 不卡av一区二区三区| 久久久久久国产a免费观看| 人妻丰满熟妇av一区二区三区| 99久久精品热视频| 99热这里只有精品一区 | а√天堂www在线а√下载| 国产伦精品一区二区三区四那| 国产午夜精品久久久久久| xxx96com| 最新中文字幕久久久久 | 99久久99久久久精品蜜桃| 久久香蕉国产精品| 国产激情欧美一区二区| 国产精品,欧美在线| www国产在线视频色| 999久久久国产精品视频| 丝袜人妻中文字幕| 欧美日韩精品网址| 91av网站免费观看| 欧美高清成人免费视频www| 天堂√8在线中文| 国产又黄又爽又无遮挡在线| 国产精品香港三级国产av潘金莲| 女人被狂操c到高潮| 波多野结衣高清无吗| 男人和女人高潮做爰伦理| 色老头精品视频在线观看| 黄色视频,在线免费观看| 天堂av国产一区二区熟女人妻| www.999成人在线观看| 精品久久久久久成人av| 成人三级黄色视频| 精品无人区乱码1区二区| 亚洲成人久久性| 欧美乱妇无乱码| 亚洲专区国产一区二区| 嫩草影院入口| h日本视频在线播放| 国产成人精品无人区| 国产午夜精品论理片| 国产一区二区三区在线臀色熟女| 精品久久久久久久人妻蜜臀av| АⅤ资源中文在线天堂| 91老司机精品| 99久久精品热视频| 日本在线视频免费播放| 国产精品av视频在线免费观看| 亚洲欧美日韩高清专用| 国产亚洲精品久久久久久毛片| 国产成人影院久久av| 亚洲欧美精品综合一区二区三区| 国产高清有码在线观看视频| 日本熟妇午夜| 岛国在线免费视频观看| 国内精品一区二区在线观看| 日韩欧美免费精品| 成人永久免费在线观看视频| 91在线精品国自产拍蜜月 | 亚洲专区字幕在线| 综合色av麻豆| 国产av麻豆久久久久久久| 亚洲精品色激情综合| 免费观看精品视频网站| 亚洲国产欧洲综合997久久,| 久久久精品大字幕| 长腿黑丝高跟| 啦啦啦观看免费观看视频高清| 波多野结衣高清作品| 国产在线精品亚洲第一网站| 免费电影在线观看免费观看| 中文字幕久久专区| 国产午夜精品久久久久久| 久久久久久久久免费视频了| 黄片小视频在线播放| 午夜久久久久精精品| 国产成年人精品一区二区| 久久久精品大字幕| 国产三级黄色录像| 亚洲av中文字字幕乱码综合| 欧美性猛交黑人性爽| 亚洲中文日韩欧美视频| a级毛片a级免费在线| 亚洲国产欧美网| 亚洲性夜色夜夜综合| 国产精品亚洲美女久久久| 免费大片18禁| 久久久久国产精品人妻aⅴ院| 国产精品 国内视频| 黑人巨大精品欧美一区二区mp4| 国模一区二区三区四区视频 | 18禁裸乳无遮挡免费网站照片| 国产视频一区二区在线看| 成年人黄色毛片网站| 国产高潮美女av| 可以在线观看的亚洲视频| 亚洲精品在线观看二区| av片东京热男人的天堂| 18禁黄网站禁片免费观看直播| 亚洲男人的天堂狠狠| 九九久久精品国产亚洲av麻豆 | 国内精品美女久久久久久| 亚洲国产色片| 村上凉子中文字幕在线| 校园春色视频在线观看| 一进一出好大好爽视频| 欧美最黄视频在线播放免费| 国产精品久久久久久人妻精品电影| 男人舔女人的私密视频| 国产人伦9x9x在线观看| 日韩欧美免费精品| 1000部很黄的大片| 午夜成年电影在线免费观看| 舔av片在线| 中文亚洲av片在线观看爽| 国产成人欧美在线观看| 九色成人免费人妻av| 国语自产精品视频在线第100页| x7x7x7水蜜桃| 国产精品,欧美在线| 亚洲自偷自拍图片 自拍| 别揉我奶头~嗯~啊~动态视频| 亚洲美女视频黄频| 国产熟女xx| 真人一进一出gif抽搐免费| 国内揄拍国产精品人妻在线| 99久久99久久久精品蜜桃| 免费观看的影片在线观看| 在线视频色国产色| 啦啦啦观看免费观看视频高清| 久久亚洲真实| 神马国产精品三级电影在线观看| 九色国产91popny在线| 欧美成人性av电影在线观看| 伊人久久大香线蕉亚洲五| 黄色丝袜av网址大全| 观看美女的网站| 老熟妇仑乱视频hdxx| 草草在线视频免费看| 国产亚洲精品一区二区www| 99在线视频只有这里精品首页| 在线观看一区二区三区| 欧美日本视频| 99在线视频只有这里精品首页| 久久精品国产亚洲av香蕉五月| 香蕉久久夜色| 99精品久久久久人妻精品| 久久久成人免费电影| 亚洲自偷自拍图片 自拍| 熟女人妻精品中文字幕| 色吧在线观看| 国产精品一区二区三区四区免费观看 | 又爽又黄无遮挡网站| 国产黄a三级三级三级人| 国产视频一区二区在线看| 特大巨黑吊av在线直播| 久久精品综合一区二区三区|