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

    海水入侵區(qū)含水層中原核微生物多樣性和群落結(jié)構(gòu)特征及其意義

    2022-06-22 08:32:38桑石磊黃柏強(qiáng)王永智劉華雪許莉佳
    環(huán)境科學(xué)研究 2022年6期
    關(guān)鍵詞:咸水菌門(mén)鹽度

    桑石磊,黃柏強(qiáng),王永智,劉華雪,許莉佳,張 晉

    1. 生態(tài)環(huán)境部華南環(huán)境科學(xué)研究所,廣東 廣州 510535

    2. 河海大學(xué)長(zhǎng)江保護(hù)與綠色發(fā)展研究院,水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210098

    3. 中國(guó)科學(xué)院新疆生態(tài)與地理研究所,新疆 烏魯木齊 830011

    地下水是一種寶貴的飲用水資源,為我國(guó)近70%的人口提供飲用,95%以上的農(nóng)村人口主要飲用地下水,40%的農(nóng)田主要用地下水來(lái)灌溉[1]. 改革開(kāi)放40 多年來(lái),社會(huì)經(jīng)濟(jì)的快速發(fā)展和各種污染物的排放引起了地下淡水資源的缺乏[2]. 海水入侵再加上地下水需求的激增以及全球氣候變化、海平面上升等自然因素,進(jìn)一步加劇了地下淡水資源的短缺[3].特別是許多沿海地下含水層正面臨著人為污染[4]、地下鹵水[5]或過(guò)量抽取地下水[6]所導(dǎo)致的海水入侵作用等因素,嚴(yán)重威脅著地下水的安全. 雖然地下含水層大多是寡營(yíng)養(yǎng)的,但微生物幾乎是地下水環(huán)境中唯一的棲息者[7]. 研究表明,含水層是受微生物群落活動(dòng)影響的特殊的生態(tài)系統(tǒng)[8],與營(yíng)養(yǎng)物質(zhì)運(yùn)移、元素的生物地球化學(xué)循環(huán)和污染物降解作用密切相關(guān)[9].

    研究[10-12]表明,微生物群落是地下水生態(tài)系統(tǒng)的重要組成部分. 近年來(lái),多種微生物參數(shù),包括微生物群落、多樣性與理化特性的相互作用等,已被用來(lái)作為水質(zhì)評(píng)價(jià)的指標(biāo). 例如,為了尋找監(jiān)測(cè)地下水水質(zhì)的生物參數(shù),Unno 等[13]分析了微生物類群與水文化學(xué)之間的關(guān)系,發(fā)現(xiàn)微生物組分析是監(jiān)測(cè)地下水健康的有效工具. Ye 等[14]比較了海底地下水不同采樣位置的微生物指紋圖譜,并選擇了一些潛在的細(xì)菌類群(如共生單胞菌等)進(jìn)行生物修復(fù). 為有效了解和保護(hù)淡水地下水,研究地下含水層微生物群落的組成和功能是十分必要的. 然而,在海水入侵地區(qū),地下水中微生物群落的多樣性和結(jié)構(gòu)組成還未得到深入刻畫(huà),且其群落結(jié)構(gòu)與水文環(huán)境的相互作用關(guān)系仍不明確.

    珠江三角洲地區(qū)供應(yīng)將近5 000 萬(wàn)人口的用水需求,但幾乎一半的含水層中地下水含鹽度超過(guò)1 g/L[15-16]. 盡管水資源豐富,但由于土壤和水污染[17]以及珠江口嚴(yán)重的咸水入侵問(wèn)題[18],導(dǎo)致珠江三角洲一些地區(qū)地下水中總?cè)芙夤腆w含量甚至超過(guò)10 g/L[19]. 研究表明,地下水鹽度的來(lái)源是由全新世海侵作用遺留的古海水[20-21]. 為了提高對(duì)地下水的管理水平,學(xué)者對(duì)珠江三角洲含水層水文地球化學(xué)特征進(jìn)行了大量研究[21-22],但對(duì)該地區(qū)地下水鹽度梯度下微生物群落結(jié)構(gòu)的研究較少. 事實(shí)上,在世界上大多數(shù)國(guó)家的水質(zhì)監(jiān)測(cè)中,對(duì)地下生態(tài)系統(tǒng)中微生物群落的研究一直都未得到充分重視[23]. 研究[24]表明,微生物是地下水生態(tài)系統(tǒng)的重要組成部分,是評(píng)價(jià)地下水水質(zhì)不可缺少的參數(shù). 為有效地了解和保護(hù)地下淡水,對(duì)含水層微生物群落的組成和功能進(jìn)行調(diào)查研究很有必要.

    對(duì)地下水微生物多樣性的研究技術(shù)包括常規(guī)的培養(yǎng)技術(shù)和基于DNA 的分子生物學(xué)研究方法[25]. 由于環(huán)境生態(tài)系統(tǒng)中絕大多數(shù)微生物無(wú)法在人工培養(yǎng)基中輕易培養(yǎng),而基于DNA 的高通量測(cè)序方法在探索復(fù)雜環(huán)境中的微生物群落方面具有更強(qiáng)的優(yōu)勢(shì). 為研究珠江三角洲地下水中微生物群落的多樣性和分類組成,采用下一代測(cè)序(next generation sequencing,NGS)技術(shù),從水文生態(tài)學(xué)的角度描述微生物群落結(jié)構(gòu)與環(huán)境變量之間的關(guān)系;通過(guò)對(duì)珠江三角洲地下水水質(zhì)空間變化和微生物群落結(jié)構(gòu)的研究,分析地下水微生物群落對(duì)鹽度變化的響應(yīng),以期為地下咸水中微生物群落結(jié)構(gòu)及其環(huán)境相互作用提供參考.

    1 材料與方法

    1.1 采樣地點(diǎn)、樣品采集

    珠江三角洲屬亞熱帶季風(fēng)氣候,氣候溫暖濕潤(rùn),年均氣溫22 °C 左右,年均降水量1 600~2 000 mm,主要集中在4—10 月. 2017 年5 月,在廣州市中部和南部地區(qū)選擇22 個(gè)監(jiān)測(cè)井,包括5 個(gè)重度咸水監(jiān)測(cè)井(S)(TDS 濃度>10 g/L,TDS 表示溶解性總固體)、5個(gè)中度咸水監(jiān)測(cè)井(M1)(3 g/L<TDS 濃度<10 g/L)、3個(gè)微咸水監(jiān)測(cè)井(M2)(1 g/L<TDS 濃度<3 g/L)和9 個(gè)淡水監(jiān)測(cè)井(TDS 濃度<1 g/L).

    為保證取得的地下水樣可以表征含水層中地下水的性質(zhì),采樣之前均對(duì)監(jiān)測(cè)井進(jìn)行洗井工作,采用專業(yè)的地下水抽水泵并以小流量抽水洗井. 根據(jù)監(jiān)測(cè)井井管尺寸等信息計(jì)算管內(nèi)水體積,保證洗井過(guò)程抽水量大于管內(nèi)水體積的3 倍,以除去殘留在監(jiān)測(cè)井中的“死水”,待水體各項(xiàng)理化指標(biāo)(溫度、鹽度和pH)穩(wěn)定后再采集樣品. 現(xiàn)場(chǎng)測(cè)試的理化指標(biāo)主要包括溫度(temperature,T)、pH、電導(dǎo)率(electrical conductivity,EC)、溶解性總固體(total dissolved solids, TDS)、溶解氧 (dissolved oxygen, DO)和氧化還原電位(oxidation-reduction potential, ORP),測(cè)量?jī)x器為電導(dǎo)率儀(Orion 320C-01A, Thermo Scientific, USA). 地下水樣本都收集在無(wú)菌的10 L 塑料容器中,用于過(guò)濾收集微生物樣本,通過(guò)抽濾法在0.22 μm 的纖維濾膜上富集微生物,將濾膜保存于無(wú)菌離心管,封口膜密封后即刻存放在—80 ℃冰箱中備用;每個(gè)樣本再過(guò)濾500 mL,一式三份,用于理化分析. 所有樣品在運(yùn)輸過(guò)程中均保存在4 ℃的冰柜中,運(yùn)送回實(shí)驗(yàn)室立即處理.

    1.2 理化性質(zhì)分析

    理化分析參考Greenberg 等[26]的方法,為了減少水樣雜質(zhì)的干擾,在測(cè)試前,將水樣用0.45 μm 的濾膜過(guò)濾. 陽(yáng)離子(K+、Ca2+、Na+和Mg2+)濃度采用電感耦合等離子體發(fā)射光譜儀(ICAP 7600 ICP-OES,Thermo Scientific, USA)測(cè)試,陰離子(NO3—、Cl—和SO42—)濃度使用離子色譜儀(883 chromatograph,Metrohm, Switzerland)測(cè)定. HCO3—和CO32—濃度采用化學(xué)滴定法(DZ/T 0064.49—1993《地下水質(zhì)檢測(cè)方法》)檢測(cè).

    1.3 DNA 提取和PCR 擴(kuò)增及測(cè)序

    將水樣按照試劑盒MOBIO PowerSoil? DNA Isolation Kit (Qiagen/MO BIO Laboratories Inc., Carlsbad,CA, USA)操作步驟進(jìn)行DNA 提取. 用1%的瓊脂糖凝膠檢測(cè)基因組DNA 的提取效果,用NanoDrop 2000(Isogen Life Science, the Netherlands)光度儀檢測(cè)DNA的純度和濃度. 利用原核微生物16S rRNA V4 區(qū)引物515FmodF(5′-GTGCCAGCMGCCGCGGTAA-3′)、806RmodR (5′-GGACTACHVGGGTWTCTAAT-3′)進(jìn)行目標(biāo)基因的PCR 擴(kuò)增[27]. 擴(kuò)增程序:95 ℃預(yù)變性3 min,27 個(gè)循環(huán)(95 ℃變性30 s,55 ℃退火30 s, 72 ℃延伸30 s),然后72 ℃穩(wěn)定延伸10 min,最后在4 ℃下保存(PCR 儀,ABI GeneAmp? 9700 型). PCR 反應(yīng)體系:5×TransStart FastPfu 緩沖液4 μL,2.5 mmol/L dNTPs 2 μL,上游引物(5 μmol/L) 0.8 μL,下游引物(5 μmol/L ) 0.8 μL, TransStart FastPfu DNA 聚合酶0.4 μL,模板DNA 10 ng,ddH2O 補(bǔ)足至20 μL. 每個(gè)樣本設(shè)3個(gè)重復(fù). 按照NEBNext?Ultra? DNA Library Prep Kit for Illumina?(New England Biolabs, Ipswich,MA, USA)標(biāo)準(zhǔn)流程建庫(kù)操作,將樣品送到上海美吉生物科技有限公司,采用Illumina HiSeq 2500 PE250 測(cè)序平臺(tái)進(jìn)行擴(kuò)增子文庫(kù)的測(cè)序.

    利用Trimmomatic 軟件(V0.33)對(duì)測(cè)序數(shù)據(jù)過(guò)濾后,通過(guò)FLASH(V1.2.11)軟件對(duì)數(shù)據(jù)進(jìn)行拼接,再利用Mothur 軟件對(duì)拼接后的序列進(jìn)行質(zhì)量控制及過(guò)濾[28],得到有效的拼接片段. 過(guò)濾reads 尾部質(zhì)量值在20 以下的堿基,設(shè)置50 bp 的窗口,如果窗口內(nèi)的平均質(zhì)量值低于20,則從窗口開(kāi)始截去后端堿基,過(guò)濾質(zhì)控后50 bp 以下的reads,去除含N 堿基的reads;根據(jù)PE reads 之間的overlap 關(guān)系,將成對(duì)reads 拼接(merge)成一條序列,最小overlap 長(zhǎng)度為10 bp;拼接序列的overlap 區(qū)允許的最大錯(cuò)配比率為0.2,篩選不符合序列;根據(jù)序列首尾兩端的barcode 和引物區(qū)分樣品,并調(diào)整序列方向,barcode 允許的錯(cuò)配數(shù)為0,最大引物錯(cuò)配數(shù)為2.

    采用 RDP Classifier algorithm 方法和 Silva 16S rRNA 數(shù)據(jù)庫(kù)對(duì)物種進(jìn)行注釋分類(設(shè)置閾值為70%),將相似度高于97%的定義為一個(gè)OTU(operational taxonomic unit,OTU). α-多樣性計(jì)算:包括群落的多樣性指數(shù)(Shannon-Wiener)、豐富度指數(shù)(Chao1)、均勻度指數(shù)(Heip)和序列的覆蓋度(Good's coverage),使用Mothur 軟件完成計(jì)算[28]. β-多樣性的計(jì)算:比較樣本之間的相似性,在OTU 水平上利用Bray-Curtis 算法做PCoA 分析,用R 語(yǔ)言進(jìn)行統(tǒng)計(jì)分析和制圖[29]. 通過(guò)Cytoscape v2.8.3 對(duì)不同樣本間的物種豐度信息進(jìn)行共現(xiàn)性網(wǎng)絡(luò)分析[30],相關(guān)性熱圖用R 語(yǔ)言制圖[29],采用VIF(Variance inflation factor,VIF)分析去掉共線性環(huán)境因子,并進(jìn)一步使用Canoco 4.5 軟件進(jìn)行典范對(duì)應(yīng)分析(CCA). 該文涉及的所有原始序列已上傳至NCBI 數(shù)據(jù)庫(kù),其登錄號(hào)為SRP118856.

    2 結(jié)果與分析

    2.1 微生物群落的豐富度和多樣性

    該研究中,從22 個(gè)地下水樣品中共獲得814 701條基因序列,對(duì)所有樣品進(jìn)行抽平,共獲得1 427 個(gè)OTUs (97%的置信度),OTUs 的數(shù)目為415~804 個(gè),表明觀測(cè)到的微生物豐度在不同水體中存在差異. 微生物群落的Chao1 豐富度指數(shù)為650~1 047,Shannon-Wiener 多樣性指數(shù)為0.32~0.60. 一般來(lái)說(shuō),重度咸水監(jiān)測(cè)井(S)樣品包含的微生物群落的豐富度和均勻度指數(shù)高于淡水監(jiān)測(cè)井(F)樣品. 然而,微生物群落的豐富度和均勻度指數(shù)的最低值出現(xiàn)在微咸水(M2)樣品中.

    Shannon-Wiener 多樣性指數(shù)表明,含水層中微生物群落多樣性隨鹽度梯度發(fā)生顯著變化. 然而,微生物群落多樣性指數(shù)最高值出現(xiàn)在重度咸水(S)樣品中,而最低值出現(xiàn)在微咸水(M2)樣品中,多樣性指數(shù)的變化并不完全與鹽度梯度一致. 微生物序列的覆蓋度都在99%以上,表現(xiàn)出了較好的序列完整性.

    2.2 微生物群落的組成

    經(jīng)過(guò)數(shù)據(jù)庫(kù)的比對(duì),該研究中基因序列可以分為細(xì)菌和古菌兩大類,大多數(shù)的序列鑒定為細(xì)菌(見(jiàn)圖1). 具體來(lái)說(shuō),共分屬于47 個(gè)原核微生物門(mén)類,大多數(shù)序列從屬于變形菌門(mén)(Proteobacteria, 占比為50%~94%)、厚壁菌門(mén)(Firmicutes, 占比為3%~45%)、奇古菌門(mén)(Thaumarchaeota, 占比為0~29%)、放線菌門(mén)(Actinobacteria, 占 比 為 0~14%)、 廣 古 菌 門(mén)(Euryarchaeota, 占 比 為 0~13%)和 擬 桿 菌 門(mén)(Bacteroidetes, 占比為0~10%),還有一些序列占比在1%~10%之 間, 包 括 Acidobacteria、 Chloroflexi、Gemmatimonadetes 和Ignavibacteriae. 從結(jié)果來(lái)看,古菌在咸水中占比更高,廣古菌群(Euryarchaeota)在重度咸水(S)樣品中占比更高,而奇古菌門(mén)(Thaumarchaeota)在高濃度咸水樣品中占比更高.

    圖 1 微生物群落在門(mén)水平上的分布Fig.1 Microbial composition analysis based on taxonomic classification at the phylum level

    占比最高的變形菌門(mén)在綱水平上主要為γ-變形菌綱(Gammaproteobacteria),在重度咸水(S)、中度咸水(M1)、微咸水(M2)和淡水(F)監(jiān)測(cè)井樣品中分別為45%、30%、46%和38% (見(jiàn)圖1). 淡水微生物群落中β-變形菌綱在總序列中的占比高達(dá)34%,而在重度咸水樣品中其占比不超過(guò)2%. 與其他樣品相比,Epsilonproteobacteria 僅在中度咸水(M1)樣品中發(fā)現(xiàn). 此外,Deltaproteobacteria 在重度咸水(S)樣品中的豐度高于其他樣品,而Zetaproteobacteria 在所有樣品中的豐度都很低.

    2.3 微生物群落在不同鹽度梯度上的變化

    利用群落隸屬關(guān)系對(duì)含水層微生物群落進(jìn)行比較,探討微生物群落分布沿含水層鹽度梯度的空間變化(見(jiàn)圖2). 樣本間OTUs 網(wǎng)絡(luò)分析表明,相對(duì)于其他鹽度的含水層樣品,中度咸水(M1)和微咸水(M2)地下水樣品微生物群落結(jié)構(gòu)更相似. 因此,為了更加清晰地顯示微生物群落在鹽度梯度上的分布,把中度咸水(M1)和微咸水(M2)地下水樣品合并起來(lái)共稱為咸水樣品(M). 維恩圖展示了共有和獨(dú)有的OTUs 數(shù)量在3 個(gè)鹽度梯度上的分布(見(jiàn)圖2). 從22 個(gè)地下水樣品中產(chǎn)生的1 427 個(gè)OTUs 中可以看出,3 個(gè)鹽度差異的微生物群落之間共有的OTUs 占比高達(dá)59%(837 個(gè)OTUs). 此外,兩個(gè)樣品的微生物群落之間共有的OTUs 約占18% (257 個(gè)OTUs). 維恩圖結(jié)果顯示,不同鹽度的樣品之間有很多共同的OTUs,在一定程度上也說(shuō)明含水層中淡水和咸水之間存在流通和交換作用.

    特別是γ-變形菌綱(Gammaproteobacteria)在所有鹽度水體中占比都很高,微生物群落的相似性反映在屬水平的高豐度類群中,如假單胞菌屬(Pseudomonas)、不動(dòng)桿菌屬 (Acinetobacter)和Unclassified Enterobacteriaceae(見(jiàn)圖3). 重度咸水(S)中的高豐度類群是甲烷球菌屬(Methanococcus)、Candidatus Nitrosoarchaeum和海洋桿菌屬(Marinobacter),分別隸屬于奇古菌門(mén)(Thaumarchaeota)、 廣古菌門(mén)(Euryarchaeota)和γ-變形菌綱(Gammaproteobacteria).然而,在淡水樣品中占優(yōu)勢(shì)的類群是β-變形菌綱(Betaproteobacteria),包括Norank Hydrogenophilaceae、Unclassified Methylophilaceae、 Denitratisoma 和Unclassified Rhodocyclaceae.

    圖 2 樣本之間OTUs 網(wǎng)絡(luò)分析和維恩圖Fig.2 Network analysis of OTUs between samples obtained Venn diagram

    圖 3 不同鹽度梯度下微生物群落(屬水平)的相對(duì)豐度Fig.3 The relative abundance of microbial taxa (genus level) along the salinity gradient

    2.4 微生物群落結(jié)構(gòu)及其與環(huán)境因子的相互關(guān)系

    對(duì)微生物群落與主要環(huán)境因子進(jìn)行典型相關(guān)分析. 結(jié)果顯示,前兩個(gè)CCA 軸中的環(huán)境參數(shù)共同解釋了微生物群落組成差異的23.7% (見(jiàn)圖4),TDS 濃度(r2=0.63,P=0.001)和TN 濃度(r2=0.8,P=0.001)對(duì)微生物分布的影響最強(qiáng),而溫度是最弱的環(huán)境變量. Spearman 相關(guān)分析(見(jiàn)表1)表明,厚壁菌相對(duì)豐 度與ORP 及HCO3—、DO 濃度均呈顯著負(fù)相關(guān)(P<0.05). 擬桿菌相對(duì)豐度與ORP、pH 均呈顯著正相關(guān),與Fe2+濃度呈顯著負(fù)相關(guān),而α-變形桿菌相對(duì)豐度與TP 濃度呈顯著正相關(guān)(P<0.05). β-桿菌相對(duì)豐度與HCO3—濃度呈顯著正相關(guān),與Cl—、TDS 濃度均呈顯著負(fù)相關(guān)(P<0.05). γ-變形菌綱相對(duì)豐度與TN濃度呈顯著正相關(guān)(P<0.01). 在屬水平上,短小桿菌屬相對(duì)豐度與HCO3—、DO 濃度均呈負(fù)相關(guān),不動(dòng)桿菌屬相對(duì)豐度與HCO3—濃度呈負(fù)相關(guān),假單胞菌屬相對(duì)豐度與ORP 及HCO3—、DO 濃度均呈負(fù)相 關(guān)(P<0.05).

    圖 4 微生物群落結(jié)構(gòu)與環(huán)境因子的典范對(duì)應(yīng)分析Fig.4 CCA analysis between microbial community structure and environmental fact

    表 1 水體樣品在不同分類水平上的微生物群落相對(duì)豐度與環(huán)境因子之間的Spearman 相關(guān)性分析Table 1 Environmental factors associated with variations of the microbial communities at the phylum/class or genus level in the four water samples

    3 討論

    3.1 微生物的豐富度和多樣性

    長(zhǎng)期以來(lái),對(duì)水體鹽度梯度下微生物多樣性的研究主要集中在地表水系統(tǒng)[31-32]. 該研究揭示了地下水中原核微生物群落的豐富度和多樣性,這些原核微生物所在的地下環(huán)境具有復(fù)雜的地球化學(xué)性質(zhì),可能發(fā)生了劇烈的生物地球化學(xué)過(guò)程. 由Chao1 和Shannon-Wiener 多樣性指數(shù)反映的微生物群落豐富度和均勻度沿鹽度梯度發(fā)生了變化,但變化并不明顯. 總體而言,在重度咸水(TDS 濃度>10 g/L)中觀察到的微生物群落豐富度和多樣性最高,其次是淡水. 該研究結(jié)果與在海灣河口[33]和原始含水層[34]中觀察到的結(jié)果相反,他們認(rèn)為鹽度的增加會(huì)導(dǎo)致細(xì)菌多樣性的減少,而在相對(duì)穩(wěn)定的波羅的海河口[35], Shannon-Wiener多樣性指數(shù)沿鹽度梯度的變化不明顯. 在微咸水(1 g/L<TDS 濃度<3 g/L)樣品中發(fā)現(xiàn)微生物群落的豐富度和均勻度最低,也就是說(shuō),雖然微生物多樣性并沒(méi)有完全遵循鹽度梯度變化規(guī)律,但在TDS 濃度>1 g/L 的咸水中,微生物群落多樣性隨著鹽度的增加而增加[23,36-37]. 因此,微生物群落的α-多樣性與其所處的環(huán)境條件關(guān)系非常密切. 在環(huán)境波動(dòng)較大的水體環(huán)境中[33-34],地下水中微生物的Shannon-Wiener 多樣性指數(shù)由秋冬季的3.22±0.28顯著變化為春夏季的1.31±0.35,這是由于營(yíng)養(yǎng)鹽、有機(jī)碳的輸入以及地下水位的強(qiáng)烈波動(dòng),引起微生物群落多樣性在不同季節(jié)具有顯著性差異. 而在該研究中,地下環(huán)境條件相對(duì)比較穩(wěn)定,微生物的多樣性的變化并不明顯.

    3.2 微生物的群落結(jié)構(gòu)

    結(jié)果顯示,地下水中的優(yōu)勢(shì)菌群是變形菌門(mén)、厚壁菌門(mén)、放線菌門(mén)和擬桿菌門(mén),與其他研究結(jié)果[38-39]相似,而古菌的群落結(jié)構(gòu)主要是奇古菌門(mén)和廣古菌門(mén). 值得注意的是,γ-變形菌綱豐富度在所有樣本中都比較高(見(jiàn)圖1). 一般來(lái)說(shuō),海水中的優(yōu)勢(shì)細(xì)菌是α-或γ-變形菌綱,β-變形細(xì)菌在淡水中占優(yōu)勢(shì)[40-41]. 該研究中,與咸水樣品相比,淡水樣品中的β-變形菌綱占比相對(duì)較高. 在咸水中檢測(cè)到豐度較高的類群主要包括脫硫弧菌屬、甲烷球菌屬和海洋桿菌屬,在其他樣品中幾乎未檢出(見(jiàn)圖3). 脫硫弧菌屬和產(chǎn)甲烷球菌屬中存在豐富的異養(yǎng)厭氧微生物,表明樣品的氧化還原電位較低. 先前的地質(zhì)調(diào)查中就發(fā)現(xiàn)地層中有大量的甲烷氣體和臭雞蛋味的硫化氫氣體,表明含水層中存在微生物硫酸鹽還原和產(chǎn)甲烷作用[20]. 該研究中一些高濃度咸水樣品中檢測(cè)到的脫硫弧菌屬和甲烷球菌厭氧微生物,可能是地層中二氧化硫和甲烷氣體生成的微生物類群. 在咸水樣品中有豐度較高的海洋桿菌屬,是嗜鹽或者耐鹽性的物種,一般生存在海洋環(huán)境中[42],而這些海洋物種序列很可能是在海水入侵的作用下被帶入海濱含水層. 在其他海濱土壤的研究中也發(fā)現(xiàn)了海洋物種序列[43],研究者認(rèn)為很可能是由于海水入侵作用被帶入海濱土壤. 這也表明了海濱含水層是海洋與內(nèi)陸物種交換的重要場(chǎng)所,對(duì)微生物的地理分布起到一定作用[44].

    較長(zhǎng)的水體停留時(shí)間和相對(duì)穩(wěn)定的環(huán)境有助于建立穩(wěn)定的土著微生物群落生態(tài)系統(tǒng)[39]. 咸水與淡水微生物群落組成差異較大(見(jiàn)圖3),然而維恩圖結(jié)果(見(jiàn)圖2)顯示,不同鹽度水體之間存在大量共有的OTUs,表明咸水與淡水之間可能存在著緊密的流動(dòng)和交換聯(lián)系[20].

    3.3 微生物群落與環(huán)境的關(guān)系

    一般認(rèn)為,水生微生物群落結(jié)構(gòu)的時(shí)空變化和分布是環(huán)境因素綜合作用的結(jié)果. 該研究中,微生物群落與主要環(huán)境因子的典型相關(guān)分析表明,前兩個(gè)CCA 軸中的環(huán)境參數(shù)共同解釋了微生物群落組成差異的23.7% (見(jiàn)圖4),TDS 濃度(r2=0.63,P=0.001)和TN 濃度(r2=0.8,P=0.001)對(duì)微生物分布的影響最強(qiáng),而溫度的影響最弱. 根據(jù)對(duì)不同生態(tài)環(huán)境的全球調(diào)查結(jié)果,鹽度被視為是微生物群落變異的主要環(huán)境影響因子[45-46]. 該研究結(jié)果表明,TDS 和TN 是微生物群落最有力的驅(qū)動(dòng)因素(見(jiàn)圖3). 鹽度是一個(gè)至關(guān)重要的環(huán)境變量,它與β-變形菌綱豐度呈負(fù)相關(guān),這與在黃海沿岸海底地下水中的研究結(jié)果[14]一致. 有研究[21]表明,珠江三角洲地下水體中存在高達(dá)8.6×109kg 的氨氮污染物,并且這些數(shù)量巨大的氨氮有機(jī)物從地下弱透水層不斷釋放,可以解釋TN 對(duì)微生物群落的重要影響. 根據(jù)廣東省水文地質(zhì)調(diào)查結(jié)果,珠江三角洲地下水中氨濃度高達(dá)560 mg/L,主要源于上覆的全新世-更新世中富含有機(jī)質(zhì)的含水層,通過(guò)地下水輸送擴(kuò)散進(jìn)入含水層[21],可能會(huì)逐漸遷移到河水和沿岸海水中,破壞生態(tài)平衡,需引起高度重視.

    4 結(jié)論

    a) 珠三角地區(qū)海水入侵區(qū)含水層中原核微生物群落中細(xì)菌以變形菌門(mén)、厚壁菌門(mén)、放線菌門(mén)和擬桿菌門(mén)為主,而古菌以奇古菌門(mén)和廣古菌門(mén)為主,且古菌在咸水中豐富度更高.

    b) 環(huán)境條件相對(duì)比較穩(wěn)定的含水層中,微生物多樣性變化并不明顯.

    c) 該研究中對(duì)地下水微生物群落影響最大的環(huán)境因子是溶解性總固體(TDS)和總氮(TN)濃度,含水層中大量的氨氮可能會(huì)對(duì)生態(tài)平衡造成一定影響,需引起高度重視.

    猜你喜歡
    咸水菌門(mén)鹽度
    野生樹(shù)鼩與人工飼養(yǎng)樹(shù)鼩消化道不同部位微生物組成的比較研究
    饑餓與重?cái)z食對(duì)河蟹腸道菌群結(jié)構(gòu)的影響
    昆蟲(chóng)體內(nèi)微生物多樣性的影響因素研究進(jìn)展
    妊娠期糖尿病腸道菌群變化及臨床價(jià)值分析
    聊城市地下咸水地質(zhì)特征與綜合開(kāi)發(fā)利用分析
    微咸水滴灌能提高紅棗果實(shí)品質(zhì)
    鹽度和pH對(duì)細(xì)角螺耗氧率和排氨率的影響
    鹽度脅迫對(duì)入侵生物福壽螺的急性毒性效應(yīng)
    適用于高鹽度和致密巖層驅(qū)油的表面活性劑
    膠州灣夏季鹽度長(zhǎng)期輸運(yùn)機(jī)制分析
    videos熟女内射| 一区二区三区四区激情视频| 免费人妻精品一区二区三区视频| 黑人巨大精品欧美一区二区蜜桃| 最黄视频免费看| 女人久久www免费人成看片| 乱人伦中国视频| 亚洲精品一区蜜桃| 午夜激情av网站| 国精品久久久久久国模美| 天堂8中文在线网| 色视频在线一区二区三区| 国产日韩一区二区三区精品不卡| 亚洲精品乱久久久久久| 91国产中文字幕| 国产精品久久久久久精品古装| 卡戴珊不雅视频在线播放| 亚洲av福利一区| 极品人妻少妇av视频| 91aial.com中文字幕在线观看| 色网站视频免费| 久久鲁丝午夜福利片| 日韩电影二区| 亚洲第一av免费看| 亚洲精品一区蜜桃| 狠狠精品人妻久久久久久综合| 国产精品国产三级国产专区5o| 一边摸一边抽搐一进一出视频| 日韩免费高清中文字幕av| 精品第一国产精品| 天堂8中文在线网| 免费黄色在线免费观看| 蜜桃在线观看..| 亚洲av欧美aⅴ国产| 最新在线观看一区二区三区 | 国产精品人妻久久久影院| av卡一久久| 久久久久久人人人人人| 伦理电影大哥的女人| a级毛片在线看网站| av网站在线播放免费| 国语对白做爰xxxⅹ性视频网站| 男人舔女人的私密视频| 天天添夜夜摸| 最新在线观看一区二区三区 | 99久久综合免费| 纯流量卡能插随身wifi吗| 亚洲在久久综合| 只有这里有精品99| 看免费av毛片| 国产欧美日韩一区二区三区在线| 久久女婷五月综合色啪小说| 十分钟在线观看高清视频www| netflix在线观看网站| 九九爱精品视频在线观看| 十八禁网站网址无遮挡| 国产精品秋霞免费鲁丝片| 国产成人免费无遮挡视频| 蜜桃国产av成人99| 成人亚洲精品一区在线观看| 人成视频在线观看免费观看| 啦啦啦在线观看免费高清www| av电影中文网址| 精品亚洲成a人片在线观看| 国产精品一区二区在线观看99| 日韩精品有码人妻一区| 亚洲综合精品二区| 秋霞伦理黄片| 亚洲美女黄色视频免费看| 制服人妻中文乱码| 婷婷成人精品国产| 精品酒店卫生间| 99热国产这里只有精品6| 最近中文字幕高清免费大全6| 国产精品久久久久久精品古装| 国产精品二区激情视频| 亚洲欧美成人综合另类久久久| 国产av国产精品国产| 女人被躁到高潮嗷嗷叫费观| 天天添夜夜摸| 一边亲一边摸免费视频| 亚洲欧美激情在线| 久久久久精品人妻al黑| 国产麻豆69| netflix在线观看网站| 日日啪夜夜爽| 国产女主播在线喷水免费视频网站| 可以免费在线观看a视频的电影网站 | 不卡av一区二区三区| 国产精品久久久av美女十八| 伊人亚洲综合成人网| 日韩av免费高清视频| 成年人午夜在线观看视频| 久久狼人影院| 汤姆久久久久久久影院中文字幕| 久久精品久久精品一区二区三区| av在线app专区| 中文字幕另类日韩欧美亚洲嫩草| 桃花免费在线播放| 一边摸一边抽搐一进一出视频| 亚洲精品久久成人aⅴ小说| 日韩人妻精品一区2区三区| 夫妻午夜视频| 老司机靠b影院| 亚洲中文av在线| 日本wwww免费看| 天天躁日日躁夜夜躁夜夜| 亚洲精品一区蜜桃| av网站免费在线观看视频| 亚洲国产中文字幕在线视频| 免费看av在线观看网站| 久久女婷五月综合色啪小说| av网站在线播放免费| 水蜜桃什么品种好| 不卡视频在线观看欧美| 国产成人91sexporn| 高清欧美精品videossex| 国产免费又黄又爽又色| 国产精品成人在线| 咕卡用的链子| 日韩视频在线欧美| 两性夫妻黄色片| 久久99热这里只频精品6学生| 女人被躁到高潮嗷嗷叫费观| 久久人妻熟女aⅴ| 精品久久蜜臀av无| 亚洲精品成人av观看孕妇| 国产不卡av网站在线观看| 国产精品无大码| 精品少妇久久久久久888优播| 天天躁夜夜躁狠狠躁躁| 国产成人精品久久二区二区91 | 欧美日韩亚洲高清精品| 精品国产乱码久久久久久男人| 啦啦啦在线观看免费高清www| 丝袜脚勾引网站| kizo精华| 极品人妻少妇av视频| 2021少妇久久久久久久久久久| 一个人免费看片子| 捣出白浆h1v1| 夫妻性生交免费视频一级片| 又大又爽又粗| 国产精品无大码| 国产成人欧美在线观看 | 亚洲国产精品国产精品| 日韩精品免费视频一区二区三区| 欧美日本中文国产一区发布| 在线天堂最新版资源| 亚洲一码二码三码区别大吗| 夫妻性生交免费视频一级片| 精品国产乱码久久久久久男人| 成人毛片60女人毛片免费| 男女国产视频网站| 熟妇人妻不卡中文字幕| 欧美在线一区亚洲| 国产免费视频播放在线视频| 欧美精品高潮呻吟av久久| 美女大奶头黄色视频| 久久影院123| 久久人人97超碰香蕉20202| 叶爱在线成人免费视频播放| 亚洲精品国产色婷婷电影| 丝袜在线中文字幕| 欧美在线一区亚洲| 国产免费视频播放在线视频| 男女无遮挡免费网站观看| 国产精品偷伦视频观看了| av卡一久久| 巨乳人妻的诱惑在线观看| 亚洲av国产av综合av卡| 高清不卡的av网站| 国产av一区二区精品久久| 一区二区三区精品91| 丝袜美足系列| 欧美激情 高清一区二区三区| 成人国产麻豆网| 日本欧美国产在线视频| 亚洲,欧美,日韩| 99热网站在线观看| 最近2019中文字幕mv第一页| 亚洲中文av在线| av又黄又爽大尺度在线免费看| 亚洲欧美激情在线| 中文字幕人妻丝袜制服| 校园人妻丝袜中文字幕| 欧美97在线视频| 日本av免费视频播放| 丰满少妇做爰视频| 99re6热这里在线精品视频| 日韩一卡2卡3卡4卡2021年| 亚洲欧美一区二区三区国产| 性高湖久久久久久久久免费观看| 国产成人啪精品午夜网站| 婷婷成人精品国产| 男女免费视频国产| 午夜福利视频精品| 熟女av电影| 精品一区二区免费观看| 欧美亚洲日本最大视频资源| 男女下面插进去视频免费观看| 天美传媒精品一区二区| 亚洲第一av免费看| 国产国语露脸激情在线看| 亚洲国产精品成人久久小说| 国产成人精品在线电影| 人人澡人人妻人| videos熟女内射| 国产精品国产三级专区第一集| 色播在线永久视频| 黄片小视频在线播放| 七月丁香在线播放| 天堂俺去俺来也www色官网| 国产熟女欧美一区二区| 人人妻人人添人人爽欧美一区卜| 在线观看免费视频网站a站| 免费观看a级毛片全部| 欧美国产精品va在线观看不卡| 亚洲久久久国产精品| 可以免费在线观看a视频的电影网站 | 天天操日日干夜夜撸| 一二三四中文在线观看免费高清| 国产一卡二卡三卡精品 | 在线亚洲精品国产二区图片欧美| 一级片免费观看大全| 日日摸夜夜添夜夜爱| 日韩一区二区视频免费看| bbb黄色大片| 免费黄色在线免费观看| 大话2 男鬼变身卡| 一区二区三区乱码不卡18| 欧美黄色片欧美黄色片| 黑人猛操日本美女一级片| 99久久人妻综合| 国产极品天堂在线| 亚洲四区av| 成年人免费黄色播放视频| 熟女少妇亚洲综合色aaa.| 男人舔女人的私密视频| 香蕉国产在线看| www日本在线高清视频| 精品人妻熟女毛片av久久网站| 亚洲综合色网址| 亚洲精品国产av蜜桃| 日韩人妻精品一区2区三区| 丝袜美足系列| 亚洲国产日韩一区二区| 亚洲国产中文字幕在线视频| 国产精品麻豆人妻色哟哟久久| 亚洲熟女精品中文字幕| 麻豆乱淫一区二区| 成人毛片60女人毛片免费| 久久韩国三级中文字幕| 亚洲,欧美精品.| 国产一区亚洲一区在线观看| 黑人巨大精品欧美一区二区蜜桃| 国产黄色免费在线视频| 成年人午夜在线观看视频| 女性生殖器流出的白浆| 国产精品女同一区二区软件| 精品国产一区二区三区四区第35| 日本vs欧美在线观看视频| 天堂俺去俺来也www色官网| 成年美女黄网站色视频大全免费| 毛片一级片免费看久久久久| 亚洲久久久国产精品| 欧美日韩亚洲国产一区二区在线观看 | 亚洲成人手机| 国产探花极品一区二区| 成人影院久久| 午夜影院在线不卡| 最新在线观看一区二区三区 | 999久久久国产精品视频| 国产成人精品无人区| 美女福利国产在线| 久久鲁丝午夜福利片| 欧美少妇被猛烈插入视频| 91成人精品电影| 欧美日韩视频精品一区| 国产精品国产三级国产专区5o| 免费观看av网站的网址| 毛片一级片免费看久久久久| 亚洲欧美成人精品一区二区| 国产激情久久老熟女| 国产野战对白在线观看| 欧美日韩一级在线毛片| 丝袜在线中文字幕| 国产精品偷伦视频观看了| 亚洲一区中文字幕在线| 男女之事视频高清在线观看 | 夜夜骑夜夜射夜夜干| 亚洲精品美女久久av网站| 午夜激情久久久久久久| 黄网站色视频无遮挡免费观看| 国产成人精品在线电影| 国产一区二区在线观看av| a 毛片基地| 亚洲男人天堂网一区| 少妇人妻久久综合中文| 在现免费观看毛片| 精品一区二区三卡| 亚洲av日韩在线播放| 日韩av免费高清视频| 亚洲精品在线美女| a级片在线免费高清观看视频| 如日韩欧美国产精品一区二区三区| 少妇猛男粗大的猛烈进出视频| 国产野战对白在线观看| 纵有疾风起免费观看全集完整版| 激情五月婷婷亚洲| 亚洲精品国产av蜜桃| 色网站视频免费| 日韩av免费高清视频| 妹子高潮喷水视频| 亚洲av电影在线进入| 1024视频免费在线观看| 国产欧美日韩综合在线一区二区| 精品久久久精品久久久| 精品一品国产午夜福利视频| 成年美女黄网站色视频大全免费| 成年女人毛片免费观看观看9 | 黄色 视频免费看| 中文字幕另类日韩欧美亚洲嫩草| 久久久久久久久免费视频了| 天天影视国产精品| 超色免费av| 一区二区三区乱码不卡18| 免费看不卡的av| 日韩大片免费观看网站| 精品少妇一区二区三区视频日本电影 | 日本欧美国产在线视频| 欧美黑人精品巨大| 看非洲黑人一级黄片| 国产精品熟女久久久久浪| 18禁国产床啪视频网站| 看非洲黑人一级黄片| 久久精品国产综合久久久| 精品亚洲成国产av| 亚洲av欧美aⅴ国产| 天堂8中文在线网| 高清在线视频一区二区三区| 国产成人精品久久久久久| 中国国产av一级| 十八禁高潮呻吟视频| 亚洲精品自拍成人| 久久青草综合色| 无遮挡黄片免费观看| 18在线观看网站| netflix在线观看网站| 久久天堂一区二区三区四区| 精品一区二区免费观看| 麻豆精品久久久久久蜜桃| 国产精品一区二区在线观看99| 国产精品 欧美亚洲| 欧美 日韩 精品 国产| 交换朋友夫妻互换小说| 欧美精品亚洲一区二区| 一级片免费观看大全| 国产精品久久久久久精品电影小说| 最近最新中文字幕大全免费视频 | 一个人免费看片子| 亚洲精品乱久久久久久| 日本91视频免费播放| 国产成人免费观看mmmm| 国产亚洲精品第一综合不卡| 午夜老司机福利片| 高清在线视频一区二区三区| 各种免费的搞黄视频| 久久精品国产a三级三级三级| 麻豆av在线久日| netflix在线观看网站| 久久精品国产亚洲av涩爱| 青春草亚洲视频在线观看| 亚洲欧美成人精品一区二区| 午夜免费鲁丝| 看十八女毛片水多多多| 一区二区三区乱码不卡18| 精品国产露脸久久av麻豆| 欧美日韩国产mv在线观看视频| 色网站视频免费| 美国免费a级毛片| 丝袜美腿诱惑在线| 不卡av一区二区三区| 欧美xxⅹ黑人| 国产精品99久久99久久久不卡 | 国产一级毛片在线| 成人国产麻豆网| 在线亚洲精品国产二区图片欧美| 日韩制服骚丝袜av| 一级爰片在线观看| √禁漫天堂资源中文www| 赤兔流量卡办理| 制服人妻中文乱码| 欧美日本中文国产一区发布| 国产免费一区二区三区四区乱码| 丝袜美腿诱惑在线| 亚洲欧美中文字幕日韩二区| 国产极品天堂在线| 18禁动态无遮挡网站| 成人国产麻豆网| 乱人伦中国视频| 一级毛片 在线播放| tube8黄色片| av网站在线播放免费| 精品人妻在线不人妻| 男女国产视频网站| 十八禁网站网址无遮挡| 国产成人免费无遮挡视频| 19禁男女啪啪无遮挡网站| 91精品伊人久久大香线蕉| 亚洲人成77777在线视频| 日韩免费高清中文字幕av| 国产成人精品久久久久久| 国产精品久久久久久精品电影小说| 欧美精品人与动牲交sv欧美| 激情视频va一区二区三区| 免费观看a级毛片全部| 国产欧美日韩一区二区三区在线| 七月丁香在线播放| 男女无遮挡免费网站观看| 激情五月婷婷亚洲| 亚洲,欧美,日韩| 最近中文字幕高清免费大全6| 成人免费观看视频高清| 老司机深夜福利视频在线观看 | 看免费成人av毛片| 精品卡一卡二卡四卡免费| 欧美日韩一区二区视频在线观看视频在线| 丝袜脚勾引网站| 欧美亚洲 丝袜 人妻 在线| 亚洲精品国产区一区二| 国产成人91sexporn| 少妇精品久久久久久久| 成年女人毛片免费观看观看9 | 国产毛片在线视频| 中文字幕高清在线视频| 男女免费视频国产| 亚洲熟女毛片儿| 亚洲一卡2卡3卡4卡5卡精品中文| videos熟女内射| av网站在线播放免费| 久久久欧美国产精品| 人人妻人人添人人爽欧美一区卜| 免费高清在线观看日韩| 曰老女人黄片| 制服诱惑二区| 涩涩av久久男人的天堂| svipshipincom国产片| 高清视频免费观看一区二区| 王馨瑶露胸无遮挡在线观看| 午夜福利网站1000一区二区三区| 久久久精品区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 免费少妇av软件| 亚洲欧洲精品一区二区精品久久久 | a级毛片在线看网站| 精品少妇内射三级| 纵有疾风起免费观看全集完整版| 天堂俺去俺来也www色官网| 91国产中文字幕| 777久久人妻少妇嫩草av网站| 日日啪夜夜爽| 99精国产麻豆久久婷婷| 亚洲,一卡二卡三卡| 午夜福利,免费看| 国产精品 欧美亚洲| 国产成人系列免费观看| e午夜精品久久久久久久| 久久精品国产a三级三级三级| 中文字幕人妻熟女乱码| 韩国高清视频一区二区三区| 亚洲av中文av极速乱| 在线 av 中文字幕| 亚洲第一区二区三区不卡| 无限看片的www在线观看| 天天添夜夜摸| av卡一久久| 伦理电影大哥的女人| 啦啦啦在线免费观看视频4| 中文字幕人妻丝袜一区二区 | 日日啪夜夜爽| 亚洲av成人精品一二三区| 三上悠亚av全集在线观看| 哪个播放器可以免费观看大片| 性色av一级| av在线观看视频网站免费| 最近2019中文字幕mv第一页| 99精品久久久久人妻精品| 丝袜喷水一区| 两个人看的免费小视频| 99国产综合亚洲精品| 最近的中文字幕免费完整| 精品国产乱码久久久久久小说| 国产精品麻豆人妻色哟哟久久| 久久ye,这里只有精品| 国产精品.久久久| 国产精品嫩草影院av在线观看| 亚洲人成电影观看| 97在线人人人人妻| 国产人伦9x9x在线观看| 国产男女超爽视频在线观看| 少妇人妻精品综合一区二区| 精品国产一区二区久久| 日韩一区二区视频免费看| av在线播放精品| 大话2 男鬼变身卡| 亚洲精品国产区一区二| 国产精品人妻久久久影院| 国产成人免费无遮挡视频| 精品国产一区二区三区久久久樱花| 国产 精品1| 99精品久久久久人妻精品| 天天影视国产精品| 我的亚洲天堂| 51午夜福利影视在线观看| 精品国产超薄肉色丝袜足j| 国产成人91sexporn| 午夜免费男女啪啪视频观看| 嫩草影院入口| 亚洲精品日韩在线中文字幕| 高清在线视频一区二区三区| 激情五月婷婷亚洲| 老司机影院毛片| 狂野欧美激情性bbbbbb| 少妇人妻 视频| 一区二区日韩欧美中文字幕| 美女福利国产在线| 国产精品 国内视频| 天天躁夜夜躁狠狠躁躁| 国产精品 欧美亚洲| 黄色视频不卡| 国产精品久久久久久精品电影小说| 免费高清在线观看视频在线观看| 亚洲伊人色综图| 国产亚洲最大av| 国产男女超爽视频在线观看| 一边摸一边做爽爽视频免费| 国产极品粉嫩免费观看在线| 精品人妻一区二区三区麻豆| 爱豆传媒免费全集在线观看| 国产色婷婷99| 欧美国产精品一级二级三级| 精品国产一区二区三区四区第35| 777米奇影视久久| 欧美97在线视频| 国产精品三级大全| 1024香蕉在线观看| 观看美女的网站| 国产伦理片在线播放av一区| 青春草亚洲视频在线观看| 熟女av电影| 青青草视频在线视频观看| 老汉色av国产亚洲站长工具| 国产免费视频播放在线视频| 大陆偷拍与自拍| 不卡av一区二区三区| 男的添女的下面高潮视频| 一边摸一边做爽爽视频免费| 成人黄色视频免费在线看| 999久久久国产精品视频| 操美女的视频在线观看| 欧美在线黄色| 欧美精品av麻豆av| 国产探花极品一区二区| 嫩草影院入口| 亚洲第一区二区三区不卡| 水蜜桃什么品种好| 国产黄频视频在线观看| 中文字幕制服av| 一级毛片 在线播放| 男人添女人高潮全过程视频| www.熟女人妻精品国产| 男女床上黄色一级片免费看| xxxhd国产人妻xxx| 波野结衣二区三区在线| 丝袜美足系列| 香蕉丝袜av| 亚洲欧美中文字幕日韩二区| 男女边摸边吃奶| 亚洲精品国产区一区二| 一区二区av电影网| 国产成人系列免费观看| 免费av中文字幕在线| 在线精品无人区一区二区三| 毛片一级片免费看久久久久| 欧美激情高清一区二区三区 | 精品免费久久久久久久清纯 | 久久久久久久久免费视频了| 五月开心婷婷网| 又黄又粗又硬又大视频| 老鸭窝网址在线观看| 婷婷色av中文字幕| 亚洲国产av新网站| 国产精品香港三级国产av潘金莲 | 99热全是精品| 国产一区亚洲一区在线观看| av有码第一页| 9热在线视频观看99| 国产精品一国产av| 精品人妻一区二区三区麻豆| 国产97色在线日韩免费| 日韩欧美一区视频在线观看| 亚洲精品av麻豆狂野| 精品免费久久久久久久清纯 | 777米奇影视久久| 青春草视频在线免费观看| 国产精品三级大全| 欧美黄色片欧美黄色片| 精品少妇内射三级| 男女下面插进去视频免费观看| 日韩中文字幕欧美一区二区 | 国产一区亚洲一区在线观看| 国产不卡av网站在线观看| 丁香六月欧美| 欧美xxⅹ黑人|