• <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新网站| 91精品伊人久久大香线蕉| 少妇人妻一区二区三区视频| 干丝袜人妻中文字幕| 久久久久精品性色| 欧美一级a爱片免费观看看| 熟妇人妻久久中文字幕3abv| 国产黄片美女视频| 男女边摸边吃奶| 纵有疾风起免费观看全集完整版 | 亚洲真实伦在线观看| 国产乱人视频| 99久国产av精品国产电影| 亚洲国产欧美在线一区| 最新中文字幕久久久久| 国产极品天堂在线| 肉色欧美久久久久久久蜜桃 | 久久鲁丝午夜福利片| h日本视频在线播放| 成年女人看的毛片在线观看| 国产乱人视频| av国产免费在线观看| 亚洲av成人精品一二三区| 亚洲精品乱久久久久久| 国国产精品蜜臀av免费| 嫩草影院新地址| 91精品伊人久久大香线蕉| 国内精品宾馆在线| 亚洲av不卡在线观看| 国产亚洲精品av在线| 男女啪啪激烈高潮av片| 高清av免费在线| 禁无遮挡网站| 精品国内亚洲2022精品成人| 久久精品人妻少妇| 国国产精品蜜臀av免费| 能在线免费看毛片的网站| 久久精品国产亚洲网站| 少妇熟女欧美另类| .国产精品久久| 亚洲国产高清在线一区二区三| 国产永久视频网站| 免费播放大片免费观看视频在线观看| 国产成人91sexporn| 亚洲自偷自拍三级| 我的女老师完整版在线观看| 免费观看a级毛片全部| 亚洲,欧美,日韩| 亚洲国产精品国产精品| 免费看不卡的av| 精品国产一区二区三区久久久樱花 | 亚洲在线自拍视频| 中国美白少妇内射xxxbb| 亚洲精品aⅴ在线观看| 亚洲怡红院男人天堂| 国产视频首页在线观看| 国语对白做爰xxxⅹ性视频网站| 欧美丝袜亚洲另类| 成人鲁丝片一二三区免费| 直男gayav资源| 日日啪夜夜爽| 一级毛片aaaaaa免费看小| 美女被艹到高潮喷水动态| 亚洲欧美日韩东京热| 寂寞人妻少妇视频99o| 免费观看a级毛片全部| 国国产精品蜜臀av免费| 久久午夜福利片| 熟女人妻精品中文字幕| 啦啦啦啦在线视频资源| 内地一区二区视频在线| 中文字幕人妻熟人妻熟丝袜美| kizo精华| h日本视频在线播放| 成人毛片60女人毛片免费| av.在线天堂| 伦理电影大哥的女人| 国产黄色免费在线视频| 国产不卡一卡二| 色吧在线观看| 日韩av免费高清视频| 麻豆国产97在线/欧美| 男女国产视频网站| 欧美成人a在线观看| 亚洲婷婷狠狠爱综合网| 国产av国产精品国产| 成人无遮挡网站| 在现免费观看毛片| av女优亚洲男人天堂| 免费黄色在线免费观看| 久久久久性生活片| 又爽又黄无遮挡网站| 欧美高清性xxxxhd video| 国产探花在线观看一区二区| 精品久久久噜噜| 欧美高清成人免费视频www| 国产一区二区三区综合在线观看 | av播播在线观看一区| 精品久久久久久久久亚洲| 亚洲av成人精品一二三区| 九九久久精品国产亚洲av麻豆| 狂野欧美激情性xxxx在线观看| 免费无遮挡裸体视频| 精品久久久久久久人妻蜜臀av| 肉色欧美久久久久久久蜜桃 | 亚洲成人精品中文字幕电影| av免费观看日本| 又粗又硬又长又爽又黄的视频| 女人十人毛片免费观看3o分钟| 大片免费播放器 马上看| 99久久精品国产国产毛片| 99久国产av精品国产电影| 久久精品国产亚洲网站| 99久久人妻综合| 日韩电影二区| 久久久久久久久久人人人人人人| 久久久亚洲精品成人影院| 亚洲国产精品专区欧美| 99久久中文字幕三级久久日本| 99九九线精品视频在线观看视频| 99久久精品国产国产毛片| 日本三级黄在线观看| 国产精品99久久久久久久久| 在线天堂最新版资源| 在线观看一区二区三区| 丰满乱子伦码专区| 六月丁香七月| 日韩一区二区三区影片| 国产伦精品一区二区三区四那| 波野结衣二区三区在线| 国产极品天堂在线| 一级毛片电影观看| 亚洲欧美清纯卡通| 久久久久久久国产电影| 日韩电影二区| 联通29元200g的流量卡| 国产单亲对白刺激| 久久亚洲国产成人精品v| 久久久亚洲精品成人影院| 日韩av免费高清视频| 麻豆成人午夜福利视频| 免费观看精品视频网站| 欧美三级亚洲精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日本一本二区三区精品| 成人高潮视频无遮挡免费网站| 国产av在哪里看| 亚洲av福利一区| 卡戴珊不雅视频在线播放| 三级毛片av免费| 一级爰片在线观看| 久久久久久久久久久免费av| 亚洲熟女精品中文字幕| 精品久久国产蜜桃| 国产精品嫩草影院av在线观看| 亚洲国产精品国产精品| 搞女人的毛片| 日韩av不卡免费在线播放| 国产国拍精品亚洲av在线观看| av国产久精品久网站免费入址| 极品少妇高潮喷水抽搐| 汤姆久久久久久久影院中文字幕 | 中文资源天堂在线| 日韩欧美三级三区| 亚洲成人久久爱视频| 内射极品少妇av片p| 色吧在线观看| 国产精品久久久久久久久免| 大话2 男鬼变身卡| 国产精品久久久久久av不卡| av专区在线播放| 国产一区二区三区av在线| 精品久久久精品久久久| 久久这里有精品视频免费| 国产精品无大码| 久久午夜福利片| 久久这里有精品视频免费| 日本午夜av视频| 一区二区三区乱码不卡18| 国产伦精品一区二区三区四那| 精华霜和精华液先用哪个| 亚洲va在线va天堂va国产| 久久这里有精品视频免费| 欧美bdsm另类| 国产精品久久久久久精品电影小说 | av女优亚洲男人天堂| 国产成人一区二区在线| 久久久久久久久久黄片| 亚洲精品456在线播放app| 精品酒店卫生间| 午夜激情久久久久久久| 国产精品久久久久久久久免| 亚洲精品乱码久久久久久按摩| 亚洲欧美日韩卡通动漫| 欧美变态另类bdsm刘玥| 精品国产三级普通话版| 国产老妇伦熟女老妇高清| 国产伦在线观看视频一区| 久久久久久伊人网av| 亚洲四区av| 一级爰片在线观看| 色尼玛亚洲综合影院| 国产美女午夜福利| 三级国产精品片| 亚洲av成人精品一区久久| 成人一区二区视频在线观看| 黄色配什么色好看| 蜜臀久久99精品久久宅男| 久久久国产一区二区| 美女cb高潮喷水在线观看| 精品久久久久久久久久久久久| 看免费成人av毛片| 3wmmmm亚洲av在线观看| 亚洲欧美日韩卡通动漫| 亚洲成人一二三区av| 久久精品国产亚洲av涩爱| 日韩精品青青久久久久久| 国产成人福利小说| 国产麻豆成人av免费视频| 国产伦精品一区二区三区四那| 2022亚洲国产成人精品| 狂野欧美激情性xxxx在线观看| 伊人久久国产一区二区| 亚洲自拍偷在线| 午夜精品在线福利| 亚洲丝袜综合中文字幕| 草草在线视频免费看| 久久久久网色| 欧美最新免费一区二区三区| 亚洲精品国产成人久久av| 伊人久久国产一区二区| 2021天堂中文幕一二区在线观| 毛片女人毛片| 亚洲av二区三区四区| 日本猛色少妇xxxxx猛交久久| 特大巨黑吊av在线直播| 日韩精品青青久久久久久| 国语对白做爰xxxⅹ性视频网站| 人妻一区二区av| 国产精品国产三级国产av玫瑰| 人妻系列 视频| 免费观看av网站的网址| 色播亚洲综合网| 99久久精品一区二区三区| 最新中文字幕久久久久| 嘟嘟电影网在线观看| 亚洲精品aⅴ在线观看| 嫩草影院入口| 日韩,欧美,国产一区二区三区| 久久久国产一区二区| 久久精品国产亚洲av天美| 自拍偷自拍亚洲精品老妇| 午夜精品在线福利| 熟女电影av网| videossex国产| 国产精品99久久久久久久久| 久久久久久国产a免费观看| 亚洲精品一区蜜桃| 又大又黄又爽视频免费| 超碰97精品在线观看| 成人毛片a级毛片在线播放| av国产久精品久网站免费入址| 91狼人影院| 久久人人爽人人爽人人片va| 午夜久久久久精精品| 日本午夜av视频| 狂野欧美白嫩少妇大欣赏| 亚洲电影在线观看av| av在线播放精品| 国产伦精品一区二区三区四那| 99热网站在线观看| 99久久人妻综合| 久久国产乱子免费精品| 国产淫片久久久久久久久| 成人亚洲欧美一区二区av| 99九九线精品视频在线观看视频| 看黄色毛片网站| 国产精品一区二区三区四区久久| av黄色大香蕉| 日韩精品有码人妻一区| 成人av在线播放网站| 最近中文字幕2019免费版| 综合色丁香网| 欧美日本视频| 男女国产视频网站| 80岁老熟妇乱子伦牲交| 夫妻午夜视频| 大又大粗又爽又黄少妇毛片口| 91av网一区二区| 国产亚洲av嫩草精品影院| 亚洲丝袜综合中文字幕| 两个人的视频大全免费| 精品久久久久久成人av| av网站免费在线观看视频 | 久久韩国三级中文字幕| 男女那种视频在线观看| 22中文网久久字幕| 久久久精品94久久精品| 建设人人有责人人尽责人人享有的 | 日产精品乱码卡一卡2卡三| 一本久久精品| 免费av观看视频| av免费在线看不卡| 嘟嘟电影网在线观看| 午夜福利高清视频| 国精品久久久久久国模美| 免费黄网站久久成人精品| 久久久久久久久久人人人人人人| 国产黄a三级三级三级人| 毛片女人毛片| 在线播放无遮挡| 久久精品夜夜夜夜夜久久蜜豆| 日本欧美国产在线视频| 中文字幕av在线有码专区| 国产精品无大码| 久久久久九九精品影院| 午夜视频国产福利| 国产一区二区在线观看日韩| 婷婷色综合www| 人人妻人人看人人澡| av国产免费在线观看| 免费观看的影片在线观看| 18+在线观看网站| 国内揄拍国产精品人妻在线| 国产精品嫩草影院av在线观看| 搞女人的毛片| 国产成年人精品一区二区| 建设人人有责人人尽责人人享有的 | 床上黄色一级片| 老女人水多毛片| 最新中文字幕久久久久| 亚洲人成网站在线播| 美女主播在线视频| 欧美成人精品欧美一级黄| 亚洲欧美精品专区久久| 欧美xxxx性猛交bbbb| 国产精品综合久久久久久久免费| 国产成人aa在线观看| 国产高潮美女av| 婷婷色综合www| 人体艺术视频欧美日本| 亚洲婷婷狠狠爱综合网| 欧美区成人在线视频| 久久韩国三级中文字幕| av福利片在线观看| 大片免费播放器 马上看| 国产精品.久久久| 一区二区三区四区激情视频| 黄色一级大片看看| 性插视频无遮挡在线免费观看| 日本熟妇午夜| 一区二区三区四区激情视频| 成年版毛片免费区| 丰满少妇做爰视频| 亚洲第一区二区三区不卡| 日韩欧美精品免费久久| 夜夜看夜夜爽夜夜摸| 国产在视频线精品| 大话2 男鬼变身卡| 久99久视频精品免费| av在线蜜桃| 国产老妇伦熟女老妇高清| 日本黄大片高清| 日本欧美国产在线视频| 男的添女的下面高潮视频| 国产人妻一区二区三区在| 91久久精品电影网| 国产中年淑女户外野战色| 22中文网久久字幕| 国产乱来视频区| 神马国产精品三级电影在线观看| 久久精品久久久久久久性| 久久久午夜欧美精品| 欧美xxxx性猛交bbbb| 男女边摸边吃奶| 免费黄色在线免费观看| 国产午夜精品久久久久久一区二区三区| 少妇熟女欧美另类| 亚洲熟女精品中文字幕| 内地一区二区视频在线| 一级毛片 在线播放| 国产在视频线在精品| 99久久中文字幕三级久久日本| 免费黄频网站在线观看国产| 少妇的逼好多水| 日韩欧美一区视频在线观看 | 国产精品一区二区三区四区免费观看| 亚洲aⅴ乱码一区二区在线播放| 99久久九九国产精品国产免费| 一区二区三区四区激情视频| 老司机影院毛片| 热99在线观看视频| 免费黄频网站在线观看国产| 亚洲国产精品专区欧美| 国产精品蜜桃在线观看| 国产av在哪里看| 久久精品久久久久久久性| 1000部很黄的大片| 嫩草影院入口| 日韩成人伦理影院| 91在线精品国自产拍蜜月| 国产精品.久久久| 欧美成人一区二区免费高清观看| 欧美日韩一区二区视频在线观看视频在线 | 超碰av人人做人人爽久久| 国产成人freesex在线| 国产高清有码在线观看视频| 看十八女毛片水多多多| videossex国产| 免费观看av网站的网址| 亚洲人与动物交配视频| 深夜a级毛片| 中文字幕制服av| 亚洲国产av新网站| 又黄又爽又刺激的免费视频.| 日韩欧美三级三区| 成人午夜精彩视频在线观看| 少妇裸体淫交视频免费看高清| 人妻制服诱惑在线中文字幕| 18禁裸乳无遮挡免费网站照片| 国产午夜福利久久久久久| 精品一区二区三区视频在线| 蜜臀久久99精品久久宅男| 成人二区视频| www.色视频.com| 男女那种视频在线观看| a级毛片免费高清观看在线播放| 亚洲无线观看免费| 亚洲欧美日韩卡通动漫| 欧美日韩一区二区视频在线观看视频在线 | 两个人视频免费观看高清| 黄色日韩在线| 日本欧美国产在线视频| 日日啪夜夜撸| 国产精品久久久久久久久免| 高清日韩中文字幕在线| 久久久精品免费免费高清| av一本久久久久| 久久久欧美国产精品| 欧美一区二区亚洲| 国产精品一区二区在线观看99 | 岛国毛片在线播放| 午夜老司机福利剧场| 人妻夜夜爽99麻豆av| 99热6这里只有精品| 亚洲精品视频女| 51国产日韩欧美| 在线天堂最新版资源| 亚洲精品国产成人久久av| 日韩强制内射视频| 久久精品国产鲁丝片午夜精品| 婷婷色麻豆天堂久久| 男人和女人高潮做爰伦理| 九九在线视频观看精品| 亚洲综合精品二区| 亚洲国产日韩欧美精品在线观看| 三级国产精品欧美在线观看| 国产一区亚洲一区在线观看| 日本三级黄在线观看| 欧美日韩精品成人综合77777| 成年免费大片在线观看| 男女啪啪激烈高潮av片| 国产片特级美女逼逼视频| 午夜福利在线观看免费完整高清在| 性色avwww在线观看| 一个人看的www免费观看视频| 一区二区三区四区激情视频| 熟妇人妻久久中文字幕3abv| 欧美性感艳星| 欧美高清成人免费视频www| 狂野欧美激情性xxxx在线观看| 一级毛片aaaaaa免费看小| 国产男人的电影天堂91| 高清视频免费观看一区二区 | 国产毛片a区久久久久| 高清欧美精品videossex| 亚洲精品自拍成人| 熟妇人妻不卡中文字幕| av又黄又爽大尺度在线免费看| 免费看不卡的av| 免费观看性生交大片5| 你懂的网址亚洲精品在线观看| 国产精品一区www在线观看| 国产三级在线视频| 欧美成人午夜免费资源| 亚洲欧洲日产国产| 岛国毛片在线播放| 男人爽女人下面视频在线观看| 久久99精品国语久久久| 亚洲精品乱码久久久v下载方式| 亚洲精品色激情综合| 91精品一卡2卡3卡4卡| 亚洲精品影视一区二区三区av| 亚洲精品aⅴ在线观看| 成人二区视频| 亚洲av免费在线观看| 小蜜桃在线观看免费完整版高清| 免费观看性生交大片5| 成年女人看的毛片在线观看| a级毛色黄片| 菩萨蛮人人尽说江南好唐韦庄| 国产欧美日韩精品一区二区| 免费观看无遮挡的男女| 欧美xxⅹ黑人| 国产精品伦人一区二区| 国产一区亚洲一区在线观看| 欧美+日韩+精品| 亚洲欧美清纯卡通| 国产黄频视频在线观看| 观看免费一级毛片| 美女xxoo啪啪120秒动态图| 国产老妇女一区| 最近2019中文字幕mv第一页| 亚洲av成人精品一二三区| 欧美潮喷喷水| 国产免费视频播放在线视频 | 特大巨黑吊av在线直播| 午夜激情久久久久久久| 黄色配什么色好看| 成人av在线播放网站| 婷婷色综合大香蕉| 欧美一级a爱片免费观看看| 18禁在线播放成人免费| av.在线天堂| 亚洲一区高清亚洲精品| 亚州av有码| 97超碰精品成人国产| 国模一区二区三区四区视频| 男女边摸边吃奶| 直男gayav资源| 国国产精品蜜臀av免费| 白带黄色成豆腐渣| 中文乱码字字幕精品一区二区三区 | 久久久久网色| 国产一区二区在线观看日韩| 日韩一区二区视频免费看| 日日干狠狠操夜夜爽| 亚洲欧美一区二区三区黑人 | 久久久国产一区二区| 久久鲁丝午夜福利片| 老司机影院成人| eeuss影院久久| 免费人成在线观看视频色| 亚洲aⅴ乱码一区二区在线播放| 91久久精品电影网| 高清毛片免费看| 男女边吃奶边做爰视频| 免费无遮挡裸体视频| 97精品久久久久久久久久精品| 在线观看一区二区三区| 国产午夜精品一二区理论片| 最后的刺客免费高清国语| 日本与韩国留学比较| 中文在线观看免费www的网站| or卡值多少钱| 亚洲电影在线观看av| 久久草成人影院| 天堂√8在线中文| 看免费成人av毛片| 日韩 亚洲 欧美在线| 热99在线观看视频| 国产精品不卡视频一区二区| 禁无遮挡网站| 精品久久久久久久久亚洲| 中文字幕久久专区| 午夜福利视频精品| 蜜臀久久99精品久久宅男| av在线老鸭窝| 国产黄色小视频在线观看| 国产免费一级a男人的天堂| 特大巨黑吊av在线直播| 久久久精品94久久精品| 九九在线视频观看精品| 国产一级毛片在线| 国产精品三级大全| 久久久久久久久久人人人人人人| 国产69精品久久久久777片| 在线观看一区二区三区| 国产高清三级在线| a级毛片免费高清观看在线播放| 99久久精品一区二区三区| 中文天堂在线官网| 国产中年淑女户外野战色| 三级国产精品欧美在线观看| 国产在视频线精品| 国产成人精品久久久久久| 久久久久精品性色| 亚洲精品,欧美精品| 超碰97精品在线观看| 久久这里只有精品中国| 亚洲在线自拍视频| 一级二级三级毛片免费看| 精品酒店卫生间| 亚洲欧美中文字幕日韩二区| 国产一区二区三区av在线| 亚洲美女搞黄在线观看| 亚洲自偷自拍三级| 成人欧美大片| 久久6这里有精品| 国内精品美女久久久久久| 一级黄片播放器| 最新中文字幕久久久久| 久久精品熟女亚洲av麻豆精品 | 国产黄色视频一区二区在线观看| 五月伊人婷婷丁香| 老师上课跳d突然被开到最大视频| 美女黄网站色视频| 午夜免费男女啪啪视频观看| 成人午夜高清在线视频| 亚洲成人av在线免费| 国产亚洲av嫩草精品影院| 亚洲av中文av极速乱| 亚洲人成网站在线播|