• <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)栽培技術
    在线观看国产h片| 国产大屁股一区二区在线视频| 中国美白少妇内射xxxbb| 国产老妇女一区| 成人毛片a级毛片在线播放| 国产精品一及| 欧美xxxx黑人xx丫x性爽| 女人被狂操c到高潮| 亚洲人成网站高清观看| 人人妻人人澡人人爽人人夜夜| 男人舔奶头视频| 一级片'在线观看视频| 天天躁日日操中文字幕| 婷婷色av中文字幕| 国产一区二区三区综合在线观看 | 日韩av不卡免费在线播放| 亚洲精品乱久久久久久| 好男人在线观看高清免费视频| 欧美成人午夜免费资源| 成人亚洲精品一区在线观看 | 久久久久久伊人网av| 免费观看无遮挡的男女| 婷婷色av中文字幕| 精品酒店卫生间| 国产精品无大码| 日韩强制内射视频| 久久久久国产精品人妻一区二区| 联通29元200g的流量卡| 午夜福利高清视频| 日韩强制内射视频| 97精品久久久久久久久久精品| 国产日韩欧美在线精品| 1000部很黄的大片| 色哟哟·www| 国产精品嫩草影院av在线观看| 乱码一卡2卡4卡精品| 精品酒店卫生间| 3wmmmm亚洲av在线观看| 欧美zozozo另类| 婷婷色av中文字幕| 精品少妇黑人巨大在线播放| 欧美日韩精品成人综合77777| 99九九线精品视频在线观看视频| 国产一级毛片在线| 女人久久www免费人成看片| 日本色播在线视频| 少妇裸体淫交视频免费看高清| 免费在线观看成人毛片| 亚洲在久久综合| 免费高清在线观看视频在线观看| 国国产精品蜜臀av免费| 狂野欧美白嫩少妇大欣赏| 青春草视频在线免费观看| 天堂中文最新版在线下载 | 在线亚洲精品国产二区图片欧美 | 国内精品宾馆在线| 欧美97在线视频| 国产午夜福利久久久久久| 五月玫瑰六月丁香| 水蜜桃什么品种好| 亚洲在久久综合| 日本色播在线视频| 能在线免费看毛片的网站| 亚洲不卡免费看| 国产精品.久久久| 97人妻精品一区二区三区麻豆| 青青草视频在线视频观看| 麻豆精品久久久久久蜜桃| 另类亚洲欧美激情| 蜜臀久久99精品久久宅男| 综合色av麻豆| 丝袜脚勾引网站| 国产淫片久久久久久久久| 国产爱豆传媒在线观看| 网址你懂的国产日韩在线| 男的添女的下面高潮视频| 欧美日韩在线观看h| 亚洲美女搞黄在线观看| 极品少妇高潮喷水抽搐| 波多野结衣巨乳人妻| 亚洲欧美一区二区三区黑人 | 久久久久久久久久人人人人人人| 亚洲第一区二区三区不卡| 午夜福利在线观看免费完整高清在| 日韩一本色道免费dvd| 婷婷色综合大香蕉| 国产国拍精品亚洲av在线观看| 亚洲色图综合在线观看| 中文乱码字字幕精品一区二区三区| 成人二区视频| 中文在线观看免费www的网站| 国产精品人妻久久久久久| 三级男女做爰猛烈吃奶摸视频| 亚洲精品视频女| 少妇人妻久久综合中文| 美女xxoo啪啪120秒动态图| 99久国产av精品国产电影| 国产v大片淫在线免费观看| 男人添女人高潮全过程视频| 国产毛片a区久久久久| 久久久亚洲精品成人影院| 国产精品国产三级国产av玫瑰| 国产综合精华液| 老司机影院成人| 国产精品久久久久久精品古装| 亚洲最大成人av| 精品久久久久久电影网| 久久久精品94久久精品| 男插女下体视频免费在线播放| 久久久久久久大尺度免费视频| 少妇的逼水好多| 99久久九九国产精品国产免费| 久久精品国产亚洲av涩爱| 中国美白少妇内射xxxbb| 高清在线视频一区二区三区| 中文资源天堂在线| 18禁动态无遮挡网站| 免费黄频网站在线观看国产| 国产精品人妻久久久久久| 色综合色国产| 国内精品美女久久久久久| 国产有黄有色有爽视频| 日韩欧美精品v在线| 亚洲欧洲国产日韩| 啦啦啦在线观看免费高清www| 深爱激情五月婷婷| 国产 精品1| 真实男女啪啪啪动态图| 在线精品无人区一区二区三 | 欧美少妇被猛烈插入视频| 午夜福利视频1000在线观看| 国产爽快片一区二区三区| 99久久人妻综合| 日日摸夜夜添夜夜添av毛片| 天天一区二区日本电影三级| 看十八女毛片水多多多| 两个人的视频大全免费| 亚洲精品日韩av片在线观看| 高清av免费在线| 黄色怎么调成土黄色| 一本一本综合久久| 成人鲁丝片一二三区免费| 日本一本二区三区精品| 人妻一区二区av| av天堂中文字幕网| 欧美xxⅹ黑人| 亚洲怡红院男人天堂| 国产高清有码在线观看视频| 国产精品国产三级国产专区5o| 人妻一区二区av| 国产精品av视频在线免费观看| 男插女下体视频免费在线播放| 人妻 亚洲 视频| 一区二区三区精品91| 日日摸夜夜添夜夜爱| 亚洲精品国产av蜜桃| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品自拍成人| 男女啪啪激烈高潮av片| 亚洲精品自拍成人| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产av国产精品国产| 插阴视频在线观看视频| 亚洲四区av| 亚洲欧美中文字幕日韩二区| 国产乱人偷精品视频| 精品久久久久久久人妻蜜臀av| 欧美+日韩+精品| 成年女人看的毛片在线观看| 3wmmmm亚洲av在线观看| 日本欧美国产在线视频| 伦理电影大哥的女人| av又黄又爽大尺度在线免费看| 97精品久久久久久久久久精品| www.av在线官网国产| 免费不卡的大黄色大毛片视频在线观看| 狂野欧美激情性bbbbbb| 国产精品一区二区在线观看99| 在线观看三级黄色| 在线亚洲精品国产二区图片欧美 | 日韩一区二区三区影片| 成年女人看的毛片在线观看| 一区二区三区四区激情视频| 美女被艹到高潮喷水动态| 中文欧美无线码| 王馨瑶露胸无遮挡在线观看| 在线观看一区二区三区激情| 又大又黄又爽视频免费| 亚洲人与动物交配视频| 色综合色国产| 美女xxoo啪啪120秒动态图| 新久久久久国产一级毛片| 22中文网久久字幕| 国语对白做爰xxxⅹ性视频网站| 国产欧美日韩精品一区二区| 国产成人免费观看mmmm| 亚洲精品中文字幕在线视频 | 免费少妇av软件| 久久99精品国语久久久| 大香蕉97超碰在线| 国产一区二区三区综合在线观看 | 国产精品久久久久久精品电影小说 | 免费av观看视频| 国内揄拍国产精品人妻在线| 在现免费观看毛片| 伊人久久精品亚洲午夜| 夫妻午夜视频| 新久久久久国产一级毛片| 国产极品天堂在线| 国产亚洲精品久久久com| 国产伦在线观看视频一区| 美女国产视频在线观看| 国产精品伦人一区二区| 22中文网久久字幕| 99久久精品国产国产毛片| 少妇的逼水好多| 菩萨蛮人人尽说江南好唐韦庄| 精品人妻熟女av久视频| 欧美日韩综合久久久久久| 久久影院123| 日日啪夜夜爽| 少妇人妻一区二区三区视频| 秋霞在线观看毛片| 国产午夜福利久久久久久| 日日啪夜夜撸| 六月丁香七月| 高清视频免费观看一区二区| 欧美精品国产亚洲| 别揉我奶头 嗯啊视频| 亚洲精品亚洲一区二区| 国产精品嫩草影院av在线观看| 国产欧美亚洲国产| 深爱激情五月婷婷| 久久久a久久爽久久v久久| 国产精品一区二区性色av| 亚洲内射少妇av| 最近2019中文字幕mv第一页| 人人妻人人爽人人添夜夜欢视频 | 亚洲va在线va天堂va国产| 日本一本二区三区精品| 99视频精品全部免费 在线| www.色视频.com| 国产精品久久久久久av不卡| 久久精品久久久久久久性| 少妇的逼水好多| 寂寞人妻少妇视频99o| 美女国产视频在线观看| 深爱激情五月婷婷| 婷婷色av中文字幕| 波多野结衣巨乳人妻| .国产精品久久| 亚洲欧美日韩无卡精品| 国产精品精品国产色婷婷| 日韩亚洲欧美综合| 久久精品熟女亚洲av麻豆精品| 日韩中字成人| 国产乱来视频区| 天堂中文最新版在线下载 | 国产精品不卡视频一区二区| 国产av不卡久久| 久久久久久久久久久丰满| 国产亚洲5aaaaa淫片| 婷婷色综合www| 黄色一级大片看看| 欧美性感艳星| 日韩在线高清观看一区二区三区| 97在线视频观看| 真实男女啪啪啪动态图| 久久久久久伊人网av| 免费看日本二区| 亚洲av欧美aⅴ国产| 亚洲高清免费不卡视频| 免费大片黄手机在线观看| 欧美另类一区| 青春草视频在线免费观看| www.色视频.com| 91精品国产九色| 天天躁夜夜躁狠狠久久av| 春色校园在线视频观看| av福利片在线观看| 国产精品久久久久久av不卡| 伦理电影大哥的女人| 欧美性感艳星| 久久亚洲国产成人精品v| 九色成人免费人妻av| 乱码一卡2卡4卡精品| 少妇丰满av| 成年免费大片在线观看| 国产黄色免费在线视频| 国产黄色视频一区二区在线观看| 在线看a的网站| 久久人人爽av亚洲精品天堂 | 老师上课跳d突然被开到最大视频| 网址你懂的国产日韩在线| 日日撸夜夜添| 最后的刺客免费高清国语| 美女高潮的动态| 岛国毛片在线播放| 午夜日本视频在线| 国产伦精品一区二区三区视频9| 夫妻午夜视频| 国产亚洲精品久久久com| 高清在线视频一区二区三区| 最近中文字幕2019免费版| 免费不卡的大黄色大毛片视频在线观看| 国语对白做爰xxxⅹ性视频网站| 五月天丁香电影| 男人添女人高潮全过程视频| 一区二区三区精品91| 国产精品一区www在线观看| 欧美区成人在线视频| 亚洲国产精品成人综合色| 亚洲精品一二三| 高清av免费在线| tube8黄色片| 亚洲国产精品专区欧美| 久久精品国产a三级三级三级| 五月伊人婷婷丁香| 亚洲精品456在线播放app| 汤姆久久久久久久影院中文字幕| 久久精品国产亚洲网站| av在线观看视频网站免费| 国产淫片久久久久久久久| 欧美三级亚洲精品| 午夜福利视频1000在线观看| 国产午夜福利久久久久久| 精品久久久精品久久久| 国产精品一区二区三区四区免费观看| 免费观看性生交大片5| 国产成人一区二区在线| av在线天堂中文字幕| 人妻 亚洲 视频| 2018国产大陆天天弄谢| 可以在线观看毛片的网站| 777米奇影视久久| 又爽又黄a免费视频| 亚洲美女视频黄频| 丝瓜视频免费看黄片| 亚洲成人中文字幕在线播放| 免费看日本二区| 2021天堂中文幕一二区在线观| 国产精品99久久99久久久不卡 | 免费看av在线观看网站| 免费观看a级毛片全部| 精品人妻熟女av久视频| 国产白丝娇喘喷水9色精品| 日韩av在线免费看完整版不卡| freevideosex欧美| 日本猛色少妇xxxxx猛交久久| 人妻系列 视频| 成人亚洲欧美一区二区av| 成人国产麻豆网| 禁无遮挡网站| 久久精品国产亚洲av涩爱| 亚洲av免费在线观看| 一级av片app| 我的老师免费观看完整版| 久久影院123| 国产日韩欧美在线精品| 又大又黄又爽视频免费| 狂野欧美白嫩少妇大欣赏| 高清午夜精品一区二区三区| av国产精品久久久久影院| 亚洲人成网站高清观看| 国产高潮美女av| 亚洲色图综合在线观看| 日韩欧美一区视频在线观看 | 视频中文字幕在线观看| 精品人妻一区二区三区麻豆| 成人毛片a级毛片在线播放| 一本一本综合久久| 欧美潮喷喷水| 自拍偷自拍亚洲精品老妇| 91精品伊人久久大香线蕉| 国产亚洲av片在线观看秒播厂| 国产老妇女一区| 日韩强制内射视频| 欧美 日韩 精品 国产| 五月开心婷婷网| 亚洲国产av新网站| 欧美xxxx性猛交bbbb| 亚洲人成网站在线观看播放| 久久99蜜桃精品久久| 久久久精品欧美日韩精品| 精品亚洲乱码少妇综合久久| 国产女主播在线喷水免费视频网站| 免费电影在线观看免费观看| 免费大片18禁| 亚洲天堂国产精品一区在线| 精华霜和精华液先用哪个| 男人爽女人下面视频在线观看| 欧美区成人在线视频| 午夜福利视频1000在线观看| 久久久久久九九精品二区国产| 国产美女午夜福利| 久久久午夜欧美精品| 久久久久网色| 久久久午夜欧美精品| 国产精品一区www在线观看| 亚洲精品自拍成人| 一级av片app| 少妇被粗大猛烈的视频| 国产日韩欧美亚洲二区| 久久久国产一区二区| 人妻 亚洲 视频| 激情 狠狠 欧美| 欧美激情国产日韩精品一区| 欧美激情在线99| 成年版毛片免费区| 99久久精品国产国产毛片| 两个人的视频大全免费| 国产乱来视频区| 深夜a级毛片| 狠狠精品人妻久久久久久综合| 日韩精品有码人妻一区| 真实男女啪啪啪动态图| 夫妻午夜视频| 内射极品少妇av片p| 久久99热这里只频精品6学生| 国产精品一区二区三区四区免费观看| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美精品自产自拍| 秋霞在线观看毛片| eeuss影院久久| 中文字幕亚洲精品专区| 久久久成人免费电影| 亚洲精品乱码久久久v下载方式| 国产成人免费观看mmmm| 亚洲丝袜综合中文字幕| 狂野欧美激情性xxxx在线观看| 熟女av电影| 精品少妇久久久久久888优播| 国产白丝娇喘喷水9色精品| 插阴视频在线观看视频| 天天躁夜夜躁狠狠久久av| 久久99热6这里只有精品| 国产 一区精品| 高清毛片免费看| 亚洲av.av天堂| 大又大粗又爽又黄少妇毛片口| 男女边吃奶边做爰视频| 搞女人的毛片| 能在线免费看毛片的网站| 人妻夜夜爽99麻豆av| 狂野欧美激情性bbbbbb| av国产免费在线观看| 久久精品国产亚洲av天美| 色吧在线观看| 另类亚洲欧美激情| 国产高清有码在线观看视频| 一级毛片黄色毛片免费观看视频| 狂野欧美白嫩少妇大欣赏| 国产真实伦视频高清在线观看| 高清毛片免费看| xxx大片免费视频| 秋霞伦理黄片| 色视频在线一区二区三区| 少妇高潮的动态图| av在线播放精品| 特级一级黄色大片| 国产一区亚洲一区在线观看| 亚洲一级一片aⅴ在线观看| 麻豆精品久久久久久蜜桃| 最近中文字幕2019免费版| 一本久久精品| 亚洲av中文av极速乱| 香蕉精品网在线| 一级毛片 在线播放| 日韩一区二区视频免费看| 欧美日韩视频高清一区二区三区二| 国产白丝娇喘喷水9色精品| 国产亚洲精品久久久com| 亚洲av中文av极速乱| 中文字幕亚洲精品专区| 五月玫瑰六月丁香| 91狼人影院| 王馨瑶露胸无遮挡在线观看| 午夜福利高清视频| 特大巨黑吊av在线直播| 国产精品人妻久久久影院| 91精品一卡2卡3卡4卡| 亚洲国产精品成人综合色| 日本-黄色视频高清免费观看| 最近2019中文字幕mv第一页| 少妇丰满av| 国产成人精品一,二区| 成人特级av手机在线观看| 欧美激情久久久久久爽电影| 春色校园在线视频观看| 只有这里有精品99| 97在线人人人人妻| a级毛片免费高清观看在线播放| 亚洲国产日韩一区二区| 禁无遮挡网站| 国产一区二区三区av在线| av在线app专区| 欧美少妇被猛烈插入视频| www.av在线官网国产| 99久国产av精品国产电影| 免费人妻精品一区二区三区视频| 国产乱人偷精品视频| 免费av中文字幕在线| 成年人免费黄色播放视频| 夜夜骑夜夜射夜夜干| 国产精品国产三级专区第一集| 日韩制服丝袜自拍偷拍| 多毛熟女@视频| 国产成人精品福利久久| 美女国产高潮福利片在线看| 69精品国产乱码久久久| 国产色婷婷99| 观看美女的网站| 99精国产麻豆久久婷婷| 日韩一区二区三区影片| e午夜精品久久久久久久| 热99久久久久精品小说推荐| 国产免费一区二区三区四区乱码| 色精品久久人妻99蜜桃| 亚洲图色成人| 免费av中文字幕在线| 性高湖久久久久久久久免费观看| 欧美日韩av久久| 中文字幕人妻熟女乱码| 日韩精品免费视频一区二区三区| 久久久久久久精品精品| av.在线天堂| 精品国产超薄肉色丝袜足j| 成人国产av品久久久| 69精品国产乱码久久久| 欧美97在线视频| 精品少妇久久久久久888优播| 国产伦人伦偷精品视频| 久久久精品区二区三区| 国产成人免费观看mmmm| 亚洲国产精品999| 欧美日韩国产mv在线观看视频| 1024视频免费在线观看| 高清在线视频一区二区三区| 中国国产av一级| 精品国产乱码久久久久久男人| 建设人人有责人人尽责人人享有的| 亚洲国产av影院在线观看| 久久毛片免费看一区二区三区| 王馨瑶露胸无遮挡在线观看| 国产极品天堂在线| 久久久国产欧美日韩av| 91老司机精品| 国产成人精品久久久久久| 另类精品久久| 亚洲综合色网址| 黄色一级大片看看| 亚洲欧美一区二区三区国产| 日韩伦理黄色片| 人妻 亚洲 视频| 色播在线永久视频| 十八禁人妻一区二区| 久久久国产欧美日韩av| 日韩成人av中文字幕在线观看| 亚洲av日韩在线播放| 在线亚洲精品国产二区图片欧美| 亚洲国产中文字幕在线视频| 国产精品一区二区在线观看99| 欧美成人午夜精品| 少妇的丰满在线观看| 91精品三级在线观看| 亚洲一区中文字幕在线| 亚洲欧洲国产日韩| 18禁国产床啪视频网站| 亚洲综合精品二区| av免费观看日本| 美女中出高潮动态图| 丁香六月欧美| 亚洲av中文av极速乱| 高清在线视频一区二区三区| 青青草视频在线视频观看| 亚洲国产精品一区二区三区在线| 天天躁夜夜躁狠狠躁躁| 久久久精品94久久精品| 国产熟女欧美一区二区| 亚洲国产欧美日韩在线播放| 国产亚洲午夜精品一区二区久久| 久久久久久久久免费视频了| 国产福利在线免费观看视频| 自线自在国产av| 国产片内射在线| 在线天堂最新版资源| 国产精品久久久久久久久免| 2018国产大陆天天弄谢| 午夜福利在线免费观看网站| 在线观看www视频免费| 又粗又硬又长又爽又黄的视频| 国产在线免费精品| 99热全是精品| 一本—道久久a久久精品蜜桃钙片| 日本黄色日本黄色录像| 97精品久久久久久久久久精品| 十分钟在线观看高清视频www| 十八禁网站网址无遮挡| 日韩 亚洲 欧美在线| 日本vs欧美在线观看视频| 亚洲成av片中文字幕在线观看| 曰老女人黄片| 精品一区二区三区四区五区乱码 | 国产男女超爽视频在线观看| 免费看av在线观看网站| www.精华液| 国产男女超爽视频在线观看| 婷婷色av中文字幕| 欧美另类一区| 亚洲第一区二区三区不卡| 久久久久久人妻| 久久人人爽av亚洲精品天堂| 丝袜美腿诱惑在线| 如何舔出高潮| 久久ye,这里只有精品|