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

    鄱陽湖7個河口水體中細菌多樣性和組成特征

    2016-09-13 00:38:44圣平于一尊田曉娟黃江麗張志紅丁建南
    關鍵詞:饒河樂安鄱陽湖

    圣平,于一尊,田曉娟,黃江麗,張志紅,2,丁建南*

    (1.江西省科學院生物資源研究所,江西 南昌330096;2. 江西省鄱陽湖重點實驗室,江西 南昌 330096)

    鄱陽湖7個河口水體中細菌多樣性和組成特征

    圣平1,于一尊1,田曉娟1,黃江麗1,張志紅1,2,丁建南1*

    (1.江西省科學院生物資源研究所,江西 南昌330096;2. 江西省鄱陽湖重點實驗室,江西 南昌 330096)

    本文對鄱陽湖7個河口(贛江、撫河、信江、樂安河、昌江、饒河和修河)水體的理化特性進行了分析,并采用細菌16S rRNA基因MiSeq測序的方法對其多樣性和組成進行了研究。結(jié)果表明,樂安河水體總氮含量最高,信江水體總磷含量最高。樂安河、修河和饒河水體細菌豐度高于其它河口,而修河、昌江和饒河水體細菌多樣性則相對較高。水體細菌群落組成分析結(jié)果表明,在7個河口水體中變形菌門細菌均占優(yōu)勢,其次為放線菌門和擬桿菌門細菌。而在綱水平,在7個河口水體中β-變形菌綱細菌占優(yōu)勢,其次為α-變形菌綱細菌和放線菌綱細菌。相關性分析結(jié)果表明,水體pH和細菌群落多樣性呈顯著負相關關系,而總氮和氨態(tài)氮含量和細菌群落多樣性呈顯著正相關性。

    鄱陽湖;7大河口;不同河口水樣;細菌多樣性;細菌群落組成

    圣平, 于一尊, 田曉娟, 黃江麗, 張志紅, 丁建南. 鄱陽湖7個河口水體中細菌多樣性和組成特征[J]. 農(nóng)業(yè)現(xiàn)代化研究,2016, 37(3): 606-612.

    Sheng P, Yu Y Z, Tian X J, Huang J L, Zhang Z H, Ding J N. Bacterial diversities and compositions in seven different estuarine water columns of Poyang Lake[J]. Research of Agricultural Modernization, 2016, 37(3): 606-612.

    河口地段是鄱陽湖的生態(tài)敏感部位,一些重要的污染參數(shù)會在河口部位得到集中體現(xiàn),對河口地段進行監(jiān)控,可以為湖區(qū)水環(huán)境狀況的預測預報提供重要依據(jù)。本研究旨在通過研究鄱陽湖七大河口水樣中微生物的組成和多樣性,并研究不同河口水樣環(huán)境因子對微生物種類組成和多樣性的影響,以期為準確、快速地評價和預測鄱陽湖整體環(huán)境狀況提供重要依據(jù),為鄱陽湖生態(tài)功能的維護提供指導。

    1 材料與方法

    1.1 樣地選擇、樣品采集和處理

    本試驗選擇贛江(116°02′E,29°19′N)、撫河(116°09′E,28°46′N)、信江(116°43′E,28°72′N)、樂 安 河(116°69′E,28°97′N)、 昌 江(116°70′E,28°97′N)、饒河(116°46′E,29°01′N)和修河(116°02′E,29°19′N)七大河口進行研究,采樣時間為2014年10月底,處于幾大河口和鄱陽湖的枯水期,樣品采集點如圖1,采用GPS進行定位。采集水樣后,分成兩部分,一部分用于理化性質(zhì)測定,另一部分用于細菌基因組DNA提取和分析。

    圖1 鄱陽湖7個河口水樣取樣點Fig. 1 Collection sites of seven estuarine water samples of Poyang Lake

    1.2 水樣理化性質(zhì)測定

    水樣的溫度采用便攜式溫度計進行測定,pH值采用校準的pH計(Sartorius PB-10)測定。水樣的總氮和總磷采用過硫酸鉀氧化法測定[12]。氨態(tài)氮采用奈氏比色法測定,硝態(tài)氮和亞硝態(tài)氮含量分別采用紫外分光光度法和重氮偶合分光光度法測定[13]。全部測定均進行三次重復。

    1.3 樣品細菌收集及其基因組DNA提取

    將部分水樣先用8 μm孔徑的定性濾紙過濾以去除顆粒物以及真核生物,然后通過聚碳酸酯膜(0.2 μm,Whatman)真空抽濾[10],然后用無菌水沖洗2次,將濾膜取出剪成碎片,裝入1.5 ml的離心管中,-80℃保存?zhèn)溆谩?/p>

    樣品細菌基因組DNA的提取采用細菌基因組DNA提取試劑盒(天根生化科技有限公司)完成,具體步驟參照試劑盒說明書,1.2%瓊脂糖凝膠電泳檢測所提取的基因組DNA。Nanodrop 1000測定所提取DNA的濃度以及純度。DNA于-80 ℃保存?zhèn)溆?,符合要求的DNA樣本進行后續(xù)MiSeq測序。

    1.4 樣本PCR擴增

    以稀釋后的DNA為模板,采用帶有Barcode的通用引物515F和806R對細菌16S rRNA基因的V4區(qū)域進行擴增,引物序列為515F:5'-GTTTCGGTGC CAGCMGCCGCGGTAA-3',806R:5'-GTGAAAGG ACTACHVGGGTWTCTAAT -3'。PCR擴增條件為:98 ℃預變性1 min,98 ℃變性10 s,50 ℃退火30 s,72 ℃復性60 s,30個循環(huán),最后72 ℃延伸l0 min,4 ℃保存。每個樣品3個重復,將同一樣品的PCR產(chǎn)物混合后用2%瓊脂糖凝膠電泳檢測,使用北京天根生化科技有限公司瓊脂糖凝膠DNA回收試劑盒(DP209)切膠回收PCR產(chǎn)物,Tris-HCL洗脫。

    1.5 文庫構(gòu)建及測序

    測序文庫的構(gòu)建采用NEB Next? Ultra? DNA Library Prep kit for Illumina(NEB,美國)建庫試劑盒完成,具體步驟參照試劑盒說明書。構(gòu)建好的文庫經(jīng)過Qubit定量和文庫檢測,合格后,使用MiSeq進行上機測序。

    1.6 生物信息學分析

    采用FLASH軟件進行DNA片段的雙末端序列拼接[14]。序列分析采用UPARSE軟件包進行(Uparse v7.0.1001,http://drive5.com/uparse/)[15],算法采用UPARSE-OTU和UPARSE- OTUref。

    OTU(Operational Taxonomic Units,操作分類單元)聚類分析:在97%的相似水平下,應用QIIME軟件中Usearch對所有序列進行OTU劃分。具體分析方法:1)使用Uparse(version 7.1 http:// drive5.com/uparse/)方法進行OTU聚類,OTU中序列相似性設為97%,得到OTU的代表序列[16];2)使用Uchime(version 4.2.40 http://drive5.com/usearch/ manual/uchime_algo.html)檢測PCR擴增中產(chǎn)生的嵌合體序列并從OTU中去除[15];3)使用Usearch_ global方法將優(yōu)化序列map比對回OTU代表序列,得到OTU各樣品序列豐度統(tǒng)計表。

    分類學分析:為了得到每個OTU對應的物種分類信息,采用RDP classifer貝葉斯算法對97%相似水平的OTU代表序列進行分類學分析,選用Silva(Release115 http://www.arb-silva.de)數(shù)據(jù)庫進行比對。軟件及算法:Qiime平臺(http://qiime.org/ scripts/assign_taxonomy.html),RDP Classifer (version 2.2 http: //sourceforge.net/projects/rdp-classifer/),置信度閾值為0.7[17]?;诜诸悓W分析結(jié)果,在Excel中編輯作群落結(jié)構(gòu)組分圖。

    多樣性指數(shù)(Alpha-diversity)計算:在97%(0.97)的相似水平下,采用Mothur軟件(version v.1.30.1 http://www.mothur. org/wiki/Schloss_ SOP#Alpha_ diversity),計算菌群豐度指數(shù)Chao1以及菌群多樣性指數(shù)Shannon指數(shù)[18]。

    群落聚類分析:采用Qiime軟件計算群落間非加權的UniFrac距離,采用的是非加權配對算術平均法。

    1.7 統(tǒng)計分析

    不同河口水樣環(huán)境因子和微生物多樣性之間的相關性分析采用SPSS 16.0軟件進行。同時,采用Canoco軟件分析環(huán)境因子對不同類群細菌分布的影響。

    2 結(jié)果與討論

    2.1 鄱陽湖7個河口水樣理化性質(zhì)

    鄱陽湖不同河口水體溫度、pH、總氮、氨態(tài)氮、亞硝態(tài)氮、硝態(tài)氮、總磷含量如表1所示。從表1可以看出,樂安河水樣pH最低,偏酸性,為6.19,其余河口水樣pH均在6.41-6.86間波動,樂安河水樣偏酸性,可能與德興銅礦廢水的排放有關。從總氮和氨態(tài)氮含量結(jié)果可以看出,樂安河水樣中含量最高,分別為8.08 mg/L和6.06 mg/L,饒河含量次之,分別為5.74 mg/L和4.12 mg/L,污染較為嚴重。其原因可能是由于樂安河下游鄱陽縣大量的水產(chǎn)養(yǎng)殖以及生活污水的排放,從而導致水體中總氮和氨態(tài)氮含量偏高[19],且樂安河匯入饒河,其水質(zhì)在一定程度上會影響?zhàn)埡拥乃|(zhì)。信江水樣中總磷含量最高,為0.51 mg/L,這可能受到上饒市朝陽磷礦生產(chǎn)廢水的影響,從而導致信江下游總磷含量偏高[20]。

    表1 鄱陽湖7個河口水樣理化特性Table 1 Physicochemical properties of seven estuarine water columns of Poyang Lake and local meteorological conditions at their sampling sites

    2.2 鄱陽湖7個河口水樣中微生物多樣性和組成

    2.2.1 水樣樣本α多樣性 在所研究的7組樣本中,修河水樣的操作分類單元(OTU)數(shù)目最多(為851個),饒河及昌江樣本所含的操作分類單元數(shù)目次之,分別為825和772個,撫河水樣樣本所含的操作分類單元數(shù)目則相對較少,為698個(表2)。進一步通過種群豐富度指數(shù)(Chao1)對各樣本微生物豐度進行了評估,研究發(fā)現(xiàn)樂安河樣本種群豐富度指數(shù)最高(為1069.50),其次是修河(為998.02)、饒河(為982.13),信江樣本種群豐富度指數(shù)最低(為782.03)(表2)。該結(jié)果說明樂安河、修河和饒河河水微生物物種豐度較高,信江河水中微生物物種豐度低。相比較而言,樂安河、修河和饒河的水質(zhì)較其它河口差,氮含量相對較高,且三者微生物豐度較高,表明其水體中可能含有大量的參與氮循環(huán)的微生物。此外,樂安河匯入饒河,其水質(zhì)會對饒河的水質(zhì)造成一定影響。

    本研究同時通過香農(nóng)指數(shù)(Shannon index)對微生物的多樣性進行了評估,研究發(fā)現(xiàn)7大河口河水樣本均有著較高的生物多樣性,其中,修河河水樣本香農(nóng)指數(shù)(9.06)最高,其次為昌江(6.88)和饒河(6.79),信江樣本(6.32)香農(nóng)指數(shù)最?。ū?)。該實驗結(jié)果表明,修河、昌江、饒河水樣中細菌群落多樣性較其他河口高,信江水樣中細菌群落多樣性最低。

    表 2 河口水樣樣本α多樣性Table 2 Alpha diversity of all samples

    2.2.2 水樣細菌群落組成 進一步對各樣本細菌群落組成進行分析,由圖2可知,7大河口樣本中,變形菌門為主要的細菌類群(各樣本所占比例為62.65%-79.03%),放線菌門(2.11%-18.32%)和擬桿菌門(5.05%-9.08%)細菌所占比例次之。其中樂安河、信江和饒河樣本中放線菌門細菌含量要顯著低于其它河口。贛江樣本中放線菌門細菌含量(18.32%)最高。之前研究結(jié)果也表明,變形菌門、放線菌門和擬桿菌門為湖泊生態(tài)系統(tǒng)的優(yōu)勢菌門,它們參與了環(huán)境中污染物的生物降解[21-23]。該結(jié)果與其他研究人員對其他淡水湖泊水體細菌多樣性研究結(jié)果相似。吳蘭[10]通過對鄱陽湖水體微生物群落組成進行研究發(fā)現(xiàn),變形菌門為優(yōu)勢菌群,擬桿菌門和放線菌門有少量分布。Wu等[7]研究了太湖梅梁湖區(qū)的水體細菌群落遺傳組成,結(jié)果顯示,梅梁湖區(qū)水體細菌群落組成隨季節(jié)而變化,其主要以變形菌門、擬桿菌門和放線菌門為主。

    圖2 7個河口水樣細菌群落組成(門水平)Fig. 2 Bacterial composition of these communities in seven estuarine water samples(Phylum level)

    從綱水平來看,7大河口水樣中β-變形菌綱細菌占據(jù)絕對優(yōu)勢(33.8%-50.4%),α-變形菌綱(11.9%-33.3%)和放線菌綱(1.7%-16.2%)細菌所占比例次之。其中,樂安河水樣中β-變形菌綱細菌相對其他河口最高(50.4%),α-變形菌綱細菌含量則相對較低(16.4%),而撫河樣本中β-變形菌綱細菌含量最低(33.8%),α-變形菌綱細菌則含量最高(33.3%)(圖3)。根據(jù)以前關于β-變形菌綱的報道指出,在大部分淡水系統(tǒng)中,其都有著大量的分布,它們的新陳代謝活動在淡水湖泊生態(tài)系統(tǒng)中可能起重要作用,其參與淡水系統(tǒng)的氮代謝,具有降解復雜有機大分子的能力,其豐富度與低分子量的營養(yǎng)物含量相關。α-變形菌則主要存在于相對寡營養(yǎng)的環(huán)境中,擬桿菌的豐富度與水體的水華程度密切相關[24]。因此,該研究結(jié)果進一步揭示了樂安河水體的主要營養(yǎng)物可能以低分子量物質(zhì)居多,而撫河含沙量大,相對提供給微生物的營養(yǎng)偏少,α-變形菌綱細菌所占比例較高。此外,各大河口中高豐度的β-變形菌綱細菌可能在各自生態(tài)系統(tǒng)中發(fā)揮著重要的功能。

    除此之外,我們還發(fā)現(xiàn)贛江樣本中放線菌綱細菌含量(16.2%)要顯著高于其它樣本(圖3)。在淡水系統(tǒng)中放線菌綱細菌是一類獨立的微生物種群,水體中溶解有機碳含量與放線菌綱細菌的含量呈負相關關系[25]。但是我們的研究中由于缺少溶解有機碳含量的數(shù)據(jù),因此無法解釋贛江中放線菌綱細菌含量顯著高于其它幾個河口的原因。

    圖3 7個河口河水樣本細菌群落組成(綱水平)Fig. 3 Bacterial composition of these communities in seven estuarine water samples(Class level)

    2.3 細菌群落結(jié)構(gòu)和多樣性比較分析

    采用加權UniFrac距離來衡量兩兩樣本之間的相異系數(shù),值越小,表示這兩個樣品在物種多樣性方面存在的差異越小。由圖4可知,修河和贛江河水樣本之間的相異系數(shù)最小(0.136),樂安河和撫河樣本之間的相異系數(shù)最大(0.351),修河和撫河的相異系數(shù)(0.306)次之,該結(jié)果與河口水體環(huán)境因子聚類結(jié)果相一致,表明水體環(huán)境因子對于其中微生物多樣性存在影響。

    對這7個河口河水樣本細菌群落組成進行聚類分析,結(jié)果表明,贛江、修河和樂安河樣本的微生物組成相似,聚在一起,其中贛江和修河樣本微生物組成更相似,昌江和撫河樣本微生物群落組成相似,聚在一起,信江和饒河樣本細菌組成相似,兩者聚在一起(圖5)。

    圖4 Beta多樣性指數(shù)熱圖Fig. 4 Heatmap analysis based on the dissimilaritycoefficient between sets of two samples

    圖5 鄱陽湖7個河口河水樣本微生物群落組成聚類圖Fig. 5 UPGMA cluster analysis of these bacterial community compositions in seven estuarine water samples of Poyang Lake(Weighted UniFrac distance)

    2.4 環(huán)境因子對細菌群落組成和多樣性的影響

    采用SPSS 16.0軟件分析河口水樣微生物多樣性和水樣理化性質(zhì)及氣候因子的相關性。結(jié)果顯示,水樣中pH、總氮和氨態(tài)氮含量與河口水樣中細菌種群豐富度存在顯著的相關性(表3),其中pH和細菌群落豐度呈顯著的負相關性(rs=-0.786),總氮和氨態(tài)氮含量和細菌群落豐度呈顯著的正相關性(rs=0.786和0.829)。此外,我們也分析了不同河口水樣微生物群落組成與水樣理化性質(zhì)及氣候因子的相關性,研究結(jié)果表明,溫度、總氮以及亞硝態(tài)氮在一定程度上影響著環(huán)境中微生物的群落組成。放線菌門細菌與溫度呈顯著的正相關性(rs=0.883),浮霉菌門與氨態(tài)氮呈顯著的負相關性(rs=-0.757),藍細菌與總氮、氨態(tài)氮和亞硝態(tài)氮呈顯著的正相關性(rs=0.893,0.847和0.929),Thermi門細菌與溫度呈顯著的正相關性(rs=0.883),綠菌門細菌與亞硝態(tài)氮呈顯著的負相關性(rs=-0.829),綠彎菌門細菌與總氮和氨態(tài)氮呈顯著的負相關性(rs=-0.857,-0.829)(表3)。

    之前大量研究也表明,生境中環(huán)境因子顯著影響該生境中微生物種類組成和多樣性。Zhang等[26]在高原淡水湖微生物種類組成和多樣性的研究中也發(fā)現(xiàn)了類似的結(jié)果。生境中氮、氨態(tài)氮等含量可以直接影響微生物的種類組成,例如,藍細菌含量與環(huán)境中氮含量呈顯著的正相關性,藍細菌含量越高,環(huán)境富營養(yǎng)化程度越嚴重[27],此外,環(huán)境中的營養(yǎng)成分也可通過影響其它浮游生物從而間接影響微生物的組成[28],當環(huán)境中營養(yǎng)成分增加,浮游植物、藻類和細菌種類也會相應增加,他們之間存在著互利共生的關系[29]。我們的研究結(jié)果似乎支持了以前的研究報告,細菌群落組成和多樣性一定程度上取決于環(huán)境因素[30],此外也進一步證明了可以通過檢測特定功能性微生物的組成和多樣性從而間接的反映相應的環(huán)境變化。

    3 結(jié)論

    通過對鄱陽湖不同河口水體環(huán)境因子進行分析,發(fā)現(xiàn)樂安河水體中總氮、氨態(tài)氮含量均最高,其次為饒河,表明兩者污染較嚴重。從總磷含量來看,信江水體含量最高,其次為樂安河。

    7大河口水體微生物多樣性和豐度均有著一定的差異,樂安河、修河、昌江和饒河水體中微生物群落多樣性和豐度較其他河口高。水體細菌群落組成分析結(jié)果表明,變形菌門、放線菌門和擬桿菌門細菌為所有河口水體中的優(yōu)勢菌群。其中,β-變形菌綱、α-變形菌綱和放線菌綱細菌為所有河口水體中的優(yōu)勢菌綱。該研究結(jié)果表明這些細菌可能在河口水生態(tài)系統(tǒng)中發(fā)揮著重要的功能。此外,河口水體微生物種群豐富度指數(shù)和水體環(huán)境因子的相關性研究結(jié)果顯示,水體pH、總氮和氨態(tài)氮含量與細菌種群豐富度存在顯著的相關性,其中pH和細菌種群豐富度呈顯著的負相關性,總氮和氨態(tài)氮含量和細菌種群豐富度呈顯著的正相關性。該研究結(jié)果進一步證明了環(huán)境因子的改變是影響環(huán)境微生物群落組成的因素之一。

    [1] Wu L, Ge G, Zhu G F, et al. Diversity and composition of the bacterial community of Poyang Lake (China) as determined by 16S rRNA gene sequence analysis[J]. World Journal of Microbiology and Biotechnology, 2012, 28(1): 233-244.

    [2] Wu L, Li M, Guo Y, et al. Influence of three gorges project on water quality of Poyang Lake[J]. Procedia Environmental Sciences, 2011, 10(1): 1496-1501.

    [3] 劉聚濤, 游文蓀, 丁惠君. 鄱陽湖流域農(nóng)村水環(huán)境污染防治對策研究[J]. 江西水利科技, 2014, 40(2): 97-100. Liu J T, You W S, Ding H J. Study on the countermeasures of rural water environment pollution in Poyang Lake basin[J]. Jiangxi Hydraulic Science and Technology, 2014, 40(2): 97-100.

    [4] 胡悅之. 試論鄱陽湖流域水環(huán)境污染問題及其改善策略[J]. 環(huán)境與生活, 2014, 79(20): 11-12. Hu Y Z. Study on the water pollution problem and its improvement strategy of Poyang Lake[J]. Green living, 2014, 79(20): 11-12.

    [5] Zhao D Y, Huang R, Zeng J, et al. Diversity analysis of bacterial community compositions in sediments of urban lakes by terminal restriction fragment length polymorphism (T-RFLP)[J]. World Journal of Microbiology and Biotechnology, 2012, 28(11): 3159-3170.

    [6] Wu Q L, Zwart G, Schauer M, et al. Bacterioplankton community composition along a salinity gradient of sixteen high mountain Lakes located on the Tibetan Plateau[J]. Applied and Environmental Microbiology, 2006, 72(8): 5478-5485.

    [7] Wu L, Kellogg L, Devol A H, et al. Microarray-based characterization of microbial community functional structure and heterogeneity in marine sediments from the Gulf of Mexico[J]. Applied and Environmental Microbiology, 2008, 74(14): 4516-4529.

    [8] Tijdens M, Hoogveld H L, Kamst-van Agterveld M P, et al. Population dynamics and diversity of viruses, bacteria and phytoplankton in a shallow eutrophic lake[J]. Microbial Ecology,2008, 56(1): 29-42.

    [9] 余德輝. 我國環(huán)境科技發(fā)展現(xiàn)狀及機遇[J]. 中國環(huán)境管理,2001, 3(1): 4-6. Yu D H. Our environmental technology development status and opportunities[J]. Chinese Journal of Environmental Management,2001, 3(1): 4-6.

    [10] 吳蘭. 鄱陽湖水體細菌物種多樣性研究[D]. 南昌: 南昌大學,2009. Wu L. Studies on the bacterial species diversity in Lake Poyang[D]. Nanchang: Nanchang University, 2009.

    [11]張鈺. 湖泊沉積物微生物多樣性和厭氧氨氧化菌鑒別研究[D].武漢: 華中農(nóng)業(yè)大學, 2012. Zhang Y. Study on microbial diversity and anammox bacteria identification in Lake sediments[D]. Wuhan: Huazhong Agricultural University, 2012.

    [12] Song Y Y. Improvement in total nitrogen determination by using potassium persulfate oxidation method[J]. Environmental Science Survey, 2012, 3(31): 80-82.

    [13] American Public Health Association, American Water Works Association, Water Environment Federation. Standard Methods for the Examination of Water and Waste Water[M]. American Public Health Association, Washington, D C, 1992: 515-516.

    [14] Mago? T, Salzberg S L. FLASH: fast length adjustment of short reads to improve genome assemblies[J]. Bioinformatics, 2011,27(21): 2957-2963.

    [15] Edgar R C. UPARSE: highly accurate OTU sequences from microbial amplicon reads[J]. Nature Methods, 2013, 10(10): 996-998.

    [16] Edgar R C, Haas B J, Clemente J C, et al. UCHIME improves sensitivity and speed of chimera detection[J]. Bioinformatics,2011, 27(16): 2194-2200.

    [17] Quast C, Pruesse E, Yilmaz P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools[J]. Nucleic Acids Research, 2013, 41(D1): 590-596.

    [18] Schloss P D, Gevers D, Westcott S L. Reducing the effects of PCR amplifcation and sequencing artifacts on 16S rRNA-based studies[J]. PLoS One, 2011, 6(12): e27310.

    [19]劉聚濤, 鐘家有, 付敏, 等. 鄱陽湖流域農(nóng)村生活區(qū)面源污染特征及其影響[J]. 長江流域資源與環(huán)境, 2014, 23(7): 1012-1018. Liu J T, Zhong J Y, Fu M, et al. Study on the characteristics and impact of rural living area pollution in Poyang Lake basin[J]. Resources and Environment in the Yangtze Basin, 2014, 23(7): 1012-1018.

    [20] 張維球. 解決用朝陽磷礦濕法生產(chǎn)過磷酸鈣水分超標的問題[J].磷肥與復肥, 2000, 15(1): 22-23. Zhang W Q. Solving the problem of excessive moisture of superphosphate production using Chaoyang phosphate wet production method[J]. Phosphate and Compound Fertilizer, 2000,15(1): 22-23.

    [21] Song H, Li Z, Du B, et al. Bacterial communities in sediments of the shallow Lake Dongping in China[J]. Journal of Applied Microbiology, 2012, 112(1): 79-89.

    [22] Cupples A M. RDX degrading microbial communities and the prediction of microorganisms responsible for RDX bioremediation[J]. International Biodeterioration and Biodegradation, 2013, 85(7): 260-270.

    [23] Fuentes S, Mendez V, Aguila P, et al. Bioremediation of petroleum hydrocarbons: catabolic genes microbial communities and applications[J]. Applied Microbiology and Biotechnology, 2014,98(11): 4781-4794.

    [24] Kirchman D L. The ecology of Cytophaga-Flavobacteria in aquatic environments[J]. FEMS Microbiology Ecology, 2002, 39(2): 91-100.

    [25] Burkert U, Warnecke F, Babenzien D, et al. Members of a readily enriched beta-proteobacterial clade are common in surface waters of a humic lake[J]. Applied and Environmental Microbiology, 2003, 69(11): 6550-6559.

    [26] Zhang J X, Yang Y Y, Zhao L, et al. Distribution of sediment bacterial and archaeal communities in plateau freshwater lakes[J]. Applied Microbiology and Biotechnology, 2015, 99(7): 3291-3302.

    [27] Reynolds C S. Cyanobacterial water blooms[J]. Advance in Botanical Research, 1987, 13(4): 67-143.

    [28] Hietala J, Vakkilainen K, Kairesalo T. Community resistance and change to nutrient enrichment and fsh manipulation in a vegetated lake littoral[J]. Freshwater Biology, 2004, 49(12): 1525-1537.

    [29] Hietala J, Vakkilainen K, Kairesalo T. Community resistance and change to nutrient enrichment and fsh manipulation in a vegetated lake littoral[J]. Freshwater Biology, 2004, 49(12): 1525-1537.

    [30] Yannarell A C, Triplett, E W. Geographic and environmental sources of variation in lake bacterial community composition[J]. Applied and Environmental Microbiology, 2005, 71(1): 227-239.

    (責任編輯:王育花)

    Bacterial diversities and compositions in seven different estuarine water columns of Poyang Lake

    SHENG Ping1, YU Yi-zun1, TIAN Xiao-juan1, HUANG Jiang-li1, ZHANG Zhi-hong1, 2, DING Jian-nan1*
    (1. Institute of Biological Resources, Jiangxi Academy of Sciences, Nanchang, Jiangxi 330096, China; 2. Poyang Lake Key Laboratory of Jiangxi, Nanchang, Jiangxi 330096, China)

    In this study, we analyzed environmental factors of seven different estuarine water columns of Poyang Lake (Gan, Fu, Xin, Le'an, Chang, Rao, and Xiu River), and the bacterial community compositions in these different estuarine sediments were analyzed using the 16S rRNA gene targeted metagenomic approach. Our results showed that Le'an River had the highest total nitrogen content, and Xin River had the highest total phosphorus content. For the bacterial richness and diversity in these seven different estuaries, we found that Le'an, Xiu and Rao rivers showed the higher bacterial richness than other estuarines. Xiu, Chang and Rao rivers showed higher bacterial diversities than others. Bacterial composition analysis indicated that Proteobacteria phylum was dominant in all water samples, followed by Actinobacteria and Bacteroidetes. At the class level, all samples were numerically dominated by members of the Betaproteobacteria, Alphaproteobacteria and Actinobacteria. Furthermore, according to the Spearman's correlation analysis, we found that there were signifcant negative relationships between the Chao1 species richness index and pH. Besides, total nitrogen and ammonia concentrations showed signifcant positive effect on bacterial community richness. Key words: Poyang Lake; seven estuaries; different estuarine water columns; bacterial diversity; bacterial community composition

    鄱陽湖是我國最大的淡水湖,地處長江中下游南岸。其上承贛、撫、信、饒、修五河之水,下接長江,其中樂安河和昌江匯入饒河,流入鄱陽湖。鄱陽湖濕地生態(tài)環(huán)境好,是世界上最大的鳥類保護區(qū)[1]。但是,近年來,由于工業(yè)、農(nóng)業(yè)和城市污水的大量排放,加上降水逐年減少,導致入湖河流以及湖泊的自凈能力下降,使得湖泊氮磷濃度持續(xù)升高、湖泊富營養(yǎng)化程度不斷加深[2-4],對湖泊周邊環(huán)境造成嚴重的威脅。

    在湖泊生態(tài)系統(tǒng)中,微生物在湖泊物質(zhì)循環(huán)和能量循環(huán)中起著重要的作用,同時也是各種新型生物活性物質(zhì)的潛在來源[5]。由于湖泊水體生境類型不同,使得生境中微生物的多樣性也不同。研究表明,湖泊鹽度、湖水深度、湖水和底泥pH、營養(yǎng)成分、有機質(zhì)等均會影響生境中微生物的組成[6]。因此,由于自然和人類活動引起自然環(huán)境的改變最終也會引起環(huán)境微生物種類組成的改變[7-8]。我國是個湖泊眾多的國家,約有兩萬多個湖泊,其中約三分之一為淡水湖泊,其水資源占中國城鎮(zhèn)飲用水的50%以上[9-10]。但是隨著資源的破壞性開發(fā),環(huán)境的惡化,世界各地湖泊富營養(yǎng)化進程加速,水質(zhì)日趨惡化[11],湖泊生物資源包括湖泊微生物種類與多樣性研究己成為各國科學家共同關注的熱點。

    The special funds of National Academy Alliance for Academy Construction (2012-1), Science and Technology Plan of Jiangxi Province (2008BB26100), The Special Funds of Jiangxi Academy of Sciences (2013-19-01).

    DING Jian-nan, E-mail: jiannanding@aliyun.com.

    26 June, 2015;Accepted 24 February, 2016.

    X172

    A

    1000-0275(2016)03-0606-07

    10.13872/j.1000-0275.2016.0050

    全國科學院聯(lián)盟專項重點項目(2012-1);江西省科技支撐計劃項目(2008BB26100);江西省科學院專項項目(2013-19-01)。

    圣平(1987-),女,湖北荊州人,博士,助理研究員,主要從事微生物學與分析生物學研究,E-mail:shengping_1014@163. com;通訊作者:丁建南(1955-),男,湖南衡陽人,博士,研究員,主要從事微生物生態(tài)學研究,E-mail: jiannanding@aliyun.com。

    2015-06-26,接受日期:2016-02-24

    猜你喜歡
    饒河樂安鄱陽湖
    ABSTRACTS
    鄱陽湖學刊(2024年3期)2024-01-01 00:00:00
    鄱陽湖水系之潦河
    空間與權力:清代鄱陽湖區(qū)饒河民間信仰神靈水神化芻議
    江漢論壇(2020年7期)2020-08-09 08:48:40
    鄱陽湖流域饒河水生態(tài)環(huán)境保護與治理對策研究
    《鄱陽湖生態(tài)系列插畫》
    江西省樂安老年大學校歌
    烏蘇里江流域三面紅旗的誕生地在饒河
    世紀橋(2017年10期)2017-11-22 17:41:59
    詩二首
    中華魂(2017年3期)2017-11-22 02:01:55
    民族英雄趙一曼之歌
    中華魂(2017年5期)2017-11-21 01:53:23
    來自銀時空的紙條
    观看免费一级毛片| 久久久久久久久久黄片| 国产又色又爽无遮挡免费看| 两个人看的免费小视频| 狂野欧美白嫩少妇大欣赏| 午夜福利免费观看在线| 日韩精品免费视频一区二区三区| 日韩免费av在线播放| 欧美激情久久久久久爽电影| 欧美久久黑人一区二区| 国产区一区二久久| 听说在线观看完整版免费高清| 久久久国产欧美日韩av| 久久精品91蜜桃| 男女那种视频在线观看| 免费在线观看成人毛片| 99久久久亚洲精品蜜臀av| 国产av一区在线观看免费| 看黄色毛片网站| 欧美成狂野欧美在线观看| 国产亚洲欧美在线一区二区| av在线播放免费不卡| 国内少妇人妻偷人精品xxx网站 | 欧美高清成人免费视频www| 18禁裸乳无遮挡免费网站照片| or卡值多少钱| 成熟少妇高潮喷水视频| 色播亚洲综合网| 亚洲男人天堂网一区| 搡老岳熟女国产| 色哟哟哟哟哟哟| 夜夜看夜夜爽夜夜摸| 搡老熟女国产l中国老女人| 欧美大码av| 又爽又黄无遮挡网站| 午夜老司机福利片| 1024视频免费在线观看| 在线视频色国产色| 国内久久婷婷六月综合欲色啪| 51午夜福利影视在线观看| 女人高潮潮喷娇喘18禁视频| 国产精品一区二区精品视频观看| 久久久久免费精品人妻一区二区| 亚洲av成人精品一区久久| 亚洲人与动物交配视频| 国产精品美女特级片免费视频播放器 | 久久午夜综合久久蜜桃| 亚洲精品美女久久久久99蜜臀| 亚洲成人久久性| 亚洲中文字幕日韩| 又黄又爽又免费观看的视频| 久久精品亚洲精品国产色婷小说| 精品久久久久久久人妻蜜臀av| 美女扒开内裤让男人捅视频| 中出人妻视频一区二区| 最新在线观看一区二区三区| 夜夜夜夜夜久久久久| 久久精品91无色码中文字幕| 国产91精品成人一区二区三区| 一区福利在线观看| 男女下面进入的视频免费午夜| 久久久国产欧美日韩av| 嫩草影视91久久| 午夜视频精品福利| 欧美一区二区精品小视频在线| 一边摸一边做爽爽视频免费| 久久香蕉精品热| 亚洲专区国产一区二区| 亚洲午夜精品一区,二区,三区| 超碰成人久久| 色老头精品视频在线观看| 国产成人精品久久二区二区免费| 欧美高清成人免费视频www| 国产99久久九九免费精品| 国产精品98久久久久久宅男小说| 国产高清videossex| 久久久久精品国产欧美久久久| 久久精品国产综合久久久| 全区人妻精品视频| 国产男靠女视频免费网站| 美女免费视频网站| 亚洲av成人av| videosex国产| 俄罗斯特黄特色一大片| 91国产中文字幕| 久久久久久亚洲精品国产蜜桃av| 色av中文字幕| 亚洲国产精品久久男人天堂| 一区二区三区国产精品乱码| 成人高潮视频无遮挡免费网站| 亚洲中文字幕一区二区三区有码在线看 | 国产av一区在线观看免费| 99久久99久久久精品蜜桃| 国产精品久久久久久久电影 | 黑人欧美特级aaaaaa片| 亚洲午夜理论影院| 好男人电影高清在线观看| 日韩成人在线观看一区二区三区| 国产精品av久久久久免费| 免费看日本二区| 老司机深夜福利视频在线观看| 18禁观看日本| 宅男免费午夜| 在线观看免费午夜福利视频| 国产人伦9x9x在线观看| 亚洲午夜精品一区,二区,三区| 两个人免费观看高清视频| 日韩大码丰满熟妇| 黄色毛片三级朝国网站| 99re在线观看精品视频| 在线观看免费日韩欧美大片| 日韩精品免费视频一区二区三区| 国产成人精品无人区| 日韩欧美在线二视频| 熟女电影av网| a级毛片a级免费在线| 欧美3d第一页| 国产精品98久久久久久宅男小说| 欧美黑人欧美精品刺激| 亚洲人成伊人成综合网2020| 国语自产精品视频在线第100页| 99久久精品热视频| 香蕉国产在线看| 精品国内亚洲2022精品成人| 亚洲精品在线观看二区| 69av精品久久久久久| 欧美日韩亚洲国产一区二区在线观看| 999久久久精品免费观看国产| 天堂动漫精品| 亚洲五月婷婷丁香| 99国产精品99久久久久| 欧美黑人精品巨大| 亚洲成av人片免费观看| 看片在线看免费视频| av免费在线观看网站| 国产精品亚洲一级av第二区| 久久这里只有精品19| 亚洲成a人片在线一区二区| 一级毛片精品| 亚洲天堂国产精品一区在线| 亚洲av中文字字幕乱码综合| 久久国产精品影院| 国产精品一区二区精品视频观看| 日韩欧美 国产精品| 欧美黑人精品巨大| 制服丝袜大香蕉在线| 99久久国产精品久久久| 久久国产精品人妻蜜桃| 亚洲va日本ⅴa欧美va伊人久久| 两个人的视频大全免费| 人人妻,人人澡人人爽秒播| 国产成人系列免费观看| 国产av一区二区精品久久| 色综合亚洲欧美另类图片| 一个人免费在线观看的高清视频| a级毛片在线看网站| 午夜老司机福利片| 一级毛片女人18水好多| 免费无遮挡裸体视频| 日本五十路高清| 99热只有精品国产| 亚洲国产精品成人综合色| 真人一进一出gif抽搐免费| 美女扒开内裤让男人捅视频| 国产成+人综合+亚洲专区| 国产精品久久久久久久电影 | 国产精品一区二区精品视频观看| 一个人免费在线观看的高清视频| 级片在线观看| 国产av在哪里看| 岛国在线免费视频观看| 国产精品野战在线观看| 久久精品国产亚洲av香蕉五月| 欧美乱妇无乱码| 久久久久久久久中文| 欧美av亚洲av综合av国产av| 久久精品成人免费网站| 久久久精品大字幕| www日本黄色视频网| 99riav亚洲国产免费| 丁香六月欧美| 久久99热这里只有精品18| 国产av又大| 国产午夜精品论理片| 午夜a级毛片| 一本久久中文字幕| 久久久久久久久中文| 国内少妇人妻偷人精品xxx网站 | 欧美性猛交黑人性爽| 国产成人精品久久二区二区91| 国产69精品久久久久777片 | 中文资源天堂在线| 国产免费av片在线观看野外av| 色精品久久人妻99蜜桃| 欧美日韩亚洲国产一区二区在线观看| 在线国产一区二区在线| 久久久久久九九精品二区国产 | 97超级碰碰碰精品色视频在线观看| 国产aⅴ精品一区二区三区波| 不卡av一区二区三区| 亚洲精品在线观看二区| 麻豆国产97在线/欧美 | 亚洲精品av麻豆狂野| 三级国产精品欧美在线观看 | 国产一区二区激情短视频| 91国产中文字幕| 99热这里只有精品一区 | 国产精品久久电影中文字幕| 欧美乱码精品一区二区三区| 99国产极品粉嫩在线观看| 亚洲男人天堂网一区| 国产精品久久视频播放| 亚洲免费av在线视频| 亚洲中文日韩欧美视频| 每晚都被弄得嗷嗷叫到高潮| 美女高潮喷水抽搐中文字幕| 成人午夜高清在线视频| 中文字幕精品亚洲无线码一区| 精品国内亚洲2022精品成人| 小说图片视频综合网站| 国产高清激情床上av| 久久人妻福利社区极品人妻图片| 操出白浆在线播放| 黑人欧美特级aaaaaa片| 久久国产精品影院| 波多野结衣高清作品| 久久久久久亚洲精品国产蜜桃av| 悠悠久久av| 欧美在线一区亚洲| 日日摸夜夜添夜夜添小说| avwww免费| 国产真实乱freesex| 国产精品乱码一区二三区的特点| 女生性感内裤真人,穿戴方法视频| 亚洲欧美精品综合久久99| 国产高清视频在线播放一区| 99riav亚洲国产免费| 久久久久免费精品人妻一区二区| 国产欧美日韩精品亚洲av| 欧美极品一区二区三区四区| 午夜免费成人在线视频| 制服诱惑二区| 国产精品98久久久久久宅男小说| 亚洲av电影在线进入| 国产v大片淫在线免费观看| 亚洲中文字幕一区二区三区有码在线看 | 露出奶头的视频| 国内久久婷婷六月综合欲色啪| 天天躁夜夜躁狠狠躁躁| 久久精品成人免费网站| 国产精品国产高清国产av| 黄色视频不卡| 制服人妻中文乱码| 男女下面进入的视频免费午夜| av福利片在线| 日韩中文字幕欧美一区二区| 久久久久久久精品吃奶| 免费在线观看视频国产中文字幕亚洲| 国产精品九九99| 成熟少妇高潮喷水视频| 人人妻人人澡欧美一区二区| 国产真实乱freesex| 婷婷亚洲欧美| 久久人妻福利社区极品人妻图片| 国产成人精品久久二区二区免费| 国产99白浆流出| 狠狠狠狠99中文字幕| 最新美女视频免费是黄的| 色在线成人网| 97超级碰碰碰精品色视频在线观看| 白带黄色成豆腐渣| 国产精品爽爽va在线观看网站| 精品久久久久久久久久久久久| 国产男靠女视频免费网站| 真人一进一出gif抽搐免费| 色老头精品视频在线观看| 国产av不卡久久| 一级片免费观看大全| 又大又爽又粗| 国产视频内射| 欧洲精品卡2卡3卡4卡5卡区| 99精品久久久久人妻精品| 精品午夜福利视频在线观看一区| 99久久精品国产亚洲精品| 欧洲精品卡2卡3卡4卡5卡区| 成人国产综合亚洲| 欧美日本亚洲视频在线播放| 久久精品夜夜夜夜夜久久蜜豆 | 麻豆国产97在线/欧美 | 久久九九热精品免费| 狂野欧美激情性xxxx| 免费看日本二区| 欧美 亚洲 国产 日韩一| av中文乱码字幕在线| 免费在线观看完整版高清| 人人妻,人人澡人人爽秒播| 国产又色又爽无遮挡免费看| 午夜免费观看网址| 欧美色视频一区免费| 久久国产精品人妻蜜桃| 欧美高清成人免费视频www| 国产又色又爽无遮挡免费看| 久久性视频一级片| 亚洲午夜理论影院| 变态另类成人亚洲欧美熟女| 波多野结衣高清作品| 国产麻豆成人av免费视频| 欧美 亚洲 国产 日韩一| 天天躁夜夜躁狠狠躁躁| 在线免费观看的www视频| 亚洲五月天丁香| 男女那种视频在线观看| 妹子高潮喷水视频| 亚洲国产欧美人成| 久久精品国产99精品国产亚洲性色| 黄色毛片三级朝国网站| 婷婷六月久久综合丁香| 美女 人体艺术 gogo| 国产精品国产高清国产av| 男人舔女人下体高潮全视频| 91av网站免费观看| 熟女少妇亚洲综合色aaa.| 免费在线观看视频国产中文字幕亚洲| 男女做爰动态图高潮gif福利片| 毛片女人毛片| 欧美高清成人免费视频www| 亚洲,欧美精品.| 两人在一起打扑克的视频| 欧美不卡视频在线免费观看 | 国产久久久一区二区三区| 精品一区二区三区视频在线观看免费| 日韩成人在线观看一区二区三区| 两个人视频免费观看高清| 亚洲男人的天堂狠狠| 日韩欧美 国产精品| 亚洲一码二码三码区别大吗| 美女 人体艺术 gogo| 高清在线国产一区| 国产黄片美女视频| 久久这里只有精品中国| 亚洲五月婷婷丁香| 99热这里只有精品一区 | 欧美乱色亚洲激情| 亚洲av美国av| 国产av又大| 欧美精品啪啪一区二区三区| 夜夜看夜夜爽夜夜摸| 久久精品国产清高在天天线| 制服丝袜大香蕉在线| 1024手机看黄色片| 别揉我奶头~嗯~啊~动态视频| 99久久精品国产亚洲精品| 两个人看的免费小视频| 成人精品一区二区免费| 成人精品一区二区免费| 午夜福利高清视频| 亚洲自拍偷在线| 在线免费观看的www视频| 90打野战视频偷拍视频| 亚洲精品久久国产高清桃花| 国产成人影院久久av| 精品熟女少妇八av免费久了| 很黄的视频免费| 正在播放国产对白刺激| 两性午夜刺激爽爽歪歪视频在线观看 | 国内少妇人妻偷人精品xxx网站 | 久久国产精品影院| 日本免费一区二区三区高清不卡| 每晚都被弄得嗷嗷叫到高潮| 搡老岳熟女国产| 国产精品一区二区三区四区久久| 变态另类丝袜制服| 成在线人永久免费视频| 美女 人体艺术 gogo| 人妻丰满熟妇av一区二区三区| 身体一侧抽搐| 久久人妻福利社区极品人妻图片| www国产在线视频色| 白带黄色成豆腐渣| 大型黄色视频在线免费观看| 欧美黄色片欧美黄色片| 日韩高清综合在线| 欧美中文综合在线视频| 熟妇人妻久久中文字幕3abv| 国产精品影院久久| 18禁观看日本| 深夜精品福利| 亚洲一区高清亚洲精品| 成人国产一区最新在线观看| 午夜精品在线福利| 欧美日韩精品网址| 美女午夜性视频免费| 国产成人av教育| 国产精品国产高清国产av| 国产91精品成人一区二区三区| 免费在线观看黄色视频的| 一级毛片高清免费大全| 在线观看一区二区三区| 久9热在线精品视频| 国产精品国产高清国产av| 女生性感内裤真人,穿戴方法视频| 欧美日韩亚洲国产一区二区在线观看| 国产高清视频在线播放一区| 99久久无色码亚洲精品果冻| 丁香六月欧美| 丰满的人妻完整版| 日本撒尿小便嘘嘘汇集6| 欧美成人午夜精品| 日本免费a在线| 高潮久久久久久久久久久不卡| 成人三级黄色视频| 香蕉国产在线看| 曰老女人黄片| 91成年电影在线观看| 黑人欧美特级aaaaaa片| 宅男免费午夜| 日日夜夜操网爽| 91成年电影在线观看| 亚洲一区二区三区不卡视频| 精品免费久久久久久久清纯| 老熟妇乱子伦视频在线观看| 91国产中文字幕| 12—13女人毛片做爰片一| 无人区码免费观看不卡| 亚洲乱码一区二区免费版| 亚洲第一欧美日韩一区二区三区| 国产免费av片在线观看野外av| 国产又黄又爽又无遮挡在线| 女生性感内裤真人,穿戴方法视频| 欧美+亚洲+日韩+国产| 桃红色精品国产亚洲av| 精品电影一区二区在线| 国产69精品久久久久777片 | 制服丝袜大香蕉在线| 国内精品一区二区在线观看| 母亲3免费完整高清在线观看| 一级黄色大片毛片| 一级作爱视频免费观看| 人人妻,人人澡人人爽秒播| 热99re8久久精品国产| 2021天堂中文幕一二区在线观| 精品欧美一区二区三区在线| 亚洲激情在线av| 免费在线观看日本一区| 露出奶头的视频| 99国产精品一区二区三区| 欧美黑人巨大hd| 动漫黄色视频在线观看| 国产成人av教育| 亚洲美女黄片视频| 国产精品久久电影中文字幕| 男人舔女人下体高潮全视频| 亚洲人成网站高清观看| 嫩草影院精品99| 天天一区二区日本电影三级| 国产精品 国内视频| 精品一区二区三区视频在线观看免费| 午夜激情福利司机影院| 国产又色又爽无遮挡免费看| 亚洲人成网站高清观看| 色综合站精品国产| 欧美激情久久久久久爽电影| 给我免费播放毛片高清在线观看| 亚洲专区中文字幕在线| 99久久综合精品五月天人人| 丰满人妻一区二区三区视频av | 成人av在线播放网站| 男人舔奶头视频| 精品日产1卡2卡| 国产私拍福利视频在线观看| 中文字幕高清在线视频| cao死你这个sao货| 亚洲精品中文字幕在线视频| 成年免费大片在线观看| 丰满的人妻完整版| 国产精品野战在线观看| 最近最新中文字幕大全免费视频| 国产视频一区二区在线看| 在线观看免费午夜福利视频| 国产精品久久久av美女十八| 最新美女视频免费是黄的| 亚洲中文av在线| 欧美性长视频在线观看| 久久性视频一级片| 国产男靠女视频免费网站| 法律面前人人平等表现在哪些方面| 99精品欧美一区二区三区四区| 免费一级毛片在线播放高清视频| 国产高清有码在线观看视频 | 色播亚洲综合网| 不卡一级毛片| 欧美中文日本在线观看视频| av欧美777| www.熟女人妻精品国产| 久久精品综合一区二区三区| 韩国av一区二区三区四区| 一级片免费观看大全| 两个人的视频大全免费| 国产v大片淫在线免费观看| 美女黄网站色视频| 国产亚洲精品av在线| 亚洲av第一区精品v没综合| 亚洲成人免费电影在线观看| 国产日本99.免费观看| 欧美又色又爽又黄视频| 狂野欧美激情性xxxx| 在线观看一区二区三区| 日本三级黄在线观看| 在线永久观看黄色视频| 级片在线观看| 亚洲激情在线av| 成人18禁高潮啪啪吃奶动态图| 在线永久观看黄色视频| √禁漫天堂资源中文www| 一级a爱片免费观看的视频| 桃色一区二区三区在线观看| 成人国产一区最新在线观看| 免费在线观看黄色视频的| 一进一出抽搐gif免费好疼| 最近在线观看免费完整版| 全区人妻精品视频| 国产亚洲av嫩草精品影院| 国产99白浆流出| 精品欧美国产一区二区三| 麻豆久久精品国产亚洲av| 国产亚洲av嫩草精品影院| 久久香蕉精品热| 国产亚洲av嫩草精品影院| 999久久久精品免费观看国产| 午夜精品在线福利| 两个人的视频大全免费| 久99久视频精品免费| 免费在线观看影片大全网站| 别揉我奶头~嗯~啊~动态视频| 日韩免费av在线播放| 午夜老司机福利片| 最好的美女福利视频网| 可以免费在线观看a视频的电影网站| 97人妻精品一区二区三区麻豆| 亚洲av电影不卡..在线观看| 他把我摸到了高潮在线观看| 亚洲天堂国产精品一区在线| 色综合站精品国产| 国产精品98久久久久久宅男小说| 热99re8久久精品国产| 黄色女人牲交| 亚洲国产日韩欧美精品在线观看 | 久久中文字幕一级| 欧美成人午夜精品| 亚洲一区高清亚洲精品| 日本精品一区二区三区蜜桃| 亚洲最大成人中文| 欧美日韩一级在线毛片| 国产精品av视频在线免费观看| 成人精品一区二区免费| 中国美女看黄片| 嫩草影视91久久| 久久久久精品国产欧美久久久| 亚洲狠狠婷婷综合久久图片| 一区二区三区激情视频| 免费看日本二区| 哪里可以看免费的av片| 亚洲午夜精品一区,二区,三区| 在线观看免费午夜福利视频| 久久婷婷人人爽人人干人人爱| 真人做人爱边吃奶动态| 中国美女看黄片| 欧美日本亚洲视频在线播放| 久久久精品大字幕| 亚洲国产精品成人综合色| 成人国产一区最新在线观看| 黄色片一级片一级黄色片| 男人舔女人下体高潮全视频| 黄色毛片三级朝国网站| 亚洲五月天丁香| 欧美+亚洲+日韩+国产| 国产av一区在线观看免费| 久久久精品欧美日韩精品| 国产精品亚洲一级av第二区| aaaaa片日本免费| 色综合亚洲欧美另类图片| 亚洲国产日韩欧美精品在线观看 | 久久中文字幕一级| 亚洲色图av天堂| 非洲黑人性xxxx精品又粗又长| 草草在线视频免费看| 精品国内亚洲2022精品成人| 免费搜索国产男女视频| 日韩欧美国产一区二区入口| 国产aⅴ精品一区二区三区波| 国产精品精品国产色婷婷| 成年版毛片免费区| 99热6这里只有精品| 91麻豆精品激情在线观看国产| 中文字幕av在线有码专区| 一边摸一边抽搐一进一小说| 久久久国产欧美日韩av| 一级毛片女人18水好多| 免费搜索国产男女视频| 国产三级中文精品| 九色国产91popny在线| 狠狠狠狠99中文字幕| 午夜福利高清视频| 成人18禁在线播放| 国产精品九九99| 成年版毛片免费区| e午夜精品久久久久久久| 免费高清视频大片| 亚洲色图av天堂| 久久这里只有精品19| 欧美黑人欧美精品刺激| 欧美日韩精品网址| 亚洲成人久久爱视频| 久久人人精品亚洲av| 91老司机精品|