• <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)栽培技術
    91大片在线观看| 91字幕亚洲| 午夜福利,免费看| 男女下面插进去视频免费观看| 中国美女看黄片| 国产成人一区二区三区免费视频网站| 丝袜美足系列| 91成人精品电影| 国产精品日韩av在线免费观看 | 欧美日韩精品网址| 午夜a级毛片| 免费女性裸体啪啪无遮挡网站| 国产精品秋霞免费鲁丝片| 夜夜躁狠狠躁天天躁| 国产视频一区二区在线看| 美女免费视频网站| 日日干狠狠操夜夜爽| 女人精品久久久久毛片| 丝袜美足系列| 国产精品亚洲一级av第二区| www.自偷自拍.com| 波多野结衣高清无吗| 亚洲国产精品999在线| 午夜福利18| 久热这里只有精品99| 18禁国产床啪视频网站| 久久中文看片网| 久久久水蜜桃国产精品网| 18禁国产床啪视频网站| 99热只有精品国产| 国产午夜精品久久久久久| 久久久水蜜桃国产精品网| 亚洲av电影在线进入| 午夜免费观看网址| 国产蜜桃级精品一区二区三区| 久久精品影院6| av视频在线观看入口| 如日韩欧美国产精品一区二区三区| 久久久久精品国产欧美久久久| 欧美性长视频在线观看| 啦啦啦免费观看视频1| 日本免费a在线| 最新在线观看一区二区三区| 岛国在线观看网站| 少妇被粗大的猛进出69影院| 桃红色精品国产亚洲av| 精品第一国产精品| 老熟妇仑乱视频hdxx| 88av欧美| 禁无遮挡网站| 亚洲精品av麻豆狂野| 每晚都被弄得嗷嗷叫到高潮| 国产成年人精品一区二区| 男人的好看免费观看在线视频 | 成年女人毛片免费观看观看9| 久久精品亚洲精品国产色婷小说| 男人舔女人下体高潮全视频| 日韩大码丰满熟妇| 成人欧美大片| 国产野战对白在线观看| 啪啪无遮挡十八禁网站| 国产日韩一区二区三区精品不卡| 悠悠久久av| 亚洲精品美女久久av网站| 制服诱惑二区| 亚洲国产欧美网| 国产av在哪里看| av在线天堂中文字幕| 欧美色欧美亚洲另类二区 | 午夜福利在线观看吧| 色尼玛亚洲综合影院| 国产精品一区二区精品视频观看| 国产精品二区激情视频| 午夜福利高清视频| 丝袜美腿诱惑在线| 人妻久久中文字幕网| 99国产精品一区二区三区| 日本欧美视频一区| 丝袜人妻中文字幕| 中文字幕av电影在线播放| 男人舔女人的私密视频| 男女之事视频高清在线观看| 校园春色视频在线观看| 中文亚洲av片在线观看爽| 亚洲,欧美精品.| 两性夫妻黄色片| 50天的宝宝边吃奶边哭怎么回事| 亚洲伊人色综图| 亚洲九九香蕉| 午夜成年电影在线免费观看| 日本黄色视频三级网站网址| 亚洲无线在线观看| 亚洲avbb在线观看| 国产一区二区三区视频了| 久久精品影院6| 日本一区二区免费在线视频| 不卡av一区二区三区| 国产又色又爽无遮挡免费看| 91国产中文字幕| 欧美最黄视频在线播放免费| 一级a爱片免费观看的视频| 最近最新中文字幕大全免费视频| 在线播放国产精品三级| 欧美一区二区精品小视频在线| 琪琪午夜伦伦电影理论片6080| 三级毛片av免费| 国产熟女xx| 美女免费视频网站| 在线观看午夜福利视频| 丝袜美足系列| 搡老岳熟女国产| 啪啪无遮挡十八禁网站| 啦啦啦韩国在线观看视频| 亚洲 国产 在线| 极品教师在线免费播放| 高清在线国产一区| 亚洲国产欧美日韩在线播放| or卡值多少钱| 高清在线国产一区| 亚洲精品在线观看二区| 一夜夜www| 亚洲国产精品久久男人天堂| 亚洲国产高清在线一区二区三 | 精品免费久久久久久久清纯| 午夜成年电影在线免费观看| 999久久久国产精品视频| 黑丝袜美女国产一区| 99久久综合精品五月天人人| 动漫黄色视频在线观看| 欧美激情高清一区二区三区| 精品第一国产精品| 午夜福利18| 久久久久久久久久久久大奶| 日日摸夜夜添夜夜添小说| 国产三级在线视频| 亚洲人成电影观看| 亚洲国产日韩欧美精品在线观看 | 夜夜爽天天搞| 国内精品久久久久精免费| 嫩草影视91久久| 国产伦人伦偷精品视频| 淫妇啪啪啪对白视频| 成人国产一区最新在线观看| 制服丝袜大香蕉在线| 国产成人啪精品午夜网站| 妹子高潮喷水视频| 日本a在线网址| 伦理电影免费视频| 国产野战对白在线观看| 又大又爽又粗| 热re99久久国产66热| 国产成人欧美| 少妇被粗大的猛进出69影院| 久久中文字幕一级| 久久精品成人免费网站| 欧美中文日本在线观看视频| svipshipincom国产片| 国产精品亚洲美女久久久| 999久久久精品免费观看国产| 亚洲自拍偷在线| 精品乱码久久久久久99久播| 久久久久久久精品吃奶| 久久这里只有精品19| 亚洲情色 制服丝袜| 熟妇人妻久久中文字幕3abv| 国产高清视频在线播放一区| 色综合亚洲欧美另类图片| 精品国产乱码久久久久久男人| 国产成人av教育| 欧美乱码精品一区二区三区| 国产精品二区激情视频| 免费不卡黄色视频| 久久久久久免费高清国产稀缺| 久久久久久大精品| 亚洲国产精品久久男人天堂| av中文乱码字幕在线| 啪啪无遮挡十八禁网站| 日韩大尺度精品在线看网址 | 成人18禁在线播放| 国产三级在线视频| 久9热在线精品视频| 一本综合久久免费| 亚洲成av片中文字幕在线观看| 性少妇av在线| 精品国产超薄肉色丝袜足j| 露出奶头的视频| 亚洲男人的天堂狠狠| 欧美亚洲日本最大视频资源| 少妇裸体淫交视频免费看高清 | 久久九九热精品免费| 亚洲av片天天在线观看| 岛国视频午夜一区免费看| 很黄的视频免费| 国产乱人伦免费视频| 黄色视频不卡| 午夜久久久久精精品| 久久久久亚洲av毛片大全| 一级片免费观看大全| 亚洲精品av麻豆狂野| 久久狼人影院| 日本撒尿小便嘘嘘汇集6| 在线观看66精品国产| 国产精品日韩av在线免费观看 | 久久久水蜜桃国产精品网| 可以免费在线观看a视频的电影网站| 国产亚洲精品久久久久久毛片| 嫩草影院精品99| 在线播放国产精品三级| 国产精品久久久av美女十八| bbb黄色大片| 老司机在亚洲福利影院| 黑丝袜美女国产一区| 国产精品一区二区在线不卡| 日本欧美视频一区| 午夜日韩欧美国产| 免费av毛片视频| 操美女的视频在线观看| 欧美亚洲日本最大视频资源| 国产成人精品无人区| 久久国产精品人妻蜜桃| 午夜福利欧美成人| 免费不卡黄色视频| 999久久久精品免费观看国产| 日本免费a在线| 高潮久久久久久久久久久不卡| 丁香欧美五月| 18禁黄网站禁片午夜丰满| 97人妻天天添夜夜摸| 亚洲精品国产精品久久久不卡| 国产单亲对白刺激| 黑人巨大精品欧美一区二区蜜桃| 如日韩欧美国产精品一区二区三区| 天天一区二区日本电影三级 | 熟妇人妻久久中文字幕3abv| 久久久国产精品麻豆| 久久久久精品国产欧美久久久| 免费看a级黄色片| 丝袜人妻中文字幕| 老司机在亚洲福利影院| 国产精品自产拍在线观看55亚洲| 国内精品久久久久精免费| 一a级毛片在线观看| 国产极品粉嫩免费观看在线| 国产午夜福利久久久久久| 亚洲av第一区精品v没综合| 日日夜夜操网爽| 久久婷婷人人爽人人干人人爱 | 一区二区日韩欧美中文字幕| 精品久久蜜臀av无| av天堂在线播放| 男人的好看免费观看在线视频 | 色播在线永久视频| 在线天堂中文资源库| 欧洲精品卡2卡3卡4卡5卡区| 国产亚洲精品第一综合不卡| 非洲黑人性xxxx精品又粗又长| 精品不卡国产一区二区三区| 男女下面插进去视频免费观看| 男男h啪啪无遮挡| 亚洲第一青青草原| 国内久久婷婷六月综合欲色啪| 涩涩av久久男人的天堂| 国产欧美日韩一区二区精品| 亚洲电影在线观看av| 两个人视频免费观看高清| 久久国产亚洲av麻豆专区| 免费av毛片视频| 亚洲第一欧美日韩一区二区三区| 人人妻,人人澡人人爽秒播| 男人的好看免费观看在线视频 | 在线观看舔阴道视频| 国产成年人精品一区二区| 亚洲国产中文字幕在线视频| 一级a爱片免费观看的视频| 国产麻豆成人av免费视频| 亚洲成人免费电影在线观看| 久久天躁狠狠躁夜夜2o2o| 中文字幕精品免费在线观看视频| 国产精品日韩av在线免费观看 | 日韩视频一区二区在线观看| 看片在线看免费视频| 久久精品国产清高在天天线| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产精品sss在线观看| av福利片在线| 麻豆久久精品国产亚洲av| 一级片免费观看大全| 久久精品亚洲精品国产色婷小说| 国产精品自产拍在线观看55亚洲| 日韩大尺度精品在线看网址 | 久久久久久久久中文| 男女床上黄色一级片免费看| 九色国产91popny在线| av欧美777| 欧美中文日本在线观看视频| 村上凉子中文字幕在线| x7x7x7水蜜桃| 欧美+亚洲+日韩+国产| 午夜精品在线福利| 国产99白浆流出| 国产精品久久视频播放| 国产精品1区2区在线观看.| 黄频高清免费视频| 亚洲自偷自拍图片 自拍| 88av欧美| 两人在一起打扑克的视频| 亚洲最大成人中文| 日韩精品免费视频一区二区三区| av视频在线观看入口| 天堂影院成人在线观看| 久久国产精品影院| 午夜福利,免费看| 午夜精品久久久久久毛片777| 成人欧美大片| 天天躁狠狠躁夜夜躁狠狠躁| 一级片免费观看大全| 人妻久久中文字幕网| 久久精品91蜜桃| 久久久久久国产a免费观看| 最新在线观看一区二区三区| 一级a爱视频在线免费观看| 亚洲精品美女久久av网站| 婷婷丁香在线五月| 99国产极品粉嫩在线观看| 在线观看日韩欧美| 纯流量卡能插随身wifi吗| 一区二区三区精品91| av免费在线观看网站| 狠狠狠狠99中文字幕| 亚洲av片天天在线观看| 久久久精品欧美日韩精品| 看片在线看免费视频| 99在线人妻在线中文字幕| 人人澡人人妻人| 国产私拍福利视频在线观看| 国产aⅴ精品一区二区三区波| 日韩精品免费视频一区二区三区| 免费久久久久久久精品成人欧美视频| 午夜a级毛片| 精品日产1卡2卡| 国产亚洲精品第一综合不卡| 欧美中文日本在线观看视频| 亚洲第一av免费看| 色av中文字幕| 国产主播在线观看一区二区| 在线播放国产精品三级| 国产一卡二卡三卡精品| 国产在线精品亚洲第一网站| 亚洲,欧美精品.| 亚洲欧洲精品一区二区精品久久久| 免费少妇av软件| 欧美中文日本在线观看视频| 啦啦啦免费观看视频1| 可以免费在线观看a视频的电影网站| 午夜影院日韩av| 极品人妻少妇av视频| 日韩免费av在线播放| 男人操女人黄网站| 午夜视频精品福利| 国产精品美女特级片免费视频播放器 | 可以在线观看毛片的网站| 岛国视频午夜一区免费看| 九色国产91popny在线| 国产精品免费一区二区三区在线| 黄色a级毛片大全视频| 日本精品一区二区三区蜜桃| 久久青草综合色| 老汉色av国产亚洲站长工具| 看片在线看免费视频| 久久精品影院6| 一级毛片女人18水好多| 色综合站精品国产| 国产熟女xx| 久久中文字幕人妻熟女| 日韩国内少妇激情av| 欧美日韩乱码在线| 成年版毛片免费区| 久久精品成人免费网站| 精品国产超薄肉色丝袜足j| 两性夫妻黄色片| 欧美亚洲日本最大视频资源| 黄色片一级片一级黄色片| 欧美+亚洲+日韩+国产| 亚洲国产欧美日韩在线播放| 久久久国产成人免费| 午夜视频精品福利| 欧美色欧美亚洲另类二区 | www.熟女人妻精品国产| 国产av一区二区精品久久| 亚洲男人的天堂狠狠| 男女做爰动态图高潮gif福利片 | 午夜a级毛片| 日本三级黄在线观看| 男人舔女人下体高潮全视频| 亚洲精品国产色婷婷电影| 他把我摸到了高潮在线观看| 我的亚洲天堂| 久久久国产精品麻豆| 91在线观看av| 亚洲七黄色美女视频| 日韩精品中文字幕看吧| cao死你这个sao货| 亚洲精品一卡2卡三卡4卡5卡| 这个男人来自地球电影免费观看| 国产午夜福利久久久久久| 波多野结衣一区麻豆| 纯流量卡能插随身wifi吗| 久久久精品国产亚洲av高清涩受| 国产野战对白在线观看| 精品久久久久久成人av| 1024香蕉在线观看| 脱女人内裤的视频| 黑丝袜美女国产一区| 欧美绝顶高潮抽搐喷水| 国产精品久久久久久人妻精品电影| 欧美日韩精品网址| 欧美乱妇无乱码| 午夜福利高清视频| 亚洲一区高清亚洲精品| 国产主播在线观看一区二区| 黄片播放在线免费| 色综合站精品国产| 亚洲va日本ⅴa欧美va伊人久久| 国产成人精品无人区| 变态另类成人亚洲欧美熟女 | 男女午夜视频在线观看| 成人亚洲精品av一区二区| 日韩欧美一区二区三区在线观看| 国产在线精品亚洲第一网站| 亚洲无线在线观看| 国产精品乱码一区二三区的特点 | 一区二区三区高清视频在线| 很黄的视频免费| 美女大奶头视频| 色综合欧美亚洲国产小说| 久久亚洲精品不卡| 欧美一级a爱片免费观看看 | 长腿黑丝高跟| 好看av亚洲va欧美ⅴa在| 高清黄色对白视频在线免费看| 一二三四在线观看免费中文在| 在线观看免费日韩欧美大片| 亚洲全国av大片| 女人精品久久久久毛片| 午夜福利一区二区在线看| 国产在线精品亚洲第一网站| or卡值多少钱| 亚洲色图av天堂| 大香蕉久久成人网| 97碰自拍视频| 亚洲avbb在线观看| 国产精品一区二区三区四区久久 | 亚洲av电影在线进入| 精品久久久久久久久久免费视频| 99精品欧美一区二区三区四区| 国产人伦9x9x在线观看| 国产又爽黄色视频| 亚洲人成电影观看| 国产亚洲精品久久久久5区| 国内久久婷婷六月综合欲色啪| 国产一区在线观看成人免费| 91麻豆av在线| 啦啦啦免费观看视频1| 欧美成人免费av一区二区三区| 99国产极品粉嫩在线观看| 久9热在线精品视频| 每晚都被弄得嗷嗷叫到高潮| 久久性视频一级片| 老司机在亚洲福利影院| а√天堂www在线а√下载| 好看av亚洲va欧美ⅴa在| 亚洲性夜色夜夜综合| 久久 成人 亚洲| 国产午夜精品久久久久久| 国产99久久九九免费精品| 国产91精品成人一区二区三区| 国内久久婷婷六月综合欲色啪| 日韩欧美一区视频在线观看| 中出人妻视频一区二区| 一区二区三区高清视频在线| 午夜影院日韩av| 可以在线观看的亚洲视频| 国产一区二区三区视频了| 亚洲国产欧美一区二区综合| 女人爽到高潮嗷嗷叫在线视频| 淫秽高清视频在线观看| 成人精品一区二区免费| 国产黄a三级三级三级人| 亚洲第一电影网av| 亚洲片人在线观看| 午夜成年电影在线免费观看| 一区福利在线观看| 精品久久久久久成人av| 午夜a级毛片| 亚洲avbb在线观看| 午夜a级毛片| 国产精品久久电影中文字幕| 一区福利在线观看| 日韩欧美在线二视频| 在线观看一区二区三区| 成熟少妇高潮喷水视频| 日韩精品免费视频一区二区三区| 亚洲av电影不卡..在线观看| 免费少妇av软件| 老司机福利观看| 中文字幕人成人乱码亚洲影| 精品欧美一区二区三区在线| 国产真人三级小视频在线观看| www国产在线视频色| 亚洲国产精品999在线| 精品欧美一区二区三区在线| 国产xxxxx性猛交| 亚洲男人的天堂狠狠| 色综合亚洲欧美另类图片| 国产精品一区二区在线不卡| 99热只有精品国产| 日韩大尺度精品在线看网址 | 欧美日韩亚洲国产一区二区在线观看| avwww免费| 欧美国产日韩亚洲一区| 亚洲激情在线av| 91老司机精品| 亚洲熟妇熟女久久| 久久久久久久久中文| 香蕉丝袜av| 岛国在线观看网站| 国产精品久久视频播放| 激情视频va一区二区三区| 免费在线观看黄色视频的| 亚洲少妇的诱惑av| 91麻豆av在线| 国产蜜桃级精品一区二区三区| 两个人看的免费小视频| 99国产精品免费福利视频| 脱女人内裤的视频| 久久久国产精品麻豆| 国产亚洲精品av在线| 最近最新中文字幕大全免费视频| 免费在线观看影片大全网站| 人妻久久中文字幕网| 国产黄a三级三级三级人| 亚洲午夜精品一区,二区,三区| 亚洲人成电影免费在线| 日韩欧美一区视频在线观看| 丝袜美足系列| 日韩高清综合在线| 久久久国产欧美日韩av| 欧美一级a爱片免费观看看 | 国产精品自产拍在线观看55亚洲| 久久香蕉精品热| 女人被狂操c到高潮| 99国产精品免费福利视频| 少妇熟女aⅴ在线视频| 99精品久久久久人妻精品| 精品一品国产午夜福利视频| 中文字幕高清在线视频| 国产成人精品无人区| 亚洲专区中文字幕在线| 色尼玛亚洲综合影院| 日韩视频一区二区在线观看| 一区二区三区激情视频| 久久久久国产一级毛片高清牌| 99在线人妻在线中文字幕| 亚洲av五月六月丁香网| 免费在线观看亚洲国产| 88av欧美| 美女国产高潮福利片在线看| 好男人在线观看高清免费视频 | 一级片免费观看大全| 久久欧美精品欧美久久欧美| 制服人妻中文乱码| 中文字幕色久视频| 19禁男女啪啪无遮挡网站| av免费在线观看网站| 99久久精品国产亚洲精品| 美女高潮喷水抽搐中文字幕| 久久精品亚洲熟妇少妇任你| 18禁黄网站禁片午夜丰满| 亚洲激情在线av| 免费在线观看视频国产中文字幕亚洲| av福利片在线| 精品国产美女av久久久久小说| 一本综合久久免费| 99精品久久久久人妻精品| 国产精品久久久人人做人人爽| 色播在线永久视频| 久久九九热精品免费| 男女下面进入的视频免费午夜 | 欧美精品亚洲一区二区| 欧美色欧美亚洲另类二区 | 一本大道久久a久久精品| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜福利一区二区在线看| 午夜福利影视在线免费观看| 久久香蕉激情| 在线永久观看黄色视频| 搞女人的毛片| 国产麻豆成人av免费视频| 亚洲成人国产一区在线观看| 天天一区二区日本电影三级 | 搡老岳熟女国产| 成人精品一区二区免费| 久久性视频一级片| 一区福利在线观看| 午夜福利影视在线免费观看| 精品熟女少妇八av免费久了| 午夜免费观看网址| 亚洲国产精品合色在线| 性少妇av在线| 99在线人妻在线中文字幕| 精品少妇一区二区三区视频日本电影| 亚洲第一欧美日韩一区二区三区| 久久久国产精品麻豆| 国产一级毛片七仙女欲春2 | 日韩中文字幕欧美一区二区|