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

    青島市不同功能區(qū)冬季空氣微生物群落代謝與多樣性特征

    2015-03-13 01:43:36吳等等宋志文徐愛(ài)玲
    生態(tài)學(xué)報(bào) 2015年7期
    關(guān)鍵詞:海濱填埋場(chǎng)水源地

    吳等等,宋志文,徐愛(ài)玲,鄭 遠(yuǎn),夏 巖

    青島理工大學(xué), 青島 266033

    青島市不同功能區(qū)冬季空氣微生物群落代謝與多樣性特征

    吳等等,宋志文*,徐愛(ài)玲,鄭 遠(yuǎn),夏 巖

    青島理工大學(xué), 青島 266033

    選取青島市5個(gè)功能區(qū)(市區(qū)街道、海濱區(qū)域、飲用水源地、垃圾填埋場(chǎng)和人工濕地污水處理系統(tǒng)),采用SAS ISO100空氣浮游菌采樣器于2013年冬季采集空氣微生物樣品,應(yīng)用BIOLOG方法分析空氣微生物群落代謝功能多樣性,闡明群落代謝與環(huán)境相關(guān)性。結(jié)果表明,不同功能區(qū)空氣微生物群落碳源代謝強(qiáng)度存在差異,代謝穩(wěn)定時(shí),海濱區(qū)域和飲用水源地樣品平均光密度值(AWCD)分別為0.302、0.210,而人工濕地、市區(qū)街道及垃圾填埋場(chǎng)分別為0.063、0.025和0.034,海濱區(qū)域和飲用水源地空氣微生物群落碳源代謝強(qiáng)度明顯高于其他功能區(qū)。不同功能區(qū)空氣微生物群落Shannon指數(shù)和Simpson指數(shù)接近,但海濱區(qū)域和飲用水源地McIntosh指數(shù)明顯高于其他功能區(qū)。海濱區(qū)域和飲用水源地空氣微生物群落碳源代謝類型豐富,代謝水平高,人工濕地、市區(qū)街道和垃圾填埋場(chǎng)碳源代謝類型單一,代謝水平低。5個(gè)功能區(qū)空氣微生物群落碳源代謝差異呈現(xiàn)區(qū)域性,分異代謝差異的主要是羧酸類碳源。風(fēng)速、溫度、濕度等非生物因素對(duì)空氣微生物群落碳源代謝具有不同程度影響,且不同功能區(qū)主導(dǎo)非生物因素存在差異。BIOLOG方法可以提供大量多維數(shù)據(jù),能夠分析樣品間微生物群落碳源代謝差異,客觀、全面表征空氣微生物群落碳源代謝多樣性特征,是研究空氣微生物群落功能多樣性較理想的方法之一。

    青島市; 城市功能區(qū); 空氣微生物; 代謝特征; 群落多樣性; BIOLOG方法

    空氣微生物是指空氣中細(xì)菌、霉菌和放線菌等有生命的活體,主要來(lái)源于土壤、水體、動(dòng)植物和人類,此外污水處理、動(dòng)物飼養(yǎng)、發(fā)酵過(guò)程和農(nóng)業(yè)活動(dòng)等也是空氣微生物的重要來(lái)源[1- 4]。已知存在空氣中的細(xì)菌及放線菌有1200種,真菌有40000種[5]??諝馕⑸锊粌H具有重要的生態(tài)系統(tǒng)功能,還與空氣污染、環(huán)境質(zhì)量和人體健康密切相關(guān)[6]。了解城市空氣微生物群落特征、物種組成及濃度,對(duì)控制城市微生物污染、改善環(huán)境質(zhì)量、維持人群健康和提高工農(nóng)業(yè)生產(chǎn)具有重要的理論和實(shí)際意義。國(guó)內(nèi)外學(xué)者已經(jīng)對(duì)城市不同功能區(qū)、垃圾填埋場(chǎng)、污水處理廠、養(yǎng)殖環(huán)境等空氣微生物區(qū)系、時(shí)空分布特征進(jìn)行研究[1,7- 14]。但這些研究大多采用傳統(tǒng)微生物培養(yǎng)方法[15- 20],由于“可培養(yǎng)類”微生物僅占總微生物的不到1%[21- 22],并且還有部分微生物處于活的不可培養(yǎng)狀態(tài),加上培養(yǎng)結(jié)果受培養(yǎng)基組分影響較大,導(dǎo)致測(cè)得的數(shù)據(jù)與實(shí)際情況有較大偏差。為了較準(zhǔn)確的反映空氣微生物群落和生態(tài)功能,需要從空氣微生物種類、群落結(jié)構(gòu)、功能多樣性以及生物量等層次來(lái)研究。新發(fā)展的基于微生物群落代謝單一碳源特征的BIOLOG方法能夠彌補(bǔ)傳統(tǒng)方法的不足,底物碳源的多樣性決定了BIOLOG微平板可以提供大量多維數(shù)據(jù),且能夠比較樣品間微生物群落差異,較全面地表征微生物群落結(jié)構(gòu)、總體代謝活性與功能信息,該方法已在微生物鑒定、污水處理工藝、污染土壤修復(fù)中得到廣泛應(yīng)用。

    本文以青島市5個(gè)城市功能區(qū)(市區(qū)街道、海濱區(qū)域、飲用水源地、城市生活垃圾填埋場(chǎng)和人工濕地污水處理系統(tǒng))為研究對(duì)象,應(yīng)用BIOLOG方法分析空氣微生物群落碳源代謝功能多樣性,闡明空氣微生物群落代謝與環(huán)境相關(guān)性。對(duì)了解沿海城市不同功能區(qū)空氣微生物特征及其對(duì)大氣環(huán)境和人體健康的影響具有重要意義。

    1 研究地區(qū)與研究方法

    1.1 實(shí)驗(yàn)地點(diǎn)概況

    選擇青島市區(qū)街道(撫順路青島理工大學(xué)段)、海濱區(qū)域(魯迅公園)、飲用水源地(棘洪灘水庫(kù))、垃圾填埋場(chǎng)(小澗西生活垃圾綜合處理廠)、人工濕地(膠南人工濕地污水處理系統(tǒng))5個(gè)功能區(qū)作為采樣地點(diǎn)。其中,(1)撫順路位于青島市市北區(qū),是青島市主要交通樞紐,人流量和車流量較大,汽車尾氣和揚(yáng)塵污染較嚴(yán)重;(2)魯迅公園位于青島市市南區(qū),南面緊鄰大海,長(zhǎng)約1 km,占地面積0.04 km2,綠化面積75%以上,是青島4A級(jí)旅游景點(diǎn),冬季人流量較少;(3)棘洪灘水庫(kù)位于青島市城陽(yáng)區(qū),是引黃濟(jì)青工程唯一調(diào)蓄水庫(kù),庫(kù)區(qū)面積14.4 km2,是青島市主要飲用水源地,供應(yīng)青島市50%以上居民用水,周圍是農(nóng)田,人流量較少;(4)小澗西生活垃圾綜合處理廠位于青島市城陽(yáng)區(qū),是青島市內(nèi)六區(qū)唯一的生活垃圾填埋場(chǎng),設(shè)計(jì)庫(kù)容7.1×106m3,占地面積0.66 km2,填埋區(qū)面積0.27 km2,綠地率41.25%,實(shí)際日填埋量3500 t/d,垃圾滲濾液經(jīng)垂直和水平收集系統(tǒng)收集后排入污水處理廠進(jìn)行處理;(5)人工濕地污水處理系統(tǒng)位于青島市黃島區(qū),東南面緊鄰黃海,采用自由表面流蘆葦濕地工藝,由99個(gè)并聯(lián)運(yùn)行濕地單元組成,每個(gè)單元大小為140 m×32 m,總占地面積76.7 km2,處理規(guī)模3×104m3/d,污水來(lái)源為生活污水和工業(yè)廢水(約1∶1)。

    1.2 研究方法

    1.2.1 樣品采集

    采樣時(shí)間為2013年1月1日—5日,采用SAS ISO 100空氣微生物采樣器采集空氣微生物樣品,空氣流量100 L/min,采樣時(shí)間30 min。采樣后取下凝膠膜,加到100 mL滅菌生理鹽水中,搖床震蕩10 min,凝膠膜溶解后作為接種液。采樣期間利用TASI- 620數(shù)字式溫濕度計(jì)測(cè)定采樣環(huán)境溫度和濕度,利用TASI- 641風(fēng)速儀測(cè)定風(fēng)速(表1)。

    表1 采樣點(diǎn)環(huán)境因子Table 1 Environmental factors of sampling sites

    1.2.2 BIOLOG-GN操作方法

    接種混勻的接種液于BIOLOG-GN板中,每孔接種量125 μL,接種后的BIOLOG-GN板置于25 ℃恒溫培養(yǎng)箱中培養(yǎng)[23],每隔24 h于BIOLOG讀數(shù)器上讀數(shù),波長(zhǎng)590 nm,直至光密度值穩(wěn)定為止,整個(gè)培養(yǎng)過(guò)程共持續(xù)8d。

    1.3 計(jì)算方法

    1.3.1 微生物碳源代謝水平

    平均光密度值(AWCD)計(jì)算方法[23]:AWCD=∑(C-R)/n,式中C為反應(yīng)孔光密度值,R為對(duì)照孔光密度值,n為底物數(shù)量,對(duì)于不同BIOLOG微平板n值不同, BIOLOG-GN板為95。

    1.3.2 微生物多樣性計(jì)算方法

    選用72 h光密度值計(jì)算5個(gè)不同功能區(qū)空氣微生物群落多樣性指數(shù)[24]: Shannon指數(shù)(H′)、McIntosh指數(shù)(U)和Simpson指數(shù)(D)。計(jì)算公式見(jiàn)表2。

    表2 空氣微生物群落多樣性指數(shù)計(jì)算方法Table 2 Computational methods for diversity index of airborne microbial community

    1.3.3 統(tǒng)計(jì)分析方法

    選取72 h平均光密度值(AWCD)利用SPSS 17.0軟件進(jìn)行主成分分析(PCA),利用CANOCA軟件進(jìn)行典范對(duì)應(yīng)分析(CCA)[25]。

    2 結(jié)果與分析

    2.1 空氣微生物群落碳源代謝強(qiáng)度變化

    平均光密度值(AWCD)能夠衡量空氣微生物利用不同碳源的整體能力,表征微生物生理活性,從代謝水平上揭示空氣微生物群落功能多樣性。青島市5個(gè)功能區(qū)AWCD隨時(shí)間延長(zhǎng)而增大,最初48 h內(nèi)曲線較平緩,72 h AWCD值進(jìn)入指數(shù)增長(zhǎng)期,192 h達(dá)到穩(wěn)定狀態(tài)(圖1)。方差分析表明,海濱區(qū)域和飲用水源地空氣微生物碳源代謝水平與其余各功能區(qū)存在顯著性差異(P<0.05)。培養(yǎng)過(guò)程中,人工濕地、市區(qū)街道和垃圾填埋場(chǎng)空氣微生物碳源代謝水平較低,在192 h時(shí)AWCD值仍小于0.1,而海濱區(qū)域和飲用水源地空氣微生物碳源代謝水平相對(duì)較高。5個(gè)功能區(qū)AWCD值變化趨勢(shì)不同,表明各功能區(qū)空氣微生物碳源利用能力、豐度等存在差異。

    2.2 空氣微生物群落碳源利用水平分析

    BIOLOG-GN板包含95種碳源,根據(jù)碳源官能團(tuán)可將其劃為6類,其中糖類(Carbohydrates)30種、羧酸類(Carboxylic acids)24種、氨基酸類(Amino acids)20種、胺/氨類(Amines/amides)6種、聚合物類(Polymers)5種、其他(Miscellaneous)10種[26]。5個(gè)功能區(qū)空氣微生物對(duì)6類碳源利用程度存在差異,海濱區(qū)域和飲用水源地空氣微生物碳源代謝水平高于其他功能區(qū)??諝馕⑸锾荚创x優(yōu)勢(shì)群落依次為羧酸類代謝群落、糖類代謝群落、氨基酸類代謝群落。5個(gè)功能區(qū)空氣微生物對(duì)羧酸類和糖類碳源利用程度均較高,對(duì)其它化合物類利用程度較低,其中,對(duì)胺/氨類代謝水平最低,市區(qū)街道、垃圾填埋場(chǎng)和飲用水源地空氣微生物基本不能代謝這類碳源。同一功能區(qū)對(duì)不同碳源的利用程度也存在差異(圖2)。

    圖1 青島市5個(gè)功能區(qū)空氣微生物AWCD變化Fig.1 Average well color development (AWCD) changes of airborne microbes at five city functional regions in Qingdao City

    圖2 青島市5個(gè)功能區(qū)空氣微生物群落對(duì)BIOLOG平板碳源的相對(duì)利用率Fig.2 Relative utilization ratios of carbon sources in BIOLOG plate by the airborne microbes at five city functional regions in Qingdao City

    2.3 空氣微生物群落種群多樣性分析

    采用Shannon指數(shù)、Simpson指數(shù)和McIntosh指數(shù)分析5個(gè)功能區(qū)空氣微生物群落功能多樣性,3種多樣性指數(shù)反映了空氣微生物群落功能多樣性不同側(cè)面,其中,Shannon指數(shù)受群落物種豐富度影響較大,反映了微生物群落物種種數(shù)及種間個(gè)體分配的均勻性;Simpson指數(shù)反映了群落中最常見(jiàn)的物種,是評(píng)估微生物群落優(yōu)勢(shì)度的指標(biāo);McIntosh指數(shù)則是度量群落物種均一度的指標(biāo)??梢钥闯觯琒hannon指數(shù)飲用水源地最低,其次為海濱區(qū)域,其余3個(gè)功能區(qū)基本相同,Simpson指數(shù)除飲用水源地較低外,其余4個(gè)功能區(qū)基本相同,McIntosh指數(shù)差異較大,飲用水源地和海濱區(qū)域明顯高于其余3個(gè)功能區(qū)(表3)。

    表3 不同功能區(qū)空氣微生物群落多樣性指數(shù)Table 3 Diversity index of airborne microbial community at five city functional regions in Qingdao City

    2.4 空氣微生物群落代謝功能主成分分析(PCA)

    利用主成分分析法(PCA)分析青島市5個(gè)功能區(qū)空氣微生物群落代謝功能,共提取2個(gè)主成分,其中主成分1(PC1)方差貢獻(xiàn)率37.6%,主成分2(PC2)方差貢獻(xiàn)率23.0%,累計(jì)方差貢獻(xiàn)率60.6%,能夠解釋原有變量大部分變異。在PC1上載荷大于0.9的基質(zhì)有15種,主要是羧酸類、氨基酸及糖類,其中D-葡萄糖二酸載荷最大(0.968);在PC2上載荷大于0.9的基質(zhì)有9種,主要是氨基酸和糖類,其中D-葡萄糖- 6-磷酸載荷最大(0.984)(表4)。PC1、PC2上載荷大于0.9的基質(zhì)類型和利用率較高的碳源類型一致。

    表4 與PC1、PC2相關(guān)系數(shù)大于0.9的基質(zhì)

    Table 4 Correlation coefficients greater than 0.9 between carbon substrates and Principal component 1(PC1) or Principal component 2 (PC2)

    基質(zhì)Carbonsource主成分1Principalcomponent1PC1(r)基質(zhì)Carbonsource主成分2Principalcomponent2PC2(r)D-葡萄糖二酸0.968D-葡萄糖-6-磷酸0.984肝糖0.960α-環(huán)式糊精0.978L-谷氨酸0.959L-蘇氨酸0.976N-乙酰-D-半乳糖胺0.958L-鳥(niǎo)氨酸0.972L-丙氨酸0.955L-絲氨酸0.966L-脯氨酸0.954D-阿拉伯醇0.965吐溫800.953麥芽糖0.963α-戊酮酸0.948尿苷0.939乙酸0.944松二糖0.900α-D-乳糖0.931D-葡萄糖0.925奎寧酸/金雞鈉酸0.914γ-羥丁酸0.904L-海藻糖0.902i-赤蘚糖醇0.903

    根據(jù)青島市5個(gè)功能區(qū)空氣微生物群落碳源代謝相似性可以將5個(gè)功能區(qū)歸為3類,其中市區(qū)街道、人工濕地及垃圾填埋場(chǎng)歸為一類;海濱區(qū)域歸為一類;飲用水源地歸為一類。在PC1上,3類區(qū)域空氣微生物碳源代謝差異顯著,碳源利用譜容易區(qū)分。考察成分載荷矩陣可知,羧酸類在PC1上載荷較大,是碳源代謝差異區(qū)域性的主要分異碳源類型(圖3)。

    2.5 非生物因素對(duì)空氣微生物群落碳源代謝水平的影響

    研究表明,非生物因素(溫度、濕度、海拔、風(fēng)速、風(fēng)向等)對(duì)空氣微生物群落碳源代謝強(qiáng)度和單一碳源利用能力有不同程度影響[27- 29]。CCA典范對(duì)應(yīng)分析能夠從統(tǒng)計(jì)學(xué)角度分析環(huán)境因子和空氣微生物群落碳源代謝相關(guān)性,揭示不同環(huán)境條件下微生物群落結(jié)構(gòu)變化原因。從青島市5個(gè)功能區(qū)空氣微生物群落碳源代謝強(qiáng)度與環(huán)境因子CCA分析排序圖中可以看出,溫度、濕度與人工濕地、垃圾填埋場(chǎng)、市區(qū)街道空氣微生物群落碳源代謝水平呈正相關(guān),與海濱區(qū)域、飲用水源地呈負(fù)相關(guān);風(fēng)速與人工濕地、垃圾填埋場(chǎng)、市區(qū)街道空氣微生物群落碳源代謝水平呈負(fù)相關(guān),與海濱區(qū)域、飲用水源地呈正相關(guān)(圖4)。

    圖3 5個(gè)功能區(qū)空氣微生物群落BIOLOG微平板中底物碳源代謝的主成分分析Fig.3 Principal Component analysis (PCA) of BIOLOG plate data profiles revealed the patterns of microbial communities at five city functional regions in Qingdao City

    圖4 5個(gè)功能區(qū)碳源代謝強(qiáng)度與環(huán)境因子CCA分析排序Fig.4 A two-dimensional graph analysis of canonical correspondence analysis (CCA) coordination for AWCD of five airborne microbial community and environment factors

    3 討論

    不同環(huán)境條件或不同環(huán)境介質(zhì)微生物群落碳源代謝強(qiáng)度、碳源代謝類型差異顯著,微生物群落碳源代謝呈現(xiàn)區(qū)域性,并且分異這種區(qū)域性的碳源類型隨環(huán)境條件和環(huán)境介質(zhì)改變而不同[23,27,29- 31]。

    3.1 空氣微生物群落碳源代謝強(qiáng)度變化

    AWCD可以評(píng)估空氣微生物利用不同碳源的整體能力,從功能代謝水平上揭示空氣微生物群落結(jié)構(gòu)多樣性,是反映空氣微生物活性、描述空氣微生物群落功能多樣性的重要指標(biāo)[29]。青島市5個(gè)功能區(qū)空氣微生物群落碳源代謝強(qiáng)度存在差異,其中,海濱區(qū)域、飲用水源地空氣微生物群落碳源代謝強(qiáng)度較高,人工濕地、市區(qū)街道和垃圾填埋場(chǎng)空氣微生物群落碳源代謝強(qiáng)度較低。研究表明,高輻射、寡營(yíng)養(yǎng)及干燥環(huán)境使得一些生長(zhǎng)緩慢、產(chǎn)孢子和芽孢的類群(如放線菌和芽孢桿菌)占優(yōu)勢(shì),導(dǎo)致碳源利用能力降低[32],而B(niǎo)IOLOG方法對(duì)生長(zhǎng)緩慢、不能很好利用BIOLOG板中碳源的微生物反映能力較差[33]。海濱區(qū)域和飲用水源地環(huán)境濕度較大,導(dǎo)致空氣微生物生理活性和碳源利用強(qiáng)度相對(duì)較高。對(duì)比BIOLOG方法在土壤微生物群落代謝中的研究結(jié)果[34- 43],發(fā)現(xiàn)空氣微生物群落碳源代謝整體水平較低,可能與空氣微生物組成及采樣季節(jié)有關(guān)[44]。

    3.2 空氣微生物群落多樣性

    Shannon指數(shù)(H′)、McIntosh指數(shù)(U)和Simpson指數(shù)(D)分別從微生物群落豐富度、均一性及常見(jiàn)種群優(yōu)勢(shì)度等方面反映空氣微生物群落多樣性。青島市不同功能區(qū)空氣微生物群落豐富度和優(yōu)勢(shì)度接近,但群落均一度存在差異,飲用水源地和海濱區(qū)域空氣微生物群落物種較均一,可能由于該功能區(qū)環(huán)境濕度大,晝夜溫差小,為一些微生物種群提供適宜生存條件[27,29],同時(shí),冬季海濱區(qū)域和飲用水源地受人為因素影響較少,揚(yáng)塵等對(duì)空氣微生物群落組成擾動(dòng)較小,從而形成相對(duì)均一穩(wěn)定的空氣微生物群落。

    3.3 空氣微生物群落碳源代謝水平差異

    青島市5個(gè)功能區(qū)空氣微生物群落碳源代謝類型和代謝水平存在差異,且呈現(xiàn)區(qū)域性。海濱區(qū)域空氣微生物群落碳源代謝類型主要是糖類、氨基酸類、羧酸類及聚合物類,對(duì)糖類代謝水平最高;飲用水源地碳源代謝類型主要是聚合物類、羧酸類、氨基酸類及糖類,對(duì)聚合物類碳源代謝水平最高;市區(qū)街道、人工濕地和垃圾填埋場(chǎng)空氣微生物群落主要碳源代謝類型是羧酸類、糖類;5個(gè)功能區(qū)空氣微生物群落對(duì)羧酸類碳源代謝水平均較高。研究發(fā)現(xiàn),空氣微生物群落碳源代謝水平與所處環(huán)境條件具有相關(guān)性,環(huán)境條件影響空氣微生物的來(lái)源,同時(shí)構(gòu)成微生物生存的特殊生境,影響其生理活性,使空氣微生物群落碳源代謝差異呈現(xiàn)區(qū)域性[17,19]。冬季人工濕地、垃圾填埋場(chǎng)和市區(qū)街道均有濕度低、溫差大的特點(diǎn),并且垃圾填埋場(chǎng)和市區(qū)街道空氣微生物群落組成受人為活動(dòng)影響較大,但研究發(fā)現(xiàn)人工濕地、市區(qū)街道和垃圾填埋場(chǎng)空氣微生物群落碳源代謝相近,推測(cè)這3個(gè)功能區(qū)空氣微生物群落組成及生理活性受溫度、濕度等環(huán)境因子影響較大。海濱區(qū)域和飲用水源地盡管有著相似的溫度、濕度條件,但由于海水環(huán)境和淡水環(huán)境差異,形成空氣微生物生境差異,同時(shí)海風(fēng)也會(huì)影響空氣微生物群落組成,導(dǎo)致海濱區(qū)域和飲用水源地空氣微生物群落碳源代謝存在差異。

    4 結(jié)論

    (1)青島市5個(gè)功能區(qū)冬季空氣微生物群落碳源代謝強(qiáng)度和代謝類型存在差異。海濱區(qū)域和飲用水源地空氣微生物群落碳源代謝強(qiáng)度明顯高于其他功能區(qū),且碳源代謝類型豐富,代謝水平相對(duì)較高。

    (2)青島市5個(gè)功能區(qū)冬季空氣微生物豐富度和種群優(yōu)勢(shì)度接近,而群落物種均一度存在差異??諝馕⑸锶郝涮荚创x差異呈現(xiàn)區(qū)域性,可以歸為3類,分別為:①海濱區(qū)域、②飲用水源地、③人工濕地、市區(qū)街道和垃圾填埋場(chǎng)。羧酸類是分異不同功能區(qū)空氣微生物碳源代謝差異區(qū)域性的主要碳源類型。

    (3)BIOLOG方法能夠從微生物群落碳源代謝強(qiáng)度、碳源代謝類型及碳源利用區(qū)域性分異等方面研究微生物群落代謝,但只能反映在BIOLOG-GN板上能夠生長(zhǎng)的微生物。要全面研究空氣微生物群落結(jié)構(gòu)特征,還需要結(jié)合基因多樣性方面的研究。

    [1] Bovallius A, Bucht B, Roffey R, An?s P. Three-Year investigation of the natural airborne bacterial flora at four localities in Sweden. Applied and Environmental Microbiology, 1978, 35(5): 847- 852.

    [2] Lindmanm J, Constantinidou H, Barchet W, Christen D. Plants as sources of airborne bacteria including ice nucleation-active bacteria. Applied and Environmental Microbiology, 1982, 44(5): 1059- 1063.

    [3] Mancinelli R, Shulls W A. Airborne bacteria in an urban environment. Applied and Environmental Microbiology, 1978, 35(6): 1095- 1011.

    [4] Martinez K F, Sheehy J W, Jones J H, Cusick L B. Microbial containment in conventional fermentation processes. Applied Industrial Hygiene, 1988, 3(6): 177- 181.

    [5] 周連玉, 喬楓. 中國(guó)綠地空氣微生物研究現(xiàn)狀與展望. 中國(guó)農(nóng)學(xué)通報(bào), 2011, 27(20): 269- 273.

    [6] Lacey J, Dutkiewicz J. Bioaerosols and occupational lung disease. Journal of Aerosol Science, 1994, 25(8): 1371- 1404.

    [7] Wright T J, Greene V W, Paulus H J. Viable microorganisms in an urban atmosphere. Journal of the Air Pollution Control Associate, 1969, 19(5): 337- 341.

    [8] 方治國(guó), 歐陽(yáng)志云, 胡利鋒, 王效科, 林學(xué)強(qiáng). 北京市夏季空氣微生物群落結(jié)構(gòu)和生態(tài)分布. 生態(tài)學(xué)報(bào), 2005, 25(1): 83- 88.

    [9] 岳冰, 夏立江, 杜文利, 鄺溯琳, 李軼倫, 許建民. 北京某垃圾填埋場(chǎng)空氣微生物污染狀況. 中國(guó)環(huán)境科學(xué), 2006, 26(S1): 102- 105.

    [10] 任啟文, 王成, 楊穎, 胡立香. 城市綠地空氣微生物濃度研究——以北京元大都公園為例. 干旱區(qū)資源與環(huán)境, 2007, 21(4): 80- 83.

    [11] 蔡旭玲, 劉堅(jiān), 姚振江, 周俊立, 陳思東. 廣州部分大型汽車客運(yùn)站候車廳空氣微生物污染狀況. 環(huán)境與健康雜志, 2010, 27(8): 725.

    [12] 王友斌, 康懷雄, 鄧俊興, 申志新, 李素紅, 邱賀民, 劉秀巖, 孫宗科, 李霞, 張偉, 丁培, 陳西平. 商場(chǎng)與超市室內(nèi)空氣微生物污染狀況調(diào)查. 中國(guó)衛(wèi)生檢驗(yàn)雜志, 2011, 21(8): 2038- 2040.

    [13] 孫宗科, 康懷雄, 申志新, 李素紅, 吳建榮, 唐建紅, 盧志堅(jiān), 張淑紅, 邱賀民, 張偉, 陳西平. 南方和北方城市住宅室內(nèi)空氣微生物調(diào)查. 環(huán)境與健康雜志, 2011, 28(2): 130- 132.

    [14] 歐陽(yáng)友生, 謝小保, 陳儀本, 黃小茉, 彭紅, 施慶珊. 廣州市空氣微生物含量及其變化規(guī)律研究. 微生物學(xué)通報(bào), 2006, 33(3): 47- 51.

    [15] Choi D S, Park Y K, Oh S K, Yoon H J, Kim J C, Seo W J, Cha S H. Distribution of airborne microorganisms in Yellow Sands of Korea. The Journal of Microbiology, 1997, 35(1): 1- 9.

    [16] Yeo H G, Kim J H. SPM and fungal spores in the ambient air of west Korea during the Asian dust (Yellow Sand) period. Atmospheric Environment, 2002, 36(35): 5437- 5442.

    [17] Griffin D W, Kubilay N, Ko?ak M, Gray M A, Borden T C, Shinn E A. Airborne desert dust and aeromicrobiology over the Turkish Mediterranean coastline. Atmospheric Environment, 2007, 41(19): 4050- 4062.

    [18] Hua N P, Kobayashi F, Iwasaka Y, Shi G Y, Naganuma T. Detailed identification of desert-originated bacteria carried by Asian dust storms to Japan. Aerobiologia, 2007, 23(4): 291- 298.

    [19] Polymenakou P N, Mandalakis M, Stephanou E G, Tselepides A. Particle size distribution of airborne microorganisms and pathogens during an intense African dust event in the eastern Mediterranean. Environmental Health Perspectives, 2008, 116(3): 292- 296.

    [20] Meyer C, Bernard N, Moskura M, Toussaint M L, Denayer F, Gilbert D. Effects of urban particulate deposition on microbial communities living in bryophytes: An experimental study. Ecotoxicology and Environmental Safety, 2010, 73(7): 1776- 1784.

    [21] Brock T D. The study of microorganismsinsitu: Progress and problems. Symposium of the Society for General Microbiology, 1987, 41: 1- 17.

    [22] Amann R I, Ludwig W, Schleifer K H. Phylogenetic identification and in situ detection of individual microbial cells without cultivation. Microbiology Reviews, 1995, 59(1): 143- 169.

    [23] 鄭華, 陳法霖, 歐陽(yáng)志云, 方治國(guó), 王效科, 苗鴻. 不同森林土壤微生物群落對(duì)BIOLOG-GN板碳源的利用. 環(huán)境科學(xué), 2007, 28(5): 1126- 1130.

    [24] Jeffrey S B, Donald D K. Microbial diversity in the rhizosphere of corn grown under conventional and low-input systems. Applied Soil Ecology, 1996, 5(1): 21- 27.

    [25] 王強(qiáng), 戴九蘭, 吳大千, 余悅, 申天琳, 王仁卿. 微生物生態(tài)研究中基于BIOLOG方法的數(shù)據(jù)分析. 生態(tài)學(xué)報(bào), 2010, 30(3): 817- 823.

    [26] Preston-Mafham J, Boddy L, Randerson P F. Analysis of microbial community functional diversity using sole-carbon-source utilisation profiles-a critique. FEMS Microbiology Ecology, 2002, 42(1): 1- 14.

    [27] 段魏魏, 婁愷, 曾軍, 胡蓉, 史應(yīng)武, 何清, 劉新春, 孫建, 晁群芳. 塔克拉瑪干沙塵暴源區(qū)空氣微生物群落的代謝特征. 環(huán)境科學(xué), 2012, 33(1): 27- 31.

    [28] 凌琪, 王晏平. 黃山風(fēng)景區(qū)夏季空氣微生物分布特征初步研究. 微生物學(xué)通報(bào), 2008, 35(9): 1379- 1384.

    [29] 凌琪, 包金梅, 李瑞, 陶勇, 鮑立寧, 毛欽焱. BIOLOG-ECO解析黃山風(fēng)景區(qū)空氣微生物碳代謝多樣性特征. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào), 2012, 20(1): 56- 63.

    [30] 鄭瓊, 崔曉陽(yáng), 邸雪穎, 金森. 不同林火強(qiáng)度對(duì)大興安嶺偃松林土壤微生物功能多樣性的影響. 林業(yè)科學(xué), 2012, 48(5): 96- 100.

    [31] 姜睿玲, 楊統(tǒng)一, 陳芳艷, 唐玉斌. 桑園土壤微生物群落功能多樣性對(duì)PAHs污染的響應(yīng). 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào), 2012, 28(4): 439- 444.

    [32] Griffin D W. Atmospheric movement of microorganisms in clouds of desert dust and implications for human health. Clinical Microbiology Reviews, 2007, 20(3): 459- 477.

    [33] 李鑫, 張會(huì)慧, 岳冰冰, 金薇薇, 許楠, 朱文旭, 孫廣玉. 桑樹(shù)-大豆間作對(duì)鹽堿土碳代謝微生物多樣性的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2012, 23(7): 1825- 1831.

    [34] 張瑞, 張貴龍, 陳冬青, 姬艷艷, 張海芳, 楊殿林. 不同施肥對(duì)農(nóng)田土壤微生物功能多樣性的影響. 中國(guó)農(nóng)學(xué)通報(bào), 2013, 29(2): 133- 139.

    [35] 邵麗, 谷潔, 張社奇, 高華, 秦清軍. 生物復(fù)混肥對(duì)土壤微生物功能多樣性及土壤酶活性的影響. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2012, 31(6): 1153- 1159.

    [36] Xue D, Yao H Y, Ge D Y, Huang C Y. Soil microbial community structure in diverse land use systems: a comparative study using Biolog, DGGE, and PLFA Analyses. Pedosphere, 2008, 18(5): 653- 663.

    [37] 曹均, 吳姬, 趙小蓉, 李貴桐, 孫明德, 曹慶昌, 林啟美. 北京9個(gè)典型板栗園土壤碳代謝微生物多樣性特征. 生態(tài)學(xué)報(bào), 2010, 30(2): 527- 532.

    [38] 孔濱, 孫波, 鄭憲清, 陳小云, 隋躍宇, 王簾里. 水熱條件和施肥對(duì)黑土中微生物群落代謝特征的影響. 土壤學(xué)報(bào), 2009, 46(1): 100- 106.

    [39] 安韶山, 李國(guó)輝, 陳利頂. 寧南山區(qū)典型植物根際與非根際土壤微生物功能多樣性. 生態(tài)學(xué)報(bào), 2011, 31(18): 5225- 5234.

    [40] 楊宇虹, 陳冬梅, 晉艷, 王海斌, 段玉琦, 郭徐魁, 何海斌, 林文雄. 不同肥料種類對(duì)連作煙草根際土壤微生物功能多樣性的影響. 作物學(xué)報(bào), 2011, 37(1): 105- 111.

    [41] 劉繼明, 黃炳軍, 徐演鵬, 盧萍, 林英華. 河南雞公山自然保護(hù)區(qū)典型森林類型土壤微生物群落功能多樣性. 林業(yè)資源管理, 2013, (1): 76- 79,85.

    [42] 陳法霖, 鄭華, 陽(yáng)柏蘇, 張凱, 肖燚, 屠乃美. 中亞熱帶幾種針、闊葉樹(shù)種凋落物混合分解對(duì)土壤微生物群落碳代謝多樣性的影響. 生態(tài)學(xué)報(bào), 2011, 31(11): 3027- 3035.

    [43] 姬艷艷, 張貴龍, 張瑞, 劉玉升, 楊殿林, 王彩靈. 耕作方式對(duì)農(nóng)田土壤微生物功能多樣性的影響. 中國(guó)農(nóng)學(xué)通報(bào), 2013, 29(6): 117- 123.

    [44] 潘立勇, 李勇, 付紅, 潘宏博. 徐州城市主要環(huán)境功能區(qū)空氣微生物污染狀況評(píng)價(jià). 中國(guó)環(huán)境監(jiān)測(cè), 2008, 24(3): 72- 74.

    Metabolic characteristics and community diversities of airborne microbes at different functional regions in Qingdao in winter

    WU Dengdeng, SONG Zhiwen*, XU Ailing, ZHENG Yuan, XIA Yan

    CollegeofEnvironmentandMunicipalEngineering,QingdaoTechnologicalUniversity,Qingdao266033,China

    Airborne microbes are major biological components of the ecosystem that have important ecosystem functions, affect human health and air pollution, and are key indicators of air quality. To determine the metabolic characteristics and community diversity of airborne microbes in different functional regions of Qingdao in winter, sampling sites were set up in five different functional regions (urban streets, a coastal area, drinking water source area, municipal landfill, and artificial wetlands). Airborne microbes were then collected using an SAS ISO100 air sampler in January 2013 and their carbon source metabolic characteristics, functional diversity, and relationship with environmental factors were systematically analyzed by the BIOLOG method. The results revealed significant differences in carbon metabolic profiles of air microbial communities from the five locations. When metabolism was stable, the average well color developments of samples from the coastal area and drinking water source area were 0.302 and 0.21, respectively, whereas those of artificial wetlands, urban streets, and municipal landfill were 0.063, 0.025, and 0.034, respectively. Therefore, the levels of carbon metabolism in the coastal area and drinking water source area were higher than those at other locations. The Shannon indexes and Simpson indexes at the five different functional regions were similar. However, the McIntosh indexes at the coastal area and drinking water source area were higher than those at other locations. Overall, the microbial diversity and dominant species differed among functional regions, but microflora was distributed more evenly in the coastal area and drinking water source area. Among the five locations, carbon catabolic types and levels in the coastal area and drinking water source area were richer and higher than those in the artificial wetlands, urban streets, and municipal landfill. Overall, the airborne microbes in various functional areas had high carbohydrate and carboxylic acid metabolism, and air microbial communities in the coastal and drinking water source areas had better ability to metabolize polymers, carbohydrates, carboxylic acids, and amino acids than those in the artificial wetlands, urban streets, and municipal landfill. The characteristics of carbon metabolism revealed regional differences that were mainly caused by carboxylic acids. Specifically, urban streets, municipal landfills, and artificial wetlands had similar metabolic characteristics, and could be classified together. Additionally, coastal and drinking water source areas had different characteristics, and could be classified respectively. Environmental factors such as wind speed, temperature, and humidity may affect carbon utilization to a certain degree; however, the dominant factors will differ among environments. The results of canonical correspondence analysis (CCA) showed that wind speed and the carbon metabolism of the coastal area and drinking water source area were positively correlated, while they were negatively correlated with temperature and humidity. However, urban streets, municipal landfill, and artificial wetlands were positively correlated with temperature and humidity and negatively correlated with wind speed. The BIOLOG method could provide a large amount of multi-dimensional data, compare carbon metabolic characteristics of the microbial community between samplers, and reflect characteristics of the carbon metabolic diversity of airborne microbes objectively and generally, making it an ideal method for studying functional diversity of airborne microbes.

    Qingdao City; city functional regions; airborne microbial community; metabolic characteristics; microbial diversity; BIOLOG method

    國(guó)家自然科學(xué)基金項(xiàng)目(31170509)

    2013- 06- 05;

    日期:2014- 05- 08

    10.5846/stxb201306051339

    *通訊作者Corresponding author.E-mail: songzhiwen@qtech.edu.cn

    吳等等,宋志文,徐愛(ài)玲,鄭遠(yuǎn),夏巖.青島市不同功能區(qū)冬季空氣微生物群落代謝與多樣性特征.生態(tài)學(xué)報(bào),2015,35(7):2277- 2284.

    Wu D D, Song Z W, Xu A L, Zheng Y, Xia Y.Metabolic characteristics and community diversities of airborne microbes at different functional regions in Qingdao in winter.Acta Ecologica Sinica,2015,35(7):2277- 2284.

    猜你喜歡
    海濱填埋場(chǎng)水源地
    夏日海濱
    多功能區(qū)固廢填埋場(chǎng)建設(shè)案例分析
    新型GCL復(fù)合垂直防滲技術(shù)及其在填埋場(chǎng)的應(yīng)用
    海濱書(shū)簡(jiǎn)
    海濱1
    基于因果分析圖法的飲用水源地保護(hù)探討
    南水北調(diào)中線水源地土壤侵蝕經(jīng)濟(jì)損失估算
    吉林省供水水源地問(wèn)題探討
    礦山疏干排水對(duì)地下水源地影響的數(shù)值模擬
    Effect of decompression in different time on hemodynamics and oxygen metabolism of porcine model with severe acute pancreatitis combined intraabdominal hypertension
    国产精品电影一区二区三区| 99在线视频只有这里精品首页| 丝袜美腿诱惑在线| 午夜亚洲福利在线播放| 欧美日韩亚洲综合一区二区三区_| 亚洲av第一区精品v没综合| а√天堂www在线а√下载| 午夜福利在线观看吧| 麻豆国产av国片精品| 国产视频内射| 欧美一级毛片孕妇| 成熟少妇高潮喷水视频| 一本精品99久久精品77| 久9热在线精品视频| 欧美日韩福利视频一区二区| 黄色a级毛片大全视频| 精品不卡国产一区二区三区| 亚洲,欧美精品.| 色精品久久人妻99蜜桃| 欧美久久黑人一区二区| 麻豆久久精品国产亚洲av| 午夜a级毛片| 最近最新中文字幕大全免费视频| www.精华液| 亚洲国产精品合色在线| 国产主播在线观看一区二区| 男插女下体视频免费在线播放| 午夜成年电影在线免费观看| 国产精品亚洲一级av第二区| 一进一出好大好爽视频| 青草久久国产| 国产久久久一区二区三区| 又黄又爽又免费观看的视频| 久久99热这里只有精品18| 69av精品久久久久久| 老司机午夜福利在线观看视频| 亚洲国产日韩欧美精品在线观看 | 亚洲一区二区三区不卡视频| 丁香欧美五月| 男女视频在线观看网站免费 | 亚洲成人国产一区在线观看| 亚洲欧美日韩高清专用| 国内精品久久久久精免费| 在线观看www视频免费| 国产蜜桃级精品一区二区三区| 99国产精品一区二区蜜桃av| 在线播放国产精品三级| 999久久久精品免费观看国产| 国产私拍福利视频在线观看| 久久久国产成人免费| 怎么达到女性高潮| 波多野结衣高清无吗| 九九热线精品视视频播放| 成人av在线播放网站| 人妻丰满熟妇av一区二区三区| 动漫黄色视频在线观看| 色综合站精品国产| av片东京热男人的天堂| 亚洲专区字幕在线| 五月伊人婷婷丁香| 日韩大码丰满熟妇| 黄色丝袜av网址大全| 欧美精品亚洲一区二区| 亚洲精品色激情综合| 久久性视频一级片| 精品久久久久久久人妻蜜臀av| 亚洲专区中文字幕在线| 亚洲性夜色夜夜综合| 97超级碰碰碰精品色视频在线观看| 精品高清国产在线一区| 国产在线观看jvid| 18美女黄网站色大片免费观看| 午夜成年电影在线免费观看| 国产精品1区2区在线观看.| 久久精品aⅴ一区二区三区四区| cao死你这个sao货| 法律面前人人平等表现在哪些方面| 午夜激情av网站| 国产精品香港三级国产av潘金莲| 久久久久九九精品影院| 中文在线观看免费www的网站 | 久久久精品国产亚洲av高清涩受| 国产成人精品久久二区二区免费| 神马国产精品三级电影在线观看 | 变态另类成人亚洲欧美熟女| 精品国产美女av久久久久小说| 国模一区二区三区四区视频 | 国产99白浆流出| www.999成人在线观看| 后天国语完整版免费观看| 日韩大尺度精品在线看网址| 国产亚洲精品第一综合不卡| 日韩av在线大香蕉| 国产精品一区二区免费欧美| 久久精品人妻少妇| а√天堂www在线а√下载| 欧美日韩亚洲国产一区二区在线观看| 亚洲中文字幕日韩| 久久香蕉国产精品| 精品欧美一区二区三区在线| 精品不卡国产一区二区三区| 亚洲精品国产精品久久久不卡| 狂野欧美白嫩少妇大欣赏| 午夜日韩欧美国产| 久久国产乱子伦精品免费另类| 国产高清有码在线观看视频 | 亚洲国产精品999在线| 亚洲aⅴ乱码一区二区在线播放 | 国产精品精品国产色婷婷| 最新在线观看一区二区三区| 国产视频一区二区在线看| 日本精品一区二区三区蜜桃| 99久久精品国产亚洲精品| 熟妇人妻久久中文字幕3abv| 亚洲,欧美精品.| 欧美日韩一级在线毛片| 在线国产一区二区在线| 午夜a级毛片| 桃红色精品国产亚洲av| 午夜日韩欧美国产| 日韩欧美一区二区三区在线观看| 国产真实乱freesex| 99久久99久久久精品蜜桃| 欧美av亚洲av综合av国产av| 俺也久久电影网| 免费在线观看日本一区| 欧美av亚洲av综合av国产av| 给我免费播放毛片高清在线观看| 国产亚洲精品av在线| 丁香欧美五月| 真人一进一出gif抽搐免费| 免费一级毛片在线播放高清视频| 最好的美女福利视频网| 18美女黄网站色大片免费观看| 午夜精品一区二区三区免费看| 又粗又爽又猛毛片免费看| 黄频高清免费视频| 午夜精品一区二区三区免费看| 波多野结衣高清无吗| 黄色视频,在线免费观看| 18禁黄网站禁片午夜丰满| x7x7x7水蜜桃| 久久久国产成人免费| 两人在一起打扑克的视频| 免费观看精品视频网站| 国产精品99久久99久久久不卡| 91老司机精品| 变态另类成人亚洲欧美熟女| 中文字幕人成人乱码亚洲影| 亚洲人成77777在线视频| 久久九九热精品免费| 色尼玛亚洲综合影院| 久久久久国产精品人妻aⅴ院| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久久久久久毛片微露脸| 毛片女人毛片| 亚洲精品国产一区二区精华液| 国产精品乱码一区二三区的特点| 男女做爰动态图高潮gif福利片| 桃红色精品国产亚洲av| 亚洲自偷自拍图片 自拍| 久久99热这里只有精品18| 亚洲一卡2卡3卡4卡5卡精品中文| 激情在线观看视频在线高清| 亚洲无线在线观看| 国产精品永久免费网站| 天天躁夜夜躁狠狠躁躁| 亚洲av第一区精品v没综合| 亚洲天堂国产精品一区在线| 亚洲精品一区av在线观看| 日韩av在线大香蕉| 9191精品国产免费久久| 99热只有精品国产| 久久这里只有精品19| 欧美日韩瑟瑟在线播放| 日韩三级视频一区二区三区| 亚洲男人的天堂狠狠| 老司机靠b影院| av在线天堂中文字幕| 国产成人啪精品午夜网站| 成人av一区二区三区在线看| 色播亚洲综合网| 久久久精品大字幕| 18禁黄网站禁片免费观看直播| 嫁个100分男人电影在线观看| 免费一级毛片在线播放高清视频| 国产欧美日韩一区二区三| 狂野欧美白嫩少妇大欣赏| 色老头精品视频在线观看| 波多野结衣高清作品| 国产单亲对白刺激| 老熟妇仑乱视频hdxx| 51午夜福利影视在线观看| 91大片在线观看| 波多野结衣高清无吗| 老鸭窝网址在线观看| 免费在线观看日本一区| 国产精品久久久人人做人人爽| 99久久国产精品久久久| 久久香蕉精品热| 舔av片在线| 免费看日本二区| 欧美激情久久久久久爽电影| 精品久久久久久,| 成在线人永久免费视频| 日韩有码中文字幕| 国产精品一区二区免费欧美| av福利片在线观看| 精品国产美女av久久久久小说| 亚洲精华国产精华精| 成人18禁高潮啪啪吃奶动态图| 变态另类成人亚洲欧美熟女| 他把我摸到了高潮在线观看| 国产精品 国内视频| 午夜成年电影在线免费观看| 欧美日韩精品网址| 一区福利在线观看| 亚洲精品粉嫩美女一区| 18禁美女被吸乳视频| 亚洲成人久久爱视频| 亚洲国产中文字幕在线视频| 在线a可以看的网站| 久久精品国产亚洲av香蕉五月| 国产精品久久视频播放| 亚洲一区二区三区不卡视频| 亚洲中文av在线| √禁漫天堂资源中文www| 久久精品人妻少妇| 在线观看免费日韩欧美大片| 成人高潮视频无遮挡免费网站| 99在线人妻在线中文字幕| 老司机福利观看| 男人舔女人的私密视频| 99riav亚洲国产免费| 搡老熟女国产l中国老女人| 亚洲成人国产一区在线观看| 美女扒开内裤让男人捅视频| 成人三级黄色视频| 欧美色视频一区免费| 国产成人av激情在线播放| 在线观看www视频免费| 一二三四社区在线视频社区8| 香蕉国产在线看| 一级毛片精品| 日韩欧美国产一区二区入口| 亚洲av日韩精品久久久久久密| 桃色一区二区三区在线观看| 18禁黄网站禁片午夜丰满| 国产精品 欧美亚洲| 妹子高潮喷水视频| 国产精品日韩av在线免费观看| 久久久国产成人精品二区| 特级一级黄色大片| 美女免费视频网站| 久久天躁狠狠躁夜夜2o2o| 男女那种视频在线观看| 国产一区二区三区视频了| 国产欧美日韩精品亚洲av| 免费在线观看黄色视频的| 欧美高清成人免费视频www| 人人妻人人澡欧美一区二区| av视频在线观看入口| 欧美黄色淫秽网站| 成人永久免费在线观看视频| 午夜亚洲福利在线播放| av视频在线观看入口| 国产片内射在线| 国产午夜福利久久久久久| 久久中文字幕人妻熟女| 精品午夜福利视频在线观看一区| 国产精品乱码一区二三区的特点| 在线观看一区二区三区| 在线观看午夜福利视频| 免费在线观看日本一区| 波多野结衣高清无吗| 亚洲欧美一区二区三区黑人| 在线播放国产精品三级| 老司机靠b影院| 精品久久久久久久毛片微露脸| 女人被狂操c到高潮| 无限看片的www在线观看| 男人舔女人下体高潮全视频| 中国美女看黄片| 天堂影院成人在线观看| 成人av一区二区三区在线看| 久久久久久大精品| 精品日产1卡2卡| 成年版毛片免费区| 国产精品一及| 欧美日韩福利视频一区二区| 午夜日韩欧美国产| 啦啦啦免费观看视频1| 色精品久久人妻99蜜桃| 黄色a级毛片大全视频| 久久久久亚洲av毛片大全| 狠狠狠狠99中文字幕| 黄色丝袜av网址大全| 国产精品综合久久久久久久免费| 亚洲性夜色夜夜综合| 久久精品综合一区二区三区| 操出白浆在线播放| 精品久久久久久久久久久久久| 精品国产亚洲在线| 亚洲av日韩精品久久久久久密| www国产在线视频色| 在线十欧美十亚洲十日本专区| 日韩欧美一区二区三区在线观看| 狂野欧美白嫩少妇大欣赏| 欧美乱码精品一区二区三区| 美女午夜性视频免费| 日本一二三区视频观看| 天堂影院成人在线观看| 欧美又色又爽又黄视频| 大型av网站在线播放| 制服人妻中文乱码| 动漫黄色视频在线观看| 国产成年人精品一区二区| 无人区码免费观看不卡| 欧美不卡视频在线免费观看 | 一本一本综合久久| 中文字幕av在线有码专区| 巨乳人妻的诱惑在线观看| 这个男人来自地球电影免费观看| 一级毛片精品| 18禁黄网站禁片免费观看直播| 午夜免费激情av| 午夜日韩欧美国产| 国产欧美日韩一区二区三| 无人区码免费观看不卡| 成年版毛片免费区| 亚洲精品国产精品久久久不卡| 欧美日韩亚洲国产一区二区在线观看| 午夜免费成人在线视频| 久久久久性生活片| 99国产极品粉嫩在线观看| 欧美日韩福利视频一区二区| 精品电影一区二区在线| 国产精品美女特级片免费视频播放器 | 日本熟妇午夜| 欧美成人性av电影在线观看| 日韩精品中文字幕看吧| 欧美在线一区亚洲| 亚洲精品一卡2卡三卡4卡5卡| 日日干狠狠操夜夜爽| 一本大道久久a久久精品| 亚洲人成77777在线视频| 99国产精品99久久久久| 国产私拍福利视频在线观看| 可以免费在线观看a视频的电影网站| 亚洲中文字幕一区二区三区有码在线看 | 午夜视频精品福利| 亚洲欧美一区二区三区黑人| 日韩欧美国产在线观看| 波多野结衣高清无吗| a级毛片在线看网站| 成人一区二区视频在线观看| 精华霜和精华液先用哪个| 制服诱惑二区| 超碰成人久久| 99在线人妻在线中文字幕| 成熟少妇高潮喷水视频| 一区二区三区高清视频在线| 51午夜福利影视在线观看| 亚洲欧美日韩无卡精品| 91国产中文字幕| 国产伦一二天堂av在线观看| 99热只有精品国产| 熟女少妇亚洲综合色aaa.| 午夜视频精品福利| 可以在线观看的亚洲视频| 999久久久国产精品视频| 可以在线观看的亚洲视频| 999久久久国产精品视频| 18禁黄网站禁片午夜丰满| 嫩草影视91久久| 三级毛片av免费| 欧美+亚洲+日韩+国产| 中文字幕精品亚洲无线码一区| 18禁观看日本| 成年免费大片在线观看| 亚洲成人中文字幕在线播放| svipshipincom国产片| 国产精品久久久人人做人人爽| 在线播放国产精品三级| 男女床上黄色一级片免费看| 黄色毛片三级朝国网站| 久久久精品欧美日韩精品| 99久久99久久久精品蜜桃| 天天添夜夜摸| 亚洲国产精品合色在线| 亚洲精品久久国产高清桃花| 久久久久久免费高清国产稀缺| 香蕉国产在线看| 成人av一区二区三区在线看| 免费在线观看亚洲国产| 精华霜和精华液先用哪个| 欧美另类亚洲清纯唯美| 18禁国产床啪视频网站| 看免费av毛片| 天天添夜夜摸| 又黄又爽又免费观看的视频| 久久久久国内视频| 三级国产精品欧美在线观看 | 欧美最黄视频在线播放免费| or卡值多少钱| 成人高潮视频无遮挡免费网站| 性色av乱码一区二区三区2| 欧美精品啪啪一区二区三区| 久久精品成人免费网站| 麻豆久久精品国产亚洲av| 一进一出抽搐动态| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人啪精品午夜网站| 香蕉av资源在线| 91老司机精品| 90打野战视频偷拍视频| 亚洲一区二区三区不卡视频| 国产精品久久久久久人妻精品电影| 欧美日韩精品网址| 成年人黄色毛片网站| 一进一出抽搐gif免费好疼| 久久这里只有精品19| 久久国产乱子伦精品免费另类| 欧美成人免费av一区二区三区| 男女做爰动态图高潮gif福利片| 欧美又色又爽又黄视频| 夜夜夜夜夜久久久久| 在线视频色国产色| 亚洲国产中文字幕在线视频| 99热只有精品国产| 成人精品一区二区免费| 久久精品国产亚洲av高清一级| 手机成人av网站| 俺也久久电影网| 日韩欧美 国产精品| 黄色片一级片一级黄色片| 午夜福利18| 欧洲精品卡2卡3卡4卡5卡区| 波多野结衣高清作品| 国产激情偷乱视频一区二区| 精品免费久久久久久久清纯| 国产亚洲av高清不卡| 久久香蕉激情| 无遮挡黄片免费观看| 国产激情久久老熟女| 久久中文看片网| 国产精品野战在线观看| 91av网站免费观看| 亚洲精品一卡2卡三卡4卡5卡| 成人手机av| 嫩草影视91久久| 老司机午夜福利在线观看视频| 99久久精品热视频| 国产激情欧美一区二区| 黄频高清免费视频| 久久香蕉精品热| 免费在线观看黄色视频的| 国产亚洲精品第一综合不卡| 蜜桃久久精品国产亚洲av| 欧美黄色片欧美黄色片| 男男h啪啪无遮挡| 亚洲欧美一区二区三区黑人| 成人特级黄色片久久久久久久| 免费在线观看黄色视频的| 无遮挡黄片免费观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲av美国av| 99久久99久久久精品蜜桃| 欧美zozozo另类| av福利片在线观看| 一进一出抽搐动态| 欧美在线黄色| 午夜精品一区二区三区免费看| 露出奶头的视频| 亚洲最大成人中文| 久久九九热精品免费| 伊人久久大香线蕉亚洲五| 91麻豆av在线| 在线免费观看的www视频| 看黄色毛片网站| 久久久久久久精品吃奶| 久久久久久免费高清国产稀缺| x7x7x7水蜜桃| 黄色女人牲交| 看片在线看免费视频| 两个人视频免费观看高清| 一边摸一边做爽爽视频免费| 日本一区二区免费在线视频| 色噜噜av男人的天堂激情| 淫妇啪啪啪对白视频| 久久中文字幕人妻熟女| 欧美成人午夜精品| 热99re8久久精品国产| 巨乳人妻的诱惑在线观看| 欧美久久黑人一区二区| 丝袜美腿诱惑在线| 国产视频一区二区在线看| 午夜老司机福利片| 在线观看一区二区三区| 99精品在免费线老司机午夜| 51午夜福利影视在线观看| 成人特级黄色片久久久久久久| 精品一区二区三区四区五区乱码| 亚洲片人在线观看| 国产1区2区3区精品| 国产成人精品久久二区二区91| 国产99白浆流出| 久久久久久亚洲精品国产蜜桃av| 国产精品,欧美在线| 色综合婷婷激情| 亚洲五月天丁香| 日韩欧美 国产精品| 特大巨黑吊av在线直播| 哪里可以看免费的av片| 亚洲va日本ⅴa欧美va伊人久久| 男插女下体视频免费在线播放| 午夜免费成人在线视频| 天堂影院成人在线观看| 国产精品日韩av在线免费观看| 亚洲av片天天在线观看| av在线天堂中文字幕| 老司机深夜福利视频在线观看| 久久香蕉国产精品| tocl精华| 妹子高潮喷水视频| 亚洲人成电影免费在线| 国产黄片美女视频| 国产精品av久久久久免费| www.精华液| 99精品欧美一区二区三区四区| 国产成人av激情在线播放| АⅤ资源中文在线天堂| 成人一区二区视频在线观看| 久久久久久国产a免费观看| 精品人妻1区二区| 国产av在哪里看| 桃色一区二区三区在线观看| www.自偷自拍.com| 首页视频小说图片口味搜索| 亚洲国产看品久久| 99久久综合精品五月天人人| 日韩 欧美 亚洲 中文字幕| 校园春色视频在线观看| 国产一区在线观看成人免费| 99热这里只有精品一区 | 女同久久另类99精品国产91| 欧美色欧美亚洲另类二区| 免费av毛片视频| 亚洲国产欧美网| 俺也久久电影网| 在线视频色国产色| 国产99久久九九免费精品| 亚洲国产精品成人综合色| 国产爱豆传媒在线观看 | 亚洲精品国产精品久久久不卡| 白带黄色成豆腐渣| 国产成人一区二区三区免费视频网站| 婷婷六月久久综合丁香| 久久精品国产综合久久久| 一本久久中文字幕| 国产在线观看jvid| 99精品在免费线老司机午夜| 亚洲五月天丁香| 国产男靠女视频免费网站| 99热这里只有是精品50| 男女之事视频高清在线观看| 在线十欧美十亚洲十日本专区| 制服诱惑二区| avwww免费| 亚洲中文日韩欧美视频| 午夜福利在线观看吧| 国产av在哪里看| 亚洲av成人精品一区久久| 精品免费久久久久久久清纯| 午夜福利高清视频| 国产人伦9x9x在线观看| 久久中文字幕人妻熟女| 老司机靠b影院| 日韩欧美在线二视频| 最近视频中文字幕2019在线8| 午夜两性在线视频| 久久精品国产99精品国产亚洲性色| 亚洲熟妇熟女久久| 国产黄片美女视频| 成年免费大片在线观看| 天堂√8在线中文| 亚洲精品一卡2卡三卡4卡5卡| 亚洲成a人片在线一区二区| www.999成人在线观看| 国产乱人伦免费视频| 免费在线观看成人毛片| 五月玫瑰六月丁香| 777久久人妻少妇嫩草av网站| 国产视频内射| 国产高清视频在线观看网站| 婷婷亚洲欧美| 亚洲成a人片在线一区二区| 国产精品一区二区三区四区久久| 又黄又爽又免费观看的视频| 亚洲人成伊人成综合网2020| 天堂动漫精品| 国产精品免费视频内射| 久久久国产成人精品二区| 亚洲国产欧美人成| 中文亚洲av片在线观看爽| 精品久久蜜臀av无| 国产一区二区激情短视频| 亚洲人成网站在线播放欧美日韩| 91麻豆av在线| 亚洲精品美女久久av网站| 777久久人妻少妇嫩草av网站| 亚洲免费av在线视频| 国产av不卡久久| 一边摸一边做爽爽视频免费| 色综合欧美亚洲国产小说|