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

    安徽升金湖越冬白頭鶴(Grusmonacha)腸道菌群組成的時間變化*

    2017-05-17 03:07:05趙廣紅周立志董元秋宋昀微
    湖泊科學 2017年3期
    關鍵詞:菌門越冬菌群

    趙廣紅,周立志**,董元秋,宋昀微,2

    (1:安徽大學資源與環(huán)境工程學院,合肥 230601) (2:安徽升金湖國家級自然保護區(qū)管理局,池州 247200)

    安徽升金湖越冬白頭鶴(Grusmonacha)腸道菌群組成的時間變化*

    趙廣紅1,周立志1**,董元秋1,宋昀微1,2

    (1:安徽大學資源與環(huán)境工程學院,合肥 230601) (2:安徽升金湖國家級自然保護區(qū)管理局,池州 247200)

    腸道微生物對于維持動物消化道內(nèi)環(huán)境的穩(wěn)定具有重要意義,而腸道微生物菌群結構也受寄主食物成份等外界環(huán)境因素的影響. 研究動物腸道微生物菌群結構可以更好地了解動物生活狀況. 以升金湖越冬白頭鶴的糞便作為研究對象,以微生物16S rRNA的 V3~V4區(qū)作為標記基因,進行高通量測序,分析升金湖越冬不同時期白頭鶴腸道菌群結構的差異性. 研究結果表明:升金湖越冬期不同時期白頭鶴糞便樣品中微生物多樣性指數(shù)如香農(nóng)威納多樣性指數(shù)和辛普森指數(shù)無顯著性差異. 但是越冬期不同階段的糞便樣品中菌群組成存在差異,主要包括:在門分類水平,越冬前期擬桿菌門含量較高,厚壁菌門含量較中后期低;越冬中期放線菌門含量較高;越冬后期變形菌門含量較高. 在屬分類水平,相對于其他兩個時期,越冬前期普氏菌屬(Prevotella)含量較高,越冬中期桿菌屬(Lysinibacillus)含量較高,越冬后期梭菌屬(Clostridium)含量較高,并且普氏菌屬含量從越冬前期到中期、再到后期有明顯的下降趨勢,而梭菌屬含量從越冬前期到中期、再到后期有明顯的升高趨勢. 這些結果表明不同越冬時期白頭鶴腸道菌群組成不同,越冬期不同階段的食物變化可能作為一個影響因素.

    白頭鶴;腸道菌群組成;越冬期;16S rRNA;升金湖

    動物腸道微生物群落具有豐富的種類和功能多樣性[1],這些微生物之間及微生物和宿主之間存在穩(wěn)定復雜的關系. 許多情況下,動物及其腸道微生物之間互利共生,協(xié)同進化[2]. 隨著16S rRNA作為微生物標記基因的應用以及高通量等分子生物學技術的發(fā)展,加之糞便取樣手段的引入,腸道微生物,特別是腸道微生物對宿主生理功能的影響,已經(jīng)成為近年來研究的熱點問題[3]. 動物胃腸道菌群對宿主的營養(yǎng)、消化和吸收方面具有重要作用,腸道微生物含有很多分解酶可以降解食物中不易消化的成分,如多糖類纖維素類物質(zhì)[4],并且可以合成氨基酸、脂肪酸及維生素等宿主必需物質(zhì)[5]. 同時腸道微生物在宿主免疫功能上也有重要作用[6].

    腸道微生物在影響宿主動物生理功能的同時,其菌群結構也受到宿主生存環(huán)境的影響[7],特別是宿主食性不僅能夠影響,而且在一定程度上能夠決定腸道微生物菌群組成[8]. 有研究表明,動物腸道微生物數(shù)目由肉食性到雜食性再到草食性呈逐步增加趨勢[9]. 在對人的腸道微生物研究時,由于飲食習慣的不同,不同地區(qū)生活的居民的腸道微生物菌落結構具有很大的差異,而相同地區(qū)生活的居民的腸道微生物菌落結構相似性較高[10-12]. 腸道微生物菌群結構對于動物的食性具有很強的適應性,大熊貓(Ailuropodamelanoleuca)以竹子為主要食性,含有豐富的纖維素,在長期的進化過程中,在其腸道微生物中含有豐富的能夠分解纖維素的梭菌屬(Clostridium)[13].

    白頭鶴(Grusmonacha) 屬于鶴形目,鶴科,為大型遷徙水鳥,屬于國家I級重點保護野生動物,世界自然保護聯(lián)盟(IUCN)紅色物種名錄中的易危物種[14],長江中下游的升金湖是白頭鶴的重要越冬地,每年在此越冬的白頭鶴大約在300只左右[15]. 升金湖作為淺水內(nèi)陸湖泊,含有豐富的生物資源,特別是每年冬季,大量的灘涂暴露,為越冬水鳥提供豐富的食物資源以及覓食棲息地[16]. 近幾年,由于湖泊漁業(yè)養(yǎng)殖、圍墾造田、水利建設、旅游開發(fā)等人為因素的干擾,白頭鶴的棲息地不斷退化,適宜的覓食生境不斷喪失[17]. 升金湖越冬白頭鶴覓食的主要生境為稻田和草灘[18],但在不同的越冬階段,它們的相對重要性發(fā)生變化[19]. 白頭鶴在越冬前期選擇稻田作為主要覓食生境,以稻谷作為主要食物;越冬中期,白頭鶴對稻谷的利用率相對前期降低,但仍以稻谷作為最主要的食物資源;越冬后期,白頭鶴逐漸選擇草灘生境作為其主要覓食生境,以朝天委陵菜(Polygonumcriopolitanu)作為主要食物資源[20]. 因此,不同越冬階段白頭鶴腸道微生物菌群組成的研究對于了解食物資源變化對腸道微生物菌群結構的影響具有重要意義.

    本文通過16S rRNA 高通量測序技術,研究升金湖越冬白頭鶴不同時期腸道菌群組成情況,揭示越冬不同時期白頭鶴腸道菌群組成的差異,結合食物資源的變化,探討越冬白頭鶴腸道微生物與外界環(huán)境之間的關系,為白頭鶴越冬生理生態(tài)和種群管理提供基礎資料.

    1 材料和方法

    1.1 樣品的采集與保存

    根據(jù)白頭鶴的遷徙規(guī)律和越冬地的自然氣候特征,將白頭鶴越冬期劃分為3個時期,即越冬前期、中期和后期[21]. 于2014年12月-2015年4月,采用非損傷采樣技術采集安徽省升金湖(30.25°~30.50°N, 116.92°~117.25°E)越冬白頭鶴的糞便樣品. 按照3個越冬時期采樣,即:2014年11月上旬-2014年12月下旬,越冬前期;2015年1月上旬-2015年2月下旬,越冬中期;2015年3月上旬-2015年4月上旬,越冬后期. 前期和中期采樣時間間隔為34天,中期和后期的采樣時間間隔為55天(表1).

    在采集糞便樣品前,用雙筒望遠鏡和單筒望遠鏡觀察白頭鶴集中覓食區(qū)域,選擇集群較大的群體,觀察確保在白頭鶴覓食期間周圍50 m范圍內(nèi)沒有其他鶴以及雁鴨類存在或者活動. 等到白頭鶴取食結束離開覓食地后,立刻前往覓食地點進行采集,根據(jù)白頭鶴留下的足跡或覓食坑,確定為白頭鶴覓食地點. 戴一次性消毒手套并用鑷子采集白頭鶴新鮮糞便樣品,去除接觸地面的部分后,迅速將樣品放入無菌50 ml離心管中,并貼好標簽,放在冰上運回實驗室,置于-80℃超低溫冰箱冷凍保存(表1).

    表1 3個越冬時期白頭鶴糞便樣采集的詳細信息

    1.2 DNA的提取

    本實驗采用德國QIAamp Fast DNA Stool Mini Kit 糞便微生物基因組DNA提取試劑盒來提取糞便樣品中的細菌基因組DNA,提取方法嚴格按照試劑盒制造商的具體操作步驟進行,DNA提取樣品保存在 -20℃超低溫冰箱中以備后續(xù)實驗使用.

    1.3 16S rRNA基因擴增及高通量測序

    選擇通用引物338F (5′-GGACTACHVGGGTWTCTAAT-3′)和806R(5′-GGACTACHVGGGTWTCTAAT-3′)作為這次實驗的擴增引物[22],對細菌16S rRNA V3~V4區(qū)進行PCR擴增. 實驗采用20 μl反應體系:10×ex Buffer,2 μl; 0.25 mmol/L dNTPs, 2 μl;正反向引物(20 mmol/L)各0.5 μl;DNA模板,10 ng;exTaq DNA Polymerase, 0.5 μl; ddH2O,12.9 μl. PCR反應循環(huán)程序設置:94℃預變性3 min;95℃變性30 s;55℃退火30 s;72℃延伸45 s,反應27個循環(huán)后72℃終端延伸10 min,終止于4℃. PCR 產(chǎn)物經(jīng)2% 瓊脂糖凝膠電泳鑒定. 將符合標準的樣本送至上海美吉生物醫(yī)藥科技有限公司,利用Illumina Miseq測序平臺進行高通量測序.

    1.4 數(shù)據(jù)分析

    高通量得到的原始數(shù)據(jù)都是以Fastq文件格式保存,首先對原始測序數(shù)據(jù)進行質(zhì)量控制 (去除含有模糊堿基、含有6個以上單堿基高重復區(qū)、長度小于200 bp的序列),獲得純凈序列;通過COPE( Connecting Overlapped Pair-Eend, V1.2.1)軟件對序列進行拼接,并去除上、下游引物,最終得到優(yōu)化數(shù)據(jù). 運用USEARCH(Vsesion 7.1http://drive5.com/uparse/)按照97%相似性對非重復序列(不含單序列)進行OTU(Operational Taxonomic Units)歸類操作并進行OTU豐度統(tǒng)計,運用Mothur進行Rarefaction分析,利用R語言工具繪制稀釋性曲線圖(Rarefaction Curve)[23]. 采用RDP Classifier貝葉斯算法對97%相似水平的OTU代表序列進行分類學分析,并在各個分類水平統(tǒng)計每個樣品的群落組成,置信度閾值為0.7[24]. 應用Mothur分析軟件,在97% OTU水平上計算Alpha多樣性指數(shù)(包括物種豐富度的兩個指數(shù)(Chao1和Ace))、辛普森(Simpson)指數(shù)和香農(nóng)威納(Shannon-Wiener)指數(shù)[25]. 運用R語言軟件繪制主坐標分析圖(Principal Co-ordinates Analysis,PCoA),研究樣本群落組成的相似性. 應用SPSS 19.0軟件對樣品數(shù)據(jù)進行統(tǒng)計學分析,P<0.05時為具有顯著性差異.

    2 結果與分析

    2.1 細菌的Alpha多樣性

    9個樣本共得到337447條有效序列,平均每個樣品測得 37494條有效序列(表2). 序列長度在428~ 447 bp之間,符合16S rRNA 基因V3~V4區(qū)域及引物設計長度. 序列經(jīng)過聚類后共得到708個OTU. 采用對序列進行隨機抽樣的方法,以抽到的序列數(shù)與它們所能代表OTU的數(shù)目構建稀釋性曲線(圖1),用來說明樣本的測序數(shù)據(jù)量是否合理. 圖1結果顯示,在測序量比較小時,隨著測序量的增加OTU數(shù)量會顯著增加,而隨著測序量的逐漸增大新增OTU的數(shù)量逐漸變少,曲線趨于平坦,說明本研究中的測序數(shù)據(jù)量合理,更多的數(shù)據(jù)量只會產(chǎn)生少量新的OTU. 通過方差分析對3個時期OTU數(shù)量進行統(tǒng)計分析,得到F2,6=0.034,P=0.967,表明基于本文的數(shù)據(jù),3個時期OTU數(shù)量未檢測出顯著性差異. 而通過取3個時期OTU數(shù)量平均值,越冬前期為323個,越冬中期和越冬后期分別為298和302個,說明在越冬前期糞便中的微生物豐度最高. 表2為越冬前中后期白頭鶴腸道微生物的Alpha多樣性指標,覆蓋率都在99%以上,為進一步了解升金湖越冬前中后期白頭鶴腸道菌群豐度和多樣性變化,通過方差分析對各個Alpha多樣性指標進行統(tǒng)計分析,得到Ace (F2,6=0.434,P=0.667)、Chao1 (F2,6=0.243,P=0.792)、香農(nóng)威納指數(shù)(F2,6=0.095,P=0.911)和辛普森指數(shù)(F2,6=0.229,P=0.802),表明基于本文的數(shù)據(jù),并未檢測出顯著性差異.

    表2 3個越冬時期白頭鶴腸道菌群的多樣性及豐富度指數(shù)

    2.2 主坐標分析

    本文運用主坐標分析形象化地展示3個時期糞便微生物組成的差異性. 結果表明,同一個越冬階段的糞便樣品微生物菌群組成相似性較高,而不同越冬階段的糞便樣品微生物菌群組成差異性較高(圖2).

    圖1 稀釋性曲線(大寫字母代表樣品編號,下同)Fig.1 Rarefaction curve

    圖2 主坐標聚類分析Fig.2 Principal co-ordinates analysis

    2.3 微生物菌落組成

    2.3.1 門分類水平分析 9個升金湖白頭鶴糞便樣品共檢測出26個門. 越冬前期共鑒定到20個門,優(yōu)勢菌門為厚壁菌門(Firmicutes, 14.2%)、藍藻細菌門(Cyanobacteria, 34.5%)、擬桿菌門(Bacteroidetes, 28.2%);越冬中期共鑒定到15個門,優(yōu)勢菌門包括厚壁菌門(40.4%)、放線菌門(Actinobacteria, 44.2%)和變形菌門(Proteobacteria, 10%);越冬后期共鑒定到18個門,優(yōu)勢菌門包括厚壁菌門(43.2%)和變形菌門(31%).

    在3個時期都檢測出的門有12個,其中厚壁菌門、變形菌門、放線菌門、擬桿菌門和藍藻細菌門存在于所有樣品中,并且作為優(yōu)勢菌,在總序列中的比例達到92.3%. 越冬前期厚壁菌門較中后期含量較低,而擬桿菌門(28.2%) 相對越冬中期(3.6%)和越冬后期(0.4%)含量較高;越冬中期放線菌門(44.2%) 較越冬前期(8.4%)和越冬后期(3.2%)含量高;越冬后期變形菌門含量(31%) 較越冬前期(11.2%)和越冬中期(10%)高,并且越冬后期未能鑒定的門含量較高,明顯高于越冬前期(0.05%)和越冬中期(0.1%),表明越冬后期樣品中存在大量未知的新類群(圖3A ).

    2.3.2 屬分類水平分析 9個升金湖越冬白頭鶴糞便樣品共檢測到252個屬,主要屬在9個樣品中的分布見表3. 從圖3B可以看出3個時期糞便樣品菌群在屬水平上的分布,越冬前期共鑒定到201個屬,優(yōu)勢菌群為普氏菌屬,平均比例為24.3%,要明顯高于越冬中期(2.1%)和越冬后期(0.01%). 并且其含量從越冬前期到越冬中期再到越冬后期有明顯的降低趨勢(圖4A),在越冬后期的3個樣品中只存在于樣品N中. 越冬中期共鑒定到175個屬,優(yōu)勢菌為桿菌屬(20.3%)、節(jié)細菌屬(18.3%)和冷細菌屬(15%). 越冬后期共鑒定到182個屬,占優(yōu)勢的菌群為腸桿菌屬(24.4%)和梭菌屬(10.6%),其中梭菌屬含量要顯著高于越冬前期(P=0.046),明顯高于越冬中期(2%),并且從越冬前期到越冬中期再到越冬后期有明顯的升高趨勢(圖4B). 另外細菌纖維菌屬(Cellulomonas)和Cellulosilyticum在白頭鶴腸道微生物中的平均含量也呈現(xiàn)出越冬后期(0.1%、0.4%)明顯大于越冬前期(0.02%、0.05%),這兩種細菌也可以降解食物中的纖維素.

    表3 升金湖越冬白頭鶴糞便微生物屬分類水平的分布

    圖3 3個越冬時期白頭鶴腸道菌群的門(A)和屬(B)組成Fig.3 Gut microbiota composition in Hooded Cranes at the phylum (A) and genus (B) taxonomical levels in the three wintering periods

    圖4 普氏菌屬(A)和梭菌屬(B)在越冬前、中、后期的含量Fig.4 Contents of Prevotella (A) and Clostridium (B) in the three wintering periods

    3 討論

    3個時期升金湖越冬白頭鶴腸道微生物菌群結構特征均表現(xiàn)出含有較高的腸道微生物多樣性,越冬前、中、后期腸道菌群OTU數(shù)量以及多樣性指數(shù)并沒有顯著性差異. 但越冬前期OTU數(shù)量平均值最高,并且檢測到的門和屬最多,說明在越冬前期檢測到的細菌豐度及細菌種類最多. 通過主坐標分析,越冬前、后期樣品分別表現(xiàn)出相似性,表明越冬前、中、后期白頭鶴腸道菌群組成可能存在差異性. 動物的生存環(huán)境,特別是食性變化影響其腸道微生物菌群結構[26],以前的研究顯示越冬期間白頭鶴覓食棲息地會發(fā)生變化,升金湖越冬白頭鶴在越冬前期主要以稻田生境為覓食棲息地,越冬后期主要以草灘生境為覓食棲息地[18],這兩個時期食物相對單一,可能造成白頭鶴腸道菌群組成分別在這兩個時期表現(xiàn)出一定的相似性,而越冬中期作為兩個時期的過渡期,生境異質(zhì)性較大,食物多樣化,可能造成白頭鶴腸道菌群構成在越冬中期相對于越冬前期和越冬后期差異性大.

    升金湖白頭鶴腸道微生物在門分類水平上主要由厚壁菌門、變形菌門、放線菌門、擬桿菌門和藍藻細菌門組成,與豬、馬、牛等哺乳動物以及麝雉、火雞、家禽等鳥類相似[27-30]. 屬分類水平上,共檢測到252個屬,以普氏菌屬、腸桿菌屬、桿菌屬、節(jié)細菌屬、梭菌屬和冷桿菌屬最為常見,這與豬的腸道微生物在屬水平上以放線桿菌屬、嗜血桿菌屬、巴氏桿菌屬、消化鏈球菌屬、梭桿菌屬、韋永氏球菌屬為主有差異,而與以草為主要食性的馬在腸道微生物分布上相似性較高,都含有較高的普氏菌屬和梭菌屬,并且都含有纖維菌屬.

    不同越冬階段下白頭鶴腸道菌群組成分析表明,在越冬前期,擬桿菌門作為其優(yōu)勢菌,其中普氏菌屬比例較高,有關研究表明這種菌群結構在人類和其他動物的腸道微生物菌群結構中普遍出現(xiàn)[31-32]. 這種菌群結構的普遍存在可能與其功能有關,根據(jù)以前的研究,擬桿菌門主要參與碳水化合物發(fā)酵以及多糖、膽汁酸和類固醇代謝[33],而普氏菌屬作為擬桿菌門中的重要細菌,主要參與腸道中碳水化合物以及植物蛋白質(zhì)的分解代謝[34],稻谷中含有大量的植物蛋白質(zhì)、淀粉和脂肪,并且在以谷物等食性為主的人群的腸道菌群中含量較高[35-38]. 3個時期越冬白頭鶴腸道菌群組成分析表明,普氏菌屬作為越冬前期的優(yōu)勢菌,高于越冬中期和后期,并且隨著越冬時間呈現(xiàn)降低趨勢;而具有降解腸道中纖維素作用的梭菌屬[13],在越冬后期含量較高,顯著高于越冬前期,明顯高于越冬中期. 這可能與升金湖越冬白頭鶴在越冬期間覓食棲息地及食性變化有關:越冬前期白頭鶴以稻谷為主要食性,越冬中期對稻谷利用率降低,在越冬后期稻田生境受到長時間的覓食消耗和春耕等原因,食物利用率最低[17];而隨著越冬時間的推移,在越冬后期,進入早春階段,植被快速生長并進入花期,白頭鶴逐漸把草灘作為主要覓食地[39],食物中可能含有較高的纖維素. 通過對2011年11月—2012年4月升金湖越冬白頭鶴糞便進行顯微鏡觀察,不同越冬期白頭鶴食物組成不同:前期為稻谷、蓼子草(51.96%、16.20%);中期為稻谷、朝天委陵菜(30.94%、17.16%);后期以朝天委陵菜、稻谷(46.74%、14.31%). 其中稻谷、蓼子草隨著越冬期的變化相對密度逐漸減小,朝天委陵菜逐漸增大[21]. 這一結果表明,白頭鶴越冬不同階段食物資源的變化可能對其腸道菌群組成造成了一定的影響,這種影響機制后期還需要進一步深入研究.

    4 結論

    本文對不同越冬階段下的白頭鶴腸道菌群組成做了初步探討,9個樣本得到337447條有效序列,經(jīng)聚類后共得到708個OTU. 白頭鶴在不同越冬階段其腸道菌群組成不同,越冬前期擬桿菌門含量較高,厚壁菌門含量較中后期低;越冬中期放線菌門含量較高;越冬后期變形菌門含量較高. 在屬分類水平,越冬前期以普氏菌屬為優(yōu)勢菌;越冬后期以梭菌屬為優(yōu)勢菌,不同越冬階段下的食物變化可能作為影響白頭鶴腸道微生物菌群結構的一個重要因素,后期還需要進一步研究其對白頭鶴腸道微生物菌群組成的影響機制.

    [1] Bj?rkstén B. The gut microbiota: A complex ecosystem.Clinical&ExperimentalAllergyJournaloftheBritishSocietyforAllergy&ClinicalImmunology, 2006, 36: 1215-1217.

    [2] Ezenwa VO, Gerardo NM, Inouye DWetal. Animal behavior and the microbiome.Science, 2012, 338: 198-199.

    [3] Cani PD, Delzenne NM. The role of the gut microbiota in energy metabolism and metabolic disease.CurrentPharmaceuticalDesign, 2009, 15(13): 1546-1558.

    [4] Young VB. The intestinal microbiota in health and disease.CurrentOpinioninGastroenterology, 2012, 28(1): 63-69.

    [5] Bentley R, Meganathan R. Biosynthesis of vitamin K (menaquinone) in bacteria.MicrobiologicalReviews, 1982, 46(3): 241-280.

    [6] Mazmanian SK, Hua LC, Tzianabos AOetal. An immunomodulatory molecule of symbiotic bacteria directs maturation of the host immune system.Cell, 2005, 122(1): 107-118.

    [7] Spor A, Koren O, Ley R. Unravelling the effects of the environment and host genotype on the gut microbiome.NatureReviewsMicrobiology, 2011, 9(4): 279-290.

    [8] Semova I, Carten J, Stombaugh Jetal. Microbiota regulate intestinal absorption and metabolism of fatty acids in the zebra fish.CellHost&Microbe, 2012, 12(3): 277-288.

    [9] Ley RE, Knight R, Gordon JIetal. Evolution of mammals and their gut microbes.Science, 2008, 320(5883): 1647-1651.

    [10] Claesson MJ, Jeffery IB, Conde Setal. Gut microbiota composition correlates with diet and health in the elderly.Nature, 2012, 488(7410): 178-184.

    [11] Zhang Jiachao. Mongolians core gut microbiota and its correlation with dietary changes[Dissertation]. Huhhot: Inner Mongolia Agricultural University, 2014.[張家超. 蒙古族腸道核心菌群及其與飲食關聯(lián)性研究[學位論文]. 呼和浩特: 內(nèi)蒙古農(nóng)業(yè)大學, 2014.]

    [12] Carlotta DF,Duccio C, Monica DPetal. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa.ProceedingsoftheNationalAcademyofSciences, 2011, 43(33): 445-446.

    [13] Zhu LF, Wu Q, Dai JYetal. Evidence of cellulose metabolism by the giant panda gut microbiome.ProceedingsoftheNationalAcademyofSciencesoftheUnitedStatesofAmerica, 2011, 108(43): 17714-17719.

    [14] Bird Life International.Grusmonacha. In: IUCN 2014. The IUCN red list of threatened species. Version 2014.3. Http://www.iucnredlist.org.

    [15] Yang L, Zhou LZ, Song YW. The effects of food abundance and disturbance on foraging flock patterns of the wintering Hooded Crane (Grusmonacha).AvianResearch, 2015, 6: 15.

    [16] Pan Yunfen, Xu Qing, Cheng Yuanqietal. Herbaceous seed elora of wetlands in Shengjin Lake Nature Reserve of Anhui Province.WetlandScience, 2008, 6: 304-308.[潘云芬, 徐慶, 程元啟等. 安徽升金湖自然保護區(qū)濕地草本種子植物區(qū)系研究. 濕地科學, 2008, 6: 304-308.]

    [17] Liu Zhengyuan, Xu Wenbin, Wang Qishanetal. Environmental carrying capacity for over-wintering Hooded Cranes in Shengjin Lake.ResourcesandEnvironmentintheYangtzeBasin, 2001, 10(5): 454-459.[劉政源,徐文斌,王岐山等. 白頭鶴在升金湖上湖越冬期環(huán)境容納量的研究.長江流域資源與環(huán)境, 2001, 10(5): 454-459.]

    [18] Zheng M, Zhou LZ, Zhao NNetal. Effects of variation in food resources on foraging habitat use by wintering Hooded Cranes (Grusmonacha).AvianResearch, 2015, 6: 11.

    [19] Bishop Mary Anne, Li Fengshan. Effects of farming practices in Tibet on wintering Black-necked Crane (Grusnigricollis) diet and food availability.ChineseBiodiversity, 2002, 10(4): 393-398.[Bishop Mary Anne, 李鳳山. 農(nóng)業(yè)耕作活動對西藏越冬黑頸鶴食性及食物可獲得性的影響. 生物多樣性, 2002, 10(4): 393-398.]

    [20] Zhou Lulu. Seasonal shifts of food habit of the Hooded Crane (Grusmonacha) wintering in the lakes of Yangtze River Floodplain in Anhui Province[Dissertation]. Hefei: Anhui University, 2013.[周璐璐. 安徽沿江湖泊越冬白頭鶴(Grusmonacha)食性的季節(jié)變化[學位論文]. 合肥: 安徽大學, 2013.]

    [21] Zhou B, Zhou L, Chen Jetal. Diurnal time-activity budgets of wintering Hooded Cranes (Grusmonacha) in Shengjin Lake, China.Waterbirds, 2010, 33(1): 110-115.

    [22] Dennis KL, Yunwei W, Blatner NRetal. Adenomatous polyps are driven by microbe-instigated focal inflammation and are controlled by il-10-producing T cells.CancerResearch, 2013, 73(19): 5905-5913.

    [23] Amato KR, Yeoman CJ, Kent Aetal. Habitat degradation impacts black howler monkey (Alouattapigra) gastrointestinal microbiomes.IsmeJournal, 2013, 7(7): 1344-1353.

    [24] Wang Q, Garrity GM, Tiedje JMetal. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy.ApplEnvironMicrob, 2007, 73(16): 5261-5267.

    [25] Schloss PD, Dirk G, Westcott SL. Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies.PLoSOne, 2011, 6(12): e27310.

    [26] Susanne M, Katiana S, Christiana Hetal. Differences in fecal microbiota in different european study populations in relation to age, gender, and country: A cross-sectional study.Applied&EnvironmentalMicrobiology, 2006, 72(2): 1027-1033.

    [27] Lowe BA, Marsh TL, Isaacs-Cosgrove Netal. Defining the "core microbiome" of the microbial communities in the tonsils of healthy pigs.BmcMicrobiology, 2012, 12(1): 20.

    [28] Donnell MMO, Harris HMB, Jeffery IBetal. The core faecal bacterial microbiome of Irish thoroughbred racehorses.LettersinAppliedMicrobiology, 2013, 57(6): 492-501.

    [29] Lu J, Domingo JS. Turkey fecal microbial community structure and functional gene diversity revealed by 16S rRNA gene and metagenomic sequences.JournalofMicrobiology, 2008, 46(5): 469-477.

    [30] Godoy-Vitorino F, Goldfarb KC, Karaoz Uetal. Comparative analyses of foregut and hindgut bacterial communities in hoatzins and cows.IsmeJournal, 2012, 6(3): 531-541.

    [31] Turnbaugh PJ, Micah H, Tanya Yetal. A core gut microbiome in obese and lean twins.Nature, 2009, 457(7228): 480-484.

    [32] Ley RE, Fredrik B, Peter Tetal. Obesity alters gut microbial ecology.ProceedingsoftheNationalAcademyofSciencesoftheUnitedStatesofAmerica, 2005, 102(31): 11070-11075.

    [33] Salyers AA. Bacteroides of the human lower intestinal tract.AnnualReviewofMicrobiology, 1984, 38(1): 293-313.

    [34] Meimandipour A, Hair-Bejo M, Shuhaimi Metal. Gastrointestinal tract morphological alteration by unpleasant physical treatment and modulating role of lactobacillus in broilers.BritishPoultryScience, 2010, 51(1): 52-59.

    [35] Maslowski KM, Vieira AT, Aylwin Netal. Regulation of inflammatory responses by gut microbiota and chemoattractant receptor gpr43.Nature, 2009, 461(7268): 1282-1286.

    [36] Matsuki T, Watanabe K, Fujimoto Jetal. Development of 16S rRNA - gene-targeted group-specific primers for the detection and identification of predominant bacteria in human feces.AppliedandEnvironmentalMicrobiology, 2002, 68(11): 5445-5451.

    [37] Matsumoto K, Takada T, Shimizu Ketal. Effects of a probiotic fermented milk beverage containing lactobacillus casei strain shirota on defecation frequency, intestinal microbiota, and the intestinal environment of healthy individuals with soft stools(brewing and food technology).JournalofBioscience&Bioengineering, 2010, 110(5): 547-552.

    [38] Mcnulty NP, Tanya Y, Ansel Hetal. The impact of a consortium of fermented milk strains on the gut microbiome of gnotobiotic mice and monozygotic twins.ScienceTranslationalMedicine, 2011, 3(3): 106ra106.

    [39] Zhang DM, Zhou LZ, Song YW. Effect of water level fluctuations on temporal-spatial patterns of foraging activities by the wintering Hooded Crane (Grusmonacha).AvianResearch, 2015, 6:16.

    The temporal variations of gut microbiota composition in overwintering Hooded Crane (Grusmonacha) at Lake Shengjin, Anhui Province

    ZHAO Guanghong1, ZHOU Lizhi1**, DONG Yuanqiu1& SONG Yunwei1,2

    (1:SchoolofResourcesandEnvironmentalEngineering,AnhuiUniversity,Hefei230601,P.R.China) (2:LakeShengjinNationalNatureReserveofAnhuiProvince,Chizhou247200,P.R.China)

    The gastrointestinal microbiome plays a significant role in maintaining the stability of gastrointestinal environment of animals, and at the same time the composition of gastrointestinal microbial communities is also influenced by external environment such as diets. Studies on the gut microbiota structure can better understand the ecological conditions of animals. In this study, a total of 9 faeces samples were collected noninvasively from overwintering Hooded Cranes (Grusmonacha) at Lake Shengjin, Anhui Province. Through high-throughput sequencing of V3-V4 regions of 16S rRNA of microbiome, there are 337447 reads and 708 Operational Taxonomic Units (OTUs). The OTU abundance reached the highest in the early wintering period, clustering 586 OTUs with average 323 OTUs. A total of 20 bacterial phyla and 201 genera were identified in the early wintering period, representing 175 genera from 15 bacterial phyla in the mid-wintering period. A total of 18 bacterial phyla with 182 genera were identified in the late wintering period. The results showed that the microbial diversity index such as Shannon-Wiener index and Simpson index had no significant differences in fecal samples of different wintering periods at Lake Shengjin. However, microbial composition had certain differences among three periods. At the microbial phylum level, the phyla Firmicutes (14.2%), Cyanobacteria (34.5%), Bacteroidetes (28.2%) were predominated in the early wintering period, and the phyla Firmicutes (40.4%), Actinobacteria (44.2%) and Proteobacteria (10%) were predominated in the mid-wintering period, while the phyla Firmicutes (43.2%) and Proteobacteria (31%) were predominated in the late wintering period. At the genus level, compared with the other two wintering periods, the amount ofClostridiumwas higher in the late wintering period and the amount ofPrevotelladecreased distinctly from the early to the middle-late wintering periods, while the amount ofClostridiumincreased distinctly from the early to the middle-late wintering periods. The variation of composition of gastrointestinal microbial may be explained by changes of diets of the Hooded Cranes in the wintering period. The study shed some light on overwintering ecology of Hooded Cranes.

    Hooded Cranes; gut microbiota composition; wintering periods; 16S rRNA; Lake Shengjin

    *國家自然科學基金項目(31472020, 31172117)和安徽大學研究生創(chuàng)新項目(YQH100288)聯(lián)合資助. 2016-05-17收稿; 2016-09-09 收修改稿. 趙廣紅(1989~),女,碩士研究生;E-mail: zhaoguanghongabc@163.com.

    ;E-mail: zhoulz@ahu.edu.cn.

    猜你喜歡
    菌門越冬菌群
    “云雀”還是“貓頭鷹”可能取決于腸道菌群
    中老年保健(2022年2期)2022-08-24 03:20:50
    奶牛安全越冬要抓好防寒保溫
    “水土不服”和腸道菌群
    科學(2020年4期)2020-11-26 08:27:06
    野生樹鼩與人工飼養(yǎng)樹鼩消化道不同部位微生物組成的比較研究
    饑餓與重攝食對河蟹腸道菌群結構的影響
    昆蟲體內(nèi)微生物多樣性的影響因素研究進展
    妊娠期糖尿病腸道菌群變化及臨床價值分析
    北方園林植物常見越冬防寒措施
    流蘇越冬種條低溫貯藏扦插繁育試驗初報
    冬季西葫蘆越冬豐產(chǎn)栽培技術
    中文字幕人妻丝袜一区二区| 中文字幕另类日韩欧美亚洲嫩草| 99精品久久久久人妻精品| 国产男女内射视频| 久久精品亚洲熟妇少妇任你| 涩涩av久久男人的天堂| 中文字幕另类日韩欧美亚洲嫩草| 日韩大片免费观看网站| 99国产精品一区二区三区| 亚洲av成人不卡在线观看播放网 | 狠狠婷婷综合久久久久久88av| 操出白浆在线播放| 亚洲精品久久午夜乱码| 飞空精品影院首页| 黄色一级大片看看| 精品熟女少妇八av免费久了| 国产免费视频播放在线视频| 黄色 视频免费看| 国产真人三级小视频在线观看| 叶爱在线成人免费视频播放| 极品人妻少妇av视频| 日日夜夜操网爽| 亚洲精品成人av观看孕妇| 成年动漫av网址| 国产福利在线免费观看视频| 午夜激情久久久久久久| 老司机亚洲免费影院| 97在线人人人人妻| 51午夜福利影视在线观看| 老司机午夜十八禁免费视频| 最近手机中文字幕大全| 日本欧美视频一区| 久久免费观看电影| 伦理电影免费视频| 国产精品久久久久久精品古装| 久久av网站| 亚洲精品久久久久久婷婷小说| 欧美国产精品一级二级三级| 亚洲精品成人av观看孕妇| 亚洲人成网站在线观看播放| 免费高清在线观看视频在线观看| 久久久久久久精品精品| 成人国产一区最新在线观看 | 不卡av一区二区三区| 日本91视频免费播放| 久久国产精品大桥未久av| 国产极品粉嫩免费观看在线| 日韩免费高清中文字幕av| 日本五十路高清| 国产男女超爽视频在线观看| 可以免费在线观看a视频的电影网站| 一本综合久久免费| 在线观看免费视频网站a站| 亚洲一码二码三码区别大吗| 精品一区在线观看国产| 黄色 视频免费看| 看免费成人av毛片| 色婷婷av一区二区三区视频| 中文字幕精品免费在线观看视频| 国产精品香港三级国产av潘金莲 | 999精品在线视频| 国产精品二区激情视频| 又大又黄又爽视频免费| 欧美 亚洲 国产 日韩一| 妹子高潮喷水视频| 999精品在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲伊人久久精品综合| xxx大片免费视频| 中国国产av一级| 精品高清国产在线一区| 又黄又粗又硬又大视频| 美女国产高潮福利片在线看| 亚洲熟女毛片儿| 极品少妇高潮喷水抽搐| 成年动漫av网址| 中文字幕人妻丝袜一区二区| 久久精品国产亚洲av高清一级| 两性夫妻黄色片| 国产精品亚洲av一区麻豆| 久久青草综合色| 一级黄色大片毛片| 在线天堂中文资源库| 欧美日韩一级在线毛片| 色婷婷久久久亚洲欧美| 亚洲,欧美精品.| 亚洲天堂av无毛| 亚洲色图综合在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 日韩一本色道免费dvd| 欧美少妇被猛烈插入视频| 国产日韩欧美视频二区| 亚洲中文av在线| 亚洲国产av新网站| 亚洲午夜精品一区,二区,三区| 日韩欧美一区视频在线观看| 午夜av观看不卡| 亚洲,一卡二卡三卡| 成年av动漫网址| 亚洲人成77777在线视频| 99久久精品国产亚洲精品| 亚洲精品国产av成人精品| 亚洲欧美色中文字幕在线| 男女下面插进去视频免费观看| 欧美中文综合在线视频| 99国产精品一区二区三区| 欧美在线黄色| 18在线观看网站| 50天的宝宝边吃奶边哭怎么回事| 午夜精品国产一区二区电影| 久久人妻熟女aⅴ| 美女扒开内裤让男人捅视频| 欧美成人精品欧美一级黄| 国产亚洲欧美在线一区二区| 色综合欧美亚洲国产小说| 精品福利观看| 大香蕉久久网| 国产成人啪精品午夜网站| 秋霞在线观看毛片| 国产精品二区激情视频| 国产一区亚洲一区在线观看| 国产精品久久久久久精品古装| 视频区欧美日本亚洲| 欧美黄色淫秽网站| 亚洲成人免费电影在线观看 | 丰满人妻熟妇乱又伦精品不卡| 国产精品国产三级国产专区5o| 欧美人与善性xxx| 国产极品粉嫩免费观看在线| 亚洲欧美一区二区三区国产| videosex国产| svipshipincom国产片| 韩国高清视频一区二区三区| 国产一区二区激情短视频 | 超碰成人久久| 深夜精品福利| 各种免费的搞黄视频| 777久久人妻少妇嫩草av网站| 亚洲,欧美精品.| 亚洲人成电影免费在线| 好男人视频免费观看在线| 天天影视国产精品| 久9热在线精品视频| 一级毛片我不卡| h视频一区二区三区| 丝袜人妻中文字幕| 女人爽到高潮嗷嗷叫在线视频| 欧美97在线视频| 18在线观看网站| 久久国产亚洲av麻豆专区| 精品亚洲成国产av| 国产精品 国内视频| 久久综合国产亚洲精品| av天堂在线播放| 交换朋友夫妻互换小说| 丝袜美足系列| 少妇被粗大的猛进出69影院| 亚洲伊人久久精品综合| 亚洲精品在线美女| 亚洲欧美一区二区三区黑人| 亚洲 欧美一区二区三区| 91成人精品电影| 伊人亚洲综合成人网| 人人澡人人妻人| 黑人巨大精品欧美一区二区蜜桃| 视频在线观看一区二区三区| 一本久久精品| 捣出白浆h1v1| 国产亚洲午夜精品一区二区久久| 99久久人妻综合| 在线看a的网站| 人人澡人人妻人| 国产亚洲精品第一综合不卡| 两性夫妻黄色片| 一区二区三区乱码不卡18| 考比视频在线观看| e午夜精品久久久久久久| 亚洲欧美精品综合一区二区三区| 人体艺术视频欧美日本| 青青草视频在线视频观看| 亚洲av在线观看美女高潮| 人成视频在线观看免费观看| 50天的宝宝边吃奶边哭怎么回事| 一区福利在线观看| 日本wwww免费看| 亚洲午夜精品一区,二区,三区| 99国产精品免费福利视频| 操美女的视频在线观看| 国产精品99久久99久久久不卡| 亚洲精品美女久久av网站| 国产精品亚洲av一区麻豆| 亚洲成人免费电影在线观看 | 一级黄色大片毛片| 免费观看av网站的网址| 久久 成人 亚洲| 国产又爽黄色视频| 一区在线观看完整版| 免费在线观看视频国产中文字幕亚洲 | 成年人午夜在线观看视频| 国产精品二区激情视频| 国产亚洲午夜精品一区二区久久| 日本五十路高清| 亚洲av国产av综合av卡| 午夜福利视频精品| 亚洲精品国产区一区二| 80岁老熟妇乱子伦牲交| 香蕉丝袜av| 久久免费观看电影| 多毛熟女@视频| 亚洲欧美中文字幕日韩二区| 欧美av亚洲av综合av国产av| 黑丝袜美女国产一区| 中文字幕人妻熟女乱码| 一级片'在线观看视频| 99久久99久久久精品蜜桃| 亚洲 国产 在线| 亚洲视频免费观看视频| 夫妻性生交免费视频一级片| 婷婷色av中文字幕| 午夜影院在线不卡| 欧美日韩亚洲综合一区二区三区_| 亚洲国产av新网站| 国产精品久久久久成人av| 男女午夜视频在线观看| 久久精品人人爽人人爽视色| 国产成人一区二区三区免费视频网站 | 国产福利在线免费观看视频| 亚洲欧美成人综合另类久久久| 看免费av毛片| 亚洲午夜精品一区,二区,三区| 亚洲国产成人一精品久久久| 日韩精品免费视频一区二区三区| 欧美日韩亚洲综合一区二区三区_| 又大又黄又爽视频免费| 啦啦啦中文免费视频观看日本| 狂野欧美激情性bbbbbb| 精品亚洲成国产av| 捣出白浆h1v1| 香蕉国产在线看| 久久午夜综合久久蜜桃| 黄色视频在线播放观看不卡| 美女午夜性视频免费| 日本a在线网址| xxx大片免费视频| 精品福利观看| 欧美日韩综合久久久久久| 中文字幕最新亚洲高清| avwww免费| 国产麻豆69| 欧美成人午夜精品| 大型av网站在线播放| 脱女人内裤的视频| 久久av网站| 高清av免费在线| www.自偷自拍.com| 婷婷色综合大香蕉| 亚洲精品在线美女| 欧美日韩成人在线一区二区| 欧美精品亚洲一区二区| 亚洲欧美一区二区三区黑人| 一本久久精品| 色网站视频免费| 青青草视频在线视频观看| 国产淫语在线视频| 精品久久蜜臀av无| av在线app专区| 免费观看a级毛片全部| 手机成人av网站| 国产精品一区二区在线观看99| 欧美大码av| 免费少妇av软件| 国产成人av教育| 一边摸一边做爽爽视频免费| 亚洲,一卡二卡三卡| 国产一区有黄有色的免费视频| 免费在线观看黄色视频的| 亚洲欧美激情在线| 最新的欧美精品一区二区| av不卡在线播放| 久久女婷五月综合色啪小说| 国产av国产精品国产| 国产午夜精品一二区理论片| 欧美日韩亚洲国产一区二区在线观看 | 欧美人与善性xxx| 夫妻午夜视频| 欧美精品高潮呻吟av久久| 日韩伦理黄色片| 亚洲熟女毛片儿| 你懂的网址亚洲精品在线观看| 中文精品一卡2卡3卡4更新| 欧美日韩视频高清一区二区三区二| 脱女人内裤的视频| 亚洲av在线观看美女高潮| 国产男女超爽视频在线观看| √禁漫天堂资源中文www| 欧美成人精品欧美一级黄| 国产伦人伦偷精品视频| √禁漫天堂资源中文www| 狂野欧美激情性bbbbbb| 99re6热这里在线精品视频| 国产成人影院久久av| 中国国产av一级| 欧美亚洲日本最大视频资源| 黄色一级大片看看| 欧美黄色淫秽网站| 亚洲欧美精品综合一区二区三区| 91精品三级在线观看| 免费不卡黄色视频| 国产黄色视频一区二区在线观看| 亚洲国产精品999| 国产成人av教育| 亚洲一区中文字幕在线| 欧美黑人欧美精品刺激| 两人在一起打扑克的视频| 午夜激情av网站| 美女高潮到喷水免费观看| 大码成人一级视频| 一边摸一边做爽爽视频免费| 女性生殖器流出的白浆| 在线av久久热| 一二三四在线观看免费中文在| 国产一区二区激情短视频 | 色播在线永久视频| 亚洲av日韩精品久久久久久密 | 考比视频在线观看| 男女国产视频网站| 亚洲精品一二三| 午夜久久久在线观看| 国产成人系列免费观看| 亚洲国产欧美在线一区| 天天影视国产精品| 2018国产大陆天天弄谢| 国产淫语在线视频| 欧美国产精品一级二级三级| 热99久久久久精品小说推荐| 亚洲av电影在线观看一区二区三区| 精品一品国产午夜福利视频| 黑人猛操日本美女一级片| 久久精品熟女亚洲av麻豆精品| 久久久久国产一级毛片高清牌| 日本91视频免费播放| 99热网站在线观看| videos熟女内射| 国产成人一区二区三区免费视频网站 | 免费看不卡的av| 夫妻午夜视频| 91麻豆精品激情在线观看国产 | 亚洲欧美日韩另类电影网站| 亚洲人成电影观看| 亚洲精品一卡2卡三卡4卡5卡 | av福利片在线| 中文字幕亚洲精品专区| 日韩伦理黄色片| 在线亚洲精品国产二区图片欧美| 在线观看免费午夜福利视频| 自线自在国产av| 亚洲人成电影免费在线| 老熟女久久久| 一本—道久久a久久精品蜜桃钙片| 十八禁高潮呻吟视频| 日韩一卡2卡3卡4卡2021年| 亚洲精品一区蜜桃| 爱豆传媒免费全集在线观看| 日韩视频在线欧美| 免费av中文字幕在线| 成人午夜精彩视频在线观看| 成人18禁高潮啪啪吃奶动态图| 国产一卡二卡三卡精品| 亚洲欧洲精品一区二区精品久久久| 国产精品国产av在线观看| 多毛熟女@视频| 国产av一区二区精品久久| 免费久久久久久久精品成人欧美视频| av欧美777| 欧美国产精品一级二级三级| 尾随美女入室| 亚洲精品av麻豆狂野| www.av在线官网国产| 一区二区三区四区激情视频| 午夜福利在线免费观看网站| 精品久久久精品久久久| 国产福利在线免费观看视频| 久久国产精品影院| 在线观看免费日韩欧美大片| 老汉色∧v一级毛片| 观看av在线不卡| 久久久久久久大尺度免费视频| 日本午夜av视频| 狠狠精品人妻久久久久久综合| 精品熟女少妇八av免费久了| 亚洲少妇的诱惑av| 狂野欧美激情性xxxx| 精品国产一区二区三区四区第35| 黑人猛操日本美女一级片| kizo精华| 亚洲精品久久午夜乱码| 国产xxxxx性猛交| 观看av在线不卡| 欧美人与性动交α欧美软件| 极品少妇高潮喷水抽搐| 叶爱在线成人免费视频播放| 国产av国产精品国产| 亚洲精品久久成人aⅴ小说| 亚洲午夜精品一区,二区,三区| 亚洲精品美女久久av网站| 2018国产大陆天天弄谢| av线在线观看网站| 欧美精品一区二区免费开放| 手机成人av网站| 色网站视频免费| 亚洲色图综合在线观看| 蜜桃在线观看..| 亚洲欧美一区二区三区国产| 两个人免费观看高清视频| 性少妇av在线| 建设人人有责人人尽责人人享有的| 国产福利在线免费观看视频| 久久热在线av| 热re99久久精品国产66热6| 免费女性裸体啪啪无遮挡网站| 啦啦啦在线观看免费高清www| a 毛片基地| 男女之事视频高清在线观看 | 涩涩av久久男人的天堂| 2018国产大陆天天弄谢| 久久亚洲国产成人精品v| 亚洲精品久久午夜乱码| 色网站视频免费| 在现免费观看毛片| 国产成人影院久久av| 国产成人欧美| 亚洲免费av在线视频| 国产又爽黄色视频| 久久ye,这里只有精品| 午夜福利,免费看| 啦啦啦视频在线资源免费观看| 久久精品人人爽人人爽视色| 免费日韩欧美在线观看| 亚洲一区二区三区欧美精品| 老司机靠b影院| bbb黄色大片| 无遮挡黄片免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 成人午夜精彩视频在线观看| 一级毛片女人18水好多 | 777久久人妻少妇嫩草av网站| 国产成人影院久久av| 午夜视频精品福利| 亚洲精品乱久久久久久| 国产欧美日韩一区二区三 | 老司机深夜福利视频在线观看 | 黑人猛操日本美女一级片| 国产成人一区二区在线| 国产免费一区二区三区四区乱码| 久久99一区二区三区| 国产日韩欧美视频二区| 国产成人av教育| 在线观看一区二区三区激情| 国产熟女午夜一区二区三区| 18禁观看日本| 色综合欧美亚洲国产小说| 蜜桃在线观看..| 丁香六月欧美| 少妇裸体淫交视频免费看高清 | 18禁观看日本| av又黄又爽大尺度在线免费看| 欧美变态另类bdsm刘玥| 亚洲精品av麻豆狂野| 国产亚洲精品久久久久5区| 亚洲欧美一区二区三区国产| 99精品久久久久人妻精品| 在线观看人妻少妇| 国产高清国产精品国产三级| 国产免费现黄频在线看| 欧美日韩av久久| 免费人妻精品一区二区三区视频| 国产爽快片一区二区三区| 亚洲av男天堂| 国产又色又爽无遮挡免| 中文字幕高清在线视频| 男女之事视频高清在线观看 | 欧美在线一区亚洲| 蜜桃国产av成人99| 亚洲精品乱久久久久久| 别揉我奶头~嗯~啊~动态视频 | 侵犯人妻中文字幕一二三四区| 亚洲少妇的诱惑av| www.999成人在线观看| 亚洲精品久久午夜乱码| 看十八女毛片水多多多| 99九九在线精品视频| 免费观看av网站的网址| 亚洲精品日本国产第一区| 亚洲欧美精品自产自拍| 亚洲国产精品成人久久小说| 欧美在线一区亚洲| 欧美日韩国产mv在线观看视频| 久久女婷五月综合色啪小说| 欧美日韩国产mv在线观看视频| 日本黄色日本黄色录像| 欧美大码av| 国产91精品成人一区二区三区 | 男女之事视频高清在线观看 | 欧美日本中文国产一区发布| 久久国产亚洲av麻豆专区| 国产成人精品在线电影| 国产免费一区二区三区四区乱码| 亚洲九九香蕉| 亚洲自偷自拍图片 自拍| 最近中文字幕2019免费版| 亚洲色图综合在线观看| 亚洲国产日韩一区二区| 高潮久久久久久久久久久不卡| 十八禁网站网址无遮挡| 成人免费观看视频高清| 亚洲国产欧美在线一区| av在线老鸭窝| 国产麻豆69| 好男人电影高清在线观看| 欧美精品一区二区免费开放| 激情视频va一区二区三区| 日本av免费视频播放| 国产黄色视频一区二区在线观看| 国产熟女欧美一区二区| 日日摸夜夜添夜夜爱| 捣出白浆h1v1| 精品欧美一区二区三区在线| videos熟女内射| 在线av久久热| 精品高清国产在线一区| 亚洲三区欧美一区| 国产成人精品在线电影| 午夜福利视频精品| 成人免费观看视频高清| 国语对白做爰xxxⅹ性视频网站| 午夜福利在线免费观看网站| 精品卡一卡二卡四卡免费| 午夜影院在线不卡| 久久久久视频综合| 日韩精品免费视频一区二区三区| 男女国产视频网站| 国产av一区二区精品久久| 中文字幕最新亚洲高清| 中文字幕人妻丝袜一区二区| 黄色毛片三级朝国网站| 在线精品无人区一区二区三| 国产成人av激情在线播放| 精品一区二区三卡| www.精华液| 久久天堂一区二区三区四区| 成人影院久久| 久久亚洲国产成人精品v| 久久国产精品人妻蜜桃| 看免费成人av毛片| 中文字幕av电影在线播放| 亚洲中文av在线| 咕卡用的链子| 欧美激情高清一区二区三区| 精品卡一卡二卡四卡免费| 国产福利在线免费观看视频| 亚洲成人免费电影在线观看 | 午夜福利影视在线免费观看| 精品一区在线观看国产| 欧美av亚洲av综合av国产av| 啦啦啦在线观看免费高清www| 成人午夜精彩视频在线观看| 免费久久久久久久精品成人欧美视频| 久久精品熟女亚洲av麻豆精品| 久久综合国产亚洲精品| 一本一本久久a久久精品综合妖精| 国产成人av激情在线播放| 国产男女超爽视频在线观看| 亚洲av片天天在线观看| 黄频高清免费视频| 曰老女人黄片| 18在线观看网站| 久久久久久免费高清国产稀缺| 成人黄色视频免费在线看| 十八禁高潮呻吟视频| 国产精品欧美亚洲77777| 国产精品久久久人人做人人爽| 一边亲一边摸免费视频| 99国产精品免费福利视频| 高清不卡的av网站| 天天添夜夜摸| 亚洲精品日本国产第一区| 亚洲精品久久久久久婷婷小说| 欧美日韩黄片免| 电影成人av| 国产无遮挡羞羞视频在线观看| 日韩一卡2卡3卡4卡2021年| 在线观看免费视频网站a站| 男女国产视频网站| 精品人妻熟女毛片av久久网站| 看免费av毛片| 久久精品人人爽人人爽视色| 91九色精品人成在线观看| 天堂中文最新版在线下载| 亚洲欧美精品综合一区二区三区| 欧美精品啪啪一区二区三区 | 999久久久国产精品视频| e午夜精品久久久久久久| 搡老岳熟女国产| 欧美另类一区| 男女无遮挡免费网站观看| 国产人伦9x9x在线观看| 黄色视频不卡| 国产1区2区3区精品| 日韩熟女老妇一区二区性免费视频| 老司机深夜福利视频在线观看 | 欧美日韩福利视频一区二区| tube8黄色片| 男女国产视频网站| 女人久久www免费人成看片| 在现免费观看毛片|