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

    常州市春季PM2.5中細(xì)菌群落特征及影響因素

    2018-09-25 00:57:12薛銀剛屠博文王利平趙亞芳江曉棟
    中國環(huán)境科學(xué) 2018年9期
    關(guān)鍵詞:顆粒物群落測序

    劉 菲,薛銀剛,*,屠博文,王利平*,趙亞芳,江曉棟,薛 柯

    常州市春季PM2.5中細(xì)菌群落特征及影響因素

    劉 菲1,薛銀剛1,2*,屠博文3,王利平1*,趙亞芳2,江曉棟1,薛 柯1

    (1.常州大學(xué)環(huán)境與安全工程學(xué)院,江蘇 常州 213164;2.江蘇省環(huán)境保護(hù)水環(huán)境生物監(jiān)測重點(diǎn)實(shí)驗(yàn)室,常州市環(huán)境監(jiān)測中心,江蘇 常州 213001;3.常州市疾病預(yù)防控制中心,江蘇 常州 213022)

    為研究常州市春季PM2.5中細(xì)菌群落特征,利用高通量測序?qū)M2.5中細(xì)菌的16S rRNA基因進(jìn)行研究.測序獲得的有效序列數(shù)為600150,以97%的相似水平劃分,各樣品的OTUs為1890~6519,同時(shí)樣本的Coverage指數(shù)較高,表明測序結(jié)果可以準(zhǔn)確地代表樣本中的空氣細(xì)菌群落.物種注釋結(jié)果表明:常州春季PM2.5中相對豐度>1%的有11個(gè)細(xì)菌門、14個(gè)細(xì)菌綱和12個(gè)細(xì)菌屬,其中變形菌門(Proteobacteria)、藍(lán)藻細(xì)菌門(Cyanobacteria)和厚壁菌門(Firmicutes)是排名前3的優(yōu)勢細(xì)菌類群,占總基因豐度的80.88%.屬水平優(yōu)勢細(xì)菌主要有擬甲色球藻屬(,6.03%)、(5.95%)、微囊藻屬(,4.86%)和鞘氨醇單胞菌屬(,3.16%),但在屬水平上未能進(jìn)行分類的基因序列比例高達(dá)81.11%.基于PM2.5中細(xì)菌群落組成進(jìn)行來源分析,發(fā)現(xiàn)常州市春季PM2.5中細(xì)菌的環(huán)境來源多變,其主要環(huán)境源可能為淡水,其次是土壤、植物和人為源.利用冗余分析探討環(huán)境因子與細(xì)菌群落的關(guān)系,結(jié)果表明,NH4+、NO3-、O3、SO42-、OC、氣壓和CO是常州市春季PM2.5中細(xì)菌群落的主要影響因子,同時(shí)不同環(huán)境因子對不同細(xì)菌類群影響不同.

    春季;PM2.5;細(xì)菌群落;來源分析;環(huán)境因子

    近年來,顆粒物已造成我國多數(shù)城市空氣環(huán)境污染嚴(yán)重,而細(xì)顆粒物(PM2.5)是近年來大中城市灰霾天氣頻發(fā)的主導(dǎo)因素[1-2].PM2.5由于粒徑小,在大氣中停留時(shí)間長,傳輸距離遠(yuǎn),并且由于其比表面積大,成為許多有毒有害物質(zhì)的載體,對人體健康產(chǎn)生危害[3-4].研究表明空氣中存在一定數(shù)量且具有代謝活性的微生物,細(xì)菌在近地面大氣中普遍存在且為大氣微生物中非常重要的一部分[5-6].目前關(guān)于PM2.5的研究大多集中于化學(xué)成分、來源分析及其對人體健康的影響[7-10],針對PM2.5中細(xì)菌群落特征及其致病性的研究相對較少,同時(shí)影響PM2.5中細(xì)菌群落結(jié)構(gòu)的因素在很大程度上未知.

    空氣中細(xì)菌的研究方法主要有傳統(tǒng)培養(yǎng)法和16S測序分析方法[5].而高通量測序技術(shù)是對傳統(tǒng)測序技術(shù)一次革命性的改變,能夠準(zhǔn)確全面地反映細(xì)菌群落結(jié)構(gòu)和多樣性特征,已被廣泛應(yīng)用于不同生境的細(xì)菌群落及其功能研究[11-12].高通量測序技術(shù)亦是分析大氣中細(xì)菌的一種有效技術(shù),已被用于研究霧霾過程中生物氣溶膠的細(xì)菌群落結(jié)構(gòu)[13]和分析空氣中細(xì)菌群落來源[14].

    已有研究表明四季內(nèi)微生物氣溶膠數(shù)量以春季較高[15],所以本次研究基于Illumina HiSeq測序平臺,利用16S rDNA測序方法分析常州市春季PM2.5中細(xì)菌群落結(jié)構(gòu),并基于細(xì)菌群落特征進(jìn)行來源解析.利用冗余分析(RDA)研究氣象因子、氣態(tài)污染物和顆粒物中化學(xué)組分對PM2.5中細(xì)菌群落結(jié)構(gòu)的影響,這對全面了解PM2.5中細(xì)菌的環(huán)境效應(yīng)具有重要意義.

    1 材料與方法

    1.1 樣品采集

    采樣地點(diǎn)位于常州市環(huán)境監(jiān)測中心東樓樓頂(119°7′47″E,31°45′44″N),是大氣顆粒物自動(dòng)監(jiān)測國控站點(diǎn),為居民文教區(qū),采樣高度距離地面10m左右.采樣時(shí)間為2017年3月28日~2017年4月20日,采樣期間溫度逐漸回升,采樣持續(xù)時(shí)間為23h (15:00~次日14:00)[16].

    采樣儀器為1臺嶗應(yīng)2031型智能空氣采樣器(青島嶗山應(yīng)用技術(shù)研究所),此采樣器為大流量采樣器,配PM2.5切割頭,流量為1050L/min.采樣所用濾膜為玻璃纖維濾膜(20.3cm×25.4cm),濾膜在采樣前均使用馬弗爐450℃煅燒2h,相關(guān)實(shí)驗(yàn)操作參照文獻(xiàn)[10],待采樣完成后將濾膜保存于-20℃,用于后續(xù)分析.

    1.2 氣象因子、氣態(tài)污染物及化學(xué)組分監(jiān)測

    采樣時(shí)記錄氣溫、氣壓、風(fēng)速和相對濕度等氣象條件.SO2、O3、CO和NO2質(zhì)量濃度由Thermo Scientific? 146i多種氣體校準(zhǔn)儀(Thermo Fisher Scientific, USA)實(shí)時(shí)監(jiān)測;水溶性離子(NO3-、SO42-、NH4+、Cl-、Na+、K+、Mg2+和Ca2+)質(zhì)量濃度由瑞士萬通中國有限公司的MARGA ADI 2080離子在線監(jiān)測儀器實(shí)時(shí)監(jiān)測;PM2.5和PM10質(zhì)量濃度由TEOM 1405系列顆粒物監(jiān)測儀(Thermo Fisher Scientific, USA)同時(shí)監(jiān)測.元素碳(EC)和有機(jī)碳(OC)質(zhì)量濃度由Sunset RT-4(Sunset Lab, USA)實(shí)時(shí)觀測.

    1.3 DNA提取

    用滅菌的剪刀將樣品膜剪成小的膜片(約為樣品膜的1/12~1/8),再用牙簽刷(滅菌)刮取樣品膜上的顆粒物,并盡量減少濾膜的摻入,接著用剪刀將刮取的樣品盡量剪碎后置于PowerBead Tube中,后續(xù)按PowerSoil DNA Isolation Kit(Mobio, USA)試劑盒操作標(biāo)準(zhǔn)方法提取顆粒物樣品DNA.用NanoDrop2000超微量蛋白質(zhì)核酸分析儀(Thermo Fisher, USA)測定提取出的DNA濃度和純度,提取成功的樣品置于-20℃中保存,用于后續(xù)的PCR擴(kuò)增.

    1.4 PCR擴(kuò)增和熒光定量

    將提取好的DNA產(chǎn)物使用16S V4區(qū)引物進(jìn)行PCR擴(kuò)增,V4區(qū)特異性引物序列正向引物: AYTGGGYDTAAAGNG,反向引物:TACNVGGGT- ATCTAATCC[17].PCR擴(kuò)增見表1,總體系共計(jì)50μL. PCR擴(kuò)增程序:95℃預(yù)變性3min,95℃變性30s,55℃退火溫度下反應(yīng)30s,72℃延伸30s,共30個(gè)循環(huán),最后在72℃延伸5min.將樣本的PCR產(chǎn)物用2%瓊脂糖凝膠電泳檢測,使用AMPure XP Beads(Beckman Coulter, USA)進(jìn)行產(chǎn)物純化.將PCR純化產(chǎn)物用Qubit?dsDNA HS Assay Kits熒光定量系統(tǒng)(Life technology, USA)進(jìn)行定量檢測.

    表1 PCR擴(kuò)增體系

    1.5 高通量測序

    使用TruSeq?DNA PCR-Free Sample Preparation Kit(Illumina, USA)建庫試劑盒進(jìn)行文庫構(gòu)建,構(gòu)建好的文庫使用Qubit @ 2.0熒光計(jì)(Thermo Scientific, USA)和Agilent Bioanalyzer 2100系統(tǒng)(Agilent, USA)評估測序文庫質(zhì)量,最后在Illumina HiSeq 2500平臺(Illumina, USA)上進(jìn)行測序.

    1.6 生物信息分析

    Hiseq測序得到的配對末端讀數(shù)首先根據(jù)重疊關(guān)系進(jìn)行拼接,同時(shí)對序列質(zhì)量進(jìn)行質(zhì)控和過濾,使用Uparse(v7.0.1001)軟件進(jìn)行序列分析,基于有效數(shù)據(jù)劃分操作分類單元(OUT),按照97%相似性對非重復(fù)序列進(jìn)行OTU聚類,在聚類過程中去除嵌合體,得到OTU的代表序列.篩選每個(gè)OTU的代表序列使用基于RDP分類器算法的GreenGene數(shù)據(jù)庫來注釋分類信息[18].為了研究不同OTU的系統(tǒng)發(fā)育關(guān)系,以及不同樣本中優(yōu)勢種的差異,使用Muscle軟件進(jìn)行多序列比對(3.8.31版本).通過Qiime軟件(Version 1.7.0)計(jì)算細(xì)菌群落Alpha多樣性指數(shù)[19].

    1.7 數(shù)據(jù)處理

    試驗(yàn)數(shù)據(jù)經(jīng)處理后使用Origin 2017繪制圖形.環(huán)境因子與細(xì)菌群落的相關(guān)性研究首先經(jīng)除趨勢對應(yīng)分析排序(DCA)分析,并根據(jù)分析結(jié)果中第一軸的梯度長度,選定冗余分析提取顯著解釋細(xì)菌群落的環(huán)境因子[20],并研究環(huán)境因子與細(xì)菌群落的相關(guān)性.

    2 結(jié)果與分析

    2.1 采樣期間環(huán)境因子監(jiān)測結(jié)果

    表2 采樣期間氣象數(shù)據(jù)

    采樣期間氣象數(shù)據(jù)見表2.采樣期間溫度為11.6~27℃;氣壓變化范圍為1002~1021hPa;采樣期間風(fēng)速為0.8~2.08m/s,軟風(fēng)及輕風(fēng)天氣偏多;相對濕度平均值為55%.

    各樣品采樣時(shí)段大氣污染物演變狀況見圖1,除O3為8h滑動(dòng)平均值外,其余污染物均為24h平均值.各污染物濃度變化范圍分別為SO2(6~47μg/m3)、NO2(18~55μg/m3)、PM10(78~251μg/m3)、O3(74~ 167μg/m3)、PM2.5(36~75μg/m3)和CO(0.627~ 2.066mg/m3).

    圖1 采樣期間氣態(tài)污染物變化特征

    圖2 采樣期間水溶性離子變化特征

    各樣品采集期間水溶性離子變化特征見圖2,可以看出本次樣品采集過程中二次顆粒物所占比重較大[21],其中NO3-、SO42-和NH4+占比重較高,NO3-變化范圍為6.324~27.870μg/m3,SO42-變化范圍為5.351~21.462μg/m3,NH4+變化范圍為3.865~ 16.570μg/m3,3種離子平均百分比之和高達(dá)97.61%.而Cl-、Na+、K+、Mg2+和Ca2+的相對百分比較小,Cl-變化范圍為1.186~2.921μg/m3,Na+變化范圍為0.150~0.770μg/m3、K+變化范圍為0.285~1.265μg/ m3、Mg2+變化范圍為0.012~0.156μg/m3,Ca2+變化范圍為0.037~0.609μg/m3.

    2.2 細(xì)菌群落多樣性與結(jié)構(gòu)特征

    經(jīng)高通量測序,PM2.5樣品在門、綱、目、科、屬和種分類水平上獲得注釋信息的基因序列分別為600150、600150、593852、577204、558179和113381.在97%分類水平下,各樣品的OTUs為1890~6519,反映各采樣點(diǎn)樣品OTU覆蓋率的Coverage指數(shù)均達(dá)到94%以上,可以看出測序效果理想(表3).

    表3 各樣品OTU聚類與多樣性分析

    通過描述群落豐富度、多樣性和均勻度的Alpha多樣性指數(shù)來解釋PM2.5細(xì)菌群落物種組成情況(表3).各樣品細(xì)菌群落豐富度和多樣性均有所差異,E樣品菌群結(jié)構(gòu)的多樣性和均勻度較高,但細(xì)菌群落的豐富度較低,F樣品菌群結(jié)構(gòu)的多樣性和均勻度較低,但細(xì)菌群落的豐富度較高,同時(shí)分析結(jié)果表明環(huán)境因子與Alpha多樣性指數(shù)間無顯著相關(guān)關(guān)系(>0.05),可能的影響因素還需進(jìn)一步研究.總體來說,Chao1指數(shù)是OTU數(shù)目的2倍左右,表明本次關(guān)于細(xì)菌豐富度的測序深度約達(dá)到了多樣性估計(jì)值的2倍[22].

    物種注釋結(jié)果表明顆粒物樣品共檢出20個(gè)細(xì)菌門,圖3為門分類水平下PM2.5的菌群結(jié)構(gòu).相對豐度最高的是變形菌門(Proteobacteria,63.87%).4類變形菌綱為PM2.5中的優(yōu)勢細(xì)菌綱,其中占比較高的是α-變形菌綱(Alphaproteobacteria,23.20%)和β-變形菌綱(Betaproteobacteria,13.37%),γ-變形菌綱(Gammaproteobacteria,7.57%)和δ-變形菌綱(Deltaproteobacteria,3.72%)平均相對豐度較低.藍(lán)細(xì)菌門(Cyanobacteria,11.21%)和厚壁菌門(Firmicutes,5.80%)是平均豐度排名第2和第3的優(yōu)勢類群;其中Cyanobacteria和Firmicutes分別在D樣品和A樣品相對豐度最高,分別為19.82%和9.31%.此外,相對豐度較高的細(xì)菌門依次為擬桿菌門(Bacteroidetes,3.27%)、放線菌門(Actinobacteria, 2.90%)、異常球菌-棲熱菌門(Deinococcus-Thermus, 2.05%)、疣微菌門(Verrucomicrobia,1.75%)和浮霉菌門(Planctomycetes,1.51%)和酸桿菌門(Acidobacteria, 1.12%),將占比較低的細(xì)菌類群歸為Others類.

    圖3 PM2.5中門分類水平細(xì)菌群落組成

    表4 優(yōu)勢細(xì)菌屬豐度與來源解析

    此外,平均相對豐度>1%的優(yōu)勢屬共有12個(gè),分別為擬甲色球藻屬(,6.03%)、(5.95%)、微囊藻屬(, 4.86%)、鞘氨醇單胞菌屬(,3.16%)、奇異球菌屬(,2.43%)、薄層桿菌屬(,1.69%)、甲基桿菌屬(,1.65%)、馬賽菌屬(, 1.61%)、副球菌屬(,1.52%)、(1.46%)、假單胞菌屬(,1.33%)和不動(dòng)細(xì)菌屬(,1.30%);優(yōu)勢細(xì)菌屬豐度變化與來源環(huán)境見表4.

    2.3 細(xì)菌群落與環(huán)境因子相關(guān)性

    圖4展示了氣象因子、氣態(tài)污染物和顆粒物化學(xué)組分對PM2.5中細(xì)菌群落的冗余分析結(jié)果,可以看出環(huán)境因子對PM2.5中細(xì)菌群落具有一定影響.經(jīng)冗余分析提取顯著解釋細(xì)菌群落的環(huán)境因子.通過篩選,NH4+、NO3-、O3、SO42-、OC、氣壓(QY)和CO對細(xì)菌群落的解釋量較高,分別為0.442、0.362、0.307、0.280、0.263、0.256和0.208.冗余分析結(jié)果表明第一主軸對門水平細(xì)菌群落方差變化的解釋量為76.9%,第二主軸的解釋量為17.6%,到軸4的累計(jì)合理解釋量為98.9%.

    圖4 環(huán)境因子與細(xì)菌群落冗余分析

    圖4中各箭頭連線夾角的余弦值表示二者的相關(guān)關(guān)系,可用于分析環(huán)境因子與不同物種間、環(huán)境因子間以及不同物種間的相關(guān)關(guān)系.其中,NH4+、NO3-與Cyanobacteria、Proteobacteria呈顯著正相關(guān)(<0.01);SO2、氣溫(QW)與Bacteroidetes、Firmicutes呈顯著負(fù)相關(guān)(<0.05);CO與Actinobacteria、Planctomycetes和Verrucomicrobia呈顯著負(fù)相關(guān)(<0.05),與Chloroflexi、Cyanobacteria和Proteobacteria呈顯著正相關(guān)(<0.05);OC、EC、氣壓與Bacteroidetes、Firmicutes呈顯著正相關(guān)(<0.05);SO42-與Firmicutes呈顯著正相關(guān)(<0.05).

    3 討論

    3.1 PM2.5中細(xì)菌群落結(jié)構(gòu)及來源分析

    大氣顆粒物中的微生物是大氣環(huán)境的重要組成部分,而細(xì)菌在PM2.5已分析出的生物組分中占比高達(dá)80%以上[21];經(jīng)高通量測序,常州春季PM2.5中細(xì)菌群落主要由11個(gè)細(xì)菌門、14個(gè)細(xì)菌綱和12個(gè)細(xì)菌屬組成(相對豐度>1%).Fang等[30]利用FA-1型撞擊式空氣微生物采樣器采集可培養(yǎng)細(xì)菌,結(jié)果表明杭州空氣中可培養(yǎng)細(xì)菌主要為微球藻屬()、芽胞桿菌屬()、葡萄球菌屬()、考克氏菌屬()和假單胞菌屬(),占總細(xì)菌種數(shù)的58%;而姚文沖[31]利用高通量測序技術(shù)研究杭州市空氣細(xì)菌群落特征,發(fā)現(xiàn)、、和是相對含量較高的菌屬;而、、和在本次研究中均有檢出,相對豐度分別為0.9%、0.3%、0.2%和0.1%,占比較低.由于空氣中只有大約1%的微生物能夠在培養(yǎng)基上形成菌落,因此運(yùn)用傳統(tǒng)培養(yǎng)法對群落結(jié)構(gòu)的判定誤差較大[5].相對來說,高通量測序技術(shù)可以更加細(xì)致全貌的反映研究樣品的菌群結(jié)構(gòu).本次PM2.5中大部分細(xì)菌是不致病的,但部分潛在的空氣傳播的病原體也被檢出,如約氏不動(dòng)桿菌(,0.16%)、魯菲不動(dòng)桿菌(,1.30%)[32]和盲腸腸球菌(,0.12%)[33].

    來源于不同環(huán)境的細(xì)菌附著在PM2.5上,以氣溶膠形式存在于城市大氣環(huán)境中,是評價(jià)環(huán)境質(zhì)量的重要指標(biāo).細(xì)菌在自然界的物質(zhì)循環(huán)中扮演著重要角色,自然界的外在條件不僅影響著細(xì)菌氣溶膠在環(huán)境空氣中的比重,而且還影響著細(xì)菌群落的組成及多樣性[15].PM2.5中細(xì)菌來源廣泛,且可在水體、土壤和生物體之間進(jìn)行交換并組成一個(gè)龐大而復(fù)雜的體系.廖旭[2]等應(yīng)用T-RFLP、克隆文庫等測序方法研究廈門冬季PM2.5中細(xì)菌和真核微型生物的群落組成,并分析其潛在的來源環(huán)境,結(jié)果表明廈門城區(qū)空氣微生物的重要環(huán)境源可能為淡水,其次是土壤、水體沉積物、污水系統(tǒng)和動(dòng)物糞便等.此次常州市春季PM2.5中細(xì)菌主要為Proteobacteria,這主要是由于Proteobacteria來源廣泛且對大氣環(huán)境有較強(qiáng)的適應(yīng)能力,在土壤、淡水水體與大氣等多種環(huán)境介質(zhì)中廣泛分布[34-36],同時(shí)對低溫和輻射具有較強(qiáng)的抵抗能力[37].除Proteobacteria外,Cyanobacteria、Firmicutes和Bacteroidetes主要是淡水環(huán)境中的優(yōu)勢細(xì)菌[38-39],同時(shí)也有研究表明淡水環(huán)境是青島市區(qū)街道四季和人工濕地夏季空氣細(xì)菌的主要來源[37],表明水體很有可能是PM2.5中細(xì)菌的最主要來源,這里的水體不僅是淡水來源,還包括天然的雨水、人為的降水以及各種污水等[40].土壤作為細(xì)菌巨大的繁殖場所和貯存體,風(fēng)會(huì)把大量細(xì)菌送入空氣中形成細(xì)菌氣溶膠,所以土壤也是PM2.5中細(xì)菌的主要發(fā)生源[22],同時(shí)、和均為土壤環(huán)境中的優(yōu)勢細(xì)菌(表4).PM2.5中還檢出了已知和植物相關(guān)的細(xì)菌類群,如、和(表4),表明植物源也是常州春季PM2.5中細(xì)菌的來源之一.經(jīng)比較,常州春季PM2.5中細(xì)菌的主要來源環(huán)境為淡水,其次為土壤和植物,同時(shí)在樣品采集過程中二次顆粒物所占的比重較大,因此人為源可能也是PM2.5中細(xì)菌的主要來源之一.

    3.2 PM2.5中細(xì)菌群落結(jié)構(gòu)影響因子研究

    細(xì)菌在大氣中存在時(shí),會(huì)受到各種環(huán)境因子的影響,比如氣候、天氣變化、當(dāng)?shù)氐纳飦碓春蜕硥m暴等人為和自然的污染源[41].空氣中細(xì)菌的新陳代謝和生長繁殖首先會(huì)受到氣溫的影響;存在于大氣中的氣態(tài)污染物(SO2、CO和O3等)會(huì)參加一系列的物理化學(xué)反應(yīng),從而引發(fā)大氣污染,會(huì)對細(xì)菌生長和存活產(chǎn)生影響;同時(shí)氣溶膠顆粒物中存在的化學(xué)成分不僅為細(xì)菌提供粘附的媒介,也可為細(xì)菌的生長和存活提供適宜的環(huán)境[42].可見,氣象因子、氣態(tài)污染物和氣溶膠顆粒物中化學(xué)組分共同影響著大氣中細(xì)菌群落結(jié)構(gòu).

    由于大氣顆粒物中存在著復(fù)雜多樣的組分,如有機(jī)物、元素、離子以及微生物,經(jīng)冗余分析(圖4)發(fā)現(xiàn)影響常州市春季PM2.5中細(xì)菌群落結(jié)構(gòu)的主要環(huán)境因子為NH4+、NO3-、O3、SO42-、OC、氣壓和CO.王步英[21]發(fā)現(xiàn)氣象因子、氣態(tài)污染物和顆粒物中化學(xué)組分都會(huì)對北京冬季大氣中細(xì)菌群落結(jié)構(gòu)產(chǎn)生影響,其中溫度、相對濕度、O3、NO2、顆粒物中SO42-和K+是影響細(xì)菌群落結(jié)構(gòu)特征的重要環(huán)境因子.Innocente等[43]利用冗余分析研究米蘭和威尼斯兩個(gè)城市空氣顆粒物中水溶性離子對細(xì)菌群落結(jié)構(gòu)的影響,結(jié)果表明,NH4+和Ca2+與細(xì)菌群落具有很強(qiáng)的相關(guān)性.可見,不同地區(qū)、不同季節(jié)的環(huán)境影響因子對細(xì)菌群落結(jié)構(gòu)產(chǎn)生不同的影響.同時(shí)也有研究[44]表明風(fēng)速較大時(shí)會(huì)引起邊界層的擾動(dòng),使顆粒物懸浮,從而增加大氣中細(xì)菌氣溶膠的濃度.經(jīng)冗余分析,O3與PM2.5中大多數(shù)細(xì)菌呈負(fù)相關(guān)關(guān)系(圖4),這主要是由于O3毒性和破壞性很強(qiáng),對生物有機(jī)體危害極大.PM2.5中優(yōu)勢細(xì)菌與不同環(huán)境因子的關(guān)系均有所差異,這不僅與細(xì)菌自身特性有關(guān),同時(shí)與環(huán)境因子的相互作用、細(xì)菌間的相關(guān)作用以及環(huán)境因子與不同細(xì)菌間的耦合關(guān)系都相關(guān).而目前關(guān)于這方面的研究還比較缺乏,今后可開展相關(guān)研究以彌補(bǔ)該方面研究的空缺.

    4 結(jié)論

    4.1 常州春季PM2.5中細(xì)菌群落多樣性和豐富度較高,相對豐度>1%的有11個(gè)細(xì)菌門、14個(gè)細(xì)菌綱和12個(gè)細(xì)菌屬,其中Proteobacteria、Cyanobacteria和Firmicutes是排名前3的優(yōu)勢細(xì)菌類群,占總基因豐度的80.88%.屬水平優(yōu)勢細(xì)菌主要有(6.03%)、(5.95%)、(4.86%)和(3.16%).

    4.2 常州市春季PM2.5中細(xì)菌的環(huán)境來源多變,其主要環(huán)境源可能為淡水,其次是土壤、植物和人為源.

    4.3 NH4+、NO3-、O3、SO42-、OC、氣壓和CO是常州市春季PM2.5中細(xì)菌群落結(jié)構(gòu)的主要影響因子,同時(shí)不同環(huán)境因子對不同細(xì)菌類群影響有所差異.

    [1] Xu Q, Li X, Wang S, et al. Fine particulate air pollution and hospital emergency room visits for respiratory disease in urban areas in Beijing, China, in 2013 [J]. Plos One, 2016,11(4):e0153099.

    [2] 廖 旭,胡安誼,楊曉永,等.廈門冬季PM2.5顆粒物中細(xì)菌和真核微型生物群落組成及其來源分析 [J]. 生態(tài)環(huán)境學(xué)報(bào), 2013,22(8): 1395-1400.

    [3] 靳璐濱,臧增亮,潘曉濱,等.PM2.5和PM2.5~10資料同化及在南京青奧會(huì)期間的應(yīng)用試驗(yàn) [J]. 中國環(huán)境科學(xué), 2016,36(2):331-341.

    [4] Ye Z L, Li Q, Liu J S, et al. Investigation of submicron aerosol characteristics in Changzhou, China: Composition, source, and comparison with co-collected PM2.5[J]. Chemosphere, 2017,183:176- 185.

    [5] 王 琳,宋志文,徐愛玲,等.青島市秋季空氣微生物群落多樣性 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 2015,26(4):1121-1129.

    [6] 王步英,郎繼東,張麗娜,等.基于16S rRNA基因測序法分析北京霾污染過程中PM2.5和PM10細(xì)菌群落特征 [J]. 環(huán)境科學(xué), 2015, 36(8):2727-2734.

    [7] Tan J H, Zhang L M, Zhou X M, et al. Chemical characteristics and source apportionment of PM2.5in Lanzhou, China [J]. Science of the Total Environment, 2017,601-602:1743-1752.

    [8] Rogula-Koz?owska W. Size-segregated urban particulate matter: mass closure, chemical composition, and primary and secondary matter content [J]. Air Quality Atmosphere & Health, 2016,9(5):533-550.

    [9] Cao S J, Kong X R, Li L Y, et al. An investigation of the PM2.5and NO2concentrations and their human health impacts in the metro subway system of Suzhou, China [J]. Environmental Science: Processes & Impacts, 2017,19(5):666-675.

    [10] Bari M A, Kindzierski W B. Ambient fine particulate matter (PM2.5) in Canadian oil sands communities: Levels, sources and potential human health risk [J]. Science of the Total Environment, 2017,595:828-838.

    [11] Sun Y J, Wang T Y, Peng X W, et al. Bacterial community compositions in sediment polluted by perfluoroalkyl acids (PFAAs) using Illumina high-throughput sequencing [J]. Environmental Science & Pollution Research, 2016,23(11):10556-10565.

    [12] Zhang W, Chen L, Zhang R, et al. High throughput sequencing analysis of the joint effects of BDE209-Pb on soil bacterial community structure [J]. Journal of Hazardous Materials, 2016,301: 1-7.

    [13] 馬曼曼,甄 毓,米鐵柱,等.青島冬季霾天不同粒徑生物氣溶膠中細(xì)菌群落特征研究 [J]. 中國環(huán)境科學(xué), 2017,37(8):2855-2865.

    [14] Adams R I, Miletto M, Lindow S E, et al. Airborne bacterial communities in residences: similarities and differences with fungi [J]. Plos One, 2014,9(3):e91283.

    [15] 黃 瓊.南昌市PM10和PM2.5中微生物群落的分布及其相關(guān)性分析 [D]. 南昌:南昌大學(xué), 2015.

    [16] Cao C, Jiang W J, Wang B Y, et al. Inhalable microorganisms in Beijing’s PM2.5and PM10pollutants during a severe smog event [J]. Environmental Science & Technology, 2014,48(3):1499-1507.

    [17] Zhang J Y, Ding X, Guan R, et al. Evaluation of different 16S rRNA gene V regions for exploring bacterial diversity in a eutrophic freshwater lake [J]. Science of the Total Environment, 2017,618:1254- 1267.

    [18] Desantis T Z, Hugenholtz P, Larsen N, et al. Greengenes: chimera- checked 16S rRNA gene database and workbenchcompatible in ARB [J]. Applied & Environmental Microbiology, 2006,72(7):5069-5072.

    [19] Caporaso J G, Kuczynski J, Stombaugh J, et al. QIIME allows analysis of high-throughput community sequencing data [J]. Nature Methods, 2010,7(5):335-336.

    [20] Ter-Braak C J F, ?milauer P. CANOCO reference manual and CanoDraw for Windows user's guide: software for canonical community ordination (version 4.5) [R]. Ithaca, New York: Microcomputer Power, 2002.

    [21] 王步英.北京冬季大氣中細(xì)菌群落結(jié)構(gòu)特征及霾污染對其影響[D]. 北京:清華大學(xué), 2015.

    [22] 茍歡歌,謝純斌,童延斌,等.石河子市春季大氣顆粒物TSP和PM10中微生物群落特征 [J]. 環(huán)境工程學(xué)報(bào), 2017,11(4):2343-2349.

    [23] 劉 洋,朱夢靈,徐 瑤,等.中國藍(lán)藻植物的新記錄屬—擬甲色球藻[J]. 水生生物學(xué)報(bào), 2013,37(3):413-417.

    [24] Figuerola E L M, Guerrero L D, Rosa S M, et al. Bacterial indicator of agricultural management for soil under no-till crop production [J]. Plos One, 2012,7(11):e51075.

    [25] Balkwill D L, Fredrickson J K, Romine M F.and related genera [M]. New York: Springer, 2006:605-629.

    [26] 陳 立,邢瑞云,董俊興.屬細(xì)菌分類和應(yīng)用的研究進(jìn)展 [J]. 微生物學(xué)雜志, 2008,28(5):73-76.

    [27] Lidstrom M E, Chistoserdova L. Plants in the pink: cytokinin production by[J]. Journal of Bacteriology, 2002, 184(7):1818.

    [28] Fahlgren C, Bratbak G, Sandaa R A, et al. Diversity of airborne bacteria in samples collected using different devices for aerosol collection [J]. Aerobiologia, 2011,27(2):107-120.

    [29] Gao C Y, Wang A J, Wu W M, et al. Enrichment of anodic biofilm inoculated with anaerobic or aerobic sludge in single chambered air- cathode microbial fuel cells [J]. Bioresource Technology, 2014,167(3): 124-132.

    [30] Fang Z G, Yao W C, Lou X Q, et al. Profile and characteristics of culturable airborne bacteria in Hangzhou, southeast of China [J]. Aerosol & Air Quality Research, 2016,16(7):1690-1700.

    [31] 姚文沖.基于高通量測序的南方典型旅游城市空氣細(xì)菌群落特征研究 [D]. 杭州:浙江工商大學(xué), 2016:72-74.

    [32] Kozińska A, Pa?dzior E, P?kala A, et al.and- the emerging fish pathogens [J]. Bulletin of the Veterinary Institute in Pulawy, 2014,58(2):193-199.

    [33] Borst L B, Suyemoto M M, Keelara S, et al. A chicken embryo lethality assay for pathogenic[J]. Avian Diseases, 2014,58(2):244-248.

    [34] Zeng Q C, An S S, Liu Y. Soil bacterial community response to vegetation succession after fencing in the grassland of China [J]. Science of the Total Environment, 2017,609:2-10.

    [35] Fan L M, Song C, Meng S L, et al. Spatial distribution of planktonic bacterial and archaeal communities in the upper section of the tidal reach in Yangtze River [J]. Scientific Reports, 2016,6:39147.

    [36] Maki T, Hara K, Kobayashi F, et al. Vertical distribution of airborne bacterial communities in an Asian-dust downwind area, Noto Peninsula [J]. Atmospheric Environment, 2015,119:282-293.

    [37] 趙亞冰.青島市市區(qū)街道和人工濕地空氣細(xì)菌群落結(jié)構(gòu)研究 [D]. 青島:青島理工大學(xué), 2015.

    [38] Su X M, Steinman A D, Xue Q J, et al. Temporal patterns of phyto- and bacterioplankton and their relationships with environmental factors in Lake Taihu, China [J]. Chemosphere, 2017,184:299-308.

    [39] Li J F, Zhang J Y, Liu L Y, et al. Annual periodicity in planktonic bacterial and archaeal community composition of eutrophic Lake Taihu [J]. Scientific Reports, 2015,5:15488.

    [40] Huffman J A, Prenni A J, Demott P J, et al. High concentrations of biological aerosol particles and ice nuclei during and after rain [J]. Atmospheric Chemistry & Physics, 2013,13(13):6151-6164.

    [41] Soleimani Z, Parhizgari N, Rad H D, et al. Normal and dusty days comparison of culturable indoor airborne bacteria in Ahvaz, Iran [J]. Aerobiologia, 2015,31(2):127-141.

    [42] Han Y P, Li L, Liu J X, et al. Microbial structure and chemical components of aerosols caused by rotating brushes in a wastewater treatment plant [J]. Environmental Science & Pollution Research, 2012,19(9):4097-4108.

    [43] Innocente E, Squizzato S, Visin F, et al. Influence of seasonality, air mass origin and particulate matter chemical composition on airborne bacterial community structure in the Po Valley, Italy [J]. Science of the Total Environment, 2017,s593-594:677-687.

    [44] 翟晴飛,金蓮姬,林振毅,等.石家莊春季大氣氣溶膠數(shù)濃度和譜的觀測特征 [J]. 中國環(huán)境科學(xué), 2011,31(6):886-891.

    Characteristics and influence factors of airborne bacterial communities in Changzhou’s PM2.5samples during spring.

    LIU Fei1, XUE Yin-gang1,2*, TU Bo-wen3, WANG Li-ping1*, ZHAO Ya-fang2, JIANG Xiao-dong1, XUE Ke1

    (1.School of Environmental and Safety Engineering, Changzhou University, Changzhou 213164, China;2.Key Laboratory of Environmental Protection of Water Environment Biological Monitoring of Jiangsu Province, Changzhou Environmental Monitoring Center, Changzhou 213001, China;3.Changzhou Centers for Disease Control and Prevention, Changzhou 213022, China)., 2018,38(9):3254~3261

    In order to investigate the characteristics of the airborne bacterial community in Changzhou’s PM2.5samplesduring spring, the 16S rRNA gene of bacterial in PM2.5sampleswas studied by high-throughput sequencing. 600150sequences obtained could be divided into 1890~6519 OTUs of each sample at 97% similarity level. The Coverage index of the samples was high and could accurately represent the airborne bacterial community in the samples. Species annotation results indicated that there were 11bacterial phyla, 14bacterial classes, 12bacterial genera with relative abundance greater than 1% in Changzhou’s PM2.5samples during spring. Proteobacteria, Bacteroidetes and Firmicutes were the top 3dominant phyla, accounting for 80.88% of the total gene abundance. At the genus level,(6.03%),(5.95%),(4.86%) and(3.16%) were predominant, but the proportion of gene sequences that failed to be classified was up to 81.11%. Based on the source analysis of bacterial community composition in PM2.5, it was found that the environmental sources of bacteria in Changzhou’s PM2.5samples during spring were varied, and the main environmental source might be fresh water, followed by soil, plants and anthropogenic sources. The results of redundancy analysis (RDA) for the relationship between environmental factors and bacterial community showed that NH4+, NO3-, O3, SO42-, OC, air pressure and CO were the main environmental factors affecting the bacterial community in Changzhou’s PM2.5samplesduring spring. At the same time, varied environmental factors had different effects on bacterial groups.

    spring;PM2.5;bacterial community;source analysis;environmental factors

    X172

    A

    1000-6923(2018)09-3254-08

    劉 菲(1994-),女,安徽宣城人,常州大學(xué)碩士研究生,主要從事氣溶膠環(huán)境效應(yīng)及抗性基因污染特征研究.發(fā)表論文8篇.

    2018-02-05

    國家自然科學(xué)基金青年科學(xué)基金項(xiàng)目(21607016);常州市科技局科技支撐(社會(huì)發(fā)展)項(xiàng)目(CE20175022);環(huán)境基準(zhǔn)與風(fēng)險(xiǎn)評估國家重點(diǎn)實(shí)驗(yàn)室開放課題(SKLECRA20160FP20);江蘇省自然科學(xué)基金青年基金項(xiàng)目(BK20150250);江蘇省預(yù)防醫(yī)學(xué)會(huì)科研課題項(xiàng)目(Y2015016)

    * 責(zé)任作者, 薛銀剛, 高級工程師, yzxyg@126.com; 王利平, 教授, wlp@cczu.edu.cn

    猜你喜歡
    顆粒物群落測序
    杰 Sir 帶你認(rèn)識宏基因二代測序(mNGS)
    新民周刊(2022年27期)2022-08-01 07:04:49
    大學(xué)生牙齦炎齦上菌斑的微生物群落
    合成微生物群落在發(fā)酵食品中的應(yīng)用研究
    二代測序協(xié)助診斷AIDS合并馬爾尼菲籃狀菌腦膜炎1例
    傳染病信息(2021年6期)2021-02-12 01:52:58
    南平市細(xì)顆粒物潛在來源分析
    錯(cuò)流旋轉(zhuǎn)填料床脫除細(xì)顆粒物研究
    多層介質(zhì)阻擋放電處理柴油機(jī)尾氣顆粒物
    春季和夏季巢湖浮游生物群落組成及其動(dòng)態(tài)分析
    基因捕獲測序診斷血癌
    單細(xì)胞測序技術(shù)研究進(jìn)展
    麻豆成人午夜福利视频| 国产精品一区二区免费欧美| 亚洲天堂国产精品一区在线| 日本一本二区三区精品| 欧美潮喷喷水| 亚洲精华国产精华液的使用体验 | 色哟哟·www| 日韩欧美免费精品| 国产亚洲欧美98| 真人做人爱边吃奶动态| 身体一侧抽搐| 国产精品人妻久久久影院| 欧美成人一区二区免费高清观看| 亚洲久久久久久中文字幕| 国产精品99久久久久久久久| 神马国产精品三级电影在线观看| 一级a爱片免费观看的视频| 精品久久国产蜜桃| 狂野欧美白嫩少妇大欣赏| 联通29元200g的流量卡| 成人毛片a级毛片在线播放| 色综合亚洲欧美另类图片| 精品免费久久久久久久清纯| 内地一区二区视频在线| 欧美日韩综合久久久久久 | 很黄的视频免费| 国产欧美日韩精品一区二区| 神马国产精品三级电影在线观看| 我要搜黄色片| 亚洲avbb在线观看| 三级男女做爰猛烈吃奶摸视频| 国产不卡一卡二| 搡老妇女老女人老熟妇| 一卡2卡三卡四卡精品乱码亚洲| 伦精品一区二区三区| 人人妻,人人澡人人爽秒播| 琪琪午夜伦伦电影理论片6080| 国产淫片久久久久久久久| 最近最新中文字幕大全电影3| 亚洲精品一卡2卡三卡4卡5卡| 最新中文字幕久久久久| 毛片一级片免费看久久久久 | 久久欧美精品欧美久久欧美| 日韩欧美免费精品| 日本三级黄在线观看| 白带黄色成豆腐渣| 美女黄网站色视频| 亚洲av免费高清在线观看| 日本一本二区三区精品| 日韩一本色道免费dvd| 精品久久久久久久久av| 精品无人区乱码1区二区| 成人精品一区二区免费| 国产精品福利在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 99久国产av精品| 不卡视频在线观看欧美| 99久久精品热视频| 麻豆国产av国片精品| 亚洲在线自拍视频| 此物有八面人人有两片| 日韩欧美免费精品| 尤物成人国产欧美一区二区三区| 亚洲第一区二区三区不卡| 精品乱码久久久久久99久播| 国产美女午夜福利| 国产亚洲精品av在线| 人人妻,人人澡人人爽秒播| 舔av片在线| 99热这里只有是精品在线观看| 国产亚洲精品av在线| 日本与韩国留学比较| 午夜福利成人在线免费观看| 亚洲,欧美,日韩| 亚洲欧美日韩无卡精品| 国产精品久久久久久av不卡| 天堂影院成人在线观看| 成人av一区二区三区在线看| 国产精品自产拍在线观看55亚洲| 国产v大片淫在线免费观看| 人妻制服诱惑在线中文字幕| www.色视频.com| 午夜福利在线观看吧| 国产精品永久免费网站| 一区二区三区高清视频在线| 国产成人福利小说| 成人鲁丝片一二三区免费| 搞女人的毛片| eeuss影院久久| 性插视频无遮挡在线免费观看| 高清在线国产一区| 在线国产一区二区在线| av天堂在线播放| 波多野结衣高清作品| 国产黄色小视频在线观看| 特大巨黑吊av在线直播| 亚洲人成网站在线播放欧美日韩| 亚洲国产欧美人成| 成年免费大片在线观看| 亚洲色图av天堂| 久久久久久伊人网av| 久久久久九九精品影院| 看黄色毛片网站| 自拍偷自拍亚洲精品老妇| 亚洲av一区综合| 人妻久久中文字幕网| 亚洲av二区三区四区| 高清日韩中文字幕在线| 51国产日韩欧美| 搡老熟女国产l中国老女人| 成人三级黄色视频| a在线观看视频网站| 午夜精品一区二区三区免费看| 亚洲成人久久性| 综合色av麻豆| 色av中文字幕| 国产午夜福利久久久久久| 亚洲国产欧美人成| 精品99又大又爽又粗少妇毛片 | 久99久视频精品免费| 赤兔流量卡办理| 久久人人精品亚洲av| а√天堂www在线а√下载| 免费黄网站久久成人精品| 变态另类成人亚洲欧美熟女| 亚洲国产欧美人成| 99热这里只有精品一区| 全区人妻精品视频| 听说在线观看完整版免费高清| 久久久国产成人精品二区| 少妇高潮的动态图| 哪里可以看免费的av片| 亚洲av不卡在线观看| 欧美区成人在线视频| 九色成人免费人妻av| xxxwww97欧美| 全区人妻精品视频| 我要搜黄色片| 天堂网av新在线| 神马国产精品三级电影在线观看| 国产精品久久视频播放| a级一级毛片免费在线观看| 精品人妻1区二区| 亚洲欧美日韩卡通动漫| 亚洲av成人精品一区久久| 日本-黄色视频高清免费观看| 色综合亚洲欧美另类图片| 久久久久久久精品吃奶| 两人在一起打扑克的视频| 精品久久久久久久久久免费视频| 熟妇人妻久久中文字幕3abv| 国产极品精品免费视频能看的| 欧美bdsm另类| 亚洲av一区综合| 国产欧美日韩一区二区精品| 日日摸夜夜添夜夜添av毛片 | 日本一本二区三区精品| 精品人妻视频免费看| 夜夜爽天天搞| 亚洲精华国产精华液的使用体验 | av在线观看视频网站免费| 亚洲欧美日韩卡通动漫| 成人精品一区二区免费| 男女下面进入的视频免费午夜| 日本一二三区视频观看| 色5月婷婷丁香| 婷婷色综合大香蕉| 亚洲av免费高清在线观看| 国产亚洲精品av在线| 国产精品av视频在线免费观看| 毛片女人毛片| 国产伦一二天堂av在线观看| 亚洲熟妇中文字幕五十中出| 欧美另类亚洲清纯唯美| 91久久精品电影网| 免费观看人在逋| av黄色大香蕉| 久久国产乱子免费精品| 麻豆成人午夜福利视频| 免费av不卡在线播放| 少妇的逼好多水| 亚洲一区高清亚洲精品| 在线天堂最新版资源| 欧美性猛交黑人性爽| 国产精品久久电影中文字幕| 国产精品久久久久久亚洲av鲁大| 午夜免费激情av| 日韩欧美在线乱码| 丰满人妻一区二区三区视频av| 精品无人区乱码1区二区| h日本视频在线播放| 天堂√8在线中文| 亚洲精品456在线播放app | 亚洲av一区综合| 51国产日韩欧美| 国产精品野战在线观看| 老熟妇仑乱视频hdxx| 久久久久久久久大av| 欧美在线一区亚洲| 两人在一起打扑克的视频| av专区在线播放| 国内精品美女久久久久久| 香蕉av资源在线| 男女视频在线观看网站免费| 国产v大片淫在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 国产精品无大码| 免费看美女性在线毛片视频| 99热6这里只有精品| 日韩在线高清观看一区二区三区 | 欧美日韩中文字幕国产精品一区二区三区| 久久久色成人| 久久久久精品国产欧美久久久| 在线观看av片永久免费下载| 中国美女看黄片| 可以在线观看毛片的网站| 美女大奶头视频| 亚洲欧美日韩高清专用| av天堂中文字幕网| 日本免费a在线| 精品一区二区三区av网在线观看| 成人三级黄色视频| 人妻少妇偷人精品九色| 变态另类成人亚洲欧美熟女| 狠狠狠狠99中文字幕| 人妻丰满熟妇av一区二区三区| 97人妻精品一区二区三区麻豆| 毛片女人毛片| 一个人看的www免费观看视频| 亚洲国产欧美人成| 精品一区二区三区av网在线观看| 黄色女人牲交| 极品教师在线视频| 神马国产精品三级电影在线观看| 午夜老司机福利剧场| 日韩欧美在线乱码| 长腿黑丝高跟| 久久国产精品人妻蜜桃| 日韩中文字幕欧美一区二区| 无遮挡黄片免费观看| 国产极品精品免费视频能看的| 少妇被粗大猛烈的视频| 少妇裸体淫交视频免费看高清| 夜夜夜夜夜久久久久| 国产精品不卡视频一区二区| 午夜免费成人在线视频| 男女边吃奶边做爰视频| 可以在线观看的亚洲视频| 国产激情偷乱视频一区二区| 国产爱豆传媒在线观看| 亚洲国产欧美人成| 亚洲成a人片在线一区二区| 日韩在线高清观看一区二区三区 | eeuss影院久久| 99精品久久久久人妻精品| 黄色一级大片看看| 一级毛片久久久久久久久女| 日韩在线高清观看一区二区三区 | 成年女人毛片免费观看观看9| 欧美激情久久久久久爽电影| 18禁黄网站禁片免费观看直播| 两个人视频免费观看高清| 尤物成人国产欧美一区二区三区| 麻豆成人av在线观看| 国产91精品成人一区二区三区| 欧美成人一区二区免费高清观看| 亚洲图色成人| 日本成人三级电影网站| 国产91精品成人一区二区三区| 日本欧美国产在线视频| 天堂av国产一区二区熟女人妻| 国产三级在线视频| 超碰av人人做人人爽久久| 久久草成人影院| 欧美不卡视频在线免费观看| 高清日韩中文字幕在线| 一本一本综合久久| 日日摸夜夜添夜夜添小说| 婷婷精品国产亚洲av在线| 岛国在线免费视频观看| 国产精品一区二区免费欧美| 欧美又色又爽又黄视频| 成人性生交大片免费视频hd| 国产熟女欧美一区二区| 久久久久国内视频| 欧美激情久久久久久爽电影| 日本a在线网址| 日韩一区二区视频免费看| 丰满的人妻完整版| 女的被弄到高潮叫床怎么办 | 人妻久久中文字幕网| 人人妻,人人澡人人爽秒播| 久久久久九九精品影院| 国产精品亚洲一级av第二区| 久久天躁狠狠躁夜夜2o2o| 国产亚洲精品综合一区在线观看| 五月玫瑰六月丁香| 欧美国产日韩亚洲一区| 亚洲美女视频黄频| 蜜桃久久精品国产亚洲av| 成人三级黄色视频| 中出人妻视频一区二区| 久久人人精品亚洲av| 精品欧美国产一区二区三| 国产一区二区亚洲精品在线观看| 欧美一区二区精品小视频在线| 丰满乱子伦码专区| 成人二区视频| 国产精品亚洲一级av第二区| 亚洲内射少妇av| 午夜精品在线福利| 国产精品人妻久久久久久| av中文乱码字幕在线| 欧美在线一区亚洲| 国产精品一区二区免费欧美| 97碰自拍视频| 欧美成人性av电影在线观看| 国产 一区 欧美 日韩| 成人毛片a级毛片在线播放| 一级毛片久久久久久久久女| 国产av一区在线观看免费| 久久久久久九九精品二区国产| 精品免费久久久久久久清纯| 乱码一卡2卡4卡精品| 久久久久久九九精品二区国产| 嫩草影院新地址| 老女人水多毛片| 日韩一本色道免费dvd| 一本一本综合久久| 国产女主播在线喷水免费视频网站 | 欧美日韩中文字幕国产精品一区二区三区| 亚洲av中文字字幕乱码综合| 女生性感内裤真人,穿戴方法视频| 久久久精品大字幕| 国产免费av片在线观看野外av| 久久精品国产99精品国产亚洲性色| 午夜免费成人在线视频| 欧美丝袜亚洲另类 | 在线免费十八禁| 在线天堂最新版资源| 一区二区三区免费毛片| 全区人妻精品视频| 在现免费观看毛片| 国内精品美女久久久久久| 99久久精品国产国产毛片| 日本成人三级电影网站| 久久久久国内视频| 精品一区二区三区视频在线观看免费| 国产亚洲精品久久久com| 在线天堂最新版资源| 18禁裸乳无遮挡免费网站照片| 三级男女做爰猛烈吃奶摸视频| 性插视频无遮挡在线免费观看| 亚洲av五月六月丁香网| 亚洲五月天丁香| 国语自产精品视频在线第100页| 69av精品久久久久久| 亚洲国产精品合色在线| 亚洲五月天丁香| 国产精品1区2区在线观看.| 免费搜索国产男女视频| 99久久九九国产精品国产免费| 欧美激情在线99| av在线蜜桃| 日本欧美国产在线视频| 一进一出抽搐gif免费好疼| 午夜精品久久久久久毛片777| 悠悠久久av| 久久久国产成人免费| 身体一侧抽搐| 18禁黄网站禁片午夜丰满| 久久精品人妻少妇| 欧美丝袜亚洲另类 | 中文在线观看免费www的网站| 香蕉av资源在线| 午夜久久久久精精品| 久久久久久久久大av| 午夜久久久久精精品| 美女黄网站色视频| 欧美最黄视频在线播放免费| 久久久久久大精品| 国产精品99久久久久久久久| 久久6这里有精品| 婷婷亚洲欧美| 18禁裸乳无遮挡免费网站照片| 亚洲中文字幕日韩| 简卡轻食公司| av在线观看视频网站免费| 桃色一区二区三区在线观看| 亚洲精品日韩av片在线观看| 婷婷丁香在线五月| 色综合亚洲欧美另类图片| 日韩欧美国产在线观看| 国产av不卡久久| 国产av麻豆久久久久久久| 成年女人永久免费观看视频| 中文资源天堂在线| 欧美激情在线99| 国产精品98久久久久久宅男小说| 日韩欧美精品免费久久| av在线蜜桃| 亚洲人与动物交配视频| 国产午夜精品久久久久久一区二区三区 | 亚洲精品久久国产高清桃花| 国产一区二区三区视频了| 麻豆成人午夜福利视频| 精品人妻1区二区| 久久久国产成人精品二区| 3wmmmm亚洲av在线观看| 99久久成人亚洲精品观看| 99热只有精品国产| 免费一级毛片在线播放高清视频| 国产精品日韩av在线免费观看| 午夜日韩欧美国产| 欧美成人性av电影在线观看| 一级黄色大片毛片| 毛片女人毛片| 欧美又色又爽又黄视频| 亚洲av不卡在线观看| 久久婷婷人人爽人人干人人爱| 国产又黄又爽又无遮挡在线| 国产一区二区激情短视频| 在线国产一区二区在线| 两个人视频免费观看高清| 亚洲经典国产精华液单| 午夜精品在线福利| 人人妻,人人澡人人爽秒播| 免费看日本二区| 一区福利在线观看| 国产色爽女视频免费观看| 九色成人免费人妻av| 亚洲精华国产精华液的使用体验 | 露出奶头的视频| 搞女人的毛片| 日韩欧美精品免费久久| 午夜福利高清视频| 亚洲精品456在线播放app | 国产伦精品一区二区三区四那| 亚洲最大成人av| 伊人久久精品亚洲午夜| 最后的刺客免费高清国语| 国产免费一级a男人的天堂| 国内精品久久久久久久电影| 精华霜和精华液先用哪个| 老司机福利观看| 99热网站在线观看| 中文字幕久久专区| 我要看日韩黄色一级片| av专区在线播放| 波多野结衣高清作品| 久久久久久久亚洲中文字幕| 在线播放国产精品三级| 国内少妇人妻偷人精品xxx网站| 毛片女人毛片| 一个人看的www免费观看视频| 噜噜噜噜噜久久久久久91| 联通29元200g的流量卡| 亚洲一级一片aⅴ在线观看| 久久精品国产鲁丝片午夜精品 | 亚洲黑人精品在线| 中国美白少妇内射xxxbb| 国产亚洲91精品色在线| 最近最新中文字幕大全电影3| 日韩精品有码人妻一区| 极品教师在线视频| 免费看a级黄色片| 狂野欧美白嫩少妇大欣赏| 色哟哟哟哟哟哟| 亚洲四区av| 久久香蕉精品热| 少妇猛男粗大的猛烈进出视频 | 嫩草影院精品99| 在线国产一区二区在线| 丰满乱子伦码专区| 可以在线观看毛片的网站| 亚洲熟妇中文字幕五十中出| 精品久久久久久久久久免费视频| 精品乱码久久久久久99久播| 中文字幕av在线有码专区| av视频在线观看入口| av.在线天堂| 免费在线观看成人毛片| 久久国产乱子免费精品| 十八禁国产超污无遮挡网站| 我要看日韩黄色一级片| 波野结衣二区三区在线| 能在线免费观看的黄片| 少妇被粗大猛烈的视频| 中文字幕久久专区| 一夜夜www| 中文亚洲av片在线观看爽| 午夜激情欧美在线| 久久亚洲精品不卡| 国产精品电影一区二区三区| 日本a在线网址| 亚洲第一区二区三区不卡| 麻豆精品久久久久久蜜桃| 久久精品国产亚洲av香蕉五月| 欧美潮喷喷水| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩黄片免| 桃色一区二区三区在线观看| 高清日韩中文字幕在线| 国产精品不卡视频一区二区| 亚洲精品久久国产高清桃花| 此物有八面人人有两片| 国内久久婷婷六月综合欲色啪| 国产成人福利小说| 婷婷亚洲欧美| 久99久视频精品免费| 久久精品国产自在天天线| 国产成人a区在线观看| 欧美一区二区亚洲| 亚洲国产高清在线一区二区三| aaaaa片日本免费| 他把我摸到了高潮在线观看| 日韩在线高清观看一区二区三区 | 亚洲最大成人中文| 国产欧美日韩一区二区精品| 日韩中字成人| 精品国内亚洲2022精品成人| 狂野欧美激情性xxxx在线观看| 看片在线看免费视频| 成年女人永久免费观看视频| 国产精品不卡视频一区二区| 国内毛片毛片毛片毛片毛片| 成人高潮视频无遮挡免费网站| 欧美黑人巨大hd| 变态另类成人亚洲欧美熟女| 春色校园在线视频观看| 国产精品一区二区性色av| 99精品久久久久人妻精品| 午夜福利高清视频| 日韩高清综合在线| 成人高潮视频无遮挡免费网站| 69人妻影院| 色综合色国产| 在线观看午夜福利视频| 99视频精品全部免费 在线| 给我免费播放毛片高清在线观看| 国产精品一区二区性色av| 精品人妻熟女av久视频| 国产精品久久久久久亚洲av鲁大| aaaaa片日本免费| 一级a爱片免费观看的视频| 精品久久久久久成人av| 丰满的人妻完整版| 亚洲在线自拍视频| 国产精品三级大全| 在线国产一区二区在线| 一级毛片久久久久久久久女| 人妻制服诱惑在线中文字幕| 亚洲精品影视一区二区三区av| 国产视频一区二区在线看| 国产亚洲91精品色在线| 国产成人一区二区在线| 亚洲av免费高清在线观看| 午夜精品在线福利| 欧美日韩精品成人综合77777| 成人一区二区视频在线观看| 免费搜索国产男女视频| 欧美又色又爽又黄视频| 国产又黄又爽又无遮挡在线| 天美传媒精品一区二区| 成人av一区二区三区在线看| 九九爱精品视频在线观看| 国产精品爽爽va在线观看网站| 亚洲人成网站在线播放欧美日韩| 久久久国产成人精品二区| 免费高清视频大片| h日本视频在线播放| 久久久久性生活片| 五月玫瑰六月丁香| 婷婷亚洲欧美| 亚洲欧美激情综合另类| 国产成年人精品一区二区| 欧美日韩瑟瑟在线播放| 三级毛片av免费| 国产高清激情床上av| 午夜精品一区二区三区免费看| 色综合站精品国产| 黄色欧美视频在线观看| 成年版毛片免费区| 欧美性猛交黑人性爽| 淫妇啪啪啪对白视频| 亚洲国产欧美人成| a级一级毛片免费在线观看| 国产男人的电影天堂91| 乱码一卡2卡4卡精品| 天美传媒精品一区二区| 久久热精品热| 国产大屁股一区二区在线视频| 搡老妇女老女人老熟妇| 国产乱人视频| 久久人妻av系列| 欧美激情在线99| 日韩精品青青久久久久久| 日本一本二区三区精品| 色播亚洲综合网| 欧美日本亚洲视频在线播放| 欧美最黄视频在线播放免费| 97热精品久久久久久| 美女大奶头视频| 中文字幕高清在线视频| 亚洲欧美清纯卡通| 黄片wwwwww| 国产精品不卡视频一区二区| 国产欧美日韩一区二区精品| 国产蜜桃级精品一区二区三区| 久久久久久久久久成人| 18禁黄网站禁片午夜丰满| 日韩欧美三级三区| 国产黄片美女视频| 女人被狂操c到高潮| 色综合站精品国产|